Unsteady shock generated vorticity

Abstract. A two-dimensional theory is developed for the vorticity just downstream of a curved, unsteady shock wave. By utilizing Crocco's equation, an explicit formula is obtained for the vorticity that does not require a perfect gas and that holds for arbitrary conditions upstream of the shock wave. The analysis is applied to the flow just downstream of the reflected shock that occurs in a single-Mach reflection pattern. Flow conditions are based on an interferometric photograph of Ben-Dor and Glass (1978). In this case, the reflected shock is weak everywhere from its upstream intersection with the wall to the triple point. The vorticity has a singularity and a change of sign near the triple point that indicates the presence of a weak shear layer downstream of this location. Received 20 December 1999 / Accepted 19 March 2000     

The flow gradients in the vicinity of a shock wave for a thermodynamically imperfect gas

Supersonic rotational planar and axisymmetric flows of a non-viscous, non-heat-conductive gas with arbitrary thermodynamic properties in the vicinity of a steady shock wave are studied. The differential equations describing the gas flow upstream and downstream of the discontinuity surface and the dynamic compatibility conditions at this discontinuity are used. The gas flow non-uniformity in the shock vicinity is described by the spatial derivatives of the gasdynamic parameters at a point on the shock surface. The parameters are the gas pressure, density, and velocity vector. The derivatives with respect to the directions of the streamline and normal to it, and of the shock surface and normal to it, are considered. Spatial derivatives of all gasdynamic parameters are expressed through the flow non-isobaric factor along the streamline, the streamline curvature, and the flow vorticity and non-isoenthalpy factors. An algorithm for determining these factors of the gas flow downstream of a shock wave is developed. Example calculations of these factors for imperfect oxygen and thermodynamically perfect gas are presented. The influence coefficients of the upstream flow factors on the downstream flow factors are calculated. The gas flow in the vicinity of the shock is described by the isolines of gasdynamic parameters. Uniform planar and axisymmetric flows at different distances from the axis of symmetry are examined; the isobars, isopycnics, isotachs and isoclines are used to characterize the downstream flow behind a curved shock in an imperfect gas.     

Curved shock theory

Curved shock theory (CST) is introduced, developed and applied to relate pressure gradients, streamline curvatures, vorticity and shock curvatures in flows with planar or axial symmetry. Explicit expressions are given, in an influence coefficient format, that relate post-shock pressure gradient, streamline curvature and vorticity to pre-shock gradients and shock curvature in steady flow. The effect of pre-shock flow divergence/convergence, on vorticity generation, is related to the transverse shock curvature. A novel derivation for the post-shock vorticity is presented that includes the effects of pre-shock flow non-uniformities. CST applicability to unsteady flows is discussed.     

Numerical studies of fluid flow through tubes with double constrictions

The flow fields in the neighbourhood of double constrictions in a circular cylindrical tube were studied numerically. The effects on the streamline, velocity and vorticity distributions as the flow passes through the constrictions in the tube were studied in the Reynolds number range 5–200. Double constrictions with dimensionless spacing ratios of 1, 2, 3 and ∞ were studied for a 50% constriction. It is noted that when the Reynolds number is below 10, no recirculation region is formed in the above constricted flow. For Reynolds numbers greater than 10, a recirculation region forms downstream of each of the constrictions. For constriction spacing ratios of 1, 2, and 3, when the Reynolds number is high, a recirculation region spreads between the valley of the constrictions. The recirculation region formed between the two constrictions has a diminishing effect on the generation of wall vorticity near the second constriction area. In general, the peak value of wall vorticity is found slightly upstream of each of the constrictions. When the Reynolds number is increased, the peak wall vorticity value increases and its location is moved upstream. Maximum wall vorticity generated by the first constriction is found to be always greater than the maximum wall vorticity generated by the second constriction. The extent of this spreading of the recirculation region from the first constriction and its effects on the second constriction depend on the constriction spacing ratio and the flow Reynolds number.     

Behavior of the vorticity vector in supersonic axisymmetric swirl flows behind a detonation wave

The behavior of the vorticity vector on a discontinuity surface arising in a supersonic nonuniform combustible gas flow with the formation of a shock or detonation wave is studied. In the general case, it is a vortex flow with prescribed distributions of parameters. It is demonstrated that the ratio of the tangential component of vorticity to density remains continuous in passing through the discontinuity surface, while the quantities proper become discontinuous. Results calculated for flow vorticity behind a steady-state detonation wave in an axisymmetric supersonic flow of a combustible mixture of gases are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 15–21, November–December, 2007     

Experimental investigation of viscoelastic drop deformation in Newtonian matrix at high capillary number under simple shear flow

The steady state deformation of a viscoelastic drop (Boger fluid) in a Newtonian liquid at high capillary number under simple shear flow is investigated by direct visualization using a specially designed Couette apparatus which enables visualization from two perpendicular directions. Two drop deformation modes are found: (1) Mode I – drop deformation in the flow direction and (2) Mode II – drop deformation in the vorticity direction. The drop deformation mode depends on the relative strength of the elastic contribution to viscous contribution. If the elastic contribution is weak compared to the viscous contribution, the drop elongates in the flow direction via Mode I. If the elastic contribution is strong, the drop elongates in the vorticity direction via Mode II. The drop size also affects the drop deformation. At the same capillary number, bigger drops have larger deformations than smaller drops.     

Baroclinic vortex sheet production by shocks and expansion waves

Vortex sheet production by shocks and expansion waves refracting at a density discontinuity was examined and compared using an analytical solution and numerical simulations. The analytical solution showed that with a small exception, vortex sheet strength is generally stronger in fast/slow shock refractions. In contrast, expansion waves generated a stronger vortex sheet in slow/fast refractions. This difference results in larger vorticity deposited by shocks in fast/slow refractions and by expansion waves in slow/fast refractions. Shock refractions become irregular and the analytical solution fails when either incident, transmitted or reflected shock, exceeded the angle limit for an attached shock. To investigate vortex sheet production outside the range of analytical solutions and to verify the applicability of the planar-interface analytical solution to a curved interface, shock refraction through a sinusoidal interface was numerically simulated in the shock frame of reference. It is found that variation in the local incidence angle along the curved interface creates pressure waves that affect the level of deposited vorticity. This contributes to the difference between predictions from local analysis and numerical computation. Furthermore, an interesting behavior of the shock and expansion wave-deposited vorticity in supersonic ramp flow was discovered. When the high- and low-density streams were swapped, while keeping the incident flow Mach numbers constant, a vortex sheet of equal magnitude but of opposite sign was generated.     

流场非均匀性对非平面激波诱导的Richtmyer-Meshkov不稳定性影响的数值研究

基于Navier-Stokes方程组,采用可压缩多介质黏性流动和湍流大涡模拟程序MVFT (multi-viscousflow and turbulence),模拟了均匀流场与初始密度呈现高斯函数分布的非均匀流场中马赫数为1.25的非平面激波加载初始扰动air/SF6界面的Richtmyer-Meshkov (RM)不稳定性现象。数值模拟结果表明,初始流场非均匀性将会影响非平面激波诱导的RM不稳定性演化过程。反射激波加载前,非平面激波导致的界面扰动振幅随着流场非均匀性增强而增大;反射激波加载后,非均匀流场与均匀流场条件下的界面扰动振幅差异有所减小。进一步,定量分析流场中环量分布及脉动速度统计量揭示了前述规律的原因。此外,还与平面激波诱导的RM不稳定性进行了简单对比,发现由于非平面激波波阵面区域的涡量与激波冲击界面时产生的涡量的共同作用,使得非平面激波与平面激波诱导的界面失稳过程存在差异。     

Vorticity produced by shock wave diffraction

Numerical simulations with a monotonicity preserving flow solver have been performed to study shock diffraction phenomena and shock wave generated vorticity. The computations were performed using the conservative Finite Element Method-Flux Corrected Transport (FEM-FCT) scheme, which has been shown to have an excellent predictive capability for various compressible flows with both strong and weak shocks. An adaptive unstructured methodology based on adapting to high density and entropy gradients was used in conjunction with a conservative shock-capturing scheme to adequately resolve strong and weak flowfield gradients. The chief interest was the formation of vorticity arising from shock wave propagation over a sharp corner and the high accuracy and resolution of the interacting compressible wave features. Numerical simulations were compared with previous experimental results and exhibited remarkably good agreement in terms of compressible wave propagation, as well as vorticity development and transport. The computations also allowed insight into the fundamental fluid dynamics, specifically shock diffraction, vortex convection and shock-vortex interactions.     

Vorticity and its rate of change just downstream of a curved shock

A theory is developed for the vorticity and its substantial derivative just downstream of a curved shock wave, the resulting formulas are exact, algebraic, and explicit. Analysis is for a cylinder-wedge or sphere-cone body, at zero incidence, whose downstream half-angle is θ_b. Derived formulas directly depend only on the ratio of specific heats, γ, the freestream Mach number, M ₁, the local slope and curvature of the shock, and the dimensionality parameter, σ, which is zero for a two-dimensional shock and unity for an axisymmetric shock. In turn, the slope and curvature depend on γ, M ₁, and θ_b. Numerical results are provided for a bow shock in which θ_b is 5°, 10°, or 15°, M ₁ is 2, 4, or 6, and γ = 1.4. There is little dependence on the half angle but a strong dependence on the freestream Mach number and on dimensionality. For vorticity and its substantial derivative, the dimensionality dependence gradually decreases with increasing Mach number. In comparison to the two-dimensional case, an axisymmetric shock generates considerable vorticity in a region relatively close to the symmetry axis. Moreover, the magnitude of the vorticity, in this region, is further enhanced in the flow downstream of the shock. This dimensionality difference in vorticity and its substantial derivative is attributed to the three-dimensional relief effect in an axisymmetric flow.

---

     

Numerical simulation and PIV study of compressible vortex ring evolution

Formation and evolution of a compressible vortex ring generated at the open end of a short driver section shock tube has been simulated numerically for pressure ratios (PR) of 3 and 7 in the present study. Numerical study of compressible vortex rings is essential to understand the complicated flow structure and acoustic characteristics of many high Mach number impulsive jets where simultaneously velocity, density and pressure fields are needed. The flow development, incident shock formation, shock diffraction, vortex ring formation and its evolution are simulated using the AUSM+ scheme. The main focus of the present study is to evaluate the time resolved vorticity field of the vortex ring and the shock/expansion waves in the starting jet for short driver section shock tubes—a scenario where little data are available in existing literature. An embedded shock and a vortex induced shock are observed for PR =  7. However the vortex ring remains shock free, compact and unaffected by the trailing jet for PR =  3. Numerical shadowgraph shows the evolution of embedded shock and shock/expansion waves along with their interactions. The velocity and vorticity fields obtained from simulation are validated with the particle image velocimetry results and these data match closely. The translational velocity of the vortex ring, velocity across the vortex and the centre line velocity of the jet obtained from simulation also agree well with the experimental results.     

Determination of the vorticity on a wing of small aspect ratio in a hypersonic stream

The theory of a thin shock layer [1–3] is used to obtain a formula for calculating the component of the vorticity in the direction of the flow on a wing of small aspect ratio in a hypersonic gas stream. It is shown that for definite shapes of the wing and flow regimes zones may occur with large local values of the vorticity, which, as is well known, have a significant influence on the structure of the flow field.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 175–178, September–October, 1980.     

Existence and Stability of Multidimensional Transonic Flows through an Infinite Nozzle of Arbitrary Cross-Sections

We establish the existence and stability of multidimensional steady transonic flows with transonic shocks through an infinite nozzle of arbitrary cross-sections, including a slowly varying de Laval nozzle. The transonic flow is governed by the inviscid potential flow equation with supersonic upstream flow at the entrance, uniform subsonic downstream flow at the exit at infinity, and the slip boundary condition on the nozzle boundary. Our results indicate that, if the supersonic upstream flow at the entrance is sufficiently close to a uniform flow, there exists a solution that consists of a C ^1,α subsonic flow in the unbounded downstream region, converging to a uniform velocity state at infinity, and a C ^1,α multidimensional transonic shock separating the subsonic flow from the supersonic upstream flow; the uniform velocity state at the exit at infinity in the downstream direction is uniquely determined by the supersonic upstream flow; and the shock is orthogonal to the nozzle boundary at every point of their intersection. In order to construct such a transonic flow, we reformulate the multidimensional transonic nozzle problem into a free boundary problem for the subsonic phase, in which the equation is elliptic and the free boundary is a transonic shock. The free boundary conditions are determined by the Rankine–Hugoniot conditions along the shock. We further develop a nonlinear iteration approach and employ its advantages to deal with such a free boundary problem in the unbounded domain. We also prove that the transonic flow with a transonic shock is unique and stable with respect to the nozzle boundary and the smooth supersonic upstream flow at the entrance.     

DIAMOND PORT JET INTERACTION WITH SUPERSONIC FLOW

Interaction flow field of the sonic air jet through diamond shaped orifices at different incidence angles （10 degrees, 27.5 degrees, 45 degrees and 90 degrees） and total pressures （0.10 MPa and 0. 46 MPa） with a Mach 5.0 freestream was studied experimentally. A 90 degrees circular injector was examined for comparison. Crosssection Mach number contours were acquired by a Pitot-cone five-hole pressure probe. The results indicate that the low Mach semicircular region close to the wall is the wake region. The boundary layer thinning is in the areas adjacent to the wake. For the detached case, the interaction shock extends further into the freestream, and the shock shape has more curvature, also the low-Mach upwash region is larger. The vortices of the plume and the height of the jet interaction shock increase with increasing incidence angle and jet pressure. 90 degrees diamond and circular injector have stronger plume vorticity, and for the circular injector low-Mach region is smaller than that for the diamond injector. Tapered ramp increases the plume vorticity, and the double ramp reduces the level of vorticity. The three-dimensional interaction shock shape was modeled from the surface shock shape, the center plane shock shape, and crosssectional shock shape. The shock total pressure was estimated with the normal component of the Mach number using normal shock theory. The shock induced total pressure losses decrease with decreasing jet incidence angle and injection pressure, where the largest losses are incurred by the 90 degrees, circular injector.     

One Singularly Perturbed Problem of Turbulent Gas Dynamics

The problem of the interaction of a Prandtl–Mayer wave with a shear layer is solved using the small parameter method for the case where the flow vorticity in the shear layer is small. A direct expansion is constructed and its inadequacy at large distances from the vortex layer is proved. The strained coordinate method is used to obtain a uniformly adequate expansion. It is shown that for certain velocity distributions in the shear layer, the characteristics in the reflected simple wave resulting from the interaction intersect each other and a shock arises in the flow. There coordinates of the shock origin and the function describing the shock shape are obtained.     

On the interaction between the shock wave attached to a wedge and freestream disturbances

We present a study of the interaction of small amplitude, unsteady, freestream disturbances with a shock wave induced by a wedge in supersonic flow. These disturbances may be acoustic waves, vorticity waves, or entropy waves (or indeed a combination of all three). Their interactions then generate behind the shock disturbances of all three classes, an aspect that is investigated in some detail. Also, the possibility of enhanced mixing owing to additional vorticity produced by the shock-body coupling is investigated. It is shown that disturbances behind the shock may either decay downstream, or alternatively experience sustained oscillations. The precise regimes under which either behaviour is found are stated.     

Particle image velocimetry study of the shock-induced single mode Richtmyer–Meshkov instability

Full field particle image velocimetry (PIV) measurements are obtained for the first time in Richtmyer–Meshkov instability shock tube experiments. The experiments are carried out in a vertical shock tube in which the light gas (air) and the heavy gas (SF₆) flow from opposite ends of the shock tube driven section and exit through narrow slots at the interface location. A sinusoidal perturbation is given to the interface by oscillating the shock tube in the horizontal direction. Richtmyer–Meshkov instability is then produced by the interaction with a weak shock wave (M _s  = 1.21). PIV measurements are obtained by seeding the flow with 0.30 μm polystyrene Latex spheres which are illuminated using a double-pulsed Nd:YAG laser. PIV measurements indicate the vorticity to be distributed in a sheet-like distribution on the interface immediately after shock interaction and that this distribution quickly rolls up into compact vortices. The integration of the vorticity distribution over one half wave length shows the circulation to increase with time in qualitative agreement with the numerical study of Peng et al. (Phys. Fluids, 15, 3730–3744, 2003).     

Application of Particle Image Velocimetry and Reference Image Topography to jet shock cells using the hydraulic analogy

This paper examines the shock cell structure, vorticity and velocity field at the exit of an underexpanded jet nozzle using a hydraulic analogy and the Reference Image Topography technique. Understanding the flow in this region is important for the mitigation of screech, an aeroacoustic problem harmful to aircraft structures. Experiments are conducted on a water table, allowing detailed quantitative investigation of this important flow regime at a greatly reduced expense. Conventional Particle Image Velocimetry is employed to determine the velocity and vorticity fields of the nozzle exit region. Applying Reference Image Topography, the wavy water surface is reconstructed and when combined with the hydraulic analogy, provides a pressure map of the region. With this approach subtraction of surfaces is used to highlight the unsteady regions of the flow, which is not as convenient or quantitative with conventional Schlieren techniques. This allows a detailed analysis of the shock cell structures and their interaction with flow instabilities in the shear layer that are the underlying cause of jet screech.     

脉冲激光等离子体与正激波相互作用的PIV实验研究

脉冲激光等离子体与超声速流场相互作用在飞行器减阻隔热、点火助燃等方面具有重要的应用价值.纹影实验方法只能定性或半定量地反映流动状态.为定量研究速度分布和旋涡结构,针对激光等离子体及其与正激波相互作用过程开展粒子图像测速PIV实验研究.在激波管实验平台上建立了纳秒脉冲激光能量沉积系统和PIV测量系统,通过定量测量,探明了激光等离子体引致的激光空气泡以及热核的流动特性,揭示了激光等离子体在正激波冲击下的流动特性与演化规律,并给出了激光能量大小和位置对相互作用过程的影响.结果表明:激光空气泡内的速度分布在激光入射方向上并不关于击穿点对称,而是在靠近激光入射方向一侧的流速略大于远离激光入射方向一侧;斜压导致热核在演化初期产生涡环,后期则由剪切主导;正激波与激光空气泡界面、热核界面相互作用时,产生斜压涡量,当激光能量为87.8 mJ、正激波马赫数1.4时,热核在正激波作用下产生的涡量比在静止空气中演化时大1个数量级;激光与正激波相互作用的关键过程是热核在正激波冲击下演化成涡环,在激波波前注入激光能量能够获得更加显著的涡环.