A HDG method with transmission variables for time-harmonic wave propagation problems with constant coefficients

Resumo Técnico: Um Método HDG com Variáveis de Transmissão para Propagação de Ondas em Regime Harmônico no Tempo

Definição do Problema
Solucionadores iterativos de elementos finitos para problemas de propagação de ondas em regime harmônico no tempo são notoriamente lentos para convergir devido às propriedades intrínsecas dos sistemas lineares subjacentes. Embora os métodos de Discontinuous Galerkin (DG) sejam bem adequados para soluções de alta fidelidade em geometrias complexas e regimes de alta frequência, os sistemas lineares esparsos e não estruturados resultantes impõem desafios significativos para solucionadores diretos em arquiteturas paralelas. Embora os solucionadores iterativos ofereçam eficiência de memória, suas taxas de convergência são frequentemente insuficientes. Métodos Hybridizable Discontinuous Galerkin (HDG) existentes melhoram o condicionamento do sistema ao eliminar incógnitas físicas dentro dos elementos, restando um sistema menor definido no esqueleto da malha. No entanto, formulações HDG padrão, que utilizam traços numéricos como variáveis híbridas, nem sempre proporcionam a convergência mais rápida para esquemas iterativos em contextos de propagação de ondas.

Metodologia
Os autores propõem uma variante do método HDG, denominada CHDG (HDG baseado em características), que substitui o traço numérico padrão por variáveis de transmissão derivadas da análise de características dos sistemas hiperbólicos governantes.

1. Estrutura Matemática: O método visa sistemas hiperbólicos simétricos lineares com coeficientes constantes no domínio da frequência. As equações governantes são analisadas via decomposição de características, dividindo a solução em variáveis de transmissão de saída ($G^+$) e de entrada ($G^-$) com base nos autovalores da matriz de fluxo.
2. Formulação DG: Um esquema DG é construído utilizando fluxos upwind expressos explicitamente em termos dessas variáveis de transmissão. O fluxo numérico em uma interface utiliza a variável de saída do elemento atual e a variável de entrada do elemento vizinho.
3. Hibridização: A inovação central reside na estratégia de hibridização. Em vez de eliminar variáveis físicas via um traço numérico, o método introduz as variáveis de transmissão de entrada ($G^-$) como as incógnitas globais no esqueleto da malha.
   • As variáveis físicas dentro de cada elemento são eliminadas ao resolver problemas locais onde $G^-$ é tratado como um dado de contorno.
   • Isso resulta em um sistema híbrido global da forma $(I - \Pi S)G^- = b$, onde $S$ é um operador de espalhamento (mapeando de entrada para saída) e $\Pi$ é um operador de troca (lidando com o acoplamento de interface e condições de contorno).
4. Propriedades Teóricas: Os autores provam que, para condições de contorno passivas, o operador $\Pi S$ é uma contração estrita. Consequentemente, o sistema hibridizado é bem posto e suscetível à iteração de ponto fixo.

Aplicação em Aeroacústica
O arcabouço é aplicado especificamente às equações de Euler linearizadas (LEE) no tempo harmônico, que modelam a propagação de ondas em um escoamento médio subsônico uniforme. Esta aplicação introduz complexidades específicas:

• O sistema suporta tanto ondas sonoras quanto ondas de vorticidade.
• O número de variáveis de transmissão de entrada e saída em uma face de elemento não é constante; depende da orientação da face em relação à velocidade do escoamento médio ($u_0 \cdot n$).
• O arcabouço geral proposto lida naturalmente com essa variabilidade no número de condições de contorno e variáveis de transmissão através do domínio.

Resultos Numéricos
O método é validado através de benchmarks numéricos bidimensionais envolvendo ondas sonoras e de vorticidade em escoamentos médios uniformes, comparando a convergência de três esquemas iterativos (Ponto Fixo, CGNR e GMRES) aplicados tanto a sistemas DG padrão quanto a sistemas CHDG:

• Propagação de Onda Plana: Testes com variações de direções de propagação em relação ao escoamento médio e diferentes condições de contorno (impedância vs. pressão).
• Fonte Pontual: Um problema de fonte singular em um domínio circular com números de Mach variados.
• Ondas de Vorticidade: Propagação em um guia de ondas rígido com modos de ordem superior.

Principais Descobertas:

• A iteração de Ponto Fixo aplicada ao sistema CHDG converge consistentemente, alinhando-se com a prova de contração teórica.
• Esquemas iterativos (CGNR e GMRES) aplicados ao sistema CHDG convergem significativamente mais rápido (menos iterações) do que quando aplicados aos sistemas DG padrão.
• No contexto CHDG, o método de Ponto Fixo performa quase tão bem quanto o GMRES e supera consistentemente o CGNR.
• A presença de ondas de vorticidade e orientações de escoamento variáveis não prejudica as propriedades de convergência da abordagem CHDG.

Significância e Alegações
O artigo alega estender a metodologia CHDG, previamente estabelecida para problemas de Helmholtz, para um arcabouço geral que cobre uma ampla classe de problemas de propagação de ondas definidos por sistemas hiperbólicos simétricos com coeficientes constantes.

• Tratamento Unificado: O arcabouço unifica o tratamento de diversos tipos de ondas (acústicas e de vorticidade) dentro de uma única estrutura matemática.
• Robustez: Gerencia efetivamente a complexidade do número variável de variáveis de transmissão causado pela orientação do escoamento médio, uma característica não tratada naturalmente por implementações específicas anteriores.
• Eficiência: A reformulação acelera significativamente a convergência de solucionadores iterativos em comparação com as abordagens DG e HDG padrão, tornando-a um candidato promissor para simulações de propagação de ondas de alta frequência e grande escala.

Os autores observam que, embora os resultados atuais sejam baseados em scripts MATLAB, trabalhos preliminares em ondas eletromagnéticas sugerem benefícios semelhantes, e trabalhos futuros visam estender este arcabouço para ondas elásticas e meios não homogêneos.