Een nieuwe methode voor het identificeren van vrije oppervlakken in de 3D MPS-methode

Waarom het volgen van bewegende wateroppervlakken ertoe doet

Van brekende golven tegen zeeweringen tot brandstofdie heen en weer schuift in ruimtetanks: ingenieurs vertrouwen steeds vaker op computersimulaties om te voorspellen hoe vloeistoffen bewegen en welke krachten ze uitoefenen. Veel van deze simulaties gebruiken een zwerm virtuele deeltjes in plaats van een vast rekenrooster om water te modelleren. Er is echter een subtiel probleem: de computer moet continu bepalen welke deeltjes zich op het buitenste oppervlak van de vloeistof bevinden, waar de druk van de omringende lucht geldt. Als die beslissing ook maar iets onnauwkeurig is, kunnen de voorspelde drukken op wanden en constructies ruisig of misleidend worden. Deze studie introduceert een nieuwe manier om die oppervlakte­deeltjes in driedimensionale simulaties te vinden, waardoor virtueel water zich meer als echt water gedraagt.

Een andere manier om de grens tussen water en lucht te trekken

In de Moving Particle Semi-Implicit (MPS)-methode wordt water gemodelleerd als talloze kleine deeltjes die bewegen en op elkaar inwerken. Traditionele benaderingen bepalen of een deeltje op het vrije oppervlak ligt door simpelweg te tellen hoeveel buren het heeft. Minder buren betekent meestal dat het zich op of nabij het oppervlak bevindt. Deze vuistregel kan echter falen wanneer de deeltjes ongelijkmatig verdeeld zijn, zoals vaak voorkomt bij turbulente stromingen of bij scherpe hoeken. De auteurs richten zich in plaats daarvan op hoe naburige deeltjes ruimtelijk zijn gerangschikt, met behulp van een grootheid genaamd relative position divergence (RPD), die meet hoezeer buren 'uitgespreid' liggen rond een gegeven deeltje. Als de buren vooral aan één zijde liggen, verandert de divergentiewaarde op een manier die aangeeft dat het waarschijnlijk een oppervlak betreft.

Richting toevoegen aan de oppervlaktetest

Hoewel de verbeterde divergentiemaat goed overeenkomt met theoretische verwachtingen, laten de auteurs zien dat in drie dimensies een eenvoudige drempel op deze waarde niet perfect binnen- en buiten­deeltjes kan scheiden, vooral bij scherpe randen en in holtes. Om dit te verhelpen introduceren zij een extra stap die de lokale oppervlaktrichting, of normaalvector, schat op basis van de omliggende deeltjes. Met deze richting definieert de methode een kegelvormig gebied dat naar buiten wijst vanaf elk kandidaatdeeltje. Als die kegel binnen een bepaalde afstand geen andere deeltjes bevat, wordt het deeltje beoordeeld als behorend tot het vrije oppervlak. Deze gecombineerde benadering, RPD+NV genoemd, gebruikt zowel de ruimtelijke ordening van buren als de aanwezigheid van lege ruimte en levert zo een betrouwbaardere weergave van de werkelijke grens tussen vloeistof en lucht.

De methode testen op complexe vormen en stilstaand water

De onderzoekers testen hun methode eerst op stilstaande, puur geometrische vormen: een kubus met een koepel bovenop, een kubus met een komvormige uitsparing en een S-vormige buis met een interne holte. Deze gevallen dagen het algoritme uit met zowel uitpuilende als uitgehouwen oppervlakken. Door de gedetecteerde oppervlakte­deeltjes te vergelijken met de bekende echte oppervlakken tonen zij aan dat de nieuwe methode zowel bolle als holle gebieden correct vastlegt en verborgen holtes nauwkeurig resolveert. Vervolgens passen ze de aanpak toe op een eenvoudige hydrostatische watercontainer en verifiëren dat, met redelijke keuzes voor deeltjesafstand en tijdstap, de voorspelde druk op de bodem zeer dicht bij de theoretische waarde ligt. Onder deze kalme omstandigheden veroorzaakt de nieuwe oppervlakte­detector minder drukoscillaties dan oudere methoden, terwijl de rekentijd vrijwel gelijk blijft.

Golven, impacten en sloshing in beeld brengen

Vervolgens dagen de onderzoekers de methode met zeer dynamische situaties uit. In een driedimensionale dambreuksimulatie stort een waterkolom in en slaat tegen een wand, waarbij snel veranderende vrije oppervlakken, opspattend water en luchtkamers ontstaan. Vergeleken met traditionele criteria op basis van enkel het aantal buren of alleen divergentie, volgt de RPD+NV-methode continuïere en gedetailleerdere oppervlakken en levert zij vloeiendere, realistischere drukvelden op. Cruciaal is dat de piekimpactdruk tegen de wand veel dichter bij laboratoriummetingen komt. In een gedempte sloshingcase, waarbij golven in een tank geleidelijk bedaren, zorgt de nieuwe methode voor de meest soepele ontwikkeling van krachten op de wanden, wat wijst op verminderde numerieke ruis. Ten slotte, in een test die brandstofsloshing in een ringvormige tank nabootst die heen en weer beweegt, komen oppervlakte­vormen en drukpieken uit de simulatie goed overeen met beelden van hoge­snelheidscamera's en drukgegevens van sensoren, wederom zonder valse drukpiekjes aan het oppervlak.

Wat dit betekent voor simulaties in de praktijk

Voor lezers buiten de computationele stromingsleer is de praktische boodschap dat dit werk verbetert hoe computers bepalen waar water eindigt en lucht begint in deeltjesgebaseerde modellen. Door een verfijnde maat voor de manier waarop buren elk deeltje omringen te combineren met een richtinggebonden "lege kegel"-test, vermindert de nieuwe RPD+NV-methode sterk het aantal foutieve classificaties van oppervlakte­deeltjes. Het resultaat is schonere, stabielere drukvoorspellingen over een reeks stromingen, van stilstaand water tot hevige impacten en complexe driedimensionale sloshingbewegingen. Omdat de methode weinig extra rekentijd vergt, biedt zij een praktische upgrade voor ingenieurs die kusterosie, scheepsbewegingen, opslagtanks en andere systemen simuleren waarbij het nauwkeurig vastleggen van bewegende wateroppervlakken essentieel is.

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

Trefwoorden: deeltjesgebaseerde vloei­gheidssimulatie, detectie van vrij oppervlak, numerieke golfimpact, sloshing in tanks, moving particle method