偏心对汇聚激波诱导的RM不稳定性影响的数值研究

计算物理 ›› 2016, Vol. 33 ›› Issue (1): 66-74.

偏心对汇聚激波诱导的RM不稳定性影响的数值研究

何惠琴, 翟志刚, 司廷, 罗喜胜

中国科学技术大学近代力学系先进推进实验室, 合肥 230027

收稿日期:2015-03-10 修回日期:2015-03-11 出版日期:2016-01-25 发布日期:2016-01-25
作者简介:何惠琴(1991-),女,安徽池州,硕士研究生,从事计算流体力学,RM不稳定性的研究,E-mail:hhq@mail.ustc.edu.cn
基金资助:
国家自然科学基金（11302219，11272308）资助项目

Numerical Investigation on Eccentricity Effect in Richtmyer-Meshkov Instability Induced by Converging Shock Wave

HE Huiqin, ZHAI Zhigang, SI Ting, LUO Xisheng

Advanced Propulsion Laboratory, Department of Modern Mechanics, University of Science and Technology of China, Hefei 230027, China

Received:2015-03-10 Revised:2015-03-11 Online:2016-01-25 Published:2016-01-25

摘要/Abstract

摘要： 数值研究汇聚激波与四种形状（圆形、小振幅单模、大振幅单模和正方形）的二维气柱界面相互作用，激波汇聚中心与界面同心和不同心（即偏心）时Richtmyer-Meshkov（RM）不稳定性的发展规律，重点考察界面中心的压力及混合区面积在两种情况下随时间的变化.数值方法使用VAS2D程序，该方法采用有限体积法结合网格自适应技术，能够达到时间和空间的二阶精度.结果表明，偏心情况下RM不稳定性是其在同心情况下的扰动和偏心小扰动叠加的结果.在本文采用的偏心程度下（20%），偏心对于圆形无扰动界面发展的影响主要表现在后期界面出现微小扰动结构；而对于单模和正方形这种原本有扰动的界面，偏心使扰动结构呈现不对称及扭曲，同时也影响了界面中心压力和混合区面积，因而加剧了不稳定性的发展.

关键词: Richtmyer-Meshkov不稳定性, 气柱, 偏心, 涡量

Abstract: Development of Richtmyer-Meshkov instability induced by a converging shock wave and its reflected shock from converging center is numerically investigated in which emphasis is on eccentricity effect. VAS2D uses finite volume method which has a secondorder precision in both temporal and spatial scales. Four types of 2-dimensional SF₆ column (circle, small-amplitude single-mode, large-amplitude single mode and square) surrounded by air are adopted. Interface morphologies show that development of RM instability in eccentric case results from perturbation growth in concentric case combined with little perturbation growth caused by eccentricity. For interface without initial perturbation, eccentricity results in appearance of tiny perturbation structures. For interfaces with initial perturbation, eccentricity makes perturbation structures asymmetric and distorted. Meanwhile, eccentricity affects pressure of interface center and area of mixed zone, which are closely related to development of RM instability.

Key words: Richtmyer-Meshkov instability, gas column, eccentricity, vorticity

中图分类号:

O534

何惠琴, 翟志刚, 司廷, 罗喜胜. 偏心对汇聚激波诱导的RM不稳定性影响的数值研究[J]. 计算物理, 2016, 33(1): 66-74.

HE Huiqin, ZHAI Zhigang, SI Ting, LUO Xisheng. Numerical Investigation on Eccentricity Effect in Richtmyer-Meshkov Instability Induced by Converging Shock Wave[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(1): 66-74.

参考文献

[1] RICHTMYER R D. Taylor instability in shock acceleration of compressible fluids[J]. Commun Pure Appl Math, 1960, 13: 297-319.
[2] MESHKOV E E. Instability of the interface of two gases accelerated by a shock wave[J]. Fluid Dyn, 1969,(4): 101-104.
[3] JACOBS J W. The dynamics of shock accelerated light and heavy gas cylinders[J]. Phys Fluids A, 1993,(5): 2239-2247.
[4] RANJAN D, OAKLEY J, BONAZZA R. Shock-bubble interactions[J]. Annu Rev Fluid Mech, 2011, 43: 117-140.
[5] ZHAI Z, SI T, LUO X, et al. On the evolution of spherical gas interface accelerated by a planar shock wave[J]. Phys Fluids, 2011, 23: 084104.
[6] HAAS J F, STUTEVANT B. Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities[J]. J Fluid Mech, 1987, 181: 41-76.
[7] 范美如, 翟志刚, 司廷, 等. 激波与不同形状界面作用的数值模拟[J]. 中国科学, 2011, 41: 862-869.
[8] SHANKAR S K, KAWAI S, LELE S K. Two-dimensional viscous flow simulation of a shock accelerated heavy gas cylinder[J]. Phys Fluids, 2011, 23: 024102.
[9] BAI J S, ZOU L Y, WANG T, et al. Experimental and numerical study of shock-accelerated elliptic heavy gas cylinders[J]. Phys Rev E, 2010, 82: 056318.
[10] BROUILLETTE M. The Richtmyer-Meshkov instability[J]. Annu Rev Fluid Mech, 2002, 34: 445-468.
[11] PERRY R W, KANTROWITZ A. The production and stability of converging shock waves[J]. J Appl Phys, 1951, 22(7): 878-886.
[12] ZHAI Z, SI T, LUO X, et al. Parametric study of cylindrical converging shock waves generated based on shock dynamics theory[J]. Phys Fluids, 2012, 24: 026101.
[13] DIMOTAKIS P E, SAMTANEY R. Planar shock cylindrical focusing by a perfect-gas lens[J]. Phys Fluids, 2006, 18(3), 031705.
[14] HOSSEINI SHR, ONODERA O, TAKAYAMA K. Characteristics of an annular vertical diaphragmless shock tube[J]. Shock Waves, 2000, 10(3): 151-158.
[15] HOSSEINI S H R, TAKAYAMA K. Experimental study of Richtmyer-Meshkov instability induced by cylindrical shock waves[J]. Phys Fluids, 2005, 17: 084101.
[16] SI T, ZHAI Z, LUO X. Experimental study of Richtmyer-Meshkov instability in a cylindrical converging shock tube[J]. Laser Part Beams, 2014, 14: 263-346.
[17] SOD G A. A numerical study of a converging cylindrical shock[J]. J. Fluid Mech, 1977, 83(04): 785-794.
[18] ZHANG Q, GRAHAM M J. A numerical study of Richtmyer-Meshkov instability driven by cylindrical shocks[J]. Phys Fluids, 1998,(10): 974-992.
[19] ZHENG J G, LEE T S, WINOTO S H. Numerical simulation of Richtmyer-Meshkov instability driven by imploding shocks[J]. Math Comput Simulat, 2008, 79(3): 749-762.
[20] SUN M, TAKAYAMA K. Conservative smoothing on an adaptive quadrilateral grid[J]. J Comput Phys, 1999, 150(1): 143-180.
[21] 杨伟航, 罗喜胜. 柱状汇聚激波冲击重气柱不稳定性的数值研究[J]. 中国科学技术大学学报, 2014,44(6): 488-495.