欢迎访问《图学学报》 分享到:

图学学报 ›› 2020, Vol. 41 ›› Issue (6): 970-979.DOI: 10.11996/JG.j.2095-302X.2020060970

• 计算机图形学与虚拟现实 • 上一篇    下一篇

修正压力与表面张力计算的两相流 自由界面运动模拟 

  

  1. (合肥工业大学数学学院,安徽 合肥 230009)
  • 出版日期:2020-12-31 发布日期:2021-01-08
  • 基金资助:
    基金项目:国家自然科学基金项目(61272024) 

Simulation of free interface motion of two-phase flow based on modified pressure and surface tension calculation 

  1. (School of Mathematics, Hefei University of Technology, Hefei Anhui 230009, China) 
  • Online:2020-12-31 Published:2021-01-08
  • Supported by:
    Foundation items:National Natural Science Foundation of China (61272024) 

摘要: 摘 要:在利用光滑粒子流体幼力学(SPH)方法进行多相流的模拟中,多相流在交界面处由 于密度不连续、粒子界面的压力计算出现误差,出现界面处压力振荡,界面破碎的问题。针对 上述问题,提出了新的压力梯度近似公式和改进的界面处人工斥力公式,使交界面更加清晰、 光滑、无穿透、模拟效果更好。此外,通过给出基于密度权重的色函数计算公式,对大密度比 多相流界面处的表面张力计算公式进行了改进,使得多相流交界面密度过渡更加光滑,模拟效 果更好。最后,通过溃坝、Rayleigh-Taylor 界面不稳定性,非 Boussinesq 锁定交换问题 3 个仿 真实验,得到了不同时刻界面粒子分布图、界面锋面距离等。结果验证了新的压力梯度近似公 式以及界面处人工斥力公式的合理性;通过空气中液滴形成仿真实验,得到了圆形液滴形成的 粒子变化图,结果验证了该改进的表面张力计算方法的有效性。

关键词: 关 键 词:SPH 方法, 多相流, 压力梯度近似公式, 人工斥力, 色函数, 表面张力

Abstract: Abstract: In the simulation of multiphase flow using the smoothed particle hydrodynamics (SPH) method, the density of the multiphase flow at the interface was discontinuous, the pressure calculation of the particle interface was erroneous, and problems occurred, such as pressure oscillation at the interface and interface fracture. In response to the above problems, a new pressure gradient approximation formula and an improved artificial repulsion formula at the interface were proposed to make the interface clearer and smoother without penetration, thus producing better simulation results. In addition, by giving a color function calculation formula based on density weights, the formula was improved for the calculation of surface tension at the interface of multiphase flow with high density ratio, leading to the smoother density transition of multiphase flow interface and better simulation effect. Finally, through three simulation experiments, such as dam break, Rayleigh-Taylor interface instability, and non-Boussinesq locked exchange problem, the particle distribution map of the interface and the distance of the interface front at different times were obtained. The results can verify the rationality of the new pressure gradient approximation formula and the artificial repulsion formula at the interface. Through the simulation experiment of droplet formation in the air, the particle change diagram formed by the circular droplet was obtained, and the result can corroborate the effectiveness of the improved surface tension calculation method in this paper

Key words: Keywords: smoothed particle hydrodynamics method, multi-phase flow, pressure gradient particle approximation, artificial repulsion, the color function, surface Tension

中图分类号: