An ultraspherical spectral element method for solving partial differential equations
Date
2023-12
Authors
Journal Title
Journal ISSN
Volume Title
Publisher
Stellenbosch : Stellenbosch University
Abstract
ENGLISH ABSTRACT: We investigate the ultraspherical spectral element method for solving second-order partial differential equations in two dimensions. Moreover, a novel coordinate transformation is introduced to broaden the scope of the method, making it applicable to rectangular domains with circular holes (square donuts), as well as certain types of curved boundaries. The presented method is an integration of two approaches, namely the ultraspherical spectral method and the hierarchical Poincaré–Steklov (HPS) scheme. The ultraspherical method is a Petrov–Galerkin scheme that presents operators in the form of sparse and almost-banded matrices, enabling both stability and computational efficiency. The HPS method is a recursive domain decomposition strategy that enables fast direct solves. It merges solution operators and Dirichlet-to-Neumann operators between subdomains, enforcing continuity of the solution and its derivative across domain boundaries. The fusion of these two methods, combined with a bilinear mapping, results in an accurate discretisation with an explicit direct solve that can be applied to problems on arbitrary polygonal domains with smooth solutions. A major advantage is the reuse of precomputed solution operators facilitated by the HPS scheme, enhancing the efficiency of elliptic solves within implicit and semi-implicit time-steppers. Additionally, the approach is highly parallelisable, allowing for efficient computation time. An implementation of the method is established as a software system, ultraSEM, which employs the HPS method to solve on rectangular and polygonal domains. We extend this implementation to allow for solving on domains with circular cavities. This extension relies on a nonlinear coordinate mapping and proves to work effectively, achieving near machine level precision accuracy. On some simple test problems, we demonstrate geometric convergence for refinement of the polynomial degree and algebraic convergence for domain refinement. Furthermore, we show that execution times scale comparably to those achieved for a rectangular domain. We demonstrate the application of the method on various time-dependent and fluid dynamics examples, including contaminant transport and reaction-diffusion systems, and underscore the practical applicability of the methodology and the new domain.
AFRIKAANSE OPSOMMING:Ons ondersoek die ultrasferiese spektrale element metode vir die oplossing van tweedeorde parsiële differensiaalvergelykings in twee dimensies. ’n Nuwe koördinaattransformasie word voorgestel om die omvang van die metode te verbreed en dit van toepassing te maak op reghoekige gebiede met sirkelvormige gate, sowel as sekere tipes gekurfde grense. Die voorgestelde metode is ’n integrasie van twee benaderings, naamlik die ultrasferiese pektrale metode en die hiërargiese Poincaré–Steklov (HPS) skema. Die ultrasferiese metode is ’n Petrov–Galerkin-skema wat operatore skep in die vorm van matrikse wat effektief is om mee te werk en goeie stabiliteit bied. Die HPS-metode is ’n rekursiewe gebied-ontbindingstrategie wat vinnige direkte oplossings moontlik maak. Dit weef oplossings-operatore en Dirichletna- Neumann-operatore tussen subgebiede saam. In die proses word die kontinuïteit van die oplossing en sy afgeleide, oor gebiedsgrense afgedwing. Die samesmelting van hierdie twee strategieë, gekombineer met ’n bilineêre afbeelding, lei tot ’n direkte strategie om akkurate oplossings te genereer vir probleme op veelhoekige gebiede met gladde oplossings. ’n Groot voordeel is die hergebruik van voorafberekende oplossings-operatore wat deur die HPS-skema gefasiliteer word. Dit verbeter veral die doeltreffendheid om elliptiese probleme op te los binne implisiete en semi-implisiete tydstappers. Daarbenewens is dit ook moontlik om die benadering in parallel toe te pas en dus ‘n doeltreffende berekeningstyd moontlik te maak. ’n Reedsbestaande implementering van die metode bestaan wel as ‘n sagtewarestelsel genaamd ultraSEM. Hierdie sagteware gebruik die HPS-metode om oplossings op reghoekige en veelhoekige gebiede te vind. Ons brei hierdie implementering uit om die oplossing op gebiede met sirkelvormige holtes moontlik te maak. Hierdie uitbreiding maak staat op ’n nie-lineêre koördinaat afbeelding en word bewys om effektief te werk. Met ’n paar eenvoudige toetsprobleme word die metode gedemonstreer. Dit vertoon geometriese konvergensie wanneer die polinoomgraad verfyn word en algebraïese konvergensie vir gebiedverfyning. Verder groei die berekeningstyd teen ‘n vergelykbare tempo as dit wat vir ’n vierhoekige gebied behaal word. Ons demonstreer die toepassing van die nuwe metode op verskeie tydafhanklike en vloeidinamika voorbeelde, insluitend kontaminant vervoer en reaksie-diffusie stelsels. Met hierdie voorbeelde word die praktiese toepaslikheid van die metodologie en die nuwe gebied onderstreep.
AFRIKAANSE OPSOMMING:Ons ondersoek die ultrasferiese spektrale element metode vir die oplossing van tweedeorde parsiële differensiaalvergelykings in twee dimensies. ’n Nuwe koördinaattransformasie word voorgestel om die omvang van die metode te verbreed en dit van toepassing te maak op reghoekige gebiede met sirkelvormige gate, sowel as sekere tipes gekurfde grense. Die voorgestelde metode is ’n integrasie van twee benaderings, naamlik die ultrasferiese pektrale metode en die hiërargiese Poincaré–Steklov (HPS) skema. Die ultrasferiese metode is ’n Petrov–Galerkin-skema wat operatore skep in die vorm van matrikse wat effektief is om mee te werk en goeie stabiliteit bied. Die HPS-metode is ’n rekursiewe gebied-ontbindingstrategie wat vinnige direkte oplossings moontlik maak. Dit weef oplossings-operatore en Dirichletna- Neumann-operatore tussen subgebiede saam. In die proses word die kontinuïteit van die oplossing en sy afgeleide, oor gebiedsgrense afgedwing. Die samesmelting van hierdie twee strategieë, gekombineer met ’n bilineêre afbeelding, lei tot ’n direkte strategie om akkurate oplossings te genereer vir probleme op veelhoekige gebiede met gladde oplossings. ’n Groot voordeel is die hergebruik van voorafberekende oplossings-operatore wat deur die HPS-skema gefasiliteer word. Dit verbeter veral die doeltreffendheid om elliptiese probleme op te los binne implisiete en semi-implisiete tydstappers. Daarbenewens is dit ook moontlik om die benadering in parallel toe te pas en dus ‘n doeltreffende berekeningstyd moontlik te maak. ’n Reedsbestaande implementering van die metode bestaan wel as ‘n sagtewarestelsel genaamd ultraSEM. Hierdie sagteware gebruik die HPS-metode om oplossings op reghoekige en veelhoekige gebiede te vind. Ons brei hierdie implementering uit om die oplossing op gebiede met sirkelvormige holtes moontlik te maak. Hierdie uitbreiding maak staat op ’n nie-lineêre koördinaat afbeelding en word bewys om effektief te werk. Met ’n paar eenvoudige toetsprobleme word die metode gedemonstreer. Dit vertoon geometriese konvergensie wanneer die polinoomgraad verfyn word en algebraïese konvergensie vir gebiedverfyning. Verder groei die berekeningstyd teen ‘n vergelykbare tempo as dit wat vir ’n vierhoekige gebied behaal word. Ons demonstreer die toepassing van die nuwe metode op verskeie tydafhanklike en vloeidinamika voorbeelde, insluitend kontaminant vervoer en reaksie-diffusie stelsels. Met hierdie voorbeelde word die praktiese toepaslikheid van die metodologie en die nuwe gebied onderstreep.
Description
Thesis (MSc)--Stellenbosch University, 2023.