Momentum and Energy Transfer in a Shock Wave

The macroparameter profiles in a strong shock wave propagating in a single-component monatomic gas are investigated. The interaction between the molecules is described by the variable-diameter sphere model. Qualitative information concerning the shock wave parameters is obtained by direct statistical simulation (Monte-Carlo) method. Using the data obtained, simple approximate gas-velocity dependences of the stress and the heat flow are found and the linear Newton and Fourier relations are generalized.     

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.     

Interaction of a shock with a sphere suspended in a vertical shock tube

Shock wave interaction with a sphere is one of the benchmark tests in shock dynamics. However, unlike wind tunnel experiments, unsteady drag force on a sphere installed in a shock tube have not been measured quantitatively. This paper presents an experimental and numerical study of the unsteady drag force acting on a 80 mm diameter sphere which was vertically suspended in a 300 mm x 300 mm vertical shock tube and loaded with a planar shock wave of M _s = 1.22 in air. The drag force history on the sphere was measured by an accelerometer installed in it. Accelerometer output signals were subjected to deconvolution data processing, producing a drag history comparable to that obtained by solving numerically the Navier-Stokes equations. A good agreement was obtained between the measured and computed drag force histories. In order to interpret the interaction of shock wave over the sphere, high speed video recordings and double exposure holographic interferometric observations were also conducted. It was found that the maximum drag force appeared not at the time instant when the shock arrived at the equator of the sphere, but at some earlier time before the transition of the reflected shock wave from regular to Mach reflection took place. A negative value of the drag force was observed, even though for a very short duration of time, when the Mach stem of the transmitted shock wave relfected and focused at the rear stagnation point of the sphere.Received: 31 March 2003, Accepted: 7 July 2003, Published online: 2 September 2003     

Acceleration of a sphere behind planar shock waves

The paper presents the results of an investigation on the motion of a spherical particle in a shock tube flow. A shock tube facility was used for studying the acceleration of a sphere by an incident shock wave. Using different optical methods and performing experiments in two different shock tubes, the trajectory and velocity of a spherical particle were measured. Based upon these results and simple one-dimensional calculations, the drag coefficient of a sphere and shading effect of sphere interaction with a shock tube flow were studied.     

The effect of an unsteady drag force on the structure of a non-equilibrium region behind a shock wave in a gas-particle mixture

For numerical analysis of shock wave propagation in gas-particle mixtures, drag coefficients of a sphere in steady flows are generally used. However, it is shown both experimentally and numerically that a shock loaded solid sphere experiences unsteady drag forces. The paper describes a model of unsteady drag force and its effect on the structure of the non-equilibrium region behind a shock front traveling in a dusty gas. The results are compared with those obtained by using a steady drag coefficient and are discussed. It is demonstrated that the large drag force at the early stage of the interaction between shock-wave induced flow and a solid particle affects the flow structure that is obtained with a steady drag force.      

Propagation of converging spherical deformation waves in a heteromodular elastic medium

The unsteady one-dimensional boundary-value problem of shock deformation of a medium bounded by a sphere is solved. The propagation of converging deformation wave fronts in an elastic material with different tensile and compressive strengths is studied. A boundary condition is obtained that provides the formation of a converging spherical shock wave with constant velocity. The impact conditions on the boundary of the heteromodular sphere are determined that can lead to the formation of a transition zone (a spherical layer of constant density) between the compression and tension regions.     

Vortex in a Supersonic Flow and its Influence on Blunt Body Flow and Heat Transfer

The dissipation of a tip vortex generated by a rectangular wing and the stagnation pressure and temperature distributions within the vortex are experimentally investigated. The interaction of the vortex with the bow shock ahead of a sphere and the heat flux distribution over the spherical surface are studied. The experiments were carried out at the Mach number M=3 and the unit Reynolds numbers Re₁=1.1·10⁷ and 3.7·10⁷ 1/m.     

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.     

Propagation of acoustic waves and shock waves radiated by a sinusoidal pulsation of a sphere during a single period

A spherical sound wave is emitted by a sphere which executes a small sinusoidal pulsation of a single period at high frequency in an inviscid fluid. Nonlinear propagation of the waves is formulated as an initial boundary value problem and is analysed in detail. The governing equation is linear near the sphere, while it is a nonlinear hyperbolic equation in a far field. The nonlinearity has a significant effect there, leading to the formation of two shocks. The exact solution to match the near field solution can easily be obtained for the far field equation. The nonlinear distortion of waveform and the shock formation distance are evaluated from the representation of the solution with strained coordinates. The evolution and nonlinear attenuation of the two shock discontinuities are also examined by making use of the equal-areas rule. In its asymptotic form the entire profile is an N wave with a long tail.     

Shock wave standoff distance for a sphere slightly above Mach one

This paper describes a numerical simulation of bow shock formation ahead of a sphere at steady supersonic flow in the Mach number range of 1.025–1.20. Turbulent viscous flow results are presented using the Spalart–Allmaras turbulence model. The purpose of this study is to determine the shock standoff distance for a spherical projectile at slightly supersonic free flight speeds. Results are compared to experimental data, including double exposure holographic interferograms obtained from a 40 mm polycarbonate sphere launched by a light gas gun. The shock standoff distance was determined from the interferograms. The present numerical simulations were found to agree with previously published data, and reached down to M = 1.025—a range where almost no previously published data exists. The computed flow structure and shock wave locations agree well with recently obtained free-flight interferograms.     

Circulation flow at the front surface of a sphere in a supersonic wake

Several theoretical and experimental studies of supersonic flow past a blunt body located in the wake behind another body have been made [1–7]. It has been shown that a reverse-circulation flow can occur in the shock layer at the front surface. The possibility of such a flow forming depends on the nonuniformity of the freestream flow and the Reynolds number. This paper presents new results of the theoretical study of the structure of the shock wave at the front surface of such a sphere, obtained on the basis of numerical solution of Navier-Stokes equations. It is shown that for a fixed nonuniformity of the freestream flow, an increase in the Reynolds number and cooling of the surface of the body lead to the formation of a secondary vortex in the region where the contour of the body intersects the axis of symmetry. A study is made of the variations of the drag and heat transfer parameters over the front surface of a cooled and thermally insulated sphere. The possibility of numerical simulation of the flow on the basis of the Euler equations is discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 143–148, May–June, 1985.     

Unsteady drag on a sphere by shock wave loading

The dynamic drag coefficient of a sphere by shock wave loading is investigated numerically and experimentally. The diameter of the sphere is varied from 8 m to 80 mm in numerical simulation. The axisymmetric Navier-Stokes equations are solved on a fine grid, and the grid convergence of the drag coefficient is achieved. The numerical result is validated by comparing the experimental data of a 80 mm sphere, measured by the accelerometer in a vertical shock tube. It is found that the sphere experiences in the early interaction one order higher drag than in the steady state. A transient negative drag, mainly resulting from the focusing of shock wave on the rear side of the sphere, is observed only for high Reynolds number flows, and the drag becomes positive because of increased skin friction for low Reynolds number flows.Received: 10 March 2004, Accepted: 24 May 2004, Published online: 20 August 2004[/PUBLISHED]M. Sun: Send offprints requests to     

球面各向同性球体内的动态热应力集中

利用有限克尔变换得到了求面各向同笥球体的热冲击作用下的动态热应力解析表达式。从表达式中,可以看出球中向同性球体的动态热应力集中现象明显不同于各向同性球体。另外,所描述的动态热应力集中现象与文献「１,３」也有一定的区别。     

Reflection of a strong shock wave from a sphere and a cylinder

We examine the flow on the axis in the vicinity of the stagnation point for reflection of a strong plane shock wave (with uniform parameters behind the wave) from a sphere and a circular cylinder whose generators are parallel to the incident wave front.The small parameter method [1, 2] is used to obtain, in closed form, relations which define the time variation of the velocity profile, pressure, enthalpy, and reflected shock wave standoff.As the time t , these relations reduce to the known formulas [3, 4] which define the steady flow on the axis for the flow behind the incident shock wave about a body, if account is taken of the leading terms containing the small parameter.The solution is extended to the case in which account for equilibrium dissociation and ionization is necessary.Comparison of the results with measurement [5] of the reflected shock wave distance from a sphere, as a function of time, shows satisfactory agreement.     

Heat transfer from a sphere in the intermediate dynamics region of a rarefied gas

The results are given of the experimental study of convective heat transfer from a sphere in a low-density subsonic stream. Generalizing the results obtained and earlier known data for sub-and supersonic velocities, we suggest approximate formulas for calculating heat transfer from a sphere under any streamline flow conditions of a rarefied gas.Moscow. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 170-172, March-April, 1972.     

Heat transfer from a heated sphere in a rarefied gas at rest

We consider the problem of heat transfer from a slightly heated sphere in a resting rarefied gas. We assume that the Krook equation is valid in this case. Two forms of the basic equations are presented, and relations are given which are obtained as a result of calculations of the heat flux and the temperature jump at the sphere surface as a function of a parameter which is inversely proportional to the Knudsen number. The results obtained are compared with results given by the known approximate theories.In conclusion the author wishes to thank M. N. Kogan for proposing the problem and for numerous discussions.     

Flow in the neighborhood of the stagnation line on a blunt body traveling at supersonic speed through a zone with elevated temperature and vibrational energy of the molecules

The unsteady flow in the neighborhood of the stagnation line on a sphere traveling at supersonic speed through a plane layer of diatomic gas with elevated temperature and nonequilibrium excitation of the molecular vibrations is investigated. (The source of the inhomogeneity could be a gas discharge [1].) The problem is solved using the viscous shock layer model which makes it possible to take molecular transport processes into account and analyze the unsteady heat transfer. Such flows were previously calculated in [2] within the framework of the inviscid gas model.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i GazNo. 3, pp. 183–185, May–June, 1990.     

Transient conjugate mixed convection from a sphere in a porous medium saturated with cold pure or saline water

In this paper, a numerical investigation of the transient conjugate mixed convection flow about a sphere embedded in a porous medium saturated with pure or saline water is carried out. The effect of density extremum is considered by using the nonlinear dependence of density on the temperature. The salinity effects are considered by assuming uniform saline concentration over the domain considered. The direction of the natural convection is changed either to aiding or to opposing the upcoming flow direction simulating the sphere is either hot or cold relative to the surrounding temperature. Results show that the initial temperature differences as well as the saline concentration alter the transient heat transfer rate in conceivable degree. It was found that the heat capacity ratio between the sphere and the surrounding media has more significant effect on the calculated heat transfer rate than the thermal conductivity ratio. The study is performed by using six nondimensional parameters and results are discussed in detail. Received on 10 November 1997     

Establishment of circumfluence (steady-state flow) when a shock wave strikes a cylinder or a sphere

The problem of the diffraction of a shock wave at a stationary sphere or cylinder is considered. The finite-difference method proposed by S. K. Godunov [1, 2] is employed Numerical solutions are obtained for the stage of the diffraction of the shock wave and for the subsequent steady state of flow around the object (circumfluence). Cases of sub-, trans-, and supersonic flow behind the shock wave are considered. When strong shock waves undergo diffraction, zones of reverse flow appear in the neighborhood of the tail part of the obstacle.Moscow. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 97–103, September–October, 1972.     

Conjugate heat transfer in the unbounded flow of a viscoelastic fluid past a sphere

This work addresses the conjugate heat transfer of a simplified PTT fluid flowing past an unbounded sphere in the Stokes regime (Re = 0.01). The problem is numerically solved with the finite-volume method assuming axisymmetry, absence of natural convection and constant physical properties. The sphere generates heat at a constant and uniform rate, and the analysis is conducted for a range of Deborah (0 ≤ De ≤ 100), Prandtl (10⁰ ≤ Pr ≤ 10⁵) and Brinkman (0 ≤ Br ≤ 100) numbers, in the presence or absence of thermal contact resistance at the solid–fluid interface and for different conductivity ratios (0.1 ≤ κ ≤ 10). The drag coefficient shows a monotonic decrease with De, whereas the normalized stresses on the sphere surface and in the wake first increase and then decrease with De. A negative wake was observed for the two solvent viscosity ratios tested (β = 0.1 and 0.5), being more intense for the more elastic fluid. In the absence of viscous dissipation, the average Nusselt number starts to decrease with De after an initial increase. Heat transfer enhancement relative to an equivalent Newtonian fluid was observed for the whole range of conditions tested. The dimensionless temperature of the sphere decreases and becomes more homogeneous when its thermal conductivity increases in relation to the conductivity of the fluid, although small changes are observed in the Nusselt number. The thermal contact resistance at the interface increases the average temperature of the sphere, without affecting significantly the shape of the temperature profiles inside the sphere. When viscous dissipation is considered, significant changes are observed in the heat transfer process as Br increases. Overall, a simplified PTT fluid can moderately enhance heat transfer compared to a Newtonian fluid, but increasing De does not necessarily improve heat exchange.