Shock-induced mixing zone characterization: an attempt by hot-wire diagnostic

An attempt to extract velocity and molar fraction from a single hot-wire trace within a turbulent mixing zone induced by a shock accelerated gaseous interface has been proposed. Experiments have been conducted for negative and positive density jumps across the interface. The hot-wire signals clearly show interfaces between mixed and unmixed regions and the locations of incident and reflected shocks. With some hypotheses on the temperature, velocity and molar fraction profiles within the turbulent mixing zone have been obtained solving an inverse problem. Results show that if the molar fraction profiles follow physically coherent evolutions, those of the local velocity are strongly correlated with the choice of its variation range. So, we reasonably think that the results obtained from single wire have to remain limited to interface and shock locations. And it is only by coupling the present technique with the laser Doppler velocimetry, which we will be able to possibly obtain reliable estimates of the variations of quantities in the turbulent mixing zone.      

An attempt to reduce the membrane effects in Richtmyer–Meshkov instability shock tube experiments

A new device has been developed to reduce the effects of the initial materialization of the gaseous interface in the context of horizontal shock tube experiments for the Richtmyer–Meshkov instability study. The thin nitrocellulosic membrane deposited on a stereolithographed grid support woven with thin wires is destroyed by thermal effect, through a powerful electrical pulse, just before the arrival of the incident shock wave. We present a first attempt realized in the light/heavy gas configuration (air/SF₆) and compared with the experiments carried out without destruction. We show that the present device allows to reduce the influence of the membrane on the instability development.      

Estimation of shock induced vorticity on irregular gaseous interfaces: a wavelet-based approach

We study the interaction of a shock with a density-stratified gaseous interface (Richtmyer–Meshkov instability) with localized jagged and irregular perturbations, with the aim of developing an analytical model of the vorticity deposition on the interface immediately after the passage of the shock. The jagged perturbations, meant to simulate machining errors on the surface of a laser fusion target, are characterized using Haar wavelets. Numerical solutions of the Euler equations show that the vortex sheet deposited on the jagged interface rolls into multiple mushroom-shaped dipolar structures which begin to merge before the interface evolves into a bubble-spike structure. The peaks in the distribution of x-integrated vorticity (vorticity integrated in the direction of the shock motion) decay in time as their bases widen, corresponding to the growth and merger of the mushrooms. However, these peaks were not seen to move significantly along the interface at early times i.e. t < 10 τ, where τ is the interface traversal time of the shock. We tested our analytical model against inviscid simulations for two test cases – a Mach 1.5 shock interacting with an interface with a density ratio of 3 and a Mach 10 shock interacting with a density ratio of 10. We find that this model captures the early time (t/τ ∼ 1) vorticity deposition (as characterized by the first and second moments of vorticity distributions) to within 5% of the numerical results. PACS 47.40.Nm; 47.20.Ma     

Overview of diagnostic methods used in shock-tube investigations of mixing induced by Richtmyer–Meshkov instability

We present an overview of the diagnostic methods used in shock-tube investigations of mixing induced by Richtmyer–Meshkov instability. The different diagnostic techniques are first briefly presented, and then reviewed in a simple single table, which lists their advantages and disadvantages, their technological characteristics and domain of validity, the physical parameters measured, the laboratory in which they were developed and an assessment of their mean cost. Received 19 November 1997 / Accepted 3 March 1998     

A Two-Time-Scale Model for Turbulent Mixing Flows Induced by Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

A two-time-scale closure model for compressible flows previously developed is extended to turbulent Rayleigh-Taylor and Richtmyer-Meshkov driven flows where mixing coexists with mean pressure gradients. Two model coefficients are calibrated with the help of Canuto-Goldman's model. For several Rayleigh-Taylor configurations, it is shown that the characteristic lengths scale as t ² while the kinetic energies and spectral transfers behave as t ² and t, respectively. The computed phenomenological coefficients of Youngs' scaling law are compared with experimental data ones. Comparisons with Youngs' three-dimensional numerical simulation (The Physics of Fluids A 3 (1991) 1312) are also performed. Finally three shock tube experiments, where the Richtmyer-Meshkov instability initiates the mixing, are simulated. The mixing thickness evolution is well reproduced while the turbulence levels seem to be overestimated with such first order models. The capability of the two-time-sale model to recover available data for different turbulent flows allows us to conclude to a more universal behavior in comparison with single-time-scale models. This revised version was published online in July 2006 with corrections to the Cover Date.     

流动肥皂膜气体界面及其在柱状汇聚激波作用下演化的实验研究

在Richtmyer-Meshkov(RM)不稳定性实验研究中,形成两种不同密度流体的初始扰动界面是前提和关健.本文提出了一种流动肥皂膜气体界面生成方法,其工作原理是由细丝构成的导流框从激波管实验段穿过,肥皂液从导流框的上端注入并在重力作用下在导流框中形成流动肥皂膜,膜的两测可以分别充入不同密度的气体从而形成稳定的气体...     

Evolution of shock waves and the primary vortex loop discharged from a square cross-sectional tube

In this paper, a numerical and experimental investigation of the evolution of a transmitting shock wave and its associated primary vortex loop, which are discharged from the open end of a square cross-sectional tube, is described. The experiments were conducted in the square tube connected to a diaphragmless shock tube and the flowfield was visualized from the axial direction with diffusive holographic interferometry. The numerical simulations were carried out by solving the three-dimensional Euler equations with a dispersion-controlled scheme. The numerical results were displayed in the form of interferograms to compare them with experimental interferograms. Good agreement between the numerical and experimental results was obtained. More detailed numerical calculations were carried out, from which the three-dimensional transition of the shock wave configuration from an initial planar to a spherical shape and the development of the primary vortex loop from a square shaped to a three-dimensional structure were clearly observed and interpreted. Received 29 January 1998 / Accepted 22 May 1998     

激波在平板上诱导的旋涡强度的测量

运动激波通过两个等攻角平板后诱导出两个同向旋转的旋涡，这两个旋涡在随当地气流向下游运动的同时，绕涡核连线中心旋转。本文通过测量涡对的转动角度速度，获得了每个旋涡的强度。实验结果表明，由此测得的旋涡强度不同用于小攻角平板起动涡公式计算了起动涡强度。     

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle flows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained. Received 14 June 1996 / Accepted 19 October 1996     

Investigation of the starting process in a Ludwieg tube

The hypersonic Ludwieg tube Braunschweig (HLB) is a valve-controlled wind tunnel that has been designed for a Mach number of Ma = 5.9 and a Reynolds number range from 2.5 × 10⁶ up to 2.5 × 10⁷ 1/m. The intermittent working principle implies an unsteady onset of flow, which leads to a delay of the time frame suitable for measurements as well as to heat loads different from steady flow conditions. This work numerically simulates the starting process. It determines whether the onset of flow leads to a significant temperature rise in the model surface which in turn impacts results gathered during measurement time. The flow field in the HLB is numerically rebuilt for two operating points including valve opening. The non-stationary flow around a hyperboloid/flare configuration in the test section is calculated for one operating point including surface heating. For laminar flow it is found that due to the short duration of the starting process no significant model heating affects results obtained during measurement time.      

Interaction of a reflected shock from a concave wall with a flame distorted by an incident shock

Observations are presented from calculations where a laminar spherical CH₄/air flame was perturbed successively by incident and reflected shock waves reflected from a planar or concave wall. The two-dimensional axi-symmetric Navier–Stokes equations with detailed chemistry were used. The computational results were qualitatively validated with experiments which were performed in a standard shock tube arrangement. Under the influence of the incident shock wave, a Richtmyer–Meshkov instability is induced in the flame, and the distorted flame finally takes the form of two separated elliptical burning bubbles in the symmetric cross plane. Then, under subsequent interactions with the shock wave reflected from the planar or the concave wall, the flame takes a mushroom-like shape. Transverse waves produced by the shock reflection from the concave wall can compress the flame towards the axis, and the focusing shock generated on the concave wall will lead to a larger mushroom-like flame than that induced by the planar reflection.      

Study on the Fuel Air Mixing Induced by a Shock Wave Propagating into a H2-Air Interface

2nd-order upwind TVD scheme was used to solve the laminar, fully Navier-Stokes equations. The numerical simulations were done on the propagation of a shock wave with Ma _S = 2 and 4 into a hydrogen and air mixture in a duct and a duct with a rearward step. The results indicate that a swirling vortex may be generated in the lopsided interface behind the moving shock. Meanwhile, the complex shock system is also formed in this shear flow region. A large swirling vortex is produced and the fuel mixing can be enhanced by a shock wave at low Mach number. But in a duct with a rearward step, the shock almost disappears in hydrogen for Ma _S = 2. The shock in hydrogen will become strong if Ma _S is large. Similar to the condition of a shock moving in a duct full of hydrogen and air, a large vortex can be formed in the shear flow region. The large swirling vortex even gets through the reflected shock and impacts on the lower wall. Then, the distribution of hydrogen behind the rearward step is divided into two regions. The transition from regular reflection to Mach reflection was observed as well in case Ma _S = 4.     

Numerical study of the oscillations induced by shock/shock interaction in hypersonic double-wedge flows

In this paper, the shock pattern oscillations induced by shock/shock interactions over double-wedge geometries in hypersonic flows were studied numerically by solving 2D inviscid Euler equations for a multi-species system. Laminar viscous effects were considered in some cases. Temperature-dependent thermodynamic properties were employed in the state and energy equations for consideration of the distinct change of the thermodynamic state. It was shown that the oscillation results in high-frequency fluctuations of heating and pressure loads over wedge surfaces. In a case with a relatively lower free-stream Mach number, the shock/shock interaction structure maintains a seven-shock configuration during the entire oscillation process. On the other hand, the oscillation is accompanied by a transition between a six-shock configuration (regular interaction) and a seven-shock configuration (Mach interaction) in a case with a higher free-stream Mach number. Numerical results also indicate that the critical wedge angle for the transition from a steady to an oscillation solution is higher compared to the corresponding value in earlier numerical research in which the perfect diatomic gas model was used.      

On the Chaotic Dynamics of a Spherical Pendulum with a Harmonically Vibrating Suspension

The equations of motion for a lightly damped spherical pendulum are considered. The suspension point is harmonically excited in both vertical and horizontal directions. The equations are approximated in the neighborhood of resonance by including the third order terms in the amplitude. The stability of equilibrium points of the modulation equations in a four-dimensional space is studied. The periodic orbits of the spherical pendulum without base excitations are revisited via the Jacobian elliptic integral to highlight the role played by homoclinic orbits. The homoclinic intersections of the stable and unstable manifolds of the perturbed spherical pendulum are investigated. The physical parameters leading to chaotic solutions in terms of the spherical angles are derived from the vanishing Melnikov–Holmes–Marsden (MHM) integral. The existence of real zeros of the MHM integral implies the possible chaotic motion of the harmonically forced spherical pendulum as a result from the transverse intersection between the stable and unstable manifolds of the weakly disturbed spherical pendulum within the regions of investigated parameters. The chaotic motion of the modulation equations is simulated via the 4th-order Runge–Kutta algorithms for certain cases to verify the analysis.     

On the pressure and stress singularities induced by steady flows of incompressible viscous fluids

Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, to date there has been relatively little explicit identification of stress singularities caused by fluid flows. In this study, stress and pressure singularities induced by steady flows of viscous incompressible fluids are asymptotically identified. This is done by taking advantage of an earlier result that the Navier-Stokes equations are locally governed by Stokes flow in angular corners. Findings for power singularities are confirmed by developing and using an analogy with solid mechanics. This analogy also facilitates the identification of flow-induced log singularities. Both types of singularity are further confirmed for two global configurations by applying convergence-divergence checks to numerical results. Even though these flow-induced stress singularities are analogous to singularities in solid mechanics, they nonetheless render a number of structural configurations singular that were not previously appreciated as such from identifications within solid mechanics alone.     

Shock induced Richtmyer-Meshkov instability in the presence of a wall boundary layer

An experimental investigation on gaseous mixing zones originated from the Richtmyer-Meshkov instability has been undertaken in a square cross section shock tube. Mass concentration fields, of one of the two mixing constituents, have been determined within the mixing zone when the shock wave passes from the heavy gas to the light one, from one gas to an other of close density, and from the light gas to the heavy one. Results have been obtained before and after the coming back of the reflected shock wave. The diagnostic method is based on the infrared absorption of one of the two constituents of the mixing zone. It is shown that the mixing zone is strongly deformed by the wall boundary layer. The consequence is the presence of strong gradients of concentration in the direction perpendicular to the shock wave propagation. Finally, it is pointed out that the mixing goes more homogeneous when the Atwood number tends to zero.     

Evolution of the Diffusion Mixing Layer of Two Gases upon Interaction with Shock Waves

A mathematical model of mechanics of a twovelocity twotemperature mixture of gases is developed. Based on this model, evolution of the mixing layer of two gases with different densities under the action of shock and compression waves is considered by methods of mathematical simulation in the onedimensional unsteady approximation. In the asymptotic approximation of the full model, a solution of an initialboundary problem is obtained, which describes the formation of a diffusion layer between two gases. Problems of interaction of shock and compression waves with the diffusion layer are solved numerically in the full formulation. It is shown that the layer is compressed as the shock wave traverses it; the magnitude of compression depends on shockwave intensity. As the shock wave passes from the heavy gas to the light gas, the mixing layer becomes overcompressed and expands after shockwave transition. The wave pattern of the flow is described in detail. The calculated evolution of the mixinglayer width is in good agreement with experimental data.     

Investigation of high-speed combustible gas ignited by a hot gas jetproduced in the shock tube

The high-speed combustible gas ignited by a hot gas jet, which is induced by shock focusing, was experimentally investigated. By use of the separation mode of shock tube, the test section of a single shock tube is split into two parts, which provide the high-speed flow of combustible gas and pilot flame of hot gas jet, respectively. In the interface of two parts of test sections the flame of jet was formed and spread to the high-speed combustible gas. Two kinds of the ignitions, 3-D “line-flame ignition” and 2-D “plane-flame ignition”, were investigated. In the condition of 3-D “line-flame ignition” of combustion, thicker hot gas jet than pure air jet, was observed in schlieren photos. In the condition of 2-D “plane-flame ignition” of combustion, the delay time of ignition and the angle of flame front in schlieren photos were measured, from which the velocity of flame propagation in the high-speed combustible gas is estimated in the range of 30–90m/s and the delay time of ignition is estimated in the range of 0.12–0.29ms. PACS 47.40.Nm; 82.40.FpPart of this paper was presented at the 5th International Workshop on Shock/Vortex Interaction, Kaohsiung, October 27–31, 2003.     

On the Long-Time Behaviour of Solutions to the Navier–Stokes–Fourier System with a Time-Dependent Driving Force

In this paper, we give a complete characterization of the asymptotic behaviour of solutions to the Navier–Stokes–Fourier system. We show that either the driving force behaves asymptotically as a gradient of a scalar function, in which case any solution tends to a static state, or the total energy goes to infinity with growing time.     

Wave Motions in a Two-Layer Fluid Driven by Small Oscillations of a Cylinder Intersecting the Interface

The plane problem of steady-state small oscillations of a horizontal cylinder located at the interface between two fluids of different densities and indefinite depth is considered in the linear formulation. Boundary integral equations for the surface source distribution are derived. The behavior of the distributed singularities at points of intersection of the body contour and the interface is investigated. The problem of oscillations of a circular cylinder is solved by the multipole expansion method. The apparent mass and damping coefficients of the radiation problem and the reflection coefficient of the problem of scattering of an impinging wave by a floating body are calculated.