Nowa metoda identyfikacji swobodnej powierzchni dla trójwymiarowej metody MPS

Dlaczego śledzenie poruszających się powierzchni wody ma znaczenie

Od uderzeń fal o falochrony po kołysanie paliwa w zbiornikach statków kosmicznych — inżynierowie coraz częściej polegają na symulacjach komputerowych, aby przewidzieć ruchy cieczy i siły, które one wywołują. Wiele z tych symulacji reprezentuje wodę za pomocą chmury wirtualnych cząstek zamiast sztywnej siatki. Istnieje jednak subtelny problem: komputer musi nieustannie decydować, które cząstki znajdują się na zewnętrznej powierzchni cieczy, gdzie działa ciśnienie powietrza. Nawet niewielki błąd w tej decyzji może sprawić, że przewidywane obciążenia ścian i konstrukcji staną się zaszumione lub mylące. W tej pracy wprowadzono nowy sposób wykrywania cząstek powierzchniowych w symulacjach trójwymiarowych, dzięki czemu wirtualna woda zachowuje się bardziej podobnie do rzeczywistej.

Inny sposób wyznaczania granicy między wodą a powietrzem

W metodzie Moving Particle Semi-Implicit (MPS) wodę modeluje się jako niezliczone drobne cząstki, które się poruszają i oddziałują ze sobą. Tradycyjne podejścia decydują, czy cząstka leży na powierzchni swobodnej, po prostu zliczając, ile ma sąsiadów. Mniej sąsiadów zwykle oznacza położenie przy lub na powierzchni. Jednak ta reguła może zawodzić, gdy cząstki są nierównomiernie rozłożone, co często zachodzi w przepływach turbulentnych lub przy ostrych krawędziach. Autorzy zamiast tego skupiają się na przestrzennym ułożeniu sąsiadów, używając wielkości zwanej względną dywergencją położeń, która mierzy, jak bardzo sąsiedzi są „rozproszeni” wokół danej cząstki. Jeśli sąsiedzi mają tendencję do leżenia głównie po jednej stronie, wartość dywergencji zmienia się w sposób sygnalizujący prawdopodobną powierzchnię.

Dodanie kierunku do testu powierzchni

Chociaż ulepszona miara dywergencji bardzo dobrze odpowiada przewidywaniom teoretycznym, autorzy pokazują, że w trzech wymiarach progi ustawione tylko na tej wartości nie są w stanie idealnie oddzielić cząstek wnętrza od cząstek powierzchniowych, szczególnie przy ostrych krawędziach i wewnątrz wnęk. Aby to naprawić, wprowadzają dodatkowy krok, który estymuje lokalny kierunek powierzchni, czyli wektor normalny, na podstawie otaczających cząstek. Korzystając z tego kierunku, metoda definiuje stożkowy obszar skierowany na zewnątrz od każdej kandydującej cząstki. Jeśli w tym stożku nie ma innych cząstek w określonej odległości, cząstka uznawana jest za należącą do powierzchni swobodnej. To połączenie, nazwane RPD+NV, używa zarówno informacji o rozmieszczeniu sąsiadów, jak i o istnieniu pustej przestrzeni, dając bardziej wiarygodny obraz rzeczywistego interfejsu między cieczą a powietrzem.

Testowanie metody na skomplikowanych kształtach i spokojnej wodzie

Badacze najpierw testują swoją metodę na nieruchomych, czysto geometrycznych kształtach: sześcianie z kopułą na górze, sześcianie z misowatym zagłębieniem oraz rurze w kształcie litery S zawierającej wewnętrzną wnękę. Przypadki te angażują algorytm zarówno na wypukłych, jak i wklęsłych powierzchniach. Porównując wykryte cząstki powierzchniowe z prawdziwymi znanymi powierzchniami, pokazują, że nowa metoda poprawnie wychwytuje regiony wypukłe i wklęsłe oraz dokładnie rozpoznaje ukryte wnęki. Następnie stosują podejście do prostego, hydrostatycznego pojemnika z wodą i weryfikują, że przy rozsądnych wyborach rozstawu cząstek i kroku czasowego przewidywane ciśnienie na dnie bardzo dobrze zgadza się z wartością teoretyczną. W tych spokojnych warunkach nowy detektor powierzchni generuje mniej oscylacji ciśnienia niż starsze metody, przy niemal niezmienionych kosztach obliczeniowych.

Rejestrowanie fal, uderzeń i kołysania

Następnie zespół wystawia metodę na próbę w bardzo dynamicznych sytuacjach. W trójwymiarowej symulacji zapory kolumna wody zapada się i uderza w ścianę, tworząc szybko zmieniające się powierzchnie swobodne, bryzgi i kieszenie powietrzne. W porównaniu z tradycyjnymi kryteriami opartymi wyłącznie na liczbie sąsiadów lub samej dywergencji, metoda RPD+NV śledzi bardziej ciągłe i szczegółowe powierzchnie oraz generuje gładsze, bardziej realistyczne pola ciśnienia. Co ważne, maksymalne ciśnienie uderzeniowe przy ścianie jest znacznie bliższe pomiarom laboratoryjnym. W tłumionym przypadku kołysania, gdzie fale wewnątrz zbiornika stopniowo się wygaszają, nowa metoda daje najpłynniejszy przebieg sił działających na ścianki, wskazując na zmniejszenie szumu numerycznego. W końcu, w teście naśladującym kołysanie paliwa wewnątrz pierścieniowego zbiornika poruszającego się tam i z powrotem, kształty powierzchni i piki ciśnienia z symulacji dobrze zgadzają się z obrazami z szybkiej kamery i danymi z czujników ciśnienia, ponownie bez fałszywych skoków ciśnienia na powierzchni.

Co to oznacza dla rzeczywistych symulacji

Dla czytelników spoza dziedziny obliczeniowej mechaniki płynów praktyczny komunikat jest taki: praca ta poprawia sposób, w jaki komputery rozstrzygają, gdzie kończy się woda, a zaczyna powietrze w modelach opartych na cząstkach. Poprzez połączenie udoskonalonej miary otoczenia cząstek z kierunkowym testem „pustego stożka”, nowa metoda RPD+NV znacznie redukuje błędne klasyfikacje cząstek powierzchniowych. Efektem są czystsze, bardziej stabilne prognozy ciśnienia w szerokim zakresie przepływów — od stojącej wody, przez gwałtowne uderzenia, po złożone trójwymiarowe kołysanie. Ponieważ metoda dodaje niewiele kosztu obliczeniowego, stanowi praktyczne ulepszenie dla inżynierów symulujących konstrukcje przybrzeżne, ruch statków, zbiorniki magazynowe i inne systemy, w których dokładne odwzorowanie poruszających się powierzchni wody jest istotne.

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

Słowa kluczowe: symulacja płynów oparta na cząstkach, detekcja powierzchni swobodnej, numeryczne uderzenia fal, kołysanie paliwa w zbiornikach, metoda poruszających się cząstek