cfdmfFTFoam: A front-tracking solver for multiphase flows on general unstructured grids in OpenFOAM

Este artigo apresenta o cfdmfFTFoam, um novo solver de código aberto baseado no método de rastreamento de interface (FTM) integrado ao OpenFOAM, capaz de simular escoamentos multifásicos com alta precisão em malhas não estruturadas gerais, superando as limitações de software existente restrito a grades cartesianas.

O essencial

Imagine que você está tentando prever como uma gota de óleo se move dentro de um copo d'água, ou como bolhas de ar sobem em uma cerveja. Para os cientistas e engenheiros, isso é um quebra-cabeça matemático enorme, porque a fronteira entre a água e o óleo (ou o ar) está sempre mudando de forma, esticando e quebrando.

Este artigo apresenta uma nova ferramenta de software chamada cfdmfFTFoam, criada por Ehsan Amani, que funciona como um "super-olho" para ver e simular essas misturas complexas.

Aqui está a explicação simplificada, usando analogias do dia a dia:

1. O Problema: A "Dança" das Fronteiras

Pense em duas equipes jogando em um campo. Uma equipe é a água, a outra é o óleo. O problema é que a linha que separa as duas equipes (a interface) não é fixa; ela se move, se deforma e interage.

• Os métodos antigos (como o VOF): Eles tratavam a fronteira como uma névoa ou uma pintura borrada na tela. É fácil de calcular, mas a "pintura" às vezes escorre ou desaparece, perdendo a precisão da forma exata da gota.
• O método deste novo software (FTM - Rastreamento de Frente): Em vez de uma névoa, imagine que a fronteira é uma malha de rede de pesca feita de triângulos (como uma rede de pesca 3D). Essa rede é "inteligente": ela se move exatamente onde a água e o óleo se encontram. É como se você tivesse um grupo de dançarinos segurando as pontas de uma rede, e eles se movem perfeitamente para seguir a forma da gota.

2. A Inovação: Conectando o Mundo Real ao Computador

O grande desafio que o autor resolveu foi fazer essa "rede de pesca" (que é muito detalhada) conversar com o "chão" do computador (a grade de cálculo).

• A Analogia: Imagine que você tem um mapa de uma cidade feito de blocos de Lego (a grade do computador) e uma rede de pesca flutuando acima dela. O software precisa saber exatamente onde cada nó da rede toca os blocos de Lego para calcular a pressão e o movimento.
• O que o cfdmfFTFoam faz: Ele é o "tradutor" perfeito. Ele pega os dados da rede de pesca (Lagrangiana) e os passa para os blocos de Lego (Euleriana) sem perder informações, mesmo que a cidade (o computador) seja feita de blocos de formatos estranhos e irregulares (grades não estruturadas), e não apenas quadrados perfeitos.

3. As Ferramentas Mágicas (Algoritmos)

O software vem com um "kit de ferramentas" para manter a rede de pesca funcionando perfeitamente:

• Correção de Volume (VC2): Às vezes, a rede de pesca pode "vazar" um pouco, fazendo a gota parecer que está perdendo água. O software tem um "bico de reparo" que ajusta a posição dos nós da rede para garantir que a gota mantenha o tamanho exato, como se estivesse enchendo um balão que está perdendo ar.
• Remeshing (Reconstrução da Rede): Se a gota se estica muito, a rede de pesca pode ficar com furos grandes ou triângulos tortos. O software é como um "alfaiate digital": ele corta os triângulos grandes em pedaços menores (refinamento) ou junta os pequenos em pedaços maiores (coarsening) para que a rede sempre fique bem ajustada à forma da gota.
• Tensão Superficial: É a "pele" da gota que a mantém unida. O software calcula exatamente quanta força essa pele exerce em cada ponto, garantindo que a gota não se desfaça artificialmente.

4. Por que isso é importante?

Antes deste software, era muito difícil simular essas misturas complexas em computadores que usam grades irregulares (como em turbinas de avião ou vasos sanguíneos) usando esse método de "rede de pesca" pura. A maioria dos softwares gratuitos só funcionava em grades retas e quadradas.

O cfdmfFTFoam é:

1. Gratuito e de Código Aberto: Qualquer pessoa pode baixar, estudar e melhorar.
2. Preciso: Ele evita que as gotas se "quebrem" sozinhas na simulação, algo comum em outros métodos.
3. Versátil: Funciona em qualquer formato de grade, permitindo simulações em geometrias complexas do mundo real.

Resumo Final

Pense no cfdmfFTFoam como um novo tipo de "óculos de realidade aumentada" para engenheiros. Em vez de ver apenas uma névoa borrada de onde está a água e o óleo, eles agora podem ver uma rede de pesca 3D precisa que segue cada movimento da gota, mesmo em ambientes complexos e irregulares. Isso permite prever com muito mais segurança como fluidos se comportam em motores, no corpo humano ou na indústria, economizando tempo e dinheiro em testes reais.

Resumo técnico

1. Problema e Contexto

A previsão confiável de escoamentos multifásicos complexos (como deformação de bolhas, colisão de gotículas, pulverização e ebulição) exige métodos de simulação totalmente resolvidos (FRS - Fully Resolved Simulation). Embora existam várias abordagens de captura de interface (como VOF e Level-Set) e rastreamento (como ALE e Imersão de Fronteira), o método de Rastreamento de Fronteira (FTM - Front-Tracking Method) destaca-se pela sua capacidade de calcular com precisão vetores normais, curvaturas e forças de tensão superficial, além de evitar erros de dissipação na advecção da interface.

No entanto, o artigo identifica uma lacuna crítica: a escassez de softwares de código aberto (Open-Source) e acesso livre (Open-Access) que implementem um FTM puro em malhas eulerianas não estruturadas gerais.

• Soluções existentes (como PARIS, FronTier++, TrioCFD) são limitadas a malhas cartesianas estruturadas.
• Soluções em OpenFOAM (como lentFoam) utilizam abordagens híbridas (Level-Set + FTM) que não possuem uma representação lagrangiana pura da superfície da interface.

O objetivo deste trabalho é preencher essa lacuna, apresentando e validando o solver cfdmfFTFoam.

2. Metodologia

O solver cfdmfFTFoam foi desenvolvido integrando o código de rastreamento de fronteira Ftc3D (originalmente em Fortran, desenvolvido por Gretar Tryggvason e Sadegh Dabiri) ao ecossistema OpenFOAM v9 (C++).

Arquitetura e Algoritmos Principais:
O método utiliza uma malha lagrangiana de superfície (triangular) que se move através de uma malha euleriana volumétrica (não estruturada) onde as equações de Navier-Stokes são resolvidas. O acoplamento entre as duas malhas e a evolução da interface são realizados através de seis etapas sequenciais em cada passo de tempo:

1. Correção de Volume (VC2): Como o FTM não conserva massa/volume inerentemente, um algoritmo corrige a posição dos vértices da fronteira resolvendo uma equação cúbica para ajustar o volume encerrado, garantindo conservação de massa.
2. Remalhamento (Remeshing): Controla a qualidade e o número de elementos da malha da fronteira através de três sub-rotinas:
   • Refinamento: Divide elementos grandes ou com alto aspecto.
   • Afinamento (Coarsening): Funde elementos pequenos.
   • Suavização (Smoothing): Remove ondulações numéricas (utilizando algoritmos como TSUR3D e VCS).
3. Cálculo da Tensão Superficial: Utiliza um método direto baseado em elementos (directElementBased) para calcular as forças de tensão superficial em cada elemento triangular.
4. Comunicação Fronte-Campo (Front-to-Field): Transfere dados da malha lagrangiana para a euleriana. O trabalho implementa o Método de Partículas com Núcleo Reprodutor (RKPM - Reproducing Kernel Particle Method) para distribuir as forças de tensão superficial e construir o campo de densidade/viscosidade, preservando momentos de força e torque.
5. Construção da Função Indicadora: Define qual fluido ocupa qual região. São implementados dois métodos:
   • Poisson (PE): Resolve uma equação de Poisson para obter uma função indicadora suave (preferido para malhas não estruturadas).
   • Transformação do Ponto Mais Próximo (CPT): Baseada na distância geométrica.
6. Advecção da Fronteira: Os vértices da malha lagrangiana são movidos interpolando a velocidade do fluido euleriano, utilizando ferramentas nativas de rastreamento de partículas do OpenFOAM.

O solver suporta execução serial e paralela (decomposição de domínio) em malhas não estruturadas gerais (hexaédricas, tetraédricas, prismas, poliedros).

3. Contribuições Chave

• Primeiro Solver FTM Puro em OpenFOAM: É a primeira ferramenta de código aberto que implementa um método de rastreamento de fronteira puro (sem Level-Set híbrido) em malhas não estruturadas gerais.
• Integração Ftc3D-OpenFOAM: Sucesso na integração de um código legado complexo (Ftc3D) na arquitetura moderna e paralela do OpenFOAM, lidando com a comunicação de dados entre malhas de topologias diferentes.
• Algoritmos Avançados: Implementação de correção de volume robusta, múltiplos esquemas de suavização de malha e o uso de RKPM para transferência de dados, superando limitações de métodos anteriores.
• Código Aberto e Acessível: Disponibilização do código sob licença GPLv3 no repositório GitHub e na biblioteca CPC, facilitando a pesquisa futura e o desenvolvimento de novos algoritmos na área.

4. Resultados e Validação

O solver foi validado através de uma série de benchmarks padrão da literatura:

• Deformação 3D de uma Gota (Campo de Velocidade Analítico):

  • O solver demonstrou alta precisão na advecção da interface, com erro de advecção de apenas 1%, comparado a 36-41% para métodos VOF (MULES e isoAdvector).
  • Evitou rupturas de interface indesejadas e não controladas que ocorrem em métodos de captura de interface sob grandes deformações.
  • A convergência da malha mostrou uma ordem de convergência aparente de ~2 e um índice de convergência de malha (GCI) abaixo de 1%.
  • A correção de volume (VC2) manteve a conservação de massa com erro < 0,2%.

• Gota Estagnada (Balanço de Pressão):

  • O solver previu a diferença de pressão interna/externa com erro de 0,046%, superando significativamente o solver VOF padrão (interFoam), que apresentou erro de 12,6%.
  • As correntes parasitas (parasitic currents) foram mantidas em níveis muito baixos ($C_{as} \approx 2.1 \times 10^{-3}$).

• Gota em Escoamento de Poiseuille (Regime de Stokes):

  • Demonstrou capacidade de operar em diversas topologias de malha (hexaédrica, prismática, poliedral, tetraédrica) com boa concordância com soluções analíticas.
  • Apresentou eficiência de paralelização próxima à do solver padrão interFoam.

• Ascensão de Bolha Taylor e Bolha Livre:

  • Simulações de bolhas Taylor em tubos verticais e bolhas livres com altos números de Reynolds e razões de densidade/viscosidade mostraram concordância qualitativa e quantitativa com dados experimentais.
  • Identificou-se que a precisão em altos números de densidade pode ser melhorada ajustando o raio do núcleo (parâmetro $\alpha$) do RKPM, embora isso exija tratamento de regularização de Tikhonov para evitar singularidades.

5. Significado e Conclusão

O cfdmfFTFoam representa um avanço significativo na simulação computacional de fluidos (CFD) de escoamentos multifásicos. Ao fornecer uma ferramenta de código aberto para FTM puro em malhas não estruturadas, ele permite:

1. Maior Precisão Física: Melhor captura da física de coalescência, ruptura e forças interfaciais em comparação com métodos de captura de interface (VOF/Level-Set) em malhas complexas.
2. Flexibilidade Geométrica: Capacidade de simular problemas complexos em geometrias que não podem ser facilmente mapeadas em malhas cartesianas.
3. Base para Pesquisa Futura: O código serve como uma plataforma sólida para o desenvolvimento de novas estratégias híbridas, dinâmica de linhas de contato, transporte de calor/massa e melhoria de algoritmos de FTM.

Em suma, o trabalho elimina uma barreira técnica importante, democratizando o acesso a métodos de rastreamento de fronteira de alta fidelidade para a comunidade científica e industrial.