Resumo Técnico: Advecção de Elementos Espectrais Miméticos

Declaração do Problema

O artigo aborda a discretização numérica da equação de advecção linear para formas diferenciais em domínios limitados. O problema clássico de advecção escalar, $\partial_t \rho + \nabla \cdot (v\rho) = 0$, é generalizado dentro de um arcabouço de geometria diferencial como:
$$\frac{\partial \alpha^{(k)}}{\partial t} + \mathcal{L}_v \alpha^{(k)} = 0$$
onde $\alpha^{(k)}$ é uma $k$-forma diferencial (representando escalares, vetores ou quantidades de dimensões superiores dependendo de $k$), $v$ é um campo de velocidade Lipschitz contínuo prescrito, e $\mathcal{L}_v$ é a derivada de Lie. O desafio reside em desenvolver uma discretização que preserve as estruturas geométricas subjacentes das leis físicas, garanta a conservação de massa local e alcance precisão espectral, distinguindo entre quantidades avaliadas em instantes de tempo específicos e aquelas integradas sobre intervalos de tempo.

Metodologia

O método proposto, Advecção de Elementos Espectrais Miméticos, estende o arcabouço mimético estabelecido em trabalhos anteriores [4] incorporando a derivada de Lie via fórmula de homotopia de Cartan. A abordagem baseia-se em três pilares centrais:

1. Fundamentação de Geometria Diferencial
O método utiliza a linguagem das formas diferenciais. A derivada de Lie $\mathcal{L}_v$ é expressa usando a fórmula de homotopia de Cartan:
$$\mathcal{L}_v \alpha^{(k)} = d\iota_v \alpha^{(k)} + \iota_v d\alpha^{(k)}$$
onde $d$ é a derivada exterior e $\iota_v$ é o produto interno. Um aspecto crucial da formulação é a relação de adjunção entre o produto interno e o produto exterior (wedge product), definida via o produto interno $L^2$:
$$(\iota_v \alpha^{(k)}, \beta^{(k-1)})_{L^2} = (\alpha^{(k)}, v^\flat \wedge \beta^{(k-1)})_{L^2}$$
Esta dualidade permite que quantidades físicas representadas por produtos internos sejam mapeadas para suas formas diferenciais duais de ordem 1.

2. Discretização Espacial (Elementos Espectrais Miméticos)
O domínio espacial é discretizado usando um complexo de células composto por pontos, segmentos de linha, superfícies e volumes. O espaço de $k$-formas suaves é aproximado por um espaço de dimensão finita $\Lambda_h^k$ gerado por formas de base $\epsilon_i^{(k)}$.

• Construção da Base: As funções de base são construídas usando expansões polinomiais por partes (produtos tensoriais) em elementos quadriláteros. Em 1D, as 0-formas são representadas por polinômios de Lagrange avaliados em nós de Gauss-Lobatto, enquanto as 1-formas são representadas por "polinômios de aresta" derivados para satisfazer propriedades integrais específicas.
• Projeção Comutativa: Um operador de projeção $\pi_h$ é definido de tal forma que comuta com a derivada exterior ($\pi_h d = d \pi_h$). Isso garante que a derivada exterior discreta $d_h$ seja representada por matrizes de incidência contendo apenas valores $\{0, 1, -1\}$, preservando a estrutura topológica dos operadores contínuos.

3. Discretização Temporal (Integrador Mimético Canônico)
A integração temporal emprega um integrador de ordem arbitrária simplético derivado de integradores de colocação de Gauss canônicos [5].

• Nós Temporais Deslocados (Staggered): O método distingue entre dois conjuntos de nós temporais:
  • Nós de Gauss-Lobatto ($t_k$): Onde as variáveis da solução (ex: $\rho$) são definidas.
  • Nós de Gauss ($\tilde{t}_q$): Onde as derivadas temporais (fluxos) são avaliadas.
• Evolução Discreta: Este deslocamento resulta em um integrador discreto onde a mudança na solução ao longo de um intervalo é igualada ao fluxo avaliado nos nós internos de Gauss. Esta estrutura espelha métodos leap-frog e a regra do ponto médio implícito, garantindo propriedades simpléticas.

4. Produto Interno Discreto
O produto interno $\iota_v$ é discretizado impondo a dualidade de emparelhamento (Eq. 16) no cenário discreto. Isso leva a um sistema onde os fluxos discretos são computados resolvendo um sistema linear envolvendo os produtos internos das formas de base e o campo de velocidade.

Principais Contribuições

• Consistência Geométrica: O esquema incorpora explicitamente a derivada de Lie usando a fórmula de Cartan, garantindo que a discretização respeite a representação livre de métrica dos operadores diferenciais e seus adjuntos de Hilbert.
• Precisão Espectral: O método utiliza funções de base polinomiais de alta ordem, permitindo convergência espectral no espaço.
• Conservação de Massa Local: O uso de matrizes de incidência e a estrutura específica da derivada exterior discreta garantem a conservação de massa local.
• Distinção de Passo de Tempo: O arcabouço separa rigorosamente as quantidades avaliadas em instantes de tempo daquelas integradas sobre intervalos, utilizando uma grade temporal deslocada para alcançar integração simplética de alta ordem.

Resultados Numéricos

Os autores apresentam experimentos numéricos em domínios 2D com condições de contorno periódicas, testando a advecção de ondas senoidais e "sine bells" em campos de velocidade constantes e de vórtice de Rudman.

• Convergência: O método demonstra convergência algébrica $h$ de ordem $(p+1)$ e convergência espectral $p$, desde que o erro de integração temporal não domine o erro espacial.
• Precisão da Integração Temporal: Mostra-se que o erro na solução depende da ordem do esquema de integração temporal ($p_t$). Se o esquema de tempo for suficientemente preciso, o erro inicial de discretização é preservado; caso contrário, a dispersão artificial aumenta com o tempo.
• Conservação de Massa: O erro de massa total permanece em zero de máquina durante os primeiros $10^3$ passos de tempo e mantém-se abaixo de $10^{-12}$ mesmo após $2 \times 10^4$ passos, demonstrando excelentes propriedades de conservação.
• Reversibilidade: Um teste envolvendo a advecção de uma onda senoidal em um vórtice de Rudman, seguido pela reversão da direção do fluxo, mostra que o método pode recuperar a solução inicial, provando a reversibilidade do método de integração.
• Dispersão Artificial: O artigo observa que a dispersão artificial é uma função da ordem da integração temporal; esquemas de tempo de ordem inferior introduzem erros de dispersão que dependem da frequência da onda advectada.

Significância e Alegações

O artigo alega que o esquema derivado combina com sucesso precisão espectral com conservação de massa local dentro de um arcabouço compatível com a física (mimético). Ao estender o arcabôço mimético para incluir a derivada de Lie via fórmula de homotopia de Cartan, os autores fornecem uma discretização que esclarece a estrutura geométrica da equação de advecção. O método é apresentado como uma abordagem robusta para resolver problemas de advecção onde preservar as propriedades geométricas e topológicas subjacentes das leis físicas é essencial. Os autores notam modestamente que, embora a discretização espacial seja altamente precisa, a precisão geral e as características de dispersão dependem da ordem do esquema de integração temporal utilizado.