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

Sintesi Tecnica: Schema Alternativo WENO basato su Invarianti di Riemann per un Modello di Film Sottile a Due Strati

Definizione del Problema
Il documento affronta la simulazione numerica di flussi di film sottili a due strati governati dalla tensione superficiale e interfacciale, specificamente in presenza di particelle anti-surfattanti solubili. Mentre il lavoro precedente di Barthwal & Rohde [1] aveva stabilito un sottosistema iperbolico del primo ordine per tali flussi, esso era limitato a un regime specifico in cui i numeri di Marangoni in entrambi gli strati erano uguali all'unità ($Ma_1 = Ma_2 = 1$). Questa restrizione limitava la capacità del modello di descrivere una gamma più ampia di dinamiche fisiche. Inoltre, gli standard metodi numerici di alto ordine per sistemi di leggi di conservazione iperboliche, come gli schemi Weighted Essentially Non-Oscillatory (WENO), soffrono spesso di oscillazioni spurie in prossimità delle discontinuità. Mitigare tali oscillazioni richiede tipicamente la Decomposizione Caratteristica Locale (LCD), che comporta il calcolo di autovettori e l'esecuzione di moltiplicazioni matrice-vettore ad ogni interfaccia, aumentando significativamente il costo computazionale. La sfida consiste nello sviluppare uno schema ad alto ordine, non oscillatorio, per il modello esteso a più dimensioni del film sottile a due strati che bilanci accuratezza ed efficienza computazionale.

Metodologia
Gli autori propongono uno schema Riemann Invariant-based Local Characteristic Decomposition WENO (RI-WENO). La metodologia procede attraverso i seguenti passaggi:

1. Estensione del Modello: Gli autori estendono il modello monodimensionale da [1] a un contesto multidimensionale rilassando la condizione sui numeri di Marangoni ($Ma_1 = \alpha, Ma_2 = \beta$). Derivano un sistema iperbolico ridotto del primo ordine che governa gli spessori del film ($f, g$) e i gradienti di concentrazione ($b, q$) sotto l'assunzione di capillarità e diffusività trascurabili.
2. Iperbolicità e Invarianti di Riemann: Gli autori dimostrano che il sistema esteso è iperbolico. Fondamentalmente, dimostrano che sotto la condizione in cui i gradienti tangenziali di concentrazione sono nulli (o sono allineati in modo tale che il sistema si riduca a un sottosistema 4x4), il sistema ammette un insieme completo di invarianti di Riemann espliciti.
3. Trasformazione delle Variabili: Invece di utilizzare l'intero set di invarianti di Riemann che forma un sistema di coordinate ma richiede la risoluzione di un'equazione quartica non lineare complessa per l'inversione (computazionalmente costosa), gli autori selezionano un set specifico di invarianti di Riemann modificati. Questa scelta produce una struttura sparsa nelle matrici di autovettori trasformate ($R$ e $L$).
4. Schema Numerico:
   • Framework AWENO: Lo schema utilizza una formulazione Alternative WENO (AWENO) operante sui valori puntuali delle variabili conservate.
   • RI-LCD: La ricostruzione viene eseguita nello spazio dei selezionati invarianti di Riemann. Le variabili conservate vengono trasformate in variabili di invariante di Riemann, viene applicata l'interpolazione WENO e i risultati vengono trasformati nuovamente nello spazio fisico utilizzando le matrici di autovettori sparse.
   • Integrazione Temporale: Per il passaggio temporale viene utilizzato un metodo Runge-Kutta di terzo ordine a stabilità forte preservata (SSP-RK3).
   • Estensione 2D: Lo schema viene esteso a due dimensioni utilizzando un approccio dimensione per dimensione, assumendo che i gradienti di concentrazione siano allineati lungo una diagonale ($y=x$) per mantenere la validità della struttura degli invarianti di Riemann.

Contributi Chiave

• Modello Esteso: La derivazione di un modello di film sottile a due strati multidimensionale con numeri di Marangoni generali ($\alpha, \beta$), andando oltre il caso restrittivo $Ma_1=Ma_2=1$.
• Struttura Analitica: La prova dell'iperbolicità e la derivazione di un insieme completo di invarianti di Riemann espliciti per il sistema ridotto.
• Schema Numerico Efficiente: Lo sviluppo dello schema RI-WENO. Selezionando un set specifico di invarianti che consente matrici di autovettori sparse, gli autori evitano i costosi solver iterativi richiesti dall'inversione completa del sistema di coordinate, mantenendo al contempo i benefici della ricostruzione basata sulle caratteristiche per la soppressione delle oscillazioni.
• Benchmarking: L'introduzione di nuovi casi di test benchmark sia in una che in due dimensioni per valutare l'effetto di distinti numeri di Marangoni sulla dinamica del flusso.

Risultati
Sono stati condotti esperimenti numerici confrontando quattro schemi: Componentwise WENO (CP-WENO), standard LCD-WENO, il proposto RI-WENO e uno schema basato sul set completo di invarianti di Riemann (CRI-WENO).

• Accuratezza: Tutti gli schemi hanno raggiunto l'attesa accuratezza del quinto ordine per soluzioni lisce. Nei problemi che coinvolgono discontinuità (shock e rarefazioni), il CP-WENO ha mostrato oscillazioni spurie, mentre sia l'LCD-WENO che l'RI-WENO le hanno efficacementamente soppresse.
• Efficienza: Lo schema RI-WENO ha dimostrato risparmi computazionali significativi. Nei problemi di Riemann 1D, l'RI-WENO ha richiesto il 15–33% in meno di tempo CPU rispetto al CRI-WENO ed è stato costantemente più veloce del classico LCD-WENO.
• Approfondimenti Fisici: Il test del "regime a due Marangoni" 1D ha rivelato che il rapporto tra i numeri di Marangoni ($k = \alpha/\beta$) determina la direzione del flusso e la struttura delle onde. Quando $k < 1$, il secondo strato domina; quando $k > 1$, il primo strato domina, portando a shock più netti e diverse direzioni di propagazione delle onde.
• Performance 2D: I test 2D (perturbazione circolare liscia e fronte solutale circolare perturbato) hanno dimostrato con successo la capacità dello schema di catturare pattern d'onda non lineari e la ridistribuzione della massa guidata dalla forza di Marangoni in un contesto multidimensionale.

Significatività e Rivendicazioni
Il documento afferma che il framework RI-WENO proposto offre un equilibrio efficace tra accuratezza ed efficienza computazionale. Fornisce una scelta pratica per risolvere il modello di film sottile a due strati riducendo il costo computazionale associato ai metodi standard di decomposizione caratteristica senza compromettere la capacità di sopprimere le oscillazioni spurie. Gli autori sottolineano che, sebbene esista un set completo di invarianti di Riemann, la complessità algebrica dell'inversione può annullare i guadagni di accuratezza ed efficienza; pertanto, la selezione accurata di un set di invarianti più rado è essenziale per progettare metodi ad alto ordine efficienti. Il lavoro pone le basi per sviluppi futuri in schemi che preservino la positività e la stabilità entropica per questi complessi modelli fluidi.