Parameter-robust well-posedness and discretisation for coupled Darcy--Forchheimer and advection-diffusion-reaction equations

Technische Samenvatting: Parameterrubuste welgebondenheid en discretisatie voor gekoppelde Darcy–Forchheimer en advectie-diffusie-reactie vergelijkingen

Probleemstelling
Dit werk onderzoekt de wiskundige analyse en numerieke discretisatie van een volledig gekoppeld systeem bestaande uit de Darcy–Forchheimer-vergelijkingen en de advectie-diffusie-reactie (ADR)-vergelijking. Het systeem modelleert nietlineaire stroming in poreuze media gekoppeld aan deeltjestransport, waarbij de stroomsnelheid $u$ en de druk $p$ voldoen aan een nietlineaire momentumbalans met een Forchheimer-weerstandsterm, en de concentratie van de deeltjes $\phi$ en de totale flux $\zeta$ voldoen aan een transportvergelijking. Een cruciale uitdaging die wordt aangepakt, is de afhankelijkheid van stabiliteit en convergentie van modelparameters, specifiek de hydraulische conductiviteit (permeabiliteit) $\kappa$, de Forchheimer-coëfficiënt $F$, de diffusiecoëfficiënt $\alpha$ en de reactiecoëfficiënt $\varrho_0$. De auteurs beogen welgebondenheid vast te stellen en foutschattingen af te leiden die robuust (onafhankelijk) zijn ten opzichte van deze parameters, die in praktische toepassingen zoals reservoirmodellering en optimalisatie van geothermische energie meerdere grootteordes kunnen variëren.

Methodologie
Het artikel maakt gebruik van een functioneel-analytische benadering binnen de setting van reflexieve Banachruimten, waarbij de standaard Hilbertruimte-framework wordt overstegen om de nietlineariteiten en de specifieke schaalvereisten van het probleem op te vangen.

1. Abstracte Theorie: De auteurs passen recente theorie voor geperturbeerde saddle-point problemen in Banachruimten aan. Zij breiden de Brezzi–Braess-theorie uit om parameteronafhankelijke stabiliteitsgrenzen vast te stellen. Dit omvat het definiëren van gepaste gewogen Banachruimten voor de snelheid, druk, flux en concentratie. De normen in deze ruimten worden geconstrueerd als sommen en intersecties van ruimten die gewogen zijn door de modelparameters (bijv. $\kappa^{-1/2}$, $F^{1/3}$, $\alpha^{-1/2}$).
2. Welgebondenheidsanalyse: Met behulp van een vaste-punt-argument wordt het gekoppelde systeem ontkoppeld in een Darcy–Forchheimer subprobleem en een ADR subprobleem. De auteurs bewijzen dat de individuele subproblemen welgebonden zijn en demonstreren vervolgens dat de compositie van hun oplossingsoperatoren een contracterende afbeelding vormt op een gesloten bal in de snelheidsruimte, mits de data voldoende klein zijn. Dit vestigt de existentie en uniciteit van de continue oplossing met parameteronafhankelijke stabiliteitsconstanten.
3. Eindige Elementen Discretisatie: Het probleem wordt gediscretiseerd met behulp van gemengde eindige elementenmethoden. De auteurs maken gebruik van Raviart–Thomas-elementen voor de snelheid- en fluxruimten en discontinue Galerkin (gebroken Sobolev) ruimten voor druk en concentratie. Cruciaal is dat de discrete ruimten zijn uitgerust met dezelfde parameter-gewogen normen die in de continue analyse worden gebruikt.
4. Foutanalyse: A priori foutschattingen worden afgeleid in deze gewogen normen. De analyse steunt op discrete inf-sup condities en een discreet vaste-punt-argument om aan te tonen dat de convergentiesnelheden optimaal zijn en robuust ten opzichte van $\kappa$, $F$, $\alpha$ en $\varrho_0$.
5. Preconditionering: Om de computationele uitdagingen bij het oplossen van de resulterende nietlineaire en niet-symmetrische lineaire systemen aan te pakken, ontwerpen de auteurs blok-diagonaal operator-preconditioners. Deze zijn geïnspireerd door Schur-complementstructuren en operator-preconditioneringsframeworks, toegesneden op de specifieke parametergewichten van het gekoppelde systeem.

Belangrijkste Bijdragen

• Uitbreiding naar Banachruimten: Het artikel breidt de Brezzi–Braess-theorie voor geperturbeerde saddle-point problemen vanuit Hilbertruimten uit naar reflexieve Banachruimten, waarbij specifiek aandacht wordt besteed aan de parameter-robuuste context.
• Parameterrubuste Stabiliteit: De auteurs bewijzen dat het gekoppelde Darcy–Forchheimer/ADR-systeem welgebonden en stabiel is met grenzen die onafhankelijk zijn van de permeabiliteit, Forchheimer, diffusie en reactiecoëfficiënten. Dit wordt bereikt door de zorgvuldige constructie van gewogen Banachruimten en normen.
• Robuste Discrete Analyse: Het artikel biedt een rigoureuze discrete analyse die aantoont dat de gemengde eindelijke elementenformulering deze robuuste stabiliteitseigenschappen behoudt. De afgeleide a priori foutschattingen zijn quasi-optimaal en onafhankelijk van de modelparameters en de maasgrootte $h$.
• Operator Preconditionering: Het ontwerp van variabele blok-diagonaal preconditioners wordt voorgesteld. Deze preconditioners blijken robuust te zijn ten opzichte van de modelparameters (met uitzondering van de diffusiecoëfficiënt in het advectie-gedomineerde regime) en de discretisatieparameter.

Numerieke Resultaten
De theoretische bevindingen worden gevalideerd door middel van een reeks numerieke experimenten in 2D en 3D:

• Convergentie: Tests met vervaardigde gladde oplossingen en heterogene permeabiliteit bevestigen de optimale convergentiesnelheden ($O(h^{k+1})$) in zowel standaard als parameter-gewogen normen. De resultaten laten zien dat terwijl standaard normen gevoelig kunnen zijn voor extreme parametere waarden, de gewogen normen het optimale convergentiegedrag herstellen.
• Preconditioner Prestaties: Numerieke tests op het gelineraliseerde systeem tonen aan dat de voorgestelde preconditioners conditienummers en MINRES-iteratieaantallen opleveren die stabiel blijven bij brede variaties in $\kappa$, $F$, $\alpha$ en $\varrho_0$. De auteurs merken op dat de robuustheid wordt waargenomen in het diffusie-gedomineerde regime; de prestaties in het advectie-gedomineerde regime (hoge Péclet-getallen) zijn niet volledig robuust, wat consistent is met de beperkingen van de gekozen stabilisatie.
• Praktische Toepassingen: De methode wordt toegepast op een lid-driven cavity flow probleem en een transiënt mengbaar verdringingsprobleem (viskeus vingervorming in een Hele–Shaw cel), wat het vermogen van de methode demonstreert om complexe flow-transport koppelingen en nietlineariteiten aan te kunnen.

Betekenis en Claims
De auteurs claimen dat de primaire betekenis van het werk ligt in het bieden van een wiskundig rigoureus framework voor de analyse en simulatie van gekoppelde stroming en transport in poreuze media, dat stabiel en accuraat blijft ongeacht de grootte van de fysieke parameters. Door de analyse te verschuiven naar Banachruimten en gebruik te maken van parameter-gewogen normen, overwinnen de auteurs de beperkingen van standaard Hilbertruimte-benaderingen die vaak niet in staat zijn om uniforme stabiliteit te bieden voor dergelijke nietlineaire, gekoppelde systemen. Het werk biedt een theoretische basis voor het ontwikkelen van efficiënte en robuuste numerieke solvers voor problemen waarbij parameters zoals permeabiliteit en reactiesnelheden aanzienlijk variëren, een veelvoorkomend scenario in de ondergrondse engineering en milieumodellering. De auteurs erkennen bescheiden dat hoewel het vaste-punt-framework robuuste theoretische grenzen oplevert, de waargenomen numerieke robuustheid in het gekoppelde systeem suggereert dat het analytische verlies van robuustheid in bepaalde gewogen normen een artefact kan zijn van de schaling of het vaste-punt-framework zelf, in plaats van een fundamentele beperking van de methode.