Métodos de Elementos Finitos Superconvergentes y Libres de Divergencia para la Ecuación de Stokes

Información Básica

• ID del Artículo: 2510.14192
• Título: Métodos de Elementos Finitos Superconvergentes y Libres de Divergencia para la Ecuación de Stokes
• Autores: Long Chen, Xuehai Huang, Chao Zhang, Xinyue Zhao
• Clasificación: math.NA cs.NA
• Fecha de Presentación: 16 de octubre de 2025 en arXiv
• Enlace del Artículo: https://arxiv.org/abs/2510.14192

Resumen

Este artículo desarrolla métodos de elementos finitos superconvergentes y libres de divergencia para resolver la ecuación de Stokes. El campo de velocidad se discretiza utilizando elementos vectoriales compatibles con H(div), mientras que el campo de presión se discretiza mediante polinomios discontinuos por partes. El esquema discreto utiliza un operador de gradiente de desviación débil construido sobre elementos finitos continuos tangencial-normal basados en tensores sin traza, sin requerir estabilización. Se establecen estimaciones de error óptimas y superconvergentes. El método se relaciona con métodos de elementos virtuales no conformes y esquemas mixtos de pseudo-esfuerzo-velocidad-presión. Experimentos numéricos validan los resultados teóricos.

Antecedentes de Investigación y Motivación

Importancia del Problema

La ecuación de Stokes es la ecuación fundamental que describe el movimiento de fluidos incompresibles y ocupa un lugar central en la dinámica de fluidos computacional. La solución numérica de esta ecuación enfrenta dos desafíos principales:

1. Restricción de divergencia nula: El campo de velocidad debe satisfacer div u = 0, que es la manifestación de la conservación de masa
2. Robustez de presión: El método numérico debe ser robusto en el tratamiento del término de presión

Limitaciones de Métodos Existentes

Los pares de elementos finitos estabilizados clásicos (como elementos Taylor-Hood, MINI, P1-P0 no conforme) presentan los siguientes problemas:

• Divergencia no puntual: No pueden satisfacer estrictamente la restricción de divergencia nula, debilitando la conservación de masa
• Robustez insuficiente de presión: Carecen de robustez de presión en las estimaciones de error
• Limitaciones en el orden de convergencia: Algunos métodos sufren pérdida de orden de convergencia

Motivación de la Investigación

Este artículo tiene como objetivo construir un método de elementos finitos que mantenga la propiedad de divergencia nula estricta y logre superconvergencia, evitando simultáneamente la estabilización requerida por los métodos tradicionales.

Contribuciones Principales

1. Propone un nuevo esquema mixto de elementos finitos: Discretiza la velocidad con elementos compatibles con H(div) y la presión con polinomios discontinuos por partes
2. Construye un operador de gradiente de desviación débil: Basado en elementos finitos continuos tangencial-normal sin traza, sin requerir estabilización
3. Establece estimaciones de error superconvergentes: Demuestra superconvergencia de orden $h^{k+1}$
4. Proporciona técnicas de postprocesamiento: Construye aproximaciones de velocidad superconvergentes
5. Establece relaciones de equivalencia: Conexiones con métodos de elementos virtuales no conformes y esquemas de pseudo-esfuerzo

Explicación Detallada del Método

Definición del Problema

Considere la ecuación de Stokes en un dominio acotado $\Omega \subset \mathbb{R}^d$: