A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries

Este artigo apresenta um método numérico 3D conservador e de interface nítida, baseado no Volume de Fluidos (VOF) e em fronteiras embebidas, que resolve dinâmicas de linhas de contato com histerese em geometrias complexas através de um esquema de reconstrução para células mistas, uma estratégia de redistribuição para eliminar restrições de passo de tempo e uma técnica inovadora de imposição de ângulo de contato.

O essencial

Imagine que você está tentando prever exatamente como uma gota de água se comporta quando cai em uma superfície estranha, como um copo de vidro com riscos, uma folha de árvore com nervuras ou até mesmo uma tela de mosquito. O problema é que, na física computacional, simular isso é como tentar montar um quebra-cabeça 3D onde as peças mudam de tamanho e forma o tempo todo, e algumas delas são cortadas ao meio por obstáculos invisíveis.

Este artigo, publicado no Journal of Computational Physics em 2026, apresenta uma nova "receita" matemática para resolver esse problema com precisão cirúrgica. Vamos descomplicar as duas grandes inovações dos autores usando analogias do dia a dia:

1. O Problema da "Gota que Vaza" (Conservação de Massa)

A Situação:
Imagine que você tem uma caixa de água (o computador) dividida em pequenos cubos. A água ocupa alguns cubos totalmente, outros parcialmente. Agora, imagine que você coloca uma pedra (o sólido) dentro dessa caixa. Alguns cubos agora são "mistos": parte é água, parte é pedra.

Quando a água se move, o computador precisa calcular quanto líquido passa de um cubo para o outro. Nos métodos antigos, quando a água tentava passar por um cubo que estava quase cheio de pedra (um "cubo cortado"), o cálculo ficava confuso. Era como tentar encher um balde furado: o computador perdia um pouquinho de água aqui, ganhava um pouquinho ali, e, com o tempo, a gota parecia evaporar ou aparecer do nada. Isso quebrava a lei da conservação de massa (a água não pode sumir nem aparecer do nada).

A Solução (O "Balde Inteligente"):
Os autores criaram um novo método de "transporte" para esses cubos mistos.

• O Truque: Eles inventaram uma forma de "compensar" o espaço ocupado pela pedra. É como se, ao calcular quanto água passa, o computador dissesse: "Ah, essa parte do cubo é pedra, então a água só pode passar por aqui". Eles ajustaram a velocidade da água para preencher exatamente o espaço vazio, sem vazar.
• O Problema do Tempo: Em cubos muito pequenos (onde a pedra ocupa quase tudo), a física exigia que o computador calculasse em passos de tempo minúsculos, como se fosse uma filmagem em câmera super-lenta. Isso tornava a simulação extremamente lenta.
• A Inovação (Redistribuição): Eles desenvolveram uma técnica de "redistribuição". Imagine que você tem um balde que transbordou um pouquinho. Em vez de parar tudo para limpar o chão, você pega o excesso e o distribui rapidamente para os baldes vizinhos que ainda têm espaço. Isso permite que o computador use passos de tempo maiores (mais rápidos) sem perder a precisão ou a quantidade total de água.

2. O Problema do "Ângulo de Contato" (O Beijo da Gota)

A Situação:
Quando uma gota de água toca uma superfície, ela não forma um ângulo aleatório. Ela forma um ângulo específico (como 90 graus ou 30 graus), dependendo se a superfície é hidrofílica (gosta de água) ou hidrofóbica (odeia água).
Em superfícies planas e lisas, calcular esse ângulo é fácil. Mas em superfícies curvas, irregulares ou com texturas complexas (como a borda de um copo quebrado), os métodos antigos tentavam "esticar" a linha da gota como se fosse uma régua reta. Isso funcionava bem em superfícies lisas, mas em curvas, a régua errava o ângulo, fazendo a gota parecer quadrada ou triangular em vez de redonda.

A Solução (O "Modelador de Argila"):
Os autores substituíram a "régua reta" por um "modelador de argila" inteligente.

• A Técnica: Em vez de apenas olhar para a linha de contato, o novo método olha para a vizinhança inteira da gota e "encaixa" uma superfície curva (um paraboloide) que se ajusta perfeitamente à forma da gota e respeita o ângulo desejado na borda da pedra.
• A Analogia: É como se você estivesse moldando uma escultura de argila. Os métodos antigos tentavam cortar a argila com uma faca reta. O novo método usa as mãos para moldar a argila, garantindo que, onde a argila toca a mesa (a pedra), o ângulo seja exatamente o que você pediu, não importa se a mesa é curva ou torta.
• Histórico de Molhabilidade (Histerese): Eles também ensinaram o computador a lembrar do "histórico" da gota. Às vezes, uma gota "gruda" em um lugar e só se move quando a força é forte o suficiente (como quando você tenta empurrar uma caixa pesada no chão). O novo método simula esse comportamento de "grudar e soltar" com precisão.

Por que isso é importante?

Até agora, simular gotas de água em objetos complexos em 3D era como tentar dirigir um carro de brinquedo em um labirinto de espelhos: o computador perdia o rumo, a gota desaparecia ou ficava com formas estranhas.

Com este novo método:

1. A água não some: A conservação de massa é perfeita, mesmo em geometrias complexas.
2. A gota é realista: Ela se espalha, sobe em cones, cai em orifícios e se move em superfícies onduladas exatamente como faria na vida real.
3. É rápido: Não precisa mais de "câmera lenta" infinita para calcular os detalhes.

Resumo da Ópera:
Os autores criaram um "super-olho" matemático que consegue ver e calcular o movimento de líquidos em qualquer superfície, por mais estranha que seja, sem cometer erros de contagem e sem ficar lento demais. Isso abre portas para melhorar desde a impressão 3D de microgotas até a eficiência de células de combustível e a manipulação de fluidos em microchips.

Resumo técnico

Resumo Técnico: Método VOF 3D Conservador e de Interface Nítida para Dinâmica de Linha de Contato em Fronteiras Complexas

1. Problema e Motivação

A dinâmica da linha de contato (a interseção entre as fases líquido, gás e sólido) é fundamental em diversos processos de engenharia, como impressão jato de tinta, manipulação de filmes finos e eletrólise. A modelagem numérica precisa desses fenômenos em geometrias complexas (superfícies curvas, irregulares ou embutidas) apresenta desafios significativos que os métodos tradicionais em grades cartesianas não conseguem resolver adequadamente:

• Conservação de Massa em Células Mistas: Em métodos de Interface Nítida (Sharp Interface) com fronteiras embutidas (Embedded Boundary Method - EBM), as células cortadas pelo sólido (células mistas) contêm frações parciais de fluido e sólido. A advecção padrão da fração volumétrica (VOF) frequentemente falha em manter a conservação de massa local nessas células devido a erros no cálculo de fluxos através das faces cortadas.
• Restrições de Passo de Tempo (CFL): Células cortadas muito pequenas ("small cut cells") impõem restrições severas ao passo de tempo (CFL), tornando simulações de longo prazo computacionalmente proibitivas.
• Imposição de Ângulo de Contato: Métodos existentes para impor condições de ângulo de contato em superfícies curvas ou irregulares frequentemente utilizam extrapolações lineares simples, o que gera erros significativos na estimativa de curvatura e na simetria da linha de contato, especialmente para ângulos extremos (< 45° ou > 135°).
• Histerese: A modelagem da histerese do ângulo de contato (diferença entre ângulos de avanço e recuo) em 3D sobre geometrias complexas é um problema não resolvido de forma robusta.

2. Metodologia Proposta

O trabalho propõe um esquema numérico 3D acoplado dentro do solver de código aberto Basilisk, combinando o método Volume-of-Fluid (VOF) Geométrico com o método de Fronteiras Embutidas (EBM).

A. Esquema de Advecção Conservador para Células Mistas:

• Reconstrução de Interface: Utiliza um algoritmo de "flood-fill" (enchimento) estendido para 3D para reconstruir a interface líquido/gas dentro de poliedros convexos irregulares formados pela interseção do sólido com a célula.
• Correção de Fluxo: Introduz uma correção na velocidade de advecção ($u_n^*$) para compensar o volume bloqueado pelo sólido nas faces das células mistas, garantindo conservação de massa local estrita.
• Estratégia de Redistribuição (Redistribution): Para eliminar a restrição de passo de tempo imposta por células cortadas pequenas, o método adota uma estratégia de redistribuição. Se uma célula interfacial ficar "superlotada" ou "vazia" devido a um passo de tempo relaxado, o excesso ou déficit de volume é redistribuído conservativamente para os vizinhos a jusante. Isso remove a dependência do passo de tempo em relação ao tamanho da fração de fluido ($s_f/c_s$), mantendo a conservação global e local.

B. Imposição de Ângulo de Contato Baseada em Função de Altura (Height Function - HF):

• Ajuste de Parabolóide: Em vez de extrapolações lineares, o método utiliza um ajuste de parabolóide (superfície quadrática) para reconstruir a interface próximo à linha de contato.
• Pré-ajuste (Pre-fitting): Um estágio inicial identifica pontos de função de altura (HF) válidos na região de fluido ao redor da célula de contato, garantindo consistência geométrica.
• Ajuste Final com Restrição: Um sistema não linear é resolvido para ajustar o parabolóide final, forçando explicitamente que a normal da interface no ponto de contato satisfaça o ângulo de contato prescrito ($\theta$) em relação à normal do sólido.
• Função de Altura Vertical (VHF): Para casos degenerados onde a superfície sólida é quase paralela ao plano de acumulação (ângulos extremos), o método introduz uma abordagem de "função de altura vertical", permitindo a imposição correta do ângulo de contato mesmo em geometrias complexas onde métodos horizontais falham.

C. Modelagem de Histerese:

• Implementa um procedimento iterativo (método de bisseção) dentro da janela de histerese $[\theta_{rec}, \theta_{adv}]$. O ângulo de contato é ajustado dinamicamente para garantir que a linha de contato permaneça "pinned" (presa) até que a força capilar supere a barreira de histerese, mantendo a posição da linha de contato fixa em relação à geometria do sólido durante a iteração.

3. Principais Contribuições

1. Esquema de Advecção 3D Conservador e Livre de Restrições de CFL: Desenvolvimento de um método de advecção geométrica completo para células mistas em 3D que garante conservação de massa local estrita e elimina a necessidade de passos de tempo extremamente pequenos devido a células cortadas, através da redistribuição de volume.
2. Algoritmo 3D de Ângulo de Contato Robusto: Introdução de uma técnica baseada em ajuste de parabolóide com pré-seleção de pontos, capaz de impor condições de ângulo de contato com alta precisão em superfícies sólidas arbitrariamente complexas e curvas, superando as limitações de métodos lineares.
3. Integração de Histerese em 3D: Primeira implementação de um modelo de histerese de ângulo de contato acoplado a um esquema VOF conservador em 3D para superfícies complexas.
4. Validação Extensiva: O método é validado em uma série de testes desafiadores, demonstrando superioridade sobre abordagens de interface nítida existentes.

4. Resultados e Validação

O método foi testado em diversos cenários benchmarks:

• Esfera Concêntrica Rotativa: Demonstrou que o esquema corrigido atinge precisão de máquina para conservação de massa local, enquanto métodos tradicionais acumulam erros significativos.
• Gota em Superfície Plana e Esférica: O método de ajuste de parabolóide produziu raios de contato de equilíbrio em excelente acordo com a teoria, independentemente da orientação ou profundidade de embutimento do sólido. Diferentemente de métodos lineares, manteve a simetria circular da linha de contato mesmo em superfícies inclinadas.
• Movimento de Gota com Histerese: Simulações de gotas sob fluxo de cisalhamento reproduziram com precisão o comportamento de "pinning" e "depinning" (fixação e soltura) da linha de contato, comparável a benchmarks experimentais e numéricos anteriores.
• Geometrias Complexas:
  • Impacto em Orifício: Simulação de impacto de gota em placas com orifícios arredondados e afiados, capturando corretamente a pinagem e a histerese nas bordas.
  • Superfície Cônica e Rugosa: Simulação do movimento espontâneo de gotas em cones e superfícies cônicas com rugosidade senoidal, capturando a dinâmica de subida e a interação com as cristas da superfície.
  • Superfície em Malha: Demonstrou estabilidade e precisão ao simular o espalhamento de gotas em estruturas de malha cilíndrica com curvaturas abruptas.

5. Significância

Este trabalho representa um avanço significativo na simulação de escoamentos multifásicos 3D. Ao resolver simultaneamente os problemas de conservação de massa em células cortadas e a imposição precisa de condições de contorno de ângulo de contato em geometrias complexas, o método proposto torna-se uma ferramenta poderosa para pesquisas fundamentais e aplicações de engenharia.

A capacidade de lidar com qualquer geometria sólida 3D sem sacrificar a conservação de massa ou a precisão da interface, e de incluir histerese realisticamente, permite a investigação de fenômenos de molhagem em microestruturas, dispositivos microfluídicos e processos de impressão que antes eram difíceis de modelar com fidelidade. A eliminação da restrição de passo de tempo CFL para células pequenas é particularmente crucial para viabilizar simulações de longo prazo em geometrias complexas.