Peakon

Last updated • 4 min readFrom Wikipedia, The Free Encyclopedia

In the theory of integrable systems, a peakon ("peaked soliton") is a soliton with discontinuous first derivative; the wave profile is shaped like the graph of the function . Some examples of non-linear partial differential equations with (multi-)peakon solutions are the Camassa–Holm shallow water wave equation, the Degasperis–Procesi equation and the Fornberg–Whitham equation. Since peakon solutions are only piecewise differentiable, they must be interpreted in a suitable weak sense. The concept was introduced in 1993 by Camassa and Holm in the short but much cited paper where they derived their shallow water equation. [1]

Contents

A family of equations with peakon solutions

The primary example of a PDE which supports peakon solutions is

where is the unknown function, and b is a parameter. [2] In terms of the auxiliary function defined by the relation , the equation takes the simpler form

This equation is integrable for exactly two values of b, namely b = 2 (the Camassa–Holm equation) and b = 3 (the Degasperis–Procesi equation).

The single peakon solution

The PDE above admits the travelling wave solution , which is a peaked solitary wave with amplitude c and speed c. This solution is called a (single) peakon solution, or simply a peakon. If c is negative, the wave moves to the left with the peak pointing downwards, and then it is sometimes called an antipeakon.

It is not immediately obvious in what sense the peakon solution satisfies the PDE. Since the derivative ux has a jump discontinuity at the peak, the second derivative uxx must be taken in the sense of distributions and will contain a Dirac delta function; in fact, . Now the product occurring in the PDE seems to be undefined, since the distribution m is supported at the very point where the derivative ux is undefined. An ad hoc interpretation is to take the value of ux at that point to equal the average of its left and right limits (zero, in this case). A more satisfactory way to make sense of the solution is to invert the relationship between u and m by writing , where , and use this to rewrite the PDE as a (nonlocal) hyperbolic conservation law:

(The star denotes convolution with respect to x.) In this formulation the function u can simply be interpreted as a weak solution in the usual sense. [3]

Multipeakon solutions

Two-peakon wave profile (solid curve) formed by adding two peakons (dashed curves):
u
=
m
1
e
-
|
x
-
x
1
|
+
m
2
e
-
|
x
-
x
2
|
{\displaystyle u=m_{1}\,e^{-|x-x_{1}|}+m_{2}\,e^{-|x-x_{2}|}} Two-peakon.svg
Two-peakon wave profile (solid curve) formed by adding two peakons (dashed curves):

Multipeakon solutions are formed by taking a linear combination of several peakons, each with its own time-dependent amplitude and position. (This is a very simple structure compared to the multisoliton solutions of most other integrable PDEs, like the Korteweg–de Vries equation for instance.) The n-peakon solution thus takes the form

where the 2n functions and must be chosen suitably in order for u to satisfy the PDE. For the "b-family" above it turns out that this ansatz indeed gives a solution, provided that the system of ODEs

is satisfied. (Here sgn denotes the sign function.) Note that the right-hand side of the equation for is obtained by substituting in the formula for u. Similarly, the equation for can be expressed in terms of , if one interprets the derivative of at x = 0 as being zero. This gives the following convenient shorthand notation for the system:

The first equation provides some useful intuition about peakon dynamics: the velocity of each peakon equals the elevation of the wave at that point.

Explicit solution formulas

In the integrable cases b = 2 and b = 3, the system of ODEs describing the peakon dynamics can be solved explicitly for arbitrary n in terms of elementary functions, using inverse spectral techniques. For example, the solution for n = 3 in the Camassa–Holm case b = 2 is given by [4]

where , and where the 2n constants and are determined from initial conditions. The general solution for arbitrary n can be expressed in terms of symmetric functions of and . The general n-peakon solution in the Degasperis–Procesi case b = 3 is similar in flavour, although the detailed structure is more complicated. [5]

Notes

  1. Camassa & Holm 1993
  2. Degasperis, Holm & Hone 2002
  3. Constantin & McKean 1999 (who treat the Camassa–Holm case b = 2; the general case is very similar)
  4. Beals, Sattinger & Szmigielski 2000 (where a different normalization and sign convention is used)
  5. Lundmark & Szmigielski 2005

Related Research Articles

<span class="mw-page-title-main">Exponential distribution</span> Probability distribution

In probability theory and statistics, the exponential distribution is the probability distribution of the time between events in a Poisson point process, i.e., a process in which events occur continuously and independently at a constant average rate. It is a particular case of the gamma distribution. It is the continuous analogue of the geometric distribution, and it has the key property of being memoryless. In addition to being used for the analysis of Poisson point processes it is found in various other contexts.

In mathematical optimization, the method of Lagrange multipliers is a strategy for finding the local maxima and minima of a function subject to equality constraints. It is named after the mathematician Joseph-Louis Lagrange. The basic idea is to convert a constrained problem into a form such that the derivative test of an unconstrained problem can still be applied. The relationship between the gradient of the function and gradients of the constraints rather naturally leads to a reformulation of the original problem, known as the Lagrangian function.

<span class="mw-page-title-main">Cayley–Hamilton theorem</span> Every square matrix over a commutative ring satisfies its own characteristic equation

In linear algebra, the Cayley–Hamilton theorem states that every square matrix over a commutative ring satisfies its own characteristic equation.

<span class="mw-page-title-main">Legendre function</span>

In physical science and mathematics, the Legendre functionsPλ, Qλ and associated Legendre functionsPμ
λ
, Qμ
λ
, and Legendre functions of the second kind, Qn, are all solutions of Legendre's differential equation. The Legendre polynomials and the associated Legendre polynomials are also solutions of the differential equation in special cases, which, by virtue of being polynomials, have a large number of additional properties, mathematical structure, and applications. For these polynomial solutions, see the separate Wikipedia articles.

The principle of maximum entropy states that the probability distribution which best represents the current state of knowledge about a system is the one with largest entropy, in the context of precisely stated prior data.

<span class="mw-page-title-main">Separation of variables</span> Technique for solving differential equations

In mathematics, separation of variables is any of several methods for solving ordinary and partial differential equations, in which algebra allows one to rewrite an equation so that each of two variables occurs on a different side of the equation.

In mathematics and its applications, classical Sturm–Liouville theory is the theory of real second-order linear ordinary differential equations of the form:

<span class="mw-page-title-main">Method of characteristics</span> Technique for solving hyperbolic partial differential equations

In mathematics, the method of characteristics is a technique for solving partial differential equations. Typically, it applies to first-order equations, although more generally the method of characteristics is valid for any hyperbolic partial differential equation. The method is to reduce a partial differential equation to a family of ordinary differential equations along which the solution can be integrated from some initial data given on a suitable hypersurface.

In the theory of stochastic processes, the Karhunen–Loève theorem, also known as the Kosambi–Karhunen–Loève theorem is a representation of a stochastic process as an infinite linear combination of orthogonal functions, analogous to a Fourier series representation of a function on a bounded interval. The transformation is also known as Hotelling transform and eigenvector transform, and is closely related to principal component analysis (PCA) technique widely used in image processing and in data analysis in many fields.

<span class="mw-page-title-main">Crank–Nicolson method</span> Finite difference method for numerically solving parabolic differential equations

In numerical analysis, the Crank–Nicolson method is a finite difference method used for numerically solving the heat equation and similar partial differential equations. It is a second-order method in time. It is implicit in time, can be written as an implicit Runge–Kutta method, and it is numerically stable. The method was developed by John Crank and Phyllis Nicolson in the mid 20th century.

In mathematics, the inverse scattering transform is a method for solving some non-linear partial differential equations. The method is a non-linear analogue, and in some sense generalization, of the Fourier transform, which itself is applied to solve many linear partial differential equations. The name "inverse scattering method" comes from the key idea of recovering the time evolution of a potential from the time evolution of its scattering data: inverse scattering refers to the problem of recovering a potential from its scattering matrix, as opposed to the direct scattering problem of finding the scattering matrix from the potential.

In physics and mathematics, the solid harmonics are solutions of the Laplace equation in spherical polar coordinates, assumed to be (smooth) functions . There are two kinds: the regular solid harmonics, which are well-defined at the origin and the irregular solid harmonics, which are singular at the origin. Both sets of functions play an important role in potential theory, and are obtained by rescaling spherical harmonics appropriately:

In mathematical physics, the Degasperis–Procesi equation

<span class="mw-page-title-main">Camassa–Holm equation</span>

In fluid dynamics, the Camassa–Holm equation is the integrable, dimensionless and non-linear partial differential equation

Dispersionless limits of integrable partial differential equations (PDE) arise in various problems of mathematics and physics and have been intensively studied in recent literature. They typically arise when considering slowly modulated long waves of an integrable dispersive PDE system.

In numerical analysis, von Neumann stability analysis is a procedure used to check the stability of finite difference schemes as applied to linear partial differential equations. The analysis is based on the Fourier decomposition of numerical error and was developed at Los Alamos National Laboratory after having been briefly described in a 1947 article by British researchers Crank and Nicolson. This method is an example of explicit time integration where the function that defines governing equation is evaluated at the current time. Later, the method was given a more rigorous treatment in an article co-authored by John von Neumann.

In mathematics, the Schauder estimates are a collection of results due to Juliusz Schauder concerning the regularity of solutions to linear, uniformly elliptic partial differential equations. The estimates say that when the equation has appropriately smooth terms and appropriately smooth solutions, then the Hölder norm of the solution can be controlled in terms of the Hölder norms for the coefficient and source terms. Since these estimates assume by hypothesis the existence of a solution, they are called a priori estimates.

Derrick's theorem is an argument by physicist G.H. Derrick which shows that stationary localized solutions to a nonlinear wave equation or nonlinear Klein–Gordon equation in spatial dimensions three and higher are unstable.

In mathematics, a linear recurrence with constant coefficients sets equal to 0 a polynomial that is linear in the various iterates of a variable—that is, in the values of the elements of a sequence. The polynomial's linearity means that each of its terms has degree 0 or 1. A linear recurrence denotes the evolution of some variable over time, with the current time period or discrete moment in time denoted as t, one period earlier denoted as t − 1, one period later as t + 1, etc.

The Fokas method, or unified transform, is an algorithmic procedure for analysing boundary value problems for linear partial differential equations and for an important class of nonlinear PDEs belonging to the so-called integrable systems. It is named after Greek mathematician Athanassios S. Fokas.

References