The Tempered Finite Element Method

Technische Samenvatting: De Tempered Finite Element Method (TFEM)

Probleemstelling
De Eindige Elementen Methode (FEM) is de dominante numerieke techniek in de computationele fysica en techniek, beroemd om haar vermogen om complexe geometrieën aan te kunnen via ongestructureerde meshes. Echter, standaard FEM-formuleringen vertonen brosheid wanneer ze worden geconfronteerd met specifieke mesh-configuraties, met name die met "caps" (driehoeken met een hoek die de 180° nadert) die in banden zijn gerangschikt. Terwijl geïsoleerde degeneratieve elementen (naalden of geïsoleerde caps) de convergentie niet verhinderen, veroorzaken banden van caps een "locking"-fenomeen. In dergelijke configuraties wordt de numerieke oplossing gedwongen lineair te zijn langs de band, ongeacht de werkelijke oplossing, wat ertoe leidt dat de effectieve mesh-grootte wordt bepaald door de lengte van de band in plaats van de dikte ervan. Dit resulteert in een verlies van convergentie voor standaard FEM. Historisch gezien werd de "maximale hoekconditie" (die vereist dat alle hoeken ver verwijderd blijven van 180°) als essentieel beschouwd voor convergentie, hoewel recent werk heeft aangetoond dat dit niet strikt noodzakelijk is voor geïsoleerde elementen. De uitdaging blijft het afhandelen van meshes met elementen van nulmaat of bijna degeneratieve elementen, die relevant zijn voor interface-tracking methoden (bijv. X-Mesh) waarbij interfaces worden gevolgd via mesh-relaying, wat potentieel banden van caps kan creëren.

Methodologie
De auteurs stellen de Tempered Finite Element Method (TFEM) voor, een modificatie van de klassieke FEM die convergentie mogelijk maakt op meshes met banden van caps. Het kernidee is geïnspireerd door metallurgische ontlaten (tempering): net zoals ontlaten interne spanningen in gehard staal vermindert om de flexibiliteit te vergroten, "verzacht" TFEM de materiaaleigenschappen lokaal in de degeneratieve band om de numerieke oplossing toe te staan te conformeren aan het continue probleem.

Wiskundig gezien bevat de standaard stijfheidsmatrix $K_{ij}$ voor elliptische problemen zoals de Poisson-vergelijking de inverse van de Jacobiaan-determinant $J$ van de elementafbeelding. In degeneratieve elementen gaat $J \to 0$, wat leidt tot deling door nul of oneindige eigenwaarden. TFEM modificeert de lokale stijfheidsmatrix door de Jacobiaan weg van nul te begrenzen met behulp van een parameter $J_{min}$:
$$\tilde{K}_{ij} = \int_{\hat{E}} \frac{1}{\max(J, J_{min})} [\dots] d\xi$$
Dit is equivalent aan het introduceren van een variabele diffusiecoëfficiënt $d(x)$ in de zwakke formulering:
$$d(x) = \begin{cases} 1 & x \in \Omega \setminus B \\ D = \frac{J}{J_{min}} & x \in B \end{cases}$$
waarbij $B$ de band van degeneratieve elementen is. De parameter $J_{min}$ wordt gekozen om afhankelijk te zijn van de mesh-grootte $h$ (specifiek $J_{min} \sim h^{d+1}$ waarbij $d$ de ruimtelijke dimensie is) om een balans te vinden tussen de reductie van de diffusiecoëfficiënt (om locking te voorkomen) en het behoud van benaderingseigenschappen.

Voor het geval van exact nulmaat-elementen ($J=0$), demonstreren de auteurs dat de limiet van het TFEM-schema overeenkomt met een mortaring-methode met penalistatie van de oplossing-sprong over de interface.

Belangrijkste Bijdragen

1. Theoretisch Bewijs van Convergentie: Het artikel levert een rigoureuze foutanalyse voor TFEM op banden van caps. Het bewijst dat met een optimale keuze van de tempering-parameter $D \sim \bar{h}/h^2$ (waarbij $\bar{h}$ de bandbreedte is), de methode een optimale $O(h)$ convergentie bereikt in de $H^1$ semi-norm, wat overeenkomt met de snelheid van standaard FEM op reguliere meshes.
2. Equivalentie aan Mortaring: Voor banden met exact nulbreedte bewijzen de auteurs dat TFEM equivalent is aan een non-conforming mortar-methode met sprong-penalistatie, wat een theoretische brug slaat tussen degeneratieve mesh-afhandeling en interface-koppelingstechnieken.
3. Extensie naar Complexe Fysica: De methode wordt succesvol uitgebreid voorbij de scalaire Laplaciaan naar:
   • Lineaire Elasticiteit: Het demonstreren van vergelijkbare nauwkeurigheid met standaard FEM voor spanningsanalyse over interfaces.
   • Non-conforming Mesh Mortaring: Het aantonen dat TFEM naadloos meshes met verschillende resoluties kan verbinden zonder transitie-elementen.
   • High-Order Elementen: Het valideren dat de methode optimale convergentiesnelheden behoudt voor polynomiale orden tot 4, mits $J_{min}$ proportioneel geschaald wordt met de elementorde.
   • Advectie: Het toepassen van TFEM op advectie-diffusie problemen met behulp van Streamline Upwind/Petrov-Galerkin (SUPG) stabilisatie, waarbij wordt opgemerkt dat de tempering alleen vereist is voor de stabilisatieterm.
4. Robuustheidsanalyse: Numerieke experimenten identificeren dat de methode robuust is ten opzichte van de constante $C$ in de relatie $J_{min} = C h^{d+1}$. De fout blijft nabij optimaal over een breed bereik van $C$-waarden (dat bijna een orde van grootte beslaat), en de conditionering van het lineaire systeem blijft beheersbaar, vooral met standaard preconditioners.

Resultaten

• Numerieke Experimenten: In 2D en 3D Poisson-problemen met geconstrueerde oplossingen, herstelt TFEM op meshes met banden van caps de convergentiesnelheden van standaard FEM op reguliere meshes. Specifiek levert het kiezen van $J_{min} \sim h^3$ in 2D en $J_{min} \sim h^4$ in 3D optimale $L^2$ en $H^1$ foutconvergentie op.
• Foutgevoeligheid: De $H^1$-fout is minder gevoelig voor de keuze van $J_{min}$ dan de $L^2$-fout, en vertoont een plateau over een breed bereik van parameters.
• Conditionering: Hoewel de condition number van de stijfheidsmatrix toeneemt naarmate de bandbreedte afneemt, zorgt het gebruik van incomplete Cholesky-preconditionering ervoor dat de convergentie van iteratieve solvers onafhankelijk blijft van de banddikte.
• Gerandomiseerde Testen: Een reeks van 100 willekeurige tests in 2D bevestigde dat zelfs met een eenvoudige keuze ($C=1$), de TFEM-fouten vrijwel identiek zijn aan de standaard FEM-fouten op reguliere meshes.

Betekenis en Claims
Het artikel claimt dat TFEM de flexibiliteit van de eindige elementenanalyse aanzienlijk vergroot door de strikte eisen aan mesh-kwaliteit te versoepelen die geassocieerd worden met degeneratieve elementen. Het maakt het gebruik van meshes met zero-measure of bijna degeneratieve elementen (zoals gegenereerd door dynamische interface tracking) mogelijk zonder te lijden onder de locking-fenomenen die standaard FEM onbruikbaar maken.

De auteurs benadrukken dat de methode trivial te implementeren is in bestaande FEM-codes, aangezien het slechts een modificatie van de Jacobiaan-evaluatie vereist om deling door nul te voorkomen. Zij positioneren TFEM niet als een vervanging voor standaard FEM op hoogwaardige meshes, maar als een robuuste extensie die het gebruik van specifieke mesh-topologieën (banden van caps) mogelijk maakt die voorheen als onbruikbaar werden beschouwd. De methode is bijzonder waardevol voor toepassingen zoals de X-Mesh interface tracking methode, waarbij het handhaven van een conforme mesh moeilijk of computationeel duur is. Het artikel concludeert dat hoewel de theoretisch optimale keuze van de parameter $C$ is afgeleid, de praktische robuustheid van de methode een breed bereik aan acceptabele waarden toestaat, wat de directe toepasbaarheid ervan maakt.