Deformable bodies in a 3-dimensional viscous flow: Vorticity-Stream vector formulation

Resumo Técnico: Corpos Deformáveis em Fluxo Viscoso 3D via Formulação de Vetor de Vorticidade-Stream

Definição do Problema
Simular a interação entre fluxos viscosos incompressíveis tridimensionais (3D) e obstáculos deformáveis ou elásticos (como vesículas, gotículas e células) apresenta desafios significativos. Métodos existentes, incluindo o Método de Lattice Boltzmann (LBM), técnicas de integral de contorno e abordagens baseadas em partículas como o Smoothed Dissipative Particle Dynamics (SDPD), frequentemente sofrem com alta complexidade computacional, implementações numéricas intrincadas ou dificuldades em lidar com mecânicas interfaciais complexas. Embora simulações 2D utilizando formulações de função de corrente-vorticidade acopladas a modelos de campo de fase tenham se mostrado robustas, faltava uma estrutura comparável, leve e classicamente fundamentada para interfaces deformáveis 3D em baixos números de Reynolds. Especificamente, há uma lacuna na literatura quanto a formulações baseadas em vorticidade para interfaces 3D governadas por energias de Canham–Helfrich (dobramento de membrana) ou Cahn–Hilliard (tensão superficial).

Metodologia
Os autores introduzem uma nova formulação numérica baseada na abordagem de vetor de vorticidade-stream, especificamente adaptada para fluxos de baixo número de Reynolds. O método generaliza a função de corrente 2D para um potencial vetorial 3D ($\vec{\psi}$), onde o campo de velocidade é definido como $\vec{v} = \nabla \times \vec{\psi}$.

O núcleo da metodologia envolve a resolução de um sistema acoplado de equações para três observáveis primários:

1. Potencial Químico da Interface ($\mu$): Derivado da derivada funcional de uma energia livre escolhida ($F$) em relação ao parâmetro de ordem do campo de fase ($\phi$). Os autores demonstram a estrutura utilizando tanto a energia de Canham–Helfrich (para dobramento de membrana) quanto a energia de Cahn–Hilliard (para tensão superficial).
2. Vorticidade do Fluido ($\vec{\omega}$): Definida como $\vec{\omega} = \nabla \times \vec{v}$.
3. Vetor de Stream ($\vec{\psi}$): O potencial vetorial que satisfaz a condição de gauge $\nabla \cdot \vec{\psi} = 0$.

As equações governantes reduzem-se a um conjunto de problemas de Poisson acoplados:

• A evolução do campo de fase $\phi$ é governada por uma equação de advecção-difusão envolvendo o potencial químico e o campo de velocidade.
• A equação da vorticidade é derivada do rotacional da equação de Stokes (ou Navier-Stokes para extensões inerciais), resultando em uma equação de Poisson onde o termo de fonte depende do produto vetorial entre o gradiente de fase e o gradiente do potencial químico: $\nabla^2 \vec{\omega} = \frac{1}{\eta} (\nabla \phi) \times (\nabla \mu)$.
• O vetor de stream é obtido resolvendo $\nabla^2 \vec{\psi} = -\vec{\omega}$.

Implementação Numérica
O sistema é resolvido utilizando um método de diferenças finitas em uma grade cúbica fixa e não ajustada (unfitted). A interface é tratada como uma região difusiva via abordagem de campo de fase, eliminando a necessidade de deformação complexa de malha ou técnicas de contorno imerso. A evolução temporal emprega integração de Euler progressivo para o campo de fase e um solver iterativo para os problemas de Poisson, utilizando um esquema estagiado (segregado sequencial). O código é implementado em Python e está disponível em um repositório de código aberto.

Principais Contribuições

• Simplificação de Solvers 3D: A formulação transforma o solver de fluido em um conjunto de problemas de Poisson acoplados, evitando o cálculo direto da pressão e garantindo a incompressibilidade estrita por construção.
• Estrutura Unificada para Energias Interfaciais: O método integra perfeitamente diferentes modelos de energia livre (Canham–Helfrich para dobramento e Cahn–Hilliard para tensão superficial) dentro da mesma estrutura numérica, permitindo a customização fácil da mecânica da membrana.
• Observáveis Diagnósticos: Os autores propõem novos observáveis baseados nas flutuações quadráticas dos campos de vorticidade e vetor de stream em relação ao fluxo não perturbado. Estas métricas mostraram-se altamente sensíveis a deformações de forma e transições morfológicas.

Resultados
O método foi validado através de simulações de vesículas e gotículas únicas em fluxos de Poiseuille e Couette Newtonianos:

• Fluxo de Poiseuille: As simulações recuperaram com sucesso formas axissimétricas canônicas para glóbulos vermelhos (RBCs), incluindo morfologias de "paraquedas" (parachute) e "saliência" (slipper), consistentes com a literatura. Os observáveis diagnósticos rastrearam eficazmente a dinâmica de relaxação, mostrando picos distintos correspondentes a mudanças de forma (ex: o desaparecimento de um sulco) e estabilizando-se quando o objeto atingiu um estado estacionário.
• Fluxo de Couette: A estrutura capturou a migração lateral, um fenômeno de sustentação (lift) não inercial onde vesículas migram perpendicularmente à direção do fluxo. As simulações demonstraram que a posição de equilíbrio e a dinâmica de migração dependem sistematicamente do contraste de viscosidade ($\zeta$) entre os fluidos interno e externo. Contrastes de viscosidade mais elevados levaram a um deslocamento mais pronunciado em direção à parede móvel.
• Análise de Tensão: O método calculou com sucesso as distribuições de tensão de cisalhamento no plano, revelando que taxas de cisalhamento mais altas induzem maior heterogeneidade de tensão e podem levar a instabilidades morfológicas.

Significância e Alegações
O artigo afirma que esta formulação oferece uma alternativa leve, flexível e computacionalmente eficiente para os solvers mesoscópicos ou de alta complexidade existentes para hidrodinâmica 3D de baixo número de Reynolds. Ao basear-se em princípios clássicos de mecânica dos fluidos em vez de teoria cinética (como no LBM), o método evita a necessidade de provas de consistência com a mecânica do contínuo.

Os autores enfatizam que, embora o trabalho atual foque na dinâmica de corpo único em fluidos Newtonianos, a estrutura foi projetada para ser extensível. Ela pode acomodar naturalmente forças de corpo, efeitos inerciais, meios viscoelásticos e interações entre múltiplos corpos. A capacidade de capturar a dinâmica de cisalhamento no plano 3D — uma característica frequentemente perdida em aproximações 2D — posiciona este método como uma ferramenta valiosa para o estudo da biomecânica celular e física da matéria mole. A natureza de código aberto do trabalho visa fomentar o benchmarking e o desenvolvimento adicional nesta área emergente.