On the Structure of the Mixing Zone at an Unstable Contact Boundary

The interface between two media of different densities (contact boundary) moving with an acceleration directed from the less dense medium to the more dense one is unstable (Rayleigh–Taylor instability) [1, 2]. The initial perturbations of the interface grow indefinitely and, as a result, a medium mixing zone growing with time is formed at the interface. The structure of such a mixing zone at gas–gas and gas–liquid interfaces is discussed on the basis of laboratory experiments on shock tubes of various types. It is concluded that the regions of turbulent and laminar flows are combined in the mixing zone.     

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

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

Understanding and modeling turbulent fluxes and entrainment in a gravity current

We present an experimental study of the mixing processes in a gravity current flowing on an inclined plane. The turbulent transport of momentum and density can be described in a very direct and compact form by a Prandtl mixing length model: the turbulent vertical fluxes of momentum and density are found to scale quadratically with the vertical mean gradients of velocity and density. The scaling coefficient, the square of the mixing length, is approximately constant over the mixing zone of the stratified shear layer. We show how, in different flow configurations, this length can be related to the shear length of the flow (ε/∂_zu³)^1/2. We also study the fluctuations of the momentum and density turbulent fluxes, showing how they relate to mixing and to the entrainment/detrainment balance. We suggest a quantitative measure of local entrainment and detrainment derived from observed conditional correlations of density flux and density or vertical velocity fluctuations.     

高速摄影在流体动力学不稳定性研究中的应用

用高速摄影技术研究了高压气体膨胀驱动空气-水界面的瑞利-泰勒不稳定性,获得了空气-水界面的不稳定性清晰图像,得到气炮尖顶运动速度及湍流混合层高度增长速度与时间关系曲线。在横式激波管上用高速纹影诊断技术研究了激波作用空气-SF6界面的里克特迈耶-梅什科夫不稳定性,初步获得了实验图像,可清晰显示混合区变化过程。     

Experimental study of initial condition dependence on turbulent mixing in shock-accelerated Richtmyer–Meshkov fluid layers

Experimental evidence is needed to verify the hypothesis that the memory of initial conditions is retained at late times in variable density flows. If true, this presents an opportunity to “design” and “control” late-time turbulence, with an improved understanding in the prediction of inertial confinement fusion and other general fluid mixing processes. In this communication, an experimental and theoretical study on the effects of initial condition parameters, namely, the amplitude δ₀ and wavenumber κ₀ , where λ₀ is the initial wavelength) of perturbations, on late-time turbulence and mixing in shock-driven Richtmyer–Meshkov (R-M) unstable fluid layers in a 2D plane is presented. Single and multi-mode membrane-free initial conditions in the form of a gas curtain having a light-heavy-light configuration (air-SF₆-air) with an Atwood number of A= 0.57 were used in our experiments. A planar shock wave with a shock Mach number M= 1.21 drives the R-M instability, and the evolution of this instability after incident shock is captured using high resolution simultaneous planar laser induced fluorescence (PLIF) and particle image velocimetry (PIV) diagnostics. Time evolution of statistics such as amplitude of the mixing layer, 2D turbulent kinetic energy, Reynolds number, rms of velocity fluctuations, probability density functions, and density-specific volume correlation were observed to quantify the amount of mixing and understand the nature of turbulence in this flow. Based on these results, it was found that the R-M mixing layer is asymmetric and non-Boussinesq. There is a correlation between initial condition parameters and large-scale, and small-scale mixing at late times, indicating an initial condition dependence on R-M mixing.     

Lock-exchange flows in inclined pipes: the relevance of the Prandtl mixing length model

We examine the applicability of the Prandtl mixing length model to transverse momentum and mass flux in strongly confined, stably stratified turbulent shear flows. These fluxes were measured in the vertical diametral plane of lock-exchange flows in an inclined pipe by the simultaneous use of planar laser-induced fluorescence and particle image velocimetry at local Reynolds numbers ranging from Re = 580 to 1770 and Richardson numbers ranging from Ri = 0.26 and 1.6. Measurements indicate that the eddy diffusivities of mass and momentum are symmetric about the pipe axis, with their maximum at the axis. The corresponding Prandtl mixing lengths decrease with increasing distance from the pipe axis within the central 60% of the pipe cross-section. Within the range of experimental conditions, the mixing lengths at the axis increase linearly with Ri so that the corresponding turbulent Prandtl number Pr_t decreases with Ri. In contrast, Pr_t and the mixing lengths do not display a systematic dependence on Re. Comparison with unbounded and semi-bound shear flows suggests that the strong confinement imposed by the pipe wall may be constraining the integral length scale and Prandtl mixing lengths.     

局域体全息光栅衍射的平面波与柱面波的波前转换研究

利用一维耦合波理论研究了入射方向相互垂直的一束平面波与一束柱面波干涉形成光折变局域体全息光栅的波前转换情况。给出了透射波和衍射波振幅的解析表达式。讨论了该局域体全息光栅的几何尺寸以及全息透镜的焦距对光栅衍射效率的影响。结果表明,随着光栅横向尺寸的增加,光栅的衍射效率也逐渐增加。然而,光栅的衍射效率却随着光栅纵向尺寸的增加而减小。全息透镜的焦距越长,光栅的衍射效率也越大。此外,分析了光栅的布拉格选择特性,该光栅具有非常好的角度选择特性。在平面波与柱面波干涉形成全息透镜的设计中,为了获得最优的衍射效率,应当根据要求合理地选择记录所用柱面波和设计光栅的几何尺寸。     

Experimental study of the stability of the interface between a liquid electrolyte and the glow discharge plasma

The stability of the interface between a liquid electrolyte and the plasma of a contracted low-pressure dc glow discharge in air is investigated by means of digital photography. Water solutions of potassium permanganate and copper sulfate were used as electrolytes. It is found that, in the case of potassium permanganate, the instability of the interface leads to ejection of the electrolyte into the plasma and extinction of the discharge. Discharge modes with different types of quasi-steady interface are observed for copper sulfate at different values of the discharge current: a smooth interface, a solitary wave perturbation, regular ripples, and a churning foamed turbulent mixing zone.     

激波冲击V形界面重气体导致的壁面与旋涡作用及其对湍流混合的影响

基于多组分混合物质量分数模型,采用色散最小耗散可控的高分辨率有限体积方法,数值模拟了弱激波冲击V形空气/SF_6界面后,界面不稳定性生成的旋涡与固体壁面作用问题.激波冲击V形界面之后,因斜压效应诱导涡量沉积在界面附近,形成沿界面规则排列的多个涡对结构.旋涡的诱导作用使界面不断变形和卷起,同时旋涡之间不断发生相互并对,诱导更多更小尺度的旋涡产生.旋涡诱导作用的叠加效应,使界面尖端处的初始涡对向上下壁面发展.随后,涡结构开始与壁面发生复杂的相互作用.旋涡与壁面作用后沿壁面加速,使得物质界面沿壁面伸展,随后,旋涡从壁面回弹,并诱导二次旋涡产生.旋涡与壁面相互作用的过程,能够明显加剧物质混合.本文从物质混合的角度研究了该过程的机理,分析了旋涡与壁面作用对物质混合的影响.     

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

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

Numerical simulation of the hydrodynamic instability experiments and flow mixing

Based on the numerical methods of volume of fluid (VOF) and piecewise parabolic method (PPM) and parallel circumstance of Message Passing Interface (MPI), a parallel multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and performed to study the hydrodynamic instability and flow mixing. Firstly, the MVPPM code is verified and validated by simulating three instability cases: The first one is a Riemann problem of viscous flow on the shock tube; the second one is the hydrodynamic instability and mixing of gaseous flows under re-shocks; the third one is a half height experiment of interfacial instability, which is conducted on the AWE’s shock tube. By comparing the numerical results with experimental data, good agreement is achieved. Then the MVPPM code is applied to simulate the two cases of the interfacial instabilities of jelly models accelerated by explosion products of a gaseous explosive mixture (GEM), which are adopted in our experiments. The first is implosive dynamic interfacial instability of cylindrical symmetry and mixing. The evolving process of inner and outer interfaces, and the late distribution of mixing mass caused by Rayleigh-Taylor (RT) instability in the center of different radius are given. The second is jelly layer experiment which is initialized with one periodic perturbation with different amplitude and wave length. It reveals the complex processes of evolution of interface, and presents the displacement of front face of jelly layer, bubble head and top of spike relative to initial equilibrium position vs. time. The numerical results are in excellent agreement with that experimental images, and show that the amplitude of initial perturbations affects the evolvement of fluid mixing zone (FMZ) growth rate extremely, especially at late times.     

爆轰波对碰驱动平面锡飞层的动力学及动载行为特征研究

强冲击下的物质变形、破坏及诱发的轻重介质混合问题,是内爆压缩科学和工程应用领域的研究重点.本文针对爆轰波对碰条件下的复杂加载动力学过程及其动载破坏形态特征,开展数值模拟研究与极曲线理论分析.设计了爆轰波对碰驱动平面锡飞层的计算模型,获得了爆轰加载动力学过程及波系相互作用物理图像,分析了锡飞层对碰区自由表面速度历史的典型特征.给出了锡飞层中折射激波对碰发生马赫反射的临界条件,解读了三波结构的传播行为,阐明了对碰区内存在"一维正冲击"区域,一维区外存在单次斜冲击向两次斜冲击过渡的复杂加载动力学过程,提出了对碰区冲击动力学模型,揭示了影响对碰区动载行为特征的机理.数值模拟结果与极曲线理论分析结果相互印证,符合较好.本文的研究成果,将为深入理解和解读对碰区特殊的物质破坏及混合现象提供重要的理论支撑.     

Optical study of parameters of the passage of a shock wave from air to water

A study is made of the penetration of shock waves from air into water. The shock wave in air is generated as a result of dielectric breakdown induced by pulsed CO₂-laser radiation. A combination of the double-exposure shadow method and holographic interferometry is used to measure the shock-wave parameters. Density and pressure profiles behind the wave front are obtained at different times after onset of breakdown. It is shown experimentally that as the wave passes through the interface from the air to the water, there is a fourfold amplification of the pressure in the shock wave front. Estimates of the width of the shock wave front formed in the water are given in the context of studies of large-scale explosion processes. It is shown that simple empirical dependences, established in the course of studies of large-scale explosions, are also valid with certain corrections for microscopic laboratory experiments. Zh. Tekh. Fiz. 68, 39–43 (August 1998)     

Numerical simulations of shock-wave interaction with a boundary layer in the plane supersonic flows with jet injection

A supersonic air flow in a plane channel with a transverse turbulent jet of hydrogen injected through a slot on the bottom wall is simulated. The algorithm for solving the Favre-averaged Navier-Stokes equations for the flow of a perfect multispecies gas on the basis of the WENO scheme is proposed. The main attention is paid to the interaction of the shock-wave structure with the boundary layers on the upper and lower duct walls under the conditions of an internal turbulent flow. Namely, a detailed study of the structure of the flow is done, and separation and mixing depending on the jet slot width are investigated. It is found that in addition to well-known shock-wave structures produced by the interaction of the free stream with the transverse jet and the bow shock interaction with the boundary layers near the walls, an additional system of shock waves and the flow separation arise on the bottom wall downstream at some distance from the jet. The comparison with the experimental data is performed.     

The impact of reduced chemistry on auto-ignition of H2 in turbulent flows

The auto-ignition behaviour of hydrogen in a turbulent flow field has been studied through a combination of detailed and systematically reduced chemistry with a transported PDF approach closed at the joint-scalar level. Radiation is accounted for through the RADCAL method and the inclusion of enthalpy into the joint-scalar PDF. Molecular mixing is closed using the modified Curl's model with the mixing frequency accounted for via two algebraic closures. The main aim of the work is to compare the impact of alternative chemical mechanisms on auto-ignition and to explore the accuracy that can be expected when reactive scalars are sequentially removed through the application of quasi-steady-state approximations (QSSAs). Two different detailed mechanisms were tested to establish the effects of intrinsic uncertainties in the detailed chemistry and to provide reference points to past work. The mechanisms feature nine solved species and 19 or 20 reversible chemical reactions. The chemical mechanisms were subsequently systematically reduced to five, four and three independent scalars through the successive introduction of QSSAs for H₂O₂, HO₂ and O. Resulting inaccuracies were quantified following each simplification step with reference to experimental data obtained in shock tubes and under turbulent flow conditions in the Cabra burner configuration. A sensitivity analysis was also performed to identify the relative impact of uncertainties in key reactions as compared to systematic simplification process. It was found that alternative recommended rates for the O + H₂ = OH + H reaction have an impact on the point of flame stabilization that is similar to that observed as a consequence of the simplification process. The work also shows that realistic results can be obtained with simplified chemistry. However, it is also concluded that the temporal evolution of the radical pool and the point of stabilization is affected by the introduction of a QSSA for the O radical. Furthermore, it is shown by comparisons with time resolved OH radical data obtained in shock tubes that the progressive elimination of species via QSSA leads to a shortening of ignition delay times and that the same effects are present, but less severe, in turbulent flow fields.     

二次曝光全息术在圆柱形冲击波会聚过程观测中的应用

本文对全息曝光技术进行了分析,设计了一种全息干涉仪,对圆柱形冲击波会聚过程的不稳定性进行了观察,并得出了几个有重要价值的结论。     

Diffraction of a two-shock configuration by a concave cylindrical surface

Diffraction of a two-shock configuration by a concave cylindrical surface with a continuously varying angle of diffraction in a perfect (inviscid and non-heat-conducting) gas is numerically investigated. Peculiarities of flow formation in the perturbed region due to a large initial angle of diffraction are revealed. The reasons for the appearance of the reverse flow over the cylindrical surface, which causes a vortex, are found. The influence of the vortex on the behavior of the loose end of the TU layer (slipstream) is studied. It is shown that the structure of the loose end of the TU layer is determined by vortex-related mixing in the gas. The nonstationary interaction of the diffracting shock wave and a sublayer shock with a concave cylindrical wall and the symmetry plane, which is accompanied by a change in the type of reflection, is examined. Published in Zhurnal Tekhnicheskoĭ Fiziki, 2007, Vol. 77, No. 10, pp. 24–33. The article was translated by the authors.     

The numerical study of shock-induced hydrodynamic instability and mixing

Based on multi-fluid volume fraction and piecewise parabolic method (PPM), a multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and applied to the problems of shock-induced hydrodynamic interfacial instability and mixing. Simulations of gas/liquid interface instability show that the influences of initial perturbations on the fluid mixing zone (FMZ) growth are significant, especially at the late stages, while grids have only a slight effect on the FMZ width, when the interface is impulsively accelerated by a shock wave passing through it. A numerical study of the hydrodynamic interfacial instability and mixing of gaseous flows impacted by re-shocks is presented. It reveals that the numerical results are in good agreement with the experimental results and the mixing growth rate strongly depends on initial conditions. Ultimately, the jelly layer experiment relevant to the instability impacted by exploding is simulated. The shape of jelly interface, position of front face of jelly layer, crest and trough of perturbation versus time are given; their simulated results are in good agreement with experimental results.     

Computations of instability and turbulent mixing by Nikiforov’s model

The results of modeling several laboratory experiments, including a large class of advanced experimental studies of turbulent flows, are presented. The results of the Meshkov’s “cylindrical” and “planar” experiments on the confluence of two zones of turbulent mixing, the experiments of Poggi, Barre, and Uberoi have been carried out using the Nikiforov’s model. The presented results attest that the Nikiforov’s model qualitatively describes the considered class of flows if the mean gas-dynamic quantities are computed with a high accuracy in the technique, and the width of the front of the finite-difference shock wave does not depend on the size of the computational grid cell.     

Flow and mixing fields of turbulent bluff-body jets and flames

The mean structure of turbulent bluff-body jets and flames is presented. Measurements of the flow and mixing fields are compared with predictions made using standard turbulence models. It is found that two vortices exist in the recirculation zone; an outer vortex close to the air coflow and an inner vortex between the outer vortex and the jet. The inner vortex is found to shift downstream with increasing jet momentum flux relative to the coflow momentum flux and gradually loses its circulation pattern. The momentum flux ratio of the jet to the coflow in isothermal flows is found to be the only scaling parameter for the flow field structure. Three mixing layers are identified in the recirculation zone. Numerical simulations using the standard k-? and Reynolds stress turbulence models underpredict the length of the recirculation zone. A simple modification to the C₁ constant in the dissipation transport equation fixes this deficiency and gives better predictions of the flow and mixing fields. The mixed-is-burnt combustion model is found to be adequate for simulating the temperature and mixing field in the recirculation zone of the bluff-body flames.