{"title":"A WCSPH-MSDEM model for combined rolling and sliding motions of complex-shaped blocks in an unsteady flow","authors":"Sen Gao , Bing Ren , Pengzhi Lin , Ping Dong","doi":"10.1016/j.jcp.2025.114083","DOIUrl":null,"url":null,"abstract":"<div><div>A coupled Weakly-Compressible Smoothed Particle Hydrodynamics (WCSPH) and Multi-Sphere Discrete Element Method (MSDEM) numerical model is employed to investigate the interaction between a discrete, complex-shaped block and a solid boundary in unsteady flows. The study focuses on enhancing the original MSDEM by developing a virtual surface algorithm that assigns a surface normal vector and a virtual surface to each boundary sphere, thereby allowing for a more accurate representation of the true boundary geometry. This approach ensures a precise description of normal contact between the block and the boundary while preventing the embedding phenomenon during the block’s sliding motion. To accurately simulate frictional behavior, both static friction theory and Coulomb’s law of kinetic friction are incorporated. A series of test cases are conducted to validate the modified MSDEM, with further validation of the coupled WCSPH-MSDEM model against experimental data concerning with the sliding motion of a hollow square under solitary waves. Finally, a numerical case study on combined rolling, sliding, and collision of a hollow square on a slope under plunging waves is presented, highlighting the stability and robustness of the proposed model.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"536 ","pages":"Article 114083"},"PeriodicalIF":3.8000,"publicationDate":"2025-05-14","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/S0021999125003663","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
A coupled Weakly-Compressible Smoothed Particle Hydrodynamics (WCSPH) and Multi-Sphere Discrete Element Method (MSDEM) numerical model is employed to investigate the interaction between a discrete, complex-shaped block and a solid boundary in unsteady flows. The study focuses on enhancing the original MSDEM by developing a virtual surface algorithm that assigns a surface normal vector and a virtual surface to each boundary sphere, thereby allowing for a more accurate representation of the true boundary geometry. This approach ensures a precise description of normal contact between the block and the boundary while preventing the embedding phenomenon during the block’s sliding motion. To accurately simulate frictional behavior, both static friction theory and Coulomb’s law of kinetic friction are incorporated. A series of test cases are conducted to validate the modified MSDEM, with further validation of the coupled WCSPH-MSDEM model against experimental data concerning with the sliding motion of a hollow square under solitary waves. Finally, a numerical case study on combined rolling, sliding, and collision of a hollow square on a slope under plunging waves is presented, highlighting the stability and robustness of the proposed model.
期刊介绍:
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.