A posteriori error bounds for finite element approximations of time-dependent mean field games

Resumo Técnico: Limites de Erro A Posteriori para Aproximações de Elementos Finitos de Jogos de Campo Médio Dependentes do Tempo

Formulação do Problema
Este trabalho aborda a aproximação numérica de sistemas de Jogos de Campo Médio (MFG) de segunda ordem e dependentes do tempo. O sistema consiste em um par acoplado de equações diferenciais parciais: uma equação de Hamilton–Jacobi–Bellman (HJB) reversa para a função de valor $u$ e uma equação de Kolmogorov–Fokker–Planck (KFP) direta para a densidade dos jogadores $m$. O sistema é definido em um domínio limitado $\Omega \subset \mathbb{R}^d$ sobre um intervalo de tempo $(0, T)$, sujeito a condições iniciais para $m$ e condições terminais para $u$. As equações apresentam um coeficiente de difusão constante $\nu > 0$, um Hamiltoniano $H$ que é $C^{1,1}$ em relação à variável do gradiente, e operadores de acoplamento $F$ e $S$ que são Lipschitz contínuos e monótonos. Um requisito crítico para a validade física da solução é a preservação da não negatividade da densidade $m$ no nível discreto.

Metodologia
Os autores analisam uma classe geral de métodos de elementos finitos (FEM) estabilizados para este sistema. A discretização espacial emprega elementos finitos afins contínuos por partes, enquanto a discretização temporal utiliza um esquema de Euler implícito. Para garantir a não negatividade da densidade e a estabilidade do esquema, termos de estabilização abstratos são introduzidos. Estas estabilizações podem incluir o agrupamento de massa (mass lumping) para derivadas temporais e estabilizações espaciais para termos advectivos.

O núcleo da metodologia envolve uma análise rigorosa de erro a posteriori. Os autores primeiro estabelecem uma equivalência no nível contínuo entre a norma do erro e a norma dual do resíduo. Eles então estendem esta análise para o cenário discreto, o que apresenta desafios específicos devido à falta de conformidade temporal nos espaços padrão do tipo Discontinuous Galerkin utilizados para a densidade. Para preencher esta lacuna, os autores definem normas estendidas em espaços de soma ($V_i + Y$) e utilizam operadores de reconstrução temporal ($I_i$) para mapear as soluções discretas para o espaço de solução contínua $Y$.

O estimador de erro $\eta$ é construído a partir de vários componentes:

1. Estimadores de Resíduo: Baseados nos resíduos de volume e nos resíduos de salto espacial das equações HJB e KFP.
2. Estimadores de Salto Temporal: Medindo a falta de conformidade entre os passos de tempo.
3. Estimadores de Estabilização: Quantificando o erro introduzido pelos termos de estabilização.
4. Estimadores de Fronteira: Contabilizando erros nas condições de tempo inicial e final.
5. Termos de Oscilação de Dados: Representando o erro de aproximação dos dados do problema.

Principais Contribuições e Resultados
O artigo apresenta três resultados teóricos primários:

1. Equivalência de Normas de Erro e Resíduo: Os autores provam que, para o problema contínuo, o erro na solução (medido em normas apropriadas de Bochner–Sobolev) é equivalente à norma dual do resíduo. Crucialmente, esta equivalência mantém-se globalmente para qualquer par $(u, m)$ com densidade não negativa, não apenas em uma vizinhança da solução exata. Isto estende resultados anteriores de problemas de estado estacionário para o caso dependente do tempo.

2. Confiabilidade e Eficiência de Estimadores Gerais: Para a classe geral de FEM estabilizados, os autores derivam limites de erro a posteriori mostrando que o erro total é limitado superiormente (confiabilidade) e inferiormente (eficiência) pelo estimador proposto $\eta$, até termos de oscilação de dados.

   • Confiabilidade: O erro é limitado pelo estimador mais a oscilação temporal de dados.
   • Eficiência: O estimador é limitado pelo erro mais a oscilação espacial de dados.
   • Uma contribuição técnica significativa é o tratamento da natureza não conforme da densidade discreta no tempo, alcançado através dos operadores de reconstrução e normas estendidas.

3. Limites Livres de Estabilização para Esquemas Estruturados: O artigo identifica uma subclasse específica de estabilizações (envolvendo mass lumping e estabilizações espaciais que preservam a afinidade por partes) onde os estimadores de estabilização podem ser limitados pelo resíduo e pelos estimadores de salto temporal. Sob a condição de que o tamanho do passo de tempo $\tau$ não seja excessivamente pequeno em relação ao quadrado do tamanho da malha ($h^2 \lesssim \tau$), os autores mostram que os estimadores de estabilização não precisam ser computados explicitamente. Isto resulta em um estimador localmente computável e localmente eficiente que depende apenas de termos padrão de resíduo e salto.

Significância
O artigo afirma fornecer os primeiros limites de erro a posteriori para uma classe geral de aproximações de elementos finitos estabilizados de sistemas de MFG dependentes do tempo. Embora trabalhos anteriores (ex: [25]) tenham estabelecido tais limites para MFGs de estado estacionário, este trabalho estende a teoria para o cenário dependente do tempo, abordando as complexidades adicionais da não conformidade temporal e do acoplamento de equações parabólicas de avanço e retrocesso.

A significância reside em permitir estratégias de refinamento de malha adaptativo para MFGs. Ao fornecer estimadores de erro confiáveis e eficientes, o trabalho permite o desenvolvimento de algoritmos que podem refinar localmente a malha e os passos de tempo com base no erro computado, oferecendo potencialmente economias computacionais significativas em comparação com o refinamento uniforme. A derivação de limites livres de estabilização para esquemas de estabilização especificamente relevantes na prática aumenta ainda mais a utilidade dos resultados, pois remove a necessidade de computar indicadores de erro complexos e específicos de estabilização na prática.