Parameter-robust well-posedness and discretisation for coupled Darcy--Forchheimer and advection-diffusion-reaction equations

Resumen Técnico: Bien posed de parámetros y discretización para ecuaciones acopladas de Darcy–Forchheimer y advección-difusión-reacción

Planteamiento del Problema
Este trabajo investiga el análisis matemático y la discretización numérica de un sistema totalmente acoplado que comprende las ecuaciones de Darcy–Forchheimer y la ecuación de advección-difusión-reacción (ADR). El sistema modela el flujo no lineal en medios porosos acoplado con el transporte de solutos, donde la velocidad del flujo $u$ y la presión $p$ satisfacen un balance de cantidad de movimiento no lineal con un término de arrastre de Forchheimer, y la concentración del soluto $\phi$ y el flujo total $\zeta$ satisfacen una ecuación de transporte. Un desafío crítico abordado es la dependencia de la estabilidad y la convergencia respecto a los parámetros del modelo, específicamente la conductividad hidráulica (permeabilidad) $\kappa$, el coeficiente de Forchheimer $F$, el coeficiente de difusividad $\alpha$ y el coeficiente de reacción $\varrho_0$. Los autores pretenden establecer la buena posición (well-posedness) y derivar estimaciones de error que sean robustas (independientes) con respecto a estos parámetros, los cuales pueden variar en varios órdenes de magnitud en aplicaciones prácticas como el modelado de yacimientos y la optimización de la energía geotérmica.

Metodología
El artículo emplea un enfoque de análisis funcional dentro del entorno de los espacios de Banach reflexivos, yendo más allá del marco estándar de los espacios de Hilbert para acomodar las no linealidades y los requisitos de escala específicos del problema.

1. Teoría Abstracta: Los autores adaptan la teoría reciente para problemas de punto de silla perturbados en espacios de Banach. Extienden la teoría de Brezzi–Braess para establecer cotas de estabilidad independientes de los parámetros. Esto implica la definición de espacios de Banach ponderados apropiados para la velocidad, la presión, el flujo y la concentración. Las normas en estos espacios se construyen como sumas e intersecciones de espacios ponderados por los parámetros del modelo (por ejemplo, $\kappa^{-1/2}$, $F^{1/3}$, $\alpha^{-1/2}$).
2. Análisis de Buena Posición: Utilizando un argumento de punto fijo, el sistema acoplado se desacopla en un subproblema de Darcy–Forchheimer y un subproblema de ADR. Los autores demuestran que los subproblemas individuales están bien condicionados y luego demuestran que la composición de sus operadores de solución forma un mapeo de contracción en una bola cerrada en el espacio de velocidades, siempre que los datos sean suficientemente pequeños. Esto establece la existencia y unicidad de la solución continua con constantes de estabilidad independientes de los parámetros.
3. Discretización por Elementos Finitos: El problema se discretiza utilizando métodos de elementos finitos mixtos. Los autores utilizan elementos de Raviart–Thomas para los espacios de velocidad y flujo, y espacios de Galerkin discontinuos (Sobolev rotos) para la presión y la concentración. Crucialmente, los espacios discretos están equipados con las mismas normas ponderadas por parámetros utilizadas en el análisis continuo.
4. Análisis de Error: Se derivan estimaciones de error a priori en estas normas ponderadas. El análisis se basa en condiciones de inf-sup discretas y un argumento de punto fijo discreto para mostrar que las tasas de convergencia son óptimas y robustas con respecto a $\kappa$, $F$, $\alpha$ y $\varrho_0$.
5. Precondicionamiento: Para abordar los desafíos computacionales de resolver los sistemas lineales no lineales y no simétricos resultantes, los autores diseñan precondicionadores de operador de bloque diagonal. Estos están inspirados en estructuras de complementos de Schur y marcos de precondicionamiento de operadores, adaptados a los pesos de parámetros específicos del sistema acoplado.

Contribuciones Clave

• Extensión a Espacios de Banach: El artículo extiende la teoría de Brezzi–Braess para problemas de punto de silla perturbados desde los espacios de Hilbert a los espacios de Banach reflexivos, abordando específicamente el contexto de robustez de parámetros.
• Estabilidad Robusta de Parámetros: Los autores demuestran que el sistema acoplado Darcy–Forchheimer/ADR es bien condicionado y estable con cotas independientes de la permeabilidad, el coeficiente de Forchheimer, la difusividad y los coeficientes de reacción. Esto se logra mediante la cuidadosa construcción de espacios y normas de Banach ponderadas.
• Análisis Discreto Robusto: El artículo proporciona un análisis discreto riguro el que muestra que la formulación de elementos finitos mixtos preserva estas propiedades de estabilidad robustas. Las estimaciones de error a priori derivadas son cuasi-óptimas e independientes de los parámetros del modelo y del tamaño de la malla $h$.
• Precondicionamiento de Operadores: Se propone el diseño de precondicionadores de bloque diagonal variables. Se demuestra que estos precondicionadores son robustos con respecto a los parámetros del modelo (excepto el coeficiente de difusividad en el régimen dominado por la advección) y al parámetro de discretización.

Resultados Numéricos
Los hallazgos teóricos se validan mediante una serie de experimentos numéricos en 2D y 3D:

• Convergencia: Pruebas con soluciones suaves manufacturadas y permeabilidad heterogénea confirman las tasas de convergencia óptimas ($O(h^{k+1})$) tanto en las normas estándar como en las ponderadas por parámetros. Los resultados demuestran que, mientras las normas estándar pueden mostrar sensibilidad a valores extremos de los parámetros, las normas ponderadas restauran el comportamiento de convergencia óptimo.
• Desempeño del Precondicionador: Las pruebas numéricas en el sistema linealizado muestran que los precondicionadores propuestos mantienen números de condición y conteos de iteración MINRES estables a través de amplias variaciones de $\kappa$, $F$, $\alpha$ y $\varrho_0$. Los autores señalan que la robustez se observa en el régimen dominado por la difusión; el desempeño en el régimen dominado por la advección (altos números de Péclet) no es totalmente robusto, lo cual es consistente con las limitaciones de la estabilización elegida.
• Aplicaciones Prácticas: El método se aplica a un problema de flujo de cavidad impulsada por tapa (lid-driven cavity) y a un problema de desplazamiento miscible transitorio (dedos viscosos en una celda de Hele–Shaw), demostrando la capacidad del método para manejar acoplamientos complejos de flujo-transporte y no linealidades.

Significación y Reivindicaciones
El artículo afirma que su principal significación radica en proporcionar un marco matemáticamente riguroso para el análisis y la simulación de flujos y transporte acoplados en medios porosos que permanece estable y preciso independientemente de la magnitud de los parámetros físicos. Al desplazar el análisis hacia los espacios de Banach y utilizar normas ponderadas por parámetros, los autores superan las limitaciones de los enfoques estándar de espacios de Hilbert, que a menudo no logran proporcionar una estabilidad uniforme para tales sistemas acoplados y no lineales. El trabajo ofrece una base teórica para el desarrollo de resolvedores numéricos eficientes y robustos para problemas donde los parámetros como la permeabilidad y las tasas de reacción varían significativamente, un escenario común en la ingeniería de subsuelo y el modelado ambiental. Los autores reconocen modestamente que, si bien el marco de punto fijo produce cotas teóricas robustas, la robustez numérica observada en el sistema acoplado sugiere que la pérdida de robustez analítica en ciertas normas ponderadas puede ser un artefacto del escalamiento o del propio marco de punto fijo, más que una limitación fundamental del método.