A hybrid $C^{0}$-interior penalty method for the nematic Helmholtz--Korteweg equation

Technische Samenvatting: Een Hybride $C^0$-Interior Penalty Methode voor de Nematicus Helmholtz–Korteweg Vergelijking

Probleemstelling
Dit werk behandelt de numerieke oplossing van de nematicus Helmholtz–Korteweg vergelijking, een vierde-orde partiële differentiaalvergelijking die de voortplanting van tijdsarmonische akoestische golven in nematicus Korteweg-vloeistoffen, zoals nematicus vloeibare kristallen, modelleert. De leidende vergelijking breidt de klassieke Helmholtz-vergelijking uit met twee aanvullende vierde-orde termen, waarvan er één anisotroop is ten opzichte van het nematicus directorveld $\mathbf{n}$. Het probleem is geformuleerd op een begrensde Lipschitz-domein $\Omega \subset \mathbb{R}^d$ ($d=2,3$) met gemengde sound-soft en impedantie randvoorwaarden.

De analyse staat voor drie primaire uitdagingen:

1. Hoge Orde en Anisotropie: De aanwezigheid van de vierde-orde anisotrope term compliceert de analyse in vergelijking met standaard Helmholtz- of biharmonische problemen.
2. Regulariteit en Stabiliteit: Eerdere analyses vertrouwden op perturbatie-argumenten rond het isotrope geval ($\beta=0$), wat vereiste dat de anisotropieparameter $\beta$ klein was ten opzichte van een stabiliteitsconstante van het isotrope probleem. Deze constante degradeert nabij resonantiefrequenties, waardoor $\beta$ kunstmatig naar nul wordt gedreven, zelfs wanneer het anisotrope probleem goed gedefinieerd zou kunnen blijven.
3. Discretisatiecomplexiteit: Standaard conforme eindige elementenmethoden voor vierde-orde problemen ($C^1$-conform) vereisen hoge polynomiale graden (bijv. Argyris-elementen vereisen $p \geq 5$ in 2D en $p \geq 9$ in 3D) en zijn berucht moeilijk te implementeren op gekromde domeinen of in drie dimensies.

Methodologie
Het artikel stelt een $C^0$-hybride interior penalty ($C^0$-HIP) discretisatie voor. Deze aanpak maakt gebruik van een $H^1$-conforme eindige elementenruimte voor de volumineffecten en discontinuous Galerkin (DG) technieken om het gebrek aan $H^2$-conformiteit te behandelen.

• Continue Analyse: De auteurs verfijnen de continue analyse door de Cordes-conditie te gebruiken om het anisotrope effect te kwantificeren. Ze stellen vast dat het probleem Fredholm is met index nul onder een verbeterde kleinheidsvoorwaarde op de anisotropieparameter $\beta$. Cruciaal is dat deze conditie onafhankelijk is van de stabiliteitsconstante van het isotrope geval, wat een breder scala aan fysiek relevante parameters mogelijk maakt.
• Discretisatie: De methode introduceert aanvullende facet-onbekenden die de normale afgeleide van de volume-onbekenden representeren. Deze facet-variabelen ontkoppelen de volume-onbekenden, wat statische condensatie mogelijk maakt om de computationele kosten te verlagen. De discrete formulering gebruikt een lifting-operator om een discrete Laplaciaan te definiëren en bevat een stabilisatieterm om sprongen in de normale afgeleide over elementinterfaces te controleren.
• Theoretisch Kader: De analyse steunt op een abstract kader voor niet-conforme discretisaties van operatoren die compacte perturbaties zijn van coërcitieve operatoren. De auteurs bewijzen een discrete Gårding-ongelijkheid, die de "coërciviteit tot compactheid"-structuur van het continue probleem reproduceert op het discrete niveau. Deze structuur, gecombineerd met asymptotische consistentie en een discrete compactheidsargument, garandeert stabiliteit en convergentie.

Belangrijkste Bijdragen

1. Verbeterde Welgedefinieerdheidsanalyse: Het artikel bewijst dat de nematicus Helmholtz–Korteweg vergelijking goed gedefinieerd is voor een aanzienlijk breder bereik van de anisotropieparameter $\beta$ dan eerder vastgesteld. Door de Cordes-conditie te gebruiken, is de afgeleide kleinheidsvoorwaarde op $\beta$ niet afhankelijk van de resonantiefrequenties van het isotrope probleem.
2. $C^0$-Hybride Interior Penalty Methode: De auteurs introduceren een $C^0$-HIP discretisatie die de implementatieproblemen van $C^1$-conforme elementen (zoals Argyris of HCT elementen) vermijdt, terwijl het minder vrijheidsgraden vereist dan volledig discontinuous Galerkin methoden.
3. Stabiliteits- en Convergentiebewijzen: De methode is bewezen stabiel voor elke polynomiale graad $p \geq 2$, onafhankelijk van de ruimtelijke dimensie, mits de anisotropie voldoende klein is. De auteurs vestigen convergentie naar de continue oplossing onder minimale regulariteitsaannames ($u \in H^2(\Omega)$) en leiden optimale convergentiegraden af onder aanvullende regulariteit.
4. Driedimensionale Simulaties: Dit werk presenteert de eerste driedimensionale numerieke simulaties van de anisotrope effecten van nematicus velden op akoestische golfvoortplanting met behulp van deze specifieke vergelijking.

Resultaten

• Theoretische Resultaten: De discrete sesquilineaire vorm voldoet aan een discrete Gårding-ongelijkheid, wat asymptotische stabiliteit garandeert. De reeks discrete oplossingen convergeert naar de continue oplossing naarmate de maasgrootte $h \to 0$.
• Numerieke Experimenten:
  • Convergentie: Numerieke tests op een eenheidsdisk bevestigen de voorspelde convergentiegraden van $O(h^{p-1})$ voor polynomiale graden $p \in \{2, 3, 4, 5\}$.
  • Stabiliteit: Experimenten verifiëren stabiliteit voor $\beta > -\alpha$ in 2D en $\beta \in (-\alpha, 3\alpha)$ in 3D, consistent met de theoretische grenzen afgeleid uit de Cordes-conditie.
  • Anisotropie: Simulaties in 2D en 3D demonstreren het anisotrope karakter van golfvoortplanting. De resultaten tonen aan dat golfsnelheid en attenuatie variëren afhankelijk van de oriëntatie van de nematicus director $\mathbf{n}$, wat leidt tot duidelijke verstrooiingspatronen vergeleken met het isotrope geval.
  • Efficiëntie: Een vergelijking van vrijheidsgraden en niet-nul elementen toont aan dat de hybridisatiestrategie de computationele kosten aanzienlijk vermindert voor hogere polynomiale graden vergeleken met standaard $C^0$-interior penalty methoden, dankzij de eliminatie van interne bubbels via statische condensatie.

Betekenis
Het artikel claimt dat de primaire betekenis ligt in het bieden van een robuust en flexibel numeriek kader voor het oplossen van de nematicus Helmholtz–Korteweg vergelijking, wat essentieel is voor het modelleren van nematoakoestiek. Door de stabiliteitsanalyse te ontkoppelen van de isotrope stabiliteitsconstante, verwijdert dit werk een kunstmatige beperking op de anisotropieparameter, wat de simulatie van sterkere nematicus effecten mogelijk maakt. Bovendien overwint de voorgestelde $C^0$-HIP methode de praktische barrières van $C^1$-conforme methoden, waardoor hoogwaardige simulaties van vierde-orde anisotrope problemen in drie dimensies en op gekromde domeinen haalbaar worden. De numerieke voorbeelden dienen als een bewijs van concept voor het vermogen van de methode om complexe anisotrope golfverschijnselen te vatten.