Existence of traveling waves in compressible Euler equations with viscosity and capillarity

Given any Lax shock of the compressible Euler dynamics equations, we show that there exists the corresponding traveling wave of the system when viscosity and capillarity are suitably added. For a traveling wave corresponding to a given Lax shock, the governing viscous–capillary system is reduced to a system of two differential equations of first-order, which admits an asymptotically stable equilibrium point and a saddle point. We then develop the method of estimating attraction domain of the asymptotically stable equilibrium point for the compressible Euler equations and show that the saddle point in fact lies on the boundary of this set. Then, we establish a saddle-to-stable connection by pointing out that there is a stable trajectory leaving the saddle point and entering the attraction domain of the asymptotically stable equilibrium point. This gives us a traveling wave of the viscous–capillary compressible Euler equations.     

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.     

Existence of traveling waves associated with Lax shocks which violate Oleinik’s entropy criterion

This paper answers to the question whether a shock wave in conservation laws satisfying the Lax shock inequalities but not Oleinik’s entropy criterion is admissible under the vanishing viscosity-capillarity effects. Such a shock appears in van der Waals fluids when a secant line meets the graph of the flux function at four distinct points, and the shock jumps between the two farthest points. The existence of the corresponding traveling waves would justify the admissibility of the shock. For this purpose, we will first show that the corresponding traveling waves satisfy a system of differential equations with two saddle points and two asymptotically stable points. Second, we estimate the domains of attraction of the asymptotically stable equilibrium points, relying on Lyapunov’s stability theory. Third, we investigate the circumstances when an unstable trajectory leaving the saddle point corresponding to the left-hand state of the shock will ever enter the domain of attraction of each of the two asymptotically stable equilibrium points. Finally, we establish the existence of traveling waves associated with a Lax shock but violating the Oleinik’s entropy criterion.     

Existence of traveling waves to any Lax shock satisfying Oleinik’s criterion in conservation laws

Given any shock wave of a conservation law where the flux function may not be convex, we want to know whether it is admissible under the criterion of vanishing viscosity/capillarity effects. In this work, we show that if the shock satisfies the Oleinik’s criterion and the Lax shock inequalities, then for an arbitrary diffusion coefficient, we can always find suitable dispersion coefficients such that the diffusive-dispersive model admits traveling waves approximating the given shock. The paper develops the method of estimating attraction domain for traveling waves we have studied before.     

Global stability of the attracting set of an enzyme-catalysed reaction system

The essential feature of enzymatic reactions is a nonlinear dependency of reaction rate on metabolite concentration taking the form of saturation kinetics. Recently, it has been shown that this feature is associated with the phenomenon of “loss of system coordination” [1]. In this paper, we study a system of ordinary differential equations representing a branched biochemical system of enzyme-mediated reactions. We show that this system can become very sensitive to changes in certain maximum enzyme activities. In particular, we show that the system exhibits three distinct responses: a unique, globally-stable steady-state, large amplitude oscillations, and asymptotically unbounded solutions, with the transition between these states being almost instantaneous. It is shown that the appearance of large amplitude, stable limit cycles occurs due to a “false” bifurcation or canard explosion. The subsequent disappearance of limit cycles corresponds to the collapse of the domain of attraction of the attracting set for the system and occurs due to a global bifurcation in the flow, namely, a saddle connection. Subsequently, almost all nonnegative data become unbounded under the action of the dynamical system and correspond exactly to loss of system coordination. We discuss the relevance of these results to the possible consequences of modulating such systems.     

Global existence of diffusive-dispersive traveling waves for general flux functions

We establish a global existence of traveling waves for diffusive-dispersive conservation laws for locally Lipschitz flux functions. Using Lyapunov stability techniques, we reduce the global problem of finding traveling waves to considering local behaviors of a stable trajectory of the saddle point.     

Nagumo方程行波系统的定性分析及其有界行波解

本文对Nagumo方程的行波系统进行了定性分析,该系统存在一端连接鞍点的有界异宿轨,进而选择奇点为鞍点的平面线性自治系统,利用该平面自治系统轨线向径的斜率,根据齐次平衡原则,构造出了Nagumo方程行波系统的行波解.其次Nagumo方程的行波系统还存在着对应中心周围闭轨的周期解,因而提出新的CPP解法,求出了对应的周期解.     

A mathematical analysis for a model arising from public goods games

We study the effects of altruistic behaviors in a public goods game model which describes the competition between the farmers and the exploiters. Corresponding to different parametric regions, we analyze in detail the stability of the equilibrium states and obtain attraction regions for stable equilibria. Then using the upper–lower solution method and monotone iterations, we further show that for a family of wave speeds, there exist traveling wave solutions connecting one of the unstable states to the stable state. This answers a conjecture made by Wakano in [J.Y. Wakano, A mathematical analysis on public goods games in the continuous space, Math. Biosci. 201 (2006) 72–89]. The results indicate that when the penalty for the altruistic behavior is small, the growth rate of the population determines its survival or extinction states in the long run. Furthermore, if the two populations have the same total growth rate, altruism in the competition leads to a wide range of co-existent states. Numerical simulations are also presented to illustrate the theoretical results.     

Connection across a Separatrix with Dissipation

Dissipative perturbations of strongly nonlinear oscillators that correspond to slowly varying double-well potentials are considered. The method of averaging, which describes the solution as nearly periodic, fails as the trajectory approaches the unperturbed separatrix, a homoclinic orbit of the saddle point, significantly before it is captured in either well. Nevertheless, perturbed initial conditions corresponding to the boundary of the basin of attraction for each well, which are the perturbed stable manifolds of the saddle point, are accurately determined using only the method of averaging modified by Melnikov energy ideas near the separatrix. To determine the amplitude and phase of the captured oscillations after crossing the separatrix, a transition region is constructed consisting of a large sequence of nearly solitary pulses along the separatrix. The amplitude and phases of the slowly varying nonlinear oscillations away from the separatrix, both before and after capture, are matched to this transition region. In this way, analytic connection formulas across the separatrix are obtained and are shown to depend on the perturbed initial conditions.     

Smooth and non-smooth traveling wave solutions of some generalized Camassa–Holm equations

In this paper we employ two recent analytical approaches to investigate the possible classes of traveling wave solutions of some members of a recently-derived integrable family of generalized Camassa–Holm (GCH) equations. A recent, novel application of phase-plane analysis is employed to analyze the singular traveling wave equations of three of the GCH NLPDEs, i.e. the possible non-smooth peakon and cuspon solutions. One of the considered GCH equations supports both solitary (peakon) and periodic (cuspon) cusp waves in different parameter regimes. The second equation does not support singular traveling waves and the last one supports four-segmented, non-smooth M-wave solutions.Moreover, smooth traveling waves of the three GCH equations are considered. Here, we use a recent technique to derive convergent multi-infinite series solutions for the homoclinic orbits of their traveling-wave equations, corresponding to pulse (kink or shock) solutions respectively of the original PDEs. We perform many numerical tests in different parameter regime to pinpoint real saddle equilibrium points of the corresponding GCH equations, as well as ensure simultaneous convergence and continuity of the multi-infinite series solutions for the homoclinic orbits anchored by these saddle points. Unlike the majority of unaccelerated convergent series, high accuracy is attained with relatively few terms. We also show the traveling wave nature of these pulse and front solutions to the GCH NLPDEs.     

A solution formula for a non-convex scalar hyperbolic conservation law with monotone initial data

In this paper we prove an explicit representation formula for the solution of a one-dimensional hyperbolic conservation law with a non-convex flux function but monotone initial data. This representation formula is similar to those of Lax [10] and Kunik [7,8] and enables us to compute the solution pointwise explicitly. This result is a generalization of a theorem given in Kunik [8] where the case of only one inflexion point for the fluxes was considered. Its proof uses the polygonal method of Dafermos [2]. The application of this method leads to a simple explicit construction of the solutions for a Kynch sedimentation process [9] and to an explicit parameter representation for the shock curves evolving during the sedimentation process.     

Nonlinear uzawa methods for solving nonsymmetric saddle point problems

In [A new nonlinear Uzawa algorithm for generalized saddle point problems, Appl. Math. Comput., 175(2006), 1432–1454], a nonlinear Uzawa algorithm for solving symmetric saddle point problems iteratively, which was defined by two nonlinear approximate inverses, was considered. In this paper, we extend it to the nonsymmetric case. For the nonsymmetric case, its convergence result is deduced. Moreover, we compare the convergence rates of three nonlinear Uzawa methods and show that our method is more efficient than other nonlinear Uzawa methods in some cases. The results of numerical experiments are presented when we apply them to Navier-Stokes equations discretized by mixed finite elements.     

A generalization of parameterized inexact Uzawa method for generalized saddle point problems

Recently, a class of parameterized inexact Uzawa methods has been proposed for generalized saddle point problems by Bai and Wang [Z.-Z. Bai, Z.-Q. Wang, On parameterized inexact Uzawa methods for generalized saddle point problems, Linear Algebra Appl. 428 (2008) 2900–2932], and a generalization of the inexact parameterized Uzawa method has been studied for augmented linear systems by Chen and Jiang [F. Chen, Y.-L. Jiang, A generalization of the inexact parameterized Uzawa methods for saddle point problems, Appl. Math. Comput. (2008)]. This paper is concerned about a generalization of the parameterized inexact Uzawa method for solving the generalized saddle point problems with nonzero (2, 2) blocks. Some new iterative methods are presented and their convergence are studied in depth. By choosing different parameter matrices, we derive a series of existing and new iterative methods, including the preconditioned Uzawa method, the inexact Uzawa method, the SOR-like method, the GSOR method, the GIAOR method, the PIU method, the APIU method and so on. Numerical experiments are used to demonstrate the feasibility and effectiveness of the generalized parameterized inexact Uzawa methods.     

A Systematic “Saddle Point Near a Pole” Asymptotic Method with Application to the Gauss Hypergeometric Function

In recent works [ 1 ] and [ 2 ], we have proposed more systematic versions of the Laplace’s and saddle point methods for asymptotic expansions of integrals. Those variants of the standard methods avoid the classical change of variables and give closed algebraic formulas for the coefficients of the expansions. In this work we apply the ideas introduced in [ 1 ] and [ 2 ] to the uniform method “saddle point near a pole.” We obtain a computationally more systematic version of that uniform asymptotic method for integrals having a saddle point near a pole that, in many interesting examples, gives a closed algebraic formula for the coefficients. The asymptotic sequence is given, in general, in terms of exponential integrals of fractional order (or incomplete gamma functions). In particular, when the order of the saddle point is two, the basic approximant is given in terms of the error function (as in the standard method). As an application, we obtain new asymptotic expansions of the Gauss Hypergeometric function ₂F₁(a, b, c; z) for large b and c with c > b + 1 .     

DIFFUSIVE—DISPERSIVE TRAVELING WAVES AND KINETIC RELATIONS IV. COMPRESSIBLE EULER EQUATIONS

The authors consider the Euler equations for a compressible fluid in one space dimension when the equation of state of the fluid does not fulfill standard convexity assumptions and viscosity and capillarity effects are taken into account. A typical example of nonconvex constitutive equation for fluids is Van der Waals' equation. The first order terms of these partial differential equations form a nonlinear system of mixed (hyperbolic-elliptic) type. For a class of nonconvex equations of state, an existence theorem of traveling waves solutions with arbitrary large amplitude is established here. The authors distinguish between classical (compressive) and nonclassical (undercompressive) traveling waves. The latter do not fulfill Lax shock inequalities, and are characterized by the so-called kinetic relation, whose properties are investigated in this paper.     

Slow Passage Through the Nonhyperbolic Homoclinic Orbit of the Saddle-Center Hamiltonian Bifurcation

Slowly varying Hamiltonian systems, for which action is a well-known adiabatic invariant, are considered in the case where the system undergoes a saddle center bifurcation. We analyze the situation in which the solution slowly passes through the nonhyperbolic homoclinic orbit created at the saddle-center bifurcation. The solution near this homoclinic orbit consists of a large sequence of homoclinic orbits surrounded by near approaches to the autonomous nonlinear nonhyperbolic saddle point. By matching this solution to the strongly nonlinear oscillations obtained by averaging before and after crossing the homoclinic orbit, we determine the change in the action. If one orbit comes sufficiently close to the nonlinear saddle point, then that one saddle approach instead satisfies the nonautonomous first Painlevé equation, whose stable manifold of the unstable saddle (created in the saddle-center bifurcation) separates solutions approaching the stable center from those involving sequences of nearly homoclinic orbits.     

Branching random walks,stable point processes and regular variation

Using the theory of regular variation, we give a sufficient condition for a point process to be in the superposition domain of attraction of a strictly stable point process. This sufficient condition is used to obtain the weak limit of a sequence of point processes induced by a branching random walk with jointly regularly varying displacements. Because of heavy tails of the step size distribution, we can invoke a one large jump principle at the level of point processes to give an explicit representation of the limiting point process. As a consequence, we extend the main result of Durrett (1983) and verify that two related predictions of Brunet and Derrida (2011) remain valid for this model.     

The gap lemma and geometric criteria for instability of viscous shock profiles

An obstacle in the use of Evans function theory for stability analysis of traveling waves occurs when the spectrum of the linearized operator about the wave accumulates at the imaginary axis, since the Evans function has in general been constructed only away from the essential spectrum. A notable case in which this difficulty occurs is in the stability analysis of viscous shock profiles. Here we prove a general theorem, the “gap lemma,” concerning the analytic continuation of the Evans function associated with the point spectrum of a traveling wave into the essential spectrum of the wave. This allows geometric stability theory to be applied in many cases where it could not be applied previously. We demonstrate the power of this method by analyzing the stability of certain undercompressive viscous shock waves. A necessary geometric condition for stability is determined in terms of the sign of a certain Melnikov integral of the associated viscous profile. This sign can easily be evaluated numerically. We also compute it analytically for solutions of several important classes of systems. In particular, we show for a wide class of systems that homoclinic (solitary) waves are linearly unstable, confirming these as the first known examples of unstable viscous shock waves. We also show that (strong) heteroclinic undercompressive waves are sometimes unstable. Similar stability conditions are also derived for Lax and overcompressive shocks and for n × n conservation laws, n ≥ 2. © 1998 John Wiley & Sons, Inc.