Riemann invariant-based alternative WENO scheme for a two-layer thin film model

Technische Samenvatting: Riemann-invariant-gebaseerd alternatief WENO-schema voor een twee-lagen dunne-filmmodel

Probleemstelling
Het artikel behandelt de numerieke simulatie van twee-lagen dunne-filmstromingen die worden beheerst door oppervlaktespanning en grensvlakspanning, specifels in de aanwezigheid van oplosbare anti-surfactant deeltjes. Hoewel eerder werk door Barthwal & Rohde [1] een eerste-orde hyperbolisch subsysteem voor dergelijke stromingen vaststelde, was dit beperkt tot een specifiek regime waarbij de Marangoni-getallen in beide lagen gelijk waren aan één ($Ma_1 = Ma_2 = 1$). Deze beperking beperkte het vermogen van het model om een breder scala aan fysische dynamica te beschrijven. Bovendien lijden standaard hogere-orde numerieke methoden voor systemen van hyperbolische behoudswetten, zoals Weighted Essentially Non-Oscillatory (WENO) schema's, vaak aan spuriële oscillaties nabij discontinuïteiten. Het mitigeren van deze oscillaties vereist doorgaans Local Characteristic Decomposition (LCD), wat het berekenen van eigenvectoren en het uitvoeren van matrix-vector vermenigvuldigingen bij elk grensvlak inhoudt, wat de computationele kosten aanzienlijk verhoogt. De uitdaging ligt in het ontwikkelen van een hoogwaardig, niet-oscillerend schema voor het uitgebreide meerdimensionale twee-lagen model dat een balans vindt tussen nauwkeurigheid en computationele efficiëntie.

Methodologie
De auteurs stellen een Riemann Invariant-gebaseerd Local Characteristic Decomposition WENO (RI-WENO) schema voor. De methodologie verloopt via de volgende stappen:

1. Modeluitbreiding: De auteurs breiden het eendimensionale model uit van [1] naar een meerdimensionale setting door de voorwaarde op de Marangoni-getallen te versoepelen ($Ma_1 = \alpha, Ma_2 = \beta$). Ze leiden een gereduceerd eerste-orde hyperbolisch systeem af dat de filmdiktes ($f, g$) en concentratiegradiënten ($b, q$) beheerst onder de aanname van verwaarloosbare capillaire druk en diffusiviteit.
2. Hyperboliciteit en Riemann-invarianten: De auteurs bewijzen dat het uitgebreide systeem hyperbolisch is. Cruciaal is dat zij aantonen dat, onder de voorwaarde dat de tangentiële concentratiegradiënten verdwijnen (of zodanig zijn uitgelijnd dat het systeem reduceert tot een 4x4 subsysteem), het systeem een volledige set expliciete Riemann-invarianten toelaat.
3. Variabele Transformatie: In plaats van de volledige set Riemann-invarianten te gebruiken die een coördinatensysteem vormt maar het oplossen van een complexe nietlineaire quartische vergelijking voor inversie vereist (wat computationeel duur is), selecteren de auteurs een specifieke, gemodificeerde set Riemann-invarianten. Deze keuze levert een ijle (sparse) structuur op in de getransformeerde eigenvectormatrices ($R$ en $L$).
4. Numeriek Schema:
   • AWENO Framework: Het schema maakt gebruik van een Alternative WENO (AWENO) formulering die opereert op puntwaarden van geconserveerde variabelen.
   • RI-LCD: De reconstructie wordt uitgevoerd in de ruimte van de geselecteerde Riemann-invarianten. De geconserveerde variabelen worden getransformeerd naar Riemann-invariant-variabelen, WENO-interpolatie wordt toegepast, en de resultaten worden teruggetransformeerd naar de fysieke ruimte met behulp van de ijle eigenvectormatrices.
   • Tijdintegratie: Een derde-orde Strong Stability Preserving Runge-Kutta (SSP-RK3) methode wordt gebruikt voor de tijdstapberekening.
   • 2D Extensie: Het schema wordt uitgebreid naar twee dimensies met behulp van een dimension-by-dimension benadering, waarbij wordt aangenomen dat de concentratiegradiënten langs een diagonale richting ($y=x$) zijn uitgelijnd om de geldigheid van de Riemann-invariant structuur te handhaven.

Belangrijkste Bijdragen

• Uitgebreid Model: De afleiding van een meerdimensionaal twee-lagen dunne-filmmodel met algemene Marangoni-getallen ($\alpha, \beta$), waarmee voorbij wordt gegaan aan de restrictieve $Ma_1=Ma_2=1$ casus.
• Analytische Structuur: Het bewijs van de hyperboliciteit en de afleiding van een volledige set expliciete Riemann-invarianten voor het gereduceerde systeem.
• Efficiënt Numeriek Schema: De ontwikkeling van het RI-WENO schema. Door een specifieke set invarianten te selecteren die ijle eigenvectormatrices mogelijk maakt, vermijden de auteurs de kostbare iteratieve solvers die vereist zijn voor de volledige coördinatensysteem-gebaseerde inversie, terwijl ze de voordelen van oscillatieonderdrukking behouden die gepaard gaat met de characteristic-gebaseerde reconstructie.
• Benchmarking: De introductie van nieuwe benchmark-testgevallen in zowel één als twee dimensies om het effect van verschillende Marangoni-getallen op de stromingsdynamica te evalueren.

Resultaten
Numerieke experimenten werden uitgevoerd waarbij vier schema's werden vergeleken: Componentwise WENO (CP-WENO), standaard LCD-WENO, het voorgestelde RI-WENO, en een schema gebaseerd op de volledige set Riemann-invarianten (CRI-WENO).

• Nauwkeurigheid: Alle schema's bereikten de verwachte vijfde-orde nauwkeurigheid voor gladde oplossingen. In problemen met discontinuïteiten (schokken en rarefacties) vertoonde CP-WENO spuriële oscillaties, terwijl zowel LCD-WENO als RI-WENO deze effectief onderdrukten.
• Efficiëntie: Het RI-WENO schema vertoonde significante computationele besparingen. In 1D Riemann-problemen vereiste RI-WENO 15–33% minder CPU-tijd vergeleken met CRI-WENO en was het consequent sneller dan standaard LCD-WENO.
• Fysische Inzichten: De 1D "Two-Marangoni regime" test onthulde dat de ratio van de Marangoni-getallen ($k = \alpha/\beta$) de stroomrichting en golfstructuur bepaalt. Wanneer $k < 1$, domineert de tweede laag; wanneer $k > 1$, domineert de eerste laag, wat leidt tot scherpere schokken en andere golfvoortplantingsrichtingen.
• 2D Prestaties: De 2D tests (gladde circulaire perturbatie en een geperturbeerde circulaire solutale front) demonstreerden succesvol het vermogen van het schema om nietlineaire golfpatronen en massaverdeling te vangen, gedreven door Marangoni-kracht in een meerdimensionale setting.

Betekenis en Claims
Het artikel claimt dat het voorgestelde RI-WENO framework een effectieve balans biedt tussen nauwkeurigheid en computationele efficiëntie. Het biedt een praktische keuze voor het oplossen van het twee-lagen dunne-filmmodel door de computationele kosten te verlagen die geassocieerd worden met standaard characteristic decomposition methoden, zonder de capaciteit tot het onderdrukken van spuriële oscillaties in gevaar te brengen. De auteurs benadrukken dat hoewel een volledige set Riemann-invarianten bestaat, de algebraïsche complexiteit van het inverteren ervan de winst in nauwkeurigheid en efficiëntie kan tenietdoen; daarom is de zorgvuldige selectie van een ijlere set invarianten essentieel voor het ontwerpen van efficiënte hogere-orde methoden. Dit werk legt de basis voor toekomstige ontwikkelingen in positiviteits-behoudende en entropie-stabiele schema's voor deze complexe fluïdummodellen.