Accelerating droplet-laden Stokes flow simulations with hierarchical surrogate modeling

Resumo Técnico: Aceleração de simulações de fluxo de Stokes com gotas utilizando modelagem substituta hierárquica

Definição do Problema
O artigo aborda os desafios computacionais inerentes à simulação de fluxos de Stokes multifásicos incompressíveis envolvendo um grande número ($M$) de gotas de líquido imiscíveis suspensas em um fluido portador. Simulações de alta fidelidade de tais sistemas são computacionalmente exigentes, particularmente quando as gotas estão altamente concentradas, aproximam-se de fronteiras ou interagem de forma próxima. Técnicas de discretização padrão, como métodos de elementos finitos (que exigem remalhamento repetitivo) ou métodos de integral de contorno (que exigem quadratura especializada para integrais quase singulares), têm dificuldade com as lacunas interfaciais estreitas e interações próximas às fronteiras. Embora métodos baseados em partículas evitem o remalhamento, eles frequentemente exigem resoluções finas para precisão e carecem da precisão de métodos de ordem elevada. Os autores visam desenvolver uma estratégia de modelagem substituta que permita a previsão eficiente da evolução das gotas em configurações densamente compactadas e variáveis no tempo, sem o custo proibitivo de simulações de múltiplos fluidos totalmente resolvidas.

Metodologia
A abordagem proposta é uma estrutura de modelagem substituta hierárquica, de múltiplas fidelidades e informada pela física. Ela se baseia em uma estratégia iterativa de correção de erro em vez de tratar o substituto (surrogate) como uma caixa preta. A metodologia procede da seguinte forma:

1. Aproximação de Baixa Fidelidade: O processo começa com uma solução de Stokes de linha de base ($u_0, p_0$), onde as gotas são tratadas como traçadores passivos. Isso assume uma viscosidade "nominal" homogênea em todo o domínio, negligenciando contrastes de viscosidade e tensão superficial.
2. Derivação da Equação de Erro: Os autores derivam uma equação diferencial parcial (PDE) que governa o campo de erro ($v, q$) entre a solução de alta fidelidade e a aproximação de baixa fidelidade. Esta equação de erro contabiliza o salto na viscosidade e nas forças de tensão superficial através das interfaces das gotas.
3. Correção Hierárquica via Iteração: Em vez de resolver a equação de erro global diretamente (o que manteria a complexidade computacional original), o método emprega um esquema iterativo alternado:
   • Correção de Gotas: A equação de erro global é decomposta em $M$ problemas locais independentes, cada um representando a perturbação do fluxo causada por uma única gota em espaço livre. Esses problemas locais são resolvidos de forma aproximada.
   • Correção de Fronteira: A superposição das correções locais viola as condições de contorno globais. Uma etapa subsequente de "correção de fronteira" resolve um problema de Stokes global (sem gotas) para restaurar a consistência com as fronteiras do domínio.
   • Iteração: Estas duas etapas (correções de gota e de fronteira) são repetidas até que um critério de convergência baseado em resíduos de fronteira seja atendido.
4. Eficiência Offline-Online: Para gotas geometricamente semelhantes (especificamente circulares nos testes apresentados), o método explora a autossimilaridade. Os problemas de gota única são mapeados para um disco unitário de referência. Uma abordagem de subespaço linear é usada para aproximar o operador de solução para esses problemas de gota única usando um pequeno conjunto de "soluções fundamentais" pré-computadas correspondentes a campos de fluxo polinomiais. Isso permite que a avaliação online das correções de gota seja realizada via combinações lineares simples, reduzindo significativamente o custo computacional.

Principais Contribuições

• Correção de Erro Baseada na Física: O artigo introduz uma estratégia iterativa de múltiplas fidelidades que recupera sistematicamente soluções de alta fidelidade ao resolver equações de erro baseadas em PDE, evitando a natureza de "caixa preta" de muitos substitutos de aprendizado de máquina.
• Substituição de Refinamento por Iteração: O método evita a alta complexidade de resolver equações de erro globais em domínios densamente compactados ao construir correções localmente através da superposição de campos de gota única. Isso reduz os requisitos de memória e o custo computacional.
• Escalabilidade para Muitas Gotas: A estrutura é projetada para lidar com sistemas de até $10^4$ gotas. Ao desacoplar o problema de fronteira global das interações locais de gotas e utilizar soluções fundamentais pré-computadas, o método alcança ganhos de velocidade substanciais em comparação com simulações totalmente resolvidas.
• Tratamento de Tensão Superficial: Para gotas circulares, os efeitos de tensão superficial são contabilizados exatamente via saltos de pressão, permitando que a correção iterativa foque nos contrastes de viscosidade e interações de fluxo.

Resultados
Experimentos numéricos em duas dimensões validam a precisão e a eficiência do substituto proposto:

• Precisão: O substituto atinge erros relativos abaixo de $\sim 0,1\%$ em comparação com simulações CutFEM de alta fidelidade totalmente resolvidas. Testes de convergência mostram convergência exponencial em relação ao grau polinomial das soluções fundamentais e convergência robusta das iterações de ponto fixo, mesmo para gotas em proximidade estreita com fronteiras ou entre si.
• Eficiência: Em uma simulação no domínio do tempo com 50 gotas, o substituto alcançou uma aceleração de 80 vezes sobre o código CutFEM de referência, mantendo uma precisão aceitável. Uma versão de "baixa fidelidade" de traçador passivo alcançou uma aceleração de $3 \times 10^5$ com um erro uma ordem de magnitude maior.
• Escalabilidade: Para problemas estáticos com até $10^4$ gotas, os custos de treinamento e avaliação do substituto escalam favoravelmente. Embora o treinamento exiba escalonamento quadrático ($O(M^2)$) devido às interações par a par, o custo de avaliação é linear ($O(M)$) quando o número de pontos de avaliação é moderado. O método permanece competitivo com o CutFEM para $M < 5000$ ao incluir custos de avaliação, e supera significativamente o CutFEM quando apenas um número limitado de avaliações (por exemplo, para advecção de marcadores em simulações no domínio do tempo) é necessário.
• Robustez: O esquema iterativo demonstrou robustez mesmo em configurações não físicas com gotas sobrepostas, convergindo onde métodos padrão poderiam falhar ou exigir refinamento excessivo de malha.

Significância e Alegações
Os autores alegam que esta estrutura de modelagem substituta hierárquica oferece um caminho prático para simular fluxos de Stokes multifásicos complexos em estações de trabalho padrão, uma tarefa que tipicamente requer recursos de supercomputação para métodos de alta fidelidade. O método é particularmente significativo para cenários envolvendo muitas gotas e pequenas separações interfaciais, onde as técnicas de discretização padrão enfrentam gargalos computacionais severos. O artigo enfatiza que a abordagem captura sistematicamente as interações gota-fluxo, gota-fronteira e gota-gota. Embora a implementação atual assuma gotas circulares, os autores observam que a estrutura é suficientemente rica para motivar extensões futuras para gotas deformáveis e potencialmente fluxos de Navier-Stokes. O trabalho demonstra que, ao combinar insights físicos (PDEs de erro) com álgebra linear numérica eficiente (aproximações de subespaço), pode-se alcançar um equilíbrio entre precisão e tratabilidade computacional anteriormente inalcançável para sistemas de gotas em grande escala.