Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems

Resumo Técnico: Estimativas de Erro de Quadratura em Malhas Não Correspondentes em um Framework de Domínio Fictício para Problemas de Interação Fluido-Estrutura

Enunciado do Problema
O artigo aborda a discretização numérica de problemas de interação fluido-estrutura (FSI) utilizando uma formulação de domínio fictício baseada em multiplicadores de Lagrange distribuídos. Neste framework, o fluido e o sólido são discretizados usando duas malhas completamente independentes e não correspondentes. O sólido é descrito em um framework lagrangiano em um domínio de referência $B$, mapeado para sua posição atual dentro de um domínio euleriano fixo de fluido $\Omega$. O acoplamento cinemático entre a velocidade do fluido e a deformação do sólido é imposto via um multiplicador de Lagrange $\lambda$ definido no domínio de referência do sólido.

Um desafio computacional crítico surge na montagem da matriz de acoplamento (interface). Esta matriz requer a integração do produto das funções de base do fluido (compostas com o mapeamento do sólido) e das funções de base do sólido sobre os elementos do sólido. Como essas funções são definidas em malhas não correspondentes, a integração exata exige o cálculo da interseção geométrica entre as malhas do fluido e do sólido mapeado para construir uma regra de quadratura composta. Este processo é computacionalmente caro, particularmente porque o termo de acoplamento deve ser montado em cada passo de tempo. Alternativamente, uma abordagem aproximada integra diretamente sobre os elementos do sólido usando regras de quadratura padrão sem a interseção de malhas, introduzindo um erro de quadratura. O artigo investiga a validade teórica e as estimativas de erro associadas a esta estratégia de integração aproximada.

Metodologia
Os autores analisam um modelo de FSI baseado em Stokes em estado estacionário, embora observem que seus resultados são representativos de casos não lineares mais gerais. A formulação variacional envolve encontrar a velocidade do fluido $u$, a pressão $p$, o mapa de deformação do sólido $X$ e o multiplicador de Lagrange $\lambda$ satisfazendo um sistema de ponto de sela.

O estudo foca em duas escolhas específicas para a forma bilinear $c(\cdot, \cdot)$ usada para definir o acoplamento:

1. $c_0$ (Pareamento dual): Baseado no produto interno $L^2(B)$, onde o espaço do multiplicador é o dual de $H^1(B)$.
2. $c_1$ (Produto escalar): Baseado no produto interno $H^1(B)$, onde o espaço do multiplicador é o próprio $H^1(B)$.

A metodologia procede através dos seguintes passos:

• Discretização: Espaços de elementos finitos são introduzidos para todas as variáveis. O fluido utiliza o elemento de Bercovier–Pironneau ($P1$-iso-$P2/P1$), enquanto o mapeamento do sólido e o multiplicador utilizam elementos contínuos de polinômio linear por partes.
• Integração Aproximada: Os autores definem formas bilineares discretas onde os termos de acoplamento são computados usando regras de quadratura padrão na malha do sólido, ignorando o fato de que as funções de base do fluido compostas são polinômios por partes com descontinuidades em relação à malha do sólido.
• Análise de Erro Abstrata: Seguindo o framework das condições de inf-sup, os autores derivam uma estimativa de erro geral (Teorema 2) que limita a diferença entre a solução discreta exata e a solução obtida com integração inexata. Este limite depende dos funcionais de erro de quadratura.
• Estimativa de Erro de Quadratura: A principal contribuição técnica envolve provar estimativas específicas para o erro de quadratura dos termos de acoplamento. Os autores utilizam o fato de que a função composta $v_h(X)$ é contínua, mas linear por partes sobre os elementos do sólido, levando a propriedades de regularidade específicas em espaços de Sobolev fracionários ($H^{1+s}$). Eles empregam teoria de interpolação e desigualdades inversas para limitar o erro.
• Prova de Estabilidade: Os autores provam que o problema discreto permanece bem posto (satisfaz as condições de inf-sup discretas) mesmo quando o termo de acoplamento é construído aproximadamente, desde que condições específicas sobre os tamanhos das malhas sejam atendidas.

Principais Contribuições e Resultados

1. Bem-posto com Acoplamento Aproximado: O artigo prova que o problema discreto é estável e bem posto mesmo quando a matriz de acoplamento é montada usando quadratura aproximada, sem exigir a interseção de malhas.
2. Estimativas de Erro de Quadratura:
   • Para o acoplamento baseado em $L^2$ ($c_0$), o erro de quadratura é limitado por $O(h_B^{3/2} |\log h_B|)$. Os autores mostram que o método aproximado produz taxas de convergência ótimas, desde que o tamanho da malha do sólido $h_B \to 0$.
   • Para o acoplamento baseado em $H^1$ ($c_1$), a análise de erro é mais complexa devido ao termo de gradiente. O limite de erro inclui termos dependentes da razão entre os tamanhos das malhas. Especificamente, a convergência ótima requer não apenas $h_B \to 0$, mas também que a razão $h_B/h_\Omega \to 0$ (onde $h_\Omega$ é o tamanho da malha do fluido). Se $h_B/h_\Omega$ for mantida constante, o método pode perder a otimalidade.
3. Desigualdades Inversas: Os autores derivam e utilizam desigualdades inversas específicas para lidar com as normas do multiplicador de Lagrange nas estimativas de erro, particularmente para o caso da norma dual.
4. Validação Numérica: Quatro testes numéricos são apresentados para validar as estimativas teóricas:
   • Teste 1 & 2: Comparam taxas de convergência para $c_0$ e $c_1$ sob diferentes estratégias de refinamento de malha (razão constante vs. razão decrescente). Os resultados confirmam que $c_0$ é robusto com integração aproximada, enquanto $c_1$ requer $h_B/h_\Omega \to 0$ para manter a precisão.
   • Teste 3: Utiliza uma solução com um salto no gradiente de velocidade através da interface (típico de FSI), confirmando a taxa de convergência subótima esperada de $O(h^{1/2})$ devido a singularidades na solução.
   • Teste 4: Investiga elementos de ordem superior ($P2/P1$), mostrando que, embora $c_0$ permaneça robusto, $c_1$ ainda requer integração exata ou razões de malha específicas para atingir taxas ótimas.

Significância e Alegações
O artigo alega fornecer uma justificativa teórica rigorosa para o uso de integração aproximada em métodos de FSI de domínio fictício, uma prática frequentemente adotada para eficiência computacional. Os autores demonstram que:

• A abordagem "inexata" não compromete a estabilidade do esquema numérico.
• O acoplamento baseado em $L^2$ ($c_0$) é particularmente adequado para integração aproximada, pois não exige que as malhas do fluido e do sólido sejam refinadas em uma razão específica para manter a precisão.
• O acoplamento baseado em $H^1$ ($c_1$), embora possa oferecer diferentes propriedades de estabilidade, impõe requisitos de malha mais estritos ($h_B/h_\Omega \to 0$) para evitar que os erros de quadratura dominem o erro da solução.

Os autores concluem que seus resultados teóricos confirmam as observações experimentais de seu trabalho anterior [12, 13], fornecendo uma base matemática sólida para a escolha entre técnicas de montagem exata (com interseção de malhas) e aproximada (sem interseção) baseada no tipo de forma de acoplamento e configuração de malha.