A HHO formulation for variable density incompressible flows where the density is purely advected

Este trabalho propõe uma formulação Hybrid High-Order (HHO) para escoamentos incompressíveis com densidade variável que garante conservação exata de volume e advecção pura da densidade, utilizando espaços híbridos e métodos temporais ESDIRK para oferecer estabilidade, precisão e eficiência na simulação de misturas de fluidos imiscíveis, incluindo a instabilidade de Rayleigh-Taylor.

O essencial

Imagine que você está tentando prever como duas tintas de cores diferentes (uma pesada e outra leve) se misturam quando você as despeja em um balde de água. Se a tinta pesada estiver em cima, ela vai descer em "túneis" e a leve vai subir em "bolhas", criando um caos bonito chamado Instabilidade de Rayleigh-Taylor.

O problema é que, para computadores fazerem essa previsão com precisão, eles precisam resolver equações matemáticas extremamente complexas. Se o computador errar um pouquinho na pressão ou na velocidade, a simulação pode ficar instável, criar "fantasmas" (valores negativos de densidade) ou simplesmente travar.

Este artigo apresenta uma nova ferramenta matemática chamada HHO (High-Order Hybrid) combinada com um método de tempo chamado ESDIRK. Vamos usar algumas analogias para entender como isso funciona:

1. O Problema: O "Quebra-Cabeça" da Mistura

Pense no fluido como uma multidão de pessoas em uma sala.

• Densidade Variável: Algumas pessoas são muito pesadas (chumbo), outras são leves (balões de hélio).
• Incompressibilidade: A sala não pode encolher nem crescer; o número de pessoas por metro quadrado deve ser constante, mesmo que elas se movam.
• O Desafio: Se você tentar calcular como cada pessoa se move, a matemática fica tão pesada que o computador "sufoca". Além disso, se o computador calcular mal a pressão, ele pode dizer que uma pessoa de chumbo se transformou em um balão de hélio magicamente, o que é fisicamente impossível.

2. A Solução: O Método HHO (O "Equipe de Construção")

Os autores criaram um método que divide o trabalho de forma inteligente, como uma equipe de construção de uma casa:

• A "Fachada" e o "Interior": Em vez de tentar calcular tudo de uma vez só em cada peça do quebra-cabeça (o computador), o método HHO olha para o interior de cada peça (o elemento da malha) e para a fachada (as bordas que tocam as outras peças).

  • Analogia: Imagine que você está pintando um mosaico. Em vez de pintar cada pedrinha inteira de uma vez, você foca em como as pedrinhas se encaixam nas bordas (a fachada) e deixa o interior ser preenchido de forma automática e eficiente. Isso economiza muita memória, como se você estivesse guardando apenas o desenho das bordas e não de cada pedrinha.

• O "Guardião da Pressão" (Robustez): Um dos maiores problemas em simulações de fluidos é que erros na pressão podem estragar a velocidade do fluido.

  • Analogia: Pense na pressão como o "chefe" e a velocidade como o "funcionário". Em métodos antigos, se o chefe gritasse errado, o funcionário trabalharia mal. Neste novo método, o "funcionário" (velocidade) é tão inteligente e bem treinado que ele ignora os gritos errados do chefe. Ele só obedece às forças reais (como a gravidade), mantendo a simulação estável mesmo se a pressão for calculada de forma aproximada. Isso é chamado de robustez à pressão.

• O "Transporte Puro": A densidade (a "cor" da tinta) é apenas transportada pelo fluxo, sem ser criada ou destruída.

  • Analogia: Imagine que a densidade é um pacote de encomendas. O método garante que o pacote nunca seja perdido, nunca seja duplicado e nunca seja esmagado. Ele viaja exatamente onde o vento (o fluido) o leva. O artigo garante que o volume total é conservado com precisão de máquina (quase perfeita).

3. O Motor do Tempo: ESDIRK (O "Relógio de Precisão")

Para simular o movimento ao longo do tempo, eles usam um método chamado ESDIRK.

• Analogia: Imagine que você está dirigindo um carro em uma estrada cheia de curvas (o fluido mudando de forma). Um método comum seria dar um "pulo" grande a cada segundo. Se a curva for muito fechada, você sai da pista.
• O ESDIRK é como um piloto de F1 que olha para frente, calcula várias curvas pequenas dentro de um único segundo e ajusta a direção instantaneamente. Isso permite que o computador dê passos de tempo maiores sem perder a precisão ou a estabilidade, economizando tempo de cálculo.

4. Os Resultados: O Teste da "Torre de Chumbo"

Os autores testaram sua ferramenta em dois cenários extremos:

1. Mistura Suave (Baixo Atwood): Duas tintas com densidades parecidas. O método funcionou perfeitamente, mostrando que é possível usar polinômios de alta ordem (cálculos muito detalhados) em malhas mais grossas (poucas peças) e obter o mesmo resultado que malhas super finas. É como conseguir ver uma imagem HD usando menos pixels, economizando bateria.
2. Mistura Extrema (Alto Atwood): Uma tinta de chumbo sobre uma de isopor. Aqui, o risco de erro é enorme. O método conseguiu manter a densidade sempre positiva (nada de "densidade negativa" ou fantasmas), mesmo com polinômios simples (ordem baixa), provando que é estável e seguro para misturas muito diferentes.

Resumo em uma frase

Os autores criaram um "super-cálculo" para fluidos que é tão eficiente que economiza memória do computador, tão inteligente que ignora erros de pressão e tão preciso que consegue simular misturas de fluidos pesados e leves sem "quebrar" a física, tudo isso usando uma estratégia de tempo que é como um piloto de corrida ajustando a direção milissegundo a milissegundo.

Por que isso importa?
Isso permite que cientistas e engenheiros simulem fenômenos reais (como explosões, mistura de combustíveis ou até o movimento de magma na Terra) com muito mais rapidez e precisão, sem precisar de supercomputadores gigantescos para tudo.

Resumo técnico

Aqui está um resumo técnico detalhado do artigo em português, cobrindo o problema, metodologia, contribuições, resultados e significância.

Título: Formulação HHO para Fluxos Incompressíveis de Densidade Variável com Densidade Puramente Advectada

Autores: Lorenzo Botti e Francesco Carlo Massa (Universidade de Bergamo, Itália)

---

1. Problema e Contexto

O artigo aborda a simulação numérica de misturas de fluidos incompressíveis imiscíveis com densidade variável. O desafio central reside em formular as equações de Navier-Stokes incompressíveis de modo a garantir:

• A conservação exata de volume (o que implica que a densidade seja uma quantidade puramente advectada, sem difusão numérica artificial).
• A estabilidade em regimes dominados pela convecção.
• A robustez em relação a forças de corpo irrotacionais (independência da pressão).
• A capacidade de lidar com altos números de Reynolds e razões de densidade elevadas (instabilidades de Rayleigh-Taylor).

Métodos existentes, como DG (Galerkin Descontínuo) ou volumes finitos, muitas vezes lutam para manter a divergência nula da velocidade ponto a ponto ou exigem estabilizações que comprometem a precisão de alta ordem.

2. Metodologia

Os autores propõem uma formulação HHO (Hybrid High-Order) híbrida, combinada com esquemas de integração temporal ESDIRK (Explicit Singly Diagonal Implicit Runge-Kutta).

Discretização Espacial (HHO)

• Variáveis Híbridas: A discretização utiliza graus de liberdade definidos tanto nas células da malha quanto nas suas fronteiras (esqueleto da malha) para velocidade, pressão e densidade.
• Conformidade H(div): A escolha de um espaço de pressão híbrido, combinado com espaços polinomiais locais apropriados, permite reconstruir uma aproximação de velocidade que é H(div)-conforme e divergente nula ponto a ponto, mesmo partindo de espaços polinomiais não conformes.
• Tratamento da Densidade: Devido à conformidade H(div) e à incompressibilidade exata, a equação de conservação de massa é formulada como a derivada temporal material da densidade. Isso garante que a densidade seja transportada puramente por advecção.
• Estabilização:
  • Termos de advecção tratados com formulações antissimétricas (skew-symmetric) para conservação de energia.
  • Termos de estabilização de upwinding para o transporte de massa e momento, essenciais para regimes dominados pela convecção.
  • Estabilização da derivada temporal para a densidade, garantindo robustez quando a velocidade normal na face se anula.
• Condições de Contorno: Imposição fraca das condições de contorno de Dirichlet e Neumann.

Discretização Temporal (ESDIRK)

• O problema é tratado como um sistema de Equações Diferenciais Algébricas (DAEs) de índice 2.
• Utilizam-se esquemas ESDIRK de alta ordem, projetados especificamente para DAEs, permitindo passos de tempo implícitos e estáveis.
• A formulação temporal lida com a transição entre variáveis primitivas (velocidade, densidade) na discretização espacial e variáveis conservativas na derivação temporal.

Estratégia Algébrica

• Condensação Estática: Os graus de liberdade internos às células são eliminados localmente, deixando apenas os graus de liberdade nas faces acoplados globalmente. Isso reduz drasticamente o tamanho e a complexidade da matriz global.
• P-multigrid: Para resolver os sistemas lineares resultantes do método de Newton, utiliza-se um pré-condicionador p-multigrid, que explora a estrutura hierárquica dos espaços polinomiais para acelerar a convergência.

3. Contribuições Principais

1. Conservação Exata de Volume: A formulação garante que a conservação de volume seja imposta célula a célula até a precisão da máquina, resultando em uma advecção pura da densidade.
2. Robustez à Pressão (Pressure-Robustness): O método é projetado para que erros na representação do campo de pressão não afetem a precisão da velocidade e da densidade. Forças de corpo irrotacionais afetam apenas a pressão.
3. Estabilidade e Limites: A formulação é provada estável (dissipativa) no nível discreto para condições de contorno de Dirichlet homogêneas. Além disso, em baixas ordens (k=0), o método preserva os limites físicos da densidade (evita densidades negativas).
4. Eficiência Computacional: O uso de condensação estática reduz a memória necessária para a matriz Jacobiana, e o uso de p-multigrid melhora o desempenho dos solvers iterativos.
5. Versatilidade: O método é aplicável tanto a fluxos com baixa quanto alta razão de densidade (Atwood numbers).

4. Resultados Numéricos

Os autores validaram o método através de três casos de teste:

• Caso Kovasznay (Fluxo Estacionário):

  • Demonstrou taxas de convergência espaciais ótimas: $k+2$ para a velocidade em norma $L^2$ e $k+1$ para o gradiente de velocidade e pressão.
  • A densidade permaneceu constante (igual a 1) com erro de máquina, confirmando a advecção pura.
  • A variante HHO-$\pi^{k+1}$ mostrou maior precisão para $k=0$ em comparação com HHO-$\pi^k$.

• Solução Manufacturada de Guermond e Quartapelle (Fluxo Transiente):

  • Validou as taxas de convergência temporal e espacial.
  • Confirmou a robustez à pressão: erros de velocidade e densidade foram independentes dos erros de pressão.
  • Esquemas ESDIRK de 3 e 6 estágios atingiram as ordens de convergência temporais esperadas (2ª e 4ª ordem, respectivamente).

• Instabilidade de Rayleigh-Taylor (RTI):

  • Baixo Número de Atwood (At=0.5): Simulações em $Re=1000$ e $Re=5000$. Resultados de alta ordem ($k=6$) em malhas grossas foram comparáveis a resultados de baixa ordem ($k=1$) em malhas finas, mas com custo computacional significativamente menor (devido à condensação estática).
  • Alto Número de Atwood (At=0.75): Cenário desafiador com razão de densidade 7:1. O método com $k=0$ preservou os limites de densidade (evitando densidades negativas), garantindo estabilidade. Embora a precisão espacial fosse limitada em $k=0$, a estrutura do pluma principal foi capturada corretamente.

5. Significância e Conclusões

Este trabalho representa um avanço significativo na modelagem computacional de fluidos multifásicos incompressíveis. A principal inovação é a combinação da conformidade H(div) (garantindo incompressibilidade exata) com a advecção pura da densidade em um esquema de alta ordem.

• Impacto Prático: O método permite simular misturas de fluidos imiscíveis com alta eficiência, reduzindo o custo computacional através da condensação estática e do p-multigrid, sem sacrificar a precisão física.
• Aplicabilidade: A capacidade de lidar com altos números de Reynolds e grandes contrastes de densidade torna a ferramenta promissora para problemas de engenharia complexos, como mistura de combustíveis, dinâmica de oceanos e atmosferas.
• Futuro: Os autores indicam que o próximo passo é estender a formulação HHO-ESDIRK para o sistema Cahn-Hilliard-Navier-Stokes (CHNS) para tratar interfaces difusas em fluxos multifásicos mais realistas.

Em resumo, a formulação proposta oferece um equilíbrio raro entre precisão de alta ordem, estabilidade numérica rigorosa e eficiência computacional, superando limitações de métodos tradicionais para fluxos de densidade variável.