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

Resumen Técnico: Esquema WENO Alternativo basado en Invariantes de Riemann para un Modelo de Película Delgada de Dos Capas

Planteamiento del Problema
El artículo aborda la simulación numérica de flujos de película delgada de dos capas gobernados por la tensión superficial e interfacial, específicamente en presencia de partículas de anti-surfactante solubles. Mientras que el trabajo previo de Barthwal & Rohde [1] estableció un subsistema hiperbólico de primer orden para tales flujos, este estaba restringido a un régimen específico donde los números de Marangoni en ambas capas eran iguales a la unidad ($Ma_1 = Ma_2 = 1$). Esta restricción limitaba la capacidad del modelo para describir una gama más amplia de dinámicas físicas. Además, los métodos numéricos de alto orden estándar para sistemas de leyes de conservación hiperbólicas, como los esquemas de Oscilación Casi Nula Ponderada (WENO, por sus siglas en inglés), suelen sufrir oscilaciones espurias cerca de las discontinuidades. Mitigar estas oscilaciones requiere típicamente una Descomposición de Características Locales (LCD), lo que implica calcular autovectores y realizar multiplicaciones matriz-vector en cada interfaz, aumentando significativamente el costo computacional. El desafío radica en desarrollar un esquema de alto orden, no oscilatorio, para el modelo extendido de dos capas multidimensional que equilibre la precisión con la eficiencia computacional.

Metodología
Los autores proponen un esquema de Descomposición de Características Locales basado en Invariantes de Riemann (RI-WENO). La metodología procede a través de los siguientes pasos:

1. Extensión del Modelo: Los autores extienden el modelo unidimensional de [1] a un entorno multidimensional relajando la condición sobre los números de Marangoni ($Ma_1 = \alpha, Ma_2 = \beta$). Derivan un sistema hiperbólico reducido de primer orden que gobierna los espesores de la película ($f, g$) y los gradientes de concentración ($b, q$) bajo la suposición de capilaridad y difusividad despreciables.
2. Hiperbolicidad e Invariantes de Riemann: Los autores demuestran que el sistema extendido es hiperbólico. Crucialmente, demuestran que bajo la condición de que los gradientes de concentración tangenciales desaparecen (o están alineados de tal manera que el sistema se reduce a un subsistema de 4x4), el sistema admite un conjunto completo de invariantes de Riemann explícitos.
3. Transformación de Variables: En lugar de utilizar el conjunto completo de invariantes de Riemann que forma un sistema de coordenadas pero requiere resolver una ecuación cuártica no lineal compleja para la inversión (computacionalmente costosa), los autores seleccionan un conjunto específico de invariantes de Riemann modificados. Esta elección produce una estructura dispersa en las matrices de autovectores transformadas ($R$ y $L$).
4. Esquema Numérico:
   • Marco AWENO: El esquema utiliza una formulación de WENO Alternativa (AWENO) que opera sobre los valores puntuales de las variables conservadas.
   • RI-LCD: La reconstrucción se realiza en el espacio de los invariantes de Riemann seleccionados. Las variables conservadas se transforman a variables de invariante de Riemann, se aplica la interpolación WENO y los resultados se transforman de vuelta al espacio físico utilizando las matrices de autovectores dispersas.
   • Integración Temporal: Se utiliza un método de Runge-Kutta de orden superior con Estabilidad de Preservación de la Fuerza (SSP-RK3) para el paso de tiempo.
   • Extensión 2D: El esquema se extiende a dos dimensiones mediante un enfoque dimensión por dimensión, asumiendo que los gradientes de concentración se alinean a lo largo de una dirección diagonal ($y=x$) para mantener la validez de la estructura de los invariantes de Riemann.

Contribuciones Clave

• Modelo Extendido: La derivación de un modelo de película delgada de dos capas multidimensional con números de Marangoni generales ($\alpha, \beta$), yendo más allá del caso restrictivo de $Ma_1=Ma_2=1$.
• Estructura Analítica: La demostración de la hiperbolicidad y la derivación de un conjunto completo de invariantes de Riemann explícitos para el sistema reducido.
• Esquema Numérico Eficiente: El desarrollo del esquema RI-WENO. Al seleccionar un conjunto específico de invariantes que permite matrices de autovectores dispersas, los autores evitan el costoso solver iterativo requerido por la inversión completa del sistema de coordenadas, manteniendo al mismo tiempo los beneficios de supresión de oscilaciones de la reconstrucción basada en características.
• Evaluación de Referencia (Benchmarking): La introducción de nuevos casos de prueba tanto en uno como en dos dimensiones para evaluar el efecto de distintos números de Marangoni en la dinámica del flujo.

Resultados
Se realizaron experimentos numéricos comparando cuatro esquemas: WENO por componentes (CP-WENO), LCD-WENO estándar, el propuesto RI-WENO, y un esquema basado en el conjunto completo de invariantes de Riemann (CRI-WENO).

• Precisión: Todos los esquemas alcanzaron la precisión de quinto orden esperada para soluciones suaves. En problemas con discontinuidades (choques y rarefacciones), el CP-WENO exhibió oscilaciones espurias, mientras que tanto el LCD-WENO como el RI-WENO las suprimieron eficazmente.
• Eficiencia: El esquema RI-WENO demostró ahorros computacionales significativos. En problemas de Riemann en 1D, el RI-WENO requirió entre un 15 y un 33% menos de tiempo de CPU en comparación con el CRI-WENO y fue consistentemente más rápido que el LCD-WENO estándar.
• Percepciones Físicas: La prueba del "régimen de dos Marangoni" en 1D reveló que la relación de los números de Marangoni ($k = \alpha/\beta$) dicta la dirección del flujo y la estructura de la onda. Cuando $k < 1$, la segunda capa domina; cuando $k > 1$, la primera capa domina, lo que conduce a choques más agudos y diferentes direcciones de propagación de ondas.
• Desempeño 2D: Las pruebas en 2D (perturbación circular suave y frente de solución circular perturbado) demostraron con éxito la capacidad del esquema para capturar patrones de ondas no lineales y la redistribución de masa impulsada por el forzamiento de Marangoni en un entorno multidimensional.

Significancia y Reivindicaciones
El artículo afirma que el marco RI-WENO propuesto ofrece un equilibrio efectivo entre precisión y eficiencia computacional. Proporciona una opción práctica para resolver el modelo de película delgada de dos capas al reducir el costo computacional asociado con los métodos estándar de descomposición de características sin comprometer la capacidad de suprimir oscilaciones espurias. Los autores enfatizan que, si bien existe un conjunto completo de invariantes de Riemann, la complejidad algebraica de invertirlos puede anular las ganancias de precisión y eficiencia; por lo tanto, la selección cuidadosa de un conjunto de invariantes más disperso es esencial para diseñar métodos de alto orden eficientes. Este trabajo sienta las bases para desarrollos futuros de esquemas de preservación de positividad y estabilidad de entropía para estos modelos de fluidos complejos.