摘要
对于近水面水下爆炸气泡的运动特性,在结合SOLA算法差分格式,利用投影算法对气泡的控制方程进行数值求解的基础上,基于VOF法进行了数值模拟,并利用MSC.DYTRAN软件进行仿真分析,通过对两者的分析结果对比研究.结果表明,基于结合SOLA算法差分格式的投影算法,利用VOF法可对水下爆炸气泡各时刻的运动界面进行较精确的捕捉和重构,从而实现精确追踪水下爆炸气泡脉动全物理过程的运动界面.同时,也研究了距离参数对气泡坍塌、射流及环状气泡的影响,掌握了距离参数对近水面水下爆炸气泡运动特性的影响.
The difference format of SOLA and projection algorithm were introduced to solve the control equation of bubble motion in underwater explosion.The dynamic behavior of bubble near free surface was simulated and analyzed with VOF and MSC.DYTRAN software. The computational results show that the VOF method can accurately track and reconstruct the motion of bubbles based on projection algorithm combined with the SOLA difference format.The interactions between bubble and free surface were studied systematically and summarized the dynamic behaviors of bubbles,including the influence of distance parameter on bubble collapse,eject and annularity bubble,that have close relation with distance parameter.
引文
[1]李健,荣吉利,项大林.近自由面水下爆炸气泡运动的数值计算研究[J].工程力学,2011,27(6):200-205.Li Jian,Rong Jili,Xiang Dalin.Numerical study of bubble motion by underwater explosion near free surface[J].Engineering Mechanics,2011,27(6):200-205.(in Chinese)
[2]Li Jian,Rong Jili.Bubble and free surface dynamics in shallow underwater explosion[J].Ocean Engineering,2011,38:1861-1868.
[3]Rajendran R.Numerical simulation of underwater explosion bulge test[J].Materials and Design,2009,30:4335-4341.
[4]刘雨,黄风雷,吕中杰.蒸汽爆炸过程中压力波冲击气泡现象的数值模拟[J].北京理工大学学报,2013,33(4):339-342.Liu Yu,Huang Fenglei,LüZhongjie.Numerical simulation of shock wave impacting bubble in vapor explosion[J].Transactions of Beijing Institute of Technology,2013,33(4):339-342.(in Chinese)
[5]Zhang A M,Yao X L,Yu X B.The dynamics of threedimensional underwater explosion bubble[J].Journal of Sound and Vibration,2008,311:1196-1212.
[6]Wang C,Khoo B C,Yeo K S.Elastic mesh technique for 3DBIM simulation with an application to underwater explosion bubble dynamics[J].Computers&Fluids,Volume,2003,32(9):1195-1212.
[7]Hirt C W,Nichols B D.Volume of fluid method for the dynamics of free boundary[J].Comput Phys,1981,93:201-225.
[8]赵大勇,李维仲.VOF方法中几种界面重构技术的比较[J].热科学与技术,2003,2(4):318-323.Zhao Dayong,Li Weizhong.Comprison of three kinds of interface reconstruction methods applied to VOF[J].Journal of Thermal Science and Technology,2003,2(4):318-323.(in Chinese)
[9]Ashgriz N,Poo J Y.FLAIR:flux line-segment model for advection and interface reconstruction[J].J Comput Phys,1991,93(2):449-468.
[10]Seifollahi M,Shirani E,Ashgriz N.An improved method for calculation of interface pressure force in PLIC-VOF methods[J].Eruopean Journal Mechanics B/Fluids,2008,27:1-23.
[11]万德成.用VOF法数值模拟孤立波迎撞及在后台阶上的演化[J].计算物理,1998,15(6):655-666.Wan Decheng.Numerical simulations of the collisions between two solitary waves and the transformation of a solitary wave over backward step by VOF method[J].Chinese Journal of Computational Physics,1998,15(6):655-666.(in Chinese)
[12]邹志利,邱大洪,王永学.VOF方法模拟波浪槽中二维非线性波[J].水动力学研究与进展:A辑,1996(1):93-103.Zou Zhili,Qiu Dahong,Wang Yongxue.Numerical simulation of nonlinear wave generated in wave flume by VOF technique[J].Journal of Hydrodynamics:Serial A,1996(1):93-103.(in Chinese)
[13]Eugen F F Botta,Marcel H M Ellenbroek.A modified SOR method for the poisson equation in unsteady freesurface flow calculations[J].Journal of Computational Physics,1985,60(1):119-134.