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

Resumen Técnico: Cuerpos Deformables en Flujo Viscoso 3D mediante la Formulación de Vector de Vorticidad-Función de Corriente

Planteamiento del Problema
La simulación de la interacción entre flujos viscosos incompresibles tridimensionales (3D) y obstáculos deformables o elásticos (como vesículas, gotas y células) presenta desafíos significativos. Los métodos existentes, incluyendo los Métodos de Lattice Boltzmann (LBM), las técnicas de integral de contorno y los enfoques basados en partículas como la Dinámica de Partículas Disipativas Suavizadas (SDPD), suelen sufrir de una alta complejidad computacional, implementaciones numéricas intrincadas o dificultades para manejar la mecánica interfacial compleja. Si bien las simulaciones en 2D utilizando formulaciones de función de corriente-vorticidad acopladas con modelos de campo de fase han demostrado ser robustas, ha faltado un marco comparable, ligero y basado en principios clásicos para interfaces deformables 3D a bajos números de Reynolds. Específicamente, existe un vacío en la literatura respecto a las formulaciones basadas en la vorticidad para interfaces 3D gobernadas por energías de Canham–Helfrich (flexión de membrana) o Cahn–Hilliard (tensión superficial).

Metodología
Los autores introducen una novedosa formulación numérica basada en el enfoque de vector de vorticidad-función de corriente, específicamente adaptada para flujos de bajo número de Reynolds. El método generaliza la función de corriente 2D a un potencial vectorial 3D ($\vec{\psi}$), donde el campo de velocidad se define como $\vec{v} = \nabla \times \vec{\psi}$.

El núcleo de la metodología consiste en resolver un sistema acoplado de ecuaciones para tres observables primarios:

1. Potencial Químico de la Interfaz ($\mu$): Derivado de la derivada funcional de una energía libre elegida ($F$) con respecto al parámetro de campo de fase ($\phi$). Los autores demuestran el marco utilizando tanto la energía de Canham–Helfrich (para la flexión de la membrana) como la energía de Cahn–Hilliard (para la tensión superficial).
2. Vorticidad del Fluido ($\vec{\omega}$): Definida como $\vec{\omega} = \nabla \times \vec{v}$.
3. Vector de Corriente ($\vec{\psi}$): El potencial vectorial que satisface la condición de calibre $\nabla \cdot \vec{\psi} = 0$.

Las ecuaciones gobernantes se reducen a un conjunto de problemas de Poisson acoplados:

• La evolución del campo de fase $\phi$ está gobernada por una ecuación de advección-difusión que involucra el potencial químico y el campo de velocidad.
• La ecuación de vorticidad se deriva del rotacional de la ecuación de Stokes (o de Navier-Stokes para extensiones inerciales), resultando en una ecuación de Poisson donde el término fuente depende del producto cruz del gradiente de fase y el gradiente del potencial químico: $\nabla^2 \vec{\omega} = \frac{1}{\eta} (\nabla \phi) \times (\nabla \mu)$.
• El vector de corriente se obtiene resolviendo $\nabla^2 \vec{\psi} = -\vec{\omega}$.

Implementación Numérica
El sistema se resuelve utilizando un método de diferencias finitas en una rejilla cúbica fija y no ajustada (unfitted). La interfaz se trata como una región difusiva mediante el enfoque de campo de fase, eliminando la necesidad de deformación compleja de la malla o técnicas de frontera inmersa. La evolución temporal emplea una integración de Euler hacia adelante para el campo de fase y un solver iterativo para las ecuaciones de Poisson, utilizando un esquema escalonado (segregado secuencial). El código está implementado en Python y está disponible en un repositorio de código abierto.

Contribuciones Clave

• Simplificación de Solvers 3D: La formulación transforma el solver de fluido en un conjunto de problemas de Poisson acoplados, evitando el cálculo directo de la presión y asegurando la incompresibilidad estricta por construcción.
• Marco Unificado para Energías Interfaciales: El método integra de forma fluida diferentes modelos de energía libre (Canham–Helfrich para flexión y Cahn–Hilliard para tensión superficial) dentro de la misma estructura numérica, permitiendo una personalización sencilla de la mecánica de la membrana.
• Observables Diagnósticos: Los autores proponen nuevos observables basados en las fluctuaciones al cuadrado de los campos de vorticidad y de vector de corriente respecto al flujo no perturbado. Se demuestra que estas métricas son altamente sensibles a las deformaciones de forma y a las transiciones morfológicas.

Resultados
El método fue validado mediante simulaciones de vesículas y gotas individuales en flujos de Poiseuille y Couette Newtonianos:

• Flujo de Poiseuille: Las simulaciones recuperaron con éxito las formas aximétricas canónicas para los glóbulos rojos (RBCs), incluyendo las morfologías de "paracaídas" y "slipper", consistentes con la literatura. Los observables diagnósticos rastrearon eficazmente la dinámica de relajación, mostrando picos distintos correspondientes a cambios de forma (por ejemplo, la desaparición de un hoyuelo/dimple) y estabilizándose una vez que el objeto alcanzó un estado estacionario.
• Flujo de Couette: El marco capturó la migración lateral, un fenómeno de elevación no inercial donde las vesículas migran perpendicularmente a la dirección del flujo. Las simulaciones demostraron que la posición de equilibrio y la dinámica de migración dependen sistemáticamente del contraste de viscosidad ($\zeta$) entre los fluidos interno y externo. Los contrastes de viscosidad más altos provocaron un desplazamiento más pronunciado hacia la pared móvil.
• Análisis de Esfuerzo: El método calculó con éxito las distribuciones de esfuerzo cortante en el plano, revelando que tasas de cizallamiento más altas inducen una mayor heterogeneidad de esfuerzo y pueden conducir a inestabilidades morfológicas.

Significancia y Reivindicaciones
El artículo afirma que esta formulación ofrece una alternativa ligera, flexible y computacionalmente eficiente a los solvers existentes de escala mesoscópica o de alta complejidad para la hidrodinámica de bajo número de Reynolds en 3D. Al apoyarse en los principios de la mecánica de fluidos clásica en lugar de la teoría cinética (como en LBM), el método evita la necesidad de pruebas de consistencia con la mecánica de medios continuos.

Los autores enfatizan que, aunque el trabajo actual se centra en la dinámica de cuerpos únicos en fluidos Newtonianos, el marco está diseñado para ser extensible. Puede acomodar naturalmente fuerzas corporales, efectos inerciales, medios viscoelásticos e interacciones de múltiples cuerpos. La capacidad de capturar la dinámica de cizalladura en el plano 3D —una característica que a menudo se pierde en las aproximaciones 2D— posiciona a este método como una herramienta valiosa para el estudio de la biomecánica celular y la física de la materia blanda en flujos confinados. La naturaleza de código abierto del programa tiene como objetivo fomentar la comparación (benchmarking) y el desarrollo posterior en esta área emergente.