Color-gradient lattice Boltzmann modeling of wetting boundary condition on curved solid boundaries

Este artigo introduz uma condição de contorno de molhabilidade para superfícies sólidas curvas no método de lattice Boltzmann de gradiente de cor ao atualizar parâmetros de ordem em nós fantasmas, um esquema validado em hardware GPU para lidar eficazmente com grandes contrastes de densidade e viscosidade, minimizando correntes espúrias e reproduzindo com precisão tanto os comportamentos de linha de contato estática quanto dinâmica.

O essencial

Imagine que você está tentando simular como uma gota de água se comporta quando atinge uma superfície curva, como uma gota de chuva caindo em uma folha ou uma bolha deslizando por um vidro curvo. Para fazer isso em um computador, cientistas usam um método chamado Método de Lattice Boltzmann. Pense neste método como uma grade gigante e invisível de pequenos azulejos cobrindo a tela do computador. Cada azulejo contém um pouco de informação de "fluido", e o computador atualiza esses azulejos passo a passo para ver como o fluido se move.

A parte complicada é a condição de contorno — especificamente, como o fluido se comporta ao tocar uma parede sólida. No mundo real, a água não apenas para abruptamente em uma parede; ela forma um ângulo específico (chamado ângulo de contato) dependendo se a superfície é molhada (como um vidro limpo) ou seca (como um carro encerado).

O Problema: O "Fantasma" na Máquina

Na simulação de computador, a parede sólida não é uma linha suave; é irregular porque é feita de azulejos de uma grade quadrada. Para fazer a matemática funcionar, o computador precisa saber o que o fluido está fazendo dentro da parede sólida, embora não haja fluido lá. Esses pontos imaginários dentro da parede são chamados de "nós fantasmas" (ghost nodes).

Métodos anteriores para dizer a esses nós fantasmas o que fazer tinham algumas falhas:

• Eles às vezes criavam "correntes fantasmas" (velocidades espúrias) onde o fluido parecia se mover sozinho, sem qualquer força.
• Eles tinham dificuldade com superfícies curvas, agindo muitas vezes como se tivessem sido projetados apenas para paredes planas.
• Eles às vezes exigiam uma matemática especial e complicada apenas para lidar com um ângulo neutro (onde a água nem se espalha nem se agrupa).

A Solução: Uma Nova Regra para os Fantasmas

Os autores deste artigo introduziram uma regra nova e mais simples para esses nós fantasmas.

A Analogia: Imagine que o fluido tem um "humor" (representado por uma cor, de 0 para gás a 1 para líquido). No mundo real, esse humor muda suavemente de gás para líquido conforme você atravessa a superfície.

• Método Antigo: Era como tentar adivinhar o humor de uma pessoa parada atrás de uma parede gritando um palpite aleatório.
• Novo Método: Os autores perceberam que, se você souber o "humor" da pessoa que está logo fora da parede (no fluido), você pode matematicamente estender essa curva de humor suave através da parede até o nó fantasma. Eles simplesmente perguntam: "Se o fluido quer formar um ângulo de 45 graus aqui, qual deve ser o humor do nó fantasma para que isso aconteça?"

Esta nova regra é como uma ponte contínua. Ela estende a forma natural da gota de fluido até chegar e avançar ligeiramente para dentro da parede sólida, garantindo que o ângulo que a gota faz com a parede seja exatamente o que o cientista solicitou.

O Que Eles Testaram

Para provar que sua nova regra funciona, eles realizaram várias simulações em um chip de computador muito poderoso (um NVIDIA A100 GPU):

1. A Gota Estática: Eles colocaram uma gota de água em uma placa plana e em um cilindro curvo. Eles verificaram se a gota se estabilizava no ângulo exato que solicitaram.
   • Resultado: A nova regra deles foi mais precisa do que o melhor método anterior, especialmente quando o ângulo era muito agudo (como uma gota que se agrupa) ou muito plano (como uma gota que se espalha).
2. A Partícula Flutuante: Eles simularam um cilindro flutuando na fronteira entre óleo e água.
   • Resultado: O método deles calculou a posição da linha d'água com mais precisão do que antes.
3. A Gota em Queda: Eles simularam uma gota caindo e atingindo um cilindro, observando-a espirrar e se espalhar.
   • Resultado: A gota se comportou de forma realista, e a nova regra não causou movimentos estranhos e falsos no fluido.

Principais Conclusões

• Precisão: O novo método lida muito melhor com superfícies curvas, mantendo o ângulo correto do fluido, quer a parede seja plana ou arredondada.
• Estabilidade: Ele cria muito pouca "interferência falsa" (correntes espúrias) na simulação, o que faz o fluido parecer mais natural.
• Simplicidade: Evita a necessidade de uma matemática especial e complicada quando o ângulo de contato é exatamente 90 graus (neutro), o que era um problema para métodos anteriores.
• Velocidade: Ao utilizar chips de computador modernos (GPUs) e um estilo de programação específico, eles fizeram as simulações rodarem muito rápido. Eles descobriram que usar um formato de número com um pouco menos de precisão (precisão simples) fez o computador rodar duas vezes mais rápido sem estragar os resultados para a maioria dos testes.

Em resumo, os autores construíram um melhor "livro de regras" para como as simulações de computador lidam com a borda onde o líquido encontra uma parede sólida, fazendo com que as gotas digitais pareçam e ajam mais como as reais, mesmo em superfícies curvas.

Resumo técnico

Resumo Técnico: Modelagem de Lattice Boltzmann com Gradiente de Cor para Condições de Contorno de Molhabilidade em Fronteiras Sólidas Curvas

Definição do Problema
A modelagem precisa de fenômenos de molhabilidade é crítica para simular processos tecnologicamente relevantes, como o fluxo em meios porosos, microfluídica e a secagem de filmes coloidais. Em estruturas de Lattice Boltzmann Method (LBM) que utilizam interfaces difusas, a implementação de condições de contorno de molhabilidade em fronteiras sólidas curvas apresenta desafios significográficos. As abordagens existentes frequentemente sofrem com limitações, incluindo a introdução de fontes ou sumidouros artificiais de massa, a geração de velocidades espúrias (parasitas) próximas às interfaces, restrições a geometrias planas e inconsistências na imposição simultânea das condições de não-escorregamento (no-slip) e de molhabilidade. Além disso, muitos métodos exigem tratamento especial para ângulos de contato neutros ($\theta = 90^\circ$) ou dependem de extrapolações lineares que introduzem erros próximos às fronteiras sólidas.

Metodologia
Os autores introduzem uma condição de contorno de molhabilidade projetada especificamente para fronteiras sólidas curvas dentro de um framework LBM de gradiente de cor bidimensional. O núcleo da metodologia envolve a atualização do parâmetro de ordem (campo de cor, $\phi$) em "nós fantasmas" (ghost nodes) localizados dentro da fase sólida.

• Regra de Atualização do Nó Fantasma: Em vez de prescrever um valor constante ou utilizar extrapolação linear, o esquema estende o perfil de equilíbrio da cor para dentro da fase sólida. O valor do campo de cor em um nó fantasma ($x_G$) é calculado com base no campo de cor em um ponto correspondente no domínio do fluido ($x_F$), localizado a uma distância $\Delta x$ ao longo da normal sólida.
• Formulação Matemática: Assumindo que a distribuição de cor está próxima do equilíbrio, a derivada do parâmetro de ordem na direção da normal sólida é derivada do perfil de equilíbrio e do ângulo de contato $\theta$ prescrito. A integração desta relação resulta em uma regra de atualização não linear:
  $$\phi_G = \frac{\phi_F}{\phi_F + (1 - \phi_F) \exp(\varepsilon)}$$
  onde $\varepsilon = -4\Delta x \cos \theta / W$, e $W$ é a largura da interface.
• Características Principais: Esta formulação garante que o campo de cor do nó fantasma permaneça limitado entre 0 e 1, evitando a transferência de massa artificial. Crucialmente, a regra não requer tratamento especial quando $\theta = 90^\circ$ (onde $\varepsilon = 0$), ao contrário de esquemas anteriores. O método é implementado utilizando um framework baseado em JAX executado em GPUs NVIDIA A100, utilizando fusão de kernels (kernel fusion) e diferenciação automática para alto desempenho.

Contribuições Principais e Resultados
O esquema proposto foi validado contra soluções analíticas e o esquema de referência de Fakhari et al. [20] através de diversos casos de teste:

1. Gota Estática sobre uma Placa Plana: O esquema produziu consistentemente erros menores nos ângulos de contato medidos em comparação com Fakh��{a}ri et al. em uma gama de $30^\circ$ a $150^\circ$. Também gerou correntes espúrias menores ($5.74 \times 10^{-7}$ a $7 \times 10^{-7}$ unidades de rede) comparado ao esquema de referência.
2. Gota sobre um Cilindro Sólido: Para uma gota situada sobre um cilindro curvo, o método proposto demonstrou melhor precisão na predição da altura do ápice da gota ($h_{max}$), particularmente em ângulos de contato extremos ($30^\circ$ e $150^\circ$), onde os erros foram reduzidos de $\sim 0,76\%$ para $\sim 0,35\%$ em comparação com o esquema de referência.
3. Partícula Fixa na Interface Líquido-Gás: Em simulações de um arranjo periódico de cilindros presos em uma interface, o esquema proposto reduziu os erros no deslocamento vertical da interface para $\sim 9\%$ em ângulos extremos, comparado a $\sim 14\%$ para o esquema de referência.
4. Impacto Dinâmico de Gota: O esquema simulou com sucesso o impacto de gotas em um cilindro em altos números de Reynolds e Bond, capturando comportamentos dinâmicos complexos, como o espalhamento da gota em superfícies hidrofílicas e a divisão em superfícies hidrofóbicas. Os resultados alinharam-se bem com as leis de escala analíticas para a espessura do filme.
5. Desempenho Computacional: A implementação em GPU mostrou que computações de precisão simples (FP32) foram pelo menos duas vezes mais rápidas que as de precisão dupla (FP64), mantendo a precisão para casos estáticos e dinâmicos moderados. No entanto, o FP32 tornou-se instável para o teste de impacto de gota de alto número de Reynolds/Weber, necessitando de FP64 para esses cenários de alta energia.

Significância e Alegações
O artigo afirma que a condição de contorno de molhabilidade proposta oferece uma alternativa robusta e precisa para simular fluxos bifásicos imiscíveis em fronteiras sólidas curvas. Ao estender o perfil de equilíbrio para dentro da fase sólida via uma regra específica de atualização de nó fantasma, o método:

• Mantém a capacidade do modelo de lidar com grandes contrastes de densidade e viscosidade.
• Produz correntes espúrias relativamente pequenas (ordens de magnitude menores que modelos de pseudopotencial populares).
• Evita as singularidades numéricas associadas a ângulos de contato neutros encontradas em outros métodos.
• Demonstra precisão superior sobre os esquemas existentes (especificamente Fakhari et al.) tanto para linhas de contato estáticas quanto dinâmicas em geometrias curvas.

Os autores observam que, embora a implementação atual seja restrita à execução em GPU única e duas dimensões, a física subjacente é diretamente aplicável a modelos de campo de fase (phase-field) e de energia livre, podendo ser estendida para três dimensões. O trabalho fornece uma base para futuras simulações envolvendo partículas sólidas em interfaces líquido-gás, particularmente quando acopladas a avanços nos cálculos de força de capilaridade da linha de contato e esquemas de histerese.