{"title":"Higher-order MPS models and higher-order Explicit Incompressible MPS (EI-MPS) method to simulate free-surface flows","authors":"Tibing Xu , Seiichi Koshizuka , Tsuyoshi Koyama , Toshihide Saka , Osamu Imazeki","doi":"10.1016/j.jcp.2025.113951","DOIUrl":null,"url":null,"abstract":"<div><div>In this study, higher-order spatial models including the gradient model and Laplacian model based on Taylor's series and using their coordinates as coefficients are evaluated by calculating some simple functions and a diffusion problem. The numerical convergence is achieved by the models as the smaller particle distance can calculate more accurate results. By using the models, when the particle distribution is significantly irregular, increasing the search radius can involve more neighboring particles which consequently improves the accuracy. Based on the proposed higher-order models, the higher-order Explicit Incompressible version of the Moving Particle Semi-implicit method (EI-MPS) is developed. The numerical scheme is validated by simulating various free surface flows including the rotation of a fluid square patch, the impact of two identical rectangular fluid patches, oscillating drop under a central force field, a hydrostatic problem, and dam-break flow. The parameters of the particle distance, search radius, and repeated time in the pressure calculation are all examined in the free surface flows. The proposed method can reproduce the free surface variations, kinetic energy, and total energy variation in the violent flows. It can also obtain the hydrostatic pressure achieving numerical convergence. Increasing the search radius can result in larger errors in simulating the hydrostatic pressure. The impacting pressure caused by the dam-break flow is reflected by the method in good agreement with the experimental measurements.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"532 ","pages":"Article 113951"},"PeriodicalIF":3.8000,"publicationDate":"2025-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125002347","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
In this study, higher-order spatial models including the gradient model and Laplacian model based on Taylor's series and using their coordinates as coefficients are evaluated by calculating some simple functions and a diffusion problem. The numerical convergence is achieved by the models as the smaller particle distance can calculate more accurate results. By using the models, when the particle distribution is significantly irregular, increasing the search radius can involve more neighboring particles which consequently improves the accuracy. Based on the proposed higher-order models, the higher-order Explicit Incompressible version of the Moving Particle Semi-implicit method (EI-MPS) is developed. The numerical scheme is validated by simulating various free surface flows including the rotation of a fluid square patch, the impact of two identical rectangular fluid patches, oscillating drop under a central force field, a hydrostatic problem, and dam-break flow. The parameters of the particle distance, search radius, and repeated time in the pressure calculation are all examined in the free surface flows. The proposed method can reproduce the free surface variations, kinetic energy, and total energy variation in the violent flows. It can also obtain the hydrostatic pressure achieving numerical convergence. Increasing the search radius can result in larger errors in simulating the hydrostatic pressure. The impacting pressure caused by the dam-break flow is reflected by the method in good agreement with the experimental measurements.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.