High order ghost-FEM for incompressible Navier-Stokes equations on moving domains

1. Problema e Motivação

O artigo aborda a simulação numérica das equações de Navier-Stokes incompressíveis em domínios móveis. O foco principal é modelar o escoamento de um fluido em torno de objetos cujo movimento é pré-definido (acoplamento unidirecional).

• Motivação: O trabalho foi motivado pelo estudo da dinâmica de surfactantes (tensoativos) ao redor de armadilhas em movimento. Surfactantes são cruciais em processos biológicos, industriais e na estabilidade de espumas, e sua interação com interfaces fluidas exige a resolução precisa de equações acopladas de transporte de massa e dinâmica de fluidos.
• Desafio Numérico: A principal dificuldade reside na necessidade de resolver as equações em geometrias complexas e em movimento sem a necessidade de remalhagem (remeshing) constante, o que é computacionalmente custoso e pode introduzir erros de conservação. Métodos tradicionais de malha ajustada (fitted) tornam-se inviáveis para movimentos complexos.

2. Metodologia Proposta

Os autores desenvolveram uma nova técnica numérica que combina três pilares principais para lidar com domínios móveis e não ajustados (unfitted meshes):

A. Discretização Espacial: Ghost-FEM (Método de Elementos Finitos Fantasma)

• Utiliza uma malha cartesiana regular que intersecta arbitrariamente o domínio físico $\Omega$.
• O método define um subconjunto de células ativas ($\Omega_{act}$) que intersectam o domínio.
• Para lidar com as células cortadas (cut cells) e evitar problemas de condicionamento (small cut cells), emprega-se uma estratégia de "snapping-back-to-grid" (retorno à grade). Pontos internos próximos à fronteira são "mapeados" para pontos fantasmas na grade, permitindo a aplicação de condições de contorno sem modificar a topologia da malha.
• Espaços de Aproximação: Utiliza elementos finitos mistos Q2/Q1 (velocidade quadrática bilinear, pressão linear bilinear), que são estáveis na condição inf-sup, garantindo precisão de 3ª ordem para velocidade e 2ª ordem para pressão.
• Condições de Contorno: As condições de Dirichlet são impostas fracamente através do Método de Nitsche.

B. Tratamento de Erros Geométricos: Shifted Boundary Method (SBM)

• Como a fronteira real $\Gamma$ é curva e a malha é poligonal, há um erro geométrico ao aplicar condições de contorno na fronteira aproximada $\Gamma_h$.
• Para manter a alta ordem de precisão, o método utiliza o Shifted Boundary Method. Ele corrige as condições de contorno em $\Gamma_h$ utilizando expansões de Taylor, mapeando pontos de $\Gamma_h$ para a fronteira física real $\Gamma$ através de uma projeção ortogonal. Isso permite que a ordem de convergência seja preservada mesmo em malhas não ajustadas.

C. Discretização Temporal: Esquemas IMEX (Implicit-Explicit)

• Para lidar com a não linearidade do termo convectivo e a rigidez do sistema, o tempo é discretizado usando esquemas IMEX (Implícito-Explícito).
• Foram implementados e testados esquemas de 2 estágios (2ª ordem) e 4 estágios (3ª ordem).
• A pressão é tratada de forma totalmente implícita, resolvendo o sistema de velocidade/pressão de maneira monolítica.
• Para domínios móveis, utiliza-se uma técnica de extrapolação quadrática (baseada no método de Aslam) para estender a solução do fluido para fora do domínio atual, permitindo que a solução seja calculada nos nós ativos do próximo passo de tempo sem perda de precisão.

D. Estabilização

• O parâmetro de penalização $\lambda$ (necessário para o método de Nitsche e para garantir a coercividade) é escolhido resolvendo um problema de autovalor generalizado. Isso garante que o sistema linear permaneça estável e bem condicionado, evitando a necessidade de escolhas heurísticas.

3. Contribuições Principais

1. Integração IMEX com Ghost-FEM: É uma das primeiras aplicações de esquemas IMEX de alta ordem em conjunto com métodos de fronteira imersa (ghost-FEM) para equações de Navier-Stokes em domínios móveis.
2. Alta Ordem de Precisão: Demonstra que é possível atingir convergência de 3ª ordem em tempo e espaço para problemas com domínios móveis e geometrias complexas, superando limitações de métodos de baixa ordem.
3. Robustez Geométrica: O uso combinado de Ghost-FEM, SBM e extrapolação de Aslam permite simulações estáveis em malhas cartesianas fixas, eliminando a necessidade de remalhagem dinâmica.
4. Seleção Automática de Parâmetros: A abordagem baseada em autovalores para o parâmetro de estabilização de Nitsche oferece uma escolha robusta e teoricamente fundamentada, em vez de empírica.

4. Resultados Numéricos

Os autores validaram o método através de uma série de testes:

• Testes de Precisão (Domínio Estático e Móvel): Utilizando soluções analíticas (problema de Kim e Moin) em geometrias como discos, elipses e formas floridas. Os resultados confirmaram as taxas de convergência esperadas:
  • Esquema de 2ª ordem: Erros $L^2$ e $H^1$ convergindo com taxa 2.
  • Esquema de 3ª ordem: Erro $L^2$ convergindo com taxa 3 e $H^1$ com taxa 2.
• Fluxo ao Redor de um Disco (Re=100): Simulação de vórtices de von Kármán. Os coeficientes de arrasto e sustentação, bem como o número de Strouhal, estiveram em excelente acordo com benchmarks da literatura.
• Cavidade Movida (Driven Cavity) com Obstáculo: Análise de vórtices em torno de um obstáculo em forma de flor para diferentes números de Reynolds (1 a 1000), mostrando a capacidade do método em capturar estruturas de fluxo complexas.
• Objeto Rotativo e Bolhas Oscilantes: Simulações de um objeto em forma de flor girando e bolhas oscilando (transversalmente e com deformação elipsoidal). Os resultados mostraram a formação periódica de vórtices e a dinâmica correta do escoamento em baixos números de Reynolds.

5. Significado e Conclusões

O trabalho representa um avanço significativo na simulação de fluidos em domínios complexos e em movimento. A metodologia proposta oferece uma alternativa eficiente e precisa aos métodos tradicionais de malha ajustada, especialmente para aplicações onde a geometria muda dinamicamente.

• Estado Atual: O método atualmente lida com acoplamento unidirecional (o movimento do objeto é pré-definido e afeta o fluido, mas o fluido não altera o movimento do objeto).
• Futuro: O objetivo final é aplicar esta técnica ao estudo da cinética de sorção de surfactantes em bolhas oscilantes, incorporando o acoplamento bidirecional completo. Isso permitirá simular cenários onde a dinâmica do fluido influencia o movimento de múltiplas bolhas interagentes, um problema complexo e relevante para a indústria e biologia.

Em resumo, o artigo apresenta uma ferramenta numérica robusta, de alta ordem e geometricamente flexível, pronta para expandir as fronteiras da simulação de escoamentos multifásicos e interações fluido-estrutura.