A consistent-splitting generalized scalar auxiliary variable scheme for the perturbed Boussinesq system

Technische Samenvatting: Een Consistent-Splitting Generalized Scalar Auxiliary Variable Schema voor het Geperturbeerde Boussinesq-systeem

Probleemstelling
Dit werk behandelt de numerieke simulatie van het tweedimensionale geperturbeerde Boussinesq-systeem, dat drijvende stromingen in stabiel gestratificeerde omgevingen modelleert (bijv. atmosferische en oceanische dynamica). Het systeem wordt afgeleid door een stabiel, lineair gestratificeerd hydrostatisch evenwicht af te trekken van de standaard Boussinesq-vergelijkingen. Het resulterende geperturbeerde systeem (Vergelijking 6) vertoont een specifieke koppeling waarbij drijfvermogen de verticale beweging in de impulsvergelijking aanstuurt, terwijl verticale beweging de lineaire achtergrondtemperatuurgradiënt in de temperatuurvergelijking advecteert. Deze koppeling genereert interne zwaartekrachtgolven en drijft de relaxatie van verstoringen naar hydrostatisch evenwicht aan. De primaire wiskundige uitdaging ligt in de onsamendrukbaarheidsrestrictie en de saddle-point structuur van het snelheid-druk-systeem, verergerd door de noodzaak van hoge-orde nauwkeurigheid en stabiliteit in aanwezigheid van nietlineaire convectie en lineaire stratificatietermen.

Methodologie
De auteurs stellen een tweede-orde consistent-splitting schema voor op basis van de Generalized Scalar Auxiliary Variable (GSAV) benadering. De methodologie integreert verschillende kerncomponenten:

1. GSAV-raamwerk: Een scalaire hulpvariabele $r(t)$ wordt geïntroduceerd om de continue energiewet van het systeem na te bootsen. Dit maakt de constructie van een gemodificeerde discrete energie mogelijk die onvoorwaardelijk stabiel blijft.
2. Consistent-Splitting Tijddiscretisatie: Het schema breidt het tijdverschoven tweede-orde BDF2 (Backward Differentiation Formula) raamwerk uit, ontwikkeld door Huang en Shen [17] voor Navier–Stokes vergelijkingen.
   • Expliciete Behandeling: Nietlineaire convectie ($u \cdot \nabla u$), advectie ($u \cdot \nabla \theta$) en lineaire koppelingen (drijfvermogen $\theta e_2$ en stratificatie $\alpha v$) worden expliciet behandeld met behulp van tweede-orde extrapolatie op het verschoven tijdniveau $t_{n+k}$.
   • Ontkoppeling: Deze expliciete behandeling reduceert elke tijdstap tot het oplossen van een klein aantal ontkoppelde lineaire systemen met constante coëfficiënten, waardoor de noodzaak voor nietlineaire solvers of monolithische gekoppelde systemen wordt vermeden.
3. Drukcorrectie: Een druk-Poisson-vergelijking wordt opgelost met behulp van een consistent-splitting benadering die een optimale tweede-orde nauwkeurigheid behoudt, waarbij de inconsistenties in randvoorwaarden die vaak bij standaard projectiemethoden voorkomen, worden vermeden.
4. Stabiliteit en Foutanalyse: De auteurs bewijzen een onvoorwaardelijk zwak stabiliteitstheorema en leiden optimale tweede-orde foutschattingen af voor snelheid, druk en temperatuur. De analyse volgt rigoureus de afhankelijkheid van foutconstanten van de kinematische viscositeit ($\nu$) en thermische diffusiviteit ($\gamma$).

Belangrijkste Bijdragen
Het artikel levert drie duidelijke bijdragen vergeleken met bestaande literatuur over SAV-type schema's voor Boussinesq-systemen:

1. Analyse van het Geperturbeerde Systeem: In tegen tegenstelling tot eerdere werken die het standaard Boussinesq-systeem analyseren, richt dit artikel zich specifiek op het geperturbeerde systeem (6) dat de lineaire stratificatiekoppelingsterm $\alpha v$ bevat in de temperatuurvergelijking.
2. Integratie van Raamwerken: Terwijl Wagner et al. [36] de GSAV-benadering toepasten op het standaard Boussinesq-systeem, neemt dit werk het specifieke tijdverschoven tweede-orde BDF consistente-splitting raamwerk van Huang en Shen [17] over, maar behoudt het de standaard GSAV-update-mechanisme, wat zorgt voor volledige tweede-orde nauwkeurigheid in zowel de snelheid- als de temperatuurvergelijkingen door consistente extrapolatie van alle koppelingstermen.
3. Expliciete Afhankelijkheid van de Foutconstante: Een cruciale bijdrage is het expliciet traceren van de afhankelijkheid van de foutconstante van $\nu^{-1}$ en $\gamma^{-1}$. De analyse laat zien dat de foutconstante negatieve machten van viscositeit en diffusiviteit bevat, wat leidt tot een viervoudig genest exponentieel verband met deze parameters.

Resultaten

• Theoretische Stabiliteit: Het schema is bewezen onvoorwaardelijk zwak stabiel. De sequenties van de hulpvariabele en de schalfactor blijven niet-negatief, en de numerieke oplossingen voldoen aan uniforme grenzen onafhankelijk van de tijdstapgrootte, mits de stabilisatieconstante $C$ een drempelwaarde overschrijdt die bepaald wordt door de forceringstermen.
• Convergentie: Optimale tweede-orde convergentiegraden worden vastgesteld voor de snelheid, druk en temperatuur in passende normen.
• Numerieke Experimenten:
  • Convergentietest: Een manufactured solution op een begrensd domein bevestigt de theoretische tweede-orde convergentie in zowel $L^2$- als $H^1$-normen.
  • Simulatie van Gestratificeerde Stroming: Een langdurige simulatie ($T=10^4$) van een gestratificeerde stroming met initiële Gaussische temperatuurverstoringen reproduceert de verwachte fysieke dynamica:
    • Een drijfvermogen-transiënt regime waarbij potentiële energie wordt omgezet in kinetische energie.
    • Een interne-golf oscillatie regime met een periode die consistent is met theoretische voorspellingen ($T \approx \pi\sqrt{2/\alpha}$).
    • Een exponentieel verval regime waarbij verstoringen relaxeren naar hydrostatisch evenwicht, met vervalsnelheden die consistent zijn met de analytische grenzen afgeleid in Lemma 2.

Betekenis en Beperkingen
Het artikel beweert dat het voorgestelde schema een computationeel efficiënte methode biedt voor het simuleren van gestratificeerde stromingen, waarbij elke tijdstap wordt gereduceerd tot ontkoppelde lineaire solvers terwijl de tweede-orde nauwkeurigheid behouden blijft. De auteurs zijn echter bescheiden over de robuustheid van het schema in het limiet van verdwijnende viscositeit of diffusiviteit.

De analyse toont aan dat de foutconstante groeit als een viervoudig geneste exponent van $\nu^{-1}$ en $\gamma^{-1}$. Gevolgelijk is het schema niet robuust wanneer $\nu, \gamma \to 0$. De numerieke experimenten bevestigen dat voor kleine viscositeit en thermische diffusiviteit, het schema steeds fijnere ruimtelijke en temporele resoluties vereist om stabiel te blijven, en dat het bij grovere resoluties kan exploderen. De auteurs concluderen dat hoewel het schema effectief is voor volledig gedissipeerde regimes, een volledig impliciete of gestabiliseerde behandeling van convectietermen noodzakelijk zou zijn om foutgrenzen te verkrijgen die uniform zijn in $\nu$ en $\gamma$. Deze beperking benadrukt een afruil tussen de computationele efficiëntie van de expliciete consistent-splitting benadering en de stabiliteit in lage-diffusie regimes.