低速下水滴撞击固体表面运动模拟

徐爽; 赵宁; 王春武; 吴杰

低速下水滴撞击固体表面运动模拟

1.
南京航空航天大学航空宇航学院,南京 210016
2.
南京航空航天大学理学院,南京 210016

基金项目: 国家自然科学基金(10972106)

计量
- 文章访问数: 1601
- HTML浏览量: 111
- PDF量: 1044
- 被引次数: 0
出版历程
- 收稿日期: 2012-03-22
- 刊出日期: 2013-03-28

Simulation of droplets impacting on solid surface at low velocity

1.
College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
2.
College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

摘要

摘要: 以水滴撞击飞机部件造成结冰问题为背景,运用level set界面追踪方法,把不可压黏性方程组中的表面张力用level set函数表示,并利用2阶变密度投影算法求解二维不可压Navier-Stokes方程组,数值模拟了低速下的水气界面追踪问题.为了能够处理密度比为1000的水气界面问题,在界面处引入了一定宽度的光滑过渡带宽,进而消除了由于界面两边大密度比引起的数值震荡,由此较为成功地模拟了大水滴在不同表面张力时撞击平板后的形态变化及其运动轨迹.数值实验结果表明本算法可以有效地处理界面处密度比和黏性很大的不可压多相流问题.考虑界面处的表面张力,本算法可模拟不同表面张力下的界面运动变化.
- 结冰 /
- level set /
- 投影算法 /
- 大密度比 /
- 水滴
Abstract: In the background of icing problem caused by the droplets which impacted on the surface of aircraft,the level set interface tracking method was used, the surface tension in the incompressible N-S equations was denoted by the level set function and the second-order projection method for variable-density flows was used to solve 2-D incompressible Navier-Stokes equations;the water and gas multiphase problems at low velocity were simulated. In order to simulate the water and gas multiphase problems of the density ratio up to 1000, the numerical oscillations caused by the large density ratio between the interface can be eliminated by smoothing the density at the interface within a few grids. So the topological changes and movements of large droplets,which are different in surface tension when impacting on the solid interface,can be successfully simulated. The numerical results show that the method can be used to model the incompressible multi-phase flow with large density and viscosity density.Considering the surface tension, the movements of the surface at different surface tensions can be simulated.
- icing /
- level set /
- projection method /
- large density ratio /
- droplet

HTML

参考文献(15)

[1]	朱春玲,李宇钦,张泉.基于旋转多圆柱的冰风洞液滴参数分析方法[J].航空动力学报,2007,22(2):180-186. ZHU Chunling,LI Yuqin,ZHANG Quan.Analysis method of droplets parameters based on rotating multicylinder in an icing wind tunnel[J].Journal of Aerospace Power,2007,22(2):180-186.(in Chinese)
[2]	Mark S,Ann S A.An adaptive level set approach for incompressible two-phase flows[J].Journal of Computational Physics,1999,148(1):81-124.
[3]	Mark S,Peter S,Osher S.A level set approach for computing solutions to incompressible two phase flow[J].Journal of Computational Physics,1994,114(1):146-159.
[4]	Mark S,Emad F,Peter S,et al.A level set approach for computing solutions to incompressible two phase flow Ⅱ[R].Lawrence Livermore National Labreport,UCRL-JC-121224,1995.
[5]	Fedkiw R P,Tariq A,Barry M,et al.A non-oscillatory Eulerian approach to interface in multimaterial flows (the ghost fluid method)[J].Journal of Computational Physics,1999,152(2):457-492.
[6]	John B,Phillip C.A second-order projection method for the incompressible Navier-Stokes equation[J].Journal of Computational Physics,1985,85(2):257-283.
[7]	Joho B,Phillip C,Louis H.An efficient second-order projection method for vicous incompressible flow[C]//AIAA 10th Computational Fluid Dynamics Conference.Honduu Hawaii:AIAA,1991:360-367.
[8]	Joho B,Daniel L.A second-order projection method for variable-density flow[J].Journal of Computational Physics,1992,101(2):334-348.
[9]	Brackbill J U,Kothe D B,Zemach C.A continuum method for modeling surface tension[J].Journal of Computational Physics,1992,100(2):335-354.
[10]	Peskin C S.Numerical analysis of blood flow in the heart[J].Journal of Computational Physics,1977,25(3):220-252.
[11]	Shu C W,Osher S.Efficient implementation of essentially non-oscillatory shock-capturing:schemes Ⅱ[J].Journal of Computational Physics,1989,83(1)32-78.
[12]	Tatebe O.The multigrid preconditioned conjugate gradient method.[EB/OL] .[2011-9-19].http://www.hpcs.cs.tsukuba.ac.jp/~tatebe/research/paper/SWoPP92.pdf.
[13]	Shu C W.Total-variation-diminishing time discretization[J].SIAM Journal on Scientific and Statical Compating,1988,9(6):1073-1084.
[14]	Mulder W,Osher S,James.A computing interface motion in compressible gas dynamics[J].Journal of Computational Physics,1992,100(2):209-228.
[15]	Mark S,Perer S.Axisymmetric free boundary problems[J].Journal of Fluid Mechanics,1997,341(1):269-294.