Analysis of a finite element method for second order uniformly elliptic PDEs in non-divergence form

A continuación se presenta un resumen técnico detallado del artículo "Análisis de un Método de Elementos Finitos para EDPs de Segundo Orden Uniformemente Elípticas en Forma No Divergente" de Weifeng Qiu.

1. Planteamiento del Problema

El artículo aborda el análisis numérico de dos clases de ecuaciones en derivadas parciales (EDP) de segundo orden definidas sobre dominios poliedrales de Lipschitz $\Omega \subset \mathbb{R}^d$ ($d=2, 3$):

1. EDP Elíptica Lineal en Forma No Divergente:
   $$A : D^2u + b \cdot \nabla u + cu = f \quad \text{en } \Omega, \quad u = 0 \text{ en } \partial\Omega$$
   donde $A$ es una matriz simétrica, uniformemente definida positiva, y $b, c$ son coeficientes acotados.

2. Ecuación de Hamilton-Jacobi-Bellman (HJB) Uniformemente Elíptica:
   $$\sup_{\alpha \in \Lambda} [A_\alpha : D^2u + b_\alpha \cdot \nabla u + c_\alpha u - f_\alpha] = 0 \quad \text{en } \Omega, \quad u = 0 \text{ en } \partial\Omega$$
   donde $\Lambda$ es un conjunto de índices, y los coeficientes satisfacen condiciones de elipticidad uniforme y acotación.

Desafíos Clave:

• Forma No Divergente: Las formulaciones débiles estándar (basadas en integración por partes) no son directamente aplicables porque el término de segundo orden $A:D^2u$ no puede transferirse fácilmente a la función de prueba sin asumir alta regularidad de $A$ o del dominio.
• Espacio de Soluciones: El análisis se centra en soluciones fuertes en el espacio de Sobolev $W^{2,p}(\Omega)$ en lugar de soluciones de viscosidad.
• Geometría del Dominio: El artículo busca manejar dominios poliedrales de Lipschitz, incluyendo polígonos no convexos en 2D, que son más generales que los dominios $C^{1,1}$ o convexos requeridos frecuentemente en la literatura existente.
• Regularidad de los Coeficientes: El método debe acomodar coeficientes discontinuos (específicamente para la ecuación HJB) manteniendo tasas de convergencia óptimas.

2. Metodología

El autor propone un Método de Elementos Finitos (MEF) unificado tanto para la EDP lineal como para la ecuación HJB.

A. Discretización y Operadores

• Malla: Una triangulación simplicial, conforme y cuasi-uniforme $\mathcal{T}_h$.
• Espacio de Elementos Finitos: $V_h = H^1_0(\Omega) \cap P_r(\mathcal{T}_h)$ con grado polinomial $r \geq 2$.
• Gradiente Discreto ($\nabla_h$): El método utiliza un operador de gradiente discreto específico definido mediante integración por partes en los elementos de la malla y términos de salto a través de las caras interiores. Este operador permite la aproximación del Hessiano $D^2u$ sin requerir continuidad $C^1$ de las funciones base.
  $$(\nabla_h v)_i = \partial_{h, x_i} v$$
• Funciones de Peso: Para manejar la forma no divergente, las ecuaciones se multiplican por una función de peso $\gamma$ (o $\gamma_\alpha$ para HJB).
  • Para el caso lineal: $\gamma(x) = \frac{\text{Tr} A(x)}{|A(x)|^2}$.
  • Esto transforma la ecuación en una forma apta para el operador discreto: $\tilde{A} : D^2 u + \dots = \tilde{f}$.

B. Los Esquemas Numéricos

• Caso Lineal (Ecuación 14): Encontrar $u_h \in V_h$ tal que para todo $v_h \in V_h$:
  $$\int_\Omega (\tilde{A} : \nabla_h(\nabla u_h) + \tilde{b} \cdot \nabla u_h + \tilde{c} u_h) v_h \, dx = \int_\Omega \tilde{f} v_h \, dx$$
• Caso HJB (Ecuación 16): Encontrar $u_h \in V_h$ tal que para todo $v_h \in V_h$:
  $$\left( \sup_{\alpha \in \Lambda} \gamma_\alpha [A_\alpha : \nabla_h(\nabla u_h) + b_\alpha \cdot \nabla u_h + c_\alpha u_h - f_\alpha], v_h \right)_\Omega = 0$$
  Nota: Si $\Lambda$ contiene un solo elemento, el esquema HJB se reduce exactamente al esquema lineal.

C. Marco Analítico

El análisis se basa en normas $W^{2,p}$ discretas definidas como:
$$\|v\|_{W^{2,p}_h(\Omega)} = \|D^2_h v\|_{L^p(\Omega)} + \left( \sum_{F \in \mathcal{F}_h^I} h_F^{1-p} \|\llbracket \nabla v \rrbracket\|_{L^p(F)}^p \right)^{1/p}$$
Las herramientas analíticas clave incluyen:

1. Estimaciones Globales $W^{2,p}$: Demostrar la buena posición de la solución fuerte para el problema continuo en poliedros de Lipschitz.
2. Compacidad Discreta: Establecer que las sucesiones acotadas en la norma discreta $W^{2,p}_h$ tienen subsucesiones convergentes en $W^{1,p}$.
3. Argumentos de Unicidad: Utilizar la unicidad de la solución fuerte para probar la estabilidad mediante contradicción (evitando argumentos de dualidad especializados utilizados en trabajos anteriores).
4. Medibilidad del Supremum: Para la ecuación HJB, probar que el supremo sobre el conjunto de índices $\Lambda$ es medible y puede aproximarse por un subconjunto numerable $\tilde{\Lambda}$ (Suposición 6).

3. Contribuciones Clave

A. Marco Unificado

El artículo demuestra que una única formulación MEF puede resolver eficazmente tanto EDPs lineales en forma no divergente como ecuaciones HJB totalmente no lineales. El esquema HJB degenera naturalmente al esquema lineal cuando el conjunto de índices es un singleton.

B. Relajación de Suposiciones sobre Dominio y Coeficientes

• Dominios: El análisis extiende los resultados de convergencia óptima desde dominios $C^{1,1}$ o estrictamente convexos hacia dominios poliedrales de Lipschitz.
  • Para $d=3$, $\Omega$ debe ser convexo.
  • Para $d=2$, $\Omega$ puede ser no convexo (con restricciones sobre $p$).
• Coeficientes: El método permite coeficientes discontinuos en la ecuación HJB, siempre que satisfagan una condición específica de aproximación (Suposición 6) que involucra un subconjunto denso numerable. Esto relaja las suposiciones de continuidad ($A \in C^0$) requeridas por muchos métodos de elementos finitos discontinuos (DG) existentes.

C. Convergencia Óptima en $W^{2,p}$

El artículo prueba tasas de convergencia óptimas en la norma discreta $W^{2,p}_h$ para un amplio rango de $p$:

• Poliedros Convexos ($d=2,3$): $p \in (1, 2]$.
• Polígonos No Convexos ($d=2$): $p \in (4/3 - \epsilon_1, 4/3 + \epsilon_2)$.
  Esto representa una mejora significativa sobre los métodos que solo garantizan convergencia en $L^2$ o $H^2$ (donde $p=2$).

D. Nueva Prueba de Estabilidad

A diferencia de trabajos anteriores (por ejemplo, Feng et al.) que dependen de argumentos de dualidad especializados condicionados a la buena posición del problema continuo, este artículo establece la estabilidad utilizando compacidad discreta y la unicidad de la solución fuerte. Este enfoque es más robusto para dominios no convexos y coeficientes discontinuos.

4. Resultados Principales

1. Buena Posición (Teoremas 3.3, 4.1):

   • Se demuestra la existencia y unicidad de soluciones fuertes en $W^{2,p}(\Omega) \cap W^{1,p}_0(\Omega)$ tanto para la EDP lineal como para la ecuación HJB bajo las condiciones especificadas sobre $A, b, c$ y el dominio.
   • Para la ecuación HJB, la condición de Cordes se generaliza para permitir coeficientes discontinuos mediante la suposición de subconjunto numerable.

2. Estabilidad (Teoremas 3.2, 4.3):

   • La solución discreta $u_h$ satisface una estimación de estabilidad: $\|u_h\|_{W^{2,p}_h} \leq C \|f\|_{L^p}$.
   • Esta estabilidad se mantiene para el rango específico de $p$ dependiente de la geometría del dominio.

3. Convergencia Óptima (Teoremas 3.3, 4.3):

   • Se deriva la estimación de error:
     $$\|u - u_h\|_{W^{1,p}} + \|u - u_h\|_{W^{2,p}_h} \leq C \left( \inf_{\chi_h \in V_h} (\|u - \chi_h\|_{W^{1,p}} + \|u - \chi_h\|_{W^{2,p}_h}) + \|P_h(D^2u) - D^2u\|_{L^p} \right)$$
   • Esto implica que si la solución exacta es suficientemente suave, el método alcanza tasas de convergencia óptimas en la norma discreta $W^{2,p}$.

5. Significado

• Puente entre Teoría y Práctica: El artículo proporciona una base teórica rigurosa para el uso de espacios estándar de elementos finitos (con continuidad $C^0$) para resolver EDPs en forma no divergente, las cuales son notoriamente difíciles de discretizar.
• Manejo de No Convexidad: Al extender el análisis a polígonos 2D no convexos, el método se vuelve aplicable a una clase más amplia de problemas prácticos de ingeniería y física donde los dominios rara vez son perfectamente convexos.
• Aplicaciones HJB: La relajación de las suposiciones de continuidad de los coeficientes hace que el método sea adecuado para problemas de control estocástico donde los coeficientes pueden ser discontinuos (por ejemplo, regímenes de conmutación), un escenario donde los métodos conformes $C^0$ anteriores fallaban.
• Enfoque Unificado: Demostrar que un único esquema numérico maneja tanto casos lineales como totalmente no lineales simplifica la implementación y el análisis teórico para sistemas complejos.

En resumen, el trabajo de Weifeng Qiu avanza el estado del arte en EDPs numéricas al establecer un marco unificado y robusto de elementos finitos para ecuaciones en forma no divergente que logra convergencia óptima en normas fuertes sobre dominios poliedrales generales con requisitos de regularidad relajados.