Numerical simulation of shock impact on material interface
-
摘要: 计算激波与物质界面相互作用问题的关键是精确捕获运动物质界面.采用level set捕获物质界面,用改进虚拟流(modified ghost fluid method)方法定义虚拟流节点参数,从而将多介质流问题分解成多个单一介质流问题求解.应用该方法对激波与物质界面相互作用系列问题进行了数值模拟,并给出了物质界面的演化过程.计算结果发现激波强度、物质界面对涡的形式有很大影响,界面两边的密度比决定了界面是否反转.
-
关键词:
- euler方程 /
- level set方法 /
- 改进虚拟流 /
- 激波 /
- richtmyer-meshkov不稳定性
Abstract: The key of numerical simulation of the problems of shock impact on interface is capturing accurately the material interfaces. The position of material interfaces were captured by level set, the parameter of ghost points were defined by modified ghost fluid method, accordingly the problem of multi-interfaces was split into several the problems of signal material interface. Based on the methods mentioned above, the numerical simulation of the series problems of shock impact on material interface were presented, and the evolvement of material interfaces were given. It is found that the strength of shock wave and the shape of material interface have great effect to vortex, whether or not the interface reverse depends on the density of two materials.-
Key words:
- euler equations /
- level set method /
- modified ghost fluid method /
- shock /
- richtmyer-meshkov instability
-
[1] Hawley J F, Zabusky N J.Vortex paradigm for shock accelerated density satisfied interfaces[J].Physical Review Letter,1989,63(12):1241-1244 [2] 严长林,孙得军,尹协远,等.激波加速的密度分层流体界面的演化[J].空气动力学学报,2001,19(2):228-233 Yan Changlin, Sun Dejun, Yin Xieyuan, et al.Evolution of shock-accelerated density-stratified interfaces [J].Acta Aerodynamica Sinica, 2001, 19(2):228-233(in Chinese) [3] Sussman M,Osher P S. A level set approach for computing solution to incompressible two-phase flow[J]. Journal of Computational Physics, 1994, 114:146-159 [4] Liu T G ,Khoo B C ,Yeo K S.Ghost fluid method for strong shock impacting on material interface[J].Journal of Computational Physics, 2003, 190:651-681 [5] Falcovitz J, Alfandary G, Hanoch G. A two-dimensional conservation laws schemes for compressible flows with moving boundaries [J].Journal of Computational Physics, 1997, 138:83-102 [6] Toro E F, Spruce M, Speares W. Restoration of the contact surface in the HLL-Riemann solver [J].Shock Waves, 1994, 4:25-34 [7] 刘儒勋,舒其望.计算流体力学的若干新方法[M].北京:科学出版社,2003 Liu Ruxun, Shu Qiwang.Several new methods of computational fluid dynamics [M].Beijing: Science Press, 2003 (in Chinese) [8] 王晨星,唐维军,程军波.应用多介质PPM方法计算斜激波与物质交界面的相互作用[J].计算物理,2004,21:532-527 Wang Chenxing,Tang Weijun,Cheng Junbo.Numerical computations of the refraction of a shock wave at interface by multi-component PPM algorithm[J].Comput Phys,2004,21:532-527 (in Chinese)
点击查看大图
计量
- 文章访问数: 3152
- HTML全文浏览量: 105
- PDF下载量: 1113
- 被引次数: 0