A new free surface identification method for 3D MPS method | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article A new free surface identification method for 3D MPS method Chong Geng, Wen-hua Wang, Mengyuan Heng, Yu Zhao, Hao Yang, Yi Huang This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-7287063/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 17 Mar, 2026 Read the published version in Scientific Reports → Version 1 posted 16 You are reading this latest preprint version Abstract In the Moving Particle Semi-Implicit (MPS) method, free surface identification in three dimensions has traditionally relied on particle number density. However, this approach often leads to misidentification of surface particles, which in turn compromises the accuracy of pressure field calculations. In this study, an improved formulation for calculating the relative position divergence (RPD) is first proposed to better align with theoretical expectations, and a surface identification criterion for three-dimensional cases is established using a model with internal cavities. Based on this, an improved surface particle identification method is introduced by combining RPD with an estimation of the surface normal vector (RPD + NV). Validation through geometry-based models demonstrates that the proposed approach can accurately capture surface features, including sharp corners, concave and convex regions, and internal cavities. Further verification using a series of MPS simulations, including static water, dam break, and damped oscillation scenarios, confirms that both the improved RPD method and the combined RPD + NV scheme can effectively suppress pressure fluctuations. Moreover, the proposed method significantly reduces false negatives and false positives in surface particle identification, resulting in a smoother and more accurate pressure field. Consequently, the proposed surface identification method is well-suited for three-dimensional MPS simulations. Physical sciences/Engineering Physical sciences/Mathematics and computing Physical sciences/Physics MPS method free surface particles relative position divergence normal vector pressure fluctuation Full Text Additional Declarations No competing interests reported. Cite Share Download PDF Status: Published Journal Publication published 17 Mar, 2026 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Revision requested 17 Nov, 2025 Reviews received at journal 13 Oct, 2025 Reviews received at journal 10 Oct, 2025 Reviews received at journal 07 Oct, 2025 Reviews received at journal 05 Oct, 2025 Reviewers agreed at journal 03 Sep, 2025 Reviewers agreed at journal 30 Aug, 2025 Reviewers agreed at journal 29 Aug, 2025 Reviewers agreed at journal 29 Aug, 2025 Reviewers agreed at journal 28 Aug, 2025 Reviewers agreed at journal 28 Aug, 2025 Reviewers agreed at journal 28 Aug, 2025 Reviewers invited by journal 28 Aug, 2025 Editor assigned by journal 14 Aug, 2025 Submission checks completed at journal 11 Aug, 2025 First submitted to journal 11 Aug, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-7287063","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":509855197,"identity":"1b3aa4d2-c40d-4297-895c-872b463fcbf0","order_by":0,"name":"Chong Geng","email":"","orcid":"","institution":"Dalian University of Technology","correspondingAuthor":false,"prefix":"","firstName":"Chong","middleName":"","lastName":"Geng","suffix":""},{"id":509855198,"identity":"c1a5f037-b24b-4e41-8f72-49b94e89c4f1","order_by":1,"name":"Wen-hua Wang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA60lEQVRIiWNgGAWjYBCDBAYG5gMSD2BsIrWwJUgkMBiQpIXHgDgtBsfPHn7NU3Enz+D4mY83Etv+MPCz5xgw/NyBR8uZvDRrnjPPig3O5G62SGwzYJDseWPA2HsGtxazAzlmxrlthxM3HMjdJgHSYnAjx4CZsQ2PlvNvgFr+AbWcf/MMrMWeoJYbOcaPcxuAWm7ksEFskSCgxf7GGzPmP8cOJ8688czYIuGcMY/EmWcFB3vxaJHszzH+OKPmcGLf+eSHNz6UycnxtydvfPATjxYgYJNA5vGAiAN4NQATygcCCkbBKBgFo2CkAwDmaVlCyPQ49QAAAABJRU5ErkJggg==","orcid":"","institution":"Dalian University of Technology","correspondingAuthor":true,"prefix":"","firstName":"Wen-hua","middleName":"","lastName":"Wang","suffix":""},{"id":509855202,"identity":"d10b1fcd-2455-426d-be0f-aa22e251225b","order_by":2,"name":"Mengyuan Heng","email":"","orcid":"","institution":"Dalian University of Technology","correspondingAuthor":false,"prefix":"","firstName":"Mengyuan","middleName":"","lastName":"Heng","suffix":""},{"id":509855204,"identity":"75f10757-2c17-45ba-9f9b-4dab9fd825bd","order_by":3,"name":"Yu Zhao","email":"","orcid":"","institution":"Dalian University of Technology","correspondingAuthor":false,"prefix":"","firstName":"Yu","middleName":"","lastName":"Zhao","suffix":""},{"id":509855205,"identity":"1f77bb1c-a961-4f8e-8215-54f1529deb24","order_by":4,"name":"Hao Yang","email":"","orcid":"","institution":"Dalian University of Technology","correspondingAuthor":false,"prefix":"","firstName":"Hao","middleName":"","lastName":"Yang","suffix":""},{"id":509855206,"identity":"72782b14-ed8e-4713-b972-c0ae15012aab","order_by":5,"name":"Yi Huang","email":"","orcid":"","institution":"Dalian University of Technology","correspondingAuthor":false,"prefix":"","firstName":"Yi","middleName":"","lastName":"Huang","suffix":""}],"badges":[],"createdAt":"2025-08-04 05:23:16","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-7287063/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-7287063/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-026-44218-9","type":"published","date":"2026-03-17T15:58:33+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":105223582,"identity":"d34fb3fd-ec2c-4738-bda1-9921f114b85e","added_by":"auto","created_at":"2026-03-23 16:08:41","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":3728766,"visible":true,"origin":"","legend":"","description":"","filename":"Anewfreesurfaceidentificationmethodfor3DMPSmethod.pdf","url":"https://assets-eu.researchsquare.com/files/rs-7287063/v1_covered_ca9dc7a4-a207-41af-8736-748ee41d4d1e.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"A new free surface identification method for 3D MPS method","fulltext":[],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":true,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":true,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"
[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"MPS method, free surface particles, relative position divergence, normal vector, pressure fluctuation","lastPublishedDoi":"10.21203/rs.3.rs-7287063/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-7287063/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eIn the Moving Particle Semi-Implicit (MPS) method, free surface identification in three dimensions has traditionally relied on particle number density. However, this approach often leads to misidentification of surface particles, which in turn compromises the accuracy of pressure field calculations. In this study, an improved formulation for calculating the relative position divergence (RPD) is first proposed to better align with theoretical expectations, and a surface identification criterion for three-dimensional cases is established using a model with internal cavities. Based on this, an improved surface particle identification method is introduced by combining RPD with an estimation of the surface normal vector (RPD\u0026thinsp;+\u0026thinsp;NV). Validation through geometry-based models demonstrates that the proposed approach can accurately capture surface features, including sharp corners, concave and convex regions, and internal cavities. Further verification using a series of MPS simulations, including static water, dam break, and damped oscillation scenarios, confirms that both the improved RPD method and the combined RPD\u0026thinsp;+\u0026thinsp;NV scheme can effectively suppress pressure fluctuations. Moreover, the proposed method significantly reduces false negatives and false positives in surface particle identification, resulting in a smoother and more accurate pressure field. Consequently, the proposed surface identification method is well-suited for three-dimensional MPS simulations.\u003c/p\u003e","manuscriptTitle":"A new free surface identification method for 3D MPS method","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-09-04 04:28:09","doi":"10.21203/rs.3.rs-7287063/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-11-17T08:55:17+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-10-13T20:23:34+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-10-10T16:30:10+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-10-07T09:03:32+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-10-05T15:37:07+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"57959664102018310128838399544404339959","date":"2025-09-03T18:31:20+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"121499622475674068873655145408225195656","date":"2025-08-30T14:21:56+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"36787502015595454736133551987213120672","date":"2025-08-29T16:07:17+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"77649645186725978964381364671566877687","date":"2025-08-29T08:58:12+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"248229988477155603597013462170210950011","date":"2025-08-28T13:47:47+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"46947597579932534644280987865722310238","date":"2025-08-28T12:59:27+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"12066870488244644245874582369624248855","date":"2025-08-28T10:18:48+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-08-28T09:49:12+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-08-14T16:34:20+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-08-11T13:46:24+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2025-08-11T13:42:22+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"
[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"19b606f4-4fa6-484f-b338-70c4eab3b27f","owner":[],"postedDate":"September 4th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":54154911,"name":"Physical sciences/Engineering"},{"id":54154912,"name":"Physical sciences/Mathematics and computing"},{"id":54154913,"name":"Physical sciences/Physics"}],"tags":[],"updatedAt":"2026-03-23T16:05:28+00:00","versionOfRecord":{"articleIdentity":"rs-7287063","link":"https://doi.org/10.1038/s41598-026-44218-9","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2026-03-17 15:58:33","publishedOnDateReadable":"March 17th, 2026"},"versionCreatedAt":"2025-09-04 04:28:09","video":"","vorDoi":"10.1038/s41598-026-44218-9","vorDoiUrl":"https://doi.org/10.1038/s41598-026-44218-9","workflowStages":[]},"version":"v1","identity":"rs-7287063","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-7287063","identity":"rs-7287063","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.