Convergence of a CutFEM for fluid--structure interaction with a deforming interface

Technische Samenvatting: Convergentie van een CutFEM voor Fluid–Structure Interaction met een Deformerend Interface

Probleemstelling
Dit werk behandelt de numerieke benadering van Fluid–Structure Interaction (FSI) problemen waarbij een incompressibele viskeuze vloeistof en een lineair elastische vaste stof betrokken zijn. Het systeem wordt gekenmerkt door een bewegend interface waarbij het vloeistofdomein wordt gedefinieerd door de deformatie van de vaste stof. De sturende vergelijkingen bestaan uit de incompressibele Navier–Stokes vergelijkingen in het vloeistofdomein (Euleriaanse beschrijving) en lineaire elasticiteitsvergelijkingen in het vaste domein (Lagrangiaanse beschrijving). De koppeling wordt afgedwongen via kinematische continuïteit (snelheidsmatching) en dynamische balans van spanningen op het evoluerende interface.

De primaire uitdaging bij het numeriek oplossen van dergelijke problemen, met name bij grote deformaties, is de beweging van het interface ten opzichte van het rekenrooster (mesh). Traditionele fitted mesh methoden (bijv. Arbitrary Lagrangian–Eulerian of ALE) lijden onder mesh-degradatie onder grote rotaties of deformaties. Unfitted methoden, zoals Immersed Boundary of Fictitious Domain methoden, vermijden mesh-regeneratie maar hebben vaak moeite met nauwkeurigheid nabij het interface door het gebrek aan resolutie van discontinuïteiten of randvoorwaarden.

Methodologie
De auteurs stellen een semi-discrete immersed CutFEM (Cut Finite Element Method) voor en analyseren deze voor dit volledig gekoppelde FSI-systeem. De methode maakt gebruik van de volgende kerncomponenten:

1. Unfitted Mesh: De vloeistof wordt gediscretiseerd op een vast achtergrondrooster dat niet voldoet aan de grens van de vaste stof. De vaste stof wordt gediscretiseerd met een isoparametrische fitted mesh op zijn referentiedomein.
2. Lagrange Multiplier Formulering: De fluid–solid koppelingsvoorwaarden worden zwak afgedwongen met behulp van een Lagrange multiplier. Deze multiplier vertegenwoordigt de vloeistofspanning (traction) die teruggetrokken is naar het referentie-interface van de vaste stof. De kinematische conditie ($u = \partial_t d \circ X^{-1}$) wordt afgedwongen op het huidige interface, terwijl de balans van spanningen wordt geformuleerd via een integraalrelatie die betrekking heeft op de referentieconfiguratie.
3. Ghost Penalty Stabilisatie: Om stabiliteit en robuustheid te garanderen tegen willekeurig kleine cut-elementen, worden globale ghost-penalty termen toegevoegd aan de snelheid- en drukformuleringen. Deze termen bestraffen sprongen in normale afgeleiden over de interne zijden van het achtergrondrooster, waardoor de stabiliteit van de eindige elementruimte effectief wordt uitgebreid naar een uitgebreid vloeistofdomein.
4. Semi-Discrete Schema: De methode is ruimtelijk semi-discreet, wat betekent dat de tijd continu is. De vloeistofsnelheid, druk, vaste stof verplaatsing en de Lagrange multiplier worden benaderd met eindige elementruimten van graad $k$ (voor snelheid en verplaatsing) en $k-1$ (voor druk).

Belangrijkste Bijdragen en Theoretisch Kader
De centrale bijdrage van dit artikel is de eerste rigoureuze a priori foutanalyse voor een unfitted eindige element methode toegepast op een volledig gekoppeld FSI-systeem met een werkelijk deformerend interface. Eerdere rigoureuze analyses waren grotendeels beperkt tot gelineariseerde modellen, benaderingen met kleine verplaatsingen of fitted mesh benaderingen.

Het bewijs van convergentie steunt op verschillende nieuwe analytische ingrediënten:

• Geïnterpoleerde Configuratie: De auteurs introduceren een geïnterpoleerde configuratie $X^*_H$ van de exacte oplossing om het exacte interface $\Gamma[X]$, het geïnterpoleerde interface $\Gamma[X^*_H]$ en het numerieke interface $\Gamma[X_H]$ met elkaar te vergelijken. Dit maakt een systematische kwantificering mogelijk van de fouten die voortvloeien uit de verplaatsing van het interface.
• Uniforme Schattingen: De analyse vestigt uniforme normequivalenties, trace-ongelijkheden, inverse schattingen en inf–sup condities op de bewegende cut-domeinen. Cruciaal is dat deze schattingen uniform zijn met betrekking tot de positie van het interface binnen het achtergrondrooster en met betrekking tot de tijd, mits de numerieke deformatie voldoende dicht bij de exacte deformatie blijft.
• Koppeling van Fouten: Een aanzienlijke analytische moeilijkheid die wordt aangepakt, is de tweezijdige afhankelijkheid tussen de vloeistoffout en de geometriefout (vaste stof verplaatsing). De fout in de verplaatsing van de vaste stof verandert het numerieke vloeistofdomein en het interface, wat op zijn beurt de vloeistofoplossing beïnvloedt. De auteurs sluiten deze lus door stabiliteitsschattingen af te leiden voor de geprojecteerde kinematische mismatch op het interface en een continuatie-in-tijd argument te gebruiken.
• Continuatie-argument: Het bewijs gaat uit van a priori bounds op de fout (specifiek $W^{1,\infty}$ en $L^\infty$ bounds) tot een maximale tijd $t^*$. Door te bewijzen dat de energie-schattingen met constanten onafhankelijk van $t^*$ werken, en te tonen dat voor voldoende kleine mesh-groottes en $k \ge 3$, de afgeleide convergentiesnelheden de aangenomen bounds impliceren, breiden de auteurs de geldigheid van de schattingen uit naar het gehele tijdsinterval $[0, T]$.

Hoofdresultaten
Onder aannames van voldoende regulariteit van de exacte oplossing en mesh-groottes die voldoen aan $c_1 H \le h \le c_2 H$ (waarbij $h$ en $H$ respectievelijk de vloeistof- en vaste stof mesh-groottes zijn), bewijzen de auteurs de volgende convergentiesnelheden voor eindige elementen van graad $k \ge 3$:

• Vloeistofsnelheid ($u_h$) en Vaste Stof Verplaatsing ($d_H$): Convergentie van orde $k$ in de natuurlijke energie-normen.
• Vloeistofdruk ($p_h$): Convergentie van orde $k-1$ in de $L^2$ norm (geassocieerd met de energie-norm).
• Lagrange Multiplier (Interface Spanning): Convergentie van orde $k-1$ in de $H^{-1/2}$ norm.

Specifiek wordt de foutbound gegeven door:
$$\max_{t \in [0,T]} \left( \|u_h - u\|_{L^2}^2 + \|d_H - d \circ \Phi_H\|_{H^1}^2 + \|\partial_t d_H - \partial_t d \circ \Phi_H\|_{L^2}^2 \right) + \int_0^T \|u_h - u\|_{H^1}^2 dt + H^2 \int_0^T (\|p_h - p\|_{L^2}^2 + \|\sigma_H - \sigma \circ \Phi_H\|_{H^{-1/2}}^2) dt \le C(h^{2k} + H^{2k})$$

Betekenis
Het artikel beweert dat dit werk de eerste rigoureuze convergentieanalyse is van een unfitted eindige element methode voor een volledig gekoppeld FSI-systeem met een deformerend interface. Hoewel eerdere werken unfitted methoden hebben geanalyseerd voor voorgeschreven bewegende domeinen of vereenvoudigde gekoppelde problemen, was de nonlineariteit en de feedbackloop tussen de vloeistofoplossing en de bewegende geometrie in een volledig gekoppeld FSI-systeem in de unfitted context nog niet rigoureus behandeld. De resultaten bieden een theoretische basis voor het gebruik van CutFEM in complexe FSI-toepassingen met grote deformaties, waarbij de nauwkeurigheid en stabiliteit van de methode worden gevalideerd zonder te vertrouwen op kleine-verplaatsings of fixed-domain benaderingen. De restrictie $k \ge 3$ wordt genoteerd als technisch, gerelateerd aan specifieke inverse schattingen gebruikt in het bewijs, en de auteurs suggereren dat het resultaat kan worden uitgebreid naar $k=2$ met toekomstige verfijningen.