A hybrid $C^{0}$-interior penalty method for the nematic Helmholtz--Korteweg equation

Resumen Técnico: Un Método de Penalización Interior Híbrido $C^0$ para la Ecuación de Helmholtz–Korteweg Nemática

Planteamiento del Problema
Este trabajo aborda la solución numérica de la ecuación de Helmholtz–Korteweg nemática, una ecuación diferencial parcial de cuarto orden que modela la propagación de ondas acústicas de frecuencia armónica en fluidos de Korteweg nemáticos, como los cristales líquidos nemáticos. La ecuación gobernante aumenta la ecuación de Helmholtz clásica con dos términos adicionales de cuarto orden, uno de los cuales es anisotrópico con respecto al campo del director nemático $\mathbf{n}$. El problema se formula en un dominio Lipschitz acotado $\Omega \subset \mathbb{R}^d$ ($d=2,3$) con condiciones de contorno mixtas de tipo suave-acústico (sound-soft) e impedancia.

El análisis enfrenta tres desafíos principales:

1. Alto Orden y Anisotropía: La presencia del término anisotrópico de cuarto orden complica el análisis en comparación con los problemas estándar de Helmholtz o biarmónicos.
2. Regularidad y Estabilidad: Los análisis previos dependían de argumentos de perturbación alrededor del caso isotrópico ($\beta=0$), lo que requería que el parámetro de anisotropía $\beta$ fuera pequeño en relación con una constante de estabilidad del problema isotrópico. Esta constante degenera cerca de las frecuencias resonantes, forzando artificialmente a $\beta \to 0$ incluso cuando el problema anisotrópico podría ser bien planteado.
3. Complejidad de la Discretización: Los métodos de elementos finitos conformes estándar para problemas de cuarto orden ($C^1$-conformes) requieren grados polinómicos altos (por ejemplo, los elementos de Argyris requieren $p \geq 5$ en 2D y $p \geq 9$ en 3D) y son notoriamente difíciles de implementar en dominios curvos o en tres dimensiones.

Metodología
El artículo propone una discretización de penalización interior híbrida $C^0$ ($C^0$-HIP). Este enfoque utiliza un espacio de elementos finitos conforme a $H^1$ para las incógnitas de volumen y técnicas de Galerkin discontinuo (DG) para manejar la falta de conformidad $H^2$.

• Análisis Continuo: Los autores refinan el análisis continuo empleando la condición de Cordes para cuantificar el efecto anisotrópico. Establecen que el problema es Fredholm de índice cero bajo una condición de pequeñez mejorada sobre el parámetro de anisotropía $\beta$. Crucialmente, esta condición es independiente de la constante de estabilidad del caso isotrópico, lo que permite un rango más amplio de parámetros físicamente relevantes.
• Discretización: El método introduce incógnitas de faceta adicionales que representan la derivada normal de las incógnitas de volumen. Estas variables de faceta desacoplan las incógnitas de volumen, permitiendo la condensación estática para reducir los costos computacionales. La formulación discreta utiliza un operador de levantamiento para definir un Laplaciano discreto e incluye un término de estabilización para controlar los saltos en la derivada normal a través de las interfaces de los elementos.
• Marco Teórico: El análisis se basa en un marco abstracto para discretizaciones no conformes de operadores que son perturbaciones compactas de operadores coercitivos. Los autores demuesten una desigualdad de Gårding discreta, que reproduce la estructura de "coercitividad hasta la compacidad" del problema continuo en el nivel discreto. Esta estructura, combinada con la consistencia asintótica y un argumento de compacidad discreta, asegura la estabilidad y la convergencia.

Contribuciones Clave

1. Análisis de Bien Planteamiento Mejorado: El artículo demuestra que la ecuación de Helmholtz–Korteweg nemática está bien planteada para un rango significativamente más amplio del parámetro de anisotropía $\beta$ de lo establecido anteriormente. Al utilizar la condición de Cordes, la condición de pequeñez derivada no degenera en las frecuencias resonantes del problema isotrópico.
2. Método de Penalización Interior Híbrido $C^0$: Los autores introducen una discretización $C^0$-HIP que evita las dificultades de implementación de los elementos conformes $C^1$ (como los elementos de Argyris o HCT), mientras requiere menos grados de libertad que los métodos de Galerkin discontinuo completos.
3. Pruebas de Estabilidad y Convergencia: Se demuestra que el método es estable para cualquier grado polinómico $p \geq 2$, independiente de la dimensión espacial, siempre que la anisotropía sea suficientemente pequeña. El artículo establece la convergencia a la solución continua bajo supuestos mínimos de regularidad ($u \in H^2(\Omega)$) y deriva tasas de convergencia de orden óptimo bajo regularidad adicional.
4. Simulaciones Tridimensionales: El trabajo presenta las primeras simulaciones numéricas en tres dimensiones de los efectos anisotrópicos de los campos nemáticos en la propagación de ondas acústicas utilizando esta ecuación específica.

Resultados

• Resultados Teóricos: La forma sesquilineal discreta satisface una desigualdad de Gårding discreta, asegurando la estabilidad asintótica. La secuencia de soluciones discretas converge a la solución continua cuando el tamaño de malla $h \to 0$.
• Experimentos Numéricos:
  • Convergencia: Las pruebas numéricas en un disco unitario confirman las tasas de convergencia predichas de $O(h^{p-1})$ para grados polinómicos $p \in \{2, 3, 4, 5\}$.
  • Estabilidad: Los experimentos verifican la estabilidad para $\beta > -\alpha$ en 2D y $\beta \in (-\alpha, 3\alpha)$ en 3D, consistente con los límites teóricos derivados de la condición de Cordes.
  • Anisotropía: Las simulaciones en 2D y 3D demuestran la naturaleza anisotrópica de la propagación de ondas. Los resultados muestran que la velocidad de onda y la atenuación varían dependiendo de la orientación del director nemático $\mathbf{n}$, lo que conduce a patrones de dispersión distintos en comparación con el caso isotrópico.
  • Eficiencia: Una comparación de grados de libertad y entradas no nulas muestra que la estrategia de hibridación reduce significativamente los costos computacionales para grados polinómicos más altos en comparación con los métodos estándar de penalización interior $C^0$, debido a la eliminación de burbujas interiores mediante condensación estática.

Significado
El artículo afirma que su principal significancia radica en proporcionar un marco numérico robusto y flexible para resolver la ecuación de Helmholtz–Korteweg nemática, la cual es esencial para modelar la nematoacústica. Al desacoplar el análisis de estabilidad de la constante de estabilidad isotrópica, este trabajo elimina una limitación artificial sobre el parámetro de anisotropía, permitiendo la simulación de efectos nemáticos más fuertes. Además, el método $C^0$-HIP propuesto supera las barreras prácticas de los métodos conformes $C^1$, haciendo que las simulaciones de alto orden de problemas anisotrópicos de cuarto orden sean factibles en tres dimensiones y en dominios curvos. Los ejemplos numéricos sirven como una prueba de concepto de la capacidad del método para capturar fenómenos complejos de ondas anisotrópicas.