A hybrid Volume of Fluid Phase-Field method for Direct Numerical Simulations of soluble surfactant-laden interfacial flows

Este artigo apresenta um método híbrido Volume de Fluido-Campo de Fase com refinamento adaptativo de malha para simulações numéricas diretas de escoamentos carregados com surfactante solúvel, que captura com precisão o acoplamento entre transporte no volume e na interface para demonstrar como tensões de Marangoni alteram significativamente a dinâmica de ascensão de bolhas em geometrias tridimensionais.

O essencial

A Visão Geral: Uma Nova Maneira de Simular Água "Sabonosa"

Imagine que você está assistindo a uma bolha subindo através de um copo de água. Se a água for pura, a bolha sobe em linha reta, oscilando um pouco, mas movendo-se rapidamente. Mas se você adicionar sabão (tensoativos) à água, a bolha se comporta de maneira diferente. Ela pode desacelerar, oscilar mais ou até mesmo mudar seu caminho.

Isso acontece porque as moléculas de sabão adoram grudar na superfície da bolha. À medida que a bolha se move, essas moléculas são empurradas, criando uma "tensão" desigual na pele da bolha. Essa tensão desigual age como mãos invisíveis empurrando a bolha em direções diferentes, alterando como ela se move.

O Problema:
Simular isso em um computador é incrivelmente difícil. É como tentar filmar uma bolha de sabão com uma câmera que possui duas configurações conflitantes:

1. A Câmera Nítida: Precisa ver a borda da bolha como uma linha fina como uma lâmina (para calcular pressão e forma).
2. A Câmera Embaçada: Precisa ver as moléculas de sabão se espalhando suavemente sobre essa borda (para calcular como o sabão se move).

A maioria dos métodos computacionais força você a escolher uma configuração de câmera, tornando a simulação ou fisicamente imprecisa ou computacionalmente impossível de executar para formas 3D complexas.

A Solução:
Os autores deste artigo construíram um método híbrido. Pense nele como uma simulação de "tela dividida" que usa o melhor dos dois mundos simultaneamente:

• A Borda Nítida (Volume de Fluido): Eles usam um método que mantém a borda da bolha nítida e conserva perfeitamente a quantidade de líquido (como um contorno de alta definição).
• O Sabão Suave (Campo de Fase): Eles usam uma segunda camada "embaçada" que atua como uma estrada suave para as moléculas de sabão viajarem. Isso permite que o sabão se mova naturalmente entre a água e a superfície da bolha sem ficar preso ou perdido.

Como Funciona: A Analogia do "Controlador de Tráfego"

Para fazer isso funcionar, os autores criaram um sistema de tráfego digital para as moléculas de sabão:

1. A Estrada (A Interface): A superfície da bolha é uma estrada movimentada. As moléculas de sabão são os carros.
2. As Rampas de Entrada e Saída (Adsorção/Desorção):
   • Adsorção: Moléculas de sabão da água (o volume) querem pular para a estrada (a superfície da bolha).
   • Desorção: Moléculas de sabão ficam cansadas e pulam da estrada de volta para a água.
   • O novo método calcula exatamente quantos carros pulam para dentro ou para fora a cada instante, garantindo que nenhum carro desapareça ou apareça do nada.
3. O Engarrafamento (Tensão de Marangoni): Quando muitos carros de sabão se acumulam em um ponto na bolha, esse ponto fica "pegajoso" (alta tensão). A pele da bolha tenta se afastar do ponto pegajoso, criando uma força que desacelera a bolha ou faz com que ela oscile. A simulação captura esse cabo de guerra perfeitamente.

O Que Eles Testaram (Os "Testes de Direção")

Antes de deixar seu novo carro dirigir na estrada, eles o levaram para uma escola de direção com três testes específicos:

1. O Teste de Estiramento (Esfera em Expansão): Eles inflaram uma bolha coberta de sabão. Verificaram se o sabão se espalhava uniformemente à medida que a bolha ficava maior. A simulação combinou perfeitamente com a matemática.
2. O Teste de Rotação (Bolha Giratória): Eles giraram uma bolha com sabão nela. Verificaram se o sabão se movia corretamente ao redor do círculo sem vazar. Novamente, a simulação foi perfeita.
3. O Teste de Troca (Parede Plana): Eles observaram o sabão se mover da água para uma parede plana e de volta. Eles testaram três cenários:
   • Apenas pulando para dentro: O sabão gruda? Sim.
   • Apenas pulando para fora: O sabão sai? Sim.
   • Ambos: Ele encontra um equilíbrio? Sim.

O Evento Principal: A Bolha Subindo

Finalmente, eles deixaram seu novo método simular uma bolha subindo em um tanque 3D de água.

• A Bolha Limpa: Subiu relativamente rápido e em linha reta.
• A Bolha de Sabão "Insolúvel": O sabão estava preso à superfície e não podia sair. Isso criou um forte "engarrafamento" na parte traseira da bolha, desacelerando-a significativamente.
• A Bolha de Sabão "Solúvel" (A Realidade): É aqui que o novo método brilha. O sabão podia pular para dentro e para fora da bolha.
  • Se o sabão pulasse para fora facilmente (alta "desorção"), a bolha se comportava quase como uma limpa.
  • Se o sabão pulasse para dentro facilmente (alta "adsorção"), a bolha se comportava como a versão de sabão preso.
  • No meio, a bolha mostrou uma dança complexa: desacelerou, mudou seu caminho e deixou um "rastro" de sabão na água atrás dela enquanto subia.

Por Que Isso Importa (Segundo o Artigo)

Os autores afirmam que este método é robusto, escalável e preciso.

• Robusto: Não falha quando a bolha assume formas estranhas ou se fragmenta.
• Escalável: Pode ser executado em supercomputadores para lidar com simulações 3D enormes e complexas de forma eficiente.
• Preciso: Prevê corretamente a velocidade com que as bolhas sobem e como elas oscilam, correspondendo à física do mundo real.

Em resumo: Eles construíram um novo motor digital que finalmente pode simular como as bolhas de sabão se comportam no espaço 3D, lidando com a dança complicada entre a forma da bolha e as moléculas de sabão movendo-se para dentro e para fora de sua pele, tudo sem perder precisão ou travar o computador.

Resumo técnico

Resumo Técnico: Um Método Híbrido Volume-of-Fluid e Phase-Field para Escoamentos Interfaciais com Surfactantes Solúveis

Enunciação do Problema
A simulação numérica precisa de escoamentos interfaciais envolvendo surfactantes solúveis apresenta um desafio formidável devido à natureza multiescala e multipísica do problema. Diferentemente dos surfactantes insolúveis, as espécies solúveis são governadas por um sistema acoplado envolvendo convecção-difusão no volume, transporte interfacial e processos de troca cinética (adsorção/desorção). Essas interações geram campos de tensão superficial que variam espacial e temporalmente, levando a tensões de Marangoni que alteram significativamente a topologia do escoamento, a estabilidade e o movimento. Os métodos numéricos existentes enfrentam limitações distintas:

• Métodos de Integral de Fronteira/Elementos: Restritos a regimes de escoamento de Stokes e com dificuldade em lidar com mudanças de topologia (ruptura/coalescência).
• Rastreamento de Fronteira (Front-Tracking): Alta precisão, mas implementação complexa em 3D, exigindo tratamento topológico explícito e enfrentando problemas de balanceamento de carga em computação paralela.
• Métodos Level-Set: Lidam naturalmente com mudanças de topologia, mas frequentemente sofrem de problemas de conservação de massa, exigindo reinicialização e correções post-hoc.
• Métodos Phase-Field: Lidam naturalmente com mudanças de topologia, mas frequentemente dependem de suposições de equilíbrio termodinâmico que podem ser violadas por convecção forte, levando à não conservação de massa.
• Métodos Volume-of-Fluid (VoF): Excelentes para conservação de massa e interfaces nítidas, mas tradicionalmente têm dificuldade em definir operadores diferenciais superficiais necessários para o transporte de surfactantes em uma interface de espessura zero.

Há uma necessidade específica de um método baseado em VoF geométrico que suporte Refinamento Adaptativo de Malha (AMR), lide com surfactantes solúveis com troca dinâmica volume-interface e opere dentro de um framework escalável e de código aberto.

Metodologia
Os autores apresentam um framework numérico híbrido implementado no código de código aberto Basilisk, que combina um método geométrico Volume-of-Fluid (VoF) com uma abordagem difusa Phase-Field.

1. Representação da Interface:

   • VoF ($c$): Um campo escalar representando a fração volumétrica de um fluido, utilizado para rastreamento conservativo da interface e solução de momento. Garante conservação exata de massa e suporta reconstrução geométrica.
   • Phase-Field ($\phi$): Uma variável auxiliar que evolui via uma formulação de Interface Difusa Conservativa Precisa baseada em Allen-Cahn (ACDI). Serve como um transportador suave para o transporte de surfactantes, fornecendo vetores normais regularizados, curvatura e gradientes superficiais sem a necessidade de cálculos explícitos de derivadas superficiais em uma interface nítida.
   • Acoplamento: O Phase-Field é periodicamente reinicializado a partir do campo VoF usando uma função de distância assinada para manter a consistência.

2. Equações Governantes:

   • Movimento do Fluido: Equações de Navier-Stokes incompressíveis com densidade e viscosidade variáveis espacialmente. As forças de tensão superficial são incorporadas como forças de corpo volumétricas: uma força capilar isotrópica ($\sigma \kappa \mathbf{n} \delta_s$) e uma força tangencial de Marangoni ($\nabla_s \sigma \delta_m$).
   • Transporte de Surfactante: Duas variáveis escalares são resolvidas: concentração no volume ($F$) e concentração interfacial ($f$, uma representação volumétrica de $\Gamma$).
     • A equação de transporte interfacial é regularizada para evitar derivadas superficiais, convertendo a concentração superficial $\Gamma$ em um campo volumétrico $f = \Gamma \delta_\phi$ localizado dentro da interface difusa.
     • Cinética: Adsorção e dessorção são modeladas via termos fonte/sumidouro regularizados ($J$) localizados na interface, ligando as concentrações no volume e na interface. O modelo suporta equações de estado gerais (por exemplo, isotermas de Langmuir) para relacionar a tensão superficial à concentração.

3. Discretização Numérica:

   • Espacial: Discretização por volumes finitos em grades quadtree/octree dinamicamente refinadas (AMR). A resolução é concentrada próximo à interface para resolver gradientes íngremes na concentração e no escoamento.
   • Temporal: Utiliza-se splitting de operadores. A advecção é tratada explicitamente (esquema geométrico não dividido Bell–Collela–Glaz), enquanto difusão, anti-difusão e termos fonte são tratados implicitamente usando um solucionador Poisson multigrid modificado.
   • Estabilidade: Uma condição Courant–Friedrichs–Lewy (CFL) é imposta com base em uma velocidade efetiva que inclui contribuições tanto convectivas quanto anti-difusivas para garantir estabilidade na presença de gradientes de interface fortes.

Contribuições Principais

• Framework Híbrido: Desenvolvimento de um método que retém as vantagens de conservação de massa e interface nítida do VoF, enquanto utiliza uma variável Phase-Field exclusivamente como transportador suave para o transporte de surfactantes, evitando as complexidades geométricas de derivadas superficiais em interfaces lineares por partes.
• Capacidade para Surfactantes Solúveis: Extensão do trabalho anterior VoF-Phase-Field (limitado a surfactantes insolúveis) para o caso totalmente solúvel, incorporando termos acoplados de advecção-difusão no volume e troca cinética.
• Escalabilidade e Adaptatividade: O framework é totalmente adaptativo (AMR) e suporta configurações 2D, axisimétricas e 3D com alta escalabilidade paralela, permitindo a resolução de gradientes multiescala tanto nos campos do volume quanto nos interfaciais.
• Implementação de Código Aberto: O método é implementado no Basilisk, fornecendo uma ferramenta reprodutível e acessível para a comunidade.

Resultados e Validação
O artigo valida o método através de uma série de testes de referência e casos de aplicação:

1. Validação sem Escoamento:

   • Esfera em Expansão (3D): Verificou a recuperação de soluções analíticas para diluição de surfactante insolúvel durante expansão isotrópica.
   • Interface Rotativa (2D): Confirmou convergência de segunda ordem para advecção-difusão de surfactante insolúvel em uma interface em movimento.
   • Adsorção/Desorção em Interface Plana: Validou o solver acoplado volume-interface contra soluções analíticas para adsorção de fluxo constante, dessorção pura e cinética completa de Langmuir. O método captura corretamente a transição para o equilíbrio e a formação de camadas de depleção/enriquecimento no volume.

2. Bolha Ascendente Axisimétrica:

   • Limite Insolúvel: Demonstrou a redução da velocidade terminal de subida devido às tensões de Marangoni, com resultados convergindo para o limite de interface limpa à medida que a concentração de surfactante diminui.
   • Dinâmica Solúvel: Variou sistematicamente o número de Biot (dessorção) e o número de Damköhler (adsorção).
     • Alto Biot (Dessorção Rápida): O sistema aproxima-se do limite de interface limpa à medida que os surfactantes são removidos rapidamente da interface.
     • Alto Damköhler (Adsorção Rápida): O sistema aproxima-se do limite de superfície incompressível (insolúvel) à medida que a interface satura rapidamente.
   • Estrutura do Escoamento: Visualizações mostraram a formação de "capuzes de estagnação" de surfactante na parte traseira da bolha, que suprimem a recirculação interna e modificam a estrutura da esteira.

3. Bolha Ascendente Tridimensional:

   • Simulações de uma bolha subindo em 3D com surfactantes solúveis capturaram fenômenos complexos incluindo deformação de forma, desvios de trajetória lateral (movimento helicoidal) e a formação de plumas de arrasto devido à dessorção.
   • O método recuperou com sucesso tanto os limites assintóticos limpos quanto insolúveis em 3D, demonstrando robustez no tratamento de padrões de escoamento não estacionários e geometrias deformáveis.

Significância e Alegações
O artigo alega fornecer um framework numérico robusto e versátil que preenche a lacuna entre modelos clássicos de surfactantes insolúveis e formulações solúveis realistas. Sua significância reside em:

• Tratamento Unificado: Lida com acoplamento não linear, interfaces em evolução e configurações totalmente tridimensionais dentro de uma única plataforma de código aberto.
• Fidelidade Física: Captura a interação complexa entre hidrodinâmica, transporte volume/interfacial e tensões de Marangoni sem depender de modelos de fechamento de turbulência ou aproximações empíricas.
• Aplicabilidade Prática: O método é adequado para escoamentos com superfície livre envolvendo reconexão de interface (por exemplo, ruptura de jatos, ondas quebrando) na presença de surfactantes solúveis relevantes industrial e biologicamente.

Os autores notam modestamente que o framework atual é limitado a números de Péclet moderados devido à difusão numérica inerente na discretização. Eles identificam a resolução precisa de gradientes interfaciais nítidos no regime de alto Péclet (dominado por convecção) como um desafio para trabalhos futuros, visando desenvolver formulações que minimizem a difusão numérica.