HIGH-ORDER WENO SIMULATIONS OF THREE-DIMENSIONAL RESHOCKED RICHTMYER-MESHKOV INSTABILITY TO LATE TIMES:DYNAMICS,DEPENDENCE ON INITIAL CONDITIONS,AND COMPARISONS TO EXPERIMENTAL DATA

The dynamics of the reshocked multi-mode Richtmyer-Meshkov instability is investigated using 513 × 257^2 three-dimensional ninth-order weighted essentially nonoscil- latory shock-capturing simulations. A two-mode initial perturbation with superposed ran- dom noise is used to model the Mach 1.5 air/SF6 Vetter-Sturtevant shock tube experiment. The mass fraction and enstrophy isosurfaces, and density cross-sections are utilized to show the detailed flow structure before, during, and after reshock. It is shown that the mixing layer growth agrees well with the experimentally measured growth rate before and after reshock. The post-reshock growth rate is also in good agreement with the prediction of the Mikaelian model. A parametric study of the sensitivity of the layer growth to the choice of amplitudes of the short and long wavelength initial interfacial perturbation is also pre- sented. Finally, the amplification effects of reshock are quantified using the evolution of the turbulent kinetic energy and turbulent enstrophy spectra, as well as the evolution of the baroclinic enstrophy production, buoyancy production, and shear production terms in the enstrophy and turbulent kinetic transport equations.

Detonation initiation developing from the Richtmyer-Meshkov instability

Detonation initiation resulting from the Richtmyer-Meshkov instability is investigated numerically in the configuration of the shock/spark-induced-deflagration interaction in a combustive gas mixture. Two-dimensional multi-species Navier-Stokes equations implemented with the detailed chemical reaction model are solved with the dispersion-controlled dissipative scheme. Numerical results show that the spark can create a blast wave and ignite deflagrations. Then, the deflagration waves are enhanced due to the Richtmyer-Meshkov instability, which provides detonation initiations with local environment conditions. By examining the deflagration fronts, two kinds of the initiation mechanisms are identified. One is referred to as the deflagration front acceleration with the help of the weak shock wave, occurring on the convex surfaces, and the other is the hot spot explosion deriving from the deflagration front focusing, occurring on the concave surfaces.

Numerical investigation of Richtmyer-Meshkov instability driven by cylindrical shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov（RM） instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

The Richtmyer-Meshkov instability of a double-layer interface in convergent geometry with magnetohydrodynamics

The interaction between a converging cylindrical shock and double density interfaces in the presence of a saddle magnetic field is numerically investigated within the framework of ideal magnetohydrodynamics.Three fluids of differing densities are initially separated by the two perturbed cylindrical interfaces.The initial incident converging shock is generated from a Riemann problem upstream of the first interface.The effect of the magnetic field on the instabilities is studied through varying the field strength.It shows that the Richtmyer-Meshkov and Rayleigh-Taylor instabilities are mitigated by the field,however,the extent of the suppression varies on the interface which leads to non-axisymmetric growth of the perturbations.The degree of asymmetry of the interfacial growth rate is increased when the seed field strength is increased.

Multiple-Relaxation-Time Lattice Boltzmann Approach to Richtmyer-Meshkov Instability

The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time （MRT） Lattice Boltzmann （LB） model proposed in [Europhys. Lett. 90 （2010） 54003]. We discuss the reason why the Gram Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka-Tsutahara model [Phys. Rev. E 69 （2004） 035701（R）] is only valid in subsonic flows. The yon Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer- Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang-Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.

The Growth of Richtmyer-Meshkov Instability under Multiple Impingements

The growth of multi-mode Richtmyer-Meshkov instability under multiple impingements and the effect of initial shock strength on the growth of RMI are numerically investigated. We obtain the time evolution of turbulent mixing zone width for initial shock with different strength. The results show that the turbulent mixing zone width grows in a different manner at different stage but strictly in a similar way for the initial shock with different strength. Also, the initial shock strength has a significant effect on the growth rate of turbulent mixing zone width, especially before reshock, but can not change the growth laws in the whole process.

激波冲击下air/SF_6界面的Richtmyer-Meshkov不稳定性

实验研究了低马赫数(1.27)激波作用air/SF6界面的RM不稳定性问题。air/SF6初始正弦界面由厚度为1～2!m的薄膜相隔得到,用阴影法测试界面演化过程。实验结果表明:由于不稳定性重流体(SF6)向轻流体(air)演化成"尖钉"结构,而轻流体演化为"气泡"结构;由于界面切向速度差的Kelvin-Helm-holtz不稳定性,"尖钉"头部翻转成蘑菇头形状;之后,蘑菇杆破碎。扰动振幅发展的实验结果与Zhang-Sohn模型和PPM数值计算的结果较吻合。

小扰动界面形态对RM不稳定性影响的数值分析

采用Navier–Stokes方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer–Meshkov不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性.研究结果发现:单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后.扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利.随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.

液柱在激波冲击下RM不稳定性和破裂过程的数值计算

对液柱受到激波冲击后在气流中的变形断裂过程进行了数值计算,研究了在该过程中RichtmyerMeshkov(RM)不稳定性的具体表现。应用Fluent软件,数值模拟了二维(2D)和三维(3D)液柱在激波马赫数为1.10、液柱初始直径为2.76mm的情况下,气/液界面上RM不稳定性的演化过程以及液体周围的流场。计算结果表明,初始扰动数目对RM不稳定性的影响显著;液柱在横截面平面(Z方向)发生变形失稳,导致沿液柱轴向出现变形失稳。数值计算的结果与已有的实验结果吻合较好。

一种用于RM不稳定性研究的竖直环形激波管的设计与验证

设计并加工了一套竖直环形同轴无膜激波管,可用于环形汇聚激波诱导下的Richtmyer-Meshkov不稳定性实验研究。与前人工作相比,本文在流体界面的形成以及流场的观测方法上做了较大的改进。通过实验和数值方法,对该竖直激波管产生的环形柱状汇聚激波的参数进行测量和分析,验证了同轴激波管形成柱状汇聚激波方法的可行性和可靠性。在界面形成方面,采用细丝约束肥皂膜技术形成正八边形气体界面,并利用数值方法考察了细丝对界面发展的影响。结果表明在界面发展的前期,细丝的影响几乎可以忽略。利用连续激光片光结合高速摄影相机对流场进行观测,获得了正八边形air/SF6气体界面在环形汇聚激波及其反射激波冲击下的演化过程,并与数值结果进行了对比,获得了较好的一致性,进一步验证了汇聚激波的对称性以及细丝约束肥皂膜技术用于形成多边形气体界面的可靠性。

RM不稳定过程中预混火焰界面演化及混合区增长预测

预混火焰界面的RM(Richtmyer-Meshkov)不稳定导致的界面混合区增长过程在自然界和工程实践中十分常见,但化学反应对其增长的影响机理仍不明确,反应性界面混合区增长速率的预测也未见报道,因此,开展RM不稳定过程中火焰界面演化和混合区预测的研究十分必要.本文采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了正弦形预混火焰界面在平面入射激波及其反射激波作用下的RM不稳定过程.结果表明,在入射激波作用后的阶段,除RM不稳定本身导致的界面演化为"钉-帽"和"泡"形结构外,化学反应一方面以预混火焰传播的方式促进了界面中"泡"结构的增长,另一方面通过与涡结构的复杂相互作用促进了"钉-帽"结构的增长.化学反应活性越强,火焰界面的"泡"结构和"钉-帽"结构的增长越快.在第一次反射激波作用后的阶段,化学反应以相同的火焰传播方式对"泡"和"钉-帽"结构产生影响,两者效应相抵,因而导致反射激波作用后的阶段中界面混合区增长不受化学反应活性的影响.根据以上分析,分别针对入射激波和第一次反射激波作用后的火焰界面混合区增长速率提出了相应的预测模型,为探索反应性RM不稳定过程的理论预测方法提供了有益参考.

A three-dimensional numerical study on instability of sinusoidal flame induced by multiple shock waves

The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier-Stokes (NS) equations with a single-step chemical reaction and a high resolution, 9th-order weighted essentially non-oscillatory scheme. The computational results were validated by the grid independence test and the experimental results in the literature. The computational results show that after the passage of incident shock wave the flame interface develops in symmetric structure accompanied by large-scale transverse vortex structures. After the interactions by successive reshock waves, the flame interface is gradually destabilized and broken up, and the large-scale vortex structures are gradually transformed into small-scale vortex structures. The small-scale vortices tend to be isotropic later. The results also reveal that the evolution of the flame interface is affected by both mixing process and chemical reaction. In order to identify the relationship between the mixing and the chemical reaction, a dimensionless parameter, , that is defined as the ratio of mixing time scale to chemical reaction time scale, is introduced. It is found that at each interaction stage the effect of chemical reaction is enhanced with time. The enhanced effect of chemical reaction at the interaction stage by incident shock wave is greater than that at the interaction stages by reshock waves. The result suggests that the parameter can reasonably character the features of flame interface development induced by the multiple shock waves.

Numerical simulations of Richtmyer-Meshkov instabilities using conservative front-tracking method

This paper presents the numerical simulations of two Richtmyer-Meshkov （RM） instability experiments using the conservative front-tracking method developed in （Mao, D. Towards front-tracking based on conservation in two space dimensions II, tracking discontinuities in capturing fashion. J. Comput. Phys., 226, 1550-1588 （2007））. The numerical results are compared with those obtained in （Holmes, R. L., Grove, J. W., and Sharp, D. H. Numerical investigation of Richtmyer-Meshkov instability using front-tracking. J. Fluid Mech., 301, 51-64 （1995））. The present simulations are generally in good agreement with those obtained by Holmes et al., and also capture the nonlinear and compessive phenomenon, i.e., the self-interactions of the transmitted and reflected wave edges, which was pointed out by Holmes et al. as the cause of the deceleration of the interfaces. However, the perturbation amplitudes and the amplitude growth rates of the interfaces obtained with the present conservative front-tracking method are a bit larger than those obtained by Holmes et al.

The Richtmyer–Meshkov instability of a 'V' shaped air/helium interface subjected to a weak shock

The Richtmyer-Meshkov instability ofa ‘V＇ shaped air/helium gaseous interface subjected to a weak shock wave is experimentally studied. A soap film technique is adopted to create a ‘V＇ shaped interface with accurate initial conditions. Five kinds of ‘V＇ shaped interfaces with different vertex angles are formed to highlight the effects of initial conditions on the flow characteristics. The results show that a spike is generated after the shock impact, and grows constantly with time. As the vertex angle increases, vortices generated on the interface become less noticeable, and the spike develops less pronouncedly. The linear growth rate of interface width after compression phase is estimated by a linear model and a revised linear model, and the latter is proven to be more effective for the interface with high initial amplitudes. The linear growth rate of interface width is, for the first time in a heavy/light interface configuration, found to be a non-monotonous function of the initial perturbation amplitude-wavelength ratio.

反应性RM不稳定混合特性研究:反应活性与介质不均匀性的影响

预混火焰界面的RM(Richtmyer-Meshkov)不稳定现象在自然界和工程实践中十分常见,但目前关于反应性RM不稳定的研究主要集中于均匀介质的情况,而实际中的预混气体往往是非均匀的,因此开展非均匀介质中火焰界面演化和混合特性的研究十分必要。采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了预混火焰界面在入射激波及反射激波作用下的RM不稳定过程,考察了化学反应活性以及介质非均匀性对RM不稳定过程中火焰界面混合特性的变化规律的影响。结果表明,在入射激波作用后的阶段,在均匀介质中的火焰界面形态呈现典型的“钉-帽-泡”结构,化学反应活性越强,界面的“泡”结构和“钉-帽”结构增长越快;而在非均匀介质中,火焰界面形态则呈现“钉-钉”结构,界面在流向速度差的诱导下被更大程度地拉伸。在第一次反射激波作用后的阶段,混合区的增长速率不依赖于反应活性和均匀性,仅与流动特性有关。时间尺度的研究表明,大尺度流动是反应性RM不稳定的主导因素,其次是化学反应,最后是小尺度混合,化学反应的强化会抑制大尺度流动,非均匀性会强化大尺度流动。

Investigation of coupling effect on the evolution of Richtmyer-Meshkov instability at double heavy square bubbles

This study investigates numerically the coupling effect on the evolution of Richtmyer-Meshkov instability at double heavy square bubbles.Five scenarios are considered,each with varying initial separations S/L(where L demotes the side length of the square)ranging from 0.125 to 1.0.Squares are filled with SF6gas,and are enclosed by N2gas.The simulations of shock-induced multispecies flow are performed by solving the two-dimensional compressible Euler equations with a higher-order explicit modal discontinuous Galerkin solver.The simulations demonstrate that the flow morphology resulting from the coupling effect is highly dependent on the separation between two squares.When the separation is large,the squares experience a weaker coupling effect and evolve independently.While,as the separation reduces,the coupling effect manifests earlier in the interaction and becomes more substantial.As a result,this phenomenon greatly intensifies the motion of inner upstream/downstream vortex rings towards the symmetry axis,leading to the emergence of multiple jets such as the twisted downward,upward,and coupled jets.A thorough exploration of the coupling effect of double squares is conducted by analyzing the vorticity production.Notably,a significant quantity of vorticity is produced along the squares interface for smaller separation.Further,these coupling effects result in various interface features(upstream/downstream movement,and height/width evolution),and temporal variations of various spatially integrated fields.Finally,the analysis of the flow structure also considers the interaction between two more flow parameters,the Mach and Atwood numbers,in order to evaluate the coupling effects.

液滴在激波冲击下的破裂过程

对液滴在入射激波作用下的变形破碎过程进行了实验研究和数值模拟,得知数值模拟结果与实验结果基本吻合,以及在什么情况下两者出现分歧。结果显示,液滴在激波的作用下要经历从压缩变形、RM不稳定性变形、细小液雾剥离到全部雾化破碎等过程。结果还表明,不同液滴直径、入射激波马赫数和液滴介质等参数下的液滴变形破碎的发展趋势是一致的,而其发展速度明显则不同。其中Weber数的增加加速了液滴的破碎,而Ohnesorge数和黏性的增加则抑制了液滴的破碎。

一个用于界面不稳定性实验研究的竖式激波管

研制了一个可用于气-气界面Rayleigh-Taylor(RT)不稳定性、气-液界面Richtmyer-Meshkov(RM)不稳定性、气-液界面RT和RM耦合不稳定性实验的竖式激波管装置,配备了相应的光学电学测量系统,给出了装置的主要技术参数,以及典型实验的压力-时间曲线。利用这套实验装置,可以精确地控制不稳定性的初始扰动界面,观测记录界面演化发展全过程。

黏性应力对液体扩张抛撒影响的实验研究

利用组合式激波管和高速摄影机,获得不同工况的气液界面不稳定性发展直至发生首次破碎的高速摄影照片。实验结果表明,RT不稳定性和RM不稳定性均随时间增长。在扰动波发展的过程中,两相混合区中气泡的前缘基本保持圆形,而尖钉的头部在发展过程中有雾化的迹象。由于黏性的存在,不稳定性的发展将受到抑制,且随着黏性系数的增大,扰动波的增长系数随之减小。

Numerical study on Rayleigh-Taylor effect on cylindrically converging Richtmyer-Meshkov instability

Evolution of a two-dimensional air/SF6 single-mode interface is numerically investigated by an upwind CE/SE method under a cylindrically converging circumstance. The Rayleigh-Taylor effect caused by the flow deceleration on the phase inversion(RTPI)is highlighted. The RTPI was firstly observed in our previous experiment, but the related mechanism remains unclear. By isolating the three-dimensional effect, it is found here that the initial amplitude(a0), the azimuthal mode number(k0) and the re-shocking moment are the three major parameters which determine the RTPI occurrence. In the variable space of(k0, a0), a critical a0 for the RTPI occurrence is solved for each k0, and there exists a threshold value of k0 below which the RTPI will not occur no matter what a0 is. There exists a special k0 corresponding to the largest critical a0, and the reduction rule of critical a0 with k0 can be well described by an exponential decay function. The results show that the occurrence of the RTPI requires a small a0 which should be less than a critical value, a large k0 which should exceed a threshold, and a right impinging moment of the re-shock which should be later than the RTPI occurrence. Finally, the effects of the incident shock strength, the density ratio and the initial position of the interface on the threshold value of k0 and on the maximum critical a0 are examined. These new findings would facilitate the understanding of the converging Richtmyer-Meshkov instability and would be helpful for designing an optimal structure of the inertia confinement fusion capsule.