Um novo método de identificação de superfície livre para o método MPS 3D

Por que rastrear superfícies móveis de água é importante

De ondas que se chocam contra quebra-mares ao combustível que se movimenta dentro de tanques de naves espaciais, engenheiros cada vez mais dependem de simulações computacionais para prever como líquidos se movem e quais forças exercem. Muitas dessas simulações usam um enxame de partículas virtuais em vez de uma grade rígida para representar a água. Mas há um problema sutil: o computador precisa decidir constantemente quais partículas estão na superfície externa do líquido, onde atua a pressão do ar ambiente. Se essa decisão estiver minimamente incorreta, as pressões previstas sobre paredes e estruturas podem ficar ruidosas ou enganosas. Este estudo introduz uma nova maneira de encontrar essas partículas de superfície em simulações tridimensionais, fazendo com que a água virtual se comporte de forma mais parecida com a real.

Uma maneira diferente de traçar a linha entre água e ar

No método Moving Particle Semi-Implicit (MPS), a água é modelada como inúmeras pequenas partículas que se movem e interagem. Abordagens tradicionais decidem se uma partícula está na superfície livre simplesmente contando quantos vizinhos ela tem. Menos vizinhos geralmente indica que está na superfície ou próxima dela. Contudo, essa regra prática pode falhar quando as partículas estão distribuídas de forma desigual, como acontece frequentemente em escoamentos turbulentos ou perto de cantos agudos. Os autores, em vez disso, concentram-se em como as partículas vizinhas estão arranjadas no espaço, usando uma grandeza chamada divergência de posição relativa, que mede o quão “espalhados” estão os vizinhos ao redor de uma dada partícula. Se os vizinhos tendem a se situar majoritariamente de um lado, o valor da divergência muda de modo a sinalizar uma provável superfície.

Adicionando direção ao teste de superfície

Embora a medida de divergência aprimorada corresponda muito bem às expectativas teóricas, os autores mostram que, em três dimensões, um simples limiar nesse valor não consegue separar perfeitamente partículas interiores e de superfície, especialmente perto de arestas afiadas e dentro de cavidades. Para corrigir isso, eles introduzem uma etapa extra que estima a direção local da superfície, ou vetor normal, a partir das partículas circundantes. Usando essa direção, o método define uma região em forma de cone que aponta para fora de cada partícula candidata. Se esse cone não contiver outras partículas dentro de certa distância, a partícula é classificada como pertencente à superfície livre. Essa abordagem combinada, chamada RPD+NV, usa tanto a disposição dos vizinhos quanto onde há espaço vazio, oferecendo uma imagem mais confiável da interface verdadeira entre líquido e ar.

Testando o método em formas complexas e água em repouso

Os pesquisadores primeiro testam seu método em formas puramente geométricas e estáticas: um cubo com uma cúpula no topo, um cubo com um rebaixo em forma de tigela e um tubo em forma de S contendo uma cavidade interna. Esses casos desafiam o algoritmo com superfícies tanto convexas quanto côncavas. Ao comparar as partículas de superfície detectadas com as superfícies verdadeiras conhecidas, mostram que o novo método captura corretamente regiões convexas e côncavas e resolve com precisão cavidades escondidas. Em seguida, aplicam a abordagem a um recipiente hidrostático simples e verificam que, com escolhas razoáveis de espaçamento de partículas e passo de tempo, a pressão prevista no fundo corresponde muito de perto ao valor teórico. Nessas condições calmas, o novo detector de superfície produz menos oscilações de pressão do que métodos antigos, mantendo o custo computacional quase inalterado.

Capturando ondas, impactos e oscilações

Em seguida, a equipe desafia o método com situações altamente dinâmicas. Em uma simulação 3D de ruptura de barragem, uma coluna de água colapsa e atinge uma parede, criando superfícies livres rapidamente variáveis, respingos e bolsões de ar. Comparado com critérios tradicionais baseados apenas na contagem de vizinhos ou na divergência, o método RPD+NV acompanha superfícies mais contínuas e detalhadas e produz campos de pressão mais suaves e realistas. Crucialmente, a pressão de impacto máxima na parede aproxima-se muito mais das medições de laboratório. Em um caso de oscilações amortecidas, onde ondas dentro de um tanque se assentam gradualmente, o novo método leva à evolução mais suave das forças sobre as paredes, indicando redução do ruído numérico. Finalmente, em um teste que simula o movimento de combustível dentro de um tanque anelar que oscila para frente e para trás, as formas de superfície e os picos de pressão da simulação concordam bem com imagens de alta velocidade e dados de sensores de pressão, novamente sem picos espúrios de pressão na superfície.

O que isso significa para simulações do mundo real

Para leitores fora da dinâmica de fluidos computacional, a mensagem prática é que este trabalho melhora a forma como computadores decidem onde a água termina e o ar começa em modelos baseados em partículas. Ao combinar uma medida refinada de como os vizinhos cercam cada partícula com um teste direcional de “cone vazio”, o novo método RPD+NV reduz drasticamente classificações equivocadas de partículas de superfície. O resultado é previsões de pressão mais limpas e estáveis em uma variedade de escoamentos, desde água parada até impactos violentos e oscilações tridimensionais complexas. Como o método acrescenta pouco custo computacional, oferece uma atualização prática para engenheiros que simulam estruturas costeiras, movimento de navios, tanques de armazenamento e outros sistemas em que capturar com precisão superfícies móveis de água é essencial.

Citação: Geng, C., Wang, Wh., Heng, My. et al. A new free surface identification method for 3D MPS method. Sci Rep 16, 13829 (2026). https://doi.org/10.1038/s41598-026-44218-9

Palavras-chave: simulação de fluidos baseada em partículas, detecção de superfície livre, impacto de ondas numérico, oscilações em tanques, método de partículas móveis