Ensemble algorithm for parametrized flow problems with energy stable open boundary conditions

Resumo Técnico: Algoritmo de Ensemble para Problemas de Fluxo Parametrizados com Condições de Contorno Abertas Estáveis em Energia

Enunciado do Problema
O artigo aborda o desafio computacional de resolver ensembles de equações de Navier-Stokes incomprimíveis sujeitas a variações nas condições iniciais, forças de corpo e coeficientes de viscosidade ($\nu_j$). Em aplicações como análise de sensibilidade, modelagem de ordem reduzida e previsão de ensemble, o número de membros do ensemble ($J$) pode ser grande. Métodos semi-implícitos padrão exigem a montagem e a resolução de um sistema linear distinto para cada membro $j$ em cada passo de tempo, devido à dependência da matriz do sistema em relação a $\nu_j$ e aos termos não lineares. Isso resulta em custos computacionais proibitivos. Abordagens totalmente explícitas permitem uma única montagem de sistema, mas impõem condições de passo de tempo restritivas, particularmente para fluxos de alto número de Reynolds. Métodos de ensemble anteriores para viscosidade constante ou variações específicas de viscosidade frequentemente exigiam suposições restritivas sobre os coeficientes de viscosidade (por exemplo, diferenças relativas limitadas) para garantir a estabilidade.

Metodologia
Os autores propõem um novo método de cálculo de ensemble que desacopla a matriz do sistema linear do índice do membro do ensemble $j$, permitindo a solução de um único sistema linear com múltiplos lados direitos em cada passo de tempo. O núcleo da metodologia envolve:

1. Estratégia de Decomposição (Splitting): O termo convectivo não linear é decomposto em um componente de média do ensemble e um componente de flutuação. O componente de média é tratado semi-implicitamente, enquanto o componente de flutuação é tratado explicitamente.
2. Tratamento da Viscosidade: Diferente de trabalhos anteriores (ex: [8]) que exigiam limites específicos na variação dos coeficientes de viscosidade para garantir estabilidade, este artigo introduz um tratamento diferente para o termo difusivo. Ao definir uma viscosidade máxima $\nu_\infty = \max_j \nu_j$ e decompor o termo de difusão, os autores derivam um esquema onde a matriz do sistema linear é independente de $j$ sem impor restrições sobre a magnitude das variações de viscosidade.
3. Condições de Contorno Abertas (OBC): O método é estendido para domínios com contornos abertos ($\Gamma_N$). O contorno é particionado em regiões de escoamento (outflow) e refluxo (backflow). Os autores utilizam uma condição de contorno estável em energia envolvendo uma escala de comprimento característica $L$ e uma função do tipo Heaviside para lidar com a natureza convectiva do escoamento de saída.
4. Discretização: O artigo apresenta esquemas de integração temporal de primeira ordem (Algoritmos 1–3) e de segunda ordem (Algoritmos 4–6). A discretização espacial utiliza espaços de elementos finitos conformes (especificamente pares Taylor-Hood $P_2/P_1$) com formas trilineares assimetrizadas para garantir a estabilidade de energia.
5. Análise de Estabilidade: Provas teóricas estabelecem a estabilidade de energia para os esquemas de primeira ordem sob condições CFL específicas envolvendo flutuações de velocidade. Para os esquemas de segunda ordem, embora uma prova de estabilidade rigorosa sem restrições de viscosidade não seja fornecida (devido à complexidade do termo $\nabla \partial_{tt} u$), experimentos numéricos sugerem estabilidade sob as mesmas condições dos esquemas de primeira ordem.

Principais Contribuições

• Estabilidade de Viscosidade Irrestrita: A principal contribuição teórica é um esquema de ensemble para coeficientes de viscosidade variáveis que garante estabilidade sem a suposição restritiva $|\nu_j - \nu|/\nu \leq \sqrt{\mu}$ encontrada na literatura anterior.
• Extensão para Condições de Contorno Abertas: Os autores estendem os métodos de ensemble para problemas com contornos abertos, derivando formulações fracas e provando a estabilidade de energia para esquemas de primeira ordem sob restrições específicas de passo de tempo.
• Eficiência Algorítmica: Os algoritmos propostos resultam em uma única matriz de coeficientes por passo de tempo, solucionável eficientemente usando solvers de múltiplos lados direitos, reduzindo significativamente o custo computacional em comparação com a resolução de $J$ sistemas independentes.

Resultos
Experimentos numéricos foram conduzidos utilizando o pacote FreeFem++:

1. Estudo de Convergência: Um teste de convergência em um domínio com soluções exatas confirmou as taxas de convergência de segunda ordem esperadas tanto para velocidade quanto para pressão para os esquemas de segunda ordem propostos. A precisão foi comparável à de simulações sequenciais independentes.
2. Escoamento ao Redor de um Cilindro: Os algoritmos foram testados em um problema de referência envolvendo o escoamento ao redor de um cilindro com três valores de viscosidade diferentes. Os resultados (arrasto, sustentação e queda de pressão) coincidiram com os valores de referência da literatura usando condições de Dirichlet e de escoamento de saída do tipo "do-nothing". A condição de contorno aberta permitiu que os vórtices saíssem do domínio naturalmente, evitando os problemas de reflexão artificial associados a perfis parabólicos de saída prescritos.
3. Escoamento em Canal com Contração: Um cenário de escoamento complexo com duas saídas e viscosidades variadas foi simulado. Os resultados do ensemble coincidiram qualitativamente com execuções independentes de Simulação Numérica Direta (DNS), apesar de as simulações de ensemble terem sido realizadas em uma malha mais grossa estabilizada por um filtro não linear adaptativo.

Significância e Alegações
O artigo afirma oferecer um arcabouço robusto e eficiente para simulações de fluxo parametrizadas. Ao remover as suposições restritivas sobre as variações de viscosidade exigidas por métodos de ensemble anteriores, o algoritmo proposto amplia a aplicabilidade das simulações de ensemble para problemas com parâmetros físicos altamente variáveis. A extensão para condições de contorno abertas com estabilidade de energia provável é um passo significativo em direção a simulações realistas de fluxos ilimitados. Os autores concluem que, embora os esquemas de ensemble apresentem bom desempenho qualitativo e quantitativo em relação às simulações independentes, pesquisas adicionais são necessárias para quantificar totalmente as vantagens computacionais, particularmente no que diz respeito à comparação de tempo de CPU em malhas adaptativamente refinadas para números de Reynolds elevados.