Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems

Resumen Técnico: Estimaciones de Error de Cuadratura en Mallas No Coincidentes en un Marco de Dominio Ficticio para Problemas de Interacción Fluido-Estructura

Planteamiento del Problema
El artículo aborda la discretización numérica de problemas de interacción fluido-estructura (FSI) utilizando una formulación de dominio ficticio basada en multiplicadores de Lagrange distribuidos. En este marco, el fluido y el sólido se discretizan utilizando dos mallas completamente independientes y no coincidentes. El sólido se describe en un marco lagrangiano sobre un dominio de referencia $B$, mapeado a su posición actual dentro de un dominio euleriano fijo de fluido $\Omega$. El acoplamiento cinemático entre la velocidad del fluido y la deformación del sólido se impone mediante un multiplicador de Lagrange $\lambda$ definido en el dominio de referencia del sólido.

Un desafío computacional crítico surge al ensamblar la matriz de acoplamiento (interfaz). Esta matriz requiere integrar el producto de las funciones de base del fluido (compuestas con el mapeo del sólido) y las funciones de base del sólido sobre los elementos del sólido. Dado que estas funciones están definidas en mallas no coincidentes, la integración exacta requiere calcular la intersección geométrica entre las mallas del fluido y del sólido mapeado para construir una regla de cuadratura compuesta. Este proceso es computacionalmente costoso, particularmente porque el término de acoplamiento debe ensamblarse en cada paso de tiempo. Alternativamente, un enfoque aproximado integra directamente sobre los elementos del sólido utilizando reglas de cuadratura estándar sin intersección de mallas, lo que introduce un error de cuadratura. El artículo investiga la validez teórica y las estimaciones de error asociadas con esta estrategia de integración aproximada.

Metodología
Los autores analizan un modelo de FSI basado en Stokes en estado estacionario, aunque señalan que sus resultados son representativos de casos no lineales más generales. La formulación variacional consiste en hallar la velocidad del fluido $u$, la presión $p$, el mapa de deformación del sólido $X$ y el multiplicador de Lagrange $\lambda$ que satisfacen un sistema de punto de silla.

El estudio se centra en dos elecciones específicas para la forma bilineal $c(\cdot, \cdot)$ utilizada para definir el acoplamiento:

1. $c_0$ (Emparejamiento dual): Basada en el producto interno $L^2(B)$, donde el espacio del multiplicador es el dual de $H^1(B)$.
2. $c_1$ (Producto escalar): Basada en el producto interno $H^1(B)$, donde el espacio del multiplicador es $H^1(B)$ mismo.

La metodología procede a través de los siguientes pasos:

• Discretización: Se introducen espacios de elementos finitos para todas las variables. El fluido utiliza el elemento de Bercovier–Pironneau ($P1$-iso-$P2/P1$), mientras que el mapeo del sólido y el multiplicador utilizan elementos continuos por trozos lineales.
• Integración Aproximada: Los autores definen formas bilineales discretas donde los términos de acoplamiento se computan utilizando reglas de cuadratura estándar en la malla del sólido, ignorando el hecho de que las funciones de base del fluido compuestas son polinomios por trozos con discontinuidades respecto a la malla del sólido.
• Análisis de Error Abstracto: Siguiendo el marco de las condiciones de inf-sup, los autores derivan una estimación de error general (Teorema 2) que acota la diferencia entre la solución discreta exacta y la solución obtenida con la integración inexacta. Este límite depende de los funcionales de error de cuadratura.
• Estimación del Error de Cuadratura: La contribución técnica central implica demostrar estimaciones específicas para el error de cuadratura de los términos de acoplamiento. Los autores utilizan el hecho de que la función compuesta $v_h(X)$ es continua pero lineal por trozos en los elementos del sólido, lo que conduce a propiedades de regularidad específicas en espacios de Sobolev fraccionarios ($H^{1+s}$). Emplean la teoría de interpolación e desigualdades inversas para acotar el error.
• Prueba de Estabilidad: Los autores demuestran que el problema discreto sigue siendo bien planteado (satisface las condiciones de inf-sup discretas) incluso cuando el término de acoplamiento se construye de forma aproximada, siempre que se cumplan ciertas condiciones sobre los tamaños de malla.

Contribuciones Clave y Resultados

1. Bien Planteamiento con Acoplamiento Aproximado: El artículo demuestra que el problema discreto es estable y bien planteado incluso cuando la matriz de acoplamiento se ensambla mediante cuadratura aproximada, sin requerir la intersección de mallas.
2. Estimaciones de Error de Cuadratura:
   • Para el acoplamiento basado en $L^2$ ($c_0$), el error de cuadratura está acotado por $O(h_B^{3/2} |\log h_B|)$. Los autores muestran que el método aproximado produce tasas de convergencia óptimas siempre que $h_B \to 0$.
   • Para el acoplamiento basado en $H^1$ ($c_1$), el análisis de error es más complejo debido al término del gradiente. El límite del error incluye términos dependientes de la relación de los tamaños de malla. Específicamente, la convergencia óptima requiere no solo que $h_B \to 0$, sino también que la relación $h_B/h_\Omega \to 0$ (donde $h_\Omega$ es el tamaño de la malla del fluido). Si $h_B/h_\Omega$ se mantiene constante, el método puede perder optimalidad.
3. Desigualdades Inversas: Los autores derivan y utilizan desigualdades inversas específicas para manejar las normas del multiplicador de Lagrange en las estimaciones de error, particularmente para el caso de la norma dual.
4. Validación Numérica: Se presentan cuatro pruebas numéricas para validar las estimaciones teóricas:
   • Prueba 1 y 2: Comparan las tasas de convergencia para $c_0$ y $c_1$ bajo diferentes estrategias de refinamiento de malla (razón constante frente a razón decreciente). Los resultados confirman que $c_0$ es robusto con integración aproximada, mientras que $c_1$ requiere que $h_B/h_\Omega \to 0$ para mantener la precisión.
   • Prueba 3: Utiliza una solución con un salto en el gradiente de la velocidad a través de la interfaz (típico de FSI), confirmando la tasa de convergencia subóptima esperada de $O(h^{1/2})$ debido a singularidades de la solución.
   • Prueba 4: Investiga elementos de orden superior ($P2/P1$), mostrando que mientras $c_0$ permanece robusto, $c_1$ aún requiere integración exacta o razones de malla específicas para alcanzar tasas óptimas.

Significado y Reivindicaciones
El artículo sostiene que proporciona una justificación teórica rigurosa para el uso de la integración aproximada en métodos de FSI de dominio ficticio, una práctica adoptada frecuentemente por eficiencia computacional. Los autores demuestran que:

• El enfoque "inexacto" no compromete la estabilidad del esquema numérico.
• El acoplamiento basado en $L^2$ ($c_0$) es particularmente adecuado para la integración aproximada, ya que no requiere que las mallas del fluido y del sólido se refinen en una relación específica para mantener la precisión.
• El acoplamiento basado en $H^1$ ($c_1$), aunque potencialmente ofrece diferentes propiedades de estabilidad, impone requisitos de malla más estrictos ($h_B/h_\Omega \to 0$) para evitar que los errores de cuadratura dominen el error de la solución.

Los autores concluyen que sus resultados teóricos confirman las observaciones experimentales de su trabajo previo [12, 13], proporcionando una base matemática sólida para elegir entre las técnicas de ensamblaje exacta (con intersección de mallas) y aproximada (sin intersección) basándose en la forma de acoplamiento específica y la configuración de la malla.