Embedded Trefftz DG method for steady Navier-Stokes flow. Part I: Oseen linearization

Resumo Técnico: Método DG de Trefftz Incorporado para Escoamento de Navier–Stokes Estacionário (Parte I: Linearização de Oseen)

Enunciado do Problema
Este artigo aborda a aproximação numérica do problema de Oseen estacionário, que serve como a linearização das equações de Navier–Stokes incompressíveis estacionárias em torno de um campo de convecção prescrito $w$. As equações governantes buscam uma velocidade $u$ e uma pressão $p$ que satisfaçam:
$$-\nu\Delta u + (w \cdot \nabla)u + \nabla p = f, \quad \text{div } u = 0$$
em um domínio limitado $\Omega \subset \mathbb{R}^d$ ($d=2,3$), com condições de contorno de Dirichlet homogêneas. O principal desafio abordado é o desenvolvimento de um método discontinuo de Galerkin (DG) de Trefftz incorporado que mantém alta qualidade de aproximação com redução de graus de liberdade (DoFs) enquanto lida com os coeficientes variáveis introduzidos pelo termo de convecção $(w \cdot \nabla)u$. Diferente do problema de Stokes (onde $w=0$), o operador de Oseen envolve coeficientes variáveis, tornando insuficientes as construções de Trefftz de coeficientes constantes padrão.

Metodologia
Os autores propõem uma abordagem DG de Trefftz incorporada que evita a construção explícita de bases de Trefftz locais. Em vez disso, o método parte de um espaço DG polinomial padrão $X$ e impõe a EDP localmente através de restrições de projeção.

1. Operador Local e Espaço de Trefftz:
   Um operador local $K_{T,w}$ é definido em cada elemento $T$, mapeando o espaço da solução para um espaço de teste dual $Q(T)'$. O espaço de Trefftz $T_w$ é definido como o subespaço de $X$ onde o resíduo local se anula fracamente contra $Q(T)$. Especificamente, as funções em $T_w$ satisfazem a EDP apenas quando testadas contra polinômios de ordem inferior (grau $k-2$ para velocidade, $k-1$ para pressão).

2. Divisão Local-Global e Espaços Complementares:
   Um passo teórico crucial é a construção de um espaço complementar local $L(T)$ tal que $X(T) = T_w(T) \oplus L(T)$.

   • Componente de Velocidade ($L_u$): Construído usando funções bolha que se anulam em uma sub-bola, enriquecidas com termos polinomiais específicos para controlar os fluxos de faceta, e restritas para garantir fluxo médio zero em $d$ facetas.
   • Componente de Pressão ($L_p$): Definido como polinômios de média zero de grau $k-1$.
   • Invertibilidade: O artigo prova que o operador local $K_{T,w}$ é limitadamente invertível em $L(T)$, desde que o campo de convecção $w$ satisfaça uma condição de resolução (especificamente, uma condição de pequenez em uma norma discreta $|w|_{h,d}$ relativa à viscosidade). Esta invertibilidade permite que a parte local da solução seja computada independentemente em cada elemento.

3. Formulação Reduzida:
   Ao eliminar os componentes de pressão de ordem superior localmente (que estão ligados à velocidade via o operador local), os autores derivam uma formulação reduzida posta unicamente em termos da incógnita de velocidade. Esta formulação:

   • Remove a estrutura de ponto de sela ao nível local.
   • Transfere a coercividade do cenário de Stokes para os espaços de Trefftz-Oseen reduzidos.
   • Retém apenas o acoplamento global para o componente de pressão constante por partes para impor a restrição de divergência discreta.

4. Análise de Estabilidade:
   A análise estabelece:

   • Continuidade das formas bilineares.
   • Coercividade da forma reduzida no subespaço de Trefftz.
   • Estabilidade inf-sup para o sistema completo acoplado em normas DG padrão.
     Estes resultados dependem do arcabouço abstrato local-global e da construção específica de $L(T)$, garantindo bem-posto e quase-otimalidade.

Principais Contribuições
O artigo estende a metodologia DG de Trefftz incorporada do problema de Stokes para as equações de Oseen estacionárias, fornecendo a base linear para a subsequente análise de Navier–Stokes não linear (Parte II). As contribuições específicas são:

• Construção de Espaços Complementares Locais: Ao contrário do caso de Stokes, que dependia de argumentos de polinômios de Taylor de coeficientes constantes, este trabalho constrói $L(T)$ para lidar com o operador de coeficientes variáveis local resultante do termo de convecção.
• Tratamento de Termos Laterais Inhomogêneos: O método incorpora termos inhomogêneos diretamente na formulação do problema local em vez de usar um argumento de correção, ajustando-se naturalmente ao arcabouço de divisão local-global abstrato.
• Formulação Reduzida e Estabilidade: Os autores derivam uma formulação reduzida apenas para a velocidade e provam a estabilidade inf-sup reduzida para o operador de Oseen completo. Isso elimina a estrutura de ponto de sela localmente e permite a transferência de coercividade sob uma condição de pequenez adequada no campo de convecção.

Resultados

• Teóricos: O artigo prova estabilidade completa e quase-otimalidade em normas DG padrão. Estabelece que o método é bem-posto e que o erro é limitado pelo melhor erro de aproximação no espaço discreto, escalonado por constantes que dependem da resolução do campo de convecção.
• Numéricos:
  • Teste de Oseen–Kovasznay: Experimentos numéricos no escoamento de Oseen–Kovasznay prescrito demonstram que o método DG de Trefftz alcança as taxas de convergência assintótica esperadas ($O(h^{k+1})$ para velocidade, $O(h^k)$ para pressão).
  • Eficiência: Embora o método DG padrão produza constantes de erro ligeiramente menores para um grau polinomial fixo, o método DG de Trefftz alcança um determinado nível de erro com significativamente menos graus de liberdade, particularmente para ordens polinomiais mais altas.
  • Custo Computacional: Resultados de tempo indicam que o método DG de Trefftz é consistentemente mais rápido que o método DG padrão, com a vantagem aumentando em malhas mais finas e para ordens mais altas.
  • Limitações: O método mostra sensibilidade no regime dominado pela convecção (baixa viscosidade), onde os erros deterioram mais cedo do que no DG padrão. Isso é atribuído à resolução e ao condicionamento da construção de Trefftz incorporada local.

Significância
O artigo afirma fornecer a base analítica linear necessária para a teoria de Navier–Stokes não linear desenvolvida na Parte II. Ao estender com sucesso a abordagem de Trefftz incorporada para o problema de Oseen com coeficientes variáveis, os autores demonstram que as técnicas de Trefftz podem ser aplicadas a operadores de coeficientes variáveis e inhomogêneos sem a construção explícita de bases. A formulação reduzida é destacada como um passo crucial para analisar o problema não linear, pois desacopla a pressão local do sistema global, simplificando a análise de estabilidade para solvers iterativos. O trabalho posiciona o problema de Oseen como uma ponte entre as discretizações clássicas de DG e a redução de modelo do tipo Trefftz, oferecendo um caminho para simulações de alto nível e baixo DoF para escoamentos viscosos incompressíveis.