Sharp-interface Simulations of Energetic Multiphase Flows with Large Density and Viscosity Ratios

Este artigo propõe um framework de Fluxos de Região de Momento de Doador Sincronizado (SynDRoM) combinado com um limitador de viscosidade para aumentar a robustez numérica e a fidelidade física de simulações de interface nítida para escoamentos multifásicos energéticos com grandes razões de densidade e viscosidade.

O essencial

A Visão Geral: Simulando um Mar Tempestuoso

Imagine que você está tentando criar uma simulação de computador de uma tempestade oceânica violenta. Você quer ver como as ondas quebram, como o ar é sugado para dentro da água e como as bolhas se formam. Isso é difícil porque a água é pesada e espessa, enquanto o ar é leve e fino. Em termos de física, eles têm uma diferença massiva de "densidade".

Quando os computadores tentam simular isso, eles costumam travar ou produzir resultados estranhos e impossíveis (como a água subitamente se transformando em um fantasma ou o ar atravessando a água como uma bala). Este artigo apresenta um novo conjunto de regras (algoritmos) para tornar essas simulações estáveis, precisas e fisicamente realistas, mesmo quando as ondas estão quebrando violentamente.

O Problema: O "Fantasma" e o "Choque"

Os autores explicam que os métodos antigos para simular esses fluxos possuem duas falhas principais:

1. O Problema do "Fantasma" (Penetração de Velocidade):
   Imagine um caminhão pesado (água) e uma pena (ar) movendo-se lado a lado. Nas simulações antigas, o "vento" da pena às vezes soprava o caminhão para trás, ou o caminhão empurrava a pena através de seu próprio corpo. Isso é chamado de "penetração de velocidade". Isso cria formas não físicas falsas na água, como um "chifre do diabo" saindo de uma onda.
2. O Problema do "Choque" (Picos de Momento):
   Para corrigir o problema do fantasma, os cientistas tentaram um novo método chamado CMOM (Momento-Massa Consistente). É como manter um registro rigoroso de quanto "impulso" (momento) cada gota de água possui. No entanto, este método tem um efeito colateral. Quando um pouco de água pesada se move para uma célula cheia de ar, a matemática fica confusa. É como dividir um número enorme por um número minúsculo, resultando em um pico de velocidade massivo e impossível. Isso cria "manchas de velocidade" — bolsões falsos de ar movendo-se a velocidades supersônicas que não deveriam existir.

A Solução: O Método "SynDRoM"

Os autores propõem uma nova correção chamada SynDRoM (Fluxo de Momento de Região de Doação Sincronizada). Veja como funciona, usando uma analogia:

A Analogia: A Esteira Rolante em Movimento
Imagine uma esteira rolante carregando caixas.

• O Jeito Antigo: Você conta as caixas (massa) e o peso das caixas (momento) separadamente. Se uma caixa se move, você pode acidentalmente contar o peso dela em um lugar onde a caixa ainda não chegou de fato. Isso causa o "choque" ou o "pico" de velocidade.
• O Jeito SynDRoM: Este método age como uma equipe sincronizada. Antes de mover o peso, você olha exatamente para qual parte da esteira rolante o peso está vindo.
  • Ele pergunta: "Se eu estou movendo este pedaço específico de ar, exatamente qual pedaço de momento está anexado a ele?"
  • Ele garante que o momento só seja movido se a massa estiver realmente lá para carregá-lo.
  • O Resultado: Sem mais picos de velocidade falsos. O ar permanece lento e a água permanece pesada, exatamente como na vida real. A simulação permanece suave e não "explode".

O Segundo Problema: A Viscosidade "Escorregadia"

O artigo também aborda um segundo problema: a Viscosidade (o quão espesso ou pegajoso um fluido é).

• O Problema: A água é pegajosa; o ar é escorregadio. Quando eles se misturam em uma fronteira nítida (como uma onda quebrando), o computador tenta adivinhar a "pegajosidade" no meio. Se ele adivinhar errado, a matemática torna-se instável, como tentar equilibrar um lápis na ponta.
• A Correção: Os autores introduzem um Limitador de Viscosidade.
  • A Analogia: Imagine uma placa de limite de velocidade. Mesmo que a matemática tente calcular uma "pegajosidade" que faria o fluido se mover de forma impossamente rápida (instável), o limitador diz: "Não, você não pode ir mais rápido do que a velocidade do fluido mais fino aqui". Ele limita o cálculo para evitar que a simulação quebre, sem alterar a física real da água ou do ar.

A Prova: Funciona?

Os autores testaram suas novas regras de três maneiras:

1. A Quebra de Represa (Dam Break): Eles simularam uma parede de água colapsando.
   • Métodos antigos: A água parecia distorcida com picos falsos.
   • SynDRoM: A água quebrou naturalmente, e o ar não foi sugado para dentro da água de maneiras estranhas.
2. A Instabilidade de Kelvin-Helmholtz: Isto é quando o vento sopra sobre a água, criando ondas que rolam (como nuvens).
   • Resultado: A simulação mostrou corretamente as ondas rolando e crescendo, sem que o computador adicionasse energia falsa ou amortecesse as ondas. Provou que o método respeita as leis da física.
3. A Onda Quebrante: Eles simularam uma onda enorme e diagonal quebrando.
   • Resultado: A onda quebrou, espirrou e criou espuma exatamente como um oceano real. A energia total do sistema permaneceu equilibrada (não desapareceu magicamente nem explodiu). Mesmo quando adicionaram "pegajosidade" (viscosidade), a simulação permaneceu estável.

A Conclusão

Este artigo apresenta um novo "policial de trânsito" para simulações de computador de água e ar.

• Ele impede que o ar atravesse a água como um fantasma.
• Ele impede que a água crie picos de velocidade impossíveis.
• Ele evita que os cálculos de "pegajosidade" quebrem a matemática.

Ao sincronizar exatamente o que está se movendo com onde está se movendo, os autores criaram uma ferramenta de simulação muito mais robusta e confiável para estudar eventos oceânicos violentos, como os que engenheiros navais precisam entender para o design de navios.

Resumo técnico

Resumo Técnico: Simulações de Interface Nítida de Fluxos Multifásicos Energéticos com Grandes Razões de Densidade e Viscosidade

Definição do Problema
Simular fluxos multifásicos energéticos com razões de densidade extremas (por exemplo, sistemas ar-água com razões $\sim 10^3$ a $10^4$) continua sendo um desafio significativo na hidrodinâmica naval. Formulações tradicionais baseadas em velocidade frequentemente sofrem de divergência numérica devido ao desalinhamento entre os gradientes de pressão e as descontinuidades de densidade, levando à transferência espúria de momento e à explosão de energia. Embora as estruturas de transporte de Momento-Momento Consistente (CMOM) melhorem a robustez ao impor a conservação de momento e a conservação de energia semidiscreta, elas introduzem novas dificuldades quando combinadas com métodos de interface nítida. Especificamente, o CMOM pode gerar choques de momento severos, resultando em oscilações não físicas de velocidade e ultrapassagens (overshoots) (violando a condição de Variação Total Limitada, ou TVD), quando interfaces nítidas cruzam as fronteiras das células. Além disso, a estimativa imprópria da viscosidade próximo às interfaces, particularmente ao dividir a viscosidade dinâmica pela densidade em grades defasadas (staggered grids), pode introduzir instabilidades numéricas em passos de tempo finitos, comprometendo a robustez do solver.

Metodologia
O artigo propõe um conjunto de esquemas numéricos de preservação de física projetados para reconciliar a fidelidade física com a robustez numérica para fluxos de alta razão de densidade.

1. Fluxo de Momento de Região de Doador Sincronizado (SynDRoM):
   Para abordar as oscilações de velocidade inerentes ao CMOM padrão, os autores introduzem o SynDRoM. Este método impõe a monotonicidade do campo de velocidade transportado sincronizando o fluxo de momento com o fluxo de massa.

• Mecanismo: Diferente dos limitadores de inclinação (slope limiters) existentes que focam apenas na reconstrução da velocidade, o SynDRoM reconstrói um par de distribuição de densidade-velocidade dentro de uma célula para satisfazer as restrições de média de célula e os limites TVD. Ele então evolui o momento usando um conceito de "região doadora" determinado estritamente pelo fluxo de massa. A velocidade advectada é calculada como uma média ponderada sobre esta região doadora específica, garantindo que a velocidade atualizada permaneça limitada.
• Implementação: O método opera dentro de uma estrutura de divisão direcional consistente com a advecção de Volume de Fluido conservativa (cVOF). Ele utiliza um esquema de integração temporal de segundo ordem de Runge-Kutta (RK2) ponderado pela densidade para lidar com a evolução não linear da velocidade em relação à densidade e ao momento.

2. Limitador de Viscosidade Cinemática Limitado:
   Para abordar as instabilidades decorrentes de saltos de viscosidade nas interfaces, o artigo introduz um limitador de viscosidade.

• Mecanismo: Os autores identificam que a interpolação padrão (aritmética ou harmônica) da viscosidade dinâmica nos vértices da grade defasada, seguida pela divisão pela densidade, pode gerar viscosidade cinemática efetiva ilimitada, violando as restrições de CFL viscoso. O limitador proposto restringe a viscosidade dinâmica efetiva usada nos termos de diferença cruzada pela densidade mínima das células de momento circundantes e por um limite superior baseado na ocupação de fase para a viscosidade cinemática. Isso garante que a viscosidade cinemática efetiva permaneça dentro de uma faixa de preservação de estabilidade.

Principais Contribuições

• Algoritmo SynDRoM: Um método generalizado para CMOM que elimina oscilações espúrias de velocidade sem sacrificar a conservação de momento ou a nitidez da interface. Ele contabiliza explicitamente a distribuição de momento e a correspondente região doadora baseada no fluxo de massa, evitando a necessidade de reconstrução total da interface em grades defasadas.
• Estabilização de Viscosidade: Um limitador de viscosidade simples, porém eficaz, que evita a divergência numérica em fluxos de alta razão de viscosidade ao limitar a viscosidade cinemática efetiva, preservando assim a estabilidade sem comprometer a fidelidade física.
• Consistência de Divisão Direcional: Uma formulação que garante que as atualizações de momento permaneçam consistentes com a advecção de massa por divisão direcional, evitando fluxos fictícios que poderiam violar a limitação da velocidade.

Resultos
Os métodos propostos foram validados através de uma série de casos de teste:

• Quebra de Barragem 2D (2D Dam Break): A formulação SynDRoM resolveu com sucesso a deformação da interface e eliminou estruturas de velocidade espúrias (como as distorções "chifre do diabo" vistas em formulações de velocidade e "blobs de velocidade" no CMOM original) enquanto mantinha a limitação de energia.
• Advecção de Gotícula 1D: Em um teste quase unidimensional com uma alta razão de densidade ($10^{-3}$), o SynDRoM manteve a monotonicidade da velocidade e as propriedades TVD, enquanto o CMOM original e variantes de média de Favre produziram ultrapassagens (overshoots) não físicas.
• Instabilidade de Kelvin-Helmholtz: O esquema capturou corretamente o crescimento linear das perturbações e a subsequente evolução não linear (pareamento de vórtices, enrolamento) sem amortecimento artificial ou explosão de energia. O momento horizontal foi conservado dentro da tolerância do solver, e as regiões de alta velocidade foram corretamente confinadas à fase de fluido leve.
• Fluxo de Poiseuille Estratificado: O limitador de viscosidade evitou com sucesso a divergência numérica observada na média aritmética padrão e corrigiu os gradientes excessivos de velocidade vistos na média harmônica, produzindo soluções estáveis e fisicamente precisas.
• Quebra de Onda Estacionária 3D: Simulações de quebras de ondas estacionárias diagonais demonstraram a capacidade do solver para lidar com extrema não linearidade, formação de jatos e entranhamento de ar. O método manteve a energia mecânica limitada mesmo na ausência de viscosidade explícita, com a dissipação de energia ocorrendo naturalmente através do mecanismo iLES (Simulação de Grandes Eddies Implícita) conforme o fluxo se tornava irregular. Estudos de convergência de malha confirmaram tendências consistentes na evolução da energia e nas taxas de dissipação.

Significância
O artigo afirma que as modificações propostas nos esquemas de advecção e viscosidade aumentam significativamente a fidelidade e a robustez das simulações de fluxo multifásico, particularmente para fluxos de superfície livre energéticos encontrados na hidrodinâmica naval. Ao garantir que os métodos CMOM satisfaçam a condição TVD próximo a células de interface de esvaziamento e ao limitar a viscosidade cinemática, a abordagem remove artefatos numéricos que poderiam desencadear eventos de quebra falsos. Os autores enfatem que o SynDRoM e o limitador de viscosidade são ortogonais ao método específico de captura de interface utilizado (demonstrado aqui com VOF geométrico), tornando a abordagem amplamente aplicável. O trabalho fornece um caminho para simular fluxos complexos de alta razão de densidade com consistência física, conservação de momento e limitação empírica de energia, abordando desafios de longa data no campo.