Mimetic Spectral Element advection

Technische Samenvatting: Mimetic Spectral Element Advectie

Probleemstelling

Het artikel behandelt de numerieke discretisatie van de lineaire advectievergelijking voor differentiaalvormen op begrensde domeinen. Het klassieke scalaire advectieprobleem, $\partial_t \rho + \nabla \cdot (v\rho) = 0$, wordt binnen een differentiaalgeometrisch kader gegeneraliseerd als:
$$\frac{\partial \alpha^{(k)}}{\partial t} + \mathcal{L}_v \alpha^{(k)} = 0$$
waarbij $\alpha^{(k)}$ een $k$-differentiaalvorm is (die scalairen, vectoren of hogere-dimensionale grootheden vertegenwoordigt, afhankelijk van $k$), $v$ een voorgeschreven Lipschitz-continue snelheidsveld is, en $\mathcal{L}_v$ de Lie-afgeleide is. De uitdaging ligt in het ontwikkelen van een discretisatie die de onderliggende geometrische structuren van de natuurkundige wetten behoudt, lokale massaconservatie waarborgt en spectrale nauwkeurigheid bereikt, terwijl er onderscheid wordt gemaakt tussen grootheden die worden geëvalueerd op specifieke tijdsmomenten en grootheden die over tijdsintervallen worden geïntegreerd.

Methodologie

De voorgestelde methode, Mimetic Spectral Element Advection, breidt het mimetic framework uit dat in eerder werk [4] is vastgesteld door de Lie-afgeleide te incorporeren via de Cartan-homotopieformule. De aanpak is gebouwd op drie kernzuilen:

1. Differentiaalgeometrische Fundering
De methode maakt gebruik van de taal van differentiaalvormen. De Lie-afgeleide $\mathcal{L}_v$ wordt uitgedrukt met behulp van de Cartan-homotopieformule:
$$\mathcal{L}_v \alpha^{(k)} = d\iota_v \alpha^{(k)} + \iota_v d\alpha^{(k)}$$
waarbij $d$ de externe afgeleide is en $\iota_v$ de inwendige product (interior product). Een cruciaal aspect van de formulering is de adjunct-relatie tussen het inwendige product en het wedge-product, gedefinieerd via de $L^2$-inwendige product:
$$(\iota_v \alpha^{(k)}, \beta^{(k-1)})_{L^2} = (\alpha^{(k)}, v^\flat \wedge \beta^{(k-1)})_{L^2}$$
Deze dualiteit maakt het mogelijk om fysieke grootheden die door inwendige producten worden gerepresenteerd, te mappen naar hun duale differentiële 1-vormen.

2. Ruimtelijke Discretisatie (Mimetic Spectral Elements)
Het ruimtelijke domein wordt gediscretiseerd met behulp van een celcomplex bestaande uit punten, lijngedeften, oppervlakken en volumes. De ruimte van gladde $k$-vormen wordt benaderd door een eindige-dimensionale ruimte $\Lambda_h^k$, gespannen door basisvormen $\epsilon_i^{(k)}$.

• Basisconstructie: De basisfuncties worden geconstrueerd met behulp van stuksgewijze polynomiale expansies (tensorproducten) op quadrilaterale elementen. In 1D worden 0-vormen gerepresenteerd door Lagrange-polynomen geëvalueerd op Gauss-Lobatto-knopen, terwijl 1-vormen worden gerepresenteerd door "rand-polynomen" (edge polynomials) die zijn afgeleid om aan specifieke integraaleigenschappen te voldoen.
• Commuterende Projectie: Een projectie-operator $\pi_h$ wordt gedefinieerd zodanig dat deze commuteert met de externe afgeleide ($\pi_h d = d \pi_h$). Dit zorgt ervoor dat de discrete externe afgeleide $d_h$ wordt gerepresenteerd door incidentiematrices die enkel waarden $\{0, 1, -1\}$ bevatten, waardoor de topologische structuur van de continue operatoren behouden blijft.

3. Temporele Discretisatie (Canonieke Mimetic Integrator)
De tijdintegratie maakt gebruik van een willekeurige orde symplectische operator afgeleid van canonieke Gauss-collocatie-integratoren [5].

• Gestaffelde Tijdknopen: De methode maakt onderscheid tussen twee sets tijdknopen:
  • Gauss-Lobatto-knopen ($t_k$): Waar de oplossingsvariabelen (bijv. $\rho$) worden gedefinieerd.
  • Gauss-knopen ($\tilde{t}_q$): Waar de tijdsafgeleiden (fluxen) worden geëvalueerd.
• Discrete Evolutie: Deze stafeling resulteert in een discrete integrator waarbij de verandering in de oplossing over een interval wordt gelijkgesteld aan de flux geëvalueerd op de interne Gauss-knopen. Deze structuur weerspiegelt leap-frog methoden en de impliciete middelpuntregel, wat de symplectische eigenschappen waarborgt.

4. Discrete Inwendige Product
Het inwendige product $\iota_v$ wordt gediscretiseerd door de dualiteitsparing (Vergelijking 16) af te dwingen in de discrete setting. Dit leidt tot een systeem waarbij de discrete fluxen worden berekend door een lineair systeem op te lossen dat de inwendige producten van de basisvormen en het snelheidsveld betreft.

Belangrijkste Bijdragen

• Geometrische Consistentie: Het schema incorporeert expliciet de Lie-afgeleide met behulp van de Cartan-formule, waardoor de discretisatie de metriek-vrije representatie van differentiële operatoren en hun Hilbert-adjuncten respecteert.
• Spectrale Nauwkeurigheid: De methode maakt gebruik van hoog-orde polynomiale basisfuncties, wat spectrale convergentie in de ruimte mogelijk maakt.
• Lokale Massaconservatie: Het gebruik van incidentiematrices en de specifieke structuur van de discrete externe afgeleide garandeert lokale massaconservatie.
• Tijdsstap-onderscheid: Het framework scheidt rigoureus grootheden die worden geëvalueerd op tijdsmomenten van grootheden die over intervallen worden geïntegreerd, waarbij een gestaffeld tijdsgrid wordt gebruikt om hoog-orde symplectische integratie te bereiken.

Numerieke Resultaten

De auteurs presenteren numerieke experimenten op 2D-domeinen met periodieke randvoorwaarden, waarbij de advectie van sinusgolven en sinusbellen in constante en Rudman-vortex snelheidsvelden wordt getest.

• Convergentie: De methode vertoont algebraïsche $h$-convergentie van orde $(p+1)$ en spectrale $p$-convergentie, mits de tijdintegratiefout niet domineert over de ruimtelijke fout.
• Nauwkeurigheid van de Tijdintegratie: De fout in de oplossing blijkt afhankelijk te zijn van de orde van het tijdintegratieschema ($p_t$). Als het tijdschema voldoende nauwkeurig is, wordt de initiële discretiefout behouden; anders neemt kunstmatige dispersie in de loop van de tijd toe.
• Massaconservatie: De totale massa-fout blijft op machine-nul voor de eerste $10^3$ tijdstappen en blijft onder de $10^{-12}$ zelfs na $2 \times 10^4$ stappen, wat uitstekende conserveringseigenschappen aantoont.
• Reversibiliteit: Een test waarbij de advectie van een sinusgolf in een Rudman-vortex wordt gevolgd door een omkering van de stroomrichting, laat zien dat de methode de initiële oplossing kan herstellen, wat de reversibiliteit van de integratiemethode bewijst.
• Kunstmatige Dispersie: Het artikel merkt op dat kunstmatige dispersie een functie is van de orde van de tijdintegratie; lagere-orde tijdschema's introduceren dispersiefouten die afhangen van de frequentie van de geadverteerde golf.

Betekenis en Claims

Het artikel claimt dat het afgeleide schema er succesvol in slaagt om spectrale nauwkeurigheid te combineren met lokale massaconservatie binnen een physics-compatible (mimetic) framework. Door het mimetic framework uit te breiden met de Lie-afgeleide via de Cartan-homotopieformule, bieden de auteurs een discretisatie die de geometrische structuur van de advectievergelijking verheldert. De methode wordt gepresenteerd als een robuuste aanpak voor het oplossen van advectieproblemen waarbij het behoud van de onderliggende geometrische en topologische eigenschappen van de natuurkundige wetten essentieel is. De auteurs merken bescheiden op dat hoewel de ruimtelijke discretisatie zeer nauwkeurig is, de algehele nauwkeurigheid en de dispersiekenmerken afhankelijk zijn van de orde van het gebruikte tijdintegratieschema.