The generalized Riemann problem for a scalar Chapman–Jouguet combustion model

The generalized Riemann problem for a scalar Chapman–Jouguet combustion model in a neighborhood of the origin (t > 0) on the (x, t) plane is studied. Under the entropy conditions, we obtain the solutions constructively. It is found that, for some cases, the perturbed Riemann solutions are essentially different from the corresponding Riemann solutions. The perturbation may transform a combustion wave CJDT into SDT in the neighborhood of the origin. Especially, it can be observed that burning happens although the corresponding Riemann solution doesn’t contain combustion waves, which exhibits the instability for unburnt states. This work is supported by NSFC 10671120     

Perturbed Riemann Problem for a Scalar Chapman-Jouguet Combustion Model

The author considers the perturbed Riemann problem for a scalar ChapmanJouguet combustion model which comes from Majda’s model with a modified, bump-type ignition function proposed in the results of Lyng and Zumbrun in 2004. Under the entropy conditions, the unique solution in a neighborhood of the origin on the(x, t) plane(t > 0) is obtained. It is found that, for some cases, the perturbed Riemann solutions are essentially different from the corresponding Riemann solutions. The perturbation may transform a strong detonation into a weak deflagration in the neighborhood of the origin. Especially, it can be observed that burning happens although the corresponding Riemann solution does not contain combustion wave, which exhibits the instability for the unburnt state.     

THE GENERALIZED RIEMANN PROBLEM FOR A SCALAR NONCONVEX COMBUSTION MODEL-THE PERTURBATION ON INITIAL BINDING ENERGY

In this article,we study the generalized Riemann problem for a scalar nonconvex Chapman-Jouguet combustion model in a neighborhood of the origin (t ＞ 0) on the (x,t) plane.We focus our attention to the...     

Interaction of oblique deflagration and shock for two-dimensional steady adiabatic combustion system

The interaction of an oblique deflagration wave and an oblique shock wave for two-dimensional steady adiabatic combustion system is analyzed. Using the shock wave polar and combustion wave polar, we exhibit the construction of the solutions. It is found that the deflagration remains if the shock is weak. However, the shock transforms the deflagration into a detonation(DDT) if it is strong or stops the deflagration if it is proper.     

The scalar Zeldovich–von Neumann–Döring combustion model (II) interactions of shock and deflagration

In this paper, the interactions of shock and deflagration for the scalar Zeldovich–von Neumann–Döring combustion model are considered. The solutions of the problem are constructed by analysing characteristics in the reaction zone. Some burning gas in the zone will be extinguished at a finite time in some cases. By studying the limits of the solutions as the reaction rate goes to infinity, we obtain that the limits are the solutions of the corresponding initial value problem for the scalar Chapman–Jouguet combustion model.     

The scalar Zeldovich–von Neumann–Döring combustion model (I) interactions of shock and detonation

In this paper, the interactions of shock and detonation for the scalar Zeldovich–von Neumann–Do¨ring combustion model are considered. The solutions of the problem are constructed by analyzing characteristics in the reaction zone. The shock speeds up the reaction front in some cases. By studying the limits of the solutions as the reaction rate goes to infinity, we obtain that the limits are the solutions of the corresponding initial value problem for the scalar Chapman–Jouguet combustion model.     

Sonic and kinetic phase transitions with applications to Chapman–Jouguet deflagrations

We consider an n × n system of hyperbolic conservation laws and focus on the case of strongly underdetermined sonic phase boundaries. We propose a Riemann solver that singles out solutions uniquely. This Riemann solver has two features: it selects phase boundaries by means of an exterior function and it allows compound waves. Then we prove the global existence of weak solutions to the Cauchy problem. Applications to Chapman–Jouguet deflagrations are given. Copyright © 2004 John Wiley & Sons, Ltd.     

Riemann problem for Chaplygin gas with combustion

We study the Riemann problem of the Chapman–Jouguet model for an ideal combustible Chaplygin gas. By analyzing the wave curves in the phase plane, we obtain constructively the unique solution of the Riemann problem under the global entropy conditions.     

Shock formation in a chemotaxis model

In this paper, we establish the existence of shock solutions for a simplified version of the Othmer–Stevens chemotaxis model (SIAM J. Appl. Math. 1997; 57 :1044–1081). The existence of these shock solutions was suggested by Levine and Sleeman (SIAM J. Appl. Math. 1997; 57 :683–730). Here, we consider the general Riemann problem and derive the shock curves in parameterized forms. By studying the travelling wave solutions, we examine the shock structure for the chemotaxis model and prove that the travelling wave speed is identical to the shock speed. Moreover, we explicitly derive an entropy–entropy flux pair to prove the uniqueness of the weak shock solutions. Some discussion is given for further study. Copyright © 2007 John Wiley & Sons, Ltd.     

UNIFORM FORMULA FOR THE RIEMANN SOLUTIONS OF A SCALAR COMBUSTION MODEL

In this paper, we construct a uniform formula for the Riemann solutions of the simplified Chapman-Jouguet model. Firstly, we define a new functional, and then, we obtain that the Riemann solutions can be expressed by the maximum value point of this functional,while Riemann solutions may contain some of strong detonation waves, Chapman-Jouguet detonation waves and contact discontinuities. Finally, Chapman-Jouguet deflagration waves are also discussed.     

The Interaction of Shocks with Dispersive Waves I. Weak Coupling Limit

We introduce and analyze a model for the interaction of shocks with a dispersive wave envelope. The model mimicks the Zakharov system from weak plasma turbulence theory but replaces the linear wave equation in that system by a nonlinear wave equation allowing the formation of shocks. This paper considers a weak coupling in which the nonlinear wave evolves independently but appears as the potential in the time-dependent Schrodinger equation governing the dispersive wave. We first solve the Riemann problem for the system by constructing solutions to the Schrodinger equation that are steady in a frame of reference moving with the shock. Then we add a viscous diffusion term to the shock equation and by explicitly constructing asymptotic expansions in the (small) diffusion coefficient, we show that these solutions are zero diffusion limits of the regularized problem. The expansions are unusual in that it is necessary to keep track of exponentially small terms to obtain algebraically small terms. The expansions are compared to numerical solutions. We then construct a family of time-dependent solutions in the case that the initial data for the nonlinear wave equation evolves to a shock as t → t* < ∞. We prove that the shock formation drives a finite time blow-up in the phase gradient of the dispersive wave. While the shock develops algebraically in time, the phase gradient blows up logarithmically in time. We construct several explicit time-dependent solutions to the system, including ones that: (a) evolve to the steady states previously constructed, (b) evolve to steady states with phase discontinuities (which we call phase kinked steady states), (c) do not evolve to steady states.     

Convergence rates toward the travelling waves for a model system of the radiating gas

The present paper is concerned with an asymptotics of a solution to the model system of radiating gas. The previous researches have shown that the solution converges to a travelling wave with a rate (1 + t)^?1/4 as time t tends to infinity provided that an initial data is given by a small perturbation from the travelling wave in the suitable Sobolev space and the perturbation is integrable. In this paper, we make more elaborate analysis under suitable assumptions on initial data in order to obtain shaper convergence rates than previous researches. The first result is that if the initial data decays at the spatial asymptotic point with a certain algebraic rate, then this rate reflects the time asymptotic convergence rate. Precisely, this convergence rate is completely same as the spatial convergence rate of the initial perturbation. The second result is that if the initial data is given by the Riemann data, an admissible weak solution, which has a discontinuity, converges to the travelling wave exponentially fast. Both of two results are proved by obtaining decay estimates in time through energy methods with suitably chosen weight functions. Copyright © 2006 John Wiley & Sons, Ltd.     

Wave interactions and stability of the Riemann solutions for the chromatography equations

We prove that the Riemann solutions are stable for the chromatography system under the local small perturbations of the Riemann initial data. The proof is based on the detailed analysis of the wave interactions by applying the method of characteristic analysis. It is noteworthy that both the propagation directions of the shock wave S and rarefaction wave R are unchanged when they interact with the contact discontinuity J. Moreover, the global structures and large time asymptotic behaviors of the perturbed Riemann solutions are constructed and analyzed case by case.     

Properties of Detonation Waves

The phenomenon of unphysical wave propagation speeds sometimes occurs in numerical computations of detonation waves on coarse grids. The strong detonation wave splits into two parts, a weak detonation which travels with the speed of one cell per time step and an ordinary shock wave. We analyse a simplified set of equations and look for travelling wave solutions. It is shown that the solution depends on the dimensionless number Kr = μK/Qρ₁. Here μ is the viscosity, K is the rate of reaction, Q is the heat release available in the process and ρ₁ is the density at the unburnt state. It is shown that the density peak of the travelling wave depends on Kr and also, that if Kr is sufficiently large there is no travelling wave solution. The erroneous behaviour above is explained as an effect of the artificial viscosity necessarily inherent in the numerical methods when coarse grids are used. To prevent this unphysical behaviour we suggest the use of an ‘artificial rate of reaction’ such that the actual value of Kr used in the numerical method retains its correct physical value.     

扰动色谱方程黎曼解的极限-delta激波的行成和转换

本文主要讨论扰动色谱方程delta激波解的行成和转换,并讨论上述方程的黎曼问题.当扰动参数趋于零时,通过研究黎曼解的极限,我们可以观察到如下两个重要现象:激波和接触间断重合行成delta激波,一类激波(一个变量含有delta函数).     

Global solutions with shock waves to the generalized Riemann problem for a class of quasilinear hyperbolic systems of balance laws II

This work is a continuation of our previous work. In the present paper, we study the existence and uniqueness of global piecewise C¹ solutions with shock waves to the generalized Riemann problem for general quasilinear hyperbolic systems of conservation laws with linear damping in the presence of a boundary. It is shown that the generalized Riemann problem for general quasilinear hyperbolic systems of conservation laws with linear damping with nonlinear boundary conditions in the half space {(t, x) | t ≥ 0, x ≥ 0} admits a unique global piecewise C¹ solution u = u (t, x) containing only shock waves with small amplitude and this solution possesses a global structure similar to that of a self‐similar solution u = U (x /t) of the corresponding homogeneous Riemann problem, if each characteristic field with positive velocity is genuinely nonlinear and the corresponding homogeneous Riemann problem has only shock waves but no rarefaction waves and contact discontinuities. This result is also applied to shock reflection for the flow equations of a model class of fluids with viscosity induced by fading memory. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

气体动力学燃烧模型的广义Riemann问题

考察了在(x,t)平面上原点(t＞0)的邻域内气体动力学燃烧模型的广义Riemann问题．在改进的熵条件下构造了此问题的唯一解．它们是自相似ZND燃烧模型的极限．发现对某些情形,广义Riemann问题的解与相应的Riemann问题的解有本质的不同．特别地,扰动会使得相应Riemann问题的强爆轰波转化为由预压激波点燃的弱爆燃波．在一些情形,尽管相应的Riemann解中不含燃烧波,扰动后燃烧波会出现．这反映了未燃气体的不稳定性．     

Shock initiation of explosives: the idealized condensed-phase model

Many current models of condensed-phase explosives employ reactionrate law models where the form of the rate has a power-law dependenceon pressure (i.e. proportional to pⁿ where n is an adjustableparameter). Here, shock-induced ignition is investigated usinga simple model of this form. In particular, the solutions arecontrasted with those from Arrhenius rate law models as studiedpreviously. A large n asymptotic analysis is first performed,which shows that in this limit the evolution begins with aninduction stage, followed by a sequence of pressure runaways,resulting in a forward propagating, decelerating, shocklesssupersonic reaction wave (a weak detonation). The theory predictssecondary shock and super-detonation formation once the weakdetonation reaches the Chapman–Jouguet speed. However,it is found that secondary shock formation does not occur untilthe weak detonation has reached a point close to the initiatingshock, whereas for Arrhenius rate laws the shock forms closerto the piston. Numerical simulations are then conducted forO(1) values of n, and it is shown that the idealized condensed-phasemodel can qualitatively describe a wide range of experimentallyobserved behaviours, from growth mainly at the shock, to smoothgrowth of a pressure pulse behind the shock, to cases wherea secondary shock and possibly a super-detonation form. Thenumerics are used to reveal the different evolutionary mechanismsfor each of these cases. However, the evolution is found tobe sensitive to n, with the whole range of behaviours coveredby varying n from about 3 to 5. The simulations also confirmthe predictions of the theory that pressure-dependent rate lawsare unable to describe homogeneous explosive scenarios wherea super-detonation forms very close to the point of initialrunaway.     

The limits of Riemann solutions to the simplified pressureless Euler system with flux approximation

The formation of vacuum state and delta shock wave in the solutions to the Riemann problem for the simplified pressureless Euler system is considered under the linear approximations of flux functions. The method is to perturb the non‐strictly hyperbolic system into a nearby strictly hyperbolic system by introducing appropriately the linear approximations of flux functions. The solutions to the Riemann problem for the approximated system can be constructed explicitly and then the formation of vacuum state and delta shock wave can be observed by taking the perturbation parameter tend to zero in the solutions.     

Finite difference method for a combustion model

We study a projection and upwind finite difference scheme for a combustion model problem. Convergence to weak solutions is proved under the Courant-Friedrichs-Lewy condition. More assumptions are given on the ignition temperature; then convergence to strong detonation wave solutions or to weak detonation wave solutions is proved.