The Tempered Finite Element Method

Resumen Técnico: El Método de Elementos Finitos Templado (TFEM)

Planteamiento del Problema
El Método de Elementos Finitos (FEM) es la técnica numérica dominante en la física computacional y la ingeniería, reconocida por su capacidad para manejar geometrías complejas mediante mallas no estructuradas. Sin embargo, las formulaciones estándar de FEM exhiben fragilidad al encontrarse con configuraciones de malla específicas, particularmente aquellas que contienen "capas" (triángulos con un ángulo que se aproxima a los 180°) dispuestas en bandas. Mientras que los elementos degenerados aislados (agujas o capas aisladas) no impiden la convergencia, las bandas de capas inducen un fenómeno de "bloqueo" (locking). En tales configuraciones, la solución numérica se ve forzada a ser lineal a lo largo de la banda, independientemente de la solución real, lo que provoca que el tamaño efectivo de la malla esté determinado por la longitud de la banda en lugar de su grosor. Esto resulta en una pérdida de convergencia para el FEM estándar. Históricamente, la "condición del ángulo máximo" (que requiere que todos los ángulos estén alejados de los 180°) se consideraba esencial para la convergencia, aunque trabajos recientes han demostrado que no es estrictamente necesaria para elementos aislados. El desafío persiste en cómo manejar mallas con elementos de medida cero o casi degenerados, los cuales son relevantes para los métodos de seguimiento de interfaces (por ejemplo, X-Mesh), donde las interfaces se rastrean mediante el relevo de la malla, lo que puede crear bandas de capas.

Metodología
Los autores proponen el Método de Elementos Finitos Templado (TFEM), una modificación del FEM clásico que permite la convergencia en mallas que contienen bandas de capas. La idea central se inspira en el templado metalúrgico: así como el templado alivia las tensiones internas en el acero endurecido para aumentar su flexibilidad, el TFEM "suaviza" localmente las propiedades del material en la banda degenerada para permitir que la solución numérica se conforme al problema continuo.

Matemáticamente, para problemas elípticos como la ecuación de Poisson, la matriz de rigidez estándar $K_{ij}$ involucra el inverso del determinante del Jacobiano $J$ del mapeo del elemento. En elementos degenerados, $J \to 0$, lo que conduce a la división por cero o a valores propios infinitos. El TFEM modifica la matriz de rigidez local acotando el Jacobiano para que no sea inferior a un valor $J_{min}$:
$$\tilde{K}_{ij} = \int_{\hat{E}} \frac{1}{\max(J, J_{min})} [\dots] d\xi$$
Esto es equivalente a introducir un coeficiente de difusión variable $d(x)$ en la formulación débil:
$$d(x) = \begin{cases} 1 & x \in \Omega \setminus B \\ D = \frac{J}{J_{min}} & x \in B \end{cases}$$
donde $B$ es la banda de elementos degenerados. El parámetro $J_{min}$ se elige de forma dependiente del tamaño de la malla $h$ (específicamente $J_{min} \sim h^{d+1}$, donde $d$ es la dimensión espacial) para equilibrar la reducción del coeficiente de difusión (para evitar el bloqueo) con la preservación de las propiedades de aproximación.

Para el caso de elementos con medida exactamente cero ($J=0$), los autores demuestran que el límite del esquema TFEM corresponde a un método de mortero (mortaring method) con penalización del salto de la solución a través de la interfaz.

Contribuciones Clave

1. Prueba Teórica de Convergencia: El artículo proporciona un análisis de error riguroso para el TFEM en bandas de capas. Demuestra que, con una elección óptima del parámetro de templado $D \sim \bar{h}/h^2$ (donde $\bar{h}$ es el ancho de la banda), el método logra una convergencia óptima de $O(h)$ en la seminorma $H^1$, igualando la tasa del FEM estándar en mallas regulares.
2. Equivalencia con el Método de Mortero: Para bandas de ancho exactamente cero, los autores prueban que el TFEM es equivalente a un método de mortero no conforme con penalización de salto, proporcionando un puente teórico entre el manejo de mallas degeneradas y las técnicas de acoplamiento de interfaces.
3. Extensión a la Física Compleja: El método se extiende con éxito más allá del Laplaciano escalar a:
   • Elasticidad Lineal: Demostrando una precisión comparable al FEM estándar para el análisis de tensiones a través de interfaces.
   • Mortero de Mallas No Conformes: Mostrando que el TFEM puede conectar sin problemas mallas de diferentes resoluciones sin necesidad de elementos de transición.
   • Elementos de Alto Orden: Validando que el método preserva las tasas de convergencia óptimas para órdenes polinómicos hasta el 4, siempre que $J_{min}$ se escale apropiadamente con el orden del elemento.
   • Advección: Aplicando el TFEM a problemas de advección-difusión utilizando la estabilización de Streamline Upwind/Petrov-Galerkin (SUPG), observando que el templado solo es requerido para el término de estabilización.
4. Análisis de Robustez: Los experimentos numéricos identifican que el método es robusto con respecto a la constante $C$ en la relación $J_{min} = C h^{d+1}$. El error permanece cercano al óptimo sobre un amplio rango de valores de $C$ (abarcando casi un orden de magnitud), y el condicionamiento del sistema lineal permanece manejable, especialmente con precondicionadores estándar.

Resultados

• Experimentos Numéricos: En problemas de Poisson 2D y 3D con soluciones fabricadas, el TFEM en mallas con bandas de capas recupera las tasas de convergencia del FEM estándar en mallas regulares. Específicamente, elegir $J_{min} \sim h^3$ en 2D y $J_{min} \sim h^4$ en 3D produce una convergencia de error óptima en $L^2$ y $H^1$.
• Sensibilidad del Error: El error $H^1$ es menos sensible a la elección de $J_{min}$ que el error $L^2$, exhibiendo una meseta sobre un amplio rango de parámetros.
• Condicionamiento: Aunque el número de condición de la matriz de rigidez aumenta a medida que el ancho de la banda disminuye, el uso de precondicionamiento de Cholesky incompleto asegura que la convergencia del solucionador iterativo sea independiente del grosor de la banda.
• Pruebas Aleatorizadas: Una serie de 100 pruebas aleatorias en 2D confirmó que, incluso con una elección simple ($C=1$), los errores del TFEM son casi idénticos a los errores del FEM estándar en mallas regulares.

Significancia y Reivindicaciones
El artículo afirma que el TFEM aumenta significativamente la flexibilidad del análisis de elementos finitos al relajar los estrictos requisitos de calidad de malla asociados con los elementos degenerados. Permite el uso de mallas con elementos de medida cero o casi degenerados (como las generadas por el seguimiento dinámico de interfaces) sin sufrir los fenómenos de bloqueo que dejan inoperante al FEM estándar.

Los autores enfatizan que el método es trivial de implementar en códigos de FEM existentes, requiriendo solo una modificación en la evaluación del Jacobiano para evitar la división por cero. Posicionan al TFEM no como un reemplazo del FEM estándar en mallas de alta calidad, sino como una extensión robusta que permite el uso de topologías de malla específicas (bandas de capas) que anteriormente se consideraban inoperativas. El método es particularmente valioso para aplicaciones como el método de seguimiento de interfaz X-Mesh, donde mantener una malla conforme es difícil o computacionalmente costoso. El artículo concluye que, si bien se deriva la elección teórica óptima del parámetro $C$, la robustez práctica del método permite una amplia gama de valores aceptables, lo que lo hace inmediatamente aplicable.