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

Resumo Técnico: Um Método de Penalidade de Interior Híbrido $C^0$ para a Equação de Helmholtz–Korteweg Nemática

Enunciado do Problema
Este trabalho aborda a solução numérica da equação de Helmholtz–Korteweg nemática, uma equação diferencial parcial de quarta ordem que modela a propagação de ondas acústicas harmônicas no tempo em fluidos de Korteweg nemáticos, como cristais líquidos nemáticos. A equação governante aumenta a equação de Helmholtz clássica com dois termos adicionais de quarta ordem, sendo um deles anisotrópico em relação ao campo de diretor nemático $\mathbf{n}$. O problema é formulado em um domínio Lipschitz limitado $\Omega \subset \mathbb{R}^d$ ($d=2,3$) com condições de contorno mistas de som suave (sound-soft) e impedância.

A análise enfrenta três desafios primários:

1. Alta Ordem e Anisotropia: A presença do termo anisotrópico de quarta ordem complica a análise em comparação com os problemas padrão de Helmholtz ou biarmônicos.
2. Regularidade e Estabilidade: Análises anteriores basearam-se em argumentos de perturbação em torno do caso isotrópico ($\beta=0$), o que exigia que o parâmetro de anisotropia $\beta$ fosse pequeno em relação a uma constante de estabilidade do problema isotrópico. Esta constante degenera próximo a frequências ressonantes, forçando artificialmente $\beta \to 0$ mesmo quando o problema anisotrópico poderia permanecer bem posto.
3. Complexidade de Discretização: Métodos de elementos finitos conformes padrão para problemas de quarta ordem ($C^1$-conformes) requerem graus polinomiais elevados (por exemplo, elementos de Argyris exigem $p \geq 5$ em 2D e $p \geq 9$ em 3D) e são notoriamente difíceis de implementar em domínios curvos ou em três dimensões.

Metodologia
O artigo propõe uma discretização $C^0$-híbrida de penalidade de interior ($C^0$-HIP). Esta abordagem utiliza um espaço de elementos finitos conforme $H^1$ para as incógnitas de volume e técnicas de Galerkin descontínuo (DG) para lidar com a falta de conformidade $H^2$.

• Análise Contínua: Os autores refinam a análise contínua empregando a condição de Cordes para quantificar o efeito anisotrópico. Eles estabelecem que o problema é Fredholm de índice zero sob uma condição de pequena magnitude melhorada no parâmetro de anisotropia $\beta$. Crucialmente, esta condição é independente da constante de estabilidade do caso isotrópico, permitindo uma gama mais ampla de parâmetros fisicamente relevantes.
• Discretização: O método introduz incógnitas de faceta adicionais que representam a derivada normal das incógnitas de volume. Essas variáveis de faceta desacoplam as incógnitas de volume, permitindo a condensação estática para reduzir os custos computacionais. A formulação discreta utiliza um operador de levantamento (lifting operator) para definir um Laplaciano discreto e inclui um termo de estabilização para controlar os saltos na derivada normal através das interfaces dos elementos.
• Estrutura Teórica: A análise baseia-se em um framework abstrato para discretizações não conformes de operadores que são perturbações compactas de operadores coercivos. Os autores provam uma desigualdade de Gårding discreta, que reproduz a estrutura de "coercividade até a compacidade" do problema contínuo no nível discreto. Esta estrutura, combinada com consistência assintótica e um argumento de compacidade discreta, garante estabilidade e convergência.

Principais Contribuições

1. Melhoria da Análise de Bem-Posticidade: O artigo prova que a equação de Helmholtz–Korteweg nemática é bem posta para uma gama significativamente mais ampla do parâmetro de anisotropia $\beta$ do que o estabelecido anteriormente. Ao utilizar a condição de Cordes, a condição de pequena magnitude derivada não degenera nas frequências ressonantes do problema isotrópico.
2. Método de Penalidade de Interior Híbrido $C^0$: Os autores introduzem uma discretização $C^0$-HIP que evita as dificuldades de implementação de elementos $C^1$-conformes (como elementos de Argyris ou HCT), enquanto requer menos graus de liberdade do que métodos de Galerkin totalmente descontínuos.
3. Provas de Estabilidade e Convergência: O método é provado como estável para qualquer grau polinomial $p \geq 2$, independente da dimensão espacial, desde que a anisotropia seja suficientemente pequena. O artigo estabelece a convergência para a solução contínua sob suposições mínimas de regularidade ($u \in H^2(\Omega)$) e deriva taxas de convergência de ordem ótima sob regularidade adicional.
4. Simulações Tridimensionais: O trabalho apresenta as primeiras simulações numéricas tridimensionais dos efeitos anisotrópicos de campos nemáticos na propagação de ondas acústicas usando esta equação específica.

Resultos

• Resultados Teóricos: A forma sesquilinear discreta satisfaz uma desigualdade de Gårding discreta, garantindo estabilidade assintótica. A sequência de soluções discretas converge para a solução contínua conforme o tamanho da malha $h \to 0$.
• Experimentos Numéricos:
  • Convergência: Testes numéricos em um disco unitário confirmam as taxas de convergência previstas de $O(h^{p-1})$ para graus polinomiais $p \in \{2, 3, 4, 5\}$.
  • Estabilidade: Experimentos verificam a estabilidade para $\beta > -\alpha$ em 2D e $\beta \in (-\alpha, 3\alpha)$ em 3D, consistente com os limites teóricos derivados da condição de Cordes.
  • Anisotropia: Simulações em 2D e 3D demonstram a natureza anisotrópica da propagação de ondas. Os resultados mostram que a velocidade da onda e a atenuação variam dependendo da orientação do diretor nemático $\mathbf{n}$, levando a padrões de espalhamento distintos em comparação com o caso isotrópico.
  • Eficiência: Uma comparação de graus de liberdade e entradas não nulas mostra que a estratégia de hibridização reduz significativamente os custos computacionais para graus polinomiais mais altos em comparação com os métodos padrão de penalidade de interior $C^0$, devido à eliminação de bolhas de interior via condensação estática.

Significância
O artigo afirma que sua principal significância reside em fornecer um framework numérico robusto e flexível para resolver a equação de Helmholtz–Korteweg nemática, que é essencial para modelar a nematoacústica. Ao desacoplar a análise de estabilidade da constante de estabilidade isotrópica, o trabalho remove uma limitação artificial no parâmetro de anisotropia, permitendo a simulação de efeitos nemáticos mais fortes. Além disso, o método $C^0$-HIP proposto supera as barreiras práticas dos métodos conformes $C^1$, tornando viáveis as simulações de alta ordem de problemas anisotrópicos de quarta ordem em três dimensões e em domínios curvos. Os exemplos numéricos servem como uma prova de conceito da capacidade do método de capturar fenômenos complexos de ondas anisotrópicas.