Korteweg–De Vries equation explained

In mathematics, the Korteweg–De Vries (KdV) equation is a partial differential equation (PDE) which serves as a mathematical model of waves on shallow water surfaces. It is particularly notable as the prototypical example of an integrable PDE and exhibits many of the expected behaviors for an integrable PDE, such as a large number of explicit solutions, in particular soliton solutions, and an infinite number of conserved quantities, despite the nonlinearity which typically renders PDEs intractable. The KdV can be solved by the inverse scattering method (ISM). In fact, Gardner, Greene, Kruskal and Miura developed the classical inverse scattering method to solve the KdV equation.

The KdV equation was first introduced by and rediscovered by Diederik Korteweg and Gustav de Vries in 1895, who found the simplest solution, the one-soliton solution. Understanding of the equation and behavior of solutions was greatly advanced by the computer simulations of Zabusky and Kruskal in 1965 and then the development of the inverse scattering transform in 1967.

Definition

The KdV equation is a partial differential equation that models (spatially) one-dimensional nonlinear dispersive nondissipative waves described by a function

\phi(x,t) adhering to:\partial _{t} \phi + \partial ^{3} _{x} \phi - 6 \phi \partial _{x} \phi =0 x \in R , t \geq 0 ,

where

\partial ^{3} _{x} \phi accounts for dispersion and the nonlinear element \phi \partial _{x} \phi is an advection term.

For modelling shallow water waves,

\phi is the height displacement of the water surface from its equilibrium height.

The constant

6 in front of the last term is conventional but of no great significance: multiplying t , x , and \phi by constants can be used to make the coefficients of any of the three terms equal to any given non-zero constants.

Soliton solutions

One-soliton solution

Consider solutions in which a fixed wave form (given by

f(X) ) maintains its shape as it travels to the right at phase speed c . Such a solution is given by \varphi(x,t) = f(x - ct - a) = f(X) . Substituting it into the KdV equation gives the ordinary differential equation- c (df)/(+ dX) (d ^{3f})/(-6 f dX ^{3}) (df)/(dX) = 0 ,

or, integrating with respect to

X ,- cf + (d ^{2} f)/(dX ^{2}) -3 f ^{2} = A

where

A is a constant of integration. Interpreting the independent variable X above as a virtual time variable, this means f satisfies Newton's equation of motion of a particle of unit mass in a cubic potentialV(f) = - \left(f ^{3+ (1)/(2)} cf ^{2} + Af\right) .IfA = 0 ,c > 0 then the potential function V(f) has local maximum at f =0 ; there is a solution in which f(X) starts at this point at 'virtual time' -inf ty , eventually slides down to the local minimum, then back up the other side, reaching an equal height, and then reverses direction, ending up at the local maximum again at time inf ty . In other words, f(X) approaches 0 as X\to -inf ty . This is the characteristic shape of the solitary wave solution.

More precisely, the solution is

\phi(x,t) =- (1)/(2) c \operatorname{sech} ^{2\left[{\sqrt{c}\over} 2 }(x - ct - a)\right]

where

\operatorname{sech} stands for the hyperbolic secant and a is an arbitrary constant. This describes a right-moving soliton with velocity c .

N-soliton solution

There is a known expression for a solution which is an

N -soliton solution, which at late times resolves into N separate single solitons. The solution depends on an decreasing positive set of parameters \chi _{1,} … , \chi _{N} > 0 and a non-zero set of parameters \beta _{1,} … , \beta _{N} . The solution is given in the form$$\backslash phi(x,t)\; =\; -2\backslash frac\; \backslash mathrm[\backslash mathrm\{det\}\; A(x,t)]$$where the components of the matrix A(x,t) are given byA _{nm} (x,t) = \delta _{nm} + (\beta _{n} e ^{8 \chi ^{3 t} _{n}} e ^{- (\chi _{n} + \chi _{m)x}})/(\chi _{n} + \chi _{m}) .

This is derived using the inverse scattering method.

Integrals of motion

The KdV equation has infinitely many integrals of motion, which do not change with time. They can be given explicitly as

\int ^{+inf ty} _{-inf ty} P _{2n-1} (\phi, \partial _{x} \phi, \partial ^{2} _{x} \phi, \ldots) d x

where the polynomials

P _{n} are defined recursively by\begin{align} P _{1&=\phi,} \\ P _{n} & =- (dP _{n-1})/(dx) + \sum ^{n -2} _{i =1} P _{i} P _{n-1-i} for n \ge 2 . \end{align} The first few integrals of motion are:

• the mass

\int \phi d x,

• the momentum

\int \phi ^{2} d x,

• the energy

\int \left[ 2 \phi ^{3} - \left( \partial _{x} \phi \right) ^{2} \right] d x .Only the odd-numbered terms P _{2n+1} result in non-trivial (meaning non-zero) integrals of motion.

Lax pairs

The KdV equation

\partial _{t\phi} = 6 \phi \partial _{x} \phi - \partial ^{3} _{x} \phi can be reformulated as the Lax equationL _{t} = [L,A] \equiv LA - AL with L a Sturm–Liouville operator:\begin{align} L & = - \partial ^{2} _{x} + \phi, \\ A & = 4 \partial ^{3} _{x} - 6 \phi \partial _{x} - 3 [\partial _{x,} \phi] \end{align} where [\partial _{x,\phi]} is the commutator such that [\partial _{x,\phi]f=f\partial} _{x\phi} . The Lax pair accounts for the infinite number of first integrals of the KdV equation.

In fact,

L is the time-independent Schrödinger operator (disregarding constants) with potential \phi(x,t) . It can be shown that due to this Lax formulation that in fact the eigenvalues do not depend on t .

Zero-curvature representation

Setting the components of the Lax connection to bethe KdV equation is equivalent to the zero-curvature equation for the Lax connection,$$\backslash partial\_t\; L\_x\; -\; \backslash partial\_x\; L\_t\; +\; [L\_x,\; L\_t]\; =\; 0.$$

Least action principle

The Korteweg–De Vries equation

\partial _{t} \phi + 6 \phi \partial _{x} \phi + \partial ^{3} _{x} \phi = 0 , is the Euler–Lagrange equation of motion derived from the Lagrangian density, l{L}

with

\phi defined by\phi : = (\partial \psi)/(\partial x) .

Since the Lagrangian (eq (1)) contains second derivatives, the Euler–Lagrange equation of motion for this field is

where

\partial is a derivative with respect to the \mu component.

A sum over

\mu is implied so eq (2) really reads,

Evaluate the five terms of eq (3) by plugging in eq (1),

\partial _{tt} \left( (\partial l{L)/() } \right) = 0 \partial _{xx} \left( (\partial l{L)/() } \right) = \partial_ \left(-\partial_ \psi \right) \partial _{t} \left( (\partial l{L)/() } \right) = \partial_t \left(\frac \partial_x \psi \right) \partial _{x} \left( (\partial l{L)/() } \right) = \partial_x \left(\frac \partial_t \psi + 3 (\partial_x \psi)^2 \right) \,(\partial l{L)/() } = 0

Remember the definition

\phi = \partial _{x} \psi , so use that to simplify the above terms,\partial _{xx} \left( - \partial _{xx} \psi \right) = - \partial _{xxx} \phi \partial _{t} \left( (1)/(2) \partial _{x} \psi \right) = (1)/(2) \partial _{t} \phi \partial _{x} \left( (1)/(2) \partial _{t} \psi + 3 (\partial _{x} \psi) ^{2} \right) = (1)/(2) \partial _{t} \phi + 3 \partial _{x} (\phi) ^{2} = (1)/(2) \partial _{t} \phi + 6 \phi \partial _{x} \phi

Finally, plug these three non-zero terms back into eq (3) to see

\left( - \partial _{xxx} \phi \right) - \left( (1)/(2) \partial _{t} \phi \right) - \left( (1)/(2) \partial _{t} \phi + 6 \phi \partial _{x} \phi \right) = 0 , which is exactly the KdV equation\partial _{t} \phi + 6 \phi \partial _{x} \phi + \partial ^{3} _{x} \phi = 0 .

Long-time asymptotics

It can be shown that any sufficiently fast decaying smooth solution will eventually split into a finite superposition of solitons travelling to the right plus a decaying dispersive part travelling to the left. This was first observed by and can be rigorously proven using the nonlinear steepest descent analysis for oscillatory Riemann–Hilbert problems.

History

The history of the KdV equation started with experiments by John Scott Russell in 1834, followed by theoretical investigations by Lord Rayleigh and Joseph Boussinesq around 1870 and, finally, Korteweg and De Vries in 1895.

The KdV equation was not studied much after this until discovered numerically that its solutions seemed to decompose at large times into a collection of "solitons": well separated solitary waves. Moreover, the solitons seems to be almost unaffected in shape by passing through each other (though this could cause a change in their position). They also made the connection to earlier numerical experiments by Fermi, Pasta, Ulam, and Tsingou by showing that the KdV equation was the continuum limit of the FPUT system. Development of the analytic solution by means of the inverse scattering transform was done in 1967 by Gardner, Greene, Kruskal and Miura.

The KdV equation is now seen to be closely connected to Huygens' principle.

Applications and connections

The KdV equation has several connections to physical problems. In addition to being the governing equation of the string in the Fermi–Pasta–Ulam–Tsingou problem in the continuum limit, it approximately describes the evolution of long, one-dimensional waves in many physical settings, including:

• shallow-water waves with weakly non-linear restoring forces,
• long internal waves in a density-stratified ocean,
• ion acoustic waves in a plasma,
• acoustic waves on a crystal lattice.

The KdV equation can also be solved using the inverse scattering transform such as those applied to the non-linear Schrödinger equation.

KdV equation and the Gross–Pitaevskii equation

Considering the simplified solutions of the form

\phi(x,t) = \phi(x \pm t) we obtain the KdV equation as\pm \partial _{x} \phi + \partial ^{3} _{x} \phi + 6 \phi \partial _{x} \phi =0 or\pm \partial _{x} \phi + \partial _{x} (\partial ^{2} _{x} \phi + 3 \phi ^{2)} =0 Integrating and taking the special case in which the integration constant is zero, we have:- \partial ^{2} _{x} \phi - 3 \phi ^{2} = \pm \phi which is the λ=1 special case of the generalized stationary Gross–Pitaevskii equation (GPE)- \partial ^{2} _{x} \phi - 3 \phi ^{λ} \phi = \pm \phi Therefore, for the certain class of solutions of generalized GPE (λ=4 for the true one-dimensional condensate andλ=2 while using the three dimensional equation in one dimension), two equations are one. Furthermore, taking the λ=3 case with the minus sign and the \phi real, one obtains an attractive self-interaction that should yield a bright soliton.

Variations

Many different variations of the KdV equations have been studied. Some are listed in the following table.

Name	Equation
Korteweg–De Vries (KdV)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u + 6 u \partial _{x} u =0
KdV (cylindrical)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u - 6 u \partial _{x} u + \tfrac{ 1 }{ 2 t}u = 0
KdV (deformed)	\displaystyle \partial _{t} u + \partial _{x} \left ( (\partial ^{2} _{x} u - 2 η u ^{3} - 3 u (\partial _{x} u) ^{2})/(2 ( η+ u ^{2)}) \right ) = 0
KdV (generalized)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u = \partial ^{5} _{x} u
KdV (generalized)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u + \partial _{x} f(u) = 0
KdV (modified)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u \pm 6 u ^{2} \partial _{x} u = 0
Gardner equation	\displaystyle \partial _{t} u + \partial ^{3} _{x} u - ( 6 \varepsilon ^{2u} ^{2} + 6 u)\partial _{x} u = 0
KdV (modified modified)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u - \tfrac{ 1 }{ 8 }(\partial _{x} u) ^{3} + (\partial _{x} u)(Ae ^{au} + B + Ce ^{-au} ) = 0
KdV (spherical)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u - 6 u \partial _{x} u + \tfrac{ 1 }{t}u = 0
KdV (super)	\displaystyle \begin{cases} \partial _{t} u = 6 u \partial _{x} u - \partial ^{3} _{x} u + 3 w \partial ^{2} _{x} w \ \partial _{t} w = 3 (\partial _{x} u) w + 6 u \partial _{x} w - 4 \partial ^{3} _{x} w \end{cases}
KdV (transitional)	\displaystyle \partial _{t} u + \partial ^{3} _{x} u - 6 f(t) u \partial _{x} u = 0
KdV (variable coefficients)	\displaystyle \partial _{t} u + \beta t ^{n} \partial ^{3} _{x} u + \alpha t ^{nu} \partial _{x} u = 0
KdV-Burgers equation	\displaystyle \partial _{t} u + \mu \partial ^{3} _{x} u + u \partial _{x} u - \nu \partial ^{2} _{x} u = 0
non-homogeneous KdV	\partial _{t} u + \alpha u + \beta \partial _{x} u + \gamma \partial ^{2} _{x} u = Ai(x), u(x, 0 ) = f(x)

See also

• Advection-diffusion equation
• Benjamin–Bona–Mahony equation
• Boussinesq approximation (water waves)
• Cnoidal wave
• Dispersion (water waves)
• Dispersionless equation
• Fifth-order Korteweg–De Vries equation
• Kadomtsev–Petviashvili equation
• KdV hierarchy
• Modified KdV–Burgers equation
• Novikov–Veselov equation
• Schamel equation
• Ursell number
• Vector soliton

References

• Berest . Yuri Y. . Loutsenko . Igor M. . Huygens' Principle in Minkowski Spaces and Soliton Solutions of the Korteweg-de Vries Equation . Communications in Mathematical Physics . 190 . 1 . 1997 . 0010-3616 . 10.1007/s002200050235 . 113–132. solv-int/9704012 .
  • Chalub . Fabio A.C.C. . Zubelli . Jorge P. . Huygens’ principle for hyperbolic operators and integrable hierarchies . Physica D: Nonlinear Phenomena . 213 . 2 . 2006 . 10.1016/j.physd.2005.11.008 . 231–245.
• Book: Darrigol, Olivier . Worlds of Flow . Oxford University Press . Oxford ; New York . 2005 . 978-0-19-856843-8.
• Book: Dauxois, Thierry . Peyrard . Michel . Physics of Solitons . Cambridge University Press . Cambridge, UK ; New York . 2006 . 0-521-85421-0 . ocm61757137.
• Book: Dingemans, M. W. . Water Wave Propagation Over Uneven Bottoms . World Scientific . River Edge, NJ . 1997 . 981-02-0427-2.
• Book: Dunajski, Maciej . Solitons, Instantons, and Twistors . OUP Oxford . Oxford ; New York . 2009 . 978-0-19-857063-9 . 320199531.
• Gardner . Clifford S. . Greene . John M. . Kruskal . Martin D. . Miura . Robert M. . Method for Solving the Korteweg-deVries Equation . Physical Review Letters . 19 . 19 . 1967 . 0031-9007 . 10.1103/PhysRevLett.19.1095 . 1095–1097.
  • Korteweg . D. J. . de Vries . G. . XLI. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves . The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science . 39 . 240 . 1895 . 1941-5982 . 10.1080/14786449508620739 . 422–443.
• Lax . Peter D. . Integrals of nonlinear equations of evolution and solitary waves . Communications on Pure and Applied Mathematics . 21 . 5 . 1968 . 0010-3640 . 10.1002/cpa.3160210503 . 467–490.
  • Book: Polyanin, Andrei D. . Zaitsev . Valentin F. . Handbook of Nonlinear Partial Differential Equations . Chapman and Hall/CRC . Boca Raton, Fla . 2003 . 978-1-58488-355-5.
• Book: Vakakis, Alexander F. . Normal Modes and Localization in Nonlinear Systems . Springer Science & Business Media . Dordrecht ; Boston . 2002 . 978-0-7923-7010-9.
• Zabusky . N. J. . Kruskal . M. D. . Interaction of "Solitons" in a Collisionless Plasma and the Recurrence of Initial States . Physical Review Letters . 15 . 6 . 1965 . 0031-9007 . 10.1103/PhysRevLett.15.240 . 240–243.

External links

• Korteweg–De Vries equation at EqWorld: The World of Mathematical Equations.
• Korteweg–De Vries equation at NEQwiki, the nonlinear equations encyclopedia.
• Cylindrical Korteweg–De Vries equation at EqWorld: The World of Mathematical Equations.
• Modified Korteweg–De Vries equation at EqWorld: The World of Mathematical Equations.
• Modified Korteweg–De Vries equation at NEQwiki, the nonlinear equations encyclopedia.
  • Derivation of the Korteweg–De Vries equation for a narrow canal.
• Three Solitons Solution of KdV Equation – https://www.youtube.com/watch?v=H4rN3Wr4ctw
• Three Solitons (unstable) Solution of KdV Equation – https://www.youtube.com/watch?v=5z5SylS2QHE
• Mathematical aspects of equations of Korteweg–De Vries type are discussed on the Dispersive PDE Wiki.
• Solitons from the Korteweg–De Vries Equation by S. M. Blinder, The Wolfram Demonstrations Project.
• Solitons & Nonlinear Wave Equations