Time integration as filtering: a space-time discretization-aware LES formulation

Resumo Técnico: Integração Temporal como Filtragem em LES Ciente da Discretização Espaço-Temporal

Definição do Problema
Na Simulação de Grandes Escalas Implícita (ILES), o filtro espacial é frequentemente identificado com a própria malha computacional. As formulações tradicionais de LES "cientes da discretização" abordam o problema do fechamento da turbulência derivando uma tensão residual exata baseada na propriedade de "troca de filtro" (filter-swap): a derivada central de diferenças finitas grosseira é reconhecida como uma derivada espacial exata filtrada por um filtro top-hat. No entanto, as formulações existentes tratam a integração temporal como uma etapa separada de pós-processamento. Consequentemente, o erro de discretização temporal situa-se fora da garantia de tensão residual exata. À medida que os passos de tempo aumentam (números de CFL mais altos), os fechamentos padrão baseados apenas no espaço degradam-se, pois falham em contabilizar o erro específico introduzido pelo integrador temporal, particularmente ao utilizar esquemas explícitos como o Euler Progressivo (Forward Euler).

Metodologia
O artigo estende o conceito de troca de filtro espacial para o domínio espaço-temporal, observando que a diferença temporal do método Euler Progressivo é, ela própria, uma derivada temporal exata filtrada (do tipo top-hat, unilateral).

1. Troca de Filtro Espaço-Temporal: Ao aplicar um filtro espaço-temporal $(\cdot)_{h,\tau} = (\cdot)_\tau \circ (\cdot)_h$ à lei de conservação escalar, os autores derivam uma lei discreta exata onde ambas as derivadas espaciais e temporais são substituídas por suas contrapartes filtradas.
2. Decomposição de Fluxo: O fluxo espaço-temporal exato é decomposto em fluxo numérico mais um fluxo residual, $r_{h,\tau}$. Este resíduo é dividido em quatro termos de erro distintos:
   • $r^{LES}_h$: O clássico comutador de LES (filtragem de malha vs. não linearidade).
   • $r^{num}_h$: O erro de fluxo numérico.
   • $r^{div}_h$: Uma correção para a divergência discreta.
   • $r^{time}_\tau$: Um novo termo que representa o erro de quadratura temporal do fluxo (a diferença entre o fluxo médio no passo de tempo e o fluxo instantâneo).
3. Modelagem: Os autores propõem um modelo de viscosidade de eddy ciente da discretização espaço-temporal. Eles expandem o termo de erro de quadratura temporal, $r^{time}_\tau$, demonstrando que este se comporta como um fluxo difusivo proporcional a $\tau (f'(u))^2 \partial_x u$. Este termo é identificado como uma difusão do tipo Lax–Wendroff, que é não negativa por construção e escala linearmente com o passo de tempo $\tau$. O fechamento final aumenta o modelo de Smagorinsky padrão com este termo temporal.

Resultados Principais
O arcabouço foi testado em turbulência de Burgers viscosa utilizando um esquema de volumes finitos Euler Progressivo com números de CFL variando de 0,06 a 0,89.

• Análise A-priori: Na decomposição do fluxo residual, o termo temporal ($r^{time}_\tau$) cresce significativamente com o número de CFL, tornando-se o maior termo de erro individual em valores de CFL práticos (ex: 0,59 do resíduo em CFL 0,89 para fluxos centrais). Por outro lado, os termos espaciais mantêm proporções fixas, independentemente do passo de tempo.
• Análise A-posteriori:
  • Fluxo Central: Os fechamentos padrão (sem modelo, Smagorinsky clássico e Smagorinsky apenas espacial) degradam-se rapidamente conforme o número de CFL aumenta, com erros subindo de 0,1 para ~0,85. O fechamento ciente do espaço-tempo proposto mantém um erro relativo baixo e constante (0,09) em toda a gama de números de CFL.
  • Fluxo Upwind: O fluxo numérico upwind é inerentemente excessivamente dissipativo. Neste caso, o erro temporal (que não é dissipativo) cancela parcialmente a sobre-dissipação do fluxo upwind. Embora o modelo espaço-temporal permaneça preciso, ele é ligeiramente excessivamente dissipativo em passos grandes, pois tenta ajustar os efeitos de cancelamento de forma independente, em vez de explorar o cancelamento fortuito.

Significância e Alegações
O artigo alega que reconhecer a diferença temporal do Euler Progressivo como um filtro top-hat permite uma decomposição de fluxo discreta exata em espaço-tempo. A principal contribuição é a identificação do erro de quadratura temporal do fluxo como um termo distinto e dominante no resíduo em passos de tempo grosseiros.

• Os autores demonstram que um fechamento aumentado com o termo de ordem principal deste erro (difusão Lax–Wendroff) permite que simulações de Euler Progressivo permaneçam precisas mesmo em passos de tempo largos, onde fechamentos baseados apenas no espaço falham.
• O trabalho estabelece que o integrador temporal contribui com exatamente um novo termo para o resíduo (o erro de quadratura) e nenhum termo de comutação temporal separado.
• O artigo nota modestamente que grande parte do erro de quadratura temporal é não dissipativo e, portanto, invisível para modelos padrão de viscosidade de eddy, sugerindo que fechamentos aprendidos por máquina (machine-learned) treinados no alvo exato ou em integradores implícitos/multi-estágios são extensões naturais para o futuro.