Parameter-robust well-posedness and discretisation for coupled Darcy--Forchheimer and advection-diffusion-reaction equations

Resumo Técnico: Bem-posedness e discretização paramétrica robusta para equações acopladas de Darcy–Forchheimer e advecção-difusão-reação

Enunciado do Problema
Este trabalho investiga a análise matemática e a discretização numérica de um sistema totalmente acoplado que compreende as equações de Darcy–Forchheimer e a equação de advecção-difusão-reação (ADR). O sistema modela o fluxo não linear em meios porosos acoplado ao transporte de solutos, onde a velocidade do fluxo $u$ e a pressão $p$ satisfazem um balanço de momento não linear com um termo de arrasto de Forchheimer, e a concentração do soluto $\phi$ e o fluxo total $\zeta$ satisfazem uma equação de transporte. Um desafio crítico abordado é a dependência da estabilidade e convergência em relação aos parâmetros do modelo, especificamente a condutividade hidráulica (permeabilidade) $\kappa$, o coeficiente de Forchheimer $F$, o coeficiente de difusividade $\alpha$ e o coeficiente de reação $\varrho_0$. Os autores visam estabelecer a bem-posedness (bem-postura) e derivar estimativas de erro que permaneçam robustas (independentes) em relação a esses parâmetros, que podem variar em várias ordens de magnitude em aplicações práticas, como modelagem de reservatórios e otimização de energia geotérmica.

Metodologia
O artigo emprega uma abordagem de análise funcional no contexto de espaços de Banach reflexivos, indo além do padrão de espaços de Hilbert para acomodar as não linearidades e os requisitos de escala específicos do problema.

1. Teoria Abstrata: Os autores adaptam a teoria recente para problemas de ponto de sela perturbados em espaços de Banach. Eles estendem a teoria de Brezzi–Braess para estabelecer limites de estabilidade independentes de parâmetros. Isso envolve a definição de espaços de Banach ponderados apropriados para a velocidade, pressão, fluxo e concentração. As normas nesses espaços são construídas como somas e interseções de espaços ponderados pelos parâmetros do modelo (por exemplo, $\kappa^{-1/2}$, $F^{1/3}$, $\alpha^{-1/2}$).
2. Análise de Bem-posedness: Utilizando um argumento de ponto fixo, o sistema acoplado é desacoplado em um subproblema de Darcy–Forchheimer e um subproblema de ADR. Os autores provam que os subproblemas individuais são bem-postos e, em seguida, demonstram que a composição de seus operadores de solução forma um mapeamento de contração em uma bola fechada no espaço de velocidade, desde que os dados sejam suficientemente pequenos. Isso estabelece a existência e unicidade da solução contínua com constantes de estabilidade independentes de parâmetros.
3. Discretização por Elementos Finitos: O problema é discretizado utilizando métodos de elementos finitos mistos. Os autores utilizam elementos de Raviart–Thomas para os espaços de velocidade e fluxo e espaços de Galerkin descontínuos (Sobolev quebrados) para pressão e concentração. Crucialmente, os espaços discretos são equipados com as mesmas normas ponderadas por parâmetros utilizadas na análise contínua.
4. Análise de Erro: Estimativas de erro a priori são derivadas nestas normas ponderadas. A análise baseia-se em condições de inf-sup discretas e em um argumento de ponto fixo discreto para mostrar que as taxas de convergência são ótimas e robustas em relação a $\kappa$, $F$, $\alpha$ e $\varrho_0$.
5. Précondicionamento: Para abordar os desafios computacionais da resolução dos sistemas lineares não lineares e não simétricos resultantes, os autores projetam précondicionadores de operador de bloco-diagonal. Estes são inspirados em estruturas de complemento de Schur e frameworks de précondicionamento de operador, adaptados aos pesos de parâmetros específicos do sistema acoplado.

Principais Contribuições

• Extensão para Espaços de Banach: O artigo estende a teoria de Brezzi–Braess para problemas de ponto de sela perturbados de espaços de Hilbert para espaços de Banach reflexivos, abordando especificamente o contexto de robustez paramétrica.
• Estabilidade Paramétrica Robusta: Os autores provam que o sistema acoplado Darcy–Forchheimer/ADR é bem-posto e estável, com limites independentes da permeabilidade, Forchheimer, difusividade e coeficientes de reação. Isso é alcançado através da construção cuidadosa de espaços e normas de Banach ponderadas.
• Análise Discreta Robusta: O artigo fornece uma análise discreta rigorosa mostrando que a formulação de elementos finitos mistos preserva essas propriedades de estabilidade robustas. As estimativas de erro a priori derivadas são quase-ótimas e independentes dos parâmetros do modelo e do tamanho da malha $h$.
• Précondicionamento de Operador: É proposto o design de précondicionadores de bloco-diagonal variáveis. Estes précondicionadores mostram-se robustos em relação aos parâmetros do modelo (exceto o coeficiente de difusão no regime dominado pela advecção) e ao parâmetro de discretização.

Resultados Numéricos
As descobertas teóricas são validadas através de uma série de experimentos numéricos em 2D e 3D:

• Convergência: Testes com soluções suaves manufaturadas e permeabilidade heterogênea confirmam as taxas de convergência ótimas ($O(h^{k+1})$) tanto em normas padrão quanto em normas ponderadas por parâmetros. Os resultados demonstram que, enquanto as normas padrão podem apresentar sensibilidade a valores extremos de parâmetros, as normas ponderadas restauram o comportamento de convergência ótima.
• Desempenho do Précondicionador: Testes numéricos no sistema linearizado mostram que os précondicionadores propostos produzem números de condição e contagens de iteração MINRES que permanecem estáveis através de amplas variações de $\kappa$, $F$, $\alpha$ e $\varrho_0$. Os autores observam que a robustez é observada no regime dominado pela difusão; o desempenho no regime dominado pela advecção (altos números de Péclet) não é totalmente robusto, consistente com as limitações da estabilização escolhida.
• Aplicações Práticas: O método é aplicado a um problema de fluxo de cavidade de tampa (lid-driven cavity) e a um problema de deslocamento miscível transiente (instabilidade de viscosidade em uma célula de Hele–Shaw), demonstrando a capacidade do método em lidar com acoplamentos complexos de fluxo-transporte e não linearidades.

Significância e Alegações
O artigo alega que sua principal significância reside em fornecer um framework matematicamente rigoroso para a análise e simulação de fluxo e transporte acoplados em meios porosos que permanece estável e preciso, independentemente da magnitude dos parâmetros físicos. Ao deslocar a análise para espaços de Banach e utilizar normas ponderadas por parâmetros, os autores superam as limitações das abordagens padrão de espaços de Hilbert, que frequentemente falham em fornecer estabilidade uniforme para tais sistemas não lineares e acoplados. O trabalho oferece uma base teórica para o desenvolvimento de solvers numéricos eficientes e robustos para problemas onde parâmetros como permeabilidade e taxas de reação variam significativamente, um cenário comum na engenharia de subsuperfície e modelagem ambiental. Os autores reconhecem modestamente que, embora o framework de ponto fixo gere limites teóricos robustos, a robustez numérica observada no sistema acoplado sugere que a perda de robustez analítica em certas normas ponderadas pode ser um artefato da escala ou do próprio framework de ponto fixo, e não uma limitação fundamental do método.