Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Resumo Técnico: Método de elementos virtuais de função de corrente – pressão para o problema de interface Stokes–Darcy

Enunciado do Problema
O artigo aborda a simulação numérica do problema acoplado de interface Stokes–Darcy, um desafio fundamental na mecânica dos fluidos que rege a interação entre regiões de escoamento livre (Stokes) e meios porosos (Darcy). O acoplamento é realizado através da continuidade das velocidades normais (conservação de massa), equilíbrio de tensões normais e a condição de escorregamento tangencial de Beavers–Joseph–Saffman. Embora existam vários métodos de Elementos Finitos (FE), Galerkin Descontínuo (dG) e Híbridos de Alta Ordem (HHO), que frequentemente dependem de pares inf–sup estáveis de velocidade–pressão, este trabalho foca em uma reformulação específica em domínios simplesmente conexos em duas dimensões. Aqui, a restrição de incompressibilidade é satisfeita identicamente através da introdução de uma função de corrente escalar, eliminando assim a variável de pressão das equações de Stokes e reduzindo o sistema a um problema biarmônico de quarta ordem acoplado a uma equação de pressão de Darcy de segunda ordem.

Metodologia
Os autores propõem uma nova discretização por Método de Elementos Virtuais (VEM) para esta formulação de função de corrente–pressão. A metodologia aproveita a capacidade do VEM de lidar com malhas poligonais gerais, o que é particularmente vantajoso para geometrias de interface complexas e malhas não coincidentes (non-matching).

1. Formulação: O problema contínuo é formulado em uma forma primal-primal. O domínio Stokes utiliza uma função de corrente $\chi \in H^2$, enquanto o domínio Darcy utiliza um potencial de pressão $\phi \in H^1$. A formulação fraca resulta em um problema de ponto de sela perturbado.
2. Discretização:
   • Domínio Stokes: É empregado um espaço VEM de ordem mínima $C^1$-conforme para a função de corrente. Isso envolve polinômios quadráticos ($P_2$) e requer conformidade $H^2$, alcançada através de espaços virtuais aumentados com graus de liberdade (DoFs) definidos por valores nodais e gradientes nodais.
   • Domínio Darcy: É utilizado um espaço VEM de ordem mínima $C^0$-conforme para a pressão, baseado em polinômios lineares ($P_1$).
   • Acoplamento: As condições de interface são impostas naturalmente dentro da formulação fraca. O método acomoda malhas não coincidentes na interface sem a necessidade de remalhamento.
3. Implementação: O esquema é implementado utilizando a biblioteca de código aberto orientada a objetos vem++. Uma característica específica, a classe vem::vemMesh2dDofFilter, é utilizada para gerenciar os DoFs através de diferentes subdomínios, permitindo a definição simultânea de diferentes espaços VEM em uma única malha e filtrando-os com base em marcadores de elementos.
4. Estabilização: As formas bilineares discretas incluem termos de estabilização ($S_S^K$ e $S_D^K$) para garantir a coercividade no núcleo dos operadores de projeção, seguindo as práticas padrão de VEM (por exemplo, receita diagonal e DOFI-DOFI).

Principais Contribuições
O artigo delineia cinco contribuições primárias:

• Análise Contínua: Desenvolvimento de uma análise contínua para o problema de interface Stokes–Darcy na formulação de função de corrente–pressão, que evita inerentemente a necessidade de condições de estabilidade inf–sup.
• Esquema Discreto: Design de uma discretização VEM de ordem mínima $C^1$–$C^0$. Os autores destacam a eficiência computacional deste esquema, observando que ele envolve apenas $3N_S + N_D + 4N_\Sigma$ graus de liberdade (onde $N_S$, $N_D$ e $N_\Sigma$ são as contagens de vértices interiores no domínio Stokes, domínio Darcy e na interface, respectivamente).
• Análise de Erro: Derivação de estimativas de erro a priori para a proposta de discretização, provando taxas de convergência ótimas sob suposições de regularidade padrão.
• Implementação de Software: Integração de novas funcionalidades na biblioteca vem++ para facilitar a implementação direta de problemas de interface via filtragem de malha.
• Validação Numérica: Execução de experimentos numéricos para confirmar os achados teóricos.

Resultados
Os experimentos numéricos validam as taxas de convergência teórica e a robustez do método:

• Soluções Fabricadas: Dois experimentos utilizando soluções fabricadas (um com valores de interface nulos, outro com valores não nulos) em vários tipos de malha (quadrilátero, não convexa, Voronoi e Voronoi perturbada) demonstram um decaimento assintótico linear no erro ($O(h)$), consistente com as previsões teóricas para o esquema de ordem mínima.
• Filtro de Fundo Cego (Dead-End Filtration): Uma simulação de um filtro de fundo cego em forma de quarto de anel recupera com sucesso resultados de estudos anteriores, demonstrando a capacidade do método de lidar com interfaces curvas e condições de contorno mistas.
• Órgãos Bioartificiais: Uma simulação do fluxo sanguíneo em um órgão bioartificial (baseado em encapsulamento) com constantes de permeabilidade variáveis ilustra a capacidade do método de modelar fluxos fisiológicos complexos e a influência da permeabilidade no fluxo através da interface.

Significância e Alegações
O artigo posiciona seu trabalho como uma ponte para uma lacuna na literatura referente à extensão das formulações de função de corrente para problemas acoplados Stokes–Darcy utilizando Elementos Virtuais. Os autores afirmam que a estrutura proposta oferece uma "estratégia completamente livre de inf–sup", simplificando o design de espaços discretos em comparação com as abordagens tradicionais de velocidade–pressão. A capacidade de lidar com geometrias de interface irregulares e malhas não coincidentes sem remalhamento é destacada como uma vantagem fundamental do framework VEM.

Os autores mantêm a modéstia quanto a desafios futuros, observando explicitamente que a extensão desta formulação de função de corrente para domínios simplesmente conexos (onde a função de corrente não é unicamente determinada por condições de contorno padrão) e o escalonamento do framework para três dimensões (exigindo potenciais vetoriais e espaços $C^1$ de alta ordem para funções vetoriais) permanecem como desafios abertos. Adicionalmente, a integração de técnicas sem estabilização e a derivação de estimadores de erro a posteriori robustos são identificadas como direções para pesquisas futuras.