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.     

Stability of reflection and refraction of shocks on interface

In this paper we study the stability of the nonlinear wave structure caused by the attack of an incident shock on an interface of two different kinds of media. The attack will produce a reflected wave and a refracted wave, and also let the interface deflected. In this paper we will mainly study the case, when the reflected wave is a shock, and the flow between the reflected wave and the refracted shock is relatively subsonic. Our result indicates that the wave structure and the flow field for the reflection-refraction problem in this case is conditionally stable.To describe the motion of the fluid we use the inviscid Euler system as the mathematical model. The reflection-refraction problem can be reduced to a free boundary value problem, where the unknown reflected shock and refracted shock are free boundaries, and the deflected interface is also to be determined. In the proof of the existence and the stability of the corresponding wave structure we apply the Lagrange transformation to fix the interface and the decoupling technique to decouple the elliptic-hyperbolic composite system in its principal part. Meanwhile, some efficient weighted Sobolev estimates are established to derive the existence for corresponding nonlinear problems.     

ON INTERACTION OF SHOCK AND SOUND WAVE （I）

This paper studies the interaction of shock and gradient wave (sound wave) of solutions to the system of inviscid isentropic gas dynamics as a model for the corresponding problems for nonlinear hyperbolic systems. The problem can be reduced to a boundary value problem in a wedged dormain, By using the method of constructing asymptotic solutions and Newton‘siteration process it is proved that if a weak shock hits a gradient wave, then the grandient wave will split into two gradient waves, while the shock continuses propagating. In this paper the author reduces the problem to a standard form and constructs asymptotic solution of the problem. The existence of the genuine solution will he given in the following paper.     

The Interaction of Shocks with Dispersive Waves II. Incompressible-Integrable Limit

This is the second in a two-part series of articles in which we analyze a system similar in structure to the well-known Zakharov equations from weak plasma turbulence theory, but with a nonlinear conservation equation allowing finite time shock formation. In this article we analyze the incompressible limit in which the shock speed is large compared to the underlying group velocity of the dispersive wave (a situation typically encountered in applications). After presenting some exact solutions of the full system, a multiscale perturbation method is used to resolve several basic wave interactions. The analysis breaks down into two categories: the nonlinear limit and the linear limit, corresponding to the form of the equations when the group velocity to shock speed ratio, denoted by ε, is zero. The former case is an integrable limit in which the model reduces to the cubic nonlinear Schrödinger equation governing the dispersive wave envelope. We focus on the interaction of a “fast” shock wave and a single hump soliton. In the latter case, the ε=0 problem reduces to the linear Schrödinger equation, and the focus is on a fast shock interacting with a dispersive wave whose amplitude is cusped and exponentially decaying. To motivate the time scales and structure of the shock-dispersive wave interactions at lowest orders, we first analyze a simpler system of ordinary differential equations structurally similar to the original system. Then we return to the fully coupled partial differential equations and develop a multiscale asymptotic method to derive the effective leading-order shock equations and the leading-order modulation equations governing the phase and amplitude of the dispersive wave envelope. The leading-order interaction equations admit a fairly complete analysis based on characteristic methods. Conditions are derived in which: (a) the shock passes through the soliton, (b) the shock is completely blocked by the soliton, or (c) the shock reverses direction. In the linear limit, a phenomenon is described in which the dispersive wave induces the formation of a second, transient shock front in the rapidly moving hyperbolic wave. In all cases, we can characterize the long-time dynamics of the shock. The influence of the shock on the dispersive wave is manifested, to leading order, in the generalized frequency of the dispersive wave: the fast-time part of the frequency is the shock wave itself. Hence, the frequency undergoes a sudden jump across the shock layer.In the last section, a sequence of numerical experiments depicting some of the interesting interactions predicted by the analysis is performed on the leading-order shock equations.     

ON INTERACTION OF SHOCK AND SOUND WAVE (Ⅰ)

1.IntroductionRecentlythestudyofdiscontinuoussolutionforthesystemofconservationlawsinhigherdimensionalspacehasbeenconsiderablydeveloped.In[1,W12]thelocalekistenceof8olutionforsuclisystemwithdiscontinuityinvolvingsingleshock,rarefactionwaveorsoundwav(gradientwave)hasbeenestablished.In[2Jand[14]theproblemsoninteractionoftwoshocksorinteractionofweaksingularitiesarealsoconsidered.Itisnaturaltoaskwhatabouttheresultwhenashockisinteractedbyawavewitliweakersillgularities,particularly,forthenbynsystem…     

The singular solutions to a nonsymmetric system of Keyfitz-Kranzer type with initial data of Riemann type

The solutions to the Riemann problem for a nonsymmetric system of Keyfitz-Kranzer type are constructed explicitly when the initial data are located in the quarter phase plane. In particular, some singular hyperbolic waves are discovered when one of the Riemann initial data is located on the boundary of the quarter phase plane, such as the delta shock wave and some composite waves in which the contact discontinuity coincides with the shock wave or the wave back of rarefaction wave. The double Riemann problem for this system with three piecewise constant states is also considered when the delta shock wave is involved. Furthermore, the global solutions to the double Riemann problem are constructed through studying the interaction between the delta shock wave and the other elementary waves by using the method of characteristics. Some interesting nonlinear phenomena are discovered during the process of constructing solutions; for example, a delta shock wave is decomposed into a delta contact discontinuity and a shock wave.     

星系螺旋结构的气体激波理论

本文用星际气体自引力星系激波来解释星系的螺旋结构、恒星的扰动引力场并非必要条件.我们首先证明,即使扰动引力场为零,也可以存在局部的星系激波解.这种解要求|ω_η0|>α,而且只要气体的密度反差比较大,就只能用激波解来解释螺旋结构.用叠代的方法求出了星际气体的自引力激波宏图.对一种特定的扰动引力场模拟气体自引力,可以在速度平面上定性分析激波解的特性.初始原星系盘中的物质分布不均匀性,通过缠卷过程、不稳定性增长和波动叠加.可以发展成星系激波宏图.这样,对星系激波的起源,演化和维持给出一个完整的图象.利用这个图象,可以解释星系螺旋结构的大量观测结果和分类特性.     

激波在异种气体中传播及诱导的剪切混合研究

利用二阶迎风TVD格式求解多组分,层流全N－S方程,针对直通道和突扩直通道,研究了马赫数为2和4的激波在H2和空气界面上的传播及诱导的燃料剪切混合,计算结果表明：（1）直通道中,剪切层中的激波阵面要发生畸变,存在对混合起主要作用的卷吸涡,激波马赫数不同,卷吸涡结构和横向混合的尺寸也不同,激波马赫数低,剪切混合效果好,（2）在突扩直通道中,马赫数为2和4的激波在H2中产生不同强度激波,在剪切层中都能产生顺时针,尺度较大的卷吸涡,后台阶增强了剪切层的混合。     

有限变形弹性杆中的几何非线性波

利用有限变形理论的Lagrange描述,借助非保守系统的Hamilton型变分原理,导出了描述弹性杆中几何非线性波的波动方程．为了使非线性波动方程有稳定的行波解,计及了粘性效应引入的耗散和横向惯性效应导致的几何弥散．运用多重尺度法将非线性波动方程简化为KdV-Bergers方程,这个方程在相平面上对应着异宿鞍-焦轨道,其解为振荡孤波解．如果略去粘性效应或横向惯性,方程将分别退化为KdV方程或Bergers方程,由此得到孤波解或冲击波解,它们在相平面上对应着同宿轨道或异宿轨道．     

Continuum models in problems of the hypersonic flow of a rarefied gas around blunt bodiest

All possible continuum (hydrodynamic) models in the case of two-dimensional problems of supersonic and hypersonic flows around blunt bodies in the two-layer model (a viscous shock layer and shock-wave structure) over the whole range of Reynolds numbers, Re, from low values (free molecular and transitional flow conditions) up to high values (flow conditions with a thin leading shock wave, a boundary layer and an external inviscid flow in the shock layer) are obtained from the Navier-Stokes equations using an asymptotic analysis. In the case of low Reynolds numbers, the shock layer is considered but the structure of the shock wave is ignored. Together with the well-known models (a boundary layer, a viscous shock layer, a thin viscous shock layer, parabolized Navier-Stokes equations (the single-layer model) for high, moderate and low Re numbers, respectively), a new hydrodynamic model, which follows from the Navier-Stokes equations and reduces to the solution of the simplified (“local”) Stokes equations in a shock layer with vanishing inertial and pressure forces and boundary conditions on the unspecified free boundary (the shock wave) is found at Reynolds numbers, and a density ratio, k, up to and immediately after the leading shock wave, which tend to zero subject to the condition that (k/Re)^1/2 → 0. Unlike in all the models which have been mentioned above, the solution of the problem of the flow around a body in this model gives the free molecular limit for the coefficients of friction, heat transfer and pressure. In particular, the Newtonian limit for the drag is thereby rigorously obtained from the Navier-Stokes equations. At the same time, the Knudsen number, which is governed by the thickness of the shock layer, which vanishes in this model, tends to zero, that is, the conditions for a continuum treatment are satisfied. The structure of the shock wave can be determined both using continuum as well as kinetic models after obtaining the solution in the viscous shock layer for the weak physicochemical processes in the shock wave structure itself. Otherwise, the problem of the shock wave structure and the equations of the viscous shock layer must be jointly solved. The equations for all the continuum models are written in Dorodnitsyn--Lees boundary layer variables, which enables one, prior to solving the problem, to obtain an approximate estimate of second-order effects in boundary-layer theory as a function of Re and the parameter k and to represent all the aerodynamic and thermal characteristic; in the form of a single dependence on Re over the whole range of its variation from zero to infinity.

An efficient numerical method of global iterations, previously developed for solving viscous shock-layer equations, can be used to solve problems of supersonic and hypersonic flows around the windward side of blunt bodies using a single hydrodynamic model of a viscous shock layer for all Re numbers, subject to the condition that the limit (k/Re)^1/2 → 0 is satisfied in the case of small Re numbers. An aerodynamic and thermal calculation using different hydrodynamic models, corresponding to different ranges of variation Re (different types of flow) can thereby, in fact, be replaced by a single calculation using one model for the whole of the trajectory for the descent (entry) of space vehicles and natural cosmic bodies (meteoroids) into the atmosphere.     

绝热流Chaplygin气体的Riemann问题

本文研究了绝热流Chaplygin气体动力学方程组,利用特征分析方法,在得到所有基本波的基础上,构造出Riemann问题的所有解.Riemann解由前向疏散波(激波)、后向疏散波(激波)、接触间断以及δ波构成.     

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.     

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.     

Traveling waves of an elliptic-hyperbolic model of phase transitions via varying viscosity-capillarity

We consider an elliptic-hyperbolic model of phase transitions and we show that any Lax shock can be approximated by a traveling wave with a suitable choice of viscosity and capillarity. By varying viscosity and capillarity coefficients, we can cover any Lax shock which either remains in the same phase, or admits a phase transition. The argument used in this paper extends the one in our earlier works. The method relies on LaSalle?s invariance principle and on estimating attraction region of the asymptotically stable of the associated autonomous system of differential equations. We will show that the saddle point of this system of differential equations lies on the boundary of the attraction region and that there is a trajectory leaving the saddle point and entering the attraction region. This gives us a traveling wave connecting the two states of the Lax shock. We also present numerical illustrations of traveling waves.     

A new averaging scheme for the riemann problem in pure water

The numerical investigation of shock phenomena in gas or liquid media where enthalpy is the preferred thermodynamic variable poses special problems. When an expression for internal energy is available, the usual procedure is to employ a splitting scheme to remove source terms from the Euler equations, then upwind-biased shock capturing algorithms are built around the Riemann problem for the conservative system which remains. However, when the governing equations arc formulated in terms of total enthalpy, treatment of a pressure time derivative as a source term leads to a Riemann problem for a system where one equation is not a conservation law. The present research establishes that successful upwind-biased shock capturing schemes can be based upon the pseudo-conservative system. A new averaging scheme for solving the associated Riemann problem is developed. The method is applied to numerical simulations of shock wave propagation in pure water.     

Weak shock waves and its interaction with characteristic shocks in polyatomic gas

We present Lie symmetry analysis for investigating the shock‐wave structure of hyperbolic differential equations of polyatomic gases. With the application of symmetry analysis, we derive particular exact group invariant solutions for the governing system of partial differential equations (PDEs). In the next step, the evolutionary behavior of weak shock along with the characteristic shock and their interaction is investigated. Finally, the amplitudes of reflected wave, transmitted wave, and the jump in shock acceleration influenced by the incident wave after interaction are evaluated for the considered system of equations.     

Vanishing pressure limits of Riemann solutions to the isentropic relativistic Euler system for Chaplygin gas

The Riemann solutions to the isentropic relativistic Euler system for Chaplygin gas with a small parameter are considered. Unlike the polytropic or barotropic gas cases, we find that firstly, as the parameter decreases to a certain critical number, the two-shock solution converges to a delta shock wave solution of the same system. Moreover, as the parameter goes to zero, that is, the pressure vanishes, the solution is nothing but the delta shock wave solution to the zero-pressure relativistic Euler system. Meanwhile, the two-rarefaction wave solution tends to the vacuum solution to the zero-pressure relativistic system, and the solution containing one rarefaction wave and one shock wave tends to the contact discontinuity solution to the zero-pressure relativistic system as pressure vanishes.     

Hyperbolic conservation laws with a moving source

The purpose of this paper is to investigate the wave behavior of hyperbolic conservation laws with a moving source. When the speed of the source is close to one of the characteristic speeds of the system, nonlinear resonance occurs and instability may result. We will study solutions with a single transonic shock wave for a general system u_t + f(u)_x = g(x, u). Suppose that the i^th characteristic speed is close to zero. We propose the following stability criteria: Here l_i and r_i are the i^th normalized left and right eigenvectors of , respectively. Through the local analysis on the evolution of the speed and strength of the transonic shock wave, the above criterion can be justified. It turns out that the speed of the transonic shock wave is monotone increasing (decreasing) most of the time in the unstable (stable) case. This is shown by introducing a global functional on nonlinear wave interactions, based on the Glimm scheme. In particular, together with the local analysis, we can study the shock speed globally. Such a global approach is absent in the previous works. Using this strategy, we prove the existence of solutions and verify the asymptotic stability (or instability). © 1999 John Wiley & Sons, Inc.     

内机械激波——海洋激流的一种解释

采用Lagrange坐标和Hamilton原理,推导了二维两层浅水系统的位移法内波方程,并在此基础上研究了二维内机械激波.通过具体的数值算例分析发现内机械激波具有流速大、持续时间短、空间范围狭小、水面存在突变的特点,指出海洋激流就是内机械激波.内机械激波同样也为海洋断崖提供了一种解释.     

Nonlinear geometric optics for contact discontinuities and shock waves in one-dimensional Euler system for gas dynamics

The aim of this paper is to study the rigorous theory of nonlinear geometric optics for a contact discontinuity and a shock wave to the Euler system for one-dimensional gas dynamics. For the problem of a contact discontinuity and a shock wave perturbed by a small amplitude, high frequency oscillatory wave train, under suitable stability assumptions, we obtain that the perturbed problem has still a shock wave and a contact discontinuity, and we give their asymptotic expansions.