Free-surface Stokes problem: stability estimates and time-step improvements

Technische Samenvatting: Vrije-oppervlakte Stokes-probleem: stabiliteitsschattingen en verbeteringen in de tijdstap

Probleemstelling
Nauwkeurige simulatie van stromingen met een zeer hoge viscositeit bij een vrij oppervlak, zoals ijskapdynamiek, mantelconvectie en lavastromen, vereist het oplossen van de gekoppelde Stokes/vrije-oppervlaktevergelijkingen. Deze problemen worden gekenmerkt door een bewegend domein waarbij de vloeistofsnelheid en druk gekoppeld zijn aan de evolutie van de hoogte van het vrije oppervlak. Hoewel standaard discretisaties met eindige elementen in de ruimte en Euler-methoden in de tijd gebruikelijk zijn, lijden deze vaak onder strikte stabiliteitsbeperkingen. Specifiek vereisen expliciete en semi-impliciete Euler-schema's vaak prohibitief kleine tijdstappen om stabiel te blijven, wat leidt tot hoge computationele kosten. Bestaande stabilisatietechnieken, zoals het Free-Surface Stabilization Algorithm (FSSA), verbeteren de tijdstapgroottes, maar missen a priori theoretische grenzen voor stabiliteit en garanderen geen onvoorwaardelijke stabiliteit of volumebewaring.

Methodologie
De auteurs presenteren een rigoureuze theoretische analyse van de stabiliteits- en conserveringseigenschappen van het gekoppelde Stokes/vrije-oppervlakte systeem voor zowel Newtoniaanse als niet-Newtoniaanse (power-law viscositeit) vloeistoffen. De analyse wordt uitgevoerd op twee niveaus:

1. Continu niveau: De auteurs leiden een totale energiebalans af die de energie van de Stokes-impulsgroepen koppelt aan de energie van de vrije-oppervlaktehoogtevergelijking. Deze afleiding berust op het vaststellen van specifieke relaties tussen de vrije-oppervlakte-advectieterm en oppervlakte-integralen, en het vervolgens koppelen hiervan aan de zwaartekrachtterm in de Stokes-vergelijkingen.
2. Discreet niveau: De studie onderzoekt volledig gediscretiseerde schema's met eindige elementen in de ruimte gecombineerd met expliciete, semi-impliciete en impliciete Euler-tijdstapmethoden.
   • Regularisatie: Om de gebrekkige inherente regulariteit van de vrije-oppervlaktehoogtefunctie aan te pakken (die vereist is voor de welgedetermineerdheid van het Stokes-probleem), introduceren de auteurs een hoogwaardige edge-stabilisatieterm, $J(s^{n+1}, w)$, die impliciet in de tijd wordt behandeld.
   • Stabilisatievoorstel: Gemotiveerd door de identificatie van "spooktermen" in de energie-inschattingen van de expliciete Euler-methode, stellen de auteurs een nieuwe stabilisatieterm voor die aan de Stokes zwakke vorm wordt toegevoegd. Deze term straft de normale snelheid aan het vrije oppervlak af en is ontworpen om de destabiliserende termen die voortvloeien uit de expliciete behandeling van de vrije-oppervlaktevergelijking precies te elimineren.

Belangrijkste Bijdragen

• Totale Energiebalans: Het artikel leidt een nieuwe totale energiebalans af voor het gekoppelde Stokes/vrije-oppervlakte systeem, waarbij de momentumenergie wordt gekoppeld aan de oppervlaktehoogte-energie. Dit wordt geïdentificeerd als een fundamentele bouwsteen voor het begrijpen van stabiliteit, die voorheen ontbrak in de literatuur.
• Theoretische Stabiliteitsbewijzen:
  • Impliciete Euler: Is bewezen onvoorwaardelijk stabiel met betrekking tot de grootte van de tijdstap (ervan uitgaande dat de oppervlaktehoogte Lipschitz-continu is) en conserveert het domeinvolume.
  • Expliciete Euler (Ongestabiliseerd): Is bewezen instabiel in de totale energienorm voor elke tijdstap vanwege spooktermen gerelateerd aan de normale oppervlaktesnelheid.
  • Semi-impliciete Euler (FSSA): Is aangetoond slechts voorwaardelijk stabiel te zijn en, cruciaal, is bewezen het domeinvolume niet te conserveren.
  • Voorgestelde Gestabiliseerde Expliciete Euler: De auteurs bewijzen dat hun nieuwe stabilisatieterm de expliciete Euler-methode onvoorwaardelijk stabiel maakt met betrekking tot de grootte van de tijdstap (onder de Lipschitz-continuïteitsveronderstelling), terwijl het domeinvolume behouden blijft.
• Volumebewaring: Het artikel biedt theoretische bewijzen en numerieke verificatie dat het voorgestelde gestabiliseerde expliciete schema het domeinvolume conserveert, terwijl het FSSA-gestabiliseerde semi-impliciete schema dat niet doet.

Resultaten
Numerieke experimenten werden uitgevoerd met de open-source bibliotheek FEniCS voor zowel Newtoniaanse als niet-Newtoniaanse gevallen, inclusere een complexe geometrie gebaseerd op gegevens van de Groenlandse ijskap.

• Stabiliteit: Het ongestabiliseerde expliciete Euler-schema was instabiel voor alle geteste tijdstappen. De FSSA-gestabiliseerde schema's (zowel semi-impliciet als expliciet) waren voorwaardelijk stabiel en faalden bij grotere tijdstappen. In tegen plaats bleef het EE-Stabilized schema stabiel voor willekeurig grote tijdstappen in zowel de Newtoniaanse als de niet-Newtoniaanse gevallen.
• Volumebewaring: Het EE-Stabilized schema conserveerde het domeinvolume tot op de machineprecisie. Het SIE-FSSA-schema slaagde er niet in het volume te conserveren, wat de theoretische bevindingen bevestigt.
• Nauwkeurigheid: Convergentietests toonden aan dat de voorgestelde stabilisatietermen consistent zijn met het gediscretiseerde probleem. Het EE-Stabilized schema behield eerste-orde nauwkeurigheid in zowel de oppervlaktehoogte als de snelheid onder gelijktijdige verfijning van het rooster en de tijdstap, met fouten die vergelijkbaar met of kleiner dan andere gestabiliseerde schema's zijn.
• Oscillaties: De impliciete edge-stabilisatieterm onderdrukte hoogfrequente oscillaties in de vrije-oppervlaktehoogte effectief zonder extra restricties op de tijdstap op te leggen.

Betekenis en Claims
Het artikel claimt dat het voorgestelde gestabiliseerde expliciete Euler-schema een robuust en efficiënt alternatief biedt voor computationeel dure impliciete oplossers. Door onvoorwaardelijke stabiliteit met betrekking tot de grootte van de tijdstap theoretisch te garanderen en volumebewaring te waarborgen, maakt de methode stabiele simulaties van vrije-oppervlaktestromingen mogelijk, zelfs met complexe geometrieën en brontermen. De auteurs benadrukken dat hun stabilisatie de eerste is die willekeurige tijdstappen theoretisch garandeert voor de expliciete Euler-discretisatie van dit gekoppelde probleem. De aanpak wordt gekenmerkt als symmetrisch (gunstig voor iteratieve oplossers), eenvoudig te implementeren en vrij van instelbare parameters. De auteurs concluderen dat hoewel hun analyse uitgaat van Lipschitz-continuïteit van de vrije-oppervlaktehoogte, de afgeleide $L^2$-norm stabiliteitsschattingen een noodzakelijke basis vormen voor welgedetermineerdheid zodra passende regularisatie wordt toegepast. Toekomstig werk wordt gesuggereerd om regularisatietermen voor welgedetermineerdheid te verkennen en om de stabiliteitsanalyse uit te breiden naar Navier-Stokes/vrije-oppervlakteproblemen en hogere-orde tijdintegratoren.