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.     

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.     

Global entropy solutions to the relativistic Euler equations for a class of large initial data

We are concerned with global entropy solutions to the relativistic Euler equations for a class of large initial data which involve the interaction of shock waves and rarefaction waves. We first carefully analyze the global behavior of the shock curves, the rarefaction wave curves, and their corresponding inverse curves in the phase plane. Based on these analyses, we use the Glimm scheme to construct global entropy solutions to the relativistic Euler equations for the class of large discontinuous initial data.     

Global entropy solutions to the relativistic Euler equations for a class of large initial data

We are concerned with global entropy solutions to the relativistic Euler equations for a class of large initial data which involve the interaction of shock waves and rarefaction waves. We first carefully analyze the global behavior of the shock curves, the rarefaction wave curves, and their corresponding inverse curves in the phase plane. Based on these analyses, we use the Glimm scheme to construct global entropy solutions to the relativistic Euler equations for the class of large discontinuous initial data.Received: May 23, 2004     

A phase decomposition approach and the Riemann problem for a model of two-phase flows

We present a phase decomposition approach to deal with the generalized Rankine–Hugoniot relations and then the Riemann problem for a model of two-phase flows. By investigating separately the jump relations for equations in conservative form in the solid phase, we show that the volume fractions can change only across contact discontinuities. Then, we prove that the generalized Rankine–Hugoniot relations are reduced to the usual form. It turns out that shock waves and rarefaction waves remain on one phase only, and the contact waves serve as a bridge between the two phases. By decomposing Riemann solutions into each phase, we show that Riemann solutions can be constructed for large initial data. Furthermore, the Riemann problem admits a unique solution for an appropriate choice of initial data.     

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

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

Bifurcation of travelling wave solutions for the modified dispersive water wave equation

The method of bifurcation of planar dynamical systems and method of numerical simulation of differential equations are employed to investigate the modified dispersive water wave equation. We obtain the parameter bifurcation sets that divide the parameter space into different regions which correspond to qualitatively different phase portraits. In different regions, different types of travelling solutions including solitary wave solutions, shock wave solutions and periodic wave solutions are simulated. Furthermore, with a generalized projective Riccati equation method, several new explicit exact solutions are obtained.     

Numerical computation of viscous profiles for hyperbolic conservation laws

Viscous profiles of shock waves in systems of conservation laws can be viewed as heteroclinic orbits in associated systems of ordinary differential equations (ODE). In the case of overcompressive shock waves, these orbits occur in multi-parameter families. We propose a numerical method to compute families of heteroclinic orbits in general systems of ODE. The key point is a special parameterization of the heteroclinic manifold which can be understood as a generalized phase condition; in the case of shock profiles, this phase condition has a natural interpretation regarding their stability. We prove that our method converges and present numerical results for several systems of conservation laws. These examples include traveling waves for the Navier-Stokes equations for compressible viscous, heat-conductive fluids and for the magnetohydrodynamics equations for viscous, heat-conductive, electrically resistive fluids that correspond to shock wave solutions of the associated ideal models, i.e., the Euler, resp. Lundquist, equations.

     

Discontinuous travelling wave solutions for certain hyperbolic systems

In an earlier paper on a malignant cell invasion model (Marchantet al., SIAM J. Appl. Math, 60, 2000) we introduced a novelform of discontinuous travelling wave solution. These solutionscould be studied easily by combining behaviour within a phaseplane with the Rankine–Hugoniot shock conditions, whichdescribe properties (such as the ratio of the jump discontinuitiesto the speed of propagation) that solutions may possess. Theseresults were new for several reasons. The shock conditions relateto hyperbolic equations (which the model is) but were appliedin a travelling wave ordinary differential equation phase planeusing techniques that usually apply to parabolic reaction–diffusionsystems. In addition the solutions possess singular behaviournear several points in the phase plane but in spite of thisthere exists a robust and stable family of physically interestingsolutions. In this paper we discuss two previously studied models, oneof detonation theory and one of angiogenesis. We show that eachof these models also possesses a family of discontinuous travellingwave solutions which was not previously discovered. Of particularinterest is the solution which has a blunt interface at thefront of the invading profile. In all three models it is thissolution that is seen to stably evolve from physically relevantinitial data, and for physically relevant parameter values. This work confirms the robustness of these novel travellingwave solutions and their applicability to a wider range of mathematicalmodelling situations.     

Formation of singularities for viscosity solutions of hamilton—jacobi equations in one space variable

In this work we study the generation and propagation of singularities (shock waves) of the solution of the Cauchy problem for Hamilton-Jacobi equations in one space variable, under no assumption on the convexity or concavity of the hamiltonian. We study the problem in the class of viscosity solutions, which is the correct class of weak solutions. We obtain the exact global structure of the shock waves by studying the way the characteristics cross. We construct the viscosity solution by either selecting a single-valued branch of the multi-valued function given as a solution by the method of characteristics or constructing explicitly the proper rarefaction waves.     

Supersonic flow onto a solid wedge

We consider the problem of two‐dimensional supersonic flow onto a solid wedge, or equivalently in a concave corner formed by two solid walls. For mild corners, there are two possible steady state solutions, one with a strong and one with a weak shock emanating from the corner. The weak shock is observed in supersonic flights. A longstanding natural conjecture is that the strong shock is unstable in some sense. We resolve this issue by showing that a sharp wedge will eventually produce weak shocks at the tip when accelerated to a supersonic speed. More precisely, we prove that for upstream state as initial data in the entire domain, the time‐dependent solution is self‐similar, with a weak shock at the tip of the wedge. We construct analytic solutions for self‐similar potential flow, both isothermal and isentropic with arbitrary γ ≥ 1. In the process of constructing the self‐similar solution, we develop a large number of theoretical tools for these elliptic regions. These tools allow us to establish large‐data results rather than a small perturbation. We show that the wave pattern persists as long as the weak shock is supersonic‐supersonic; when this is no longer true, numerics show a physical change of behavior. In addition, we obtain rather detailed information about the elliptic region, including analyticity as well as bounds for velocity components and shock tangents. © 2007 Wiley Periodicals, Inc.     

The Geometry of 2 × 2 Systems of Conservation Laws

We consider one typical two-parameter family of quadratic systems of 2 × 2 conservation laws, and study the geometry of the behaviour of the possible solutions of the Riemann problem near an umbilic point, following the geometric approach presented by Isaacson, Marchesin, Palmeira, Plohr, in A global formalism for nonlinear waves in conservation laws, Commun. Math. Phys. (1992). The corresponding phase portraits for the rarefaction curves, shock curves and composite curves are discussed. Financial support from FCT and Calouste Gulbenkian Foundation.     

Converging strong shock waves in magnetogasdynamics under isothermal condition

The similarity solutions to the problem of cylindrically symmetric strong shock waves in an ideal gas with a constant azimuthal magnetic field are presented. The flow behind the shock wave is assumed to spatially isothermal rather than adiabatic. We use the method of Lie group invariance to determine the possible class of self-similar solutions. Infinitesimal generators of Lie group transformations are determined by using the invariance surface conditions to the system and on the basis of arbitrary constants occurring in the expressions for the generators, four different possible cases of the solutions are reckoned and we observed that only two out of all possibilities hold self-similar solutions, one of which follows the power law and another follows the exponential law. To obtain the similarity exponents numerical calculations have been performed and comparison is made with the existing results in the literature. The flow patterns behind the shock are analyzed graphically.

     

广义Chaplygin气体磁流体力学方程组的Riemann问题

利用特征分析和相平面分析的方法,由Rankine-Hugoniot条件和稳定性条件,构造性地得到了一维等熵广义Chaplygin气体磁流体力学方程组的Riemann解的存在唯一性.同时,详细研究了疏散波曲线和激波曲线的性质.     

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.     

Separation of solutions of quasilinear hyperbolic systems in the neighborhood of a large discontinuity

We consider solutions of quasilinear hyperbolic systems of any dimension in one space variable. Locally, our solutions are differentiable excapt for a single jump discontinuity, either an entropy shock or a contacr discontinuity, which may be of any strengrth. With some additional assumptions. we show that the map of the initrial into the solutions at some later time is continuous in L¹     

Uniqueness of transonic shock solutions in a duct for steady potential flow

We study the uniqueness of solutions with a transonic shock in a duct in a class of transonic shock solutions, which are not necessarily small perturbations of the background solution, for steady potential flow. We prove that, for given uniform supersonic upstream flow in a straight duct, there exists a unique uniform pressure at the exit of the duct such that a transonic shock solution exists in the duct, which is unique modulo translation. For any other given uniform pressure at the exit, there exists no transonic shock solution in the duct. This is equivalent to establishing a uniqueness theorem for a free boundary problem of a partial differential equation of second order in a bounded or unbounded duct. The proof is based on the maximum/comparison principle and a judicious choice of special transonic shock solutions as a comparison solution.     

Reflection of Shock Fronts in a van der Waals Fluid

In this paper, the reflection phenomenon of a vapor shock front （both sides of the front are in the vapor phase） in a van der Waals fluid is considered. Both the 1-dimensional case and the multidimensional case are investigated. The authors find that under certain conditions, the reflected wave can be a single shock, or a single subsonic phase boundary, or one weak shock together with one subsonic phase boundary, which depends on the strength of the incident shock. This is different from the known result for the reflection of shock fronts in a gas dynamical system due to Chen in 1989.