Free-surface Stokes problem: stability estimates and time-step improvements

Resumo Técnico: Problema de Stokes de superfície livre: estimativas de estabilidade e melhorias no passo de tempo

Enunciado do Problema
A simulação precisa de fluxos de superfície livre altamente viscosos, como a dinâmica de camadas de gelo, convecção de manto e fluxo de lava, requer a resolução das equações acopladas de Stokes/superfície livre. Estes problemas são caracterizados por um domínio móvel onde a velocidade do fluido e a pressão estão acopladas à evolução da altura da superfície livre. Embora as discretizações padrão usando elementos finitos no espaço e métodos de Euler no tempo sejam comuns, elas frequentemente sofrem de restrições de estabilidade severas. Especificamente, esquemas de Euler explícitos e semi-implícitos frequentemente exigem passos de tempo proibitivamente pequenos para permanecerem estáveis, levando a altos custos computacionais. As técnicas de estabilização existentes, como o Algoritmo de Estabilização de Superfície Livre (FSSA), melhoram os tamanhos dos passos de tempo, mas carecem de limites teóricos a priori sobre estabilidade e não garantem estabilidade incondicional ou conservação de volume.

Metodologia
Os autores apresentam uma análise teórica rigorosa das propriedades de estabilidade e conservação do sistema acoplado Stokes/superfície livre para fluidos Newtonianos e não Newtonianos (viscosidade de lei de potência). A análise é conduzida em dois níveis:

1. Nível Contínuo: Os autores derivam um balanço de energia total que acopla a energia das equações de momento de Stokes com a energia da equação da altura da superfície livre. Esta derivação baseia-se no estabelecimento de relações específicas entre o termo de advecção da superfície livre e integrais de superfície, vinculando-os posteriormente ao termo de gravidade nas equações de Stokes.
2. Nível Discreto: O estudo examina esquemas totalmente discretos usando elementos finitos no espaço combinados com integração temporal de Euler explícita, semi-implícita e implícita.
   • Regularização: Para abordar a falta de regularidade inerente à função de altura da superfície livre (necessária para a correção do problema de Stokes), os autores introduzem um termo de estabilização de aresta de alta ordem, $J(s^{n+1}, w)$, tratado implicitamente no tempo.
   • Proposta de Estabilização: Motivados pela identificação de termos "espúrios" nas estimativas de energia do método de Euler explícito, os autores propõem um novo termo de estabilização adicionado à forma fraca de Stokes. Este termo penaliza a velocidade normal na superfície livre e é projetado para cancelar precisamente os termos desestabilizadores decorrentes do tratamento explícito da equação da superfície livre.

Principais Contribuições

• Balanço de Energia Total: O artigo deriva um novo balanço de energia total para o sistema acoplado Stokes/superfície livre, acoplando a energia do momento com a energia da altura da superfície. Isto foi identificado como um bloco fundamental para a compreensão da estabilidade, que estava ausente na literatura.
• Provas de Estabilidade Teórica:
  • Euler Implícito: Provado ser incondicionalmente estável em relação ao tamanho do passo de tempo (assumindo continuidade de Lipschitz da altura da superfície) e conserva o volume do domínio.
  • Euler Explícito (Não Estabilizado): Provado ser instável na norma de energia total para qualquer passo de tempo devido a termos espúrios relacionados à velocidade normal da superfície.
  • Euler Semi-implícito (FSSA): Mostrado como sendo apenas condicionalmente estável e, criticamente, provado não conservar o volume do domínio.
  • Euler Explícito Estabilizado Proposto: Os autores provam que o seu novo termo de estabilização torna o método de Euler explícito incondicionalmente estável em relação ao tamanho do passo de tempo (sob a hipótese de continuidade de Lipschitz), preservando simultaneamente o volume do domínio.
• Conservação de Volume: O artigo fornece provas teóricas e verificação numérica de que o esquema explícito estabilizado proposto conserva o volume do domínio, enquanto o esquema semi-implícito estabilizado pelo FSSA não o faz.

Resultos
Experimentos numéricos foram conduzidos utilizando a biblioteca de código aberto FEniCS para casos de fluidos Newtonianos e não Newtonianos, incluindo uma geometria complexa baseada em dados da camada de gelo da Gronelândia.

• Estabilidade: O esquema de Euler explícito não estabilizado foi instável para todos os passos de tempo testados. Os esquemas estabilizados pelo FSSA (tanto semi-implícitos quanto explícitos) foram condicionalmente estáveis, falhando em passos de tempo maiores. Em contraste, o esquema EE-Stabilized permaneceu estável para passos de tempo arbitrariamente grandes, tanto em casos Newtonianos quanto não Newtonianos.
• Conservação de Volume: O esquema EE-Stabilized conservou o volume do domínio até à precisão de máquina. O esquema SIE-FSSA falhou em conservar o volume, confirmando os achados teóricos.
• Precisão: Testes de convergência mostraram que os termos de estabilização propostos são consistentes com o problema discretizado. O esquema EE-Stabilized manteve a precisão de primeira ordem na altura da superfície e na velocidade sob refinamento simultâneo de malha e passo de tempo, com erros comparáveis ou menores do que outros esquemas estabilizados.
• Oscilações: O termo de estabilização de aresta implícito suprimiu eficazmente as oscilações de alta frequência na altura da superfície livre sem introduzir restrições adicionais ao passo de tempo.

Significância e Alegações
O artigo afirma que o esquema de Euler explícito estabilizado proposto oferece uma alternativa robusta e eficiente aos solvers implícitos computacionalmente caros. Ao garantir teoricamente a estabilidade incondicional em relação ao tamanho do passo de tempo e assegurar a conservação do volume, o método permite simulações estáveis de fluxos de superfície livre mesmo com geometrias complexas e termos de fonte. Os autores enfatizam que a sua estabilização é a primeira a garantir teoricamente passos de tempo arbitrários para a discretização de Euler explícito deste problema acoplado. A abordagem é notada por ser simétrica (benéfica para solvers iterativos), simples de implementar e livre de parâmetros ajustáveis. Os autores concluem que, embora a sua análise assuma a continuidade de Lipschitz, as estimativas de estabilidade de norma $L^2$ derivadas são uma base necessária para a correção uma vez aplicada a regularização apropriada. Sugere-se como trabalho futuro explorar termos de regularização para a correção e estender esta análise de estabilidade para problemas de Navier-Stokes/superfície livre e integradores de tempo de ordem superior.