A three-field Multiscale Method

Resumo Técnico: Um Método Multiescala de Três Campos (MH2M)

Definição do Problema
O artigo aborda a simulação numérica do fluxo de fluidos em meios porosos heterogêneos, modelados pela equação de Darcy (ou, equivalentemente, pela equação de Poisson com um coeficiente multiescala). O desafio principal reside em capturar soluções multiescala sem recorrer a malhas finas computacionalmente proibitivas. Métodos de elementos finitos mistos (utilizando elementos de Raviart-Thomas ou Brezzi-Douglas-Marini) oferecem vantagens como conservação de massa local e continuidade de fluxo, mas resultam em sistemas lineares de ponto de sela que são difíceis de resolver eficientemente. Técnicas de hibridização têm sido empregadas para convertê-los em sistemas definidos positivos, mas os métodos multiescala híbrido-mistos (MHM) existentes frequentemente retêm uma estrutura de ponto de sela ou dependem de termos de estabilização que impõem restrições à malha.

Metodologia: O Método Híbrido-Híbrido Multiescala (MH2M)
Inspirado pelo seminal método de decomposição de domínio de três campos de Brezzi e Marini, os autores propõem o Método Híbrido-Híbrido Multiescala (MH2M). O método opera em uma partição grossa do domínio e introduz uma formulação de três campos envolvendo:

1. A variável primal (pressão) $u$.
2. A variável dual (fluxo) $\lambda$.
3. Um segundo multiplicador de Lagrange $\rho$, representando o traço da pressão.

A inovação central envolve o relaxamento da continuidade da variável de fluxo $\lambda$ e a imposição fraca desta através da ação do segundo multiplicador de Lagrange $\rho$. A solução exata é decomposta em contribuições locais e globais. O problema global é definido exclusivamente no esqueleto da partição grossa (as interfaces entre os elementos), enquanto os problemas locais em cada elemento fornecem funções de base multiescala.

Características principais da formulação incluem:

• Sistema Global Simétrico Positivo Definido (SPD): Diferente do método MHM original que induz um sistema de ponto de sela, o MH2M gera um sistema linear global simétrico positivo definido para computar os graus de liberdade (o traço $\rho$).
• Base de Fluxo Multiescala: A variável de fluxo $\lambda$ é aproximada usando uma base multiescala derivada de um operador de Dirichlet-para-Neumann definido por novos problemas de Neumann locais, impulsionados pelas funções de base de $\rho$. Isso contrasta com métodos anteriores onde a base de fluxo era frequentemente escolhida de forma ad-hoc.
• Continuidade Fraca: O método impõe continuidade fraca tanto da variável primal (pressão) quanto da variável dual (fluxo) no esqueleto grosso, garantindo que o fluxo discreto esteja em equilíbrio local com as forças externas.

Principais Contribuições
O artigo estabelece as seguintes contribuições teóricas e práticas:

1. Bem-posto e Melhor Aproximação: Os autores demonstram a natureza bem-posta e as propriedades de melhor aproximação do MH2M sob condições de compatibilidade abstratas entre os espaços de interpolação.
2. Análise de Convergência: Eles propõem famílias específicas de espaços de elementos finitos (por exemplo, $P_{k+1}(F) \times P_k(F) \times P_{k+1}(\tau)$) que satisfazem as condições de compatibilidade exigidas. Taxas de convergência ótimas são provadas sob suposições de regularidade local, destacando a influência de diferentes parâmetros de malha (tamanho da malha grossa, tamanho da malha da face e tamanho da submalha local).
3. Conexão com MsFEM: O trabalho estabelece uma ponte entre o MH2M e o Método de Elementos Finitos Multiescala (MsFEM), mostrando que, sob condições específicas (operadores contínuos e espaços polinomiais específicos), o traço do MH2M coincide com o traço do MsFEM.
4. Implementação Algorítmica: O método é estruturado para permitir computação "paralelamente trivial" (embarrassingly parallel). Os problemas locais (computação das funções de base multiescala) são independentes e podem ser resolvidos offline, enquanto o problema global envolve apenas os graus de liberdade no esqueleto grosso.

Resultos Numéricos
O artigo apresenta validações numéricas comparando o MH2M com MsFEM, MHM e FEM padrão:

• Precisão vs. Graus de Liberdade: Em testes com soluções explícitas, o MH2M atinge níveis de erro comparáveis ao MHM e MsFEM, mas com significativamente menos graus de liberdade globais.
• Superconvergência: Quando as malhas das faces associadas às variáveis de traço e fluxo são refinadas enquanto a malha grossa permanece fixa, o método exibe comportamento de superconvergência (uma taxa adicional de O(H^{1/2}_\Gamma}), um fenômeno também observado no MHM, mas alcançado aqui sem aumentar o tamanho do sistema global.
• Coeficientes Oscilatórios: Em problemas com coeficientes altamente oscilatórios, o MH2M demonstra desempenho superior ao MHM. Especificamente, o MH2M evita os "efeitos de ressonância" observados no MHM quando o tamanho da malha grossa $H$ é comparável à escala de oscilação $\sqrt{\epsilon}$. Isso é atribuído ao refinamento local da malha de fluxo dentro da estrutura do MH2M, que não infla o tamanho do sistema global.

Significância
O artigo afirma que o MH2M oferece uma alternativa robusta para problemas de Darcy multiescala ao combinar os benefícios da hibridização (sistemas SPD, conservação de massa local) com uma nova decomposição espacial que elimina a necessidade de termos de estabilização frequentemente exigidos em métodos de três campos similares. Ao desacoplar a computação da base multiescala local da resolução global e garantir que o sistema global seja SPD, o método é bem adequado para arquiteturas paralelas e solvers eficientes. Os resultados sugerem que o MH2M pode alcançar alta precisão com menor custo computacional (em termos de graus de liberdade globais) em comparação com métodos multiescala existentes, particularmente em cenários envolvendo coeficientes oscilatórios onde a ressonância é uma preocupação.