摘要
传统连续表面力(continuum surface force,CSF)模型模拟流体表面张力时,流体表面的粒子受到的表面张力都是法向力,很难保证流体表面的光滑性,从而使得表面张力模拟失真;此外,流体表面的粒子不足还导致流体表面密度计算精度降低,模拟稳定性差。文章在CSF模型的基础上,对流体表面施加切向力,使得流体表面更加光滑;对边界粒子进行密度修正,提高了密度计算精度和模拟的稳定性。仿真实验结果表明,该方法模拟的流体表面张力效果更好。
When the traditional continuum surface force(CSF)model simulates the fluid surface tension,the surface tension of particles on the fluid surface is normal force,and it is difficult to guarantee the smoothness of the fluid surface,thus making the surface tension simulation distorted.In addition,the lack of particles on the fluid surface leads to lower computational accuracy of the surface density and poor stability of the simulation.In this paper,based on the CSF model,the tangential force is applied to the fluid surface to make the fluid surface smoother.Then the density correction of boundary particles is conducted to improve the computational accuracy of density and the stability of the simulation.Simulation results show that the proposed method performs better in simulating the fluid surface tension.
引文
[1] NUGENT S,POSCH H A.Liquid drops and surface tension with smoothed particle applied mechanics[J].Physical Review E,2000,62:4968-4975.
[2] MELEAN Y,SIGALOTTI L D G.Coalescence of colliding van der Waals liquid drops[J].International Journal of Heat and Mass Transfer,2005,48:4041-4061.
[3] TARTAKOVSKY A,MEAKIN P.Modeling of surface tension and contact angles with smoothed particle hydrodynamics[J].Physical Review E,2005,72(2):026301
[4] ZHOU G Z,GE W,LI J H.A revised surface tension model for macro-scale particle methods[J].Powder Technology,2008,183:21-26.
[5] AKINCI N,AKINCI G,TESCHNER M.Versatile surface tension and adhension for SPH fluids[J].ACM Transactions on Graphics,2013,32(6):1-8.
[6]董兰芳,章恒,陈骏雄.基于IIF模型的流体表面张力模拟[J].图学学报,2016,37(3):302-307.
[7]MORRIS J P.Simulating surface tension with smoothed particle hydrodynamics[J].International Journal for Numerical Methods in Fluids,2000,33:333-353.
[8] HU X Y,ADAMS N.Multi-phase SPH method for macroscopic and mesoscopic flows[J].Journal of Computational Physics,2006,213(2):844-861.
[9] FANG H S,BAO K,WEI J A,et al.Simulations of droplet spreading and solidification using an improved SPH model[J].Numerical Heat Transfer,Part A:Applications,2009,55:124-143.
[10] BREINLINGER T,POLFER P,HASHIBON A,et al.Surface tension and wetting effects with smoothed particle hydrodynamics[J].Journal of Computational Physics,2013,243(12):14-27.
[11]强洪夫,陈福振,高巍然.修正表面张力算法的SPH方法及其实现[J].计算物理,2011,28(3):375-384.
[12] LUCY L B.A numerical approach to the testing of the fission hypothesis[J].The Astronomical Journal,1977,82:1013-1024.
[13]GINGOLD R A,MONAGHAN J J.Smoothed particle hydro-dynamics:theory and application to non-spherical stars[J].Monthly Notices of the Royal Astronomical Society,1977,181(3):375-389.
[14] DESBRUN M,GASCUEL M P.Smoothed particles:a new paradigm for animating highly deformable bodies[C]//EGCA:Proceedings of the 1996Eurographics Workshop on Computer Animation and Simulation.Berlin:Springer,1996:61-76.
[15] MU··LLER M,CHARYPAR D,GROSS M.Particle-based fluid simulation for interactive applications[C]//Proceedings of the ACM SIGGRAPH/Eurographics Symposium on Computer Animation.Aire-la-Ville:Eurographics Association Press,2003:154-159.
[16] RANDLES P W,LIBERSKY L D.Smoothed particle hydrodynamics:some recent improvements and applications[J].Computer Methods in Applied Mechanics and Engineering,1996,139(1/2/3/4):375-408.