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

Resumo Técnico: Esquema Alternativo WENO Baseado em Invariantes de Riemann para um Modelo de Filme Fino de Duas Camadas

Definição do Problema
O artigo aborda a simulação numérica de fluxos de filmes finos de duas camadas governados pela tensão superficial e interfacial, especificamente na presença de partículas anti-surfactantes solúveis. Embora o trabalho anterior de Barthwal & Rohde [1] tenha estabelecido um subsistema hiperbólico de primeira ordem para tais fluxos, este estava restrito a um regime específico onde os números de Marangoni em ambas as camadas eram iguais a um ($Ma_1 = Ma_2 = 1$). Essa restrição limitou a capacidade do modelo de descrever uma gama mais ampla de dinâmicas físicas. Além disso, métodos numéricos de alta ordem padrão para sistemas de leis de conservação hiperbólicas, como os esquemas Weighted Essentially Non-Oscillatory (WENO), frequentemente sofrem com oscilações espúrias próximas a descontinuidades. Mitigar essas oscilações tipicamente requer Decomposição de Características Locais (LCD), que envolve o cálculo de autovetores e a realização de multiplicações matriz-vetor em cada interface, aumentando significativamente o custo computacional. O desafio reside no desenvolvimento de um esquema de alta ordem, não oscilatório, para o modelo estendido de duas camadas multidimensional que equilibre precisão com eficiência computacional.

Metodologia
Os autores propõem um esquema de Riemann Invariant-based Local Characteristic Decomposition WENO (RI-WENO). A metodologia procede através das seguintes etapas:

1. Extensão do Modelo: Os autores estendem o modelo unidimensional de [1] para um cenário multidimensional, relaxando a condição sobre os números de Marangoni ($Ma_1 = \alpha, Ma_2 = \beta$). Eles derivam um sistema hiperbólico reduzido de primeira ordem que governa as espessuras do filme ($f, g$) e os gradientes de concentração ($b, q$) sob a premissa de capilaridade e difusividade desprezíveis.
2. Hiperbolicidade e Invariantes de Riemann: Os autores provam que o sistema estendido é hiperbólico. Crucialmente, eles demonstram que, sob a condição onde os gradientes de concentração tangenciais desaparecem (ou estão alinhados de tal forma que o sistema se reduz a um subsistema 4x4), o sistema admite um conjunto completo de invariantes de Riemann explícitos.
3. Transformação de Variáveis: Em vez de usar o conjunto completo de invariantes de Riemann que forma um sistema de coordenadas mas exige a resolução de uma equação quártica não linear complexa para inversão (computacionalmente dispendiosa), os autores selecionam um conjunto específico de invariantes de Riemann modificados. Esta escolha resulta em uma estrutura esparsa nas matrizes de autovetores transformadas ($R$ e $L$).
4. Esquema Numérico:
   • Estrutura AWENO: O esquema utiliza uma formulação Alternative WENO (AWENO) operando sobre valores pontuais de variáveis conservadas.
   • RI-LCD: A reconstrução é realizada no espaço dos invariantes de Riemann selecionados. As variáveis conservadas são transformadas para variáveis de invariante de Riemann, a interpolação WENO é aplicada e os resultados são transformados de volta para o espaço físico usando as matrizes de autovetores esparsas.
   • Integração Temporal: Um método de integração temporal de terceira ordem, Strong Stability Preserving Runge-Kutta (SSP-RK3), é utilizado para o passo de tempo.
   • Extensão 2D: O esquema é estendido para duas dimensões usando uma abordagem dimensão por dimensão, assumindo que os gradientes de concentração se alinham ao longo de uma direção diagonal ($y=x$) para manter a validade da estrutura de invariante de Riemann.

Principais Contribuições

• Modelo Estendido: A derivação de um modelo de filme fino de duas camadas multidimensional com números de Marangoni gerais ($\alpha, \beta$), indo além do caso restritivo $Ma_1=Ma_2=1$.
• Estrutura Analítica: A prova de hiperbolicidade e a derivação de um conjunto completo de invariantes de Riemann explícitos para o sistema reduzido.
• Esquema Numérico Eficiente: O desenvolvimento do esquema RI-WENO. Ao selecionar um conjunto específico de invariantes que permite matrizes de autovetores esparsas, os autores evitam os solvers iterativos custosos exigidos pela inversão completa baseada em sistema de coordenadas, mantendo simultaneamente os benefícios de supressão de oscilações da reconstrução baseada em características.
• Benchmarking: A introdução de novos casos de teste de referência tanto em uma quanto em duas dimensões para avaliar o efeito de distintos números de Marangoni na dinâmica do fluxo.

Resultos
Experimentos numéricos foram conduzidos comparando quatro esquemas: Componentwise WENO (CP-WENO), LCD-WENO padrão, o proposto RI-WENO e um esquema baseado no conjunto completo de invariantes de Riemann (CRI-WENO).

• Precisão: Todos os esquemas alcançaram a precisão de quinta ordem esperada para soluções suaves. Em problemas envolvendo descontinuidades (choques e rarefações), o CP-WENO exibiu oscilações espúrias, enquanto o LCD-WENO e o RI-WENO suprimiram-nas efetivamente.
• Eficiência: O esquema RI-WENO demonstrou economias computacionais significativas. Em problemas de Riemann 1D, o RI-WENO exigiu de 15 a 33% menos tempo de CPU comparado ao CRI-WENO e foi consistentemente mais rápido que o LCD-WENO padrão.
• Insights Físicos: O teste de "Regime de Dois Marangoni" 1D revelou que a razão entre os números de Marangoni ($k = \alpha/\beta$) dita a direção do fluxo e a estrutura da onda. Quando $k < 1$, a segunda camada domina; quando $k > 1$, a primeira camada domina, levando a choques mais nítidos e diferentes direções de propagação de ondas.
• Desempenho 2D: Os testes 2D (perturbação circular suave e frente solutal circular perturbada) demonstraram com sucesso a capacidade do esquema de capturar padrões de ondas não lineares e redistribuição de massa impulsionadas pelo forçamento de Marangoni em um ambiente multidimensional.

Significância e Alegações
O artigo afirma que o framework RI-WENO proposto oferece um equilíbrio eficaz entre precisão e eficiência computacional. Ele fornece uma escolha prática para resolver o modelo de filme fino de duas camadas ao reduzir o custo computacional associado aos métodos padrão de decomposição de características sem comprometer a capacidade de suprimir oscilações espúrias. Os autores enfatizam que, embora um conjunto completo de invariantes de Riemann exista, a complexidade algébrica de invertê-los pode anular os ganhos de precisão e eficiência; portanto, a seleção cuidadosa de um conjunto de invariantes mais esparso é essencial para o design de métodos de alta ordem eficientes. O trabalho estabelece as bases para desenvolvimentos futuros de esquemas preservadores de positividade e de estabilidade de entropia para estes modelos de fluidos complexos.