Supersonic low-concentration dusty-gas flow past a plane cylinder in the presence of an oblique shock wave interacting with the bow shock

The motion of an inertial dispersed admixture near a plane cylinder immersed in a steady-state hypersonic dusty flow in the presence of an oblique shock wave interacting with the bow shock is considered. It is assumed that the free-stream particle mass concentration is small and the particles do not affect the carrier flow. The III and IV shock wave interaction regimes are considered. The gas flow parameters in the shock layer are calculated from the numerical solution of the full Navier-Stokes equations for the perfect gas. A TVD second-order finite-difference scheme constructed on the basis of a finite volume method is used. For calculating the dispersed-phase parameters, including the concentration, the full Lagrangian method is used. On a wide range of variation of the particle inertia parameters, the patterns of the particle trajectories, velocity, concentration, and temperature in the shock layer are studied. The possibility of aerodynamic focusing of the particles behind the shock wave intersection point and the formation of narrow beams with a high particle concentration is revealed. These beams impinge on the cylinder surface and result in a sharp increase in the local heat fluxes. The maximal possible increase in the heat fluxes caused by the particles colliding with the cylinder surface is estimated for the flows with and without the incident oblique shock wave.     

Three-dimensional thin viscous shock layer in the absence of planes of symmetry in the flow

The basic laws of viscous homogeneous gas flow at high supersonic speeds past smooth blunt bodies with a permeable surface are investigated within the framework of the thin viscous shock layer model. An efficient numerical method of solving these equations, which makes it possible to consider cases of flow past bodies at angles of attack and slip, when there are no planes of symmetry in the flow, is proposed. Some results of calculating the flow past a triaxial ellipsoid with an axial ratio of 103n73 at angles of attack =0–45° and slip angles =0–45° over a broad interval of Reynolds numbers are presented as an example. The effect of the principal determining parameters of the problem on the flow structure in the shock layer and the surface friction and heat transfer coefficients is analyzed. An expression for calculating the heat fluxes to the impermeable surface of smooth blunt bodies in a supersonic homogeneous viscous gas flow over a broad interval of Reynolds numbers is proposed on the basis of the solutions obtained and the results of other authors.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 150–158, March–April, 1989.     

Interference between a blunt edge shock layer and an impinging inclined shock at low Reynolds numbers

The interference between the shock layer on a cylinder modeling the leading edge of an air intake and an impinging plane inclined shock is investigated experimentally and numerically for a Reynolds number Re₀=32. The low-pressure wind tunnel experiments made it possible to visualize the flow and determine the local heat transfer in the presence of interference. The corresponding flow regimes were calculated numerically within the framework of the system of Navier-Stokes equations by the through-calculation method. The principal properties of the distribution of the flow characteristics for a low value of the Reynolds number were obtained for various types of interference and the differences with respect to the previously investigated interference regimes for high Reynolds numbers were examined.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 166–171, January–February, 1993.     

Unsteady convective surface heat flux measurements on cylinder for CFD code validation

The surface convective heat transfer rates to a cylinder have been measured using platinum thin film gauges in a shock tube and the results have been used to validate two numerical codes. The investigations have been carried out at different incident shock Mach numbers and flow Reynolds numbers. The measured and simulated results give an insight into the transient flow fields around the model in the shock tube and are a valuable means of complementing the numerical and experimental techniques used in this study.Received: 27 August 2003, Accepted: 11 November 2003, Published online: 17 February 2004Correspondence to: T. Saito     

Interaction of a centered isentropic compression wave with the bow shock ahead of a body in a rarefied gas

Rarefied gas flow with a centered isentropic compression wave is investigated using direct Monte Carlo simulation of the solution of the Boltzmann equation. For monatomic gas flow the pattern of formation of a suspended compression shock near the geometric center of the compression wave is considered. The flow pattern is compared with the results obtained within the framework of gas dynamics. For a diatomic gas the interference of a centered compression wave with the bow shock ahead of a cylinder is investigated. The dependence of the pressure and the heat transfer to the surface on the Reynolds number and the wave center position relative to the cylinder center is analyzed. The results are compared with those of numerical simulation of the Euler and boundary-layer equations.     

Multicomponent chemically-reacting turbulent viscous shock layer on a catalytic surface

Turbulent flows past blunt bodies at high supersonic speeds are mainly investigated within the framework of the boundary layer model. However, even at large Reynolds numbers owing to the strong entropy gradient on the lateral surface it becomes necessary to take boundary layer corrections into account in the higher approximations [1]. The use of viscous shock layer theory makes it possible to obtain fairly accurate results over a broad interval of variation of the Reynolds numbers without organizing iterations with respect to vorticity and displacement thickness. The nonequilibrium nature of both homogeneous and heterogeneous catalytic reactions is taken into account. The results obtained are compared with the experimental data [2, 3]. Previously, in [4, 5] turbulent flow was investigated within the framework of viscous shock layer theory in the case of equilibrium homogeneous reactions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 144–149, March–April, 1989.     

PLIF thermometry in shock tunnel flows using a Raman-shifted tunable excimer laser

Planar laser-induced fluorescence is performed in a free-piston shock tunnel by using a Raman-shifted tunable excimer laser to excite nitric oxide molecules in the flow. Two different flowfields are examined to test the difficulties associated with applying the technique to shock tunnels: the bluff body flow produced by a 25 mm diameter cylinder; and the oblique shock and expansion fan produced by a 35° half-angle wedge. For the cylinder, the maximum flow enthalpy was limited to 4.1 MJ kg due to high flow luminosity which is produced by metallic contaminants in the flow. A reflective filter is used to reduce the influence of flow luminosity making these measurements feasible. Freestream temperature measurements are in excellent agreement with those predicted from numerical flow calculations. Large uncertainties were observed for the high-temperature post-shock results. Several higher enthalpy shots (14 MJ kg) were also performed with the wedge and showed an insignificant amount of contaminant emission. Received 5 June 1996 / Accepted 8 February 1998     

Investigation of a hypersonic viscous shock layer on rotating axisymmetric bodies in the presence of blowing

The hypersonic flow of a laminar stream of viscous compressible gas past blunt axisyrametric bodies rotating about the longitudinal axis is considered. It is assumed that gas blows from the surface of the body. The solution of the problem is obtained by a finite-difference method in a wide range of Reynolds numbers and blowing and rotation parameters. Some results of the calculations characterizing the effect of the rotation on the velocity and temperature profiles across the shock layer, on the friction and heat transfer coefficients, and the shock wave separation are given for the neighborhood of the stagnation point. For large Reynolds numbers and strong blowing an analytic solution of the problem is found in an approximation of two inviscid layers separated by a contact surface. The calculations are made for the flow past a sphere and a paraboloid and it is shown that in the presence of rotation the maximum of the heat flux is shifted from the stagnation point onto the side surface of the body. The dependence of the pressure distribution, the heat flux, and the friction coefficient is investigated for cases of constant and variable blowing over the contour of the body.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 106–114, January–February, 1986.     

Interference between a centered rarefaction wave and an inclined shock wave

The gas flow in the zone of interaction between an oblique shock and a centered isentropic rarefaction wave is studied using the direct statistical simulation method for solving the Boltzmann equation. The data of calculations of the shock and rarefaction wave structures, flow fields, and streamlines are given for the free-stream Mach number M_∞ = 6, 4 and 2. The formation of the interaction zone is simulated by a gas flow past a double-plane wedge in which the break of the generating line leads to formation of the centered isentropic rarefaction wave. The results of calculations of this flow in solving the Boltzmann equation are given in the Euler approximation.     

Investigation of the hypersonic viscous shock layer around blunt bodies in a nonuniform flow

Unseparated viscous gas flow past a body is numerically investigated within the framework of the theory of a thin viscous shock layer [13–15]. The equations of the hypersonic viscous shock layer with generalized Rankine-Hugoniot conditions at the shock wave are solved by a finite-difference method [16] over a broad interval of Reynolds numbers and values of the temperature factor and nonuniformity parameters. Calculation results characterizing the effect of free-stream nonuniformity on the velocity and temperature profiles across the shock layer, the friction and heat transfer coefficients and the shock wave standoff distance are presented. The unseparated flow conditions are investigated and the critical values of the nonuniformity parameter a_k [10] at which reverse-circulatory zones develop on the front of the body are obtained as a function of the Reynolds number. The calculations are compared with the asymptotic solutions [10, 12].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 154–159, May–June, 1987.     

Hypersonic viscous shock layer in a swirling gas flow on a permeable surface

A hypersonic swirling flow of viscous compressible gas past rotating axisymmetric blunt bodies is considered, its velocity vector being parallel to the axis of rotation of the body. The body surface is assumed permeable, while, in the general case, the gas is not injected (drawn off) along the normal to the body surface. An analytic solution of the problem, valid at small Reynolds numbers, is found in the first approximation of the integral method of successive approximations. On the basis of the results of the numerical solution, obtained in a wide range of variation of the determining parameters of the problem, we investigate the influence of the swirling of the free-stream flow, the angular velocity of rotation of the body, the Reynolds number and the injection (suction) parameter on the structure of the compressed layer, and the coefficients of friction and heat transfer on the body surface. The influence of the swirling of the flow on the nature of the asymptotic behavior of the viscous shock layer equations at large Reynolds numbers is studied. It is shown that the presence of a nonzero peripheral component for the velocity vector of the gas in the shock layer can lead to a qualitative change in the nature of the flow. Deceased Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 27–37, November–December, 1986. The authors thank G. G. Chernyi for his useful discussion of the results of the work.     

Heat flux measurement over a cone in a shock tube flow

This paper is the part 2 of our previous thin film heat transfer measurements. In the first report we measured time variations of heat flux over a cylinder placed in a shock tube flow and compared experimental results with CFD results, Saito et al. (Shock Waves 14:327–333, 2004). We report a result of heat transfer measurements over an 86° apex angle cone surface impinged by a Ms = 2.38 shock wave in air with distributed thin film transfer gauges along cone surface and its comparison with results of numerical simulations. We performed double exposure holographic interferometric observation, and also from the heat transfer measurement and numerical simulation, confirmed the presence of delayed transition from regular to Mach reflection over the cone. The numerical estimation of delayed transition distance from the apex agreed very well with experimental one.      

Vibrations of a cylinder in a uniform flow in the presence of a no-slip side-wall

A circular cylinder placed in a uniform flow, and that spans the entire length between two side walls, may experience either parallel or oblique vortex shedding depending on the end conditions. It was shown by Mittal and Sidharth (2014) that the spatio-temporal periodicity of the oblique vortex shedding results in constant-in-time force experienced by the cylinder. On the contrary, parallel vortex shedding leads to fluid force that fluctuates with time. The free vibrations of a circular cylinder, in the presence of a wall, are investigated. For comparison, computations with end walls, where a slip condition on velocity is specified, are also carried out. The Reynolds number, based on the diameter of the cylinder and free-stream speed of the flow, is Re=100. The initial condition for the free vibrations is the fully developed unsteady flow past a stationary cylinder with oblique shedding. It is found that as the amplitude of vibration of the cylinder builds up, the vortices shed from the cylinder align with its axis leading to parallel shedding. The response of the cylinder is associated with two branches: initial and lower. On the lower branch, the response of the cylinder is virtually identical from two- and three-dimensional computations. The flow as well as the response is different on the initial branch and outside the synchronization regime. Forced vibrations confirm the phenomena.     

Spontaneous symmetry-breaking in the 3-D wake of a long cylinder simulated by the lattice gas method and drag coefficient measurements

We present a direct numerical simulation by the lattice gas method of the three-dimensional non-stationary incompressible flow at a Reynolds number of 74, past a circular cylinder, with a uniform incident flow. We describe the three-dimensional structure and the time-evolution of the wake, which leads to an oblique vortex shedding situation.We also present early results on measurements of drag coefficients for spheres and cylinders at a Reynolds number of 20.     

Viscous shock layer on the surface of a blunt body in a diverging supersonic flow

The results are given of calculations of supersonic diverging flow past the spherical front of a body obtained using the Navier-Stokes equations. Flows at moderate Reynolds numbers are considered. A study is made of the influence of the nonuniformity of the oncoming flow on the flow field in the shock layer and on the distribution of the pressure, the friction coefficient, and the heat flux over the surface of the sphere.     

Asymptotic solution of the two-dimensional equations for a thin viscous rarefied shock layer on a cold surface

Hypersonic rarefied flow past blunt bodies is studied in the continuum-free-molecular transition regime. On the basis of an asymptotic analysis three rarefied gas flow patterns are established depending on the relation between the relevant parameters of the problem. In the first regime corresponding to a cold surface asymptotic solutions of the equations of a thin viscous shock layer are derived at low Reynolds numbers in the axisymmetric and plane cases. Simple analytical expressions for the pressure and the heat transfer and friction coefficients are obtained as functions of the freestream parameters and the body geometry. With decrease in the Reynolds number the coefficients approach the values corresponding to free-molecular flow. In this regime a similarity parameter for the hypersonic rarefied flow past bodies is determined. The asymptotic solutions are compared with numerical solutions and the results of direct statistical simulation by the Monte Carlo method.     

Shock tensors and shock polars

A refraction law for the velocity at an oblique shock in a compressible fluid is derived in dyadic form similar to that for refraction of light rays at an interface. The shock tensor embodies only the assumptions of conservation of mass and equality of tangential velocity components. Given the shock inclination and density ratio, a quadratic equation in the ratio of the flow speeds can be found with flow turning angle as a parameter. Analysis of the two solutions shows that they lie on a circle in the polar plane, a result independent of the equation of state or other conservation laws. If the density ratio is allowed to vary, a pencil of circles is generated in the hodograph plane ; or, equivalently a right, elliptic cone with two nappes appears in the three-space formed when the density ratio coordinate is added at right angles to the hodograph plane. The further requirements that momentum and energy be conserved taken together with weak restrictions on the functional form of the equation of state are sufficient to permit the development of a general theory of shock polars. The allowed shock states are seen to lie on the space curve formed by intersection of a surface called the Hugoniot cylinder with the elliptic cone. The projection of this space curve on the hodograph plane is the shock polar. The theory is applied to the special case of a polytropic gas by way of illustration.     

Study of unsteady flow past bodies behind shock waves

Several theoretical and experimental studies have been devoted to the problem of the nonstationary action of the stream behind a shock wave on bodies of varied shape. In particular, in [1], the pressure and density are calculated for flow about bodies of the more typical shapes in the initial stage of the process. The basic relations which accompany the interaction of shock waves are considered in [2, 3]. The analysis of the phenomena of diffraction of shock waves on the sphere, cylinder, and cone is presented in [4]. Problems of unsteady flow about a wing are examined in [5, 6]. A detailed review of the foreign studies on unsteady flow is given in [7]. Of great practical interest is the question of the time for flow formation and the magnitudes of the unsteady loads during this period. Experimental investigations have been made recently [8, 9] in which some criteria are presented for estimating the bow shock formation time for supersonic flow about the sphere and cylinder with flat blunting. However the question of the formation time of the stationary pressure on the body surface is not referred to in these studies and no relationship is shown between the transient position of the reflected wave and the corresponding unsteady pressure on the surface. Moreover, in [8] the dimensionless time criterion is determined very approximately, independently of the Mach number of the shock wave. The present study was undertaken with the object of determining the basic criteria which characterize unsteady flow about bodies behind a plane shock wave which has time-independent parameters, and clarification of the shock wave reflected from the body and the pressure on the surface of the body during the transient period. The most typical body shapes were studied: 1) a cylinder with flat face aligned with the stream; 2) a spherically-blunted cylinder; and 3) a cylinder transverse to the stream. The experiments were conducted in a conventional shock tube using the single-diaphragm scheme. The measurements of the pressure on the models and the velocity of the incident shock wave were made using the technique analogous to that of [10, 11]. A highspeed movie camera was used to record the pattern of the wave diffraction on the body. The Mach number of the incident shock wave varied in the range from M=1.5 to M≈6.0, which corresponded to a range of Mach numbers M_∞ of the stream behind the shock wave from 0.6 to 2.1. The calculations of the required gas dynamic parameters for high temperatures were made with account for equilibrium dissociation of the air on the basis of the data of [10, 12, 13]. The magnitude of the relative maximal shock wave standoff Δ at the stagnation point obtained in the present experiments was compared with the values of Δ from other studies. In the case of the flat-blunted cylinder it was in good agreement with the results of [8–14], and in the case of the spherically-blunted cylinder and the transverse cylinder it was in agreement with the results of [15].     

Coherent and turbulent processes in the bistable regime around a tandem of cylinders including reattached flow dynamics by means of high-speed PIV

The turbulent flow around two cylinders in tandem at the sub-critical Reynolds number range of order of 10⁵ and pitch to diameter ratio of 3.7 is investigated by using time-resolved Particle Image Velocimetry (TRPIV) of 1 kHz and 8 kHz. The bi-stable flow regimes including a flow pattern I with a strong vortex shedding past the upstream and the downstream cylinder, as well as a flow pattern II corresponding to a weak alternating vortex shedding with reattachment past the upstream cylinder are investigated. The structure of this “reattachment regime” has been analyzed in association with the vortex dynamics past the downstream cylinder, by means of POD and phase-average decomposition. These elements allowed interconnection among all the measured PIV planes and hence analysis of the reattachment structure and the flow dynamics past both cylinders. The results highlight fundamental differences of the flow structure and dynamics around each cylinder and provide the ‘gap’ flow nature between the cylinders. Thanks to a high-speed camera of 8 kHz, the shear-layer vortices tracking has been possible downstream of the separation point and the quantification of their shedding frequency at the present high Reynolds number range has been achieved. This issue is important regarding fluid instabilities involved in the fluid–structure interaction of cylinder arrays in nuclear reactor systems, as well as acoustic noise generated from the tandem cylinders of a landing gear in aeronautics.     

Asymptotic solution of the thin viscous shock layer equations at low Reynolds numbers for a cold surface

Hypersonic three-dimensional viscous rarefied gas flow past blunt bodies in the neighborhood of the stagnation line is considered. The question of the applicability of the gasdynamic thin viscous shock layer model [1] is investigated for the transition flow regime from continuum to free-molecular flow. It is shown that for a power-law temperature dependence of the viscosity coefficient T the quantity (Re)^1/(1+), where = ( – 1)/2 and is the specific heat ratio, is an important determining parameter of the hypersonic flow at low Reynolds numbers. In the case of a cold surface approximate asymptotic solutions of the thin viscous shock layer equations are obtained for noslip conditions on the surface and generalized Rankine-Hugoniot relations on the shock wave at low Reynolds numbers. These solutions give simple analytic expressions for the thermal conductivity and friction coefficients as functions of the determining flow parameters. As the Reynolds number tends to zero, the values of the thermal conductivity and friction coefficients determined by this solution tend to their values in free-molecular flow for an accommodation coefficient equal to unity. This tending of the thermal conductivity and friction coefficients to the free-molecular limit takes place for both two-and three-dimensional flows. The asymptotic solutions are compared with numerical calculations and experimental data.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 159–170. Original Russian Text Copyright © 2004 by Brykina.