Distributed optimal control problems governed by poroelasticity equations

Resumen Técnico: Problemas de Control Óptimo Distribuido Gobernados por Ecuaciones de Poroelasticidad

Planteamiento del Problema
Este artículo aborda la aproximación numérica de problemas de control óptimo distribuido gobernados por el modelo de consolidación de Biot en poroelasticidad. El sistema físico describe la interacción entre el flujo de fluido intersticial y la deformación de un medio poroso, un modelo crítico para aplicaciones en geofísica, mecánica de tejidos biológicos y secuestro de carbono. El problema de control específico busca minimizar un funcional de costo que mide la desviación del desplazamiento sólido ($u$) y la presión del fluido ($p$) respecto a valores objetivo deseados ($u_C, p_C$), penalizando al mismo tiempo la magnitud de la variable de control. El control se identifica como el término de fuente/sumidero de fluido ($m$) distribuido dentro del dominio. El estudio se centra en la consolidación cuasi-estática donde se descuidan los efectos de aceleración, y aborda específicamente los desafíos matemáticos y numéricos asociados con el "bloqueo poroelástico", que ocurre cuando el coeficiente de almacenamiento específico ($s_0$) es cercano a cero, volviendo al fluido y al medio elástico casi incompresibles.

Metodología
Los autores proponen y analizan una nueva formulación simétrica de dos campos donde las incógnitas primarias son el desplazamiento sólido y la presión del fluido. La metodología procede a través de las siguientes etapas:

1. Formulación Continua y Bien Posebilidad:

   • Los autores reformulan las ecuaciones gobernantes tomando la derivada temporal del equilibrio de cantidad de movimiento, convirtiendo el sistema en un conjunto de ecuaciones diferenciales.
   • Definen los espacios de funciones apropiados ($V$ para el desplazamiento, $Q$ para la presión) y establecen la bien posebilidad de la solución débil utilizando la teoría de problemas de punto de silla y la condición de Ladyzhenskaya–Babuška–Brezzi (LBB).
   • La existencia y unicidad del control óptimo se demuestran mediante el método directo en el cálculo de variaciones, aprovechando la convexidad del funcional de costo y la linealidad del operador de estado.
   • Se derivan las condiciones de optimalidad necesarias y suficientes de primer orden mediante un enfoque de estado adjunto, lo que conduce a una desigualdad variacional para el control.

2. Esquema Discreto:

   • Discretización Espacial: El dominio se discretiza utilizando elementos finitos ($V_h \times Q_h$) que satisfacen la condición inf-sup (por ejemplo, elementos MINI o elementos de Taylor-Hood).
   • Discretización Temporal: Se emplea un esquema de Euler hacia atrás (específicamente un método dG(0) en el tiempo), lo que resulta en un sistema totalmente discreto.
   • Discretización del Control: El artículo utiliza un enfoque de discretización variacional (propuesto por Hinze) para la variable de control. Esto significa que el espacio de control no se discretiza explícitamente; en su lugar, el control discreto se obtiene proyectando el estado adjunto discreto sobre el conjunto de control admisible. Esto evita los errores de aproximación que típicamente se introducen al discretizar el espacio de control.

3. Análisis de Error:

   • Los autores derivan estimaciones de error a priori para el esquema totalmente discreto.
   • El análisis separa el error total en errores de discretización (para un control fijo) y errores de perturbación debidos a la diferencia entre los controles continuos y discretos.
   • Se establecen estimaciones de estabilidad tanto para el sistema de estado hacia adelante como para el sistema adjunto hacia atrás para acotar estos errores.

Contribuciones Clave
El artículo realiza tres contribuciones principales a la literatura:

• Nueva Formulación Simétrica: Introduce una formulación simétrica de dos campos (desplazamiento y presión) diseñada específicamente para problemas de control óptimo dentro del modelo de Biot, distinta de las formulaciones previas de cuatro campos o mixtas.
• Existencia y Unicidad: Demuestra rigurosamente la bien posebilidad de la nueva formulación y establece la existencia y unicidad del control óptimo, junto con la derivación de las condiciones de optimalidad de primer orden.
• Estimaciones de Error A Priori: Proporciona un análisis de error a priori completo para el esquema totalmente discreto. Notablemente, deriva estimaciones de error para el desplazamiento y la presión en las normas $L^\infty(0, T; H^1(\Omega))$ y $L^2(0, T; L^2(\Omega))$, y para la variable de control en la norma $L^2(0, T; L^2(\Omega))$, logrando tasas de convergencia óptimas bajo supuestos de regularidad estándar.

Resultados

• Teóricos: El artículo demuestra que el esquema numérico converge con un orden de $O(h^s + \Delta t)$, donde $h$ es el tamaño de la malla espacial y $\Delta t$ es el paso de tiempo, y $s$ se relaciona con la regularidad de la solución. Específicamente, si la solución es suficientemente regular ($s=1$), la convergencia es lineal tanto en el espacio como en el tiempo ($O(h + \Delta t)$).
• Numéricos: Un experimento numérico en un dominio de cuadrado unitario valida los hallazgos teóricos. Utilizando la biblioteca de código abierto FEniCS, los autores resolvieron el sistema de optimalidad mediante un algoritmo de ida y vuelta (forward-backward). Los resultados demuestran que el esquema propuesto alcanza las tasas de convergencia óptimas predichas para las variables de desplazamiento, presión y control utilizando elementos de Taylor-Hood y MINI de bajo orden.

Significancia
Los autores posicionan este trabajo como un paso fundacional en el análisis numérico de problemas de control óptimo gobernados por sistemas poroelásticos. Señalan que, si bien el control óptimo para el flujo de fluidos en medios porosos deformables ha sido estudiado (por ejemplo, por Bociu et al. bajo la suposición de un coeficiente de almacenamiento evanescente), la aproximación numérica de tales problemas para el modelo de Biot general no había sido investigada previamente. Al establecer la bien posebilidad de una nueva formulación simétrica y proporcionar estimaciones de error rigurosas para un esquema totalmente discreto, este artículo llena un vacío en la literatura, ofreciendo un marco teórico robusto para controlar procesos poroelásticos en aplicaciones como el modelado de tejidos biológicos y simulaciones geofísicas.