Asymptotic numerical hypocoercivity of the space-time discontinuous Galerkin method for the Kolmogorov equation

Resumo Técnico: Hipocoercividade Numérica Assintótica do Método de Galerkin Descontínuo Espaço-Tempo para a Equação de Kolmogorov

Enunciado do Problema
O artigo aborda a discretização numérica da clássica equação de Kolmogorov, um modelo representativo para equações do tipo cinética que exibem difusão degenerada. A equação é definida em um domínio espacial limitado $\Omega \subset \mathbb{R}^2$ e assume a forma $u_t + Lu = 0$, onde $L = -\partial_{xx} + x\partial_y$. Uma característica definidora deste operador é que a difusão ocorre apenas na direção $x$, enquanto o termo de transporte $x\partial_y$ atua na direção $y$. Apesar da ausência de difusão explícita na direção $y$, o sistema exibe "hipocoercividade", uma propriedade onde a interação entre transporte e dissipação degenerada propaga o decaimento por todo o domínio espacial, levando à convergência exponencial ao equilíbrio.

O desafio central é demonstrar que os esquemas numéricos padrão, especificamente o método de Galerkin descontínuo (dG) clássico espaço-tempo, preservam essa estrutura hipocoerciva no nível discreto. Embora a hipocoercividade seja bem compreendida no cenário contínuo (notadamente através do trabalho de Villani, Hérau, Nier e outros), sua preservação em métodos de Galerkin tem permanecido amplamente inexplorada, com resultados numéricos anteriores amplamente confinados a métodos de diferenças finitas, volumes finitos ou elementos finitos não conformes que exigem formulações especializadas e não padronizadas.

Metodologia
Os autores analisam o método dG espaço-tempo padrão aplicado à equação de Kolmogorov. A discretização espacial emprega um esquema dG upwind para o termo de transporte e um esquema dG de penalidade interior simétrica para o termo difusivo. A discretização temporal utiliza uma abordagem dG padrão com grau polinomial $q$.

O núcleo da metodologia envolve a construção de uma análise de estabilidade modificada que mimetiza o arcabouço de hipocoercividade contínua. No cenário contínuo, a hipocoercividade é frequentemente provada modificando o produto interno para incluir derivadas mistas ou termos de gradiente específicos, efetivamente criando uma norma "mais forte" que admite um hiato espectral (spectral gap).

Para alcançar um resultado semelhante no nível discreto sem alterar a forma bilinear padrão do método dG, os autores empregam uma estratégia específica de função de teste. Em vez de testar a equação discreta apenas contra $U$, eles testam contra uma função composta $V = V_1 + V_2 + V_3$:

1. $V_1 = U$: O teste de energia padrão.
2. $V_2 = \tau(xU_y + U_t)$: Um termo que mimetiza o operador de transporte, escalonado por um parâmetro $\tau$ dependente do tamanho da malha $h$ e do grau polinomial $p$.
3. $V_3 = -\text{div}(A\nabla U)$: Um termo envolvendo uma matriz modificada $A$ (dependente de $h$ e $p$) e o divergente do gradiente, projetado para recuperar o controle sobre a derivada $y$ ausente.

Crucialmente, os autores observam que, para os espaços polinomiais descontínuos usados em dG, a função $V$ permanece dentro do mesmo espaço discreto que $U$. Isso permite o uso de $V$ como uma função de teste admissível sem introduzir os termos de Petrov-Galerkin não padronizados exigidos em trabalhos anteriores que utilizavam elementos contínuos.

Contribuições Principais e Resultados

1. Hipocoercividade Discreta para dG Padrão: O artigo estabelece que o método dG espaço-tempo clássico admite uma propriedade de hipocoercividade numérica assintótica. Ao provar uma condição inf-sup de um lado em uma norma "aumentada por hipocoercividade", os autores demonstram que a solução numérica decai exponencialmente ao equilíbrio para tempos longos. Este é o primeiro resultado para um esquema de Galerkin padrão.

2. Norma Aumentada e Hiato Espectral: A análise introduz uma família de normas contendo o gradiente completo da solução numérica. Os autores provam um hiato espectral dependente da malha (desigualdade do tipo Poincaré) nesta norma. A taxa de decaimento $\kappa$ depende dos parâmetros de discretização, escalonando como $\kappa \sim h_{\min}^4 p^{-8}$. Embora a taxa de decaimento degenerre conforme a malha é refinada ($h \to 0$) ou o grau polinomial aumenta ($p \to \infty$), a estrutura garante estabilidade e decaimento para qualquer discretização fixa.

3. Estabilidade Inf-Sup: Além da hipocoercividade, o artigo prova que o método dG espaço-tempo é inf-sup estável na norma aumentada. Este resultado é mais forte do que os resultados de estabilidade existentes para métodos dG espaço-tempo para problemas de advecção-difusão, pois fornece controle sobre o gradiente completo em vez de apenas a direção da advecção.

4. Limites de Erro A Priori: Utilizando a estabilidade inf-sup estabelecida, os autores derivam limites de erro a priori ótimos para soluções suficientemente regulares na norma espaço-tempo aumentada. A estimativa de erro é ótima em relação tanto ao tamanho da malha espacial $h$ quanto ao passo de tempo $k$ (sob a condição $h \approx O(k)$).

5. Validação Numérica: Experimentos numéricos confirmam as taxas de convergência teóricas e demonstram o decaimento exponencial da solução ao equilíbrio. Os experimentos mostram que, embora a taxa de decaimento teórica dependa de $h$ e $p$, a convergência observada ao equilíbrio é robusta.

Significância e Alegações
Os autores alegam que este trabalho representa a primeira demonstração de hipocoercividade numérica (assintótica) para um método de Galerkin clássico e bem conhecido. A significância reside no fato de que o método preserva a estrutura hipocoerciva sem exigir a modificação da forma bilinear com termos adicionais de Petrov-Galerkin, que foram necessários em abordagens anteriores de elementos finitos não conformes.

O artigo enfatiza que a técnica de prova é genérica. Ao construir funções de teste que mimetizam o efeito dos produtos internos modificados usados na teoria de hipocoercividade contínua, a abordagem pode potencialmente ser estendida a outras equações do tipo cinética e problemas parabólicos degenerados. Os resultados fornecem uma base rigorosa para o uso de métodos dG espaço-tempo de alta ordem para simulações de longo prazo de problemas de difusão degenerada, garantindo que os esquemas numéricos não suprimam artificialmente os mecanismos de dissipação física inerentes ao modelo contínuo.