Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Resumen Técnico: Método de elementos virtuales de función de corriente – presión para el problema de la interfaz Stokes–Darcy

Planteamiento del Problema
El artículo aborda la simulación numérica del problema acoplado de la interfaz Stokes–Darcy, un desafío fundamental en la mecánica de fluidos que rige la interacción entre regiones de flujo libre (Stokes) y medios porosos (Darcy). El acoplamiento se logra mediante la continuidad de las velocidades normales (conservación de la masa), el equilibrio de tensiones normales y la condición de deslizamiento tangencial de Beavers–Joseph–Saffman. Aunque existen diversos métodos de Elementos Finitos (FE), de Galerkin discontinuo (dG) e Híbridos de Alto Orden (HHO), que a menudo dependen de pares de velocidad–presión que requieren estabilidad inf–sup, este trabajo se centra en una reformulación específica en dominios simplemente conexos en dos dimensiones. Aquí, la restricción de incompresibilidad se satisface de forma idéntica mediante la introducción de una función de corriente escalar, eliminando así la variable de presión de las ecuaciones de Stokes y reduciendo el sistema a un problema biarmónico de cuarto orden acoplado con una ecuación de presión de Darcy de segundo orden.

Metodología
Los autores proponen una novedosa discretización mediante el Método de Elementos Virtuales (VEM) para esta formulación de función de corriente–presión. La metodología aprovecha la capacidad del VEM para manejar mallas poligonales generales, lo cual es particularmente ventajoso para geometrías de interfaz complejas y mallas no coincidentes.

1. Formulación: El problema continuo se plantea en una forma primal-primal. El dominio de Stokes utiliza una función de corriente $\chi \in H^2$, mientras que el dominio de Darcy utiliza un potencial de presión $\phi \in H^1$. La formulación débil resulta en un problema de punto de silla perturbado.
2. Discretización:
   • Dominio de Stokes: Se emplea un espacio VEM de orden inferior con conformidad $C^1$. Esto implica polinomios cuadráticos ($P_2$) y requiere conformidad $H^2$, lograda mediante espacios virtuales mejorados con grados de libertad (DoFs) definidos por valores nodales y gradientes nodales.
   • Dominio de Darcy: Se utiliza un espacio VEM de orden inferior con conformidad $C^0$ para la presión, basado en polinomios lineales ($P_1$).
   • Acoplamiento: Las condiciones de interfaz se imponen naturalmente dentro de la formulación débil. El método acomoda mallas no coincidentes en la interfaz sin necesidad de remallar.
3. Implementación: El esquema se implementa utilizando la biblioteca de código abierto orientada a objetos vem++. Se utiliza una característica específica, la clase vem::vemMesh2dDofFilter, para gestionar los DoFs a través de diferentes subdominios, permitiendo la definición simultánea de diferentes espacios VEM en una sola malla y el filtrado de los mismos basados en marcadores de elementos.
4. Estabilización: Las formas bilineales discretas incluyen términos de estabilización ($S_S^K$ y $S_D^K$) para asegurar la coercitividad en el núcleo de los operadores de proyección, siguiendo las prácticas estándar de VEM (por ejemplo, receta diagonal y DOFI-DOFI).

Contribuciones Clave
El artículo describe cinco contribuciones primarias:

• Análisis Continuo: Desarrollo de un análisis continuo para el problema de la interfaz Stokes–Darcy en la formulación de función de corriente–presión, la cual evita inherentemente la necesidad de condiciones de estabilidad inf–sup.
• Esquema Discreto: Diseño de una discretización VEM de orden inferior $C^1$–$C^0$. Los autores destacan la eficiencia computacional de este esquema, señalando que involucra solo $3N_S + N_D + 4N_\Sigma$ grados de libertad (donde $N_S$, $N_D$ y $N_\Sigma$ son los recuentos de vértices interiores en el dominio de Stokes, el dominio de Darcy y en la interfaz, respectivamente).
• Análisis de Error: Derivación de estimaciones de error a priori para la propuesta de discretización, demostrando tasas de convergencia óptimas bajo supuestos de regularidad estándar.
• Implementación de Software: Integración de nueva funcionalidad en la biblioteca vem++ para facilitar la implementación directa de problemas de interfaz mediante el filtrado de mallas.
• Validación Numérica: Ejecución de experimentos numéricos para confirmar los hallazgos teóricos.

Resultados
Los experimentos numéricos validan las tasas de convergencia teórica y la robustez del método:

• Soluciones Manufacturadas: Dos experimentos utilizando soluciones manufacturadas (uno con valores de interfaz cero, otro con valores no nulos) en varios tipos de mallas (cuadriláteros, no convexas, Voronoi y Voronoi perturbadas) demuestran un decaimiento lineal asintótico del error ($O(h)$), consistente con las predicciones teóricas para el esquema de orden inferior.
• Filtración de Calle Sin Salida (Dead-End): Una simulación de un filtro de calle sin salida en forma de cuarto de anillo recupera con éxito los resultados de estudios previos, demostrando la capacidad del método para manejar interfaces curvas y condiciones de contorno mixtas.
• Órganos Bioartificiales: Una simulación del flujo sanguíneo en un órgano bioartificial (basado en encapsulación) con constantes de permeabilidad variables ilustra la capacidad del método para modelar flujos fisiológicos complejos y la influencia de la permeabilidad en el flujo a través de la interfaz.

Significancia y Reivindicaciones
El artículo posiciona su trabajo como un puente para cerrar la brecha en la literatura respecto a la extensión de las formulaciones de función de corriente a problemas acoplados de Stokes–Darcy utilizando Elementos Virtuales. Los autores afirman que la estructura propuesta ofrece una "estrategia completamente libre de inf–sup", simplificando el diseño de espacios discretos en comparación con los enfoques tradicionales de velocidad–presión. Se destaca la capacidad de manejar geometrías de interfaz irregulares y mallas no coincidentes sin remallar como una ventaja clave del marco de trabajo VEM.

Los autores mantienen la modestia respecto a desafíos futuros, señalando explícitamente que la extensión de esta formulación de función de corriente a dominios múltiplemente conexos (donde la función de corriente no está unívocamente determinada por condiciones de contorno estándar) y el escalado de este marco a tres dimensiones (que requeriría potenciales vectoriales y espacios $C^1$ de alto orden para funciones con valores vectoriales) siguen siendo desafíos abiertos. Además, la integración de técnicas sin estabilización y la derivación de estimadores de error a posteriori robustos se identifican como direcciones para investigaciones futuras.