Achieving Material Robustness via Symmetric Stress Finite Element Discretizations

Resumen Técnico: Logro de Robustez Material mediante Discretizaciones de Elementos Finitos de Tensión Simétrica

Enunciado del Problema
En mecánica de medios continuos, formulaciones variacionales como la elasticidad de Hellinger–Reissner y las formulaciones de velocidad-tensión para flujo incompresible involucran tensores de tensión simétricos. La simetría del tensor de tensión de Cauchy está físicamente ligada a la conservación del momento angular. Al discretizar estos problemas utilizando elementos finitos conformes en $H(\text{div})$, los practicantes deben elegir entre imponer la simetría fuertemente (requiriendo que los tensores discretos sean simétricos punto a punto) o débilmente (imponiendo la simetría mediante un multiplicador de Lagrange). Aunque la simetría débil suele preferirse en la literatura debido a la complejidad percibida de construir elementos fuertemente simétricos, este artículo investiga el impacto de esta elección en la precisión numérica, particularmente para materiales con leyes constitutivas anisotrópicas.

Metodología
Los autores desarrollan un marco teórico y una serie de benchmarks numéricos para comparar esquemas simétricos fuertemente y débilmente.

1. Formulación Matemática: El estudio considera un problema estático linealizado general tanto para fluidos como para sólidos, donde la tensión $\sigma$ se relaciona con la deformación $\varepsilon(u)$ y un término de campo externo $\tilde{F}$ (que representa la anisotropía). El problema se plantea como un sistema de punto de silla.
2. Esquemas de Discretización:
   • Fuertemente Simétricos: Se utilizan el esquema de Johnson-Mercier-Křížek (JMK) y el esquema de Hu-Zhang (HZk). Estos imponen $\sigma_h \in \Sigma^{\text{sym}}_h$ punto a punto.
   • Débilmente Simétricos: Se utilizan los esquemas PEERS y Arnold-Falk-Winther (AFWk). Estos introducen un multiplicador de Lagrange $\omega_h$ (que representa la parte antisimétrica del gradiente) para imponer la simetría en sentido débil.
3. Construcción de Benchmarks: Los autores construyen "soluciones manufacturadas" donde la tensión exacta es idénticamente cero ($\sigma \equiv 0$) pero los campos de desplazamiento y rotación son no triviales. Estas soluciones se derivan de leyes constitutivas inspiradas en:
   • Sólidos isotrópicos lineales.
   • Sólidos transversalmente isotrópicos (materiales reforzados con fibras).
   • Fluidos polares (redes de polímeros de cristal líquido).
     La magnitud del término de forzamiento anisotrópico se controla mediante un parámetro de escala $\delta$.
4. Análisis Teórico: Se desarrolla una teoría abstracta de espacios de Hilbert para problemas genéricos de punto de silla. Los autores definen discretizaciones "preservadoras de estructura" y derivan estimaciones de error a priori para explicar el comportamiento numérico observado.

Contribuciones Clave

• Definición de Robustez Material: El artículo introduce el concepto de "robustez material", definido como la propiedad por la cual un esquema numérico produce una aproximación de tensión cero para cualquier configuración exacta de tensión cero, independientemente de la ley constitutiva o de la magnitud del forzamiento anisotrópico.
• Identificación de Modos de Fallo: Los autores demuestran que los esquemas débilmente simétricos pueden producir aproximaciones de tensión arbitrariamente pobres (errores que escalan con el parámetro de forzamiento $\delta$) incluso cuando la tensión exacta es cero. Esto ocurre específicamente cuando el tensor de rotación exacto $\omega = \text{skw}(\nabla u)$ no pertenece al espacio discreto de multiplicadores de Lagrange $\Xi_h$.
• Teoría Unificadora: Se presenta una teoría rigurosa que muestra que los esquemas fuertemente simétricos que satisfacen propiedades específicas de inclusión de núcleo ($\ker B_h \subseteq \ker B$) son inherentemente robustos materialmente. Por el contrario, los esquemas débilmente simétricos carecen de esta invariancia a menos que el espacio de rotación discreto contenga el tensor de rotación exacto.
• Estimaciones de Error: Se derivan nuevas estimaciones de error a priori para discretizaciones parcialmente preservadoras de estructura, mostrando explícitamente cómo el error de tensión depende de la calidad de aproximación del tensor de rotación en el caso débilmente simétrico.

Resultados
Los experimentos numéricos en mallas no estructuradas confirman las predicciones teóricas:

• Caso Isotrópico: Para sólidos isotrópicos lineales donde la rotación exacta es constante (y por tanto fácilmente representable), tanto los esquemas fuertes como los débiles se comportan de manera similar, con errores dominados por el redondeo.
• Casos Anisotrópicos: Para sólidos transversalmente isotrópicos y fluidos polares, donde el tensor de rotación es complejo:
  • Esquemas Fuertemente Simétricos (JMK, HZ3): Producen consistentemente errores de tensión al nivel de la precisión de la máquina (estado de tensión cero), independientemente del parámetro de escala $\delta$.
  • Esquemas Débilmente Simétricos (PEERS, AFWk): Producen errores de tensión que escalan linealmente con $\delta$. Incluso para esquemas de muy alto orden (por ejemplo, AFW3), si el tensor de rotación exacto no está contenido en el espacio discreto $\Xi_h$, el error de tensión permanece grande.
• Comportamiento Transitorio: Un ejemplo de fluido polar transitorio muestra que, aunque ambos esquemas coinciden inicialmente, el esquema débilmente simétrico no robusto falla al relajarse al estado estacionario de tensión cero correcto, mientras que el esquema fuertemente simétrico robusto sí lo hace.

Significado y Afirmaciones
El artículo afirma que la elección entre la imposición de simetría fuerte y débil no es meramente una cuestión de complejidad de implementación, sino que afecta fundamentalmente la fidelidad física de la solución numérica. Los autores argumentan que imponer la simetría fuertemente es crucial para la robustez material, asegurando que los esquemas numéricos no generen tensiones espurias en configuraciones libres de tensión, particularmente en aplicaciones que involucran invariantes rotacionales o materiales anisotrópicos.

El trabajo aboga por una adopción más amplia de elementos finitos fuertemente simétricos (como JMK y Hu-Zhang), señalando que estos elementos están cada vez más disponibles en software de código abierto como Firedrake y NGSolve. Los autores concluyen que, para problemas sensibles a la conservación del momento angular, las discretizaciones fuertemente simétricas proporcionan una garantía necesaria de precisión que los esquemas débilmente simétricos no pueden ofrecer sin condiciones restrictivas sobre los espacios discretos.