Evolving finite elements for advection diffusion with an evolving interface

Resumo Técnico: Elementos Finitos Evolutivos para Advecção-Difusão com uma Interface Evolutiva

Enunciação do Problema
Este artigo aborda a aproximação numérica de equações parabólicas definidas em um domínio estacionário $\Omega$ contendo uma interface móvel $\Gamma(t)$. A interface divide o domínio em dois subdomínios, $\Omega_1(t)$ e $\Omega_2(t)$. As equações governantes modelam o transporte e a difusão de campos escalares $u_1$ e $u_2$ dentro desses subdomínios, sujeitos a:

• Dinâmicas de advecção-difusão-reação dentro de cada subdomínio com coeficientes dependentes do tempo ($A_i, B_i, C_i$).
• Uma condição de contorno de Dirichlet na fronteira externa $\partial\Omega$.
• Continuidade da solução através da interface ($u_1 = u_2$).
• Um salto de fluxo prescrito através da interface ($A_1 \partial_{\nu} u_1 - A_2 \partial_{\nu} u_2 = g$).

Os principais desafios identificados são a evolução dos subdomínios ao longo do tempo e a presença do salto de fluxo através da interface móvel. Embora existam aplicações em dinâmica de fluidos, ciência dos materiais e biologia celular, o artigo foca na formulação matemática e numérica do problema.

Metodologia
Os autores propõem uma abordagem Lagrangiana-Euleriana Arbitrária (ALE) utilizando Espaços de Elementos Finitos Evolutivos (ESFEM) baseados na estrutura abstrata estabelecida em [21].

1. Contexto Analítico Funcional: O problema contínuo é formulado no contexto de espaços de Hilbert evolutivos. Os autores definem pares compatíveis de espaços $(X(t), \phi_t)$, onde $\phi_t$ é um mapa de fluxo transportando funções de um domínio de referência para o domínio evolutivo. Isso permite a definição de derivadas materiais fracas e a derivação de uma formulação fraca bem-posta (Teorema 2.16).
2. Construção da Malha:
   • Malha Inicial: O domínio é inicialmente triangulado usando símplices lineares por partes. Crucialmente, os vértices das facetas da interface são restritos a cair exatamente sobre a interface inicial $\Gamma(0)$ (Hipótese M5).
   • Mapeamento Isoparamétrico: Para alcançar maior precisão de ordem, a malha linear é mapeada para uma malha isoparamétrica usando um difeomorfismo $\Psi_h$. Este mapa empurra as facetas da interface para a interface curva exata $\Gamma(0)$, preservando a continuidade da malha.
   • Evolução Temporal: A malha evolui movendo os nós Lagrangianos de acordo com um campo de velocidade global conhecido e suave $w$. O fluxo discreto $\Phi^h$ é construído como uma interpolação polinomial do fluxo contínuo, garantindo que a malha permaneça uma subdivisão admissível conformante em todos os tempos.
3. Mapeamento de Elevação (Lift): Um operador de elevação $\Lambda_h$ é introduzido para mapear funções do domínio computacional discreto (que aproxima a interface curva) de volta para o domínio contínuo exato. Isso é essencial para comparar a solução discreta com a solução exata na análise de erro.
4. Esquema Discreto: O método emprega elementos finitos isoparamétricos de ordem $k$. O problema variacional discreto é formulado usando a derivada material discreta e os coeficientes elevados. O esquema é resolvido usando um método de passo temporal (especificamente Fórmulas de Diferenças Regressivas, BDF, nos experimentos numéricos).

Principais Contribuições

• Bem-postura: O artigo fornece uma estrutura analítica funcional rigorosa provando a existência e unicidade da solução para o problema contínuo (Teorema 2.16), incluindo resultados adicionais de regularidade sob hipóteses específicas sobre os dados.
• Estrutura ESFEM Isoparamétrica: Os autores estendem a estrutura abstrata de elementos finitos evolutivos para lidar com problemas de interface de dois lados com saltos de fluxo. Eles detalham a construção de elementos isoparamétricos que aproximam com precisão a fronteira curva, permitindo convergência de ordem superior.
• Limites de Erro Ótimos: O principal resultado teórico (Teorema 5.1) estabelece limites de erro de ordem ótima. Especificamente, o erro na norma $L^2$ é limitado por $O(h^{k+1})$, e o erro na semi-norma $H^1$ é limitado por $O(h^k)$, onde $h$ é o tamanho da malha e $k$ é o grau polinomial. Esses limites valem até uma constante e correspondem ao erro de interpolação de uma solução suave conhecida.
• Implementação e Verificação: O artigo inclui uma implementação robusta usando o pacote Firedrake e o GMSH. Experimentos numéricos em 2D e 3D verificam as taxas de convergência teóricas para várias ordens polinomiais.

Resultados
Os experimentos numéricos demonstram a eficácia do método em um problema modelo com uma solução exata conhecida que é contínua, mas não diferenciável através da interface.

• Convergência: Para discretização espacial de ordem $k$, os resultados numéricos confirmam convergência de ordem $k+1$ na norma $L^2$ e ordem $k$ na semi-norma $H^1$.
• Robustez: Os resultados valem para domínios tanto 2D quanto 3D, mesmo quando as malhas são geradas independentemente (em vez de via refinamento uniforme), indicando a estabilidade do método.
• Disponibilidade de Código: Os autores fornecem o código de simulação para permitir a reprodutibilidade e o uso adicional pela comunidade.

Significado e Escopo
O artigo afirma fornecer uma estrutura teórica e prática completa para resolver problemas de advecção-difusão com interfaces móveis usando elementos finitos isoparamétricos de alta ordem. O significado reside em preencher a lacuna entre a teoria abstrata de espaços evolutivos e a implementação prática para problemas de interface com saltos de fluxo.

Os autores mantêm um escopo modesto em relação ao campo de velocidade: eles assumem que um campo de velocidade global e suave é dado a priori. Eles afirmam explicitamente que o trabalho não aborda como computar tal campo de velocidade (por exemplo, via re-parametrização para evitar emaranhamento da malha), observando que tais métodos existem na literatura [8, 18, 19]. Além disso, embora o método pudesse teoricamente aplicar-se a sistemas de EDPs acoplados onde a evolução do domínio faz parte da solução, a análise atual não cobre tais problemas acoplados. O foco permanece estritamente na aproximação do problema de transmissão parabólico dada uma evolução de domínio conhecida.