Role of the momentum interpolation mechanism of the Roe scheme in shock instability

The shock instability phenomenon is a well‐known problem for hypersonic flow computation by the shock‐capturing Roe scheme. The pressure checkerboard is another well‐known problem for low‐Mach‐number flow computation. The momentum interpolation method (MIM) is necessary for low‐Mach‐number flows to suppress the pressure checkerboard problem, and the pressure‐difference‐driven modification for cell face velocity can be regarded as a version of the MIM by subdividing the numerical dissipation of the Roe scheme. In this paper, MIM has been discovered through analysis and numerical tests to have the most important function in shock instability. MIM should be completely removed for nonlinear flows. However, the unexpected MIM is activated on the cell face nearly parallel to the flow for the high‐Mach‐number flows or low‐Mach‐number cells in numerical shock. Therefore, MIM should be retained for low‐Mach‐number flows and be completely removed for high‐Mach‐number flows and low‐Mach‐number cells in numerical shock. For such conditions, two coefficients are designed on the basis of the local Mach number and a shock detector. Thereafter, the improved Roe scheme is proposed. This scheme considers the requirement of MIM for incompressible and compressible flows, and is validated for good performance of numerical tests. An acceptable result can also be obtained with only the Mach number coefficient for general practical computation. Therefore, the objective of decreasing rather than increasing numerical dissipation to cure shock instability can be achieved with simple modification. Moreover, the mechanism of shock instability has been profoundly understood, in which MIM plays the most important role, although it is not the only factor. Copyright © 2016 John Wiley & Sons, Ltd.     

Unstructured adaptive remeshing finite element method for dusty shock flow

The passage of planar shocks in a dusty gas was investigated to note effects due to particle loading and initial shock Mach number. Two-phase flow equations have been added to a conservative, monotonic flow solver to allow study of compressible particle and droplet flows, which are of importance for shock propagation in two-phase flows and spray propulsion systems. The formulation developed herein employed a conservative Eulerian treatment for the gas and particle phases. The computations were performed using the finite element method-flux corrected transport (FEM-FCT) scheme, which has shown excellent predictive capability of various compressible flows which include both strong and weak shocks. The flux limiting technique was modified to provide monotonic particle velocity fields to increase the scheme's computational stability. Adaptive unstructured methodology based on adapting to high gradients of both the fluid and particle densities was used in conjunction with the conservative shock-capturing scheme to adequately resolve strong flowfield gradients. The shock attenuation of this scheme was then compared with previous experimental and numerical results and was found to yield robust predictions. Various interphase coupling terms were also considered to note their effect on the shock attenuation.This article was processed using Springer-Verlag t_ex Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Simultaneous concentration and velocity field measurements in a shock-accelerated mixing layer

A novel technique to obtain simultaneous velocity and concentration measurements is applied to the Richtmyer–Meshkov instability. After acceleration by a Mach 2.2 shock wave, the interface between the two gases develops into a turbulent mixing layer. A time-separated pair of acetone planar laser-induced fluorescence images are processed to yield concentration and, through application of the Advection-Corrected Correlation Image Velocimetry technique, velocity fields. This is the first application of this technique to shock-accelerated flows. We show that when applied to numerical simulations, this technique reproduces the velocity field to a similar quality as particle image velocimetry. When applied to the turbulent mixing layer of the experiments, information about the Reynolds number and anisotropy of the flow is obtained.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of Compressibility and Shock-Wave Interactions on Turbulent Shear Flows

Compressibility effects are present in many practical turbulent flows, ranging from shock-wave/boundary-layer interactions on the wings of aircraft operating in the transonic flight regime to supersonic and hypersonic engine intake flows. Besides shock wave interactions, compressible flows have additional dilatational effects and, due to the finite sound speed, pressure fluctuations are localized and modified relative to incompressible turbulent flows. Such changes can be highly significant, for example the growth rates of mixing layers and turbulent spots are reduced by factors of more than three at high Mach number. The present contribution contains a combination of review and original material. We first review some of the basic effects of compressibility on canonical turbulent flows and attempt to rationalise the differing effects of Mach number in different flows using a flow instability concept. We then turn our attention to shock-wave/boundary-layer interactions, reviewing recent progress for cases where strong interactions lead to separated flow zones and where a simplified spanwise-homogeneous problem is amenable to numerical simulation. This has led to improved understanding, in particular of the origin of low-frequency behaviour of the shock wave and shown how this is coupled to the separation bubble. Finally, we consider a class of problems including side walls that is becoming amenable to simulation. Direct effects of shock waves, due to their penetration into the outer part of the boundary layer, are observed, as well as indirect effects due to the high convective Mach number of the shock-induced separation zone. It is noted in particular how shock-induced turning of the detached shear layer results in strong localized damping of turbulence kinetic energy.     

HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

Study on Mach stems induced by interaction of planar shock waves on two intersecting wedges

The properties of Mach stems in hypersonic corner flow induced by Mach interaction over 3D inter-secting wedges were studied theoretically and numerically. A new method called “spatial dimension reduction” was used to analyze theoretically the location and Mach num-ber behind Mach stems. By using this approach, the problem of 3D steady shock/shock interaction over 3D intersecting wedges was transformed into a 2D moving one on cross sec-tions, which can be solved by shock-polar theory and shock dynamics theory. The properties of Mach interaction over 3D intersecting wedges can be analyzed with the new method, including pressure, temperature, density in the vicinity of triple points, location, and Mach number behind Mach stems. Theoretical results were compared with numerical results, and good agreement was obtained. Also, the influence of Mach number and wedge angle on the properties of a 3D Mach stem was studied.     

Application of particle image velocimetry and the background-oriented schlieren technique in the high-enthalpy shock tunnel G?ttingen

The applicability of the particle image velocimetry (PIV) and the background-oriented schlieren (BOS) techniques in the high-enthalpy shock tunnel G?ttingen of the German Aerospace Center, DLR is demonstrated. As a part of this feasibility study two different experiments are performed. The velocity field past a wedge in a Mach 6 flow at a total specific enthalpy of 1.5 MJ/kg is determined by means of PIV and the results are compared to numerical predictions. The BOS technique is applied to investigate the density field in the shock layer of a sphere at 12 and 22 MJ/kg total specific enthalpies. Using a ray tracer method, the BOS results are compared to the data obtained by corresponding numerical computations.     

Analysis of unsteady wave processes in a rotating channel

The impact of passage rotation on the gasdynamic wave processes is analyzed through a numerical simulation of ideal shock-tube flow in a closed rotating-channel containing a gas in an initial state of homentropic solid-body rotation. Relevant parameters of the problem such as wheel Mach number, hub-to-tip radius ratio, length-to-tip radius ratio, diaphragm temperature ratio, and diaphragm pressure ratio are varied. It is shown that for a fixed geometry and initial conditions, the contact interface acquires a distorted three-dimensional time-dependent orientation at non-zero wheel Mach numbers. At a fixed wheel Mach number, the level of distortion depends primarily on the density ratio across the interface and also the hub-to-tip radius ratio. The nature of the rarefaction and shock wave propagation is one-dimensional, although the acoustic waves are diffracted due to the radially varying propagation speed. Under conditions of initially homentropic solid-body rotation, a degree of similarity exists between rotating and stationary shock-tube flows. This similarity is exploited to arrive at an approximate analytical solution to the Riemann problem in a rotating shock-tube.     

The formation of a secondary shock wave behind a shock wave diffracting at a convex corner

This paper deals with the formation of a secondary shock wave behind the shock wave diffracting at a two-dimensional convex corner for incident shock Mach numbers ranging from 1.03 to 1.74 in air. Experiments were carried out using a 60 mm 150 mm shock tube equipped with holographic interferometry. The threshold incident shock wave Mach number () at which a secondary shock wave appeared was found to be = 1.32 at an 81° corner and = 1.33 at a 120° corner. These secondary shock waves are formed due to the existence of a locally supersonic flow behind the diffracting shock wave. Behind the diffracting shock wave, the subsonic flow is accelerated and eventually becomes locally supersonic. A simple unsteady flow analysis revealed that for gases with specific heats ratio the threshold shock wave Mach number was = 1.346. When the value of is less than this, the vortex is formed at the corner without any discontinuous waves accompanying above the slip line. The viscosity was found to be less effective on the threshold of the secondary shock wave, although it attenuated the pressure jump at the secondary shock wave. This is well understood by the consideration of the effect of the wall friction in one-dimensional duct flows. In order to interpret the experimental results a numerical simulation using a shock adaptive unstructured grid Eulerian solver was also carried out. Received 1 May 1996 / Accepted 12 September 1996     

Submerged shocks in conical gas flows in the presence of a mach shock-wave configuration

The results of a numerical study of a new type of singularities in the Mach shock-wave structure realized in supersonic nonsymmetric conical flows over V-wings with a bow shock attached to the leading edges are presented. Within the framework of the ideal gas model we study the changes in the shock system on transition, with increase in the sweep angle, from the region of nonsymmetric Mach interaction of the shocks attached to the leading edges of the wing to the region of special flow patterns, where on the windward cantilever surface a rarefaction flow is realized rather than a flow with an internal shock. It is shown, in particular, that in the region with special wing flow patterns a Mach system of shocks with a submerged shock proceeding from the branch point above the windward cantilever may exist.     

Numerical simulation of fully-developed compressible flows over wavy surfaces

Rough surfaces are common on high-speed vehicles, for example on heat shields, but compressibility is not usually taken into account in the flow modelling other than through the mean density. In the present study, supersonic fully-developed turbulent rough wall channel flows are simulated using direct numerical simulation to investigate whether strong compressibility effects significantly alter the mean flow and turbulence properties across the channel. The simulations were run for three different Mach numbers M = 0.3, 1.5 and 3.0 over a range of wall amplitude-to-wavelength ratios from 0.01 to 0.08, corresponding to transitionally and fully rough cases respectively. The velocity deficit values are found to decrease with increasing Mach number. It is also found that at Mach 3.0 significant differences occur in the mean flow and turbulence statistics throughout the channel and not just in a roughness sublayer. These differences are found to be due to the presence of strong shock waves created by the peaks of the roughness elements.     

Investigation on onset of shock-induced separation

A great number of experimental data indicating shock wave/boundary layer interactions in internal or external supersonic flows were reviewed to make clear the mechanism of the interaction and to decide the onset of shock-induced separation. The interesting conclusions were obtained for the considerably wide range of flow geometries that the onset of separation is independent of the flow geometries and the boundary layer Reynolds number. It is found that the pressure rise necessary to separate the boundary layer in supersonic external flows could be applied to such internal flows as overexpanded nozzles or diffusers. This is due to the fact that the separation phenomenon caused by shock wave/boundary layer interactions is processed through a supersonic deceleration. The shock-induced separation in almost all of interacting flow fields is governed by the concept of free interaction, and the onset of shock-induced separation is only a function of the Mach number just upstream of shock wave. However, physical scales of the produced separation are not independent of the downstream flow fields.     

Existence criteria for reflection configurations in shock–vortex interactions

The length scale criteria is widely accepted as an explanation for transition and hence existence of different shock wave reflection configurations in pseudo-steady flows. However, there has not been any attempt to validate this criteria using information obtained from a time-dependent numerical simulation. A high resolution time-dependent numerical simulation in pseudo-steady flow is carried out in the present work. Time-dependent numerical data is used to calculate flow features in a laboratory frame of reference to verify validity of the length scale criteria for existence of different shock wave reflection configurations in pseudo-steady flow. This analysis is then extended to the study of unsteady shock wave reflection configurations in shock–vortex interactions. It is shown that the existence of regular reflection (RR) and Mach reflection (MR) configurations in an unsteady flowfield resulting from shock–vortex interactions can also be explained locally based on limiting conditions similar to that prescribed by the length scale criteria for pseudo-steady flow.

---

     

An aeroacoustic hybrid approach for non‐isothermal flows at low Mach number

This paper presents an aeroacoustic hybrid technique for the study of non‐isothermal flows at low Mach number. The flow dynamics and the acoustic production and propagation are computed separately. The fully compressible Navier–Stokes equations are modified through an expansion of the physical quantities using a low Mach number approximation. Compressibility effects are thus removed in the CFD while inhomogeneities of the flow related to heat transfer are preserved. One advantage is a reduction of the time step constraint. Another advantage is that the Mach number does not appear explicitly and a simple rescaling allows a study over a relatively wide band of subsonic Mach number flows with a single dynamic simulation. Compatible acoustic source terms for LEE based propagation have been defined and the procedure is implemented in the case of a temporal mixing layer. Compressible simulations for Mach numbers of 0.2, 0.3 and 0.4 are compared with the numerical results obtained using the proposed method. Very good agreement is obtained even at relatively high subsonic Mach number demonstrating the efficiency of the proposed technique. Copyright © 2004 John Wiley & Sons, Ltd.     

Gas-kinetic description of shock wave structures by solving Boltzmann model equation

On the basis of the mesoscopic theory of Boltzmann-type velocity distribution function, the modified Boltzmann model equation describing the one-dimensional gas flows from various flow regimes is presented by incorporating the molecular interaction models relating to the viscosity and diffusion cross-sections, density, temperature and the dependent exponent of viscosity into the molecular collision frequency. The gas-kinetic numerical method for directly solving the molecular velocity distribution function is studied by introducing the reduced distribution functions and the discrete velocity ordinate method, in which the unsteady time-splitting method and the NND finite difference scheme are applied. To study the inner flows of non-equilibrium shock wave structures, the one-dimensional unsteady shock-tube problems with various Knudsen numbers and the steady shock wave problems at different Mach numbers are numerically simulated. The computed results are found to give good agreement with the theoretical, DSMC and experimental results. The computing practice has confirmed the good precision and reliability of the gas-kinetic numerical algorithm in solving the highly nonequilibrium shock wave disturbances from various flow regimes.     

Numerical investigation of axisymmetric shock wave focusing over paraboloidal reflectors

The problem of a plane shock wave incident to a paraboloidal reflector is numerically investigated. The numerical solver used is developed by an improved, implicit, upwind total variation diminishing scheme in a finite-volume approach. The real-gas effect is taken into account if high temperature occurs. The solver is validated on four test problems. The complicated flow fields of axisymmetric shock wave focusing for different-depth reflectors at various incident shock Mach numbers are studied. An interesting result of a maximum pressure happening at the reflector center is found. This is due to the occurrence of an implosion phenomenon. A maximum temperature might occur at the reflector center or at other locations, depending on the incident shock Mach number and the reflector depth. Moreover, vortical flows induced by shock wave focusing and their formation mechanism are explored. It was found that the vortices near the reflector are caused by a ring-shaped shock/slipline interaction. Owing to the slipline on the symmetry axis, a jet flow is induced, resulting in the formation of vortices near the symmetry axis. Received 13 January 1998 / Accepted 10 November 1998     

Aerodynamic focusing of inertial particles in the shock-wave intersection region

The flow structure in the region of intersection of two stationary plane shock waves in a dusty gas is investigated. The regular symmetric and asymmetric shock interaction and the symmetric Mach interaction are considered. For a small particle mass concentration, on the basis of numerical calculations using the full Lagrangian approach it is shown that behind the shock wave intersection point a long thin region is formed, in which the particle trajectories intersect and the particle concentration sharply increases. A parametric study of the particle concentration distributions in this region is performed and the range of governing parameters on which the particle focusing effect is maximal is found.     

Low-diffusion preconditioning scheme for numerical simulation of low-speed flows past airfoil

The preconditioning technique can address the stiffness of a low Mach number flow, while its stability is poor. Based on the conventional preconditioning method of Roe's scheme, a low-diffusion scheme is proposed. An adjustable parameter is introduced to control numerical dissipation, especially over the dissipation in the boundary layer and extremely in a low speed region. Numerical simulations of the low Mach number and low Reynolds number flows past a cylinder and the low Mach number and high Reynolds number flows past NACA0012 and S809 airfoils are performed to validate the new scheme. Results of the three tests well agree with experimental data, showing the applicability of the proposed scheme to low Mach number flow simulations.     

Two-phase reactive particle flow through normal shock waves

An analysis is presented for the steady one-dimensional flow behind a normal shock wave of a compressible gas containing small spherical particles of solid propellant. The solids mass fraction is assumed large enough to require that the void volume fraction be retained as a variable in the governing conservation equations. The particles are ignited by the shocked air and by viscous interaction. Propellant gases are then generated which depend on the instantaneous size of the particles and on the linear burning rate. The latter is assumed dependent upon the local pressure and the particle temperature. These calculations are of interest because of the potential hazards of such particle flows, in that extreme pressures are predicted within the relaxation zone, pressures even greater than those calculated for the final equilibrium conditions. The results stress the importance of the Mach number of the normal (strong shock) and the energy content of the propellant (J/kg).