Une nouvelle méthode d’identification de surface libre pour la méthode MPS 3D

Pourquoi il est important de suivre les surfaces d’eau en mouvement

Des vagues qui s’écrasent contre des murs de mer au balancement de carburant à l’intérieur des réservoirs spatiaux, les ingénieurs s’appuient de plus en plus sur des simulations numériques pour prédire les mouvements des liquides et les forces qu’ils exercent. Nombre de ces simulations représentent l’eau par un essaim de particules virtuelles plutôt que par une grille fixe. Mais il existe un problème subtil : l’ordinateur doit constamment déterminer quelles particules se trouvent à la surface externe du liquide, où la pression de l’air ambiant s’applique. Si cette décision est même légèrement erronée, les pressions prédites sur les parois et les structures peuvent devenir bruyantes ou trompeuses. Cette étude introduit une nouvelle manière d’identifier ces particules de surface dans des simulations tridimensionnelles, rendant l’eau virtuelle plus conforme au comportement réel.

Une autre manière de tracer la frontière entre eau et air

Dans la méthode Moving Particle Semi-Implicit (MPS), l’eau est modélisée par d’innombrables petites particules qui se déplacent et interagissent. Les approches traditionnelles décident si une particule est à la surface libre en comptant simplement son nombre de voisines. Moins de voisines signifie généralement qu’elle est à la surface ou à proximité. Cependant, cette règle empirique peut échouer lorsque les particules sont inégalement espacées, comme c’est souvent le cas dans les écoulements turbulents ou près d’angles vifs. Les auteurs se concentrent plutôt sur la façon dont les particules voisines sont disposées dans l’espace, en utilisant une quantité appelée divergence de position relative, qui mesure dans quelle mesure les voisines sont « réparties » autour d’une particule donnée. Si les voisines ont tendance à se situer majoritairement d’un seul côté, la valeur de divergence varie d’une manière qui signale une probable surface.

Ajouter une direction au test de surface

Bien que la mesure de divergence améliorée concorde très bien avec les attentes théoriques, les auteurs montrent qu’en trois dimensions un simple seuil sur cette valeur ne peut pas séparer parfaitement les particules internes et de surface, notamment près d’angles vifs et à l’intérieur de cavités. Pour remédier à cela, ils introduisent une étape supplémentaire qui estime la direction locale de la surface, ou vecteur normal, à partir des particules environnantes. En utilisant cette direction, la méthode définit une région conique qui pointe vers l’extérieur à partir de chaque particule candidate. Si ce cône ne contient aucune autre particule dans une certaine distance, la particule est jugée située sur la surface libre. Cette approche combinée, appelée RPD+NV, utilise à la fois la disposition des voisines et l’emplacement de l’espace vide, offrant une image plus fiable de l’interface réelle entre le liquide et l’air.

Tester la méthode sur des formes complexes et de l’eau au repos

Les chercheurs testent d’abord leur méthode sur des formes purement géométriques et immobiles : un cube surmonté d’un dôme, un cube avec un creux en forme de coupe, et un tube en S contenant une cavité interne. Ces cas mettent l’algorithme au défi avec des surfaces à la fois bombées et creusées. En comparant les particules de surface détectées avec les surfaces vraies connues, ils montrent que la nouvelle méthode capture correctement les régions convexes et concaves et résout avec précision les cavités cachées. Ils appliquent ensuite l’approche à un conteneur hydrostatique simple et vérifient qu’avec des choix raisonnables d’espacement des particules et de pas de temps, la pression prédite au fond correspond très étroitement à la valeur théorique. Dans ces conditions calmes, le nouveau détecteur de surface produit moins d’oscillations de pression que les méthodes plus anciennes, tout en conservant un coût de calcul quasi identique.

Capturer vagues, impacts et sloshing

Ensuite, l’équipe met la méthode à l’épreuve dans des situations très dynamiques. Dans une simulation tridimensionnelle de rupture de digue, une colonne d’eau s’effondre et percute un mur, générant des surfaces libres qui évoluent rapidement, des éclaboussures et des poches d’air. Comparée aux critères traditionnels basés uniquement sur le nombre de voisines ou la divergence, la méthode RPD+NV suit des surfaces plus continues et détaillées et produit des champs de pression plus lisses et plus réalistes. Surtout, la pression d’impact maximale sur le mur se rapproche beaucoup des mesures en laboratoire. Dans un cas de sloshing amorti, où les vagues à l’intérieur d’un réservoir se calment progressivement, la nouvelle méthode conduit à l’évolution la plus régulière des forces sur les parois, indiquant une réduction du bruit numérique. Enfin, dans un test imitant le balancement de carburant à l’intérieur d’un réservoir annulaire qui oscille, les formes de surface et les pics de pression obtenus en simulation concordent bien avec des images haute vitesse et des données de capteurs de pression, là encore sans pics de pression parasites à la surface.

Ce que cela signifie pour les simulations du monde réel

Pour les lecteurs extérieurs à la dynamique numérique des fluides, le message pratique est que ce travail améliore la façon dont les ordinateurs décident où l’eau s’arrête et où commence l’air dans les modèles basés sur des particules. En combinant une mesure raffinée de la façon dont les voisines entourent chaque particule avec un test directionnel du « cône vide », la nouvelle méthode RPD+NV réduit fortement les classifications erronées des particules de surface. Le résultat est des prédictions de pression plus propres et plus stables pour une gamme d’écoulements, de l’eau au repos aux impacts violents et aux mouvements complexes tridimensionnels. Comme la méthode ajoute peu de coût informatique, elle constitue une amélioration pratique pour les ingénieurs simulant des structures côtières, le comportement des navires, des réservoirs de stockage et d’autres systèmes où la capture précise des surfaces d’eau en mouvement est essentielle.

Citation: 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

Mots-clés: simulation de fluides par particules, détection de surface libre, impact de vagues numérique, sloshing dans les réservoirs, méthode des particules mobiles