Fully discrete error analysis of finite element discretizations of time-dependent Stokes equations in a stream-function formulation

Enunciado del Problema
Este artículo aborda el análisis de error totalmente discreto de las discretizaciones por elementos finitos para las ecuaciones de Stokes dependientes del tiempo en dos dimensiones, utilizando específicamente una formulación de función de corriente. El problema clásico de Stokes consiste en encontrar la velocidad $\mathbf{u}$ y la presión $p$ que satisfacen las ecuaciones de momento y continuidad. En el enfoque de función de corriente, la velocidad se representa como el rotacional vectorial de una función de corriente escalar $\psi$ ($\mathbf{u} = \text{rot}(\psi)$), lo cual satisface automáticamente la condición de divergencia nula. Esta transformación reduce el sistema a una única ecuación parabólica de cuarto orden que involucra el operador biarmónico espacial ($\Delta^2$) y la derivada temporal del laplaciano ($\partial_t \Delta$).

Los autores destacan una limitación crítica en los métodos estándar de elementos finitos (tales como los elementos Taylor-Hood o MINI) para problemas de Stokes: las estimaciones de error para la velocidad a menudo dependen de la aproximación de la presión. En consecuencia, las perturbaciones en el término fuente por un campo gradiente pueden no ser absorbidas por la presión discreta, lo que conduce a errores de velocidad subóptimos. La formulación de función de corriente ofrece una vía hacia métodos "robustos frente a la presión", pero ha faltado una teoría rigurosa para estimaciones de error totalmente discretas (combinando discretización espacial y temporal) para este problema específico de cuarto orden dependiente del tiempo.

Metodología
El artículo emplea un marco Galerkin totalmente discreto con los siguientes componentes:

1. Discretización Temporal: Los autores utilizan el método de Galerkin discontinuo (dG) de grado arbitrario $r$. Esto implica dividir el intervalo de tiempo en subintervalos y utilizar polinomios a trozos, permitiendo saltos en los nodos temporales.
2. Discretización Espacial: La discretización espacial se presenta en un marco general en lugar de estar vinculada a un tipo de elemento específico. El análisis se basa en una forma bilineal discreta simétrica y elíptica $a_h(\cdot, \cdot)$ y un operador discreto $A_h$. El marco está diseñado para acomodar diversos métodos, siempre que satisfagan condiciones específicas de ortogonalidad de Galerkin. El artículo cubre explícitamente:
   • Elementos finitos conformes $C^1$ (por ejemplo, elementos de Argyris).
   • Métodos de penalización interior $C^0$ (utilizando elementos de Lagrange estándar de grado $\ge 2$).
3. Análisis de Regularidad: Una parte significativa del trabajo se dedica a establecer la regularidad de la solución continua $\psi$. Los autores prueban resultados de existencia, unicidad y estabilidad para el problema continuo bajo supuestos mínimos de regularidad en los datos ($f \in L^2(I; H^{-1}(\Omega))$ y $\psi_0 \in H^2_0(\Omega)$). Crucialmente, demuestran que, aunque la solución posee ciertas propiedades de estabilidad, no se puede garantizar una regularidad superior de la derivada temporal $\partial_t \psi$ (específicamente en $L^2$ o $H^{-1}$) bajo estos supuestos mínimos, un hecho probado mediante un contraejemplo en el Apéndice A.

Contribuciones y Resultados Clave
La contribución principal es el establecimiento de estimaciones de error de tipo mejor aproximación para la solución totalmente discreta $\psi_{kh}$. El resultado principal (Teorema 4.5) establece que el error en el gradiente de la función de corriente está acotado por la suma de tres errores de aproximación:
$$\|\nabla(\psi - \psi_{kh})\|_{L^2(I \times \Omega)} \le C \left( \|\nabla(\psi - \chi_{kh})\|_{L^2(I \times \Omega)} + \|\nabla(\psi - R_h\psi)\|_{L^2(I \times \Omega)} + \|\nabla(\psi - \pi_k\psi)\|_{L^2(I \times \Omega)} \right)$$
donde:

• $\chi_{kh}$ es un elemento arbitrario del espacio discreto.
• $R_h$ es la proyección de Galerkin asociada al problema biarmónico estacionario.
• $\pi_k$ es una proyección temporal específica.

Al aplicar las propiedades de aproximación estándar de estas proyecciones bajo los supuestos de regularidad establecidos, los autores derivan tasas de convergencia óptimas en dominios poligonales convexos:
$$\|\nabla(\psi - \psi_{kh})\|_{L^2(I \times \Omega)} \le C(k + h^2) \left( \|f\|_{L^2(I; H^{-1}(\Omega))} + \|\psi_0\|_{H^2_0(\Omega)} \right)$$
Esto implica una convergencia de orden $O(k)$ en el tiempo y $O(h^2)$ en el espacio para el error de velocidad (dado que $\|\mathbf{u} - \mathbf{u}_{kh}\| = \|\nabla(\psi - \psi_{kh})\|$).

Significado y Afirmaciones
El artículo se posiciona como un paso fundamental en el análisis numérico de formulaciones de función de corriente para problemas dependientes del tiempo.

• Robustez frente a la presión: Las estimaciones derivadas no contienen términos de presión, abordando el problema conocido de que los métodos estándar sufren errores de velocidad dependientes de las aproximaciones de presión. Los experimentos numéricos confirman que el método de función de corriente permanece robusto incluso cuando el término fuente se modifica mediante un campo gradiente, mientras que los métodos mixtos estándar (elementos MINI) muestran un rendimiento degradado.
• Supuestos mínimos de regularidad: El análisis no requiere supuestos adicionales de regularidad más allá de los proporcionados naturalmente por el dominio y los datos. Los autores argumentan explícitamente que trabajos anteriores a menudo imponían supuestos excesivamente estrictos (por ejemplo, requerir que $\partial_{tt}\psi$ estuviera acotado) para lograr estimaciones similares.
• Generalidad: El marco es aplicable a una amplia variedad de métodos de discretización espacial, siempre que satisfagan las condiciones de ortogonalidad necesarias, lo que hace que los resultados sean versátiles para métodos conformes $C^1$ y métodos de penalización interior no conformes $C^0$.
• Control Óptimo: Los autores señalan que los supuestos de baja regularidad hacen que estos resultados sean adecuados para su aplicación en problemas de control óptimo.

El artículo concluye con experimentos numéricos que ilustran las tasas de convergencia teóricas y demuestran la robustez frente a la presión del método propuesto en comparación con las formulaciones estándar de elementos finitos mixtos. Los autores reconocen esto como un paso inicial, con trabajos futuros planeados para extender estos resultados a flujos superficiales, dominios múltiplemente conexos y las ecuaciones de Navier-Stokes.