Accelerating droplet-laden Stokes flow simulations with hierarchical surrogate modeling

Resumen Técnico: Aceleración de simulaciones de flujo de Stokes con gotas mediante modelado sustituto jerárquico

Planteamiento del problema
El artículo aborda los desafíos computacionales inherentes a la simulación de flujos de Stokes multifásicos incompresibles que involucran un gran número ($M$) de gotas líquidas inmiscibles suspendidas en un fluido portador. Las simulaciones precisas de alta fidelidad de tales sistemas son computacionalmente exigentes, particularmente cuando las gotas están altamente concentradas, se acercan a los límites o interactúan estrechamente. Las técnicas de discretización estándar, como los métodos de elementos finitos (que requieren remallado repetido) o los métodos de integral de contorno (que requieren una cuadratura especializada para integrales casi singulares), tienen dificultades con los espacios intersticiales estrechos y las interacciones cercanas a las fronteras. Aunque los métodos basados en partículas evitan el remallado, a menudo requieren resoluciones finas para lograr precisión y carecen de la exactitud de los métodos de alto orden. Los autores pretenden desarrollar una estrategia de modelado sustituto que permita la predicción eficiente de la evolución de las gotas en configuraciones densamente empaquetadas y variables en el tiempo sin el costo prohibitivo de las simulaciones de multifluidos totalmente resueltas.

Metodología
El enfoque propuesto es un marco de modelado sustituto jerárquico, de múltiples fidelidades e informado por la física. Se basa en una estrategia iterativa de corrección de errores en lugar de tratar al sustituto como una "caja negra". La metodología procede de la siguiente manera:

1. Aproximación de la fidelidad más baja: El proceso comienza con una solución de Stokes base ($u_0, p_0$) donde las gotas se tratan como trazadores pasivos. Esto asume una viscosidad "nominal" homogénea en todo el dominio, ignorando los contrastes de viscosidad y la tensión superficial.
2. Derivación de la ecuación de error: Los autores derivan una ecuación de derivadas parciales (PDE) que gobierna el campo de error ($v, q$) entre la solución de menor fidelidad y la aproximación real. Esta ecuación de error tiene en cuenta el salto en la viscosidad y las fuerzas de tensión superficial a través de las interfaces de las gotas.
3. Corrección jerárquica mediante iteración: En lugar de resolver la ecuación de error global directamente (lo que retiene la complejidad computacional original), el método emplea un esquema de iteración alternante:
   • Corrección de gotas: La ecuación de error global se descompone en $M$ problemas locales independientes, cada uno de los cuales representa la perturbación del flujo causada por una sola gota en el espacio libre. Estos problemas locales se resuelven de forma aproximada.
   • Corrección de frontera: La superposición de las correcciones locales viola las condiciones de contorno globales. Un paso de "corrección de frontera" posterior resuelve un problema de Stokes global (sin gotas) para restaurar la consistencia con los límites del dominio.
   • Iteración: Estos dos pasos (correcciones de gota y de frontera) se repiten hasta que se cumple un criterio de convergencia basado en los residuales de la frontera.
4. Eficiencia Offline-Online: Para gotas geométricamente similares (específicamente circulares en las pruebas presentadas), el método aprovecha la autosimilitud. Los problemas de una sola gota se mapean a un disco unitario de referencia. Se utiliza un enfoque de subespacio lineal para aproximar el operador de solución de estos problemas de una sola gota utilizando un pequeño conjunto de "soluciones fundamentales" precomputadas correspondientes a campos de flujo polinómicos. Esto permite que la evaluación en línea de las correcciones de las gotas se realice mediante combinaciones lineales simples, reduciendo significativamente el costo computacional.

Contribuciones clave

• Corrección de error basada en la física: El artículo introduce una estrategia iterativa de múltiples fidelidades que recupera sistemáticamente soluciones de alta fidelidad mediante la resolución de ecuaciones de error de PDE, evitando la naturaleza de "caja negra" de muchos sustitutos de aprendizaje automático.
• Reemplazo del refinamiento por iteración: El método evita la alta complejidad de resolver ecuaciones de error globales en domios estrechamente empaquetados mediante la construcción de correcciones locales a través de la superposición de campos de una sola gota. Esto reduce los requisitos de memoria y el costo computacional.
• Escalabilidad para muchas gotas: El marco está diseñado para manejar sistemas con hasta $10^4$ gotas. Al desacoplar el problema de frontera global de las interacciones locales de las gotas y utilizar soluciones fundamentales precomputadas, el método logra aceleraciones sustanciales en comparación con las simulaciones totalmente resueltas.
• Tratamiento de la tensión superficial: Para gotas circulares, los efectos de la tensión superficial se contabilizan exactamente mediante saltos de presión, lo que permite que la iteración de corrección se centre en los contrastes de viscosidad y las interacciones de flujo.

Resultados
Experimentos numéricos en dos dimensiones validan la precisión y eficiencia del sustituto propuesto:

• Precisión: El sustituto logra errores relativos inferiores a $\sim 0.1\%$ en comparación con las simulaciones CutFEM de alta fidelidad totalmente resueltas. Las pruebas de convergencia muestran una convergencia exponencial con respecto al grado polinómico de las soluciones fundamentales y una convergencia robusta de las iteraciones de punto fijo, incluso para gotas en proximidad cercana a las fronteras o entre sí.
• Eficiencia: En una simulación en el dominio del tiempo con 50 gotas, el sustituto logró una aceleración de 80 veces sobre el código de referencia CutFEM manteniendo una precisión aceptable. Una versión de "menor fidelidad" de trazador pasivo logró una aceleración de $3 \times 10^5$ con un error un orden de magnitud mayor.
• Escalabilidad: Para problemas estáticos con hasta $10^4$ gotas, los costos de entrenamiento y evaluación del sustituto escalan favorablemente. Aunque el entrenamiento exhibe un escalamiento cuadrático ($O(M^2)$) debido a las interacciones por pares, el costo de evaluación es lineal ($O(M)$) cuando el número de puntos de evaluación es moderado. El método sigue siendo competitivo con CutFEM para $M < 5000$ al incluir los costos de evaluación, y supera significativamente a CutFEM cuando solo se requiere un número limitado de evaluaciones (por ejemplo, para la advección de marcadores en simulaciones de dominio temporal).
• Robustez: El esquema iterativo demostró robustez incluso en configuraciones no físicas con gotas superpuestas, convergiendo donde los métodos estándar podrían fallar o requerir un refinamiento de malla excesivo.

Significancia y afirmaciones
Los autores afirman que este marco de modelado sustituto jerárquico ofrece una vía práctica para simular flujos de Stokes multifásicos complejos en estaciones de trabajo estándar, una tarea que típicamente requiere recursos de supercomputación para métodos de alta fidelidad. El método es particularmente significativo para escenarios que involucran muchas gotas y pequeñas separaciones interfaciales, donde las técnicas de discretización estándar enfrentan cuellos de botella computacionales severos. El artículo enfatiza que el enfoque captura sistemáticamente las interacciones gota-flujo, gota-frontera y gota-gota. Aunque la implementación actual asume gotas circulares, los autores señalan que el marco es suficientemente rico como para motivar futuras extensiones a gotas deformables y potencialmente a flujos de Navier-Stokes. El trabajo demuestra que, al combinar conocimientos físicos (PDE de error) con álgebra lineal numérica eficiente (aproximaciones de subespacio), se puede lograr un equilibrio entre precisión y tractabilidad computacional previamente inalcanzable para sistemas de gotas a gran escala.