A Nitsche method for Navier--Stokes/generalized poroelasticity interface problems

Resumo Técnico: Um Método de Nitsche para Problemas de Interface Navier–Stokes/Poroelasticidade Generalizada

Enunciação do Problema
Este trabalho aborda a modelagem matemática e numérica da interação entre um fluido livre incompressível e um meio poroso deformável, um fenômeno conhecido como Interação Fluido–Estrutura Poroelástica (FPSI). A região do fluido livre é governada pelas equações de Navier–Stokes dependentes do tempo, enquanto o meio poroso é modelado usando um sistema poro-hiperelástico generalizado e linearizado (incorporando o modelo de Brinkman para efeitos viscosos dentro da matriz porosa). Os dois domínios são acoplados através de um conjunto de condições físicas de interface: continuidade da velocidade normal, equilíbrio das tensões normais, equilíbrio das forças de contato e a condição de deslizamento com atrito de Beavers–Joseph–Saffman (BJS) para o fluxo tangencial. Os autores concentram-se no desafio de impor essas restrições de interface dentro de uma estrutura de elementos finitos (EF), particularmente ao lidar com malhas não coincidentes na interface.

Metodologia
O artigo propõe uma discretização de elementos finitos unificada e monolítica baseada em integração temporal implícita (especificamente o método de Euler implícito). Para lidar com a interface fluido-estrutura sem introduzir incógnitas adicionais (como multiplicadores de Lagrange), os autores empregam uma formulação do tipo Nitsche. Esta abordagem impõe fracamente a continuidade da velocidade normal através da interface usando um parâmetro de penalidade $\gamma$.

Componentes metodológicos-chave incluem:

• Formulação Fraca: O problema contínuo é convertido em uma forma variacional envolvendo formas bilineares para os operadores de Navier–Stokes, Brinkman e elasticidade, juntamente com termos de interface derivados das condições BJS e da penalidade de Nitsche.
• Discretização: A discretização espacial utiliza espaços padrão de elementos finitos (polinômios de grau $k$ e $k+1$) para velocidades, pressões e deslocamentos. A discretização temporal utiliza o esquema de Euler implícito.
• Análise Teórica:
  • Bem-postura: A existência e unicidade da solução semidiscreta são estabelecidas usando a teoria de Equações Diferenciais-Algebraicas (EDAs) e o teorema do ponto fixo de Banach. Esta análise depende de condições específicas relativas ao parâmetro de penalidade $\gamma$ e à magnitude das forças externas.
  • Estabilidade: Uma análise de estabilidade abrangente é realizada para o esquema totalmente discreto, derivando estimativas de energia que garantem que a solução permaneça limitada.
  • Estimativas de Erro: Estimativas de erro a priori são derivadas para o esquema totalmente discreto, demonstrando taxas de convergência dependentes do tamanho da malha $h$ e do passo de tempo $\tau$. A análise detalha explicitamente a dependência dessas estimativas em relação ao parâmetro de estabilização $\gamma$.

Principais Contribuições
Os autores identificam seu trabalho como o primeiro a abordar a discretização monolítica de problemas de interface Navier–Stokes/poroelasticidade generalizada especificamente usando o método de Nitsche. Suas principais contribuições são:

1. Extensão da Análise: O trabalho estende estudos anteriores sobre Stokes–poroelasticidade generalizada, incorporando o termo convectivo não linear das equações de Navier–Stokes e aplicando o método de Nitsche para impor fracamente a continuidade da velocidade normal.
2. Fundamentação Teórica: Eles estabelecem a bem-postura do problema semidiscreto e fornecem uma análise de convergência completa para o esquema totalmente discreto, incluindo estimativas de estabilidade e erro.
3. Validação Numérica: O artigo apresenta experimentos numéricos para validar as taxas de convergência teóricas e demonstrar a capacidade do método de capturar dinâmicas acopladas.

Resultados Numéricos
Os autores apresentam três conjuntos de testes numéricos:

1. Testes de Convergência: Usando soluções fabricadas, o método é testado quanto à convergência espacial e temporal. Os resultados confirmam as taxas de convergência teóricas, mostrando convergência ótima para a maioria das variáveis e taxas subótimas (mas aceitáveis) para a velocidade relativa em configurações específicas.
2. Escoamento em Canal 2D: Uma simulação de escoamento em um canal com obstáculos rígidos situados entre camadas porosas. Os resultados ilustram a formação de zonas de recirculação, quedas de pressão através da interface e deformação sólida, destacando o forte acoplamento entre o fluido e o meio poroelástico.
3. Simulação Microfluídica 3D: Uma simulação tridimensional do fluxo sanguíneo através de um chip microfluídico contendo obstáculos poroelásticos cilíndricos. Este teste demonstra a estabilidade e a precisão do método em geometrias 3D complexas, capturando a deformação dos obstáculos e os campos de pressão e velocidade associados.

Significado e Afirmações
O artigo afirma que a abordagem baseada em Nitsche proposta oferece uma alternativa robusta e computacionalmente eficiente aos métodos de multiplicadores de Lagrange. Ao evitar a introdução de incógnitas extras, o método mantém a simetria (ou antissimetria) do sistema algébrico e simplifica a implementação, particularmente para malhas não coincidentes. Os autores afirmam que a estabilidade do esquema é garantida por um parâmetro de penalidade que é independente do tamanho da malha, evitando assim a necessidade de condições adicionais de inf-sup de estabilidade frequentemente requeridas em abordagens baseadas em multiplicadores.

O trabalho conclui que o método captura com sucesso a interação complexa entre a dinâmica do fluido livre e as respostas poroelásticas. Os autores observam que, embora o estudo atual se concentre no modelo poro-hiperelástico linearizado, trabalhos futuros poderiam estender essa estrutura para regimes totalmente não lineares e outras condições de transmissão para abordar uma classe mais ampla de problemas de poromecânica.