Modulation theory solution for nonlinearly resonant,fifth‐order Korteweg–de Vries,nonclassical, traveling dispersive shock waves

A new class of resonant dispersive shock waves was recently identified as solutions of the Kawahara equation— a Korteweg–de Vries (KdV) type nonlinear wave equation with third‐ and fifth‐order spatial derivatives— in the regime of nonconvex, linear dispersion. Linear resonance resulting from the third‐ and fifth‐order terms in the Kawahara equation was identified as the key ingredient for nonclassical dispersive shock wave solutions. Here, nonlinear wave (Whitham) modulation theory is used to construct approximate nonclassical traveling dispersive shock wave (TDSW) solutions of the fifth‐ order KdV equation without the third derivative term, hence without any linear resonance. A self‐similar, simple wave modulation solution of the fifth order, weakly nonlinear KdV–Whitham equations is obtained that matches a constant to a heteroclinic traveling wave via a partial dispersive shock wave so that the TDSW is interpreted as a nonlinear resonance. The modulation solution is compared with full numerical solutions, exhibiting excellent agreement. The TDSW is shown to be modulationally stable in the presence of sufficiently small third‐order dispersion. The Kawahara–Whitham modulation equations transition from hyperbolic to elliptic type for sufficiently large third‐order dispersion, which provides a possible route for the TDSW to exhibit modulational instability.     

Zero-dispersion limit to the Korteweg-de Vries equation: a dressing chain approach

An asymptotic solution of the KdV equation with small dispersion is studied for the case of smooth hump-like initial condition with monotonically decreasing slopes. Despite the well-known approaches by Lax-Levermore and Gurevich-Pitaevskii, a new way of constructing the asymptotics is proposed using the inverse scattering transform together with the dressing chain technique developed by A. Shabat [1]. It provides the Whitham-type approximaton of the leading term by solving the dressing chain through a finite-gap asymptotic ansatz. This yields the Whitham equations on the Riemann invariants together with hodograph transform which solves these equations explicitly. Thus we reproduce an uniform in x asymptotics consisting of smooth solution of the Hopf equation outside the oscillating domain and a slowly modulated cnoidal wave within the domain. Finally, the dressing chain technique provides the proof of an asymptotic estimate for the leading term.      

Oscillations of the zero dispersion limit of the korteweg-de vries equation

We attack the multiphase averaged systems for the zero dispersion limit of the KdV equation. Attention is paid to the most important case—the single phase oscillations. A scheme is developed to solve the Whitham averaged system (single phase averaged system). This system, under our scheme, is transformed to a linear over-determined system of Euler-Poisson-Darboux type whose solution can be written down explicitly. We show that, for any smooth initial data which has only one hump or is a nontrivial monotone function, the weak limit has single-phase oscillations within a cusp in the x-t plane for a short time after the breaking time for the corresponding Burgers equation. Outside the cusp, the limit satisfies the Burgers equation. More surprisingly, we also show that the weak limit has global single-phase oscillations within a cusp for any smooth nontrivial monotone initial data with only one inflection point. © 1993 John Wiley & Sons, Inc.     

Whitham Equations, Bergman Kernel and Lax—Levermore Minimizer

We study multiphase solutions of the Whitham equations. The Whitham equations describe the zero dispersion limit of the Cauchy problem for the Korteweg—de Vries (KdV) equation. The zero dispersion solution of the KdV equation is determined by the Lax—Levermore minimization problem. The minimizer is a measurable function on the real line. When the support of the minimizer consists of a finite number of disjoint intervals to be determined, the minimization problem can be reduced to a scalar Riemann Hilbert (RH) problem. For each fixed x and t 0, the end-points of the contour are determined by the solution of the Whitham equations. The Lax—Levermore minimizer and the solution of the Whitham equations are described in terms of a kernel related to the Bergman kernel. At t = 0 the support of the minimizer consists of one interval for any value of x, while for t > 0, the number of intervals is larger than one in some regions of the (x,t) plane where the multiphase solutions of the Whitham equations develop. The increase of the number of intervals happens whenever the solution of the Whitham equations has a point of gradient catastrophe. For a class of smooth monotonically increasing initial data, we show that the support of the Lax—Levermore minimizer increases or decreases the number of its intervals by one near each point of gradient catastrophe. This result justifies the formation and extinction of the multiphase solutions of the Whitham equations. Furthermore we characterize a class of initial data for which all the points of gradient catastrophe occur only a finite number of times and therefore the support of the Lax—Levermore minimizer consists of a finite number of disjoint intervals for any x and t 0. This corresponds to give an upper bound to the genus of the solution of the Whitham equations. Similar results are obtained for the semi-classical limit of the defocusing nonlinear Schrödinger equation.     

The global attractor of the viscous Fornberg–Whitham equation

This paper aims to present a proof of the existence of the attractor for the one-dimensional viscous Fornberg–Whitham equation. In this paper, the global existence of solution to the viscous Fornberg–Whitham equation in L² under the periodic boundary conditions is studied. By using the time estimate of the Fornberg–Whitham equation, we get the compact and bounded absorbing set and the existence of the global attractor for the viscous Fornberg–Whitham equation.     

Breakdown of the Whitham Modulation Theory and the Emergence of Dispersion

The Whitham modulation theory for periodic traveling waves of PDEs generated by a Lagrangian produces first‐order dispersionless PDEs that are, generically, either hyperbolic or elliptic. In this paper, degeneracy of the Whitham equations is considered where one of the characteristic speeds is zero. In this case, the Whitham equations are no longer valid. Reformulation and rescaling show that conservation of wave action morphs into the Korteweg–de Vries (KdV) equation on a longer time scale thereby generating dispersion in the Whitham modulation equations even for finite amplitude waves.     

Temporal asymptotics for soliton equations in problems with step initial conditions

In this survey, we present modern approaches to the construction and justification of large time asymptotics for solutions of main soliton equations with step-like initial condition whose boundary conditions as x ± are finite-gap, quasi-periodic solutions. The principal term of the asymptotic is also a finite-gap, quasi-periodic solution whose phase vectors are modulated with respect to the slow space-like variable. The Whitham equations describing this modulation are studied in detail. For the KdV equations, we construct and justify the principal term of the asymptotic for arbitrary finite-gap boundary conditions. By examining the sine-Gordon equation, we study the case of boundary conditions with complex-valued, self-conjugated quasi-periods. We prove the existence and uniqueness theorems in the case of the complex-valued group velocities that appear in the Whitham equations in the case considered. We present a complete picture (uniform in x) of the Whitham deformation for the case of 1-gap boundary conditions in the sine-Gordon equation.Translated from Sovremennaya Matematika i Ee Prilozheniya (Contemporary Mathematics and Its Applications), Vol. 5, Asymptotic Methods, 2003.     

Unified approach to KdV modulations

We develop a unified approach to integrating the Whitham modulation equations. Our approach is based on the formulation of the initial‐value problem for the zero‐dispersion KdV as the steepest descent for the scalar Riemann‐Hilbert problem [6] and on the method of generating differentials for the KdV‐Whitham hierarchy [9]. By assuming the hyperbolicity of the zero‐dispersion limit for the KdV with general initial data, we bypass the inverse scattering transform and produce the symmetric system of algebraic equations describing motion of the modulation parameters plus the system of inequalities determining the number the oscillating phases at any fixed point on the (x, t)‐plane. The resulting system effectively solves the zero‐dispersion KdV with an arbitrary initial datum. © 2001 John Wiley & Sons, Inc.     

Construction and study of high-order accurate schemes for solving the one-dimensional heat equation

The method of undetermined coefficients on multipoint stencils with two time levels was used to construct compact difference schemes of O(τ³, h ⁶) accuracy intended for solving boundary value problems for the one-dimensional heat equation. The schemes were examined for von Neumann stability, and numerical experiments were conducted on a sequence of grids with mesh sizes tending to zero. One of the schemes was proved to be absolutely stable. It was shown that, for smooth solutions, the high order of convergence of the numerical solution agrees with the order of accuracy; moreover, solutions accurate up to ～10^?12 are obtained on grids with spatial mesh sizes of ～10^?2. The formulas for the schemes are rather simple and easy to implement on a computer.     

Uniform grid approximation of nonsmooth solutions to the mixed boundary value problem for a singularly perturbed reaction-diffusion equation in a rectangle

A mixed boundary value problem for a singularly perturbed reaction-diffusion equation in a square is considered. A Neumann condition is specified on one side of the square, and a Dirichlet condition is set on the other three. It is assumed that the coefficient of the equation, its right-hand side, and the boundary values of the desired solution or its normal derivative on the sides of the square are smooth enough to ensure the required smoothness of the solution in a closed domain outside the neighborhoods of the corner points. No compatibility conditions are assumed to hold at the corner points. Under these assumptions, the desired solution in the entire closed domain is of limited smoothness: it belongs only to the Hölder class C ^μ, where μ ∈ (0, 1) is arbitrary. In the domain, a nonuniform rectangular mesh is introduced that is refined in the boundary domain and depends on a small parameter. The numerical solution to the problem is based on the classical five-point approximation of the equation and a four-point approximation of the Neumann boundary condition. A mesh refinement rule is described under which the approximate solution converges to the exact one uniformly with respect to the small parameter in the L _∞ ^h norm. The convergence rate is O(N ^?2ln² N), where N is the number of mesh nodes in each coordinate direction. The parameter-uniform convergence of difference schemes for mixed problems without compatibility conditions at corner points was not previously analyzed.     

Spatial Analyticity of Solutions to Integrable Systems. I. The KdV Case

We are concerned with the Cauchy problem for the KdV equation for nonsmooth locally integrable initial profiles q's which are, in a certain sense, essentially bounded from below and q(x) = O(e^?cx? ), x → + ∞, with some positive c and ?. Using the inverse scattering transform, we show that the KdV flow turns such initial data into a function which is (1) meromorphic (in the space variable) on the whole complex plane if ? > 1/2, (2) meromorphic on a strip around the real line if ? = 1/2, and (3) Gevrey regular if ? <1/2. Note that q's need not have any decay or pattern of behavior at ? ∞.     

The generation,propagation, and extinction of multiphases in the KdV zero‐dispersion limit

We study the multiphases in the KdV zero‐dispersion limit. These phases are governed by the Whitham equations, which are 2g + 1 quasi‐linear hyperbolic equations where g is the number of phases. We are interested in both the interaction of two single phases and the breaking of a single phase for general initial data. We analyze in detail how a double phase is generated from the interaction or breaking, how it propagates in space‐time, and how it collapses to a single phase in a finite time. The Whitham equations are known to be integrable via a hodograph transform. The crucial step in our approach is to formulate the hodograph transform in terms of the Euler‐Poisson‐Darboux solutions. Under our scheme, the zeros of the Jacobian of the transform are given by the zeros of the Euler‐Poisson‐Darboux solution. Hence, the problem of inverting the hodograph transform to give the Whitham solution reduces to that of counting the zeros of the Euler‐Poisson‐Darboux solution. © 2002 Wiley Periodicals, Inc.     

Computable analysis of a boundary‐value problem for the generalized KdV–Burgers equation

In this paper, we investigate the computability of the solution operator of the generalized KdV‐Burgers equation with initial‐boundary value problem. Here, the solution operator is a nonlinear map H^{3m ? 1}(R⁺) × H^m(0,T)→C([0,T];H^{3m ? 1}(R⁺)) from the initial‐boundary value data to the solution of the equation. By a technique that is widely used for the study of nonlinear dispersive equation, and using the type 2 theory of effectivity as computable model, we prove that the solution map is Turing computable, for any integer m ≥ 2, and computable real number T > 0. Copyright © 2014 John Wiley & Sons, Ltd.     

Strong stability of a scheme on locally uniform meshes for a singularly perturbed ordinary differential convection-diffusion equation

The Dirichlet problem for a singularly perturbed ordinary differential convection-diffusion equation with a small parameter ? (? ?? (0, 1]) multiplying the higher order derivative is considered. For the problem, a difference scheme on locally uniform meshes is constructed that converges in the maximum norm conditionally, i.e., depending on the relation between the parameter ? and the value N defining the number of nodes in the mesh used; in particular, the scheme converges almost ?-uniformly (i.e., its accuracy depends weakly on ?). The stability of the scheme with respect to perturbations in the data and its conditioning are analyzed. The scheme is constructed using classical monotone approximations of the boundary value problem on a priori adapted grids, which are uniform on subdomains where the solution is improved. The boundaries of these subdomains are determined by a majorant of the singular component of the discrete solution. On locally uniform meshes, the difference scheme converges at a rate of O(min[?^?1 N ^?K lnN, 1] + N ^?1lnN), where K is a prescribed number of iterations for refining the discrete solution. The scheme converges almost ?-uniformly at a rate of O(N ^?1lnN) if N ^?1 ?? ?^??, where ?? (the defect of ?-uniform convergence) determines the required number K of iterations (K = K(??) ?? ??^?1) and can be chosen arbitrarily small from the half-open interval (0, 1]. The condition number of the difference scheme satisfies the bound ?? _P = O(?^?1/K ln^1/K ?^?1??^{?(K + 1)/K} ), where ?? is the accuracy of the solution of the scheme in the maximum norm in the absence of perturbations. For sufficiently large K, the scheme is almost ?-uniformly strongly stable.     

Generalized solitary waves in a finite‐difference Korteweg‐de Vries equation

Generalized solitary waves with exponentially small nondecaying far field oscillations have been studied in a range of singularly perturbed differential equations, including higher order Korteweg‐de Vries (KdV) equations. Many of these studies used exponential asymptotics to compute the behavior of the oscillations, revealing that they appear in the solution as special curves known as Stokes lines are crossed. Recent studies have identified similar behavior in solutions to difference equations. Motivated by these studies, the seventh‐order KdV and a hierarchy of higher order KdV equations are investigated, identifying conditions which produce generalized solitary wave solutions. These results form a foundation for the study of infinite‐order differential equations, which are used as a model for studying lattice equations. Finally, a lattice KdV equation is generated using finite‐difference discretization, in which a lattice generalized solitary wave solution is found.     

Quasi-classical approach to the inverse scattering problem for the KdV equation,and solution of the Whitham modulation equations

We consider an initial value problem for the KdV equation in the limit of weak dispersion. This model describes the formation and evolution in time of a nondissipative shock wave in plasma. Using the perturbation theory in power series of a small dispersion parameter, we arrive at the Riemann simple wave equation. Once the simple wave is overturned, we arrive at the system of Whitham modulation equations that describes the evolution of the resulting nondissipative shock wave. The idea of the approach developed in this paper is to study the asymptotic behavior of the exact solution in the limit of weak dispersion, using the solution given by the inverse scattering problem technique. In the study of the problem, we use the WKB approach to the direct scattering problem and use the formulas for the exact multisoliton solution of the inverse scattering problem. By passing to the limit, we obtain a finite set of relations that connects the space-time parameters x, t and the modulation parameters of the nondissipative shock wave.Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 106, No. 1, pp. 44–61, January, 1996.     

Finite difference schemes for the singularly perturbed reaction-diffusion equation in the case of spherical symmetry

The boundary value problem for the singularly perturbed reaction-diffusion parabolic equation in a ball in the case of spherical symmetry is considered. The derivatives with respect to the radial variable appearing in the equation are written in divergent form. The third kind boundary condition, which admits the Dirichlet and Neumann conditions, is specified on the boundary of the domain. The Laplace operator in the differential equation involves a perturbation parameter ?², where ? takes arbitrary values in the half-open interval (0, 1]. When ? → 0, the solution of such a problem has a parabolic boundary layer in a neighborhood of the boundary. Using the integro-interpolational method and the condensing grid technique, conservative finite difference schemes on flux grids are constructed that converge ?-uniformly at a rate of O(N ^?2ln² N + N ₀ ^?1 ), where N + 1 and N ₀ + 1 are the numbers of the mesh points in the radial and time variables, respectively.     

Finite element method for Grwünwald–Letnikov time-fractional partial differential equation

In this article, we consider the finite element methods (FEM) for Grwünwald–Letnikov time-fractional diffusion equation, which is obtained from the standard two-dimensional diffusion equation by replacing the first-order time derivative with a fractional derivative (of order α, with 0?h ^r+1?+?τ^2-α), where h, τ and r are the space step size, time step size and polynomial degree, respectively. A numerical example is presented to verify the order of convergence.     

On numerical study of periodic solutions of a delay equation in biological models

Some results are presented of the numerical study of periodic solutions of a nonlinear equation with a delayed argument in connection with themathematical models having real biological prototypes. The problem is formulated as a boundary value problem for a delay equation with the conditions of periodicity and transversality. A spline-collocation finite-difference scheme of the boundary value problem using a Hermitian interpolation cubic spline of the class C ¹ with fourth order error is proposed. For the numerical study of the system of nonlinear equations of the finitedifference scheme, the parameter continuation method is used, which allows us to identify possible nonuniqueness of the solution of the boundary value problem and, hence, the nonuniqueness of periodic solutions regardless of their stability. By examples it is shown that the periodic oscillations occur for the parameter values specific to the real molecular-genetic systems of higher species, for which the principle of delay is quite easy to implement.     

Nonexistence for the Laplace equation with a dynamical boundary condition of fractional type

We consider the Laplace equation in ? ^d?1 × ?⁺ × (0,+∞) with a dynamical nonlinear boundary condition of order between 1 and 2. Namely, the boundary condition is a fractional differential inequality involving derivatives of noninteger order as well as a nonlinear source. Nonexistence results and necessary conditions are established for local and global existence. In particular, we show that the critical exponent depends only on the fractional derivative of the least order.