Embedded Trefftz DG method for reaction-diffusion problems on anisotropic meshes

Resumo Técnico: Método DG de Trefftz Incorporado para Problemas de Reação–Difusão em Malhas Anisotrópicas

Enunciado do Problema
O artigo aborda a solução numérica de problemas de reação–difusão singularmente perturbados da forma $-\varepsilon^2 \Delta u + u = f$ em um domínio limitado $\Omega \subset \mathbb{R}^2$, com condições de contorno de Dirichlet. O foco está no regime onde $0 < \varepsilon \ll 1$, um cenário onde as soluções exibem camadas de contorno agudas e características fortemente anisotrópicas. Malhas isotrópicas padrão são frequentemente ineficientes para esses problemas, necessitando do uso de refinamento de malha anisotrópica (elementos quadrilaterais esticados em direções específicas) para resolver características localizadas com menos graus de liberdade. Os métodos de Discontinuous Galerkin (DG) são flexíveis para discretizações de alta ordem em malhas gerais, mas as formulações DG de produto tensorial padrão resultam em sistemas algébricos globais extensos. Métodos de Trefftz, que utilizam espaços de aproximação locais contidos no núcleo do operador diferencial, oferecem uma redução nos graus de liberdade (DoFs); no entanto, construir espaços de Trefftz adequados para PDEs elípticas gerais é não trivial, particularmente em elementos quadrilaterais anisotrópicos e curvos.

Metodologia
Os autores propõem um método Discontinuous Galerkin (DG) de Trefftz Incorporado adaptado para malhas quadrilaterais anisotrópicas. O método é construído através dos seguintes passos:

1. Malha e Espaços: O domínio é particionado em elementos quadrilaterais curvilíneos. Cada elemento $K$ é mapeado de um quadrado unitário via uma transformação afim (lidando com a anisotropia) e um difeomorfismo $C^1$ (lidando com a curvatura). O espaço discreto consiste em polinômios de produto tensorial nas coordenadas auxiliares.
2. Construção de Trefftz Incorporada: Em vez de construir explicitamente o espaço de Trefftz global, o método emprega uma estratégia de incorporação. Um espaço de teste local $L_h(K)$ é definido em cada elemento. Este espaço consiste em polinômios de produto tensorial que se anulam nas duas arestas mais longas do elemento de pullback.
3. Operador Local e Restrições: Um operador de reação–difusão local $A_K$ é definido. O espaço de Trefftz incorporado $T_h(K)$ é definido como o subespaço do espaço DG padrão onde o resíduo da PDE, testado contra funções em $L_h(K)$, é nulo elemento a elemento: $(A_K v_h, q_h)_K = 0$ para todo $q_h \in L_h(K)$.
4. Formulação Global: O método global busca uma solução $u_h$ satisfazendo duas condições:
   • Restrições locais: $(A_K u_h, v_L)_K = \ell_K(v_L)$ para todas as funções de teste locais $v_L \in L_h(K)$.
   • Consistência global: A forma bilinear DG de penalidade interior padrão $a_h(u_h, v_T) = \ell_h(v_T)$ é válida para todas as funções de teste globais $v_T$ no espaço de Trefftz.
5. Implementação: Na prática, o espaço de Trefftz não é explicitamente montado. Em vez disso, o método utiliza a Decomposição de Valores Singulares (SVD) das matrizes do operador local para construir uma incorporação, retendo apenas os DoFs globais associados ao espaço de Trefftz.

Principais Contribuições
O artigo apresenta várias contribuições técnicas específicas:

• Design de Espaços de Teste Locais: A principal novidade é o design de um espaço de teste local compatível com bases de polinômios de produto tensorial. Diferente de trabalhos anteriores que utilizavam bases simpliciais, esta abordagem permite a construção sistemática em malhas quadrilaterais e hexaédricas.
• Estabilidade em Elementos Anisotrópicos: Os autores provam estabilidade e quase-otimalidade em elementos quadrilaterais anisotrópicos, possivelmente curvos. Um aspecto crucial é a escolha específica do espaço de teste (anulando-se nas arestas mais longas), o que permite a derivação de desigualdades inversas anisotrópicas e desigualdades de Poincaré–Steklov essenciais para a estabilidade.
• Análise de Erro: O artigo deriva estimativas de erro a priori anisotrópicas na norma de energia. Estas estimativas capturam o comportamento da solução nas direções reais de anisotropia, fornecendo limites mais agudos do que as estimativas isotrópicas.
• Redução de Dimensão: O método reduz a dimensão do espaço de aproximação local de $(p+1)^2$ (produto tensorial padrão) para $2p+2$, diminuindo significativamente o número de graus de liberdade globalmente acoplados em comparação com os métodos DG e HDG padrão, particularmente em 3D.

Resultados
Experimentos numéricos validam as descobertas teóricas:

• Convergência $h$: Em domínios quadrados e circulares com camadas de contorno, o método de Trefftz incorporado atinge taxas de convergência ótimas na norma de energia. Ele requer significativamente menos graus de liberdade para atingir a mesma precisão que os métodos DG padrão.
• Domínios Curvos: O método opera de forma robusta em malhas curvas, incluindo casos de elementos mistos triangulares e quadrilaterais (onde os elementos triangulares são tratados via técnicas de incorporação padrão).
• Convergência $hp$: Em uma malha de nove elementos adaptada a camadas de contorno, o método de Trefftz incorporado exibe convergência exponencial. No entanto, os autores observam que para graus polinomiais elevados ($p \gtrsim 5$) e grandes tamanhos de camada de contorno, a incorporação padrão pode tornar-se instável. Isso é mitigado modificando o espaço de teste nos elementos de canto para relaxar a condição de Trefftz, restaurando a estabilidade e a convergência exponencial.
• Difusão Anisotrópica: O método é testado em um problema de difusão puramente anisotrópica (sem termo de reação), que está fora do framework teórico estrito. Os resultados mostram que o método de Trefftz incorporado consistentemente produz erros menores do que o DG padrão sob estratégias de refinamento $h$ e $p$.

Significância e Alegações
O artigo afirma fornecer um framework geral para métodos DG do tipo Trefftz que permanece eficaz mesmo quando o núcleo da PDE ou o seu kernel é não conhecido explicitamente ou é trivial. Ao estender com sucesso a abordagem de Trefftz incorporado para malhas quadrilaterais anisotrópicas, este trabalho abre as portas para discretizações de alta ordem em malhas hexaédricas para problemas singularmente perturbados. Os autores enfatizam que o método oferece uma alternativa competitiva ao HDG, gerando os menores sistemas globais em três dimensões enquanto mantém as propriedades de aproximação de discretizações de alta ordem. A análise é modesta quanto à robustez em $p$, reconhecendo que a estabilidade para graus polinomiais muito elevados em malhas grosseiras com grandes camadas pode exigir modificações específicas no espaço de teste, conforme demonstrado nos experimentos numéricos.