En ny metod för att identifiera fria ytor i 3D MPS-metoden

Varför det är viktigt att följa rörliga vattensytor

Från vågor som slår mot sjöförsvar till bränsleskvall inne i rymdfarkosters tankar, förlitar sig ingenjörer i allt större utsträckning på datorsimuleringar för att förutsäga hur vätskor rör sig och vilka krafter de utövar. Många av dessa simuleringar använder en svärm av virtuella partiklar istället för ett styvt rutnät för att representera vatten. Men det finns ett subtilt problem: datorn måste ständigt avgöra vilka partiklar som ligger på vätskans yttre yta, där omgivande luftens tryck verkar. Om det beslutet är ens lite felaktigt kan de beräknade trycken på väggar och konstruktioner bli brusiga eller missvisande. Denna studie presenterar ett nytt sätt att hitta dessa ytplaster i tredimensionella simuleringar, vilket får det virtuella vattnet att bete sig mer som det verkliga.

Ett annat sätt att dra gränsen mellan vatten och luft

I Moving Particle Semi-Implicit (MPS)-metoden modelleras vatten som otaliga små partiklar som rör sig och interagerar. Traditionella tillvägagångssätt avgör om en partikel sitter på den fria ytan genom att helt enkelt räkna hur många grannar den har. Färre grannar betyder vanligen att den ligger vid eller nära ytan. Denna tumregel kan dock misslyckas när partiklar är ojämnt fördelade, vilket ofta sker i turbulenta flöden eller nära skarpa hörn. Författarna fokuserar istället på hur närliggande partiklar är ordnade i rummet, med hjälp av en kvantitet kallad relativ positionsdivergens, som mäter hur "utbredda" grannarna är kring en given partikel. Om grannarna tenderar att ligga mest på ena sidan förändras divergensevärdet på ett sätt som signalerar en sannolik yta.

Lägga in riktning i yttestet

Även om den förbättrade divergensmåttet väl överensstämmer med teoretiska förväntningar visar författarna att i tre dimensioner kan en enkel tröskel på detta värde inte perfekt skilja inre och yttre partiklar, särskilt nära skarpa kanter och inuti håligheter. För att åtgärda detta introducerar de ett extra steg som uppskattar den lokala ytriktningen, eller normalvektorn, från de omgivande partiklarna. Med denna riktning definierar metoden ett konformat område som pekar utåt från varje kandidatpartikel. Om den konen inte innehåller några andra partiklar inom ett visst avstånd bedöms partikeln ligga på den fria ytan. Detta kombinerade tillvägagångssätt, kallat RPD+NV, använder både hur grannarna är ordnade och var tomrum finns, vilket ger en mer pålitlig bild av den verkliga gränsytan mellan vätska och luft.

Testning av metoden på komplexa former och lugnt vatten

Forskarlaget testar först sin metod på stillastående, rent geometriska former: en kub krönt av en kupol, en kub med en skålformad fördjupning och ett S-format rör som innehåller en inre hålighet. Dessa fall utmanar algoritmen med både utbuktande och ihåliga ytor. Genom att jämföra de upptäckta ytplasterna med de kända verkliga ytorna visar de att den nya metoden korrekt fångar både konvexa och konkava regioner och noggrant löser ut dolda håligheter. De tillämpar sedan metoden på en enkel hydrostatisk vattentank och verifierar att, med rimliga val av partikelavstånd och tidssteg, det beräknade trycket vid botten mycket väl överensstämmer med det teoretiska värdet. Under dessa lugna förhållanden ger den nya ytdetektorn färre tryckoscillationer än äldre metoder, samtidigt som den behåller nästan samma beräkningskostnad.

Fånga vågor, kollisioner och skvalpning

Därefter utsätter teamet metoden för högdynamiska situationer. I en tredimensionell dammbrottssimulering kollapsar en vattenpelare och slår mot en vägg, vilket skapar snabbt föränderliga fria ytor, stänk och luftfickor. Jämfört med traditionella kriterier baserade enbart på grannantal eller divergens spårar RPD+NV-metoden mer kontinuerliga och detaljerade ytor och ger jämnare, mer realistiska tryckfält. Avgörande är att toppen av kollisionstrycket mot väggen kommer mycket närmare laboratoriemätningarna. I ett dämpat skvalpningsexempel, där vågor inuti en tank gradvis lägger sig, leder den nya metoden till den mest jämna utvecklingen av krafter på väggarna, vilket indikerar minskat numeriskt brus. Slutligen, i ett test som efterliknar bränsleskvallning i en ringformad tank som gungar fram och tillbaka, överensstämmer ytformer och trycktoppar i simuleringen väl med bilder från höghastighetskamera och trycksensordata, återigen utan falska tryckspikar vid ytan.

Vad detta innebär för verkliga simuleringar

För läsare utanför beräkningsfluiddynamik är den praktiska slutsatsen att detta arbete förbättrar hur datorer avgör var vatten slutar och luft börjar i partikelbaserade modeller. Genom att kombinera ett förfinat mått på hur grannar omger varje partikel med ett riktat "tomkon"-test, reducerar den nya RPD+NV-metoden kraftigt felaktiga klassificeringar av ytplaster. Resultatet är renare, mer stabila tryckprediktioner över en rad flöden, från stilla vatten till våldsamma kollisioner och komplex tredimensionell skvalpning. Eftersom metoden tillför liten beräkningskostnad erbjuder den en praktisk uppgradering för ingenjörer som simulerar kuststrukturer, fartygsrörelser, lagringstankar och andra system där en noggrann fångst av rörliga vattensytor är avgörande.

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

Nyckelord: partikelbaserad vätskesimulering, detektion av fri yta, numerisk vågpåverkan, skvalpning i tankar, moving particle method