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     

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.      

Lift force on rotating sphere at low Reynolds numbers and high rotational speeds

The lift force on an isolated rotating sphere in a uniform flow was investigated by means of a three-dimensional numerical simulation for low Reynolds numbers (based on the sphere diameter) (Re<68.4) and high dimensionless rotational speeds (Г5). The Navier-Stokes equations in Cartesian coordinate system were solved using a finite volume formulation based on SIMPLE procedure. The accuracy of the numerical simulation was tested through a comparison with available theoretical, numerical and experimental results at low Reynolds numbers, and it was found that they were in close agreement under the above mentioned ranges of the Reynolds number and rotational speed. From a detailed computation of the flow field around a rotational sphere in extended ranges of the Reynolds number and rotational speed, the results show that, with increasing the rotational speed or decreasing the Reynolds number, the lift coefficient increases. An empirical equation more accurate than those obtained by previous studies was obtained to describe both effects of the rotational speed and Reynolds number on the lift force on a sphere. It was found in calcttlations that the drag coefficient is not significantly affected by the rotation of the sphere. The ratio of the lift force to the drag force, both of which act on a sphere in a uniform flow at the same time, was investigated. For a small spherical particle such as one of about 100μm in diameter, even if the rotational speed reaches about 10^6 revolutions per minute, the lift force can be neglected as compared with the drag force.     

Computation of wave interference and relaxation of particles after passing of a shock wave

Interaction of a shock wave with a system of motionless or relaxing particles is numerically simulated. Regimes of the gas flow around these particles are described, and the influence of the initial parameters of the examined phenomenon on the flow pattern is analyzed. The drag coefficient of particles is calculated as a function of the Mach number behind the shock wave at a fixed Reynolds number. The dynamics of heat exchange for particles of different sizes (10 μm–1 mm) is determined, and the laws of thermal relaxation after passing of a shock wave over the system of particles are found. The times of thermal and velocity relaxation of particles are estimated as functions of the Reynolds number, and the predicted relaxation time is compared with the corresponding empirical dependences.     

Some further observations on the creeping motion of spheres through Ellis model fluids

Using the formulation of Hopke and Slattery, upper and lower bounds on the drag coefficient of a sphere moving slowly in Ellis model fluids have been calculated, over wide range of conditions, and compared with the suitable experimental data available in the literature. C _D drag coefficient - d diameter of sphere - El Ellis number - Re Reynolds number based on zero-shear viscosity - V terminal falling velocity of a sphere - X drag correction factor - Ellis model parameter - ₀ zero-shear viscosity - _1/2 Ellis model parameter     

Unsteady flow with separation behind a shock wave diffracting over curved walls

The unsteady separation of the compressible flow field behind a diffracting shock wave was investigated along convex curved walls, using shock tube experimentation at large length and time scales, complemented by numerical computation. Tests were conducted at incident shock Mach numbers of $M_{\hbox {s}} =$ 1.5 and 1.6 over a 100 mm radius wall over a dimensionless time range up to $\tau \le $ 6.45. The development of the near wall flow at $M_{\hbox {s}} =$ 1.5 has been described in detail and is very similar to that observed for slightly lower $\tau $ ’s at $M_{\hbox {s}} =$ 1.6. Computations were performed at wall radii of 100 and 200 mm and for incident shock Mach numbers from 1.5 up to and including Mach 2.0. Comparing dimensionless times for different size walls shows that for a given value of $\tau $ the flow field is very similar for the various wall radii published to date and tested in this study. Previously published results that were examined alongside the results from this study had typical values of $1.6 < \tau < 3.2$ . At the later times presented here, flow features were observed that previously had only been observed at higher Mach numbers. The larger length scales allowed for a degree of Reynolds number independence in the results published here. The effect of turbulence on the numerical and experimental results could not be adequately examined due to limitations of the flow imaging system used and a number of questions remain unanswered.     

Fluid Flows past Spheroids at Moderate Reynolds Numbers

Axisymmetric viscous fluid flows past spheroids are considered. The time-independent complete Navier-Stokes equations written in a spherical coordinate system are used for describing the flow. The problem is solved by the stabilization method on the basis of a variable direction scheme. The input coordinate system is transformed in order to construct a regular computational grid. As a result of the numerical investigation, the stream patterns of flow past elongated and oblate spheroids are obtained for various values of the determining parameters. Numerical values of the dimensions of the circulation zone and the drag coefficient are given for various values of the spheroid semi-axis ratio in the domain of moderate Reynolds numbers.     

Hyperbolic Approximation of the Navier-Stokes Equations for Viscous Mixed Flows

Simplified two-dimensional Navier-Stokes equations of the hyperbolic type are derived for viscous mixed (with transition through the sonic velocity) internal and external flows as a result of a special splitting of the pressure gradient in the predominant flow direction into hyperbolic and elliptic components. The application of these equations is illustrated with reference to the calculation of Laval nozzle flows and the problem of supersonic flow past blunt bodies. The hyperbolic approximation obtained adequately describes the interaction between the stream and surfaces for internal and external flows and can be used over a wide Mach number range at moderate and high Reynolds numbers. Examples of the calculation of viscous mixed flows in a Laval nozzle with large longitudinal throat curvature and in a shock layer in the neighborhood of a sphere and a large-aspect-ratio hemisphere-cylinder are given. The problem of determining the drag coefficient of cold and hot spheres is solved in a new formulation for supersonic air flow over a wide range of Reynolds numbers. In the case of low and moderate Reynolds numbers a drag reduction effect is detected when the surface of the sphere is cooled.     

脉冲激光与正激波相互作用过程和减阻机理的实验研究

纳秒脉冲激光具有峰值功率密度高、易于击穿空气形成等离子体这一突出优势,在降低超声速波阻方面具有重要应用价值.以深刻揭示减阻机理为目的,针对激光与正激波相互作用这一基本物理现象开展实验研究.发展高精度纹影技术以测量复杂激波结构,时间分辨率达到 30ns,空间分辨率达到 1mm;搭建快速~PIV 实验系统以定量测量流场速度和涡量,时间分辨率达到 500ns.探明了激光等离子体引致的球面激波和高温低密度区域特性,揭示了激光等离子体在正激波冲击下的流动特性与演化规律,并结合数值模拟结果阐明了脉冲激光等离子体降低超声速波阻的根本原因.研究表明:激光等离子体引致激波的初始马赫数随着激光能量而增大,形状由水滴形逐渐发展为球面形,传播速度随着时间降低,在50$\mu$s 后接近于声速;高温低密度区域初始近似于球形,而后从激光入射方向的下游开始失稳,形成尖刺结构;在正激波冲击下,高温低密度区域演化为上下对称的双涡环结构,尺寸随着激光能量而增大.涡的卷吸和逆流可改变飞行器头部激波结构,是流场重构的重要形式,引起飞行器表面压力的大幅降低,是引起超声速飞行器波阻降低的重要机理.      

Hypervelocity shock standoff on spheres in air

To provide data for the validation of computational fluid dynamics models, measurements of the shock standoff distance on spheres in hypervelocity flows have been made. Test flows of air at 8.7 and 9.7 km/s were generated in the X2 expansion tunnel fitted with a Mach 10 nozzle. High-speed video images were analysed with a least-squares shape-fitting algorithm. Assuming a spherical shock shape near the nose enabled increased resolution measurements beyond the native pixel size. Normalised shock standoff distances, $\Delta $ / $D$ , in the range 0.03–0.04 were measured, with sphere diameters, $D$ , of 40, 60 and 80 mm.