Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems

Technische Samenvatting: Kwadratuurfoutschattingen op niet-aansluitende roosters in een fictief domeinraamwerk voor vloeistof-structuurinteractieproblemen

Probleemstelling
Het artikel behandelt de numerieke discretisatie van vloeistof-structuurinteractieproblemen (FSI) met behulp van een fictief domeinformulering gebaseerd op verdeelde Lagrange-multiplicatoren. In dit raamwerk worden de vloeistof en de vaste stof gediscretiseerd met behulp van twee volledig onafhankelijke, niet-aansluitende roosters. De vaste stof wordt beschreven in een Lagrangiaans raamwerk op een referentiedomein $B$, dat wordt afgebeeld op de huidige positie binnen een vast Euleriaans vloeistofdomein $\Omega$. De kinematische koppeling tussen de vloeistofsnelheid en de vervorming van de vaste stof wordt afgedwongen via een Lagrange-multiplier $\lambda$ die gedefinieerd is op het referentiedomein van de vaste stof.

Een kritieke computationele uitdaging bij het assembleren van de koppelingsmatrix (interface-matrix). Deze matrix vereist de integratie van het product van vloeistofbasisfuncties (samengesteld met de mapping van de vaste stof) en basisfuncties van de vaste stof over de elementen van de vaste stof. Omdat deze functies gedefinieerd zijn op niet-aansluitende roosters, vereist exacte integratie het berekenen van de geometrische intersectie tussen de vloeistof en de gemapte vaste stof om een samengestelde kwadratuurregel te construeren. Dit proces is computationeel duur, vooral omdat de koppelingsterm bij elke tijdstap opnieuw moet worden geassembleerd. Alternatief kan een benaderende aanpak direct over de elementen van de vaste stof integreren met behulp van standaard kwadratuurregels zonder roosterintersectie, wat een kwadratuurfout introduceert. Het artikel onderzoekt de theoretische geldigheid en de foutschattingen die verbonden zijn aan deze benaderende integratiestrategie.

Methodologie
De auteurs analyseren een stationair Stokes-gebaseerd FSI-model, maar merken op dat hun resultaten representatief zijn voor meer algemene niet-lineaire gevallen. De variationele formulering behelst het vinden van de vloeistofsnelheid $u$, druk $p$, de vervormingskaart $X$ van de vaste stof, en de Lagrange-multiplier $\lambda$ die een saddle-point systeem bevredigen.

De studie richt zich op twee specifieke keuzes voor de bilineaire vorm $c(\cdot, \cdot)$ die wordt gebruikt om de koppeling te definiëren:

1. $c_0$ (Dual pairing): Gebaseerd op de $L^2(B)$ inwendige product, waarbij de multiplicatieruimte de duale ruimte van $H^1(B)$ is.
2. $c_1$ (Scalar product): Gebaseerd op het $H^1(B)$ inwendige product, waarbij de multiplicatieruimte zelf $H^1(B)$ is.

De methodologie verloopt via de volgende stappen:

• Discretisatie: Eindige-elementeruimten worden geïntroduceerd voor alle variabelen. Voor de vloeistof wordt het Bercovier–Pironneau-element ($P1$-iso-$P2/P1$) gebruikt, terwijl de mapping van de vaste stof en de multiplier gebruikmaken van continue stuksgewijze lineaire elementen.
• Benaderende Integratie: De auteurs definiëren discrete bilineaire vormen waarbij de koppelingstermen worden berekend met standaard kwadratuurregels op het rooster van de vaste stof, waarbij wordt genegeerd dat de samengestelde vloeistofbasisfuncties stuksgewijze polynomen zijn met discontinuïteiten ten opzichte van het rooster van de vaste stof.
• Abstracte Foutanalyse: Volgens het raamwerk van inf-sup condities leiden de auteurs een algemene foutschatting (Stelling 2) af die het verschil tussen de exacte discrete oplossing en de oplossing verkregen met inexacte integratie begrenst. Deze grens is afhankelijk van de kwadratuurfoutfunctionalen.
• Kwadratuurfoutschatting: De kern van de technische bijdrage betreft het bewijzen van specifieke schattingen voor de kwadratuurfout van de koppelingstermen. De auteurs maken gebruik van het feit dat de samengestelde functie $v_h(X)$ continu maar stuksgewijs lineair is op de elementen van de vaste stof, wat leidt tot specifieke regulariteitseigenschappen in fractionele Sobolev-ruimten ($H^{1+s}$). Zij maken gebruik van interpolatietheorie en inverse ongelijkheden om de fout te begrenzen.
• Stabiliteitsbewijs: De auteurs bewijzen dat het discrete probleem welgesteld blijft (voldoet aan de discrete inf-sup condities), zelfs wanneer de koppelingsterm benaderd wordt geconstrueerd, mits aan specifieke voorwaarden met betrekking tot de roostergroottes wordt voldaan.

Belangrijkste Bijdragen en Resultaten

1. Welgesteldheid met Benaderende Koppeling: Het artikel bewijst dat het discrete probleem stabiel en welgesteld is, zelfs wanneer de koppelingsmatrix wordt geassembleerd met behulp van benaderende kwadratuur, zonder dat hiervoor roosterintersectie vereist is.
2. Kwadratuurfoutschattingen:
   • Voor de $L^2$-gebaseerde koppeling ($c_0$) is de kwadratuurfout begrensd door $O(h_B^{3/2} |\log h_B|)$. De auteurs tonen aan dat de benaderende methode optimale convergentiesnelheden oplevert, mits $h_B \to 0$ voor de roostermaat van de vaste stof.
   • Voor de $H^1$-gebaseerde koppeling ($c_1$) is de foutanalyse complexer vanwege de gradiëntterm. De foutgrens bevat termen die afhankelijk zijn van de verhouding tussen de roostermaataanpassingen. Specifiek vereist optimale convergentie niet alleen dat $h_B \to 0$, maar ook dat de verhouding $h_B/h_\Omega \to 0$ (waarbij $h_\Omega$ de vloeistofroostermaat is). Indien $h_B/h_\Omega$ constant wordt gehouden, kan de methode haar optimaliteit verliezen.
3. Inverse Ongelijkheden: De auteurs leiden specifieke inverse ongelijkheden af en gebruiken deze om de normen van de Lagrange-multiplier in de foutschattingen te behandteren, met name voor het duale normgeval.
4. Numerieke Validatie: Er worden vier numerieke tests gepresenteerd om de theoretische schattingen te valideren:
   • Test 1 & 2: Vergelijken de convergentiesnelheden voor $c_0$ en $c_1$ onder verschillende roosterverfijningsstrategieën (constante ratio versus afnemende ratio). De resultaten bevestigen dat $c_0$ robuust is met benaderende integratie, terwijl $c_1$ een $h_B/h_\Omega \to 0$ vereist om nauwkeurigheid te behouden.
   • Test 3: Gebruikt een oplossing met een sprong in de snelheidgradiënt over de interface (typisch voor FSI), wat de verwachte suboptimale convergentiesnelheid van $O(h^{1/2})$ bevestigt als gevolg van oplossingssingulariteiten.
   • Test 4: Onderzoekt hogere-orde elementen ($P2/P1$), waarbij wordt getoond dat hoewel $c_0$ robuust blijft, $c_1$ nog steeds specifieke roosterverhoudingen of exacte integratie vereist om optimale snelheden te bereiken.

Betekenis en Claims
Het artikel claimt een rigoureuze theoretische rechtvaardiging te bieden voor het gebruik van benaderende integratie in fictieve domein FSI-methoden, een praktijk die vaak wordt toegepast voor computationele efficiëntie. De auteurs demonstreren dat:

• De "inexacte" aanpak de stabiliteit van het numerieke schema niet in gevaar brengt.
• De $L^2$-gebaseerde koppeling ($c_0$) bijzonder geschikt is voor benaderende integratie, omdat deze geen specifieke verhouding tussen de verfijning van de vloeistof- en vaste stofroosters vereist om nauwkeurigheid te behouden.
• De $H^1$-gebaseerde koppeling ($c_1$), hoewel deze potentieel andere stabiliteitseigenschappen biedt, striktere roostervereisten oplegt ($h_B/h_\Omega \to 0$) om te voorkomen dat kwadratuurfouten de oplossingsfout domineren.

De auteurs concluderen dat hun theoretische resultaten de experimentele observaties uit hun eerdere werk [12, 13] bevestigen, en een solide wiskundige basis bieden voor de keuze tussen exacte (met roosterintersectie) en benaderende (niet-intersecterende) assemblage technieken op basis van de specifieke koppelingsvorm en roosterconfiguratie.