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.     

Solution of Riemann problem for dusty gas flow

A direct approach is used to solve the Riemann problem for a quasilinear hyperbolic system of equations governing the one dimensional unsteady planar flow of an isentropic, inviscid compressible fluid in the presence of dust particles. The elementary wave solutions of the Riemann problem, that is, shock waves, rarefaction waves and contact discontinuities are derived and their properties are discussed for a dusty gas. The generalised Riemann invariants are used to find the solution between rarefaction wave and the contact discontinuity and also inside rarefaction fan. Unlike the ordinary gasdynamic case, the solution inside the rarefaction waves in dusty gas cannot be obtained directly and explicitly; indeed, it requires an extra iteration procedure. Although the case of dusty gas is more complex than the ordinary gas dynamics case, all the parallel results for compressive waves remain identical. We also compare/contrast the nature of the solution in an ordinary gasdynamics and the dusty gas flow case.     

Generalized Riemann problem for gas dynamic combustion

The generalized Riemann problem for gas dynamic combustion in a neighborhood of the origin t > 0 in the (x, t) plane is considered. Under the modified entropy conditions, the unique solutions are constructed, which are the limits of the selfsimilar Zeldovich-von Neumann-Dring (ZND) combustion model. The results show that, for some cases, there are intrinsical differences between the structures of the perturbed Riemann solutions and the corresponding Riemann solutions. Especially, a strong detonation in the...     

Analytical solution of a problem of a collision of two hypersonic gas flows from symmetric sources

An asymptotic solution of the Euler equations that describe stationary interaction of two hypersonic gas flows from two identical spherically symmetric sources and an integral equation determining the shock wave shape are obtained with the use of a modified method of expansion of the sought functions with respect to a small parameter, which is the ratio of gas densities in the incoming flow and behind the shock wave. The solution of this equation near the axis of symmetry allows the shock wave stand-off distance from the contact plane and the radius of its curvature to be found. It is shown that the solution obtained agrees well with the known numerical solutions.     

Vanishing Viscosity Limit of the Compressible Navier–Stokes Equations for Solutions to a Riemann Problem

We study the vanishing viscosity limit of the compressible Navier–Stokes equations to the Riemann solution of the Euler equations that consists of the superposition of a shock wave and a rarefaction wave. In particular, it is shown that there exists a family of smooth solutions to the compressible Navier–Stokes equations that converges to the Riemann solution away from the initial and shock layers at a rate in terms of the viscosity and the heat conductivity coefficients. This gives the first mathematical justification of this limit for the Navier–Stokes equations to the Riemann solution that contains these two typical nonlinear hyperbolic waves.     

Contributions to numerical developments in shock waves attenuation in porous filters

The paper deals with the numerical method of the compressible gas flow through a porous filter emphasizing the treatment of the interface between a pure gaseous phase and a solid phase. An incident shock wave is initiated in the gaseous phase interacting with a porous filter inducing a transmitted and a reflected wave. To take into account the discontinuity jump in the porosity between the gaseous phase and the porous filter, an approximate Riemann solver is used to compute homogeneous non-conservative Euler equations in porous media using ideal gas state law. The discretization of this problem is based on a finite volume method where the fluxes are evaluated by a “volumes finis Roe” (VFRoe) scheme. A stationary solution is determined with a continuous variable porosity in order to test the numerical scheme. Numerical results are compared with the two-phase shock tube experiments and simulations of a shock wave attenuation and gas filtration in porous filters are presented.      

Generalized and exact solutions for oblique shock waves of real gases with application to real air

Calculation of the oblique shock wave of real gases is a difficult and time consuming problem because it involves numerical solution of a set of 10 equations, two of which (i.e., the equation of state and enthalpy function)—if available—are of a very complicated algebraic form. The present work presents a generalized method for calculating oblique shock waves of real gases, based on the Redlich-Kwong equation of state. Also described is an exact method applicable when the exact equation of state and enthalpy function of a real gas are available. Application of the generalized and the exact methods in the case of real air showed that the former is very accurate and at least twenty times faster than the latter. An additional contribution of the study is the derivation of real gas oblique shock wave equations, which are of the same algebraic form as the well known ideal gas normal shock wave relations.     

On the Riemann problem for liquid or gas-liquid media

Self-similar solutions to the Riemann problem for water with the modified Tait equation of state are presented. The methods of Smoller for gas dynamics are employed to reduce the problem to the solution of a single non-linear equation. The same methods are used for solving the Riemann problem at a gas-water interface. In both cases the method of interval bisections affords a solution technique free of problems with convergence.     

Multidimensional Shock Interaction for a Chaplygin Gas

A Chaplygin gas is an inviscid, compressible fluid in which the acoustic fields are linearly degenerate. We analyze the multidimensional shocks in such a fluid, which turn out to be sonic. Two shocks in general position interact rather simply. We investigate several two-dimensional Riemann problems and prove the existence of a unique solution. Among them is the supersonic reflection of a planar shock against a wedge; we remark that the solution cannot be a Mach reflection, contrary to what happens for other gases, and that there always exists a solution in the form of a regular reflection.     

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.     

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.     

Investigation of the influence of blowing on the flow in a hypersonic viscous shock layer near the stagnation streamline on a blunt body

In the framework of the approximation of local similarity to the Navier-Stokes equations, an investigation is made of the axisymmetric flow of homogeneous gas in a hypersonic shock layer, this including the region of transition through the shock wave. Boundary conditions, which take into account blowing of gas, are specified on the surface of the body and in the undisturbed flow. A numerical solution to the problem is obtained in a wide range of variation of the Reynolds number and the blowing parameter. Expressions are found for the dependences on the blowing parameter usually employed in boundary layer theory of the coefficients of friction and heat transfer on the surface of the body, which are divided by their values obtained for blowing parameter equal to zero. It is shown that these dependences are universal and the same as the dependences obtained from the solution of the equations of a hypersonic viscous shock layer with modified Rankin-Hugoniot relations across the shock wave and from the solution of the boundary layer equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 199–202, January–February, 1980.     

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.     

A CWENO ghost fluid method for compressible multimaterial flow

In this work a new ghost fluid method (GFM) is introduced for multimaterial compressible flow with arbitrary equation of states. In previous researches, it has been shown that accurate wave decomposition at the interface by solving a Riemann problem alleviates the shortcomings of the standard GFM in dealing with the impingement of strong waves onto the interface but these Riemann‐based GFM are not consistent with the framework of the central WENO scheme in which the emphasis is to avoid solving Riemann problems at control volume faces and enjoy the black box property (being independent of equation of state). The aim of this work is to develop a new GFM that is completely consistent with the methodology behind central schemes; that is, it enjoys a black box property. The capabilities of the proposed GFM method is shown by solving various types of multimaterial compressible flows including gas–gas, gas–water and fluid–solid interfaces interacting with strong shock waves in one and two space dimensions. Copyright © 2013 John Wiley & Sons, Ltd.     

CHAPLYGIN EQUATION IN THREE-DIMENSIONAL NON-CONSTANT ISENTROPIC FLOW-THE THEORY OF FUNCTIONS OF A COMPLEX VARIABLE UNDER DIRAC-PAULI REPRESENTATION AND ITS APPLICATION IN FLUID DYNAMICS(Ⅲ)

This work is the continuation of the discussion of Ref. [1]. In this paper we resolve the equations of isentropic gas dynamics into two problems: the three-dimensional non-constant irrotational flow (thus the isentropic flow, too), and the three-dimensional non-constant indivergent flow (i. c. the in compressible isentropic flow). We apply the theory of functions of a complex variable under Dirac-Pauli representation and the Legendre transformation, transform these equations of two problems from physical space into velocity space, and obtain two general Chaplygin equations in this paper. The general Chaplygin equation is a linear difference equation, and its general solution can be expressed at most by the hypergeometric functions. Thus we can obtain the general solution of general problems for the three-dimensional non-constant isentropic flow of gas dynamics.     

A particular solution for the chaplygin equation

In the theoretical studies of several gasdynamic problems a major role is played by the hodograph plane, where the equations in terms of velocity component variables are linear. In these studies a primary role is played by the Chaplygin equation for the stream function . Chaplygin [1] obtained a general solution for the equation of motion in the hodograph plane. Particular exact solutions of the hodograph are also known [2]: radial flow, spiral flow, etc. Below we consider a particular solution of the Chaplygin equation.     

Constitutive behaviour of anisotropic materials under shock loading

Thermodynamically and mathematically consistent constitutive equations suitable for shock wave propagation in an anisotropic material are presented in this paper. Two fundamental tensors α_ij and β_ij which represent anisotropic material properties are defined and can be considered as generalisations of the Kronecker delta symbol, which plays the main role in the theory of isotropic materials. Using two fundamental tensors α_ij and β_ij, the concept of total generalised “pressure” and pressure corresponding to the thermodynamic (equation of state) response are redefined. The equation of state represents mathematical and physical generalisation of the classical Mie–Grüneisen equation of state for isotropic material and reduces to the Mie–Grüneisen equation of state in the limit of isotropy. Based on the generalised decomposition of the stress tensor, the modified equation of state for anisotropic materials, and the modified Hill criteria, combined with the associated flow rule, a system of constitutive equations suitable for shock wave propagation is formulated. The behaviour of aluminium alloy 7010-T6 under shock loading conditions is considered. A comparison of numerical simulations with existing experimental data shows good agreement of the general pulse shape, Hugoniot Elastic Limits (HELs), and Hugoniot stress levels, and suggests that the constitutive equations are performing satisfactorily. The results are presented and discussed, and future studies are outlined.     

Quasi‐One‐Dimensional Riemann Problems and Their Role in Self‐Similar Two‐Dimensional Problems

We study two‐dimensional Riemann problems with piecewise constant data. We identify a class of two‐dimensional systems, including many standard equations of compressible flow, which are simplified by a transformation to similarity variables. For equations in this class, a two‐dimensional Riemann problem with sectorially constant data becomes a boundary‐value problem in the finite plane. For data leading to shock interactions, this problem separates into two parts: a quasi‐one‐dimensional problem in supersonic regions, and an equation of mixed type in subsonic regions. We prove a theorem on local existence of solutions of quasi‐one‐dimensional Riemann problems. For 2 × 2 systems, we generalize a theorem of Courant & Friedrichs, that any hyperbolic state adjacent to a constant state must be a simple wave. In the subsonic regions, where the governing equation is of mixed hyperbolic‐elliptic type, we show that the elliptic part is degenerate at the boundary, with a nonlinear variant of a degeneracy first described by Keldysh. (Accepted December 4, 1997)     

Numerical simulation of compressible multifluid flows using an adaptive positivity-preserving RKDG-GFM approach

The Runge-Kutta discontinuous Galerkin method together with a refined real-ghost fluid method is incorporated into an adaptive mesh refinement environment for solving compressible multifluid flows, where the level set method is used to capture the moving material interface. To ensure that the Riemann problem is exactly along the normal direction of the material interface, a simple and efficient modification is introduced into the original real-ghost fluid method for constructing the interfacial Riemann problem, and the initial conditions of the Riemann problem are obtained directly from the solution polynomials of the discontinuous Galerkin finite element space. In addition, a positivity-preserving limiter is introduced into the Runge-Kutta discontinuous Galerkin method to suppress the failure of preserving positivity of density or pressure for the problems involving strong shock wave or shock interaction with material interface. For interfacial cells in adaptive mesh refinement, the data transfer between different grid levels is achieved by using a L² projection approach along with the least squares fitting. Various numerical cases, including multifluid shock tubes, underwater explosions, and shock-induced collapse of a underwater air bubble, are computed to assess the capability of the present adaptive positivity-preserving RKDG-GFM approach, and the simulated results show that the present approach is quite robust and can provide relatively reasonable results across a wide variety of flow regimes, even for problems involving strong shock wave or shock wave impacting high acoustic impedance mismatch material interface.     

<Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> Stability Estimates for Shock Solutions of Scalar Conservation Laws Using the Relative Entropy Method

We consider scalar nonviscous conservation laws with strictly convex flux in one spatial dimension, and we investigate the behavior of bounded L ² perturbations of shock wave solutions to the Riemann problem using the relative entropy method. We show that up to a time-dependent translation of the shock, the L ² norm of a perturbed solution relative to the shock wave is bounded above by the L ² norm of the initial perturbation.