Diffuse holographic interferometric observation of shock wave reflection from a skewed wedge

The pattern of shock wave reflection over a wedge is, in general, either a regular reflection or a Mach reflection, depending on wedge angles, shock wave Mach numbers, and specific heat ratios of gases. However, regular and Mach reflections can coexist, in particular, over a three-dimensional wedge surface, whose inclination angles locally vary normal to the direction of shock propagation. This paper reports a result of diffuse double exposure holographic interferometric observations of shock wave reflections over a skewed wedge surface placed in a 100 × 180 mm shock tube. The wedge consists of a straight generating line whose local inclination angle varies continuously from 30° to 60°. Painting its surface with fluorescent spray paint and irradiating its surface with a collimated object beam at a time interval of a few microseconds, we succeeded in visualizing three-dimensional shock reflection over the skewed wedge surface. Experiments were performed at shock Mach numbers, 1.55, 2.02, and 2.53 in air. From reconstructed holographic images, we estimated critical transition angles at these shock wave Mach numbers and found that these were very close to those over straight wedges. This is attributable to the flow three-dimensionality.      

Numerical simulation for nonstationary Mach reflection of a shock wave: A kinetic-model approach

A numerical simulation was performed for the process of formation of single Mach reflection on a wedge by solving a BGK type kinetic equation for the reduced distribution function with a finite difference scheme. The calculations were carried out for a shock Mach number 2.75 and wedge angle 25° in a monatomic gas, which corresponds to the conditions of single Mach reflection in the classical von Neumann theory. The calculations were performed for both diffuse and specular reflection of molecules at the wall surface. It is concluded that the diffuse reflection of molecules at the wall surface or the existence of the viscous or thermal layer is an essential factor for a nonstationary process at the initial stage of Mach reflection. Furthermore, the numerical results for diffuse reflection are found to simulate the experimental results very well, such as a transient process from regular reflection to Mach reflection along with shock propagation.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Analysis of the wave configuration resulting from the termination of an inverse Mach reflection

An analytical model for solving the wave configuration which is formed when an inverse Mach reflection terminates after its triple point collides with the reflecting surface has been developed. The predictions of the model were compared with available experimental results and good agreement was obtained.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Numerical simulations of shock wave reflection phenomena in non-stationary flows using regularized smoothed particle hydrodynamics (RSPH)

In this paper we wish to demonstrate to what extent the numerical method regularized smoothed particle hydrodynamics (RSPH) is capable of modelling shocks and shock reflection patterns in a satisfactory manner. The use of SPH based methods to model shock wave problems has been relatively sparse, both due to historical reasons, as the method was originally developed for studies of astrophysical gas dynamics, but also due to the fact that boundary treatment in Lagrangian methods may be a difficult task. The boundary conditions have therefore been given special attention in this paper. Results presented for one quasi-stationary and three non-stationary flow tests reveal a high degree of similarity, when compared to published numerical and experimental data. The difference is found to be below 5, in the case where experimental data was found tabulated. The transition from regular reflection (RR) to Mach reflection (MR) and the opposite transition from MR to RR are studied. The results are found to be in close agreement with the results obtained from various empirical and semi-empirical formulas published in the literature. A convergence test shows a convergence rate slightly steeper than linear, comparable to what is found for other numerical methods when shocks are involved.     

Investigation of the acceleration of aluminum particles behind a shock wave using instantaneous Laser Doppler Velocimetry

The acceleration of aluminum particles with a 5μm diameter in the flow field behind an incident shock wave was investigated experimentally in a 10-m long and 70 mm inner diameter shock tube. By means of instantaneous Laser Doppler Velocimetry (LDV) the velocity of the particles was observed directly. The light scattered by the moving particles is Doppler shifted and sent to the laser Doppler velocimeter. The velocimeter essentially consists of a phase-stabilized Michelson interferometer used as a sensitive spectrometer. An electro-optical circuit ensures the phase stabilization that results in a voltage signal independent of the scattered light intensity and proportional to the mean velocity of the particles at the measurement point. Because of the very short response time (1μs) of the LDV system used here, the latter gives a continuous real-time signal of the particle acceleration. To avoid particle oxidation the particles were accelerated by a high-speed nitrogen gas flow. From the measured velocity the dimensionless drag coefficient was calculated. The drag coefficient is related to the fluid dynamic force exerted by the gas on the particles. The experimental data were compared to theoretical models from the literature. A significant deviation between the model and the experimental data was observed. This deviation is supposed to be induced by the shock wave, which hits the particles and breaks them into pieces of a smaller diameter. Further experiments will be carried out in the future to check the size distribution of the particles after the shock has gone past them.      

Numerical analysis of breaking waves using the moving particle semi-implicit method

The numerical method used in this study is the moving particle semi-implicit (MPS) method, which is based on particles and their interactions. The particle number density is implicitly required to be constant to satisfy incompressibility. A semi-implicit algorithm is used for two-dimensional incompressible non-viscous flow analysis. The particles whose particle number densities are below a set point are considered as on the free surface. Grids are not necessary in any calculation steps. It is estimated that most of computation time is used in generation of the list of neighboring particles in a large problem. An algorithm to enhance the computation speed is proposed. The MPS method is applied to numerical simulation of breaking waves on slopes. Two types of breaking waves, plunging and spilling breakers, are observed in the calculation results. The breaker types are classified by using the minimum angular momentum at the wave front. The surf similarity parameter which separates the types agrees well with references. Breaking waves are also calculated with a passively moving float which is modelled by particles. Artificial friction due to the disturbed motion of particles causes errors in the flow velocity distribution which is shown in comparison with the theoretical solution of a cnoidal wave. © 1998 John Wiley & Sons, Ltd.     

Formation of a jet and vortices behind a shock wave reflected from a concave body

A jet and vortices have been observed when a plane shock wave reflects from a concave body in a shock tube. If the cavity is deep enough then two reflected shocks appear near its edges. Air, carbon tetrafluoride (CF) and dichlorodifluoromethane (CClF) were chosen as test gases. The flow was visualized with the aid of a conventional shadow technique. Pressure measurements at the body surface were also obtained. Numerical studies have been conducted using a two-dimensional inviscid model. There is a good qualitative agreement between the experimental and numerical results. Received 8 February 1996 / Accepted 30 June 1997     

The Mach reflection triple-point locus for internal and external conical diffraction of a moving shock wave

This paper examines the different behavior that occurs for the Mach reflection triple-point loci between the two fundamental axisymmetric cases, these being the external diffraction by a cone and the internal diffraction within a conically contracting channel. From equations derived in this paper using a shock dynamics approach, it has been shown that, for external diffraction over a cone, a possible solution is that the triple-point locus is a straight line which corresponds to the experimental results available, while for internal diffraction along a conically converging channel, it cannot be straight and is, in fact, a convex curve. In the latter case, a transition point is noted on the triple-point locus before which the locus is nearly straight but after which the curvature becomes marked. The second region diminishes as a proportion of the total locus with decreasing half cone angle.For the external case, a set of simple, axisymmetric equations are derived which allow a rapid estimation of the triple point locus angle and the Mach stem strength for any incident shock Mach number and cone angle combination. The equations for internal diffraction are similar and allow a quick computation of both the curved triple-point locus and the strength of the diffracting front of the shock wave. A comparison with experiment has been carried out and agreement is good.     

The method of mixed boundary condition for a kion of linear and nonlinear composite structure

This paper deals with the axisymmetrical deformation of shallow shells in large deflection which are in conjunction with linear elastic structures at the boundary: A method of mixed boundary condition for this problem is introduced, then the problem of a composite structure is transformed into a problem of a single structure and the integral equations are given. The perturbation method is used to obtain the solutions and an example of composite structure consisting of a shallow spherical and a cylindrical shell is presented.Communicated by Yeh Kai-yuan     

Thermal shielding of a sphere from the action of a strong shock wave

The time-dependent interaction of an incident shock wave with a sphere is considered in the presence of a heat supply region ahead of the body. The reflected shock configuration and the flow pattern are numerically investigated. The efficiencies of heat shields of different shapes are compared with respect to the longitudinal force acting on the sphere.     

Development of a multi‐scale ocean model by using particle Laplacian method for anisotropic mass transfer

The direct injection of CO₂ into the deep ocean is one of the feasible ways for the mitigation of the global warming, although there is a concern about its environmental impact near the injection point. To minimize its biological impact, it is necessary to make CO₂ disperse as quickly as possible, and it is said that injection with a pipe towed by a moving ship is effective for this purpose. Because the injection ship moves over a spatial scale of O(10²km), a mesoscale model is necessary to analyse the dispersion of CO₂. At the same time, since it is important to investigate high CO₂ concentration near the injection point, a small‐scale model is also required. Therefore, in this study, a numerical model was developed to analyse CO₂ dispersion in the deep ocean by using a fixed mesoscale and a moving small‐scale grid systems, the latter of which is nested and moves in the former along the trajectory of the moving ship. To overcome the artificial diffusion of mass concentration at the interface of the two different grid systems and to keep its spatial accuracy almost the same as that in the small‐scale, a particle Laplacian method was adopted and newly modified for anisotropic diffusion in the ocean. Copyright © 2010 John Wiley & Sons, Ltd.     

An analysis of the modeling of stress and strain fields in a real microstructure using a hybrid method

The distribution of stress and strain fields in a micro-structural area of a particle reinforced composite is studied by a combination of experimental and numerical method (hybrid method). With the experimental values of displacements in a micro-region as the boundary loading condition, strain and stress fields inside the micro-region are calculated by the finite element method under two different kinds of modeling, namely, as plane stress and plane strain condition. The differences between the two kinds of modeling conditions as applied to micro-structural areas are discussed. Project Supported by the National Natural Science Foundation of China (19972046) and National Overseas Study Foundation.