An improved weakly compressible SPH method for simulating 2D multiphase flows with complex interface and large density ratios

IF 4.2 2区工程技术 Q1 ENGINEERING, MULTIDISCIPLINARY

Engineering Analysis with Boundary Elements Pub Date : 2024-11-26 DOI:10.1016/j.enganabound.2024.106043

Xiaoyang Xu , Erdi Wang , Xiang Li , Peng Yu

{"title":"An improved weakly compressible SPH method for simulating 2D multiphase flows with complex interface and large density ratios","authors":"Xiaoyang Xu , Erdi Wang , Xiang Li , Peng Yu","doi":"10.1016/j.enganabound.2024.106043","DOIUrl":null,"url":null,"abstract":"<div><div>This paper develops an improved weakly compressible smoothed particle hydrodynamics (SPH) method for simulating multiphase flows with complex interface and large density ratios. Surface tension is computed using a continuum surface force method along with a kernel gradient correction algorithm, thereby improving the numerical precision of normal vectors and curvatures. To maintain a uniform particle distribution and prevent instabilities resulting from particle stretching, a particle shifting technique is implemented. Additionally, an anisotropic interfacial repulsion force is introduced at the multiphase interface to create smooth phase boundaries and mitigate issues related to particle penetration. The efficacy of the method is demonstrated through simulations of various scenarios including the deformation of a square droplet, the Rayleigh-Taylor instability, the lock-exchange phenomena, a single rising bubble, and two rising bubbles. By comparing the results with analytical solutions and existing literature data, the accuracy and consistency of the numerical results are confirmed. These comparisons demonstrate that the improved multiphase SPH method developed in this study can effectively track deformable phase interfaces and simulate multiphase flows with complex interfaces and large density ratios in a stable way.</div></div>","PeriodicalId":51039,"journal":{"name":"Engineering Analysis with Boundary Elements","volume":"169 ","pages":"Article 106043"},"PeriodicalIF":4.2000,"publicationDate":"2024-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Engineering Analysis with Boundary Elements","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0955799724005162","RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"ENGINEERING, MULTIDISCIPLINARY","Score":null,"Total":0}

引用次数: 0

Abstract

This paper develops an improved weakly compressible smoothed particle hydrodynamics (SPH) method for simulating multiphase flows with complex interface and large density ratios. Surface tension is computed using a continuum surface force method along with a kernel gradient correction algorithm, thereby improving the numerical precision of normal vectors and curvatures. To maintain a uniform particle distribution and prevent instabilities resulting from particle stretching, a particle shifting technique is implemented. Additionally, an anisotropic interfacial repulsion force is introduced at the multiphase interface to create smooth phase boundaries and mitigate issues related to particle penetration. The efficacy of the method is demonstrated through simulations of various scenarios including the deformation of a square droplet, the Rayleigh-Taylor instability, the lock-exchange phenomena, a single rising bubble, and two rising bubbles. By comparing the results with analytical solutions and existing literature data, the accuracy and consistency of the numerical results are confirmed. These comparisons demonstrate that the improved multiphase SPH method developed in this study can effectively track deformable phase interfaces and simulate multiphase flows with complex interfaces and large density ratios in a stable way.

查看原文本刊更多论文

用于模拟具有复杂界面和大密度比的二维多相流的改进型弱可压缩 SPH 方法

本文开发了一种改进的弱可压缩平滑粒子流体力学（SPH）方法，用于模拟具有复杂界面和大密度比的多相流。表面张力的计算采用连续表面力法和核梯度修正算法，从而提高了法向量和曲率的数值精度。为了保持颗粒分布均匀，防止颗粒拉伸造成不稳定，采用了颗粒移动技术。此外，还在多相界面上引入了各向异性的界面斥力，以创建平滑的相界并缓解与粒子穿透有关的问题。通过模拟各种情况，包括方形液滴的变形、瑞利-泰勒不稳定性、锁定交换现象、单个上升气泡和两个上升气泡，证明了该方法的有效性。通过将结果与分析解和现有文献数据进行比较，证实了数值结果的准确性和一致性。这些比较结果表明，本研究开发的改进型多相 SPH 方法可以有效地跟踪可变形的相界面，稳定地模拟具有复杂界面和大密度比的多相流。

本文章由计算机程序翻译，如有差异，请以英文原文为准。

求助全文

约1分钟内获得全文求助全文

来源期刊

Engineering Analysis with Boundary Elements 工程技术-工程：综合

CiteScore

5.50

自引率

18.20%

发文量

368

审稿时长

56 days

期刊介绍： This journal is specifically dedicated to the dissemination of the latest developments of new engineering analysis techniques using boundary elements and other mesh reduction methods. Boundary element (BEM) and mesh reduction methods (MRM) are very active areas of research with the techniques being applied to solve increasingly complex problems. The journal stresses the importance of these applications as well as their computational aspects, reliability and robustness. The main criteria for publication will be the originality of the work being reported, its potential usefulness and applications of the methods to new fields. In addition to regular issues, the journal publishes a series of special issues dealing with specific areas of current research. The journal has, for many years, provided a channel of communication between academics and industrial researchers working in mesh reduction methods Fields Covered: • Boundary Element Methods (BEM) • Mesh Reduction Methods (MRM) • Meshless Methods • Integral Equations • Applications of BEM/MRM in Engineering • Numerical Methods related to BEM/MRM • Computational Techniques • Combination of Different Methods • Advanced Formulations.