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.     

Pairwise Wave Interactions in Ideal Polytropic Gases

We consider the problem of resolving all pairwise interactions of shock waves, contact waves, and rarefaction waves in the one-dimensional flow of an ideal polytropic gas. Here, resolving an interaction means to determine the types of the three outgoing (backward, contact, and forward) waves in the Riemann problem defined by the extreme left and right states of the two incoming waves, together with possible vacuum formation. This problem has been considered by several authors and turns out to be surprisingly involved. For each type of interaction (head-on, involving a contact, or overtaking) the outcome depends on the strengths of the incoming waves. In the case of overtaking waves the type of the reflected wave also depends on the value of the adiabatic constant. Our analysis provides a complete breakdown and gives the exact outcome of each interaction.     

Solution to the Riemann problem for drift-flux model with modified Chaplygin two-phase flows

In this paper, we concern about the Riemann problem for compressible no-slip drift-flux model which represents a system of quasi-linear partial differential equations derived by averaging the mass and momentum conservation laws with modified Chaplygin two-phase flows. We obtain the exact solution of Riemann problem by elaborately analyzing characteristic fields and discuss the elementary waves namely, shock wave, rarefaction wave and contact discontinuity wave. By employing the equality of pressure and velocity across the middle characteristic field, two nonlinear algebraic equations with two unknowns as gas density ahead and behind the middle wave are formed. The Newton–Raphson method of two variables is applied to find the unknowns with a series of initial data from the literature. Finally, the exact solution for the physical quantities such as gas density, liquid density, velocity, and pressure are illustrated graphically.     

Asymptotic Stability of Combination of Viscous Contact Wave with Rarefaction Waves for One-Dimensional Compressible Navier–Stokes System

We are concerned with the large-time behavior of solutions of the Cauchy problem to the one-dimensional compressible Navier–Stokes system for ideal polytropic fluids, where the far field states are prescribed. When the corresponding Riemann problem for the compressible Euler system admits the solution consisting of contact discontinuity and rarefaction waves, it is proved that for the one-dimensional compressible Navier–Stokes system, the combination wave of a “viscous contact wave”, which corresponds to the contact discontinuity, with rarefaction waves is asymptotically stable, provided the strength of the combination wave is suitably small. This result is proved by using elementary energy methods.     

The Riemann problem for a hyperbolic model of two-phase flow in conservative form

This article is to continue the present author's work (International Journal of Computational Fluid Dynamics (2009) 23 (9), 623–641) on studying the structure of solutions of the Riemann problem for a system of three conservation laws governing two-phase flows. While existing solutions are limited and found quite recently for the Baer and Nunziato equations, this article presents the first instance of an exact solution of the Riemann problem for two-phase flow in gas–liquid mixture. To demonstrate the structure of the solution, we use a hyperbolic conservative model with mechanical equilibrium and without velocity equilibrium. The Riemann problem solution for the model equations comprises a set of elementary waves, rarefaction and discontinuous waves of various types. In particular, such a solution treats both the wave structure and the intermediate states of the two-phase gas–liquid mixture. The resulting exact Riemann solver is fully non-linear, direct and complete. On this basis then, we use locally the exact Riemann solver for the two-phase flow in gas–liquid mixture within the framework of finite volume upwind Godunov methods. In order to demonstrate the effectiveness and accuracy of the proposed solver, we consider a series of test problems selected from the open literature and compare the exact and numerical results by using upwind Godunov methods, showing excellent oscillation-free results in two-phase fluid flow problems.     

An Exact Riemann Solver for Multicomponent Turbulent Flow

The main objective of this paper is to provide some adequate way to compute the non-conservative hyperbolic system which describes a multicomponent turbulent flow. The model is written for an isentropic gas. The exact solution of the Riemann Problem (RP) associated to the hyperbolic system is exhibited. It is composed of constant states separated by rarefaction waves, or shock waves and a contact discontinuity.

The selection of the admissible part of the shock curve is obtained using an entropy criterion. This entropy is the total energy of the system. Thanks to the latter, one may compute the exact solution of the Riemann problem, assuming genuinely non linear fields contain sufficiently weak shocks.     

The structure and asymptotic behaviour of compression waves

In the study of weak solutions to nonlinear hyperbolic partial differential equations both rarefaction waves and compression waves arise. Although the behavior of rarefaction waves is known for all time, the characteristics that determine a compression wave intersect and hence the development of the wave is not easily determined. The purpose of this paper is to study compression waves. As a first step we consider the Cauchy problem for the nonlinear wave equation. We show that if the data outside some finite interval consist of constant states, then after finite time the solution involves the same states as does the solution to the Riemann problem determined by these constant states. This result is then applied to compression waves to obtain information on the shock that arises and on the steady-state solution. The region of interaction is also described. This information is obtained via a constructive procedure.     

Lifting of Disperse Particles from a Cavity Behind the Front of an Unsteady Shock Wave with a Triangular Velocity Profile

The gas flow in plane shock waves slipping along an impermeable surface with a rectangular cavity where solid disperse particles are suspended is considered numerically. The motion of the gas and particles (gas suspension) is modeled by equations of mechanics of multiphase media. Some laws of the behavior of the dusty cloud in the cavity are established for the case of wave interaction with the cavity.     

Impact of the interplanetary magnetic field on thewave flow pattern in the neighborhood of the Earth’s bow shock under sharp variations in the solar wind dynamic pressure

The impact of the interplanetary magnetic field on transformation and disintegration of the Earth’s bow shock into a system of magnetohydrodynamic (MHD) shock waves, rotational discontinuities and rarefaction waves under the action of abrupt variations in the solar wind dynamic pressure is simulated in the three-dimensional non-plane-polarized formulation within the framework of the ideal magnetohydrodynamic model using the solution of the MHD Riemann problem of breakdown of an arbitrary discontinuity. This discontinuity arises when a contact discontinuity, on which the solar wind density increases or decreases suddenly and which travels together with the solar wind, impinges on the Earth’s bow shock and propagates along its surface. The interaction pattern is constructed in the quasisteady- state formulation as a mosaic of exact solutions obtained on computer using an original MHD Riemann solver. The wave flow patterns are found for all elements of the surface of the bow shock as functions of their latitude and longitude for various jumps in the density on the contact discontinuity and characteristics parameters of the solar wind and interplanetary magnetic field at the Earth’s orbit. It is found that when the solar wind dynamic pressure increases, a fast MHD shock wave, that first penetrates into the magnetosheath, is always formed. When the solar wind dynamic pressure decreases, the influence of the interplanetary magnetic field can lead to the development of the leading fast MHD shock wave in certain zones on the surface of the Earth’s bow shock. The solution obtained can be used to interpret measurements on spacecraft in the solar wind at the libration point and in the neighborhood of the Earth’s magnetosphere.     

The Riemann problem for non-ideal isentropic compressible two phase flows

We consider the Riemann problem for a five-equation, two-pressure (5E2P) model of non-ideal isentropic compressible gas–liquid two-phase flows. This system is more complex due to the extended thermodynamics model for van der Waals gases, that is, typical real gases for gas phase and Tait׳s equation of state for liquid phase. The overall model is strictly hyperbolic and non-conservative form. We investigate the structure of Riemann problem and construct the solution for it. To construct solution of Riemann problem approximately assuming that all waves corresponding to the genuinely non-linear characteristic fields are rarefaction and then we discuss their properties. Lastly, we discuss numerical examples and study the solution influenced by the van der Waals excluded volume.     

Chapman-Jouguet waves in a dusty gas

An analytic solution is obtained of the problem of flow of a two-phase medium, representing a mixture of gas and solid or liquid particles behind plane, cylindrical, and spherical Chapman-Jouguet detonation waves. It is assumed that all the particles are the same, are chemically inert, have a true density much greater than the density of the gas, and that their volume concentration a is low. The interaction of the particles and the influence of Brownian motion on them are disregarded. The gas is assumed to be perfect. On the detonation wave, the particle parameters are assumed to be continuous, and the usual gas-dynamical relations on the detonation wave have been applied for the gas parameters because is low. Behind the detonation front, the phases interact through interphase forces and heat transfer. It has been found that the dust content of the combustible gas qualitatively changes the character of flows with Chapman-Jouguet (C-J) waves. It is shown that a plane C-J wave is an envelope of one of the acoustic families of characteristics, and not a characteristic, as occurs in a pure gas [1]. In view of this, only two solutions of the problem of flow behind a plane C-J wave are possible: one solution corresponds to a rarefaction flow and the other to a compression flow. In a pure gas such a problem has a nondenumerable set of solutions: an arbitrary Riemann rarefaction wave can adjoin the plane C-J wave. It is found that in a dusty gas there are converging cylindrical and spherical C-J waves. In a pure gas, there are no converging C-J waves [2, 3]. An expression is found for the distance r_* from the axis (center) of symmetry on which the converging cylindrical (spherical) C-J wave changes into a supercompressed detonation wave. It has been found that r_* d/₀, = 1, 2 for the cylindrical and spherical waves, respectively, d is the particle diameter, ₀ is their initial volume concentration, and the proportionality factor decreases together with d. For the detonating mixture 2H₂ + O₂ the calculations of r_* are given in a number of cases.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 111–118, July–August, 1985.The author wishes to thank V. A. Levin for his interest in the work and his useful discussion of the results.     

MHD simulation of the distribution of the gasdynamic parameters and magnetic field behind the Earth’s bow shock under sharp variations in the solar wind dynamic pressure

The distributions of the gasdynamic parameters (density, pressure, and velocity) and the magnetic field behind the Earth’s bow shock (on the outer boundary of the magnetosheath) generated under sharp variations in the solar wind dynamic pressure are found in the three-dimensional non-planepolarized formulation with allowance for the interplanetary magnetic field within the framework of the ideal magnetohydrodynamic model using the solution to the MHD Riemann problem of breakdown of an arbitrary discontinuity. Such a discontinuity which depends on the inclination of an element of the bow shock surface arises when a contact discontinuity traveling together with the solar wind and on which the solar wind density and, consequently, the dynamic pressure, increases or decreases suddenly impinges on the Earth’s bow shock and propagates along its surface initiating the development of to six waves or discontinuities (shocks). The general interaction pattern is constructed for the entire bow shock surface as a mosaic of exact solutions to the MHD Riemann problem obtained on computer using an original software (MHD Riemann solver) so that the flow pattern is a function of the angular surface coordinates (latitude and longitude). The calculations are carried out for various jumps in density on the contact discontinuity and characteristics parameters of the solar wind and interplanetary magnetic field at the Earth’s orbit. It is found that there exist horseshoe zones on the bow shock in which the increase in the density and the magnetic field strength in the fast shock waves or their reduced decrease in the fast rarefaction waves penetrating into the magnetosheath and arising as a result of sharp variation in the solar wind dynamic pressure is superposed on significant drop in the density and growth in the magnetic field strength in slow rarefaction waves. The distributions of the hydrodynamic parameters and the magnetic field can be used to interpret measurements carried out on spacecraft in the solar wind at the libration point and orbiters in the neighborhood of the Earth’s magnetosphere.     

Hyperbolicity singularities in Rarefaction Waves

For mixed-type systems of conservation laws, rarefaction waves may contain states at the boundary of the elliptic region, where two characteristic speeds coincide, and the Lax family of the wave changes. Such contiguous rarefaction waves form a single fan with a continuous profile. Different pairs of families may appear in such rarefactions, giving rise to novel Riemann solution structures. We study the structure of such rarefaction waves near regular and exceptional points of the elliptic boundary and describe their effect on Riemann solutions.     

WAVE INTERACTIONS IN SULICIU MODEL FOR DYNAMIC PHASE TRANSITIONS

Elementary waves in Suliciu model for dynamic phase transitions are obtained through traveling wave analysis.For any given initial data with two pieces of constant states,the Riemann solutions are constructed as a combination of elementary waves. When the initial profile contains three pieces of constant states,the solution may be constructed from the Riemann solutions,with each two adjacent states connected by elementary waves.A new Riemann problem forms when these two waves collide.Through the exploration of these Riemann problems,the outcome of wave interactions may be classified in a suitable parametric space.     

Over forty years of continuous research at UTIAS on nonstationary flows and shock waves

Analytical and experimental research on non-stationary shock waves, rarefaction waves and contact surfaces has been conducted continuously at UTIAS since its inception in 1948. Some unique facilities were used to study the properties of planar, cylindrical and spherical shock waves and their interactions. Investigations were also performed on shock-wave structure and boundary layers in ionizing argon, water-vapour condensation in rarefaction waves, magnetogasdynamic flows, and the regions of regular and various types of Mach reflections of oblique shock waves. Explosively-driven implosions have been employed as drivers for projectile launchers and shock tubes, and as a means of producing industrial-type diamonds from graphite, and fusion plasmas in deuterium. The effects of sonic-boom on humans, animals and structures have also formed an important part of the investigations. More recently, interest has focussed on shock waves in dusty gases, the viscous and vibrational structure of weak spherical blast waves in air, and oblique shock-wave reflections. In all of these studies instrumentation and computational methods have played a very important role. A brief survey of this work is given herein and in more detail in the relevant references.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Wave interactions and variation estimates for self-similar zero-viscosity limits in systems of conservation laws

We consider the problem of self-similar zero-viscosity limits for systems ofN conservation laws. First, we give general conditions so that the resulting boundary-value problem admits solutions. The obtained existence theory covers a large class of systems, in particular the class of symmetric hyperbolic systems. Second, we show that if the system is strictly hyperbolic and the Riemann data are sufficiently close, then the resulting family of solutions is of uniformly bounded variation and oscillation. Third, we construct solutions of the Riemann problem via self-similar zero-viscosity limits and study the structure of the emerging solution and the relation of self-similar zero-viscosity limits and shock profiles. The emerging solution consists ofN wave fans separated by constant states. Each wave fan is associated with one of the characteristic fields and consists of a rarefaction, a shock, or an alternating sequence of shocks and rarefactions so that each shock adjacent to a rarefaction on one side is a contact discontinuity on that side. At shocks, the solutions of the self-similar zero-viscosity problem have the internal structure of a traveling wave.     

A Free-Lagrange method for unsteady compressible flow: simulation of a confined cylindrical blast wave

A Free-Lagrange numerical procedure for the simulation of two-dimensional inviscid compressible flow is described in detail. The unsteady Euler equations are solved on an unstructured Lagrangian grid based on a density-weighted Voronoi mesh. The flow solver is of the Godunov type, utilising either the HLLE (2 wave) approximate Riemann solver or the more recent HLLC (3 wave) variant, each adapted to the Lagrangian frame. Within each mesh cell, conserved properties are treated as piece-wise linear, and a slope limiter of the MUSCL type is used to give non-oscillatory behaviour with nominal second order accuracy in space. The solver is first order accurate in time. Modifications to the slope limiter to minimise grid and coordinate dependent effects are described. The performances of the HLLE and HLLC solvers are compared for two test problems; a one-dimensional shock tube and a two-dimensional blast wave confined within a rigid cylinder. The blast wave is initiated by impulsive heating of a gas column whose centreline is parallel to, and one half of the cylinder radius from, the axis of the cylinder. For the shock tube problem, both solvers predict shock and expansion waves in good agreement with theory. For the HLLE solver, contact resolution is poor, especially in the blast wave problem. The HLLC solver achieves near-exact contact capture in both problems. Received May 25, 1995 / Accepted September 11, 1995     

Simple waves for two-dimensional compressible pseudo-steady Euler system

A simple wave is defined as a flow in a region whose image is a curve in the phase space. It is well known that ＂the theory of simple waves is fundamental in building up the solutions of flow problems out of elementary flow patterns＂ see Courant and Friedrichs＇s chassical book ＂Supersonic Flow and Shock Waves＂. This paper mainly concerned with the geometric construction of simple waves for the 2D pseudo-steady compressible Euler system. Based on the geometric interpretation, the expansion or compression simple wave flow construction around a pseudo-stream line with a bend part are constructed. It is a building block that appears in the global solution to four contact discontinuities Riemann problems.     

Solution of variational problems of gasdynamics of supersonic nonequilibrium flows

The problem of the optimization of the supersonic portion of a nozzle for gas flow in the case of certain nonequilibrium processes was solved in [1, 2]. The authors examined the flow scheme in which the closing Mach line of the first family arrives at the initial rarefaction wave fan. At the same time, in [3] in the solution of the analogous problem for the case of gas flow with foreign particles it was shown that it is advisable to consider also a different scheme, in which the closing characteristic arrives at the axis of symmetry outside the initial rarefaction wave fan. In the following we present results of a study of such a scheme for gas flow with nonequilibrium processes taking place. The necessary conditions which define the optimum contour are obtained and, in particular, the conditions which define the coordinate x and the magnitude of the angle at the corner points.