Korteweg - de Vries equation


The Korteweg - de Vries equation is a partial differential equationMathworldPlanetmath defined as

ut+6uux+uxxx=0 (1)

where u=u(x,t) and the subscripts indicate derivatives. This equation arises in hydrodynamics and was originally proposed to model waves in a canal. In addition to its practical applications, this equation is quite interesting as an object of mathematical study. It exhibits interesting solitonMathworldPlanetmath solutions, has a large algebra of conserved quantities, and can be solved using methods of inverse scattering.

1 Travelling Wave Solution

It is easy to exhibit a solution which describes a traveling wave. To find this solution, one substitutes the ansatz u(x,t)=f(x-ct) 11The most general form of a solution to the two-dimensional wave equationMathworldPlanetmath is u(x,t)=f(x-ct)+f(x+ct) which describes waves propagating in both directions with velocity c. Here, we will look for a solution which only describes a wave propagatng in one direction. into the equation to obtain the following equation for f:

-cf+6ff+f′′′=0

This equation can be written as

(-cf+3f2+f′′)=0

and, hence,

-cf+3f2+f′′=k

for some constant k. Multiplying by f, we can repeat the same trick:

-cff+3f2f+ff′′-kf=(-cf2/2+f3+(f)2-kf)=0

hence

-cf2/2+f3+(f)2-kf=h

or

f=kf+cf2/2-f3

This can be solved implicitly by an integralDlmfPlanetmath:

y=f(y)kf+cf2/2-f3𝑑y+C

Since this is an elliptic integralMathworldPlanetmath, the result is an elliptic functionMathworldPlanetmath.

The solution obtained above is know as a solitary wave, or soliton. This term “solitary wave” refers to the fact that this solution describes a single wave pulse traveling with velocity c. Note that the amplitude of the pulse is determined by its velocity, unlike in the case of linear wave equations where the velocity of propagation does not depend upon the amplitude. There are also solutions which describe more than one solitary wave. In particular, there are solutions in which two of these waves collide and then re-emerge from the collision.

2 Conserved Currents

It is possible to exhibit integrals whose value is conserved under the evolution. To construct these quantities, we begin with vector fields. Suppose that u is a thrice differentiable function of x and t and consider the following vector fields:

v(1) =(u,3u2+uxx)
v(2) =(u2,4u3-ux2+2uuxx)

Computing their divergences and rearranging the result,

vt(1)t+vt(1)x =ut+6uux+uxxx
vt(2)t+vt(2)x =2u(ut+6uux+uxxx).

If u happens to satisfy our differential equation, then the divergence of these vector fields will go zero.

To obtain the conserved quantities, we will integrate them over a rectangle in the t,x plane with sides parallel to the x and t axes use Green’s theorem .

t1t2𝑑tx1x2𝑑x (ut(t,x)+6u(t,x)ux(t,x)+uxxx(t,x))
= x1x2𝑑x(3u2(x,t2)+uxx(x,t2))+t2t1𝑑tu(t,x2)+
x1x2𝑑x(3u2(x,t1)+uxx(x,t1))+t1t2𝑑tu(t,x1)
= 3x1x2𝑑xu2(x,t2)-3x1x2𝑑xu2(x,t1)
+ t1t2𝑑tu(t,x1)-t1t2𝑑tu(t,x2)+ux(x2,t2)-ux(x1,t2)-ux(x2,t1)+ux(x1,t1)
t1t2𝑑tx1x2𝑑x 2u(ut(t,x)+6u(t,x)ux(t,x)+uxxx(t,x))
= x1x2𝑑x(4u3(x,t2)-ux2(x,t2)+2u(x,t2)uxx(x,t2))+t2t1𝑑tu2(t,x2)+
x2x2𝑑x(4u3(x,t1)-ux2(x,t1)+2u(x,t1)uxx(x,t1))+t1t2𝑑tu2(t,x1)
= x1x2𝑑x(4u3(x,t2)-3ux2(x,t2))-x1x2𝑑x(4u3(x,t1)-3ux2(x,t1))
+ t1t2𝑑tu2(t,x1)-t1t2𝑑tu2(t,x2)
+ u(x2,t2)ux(x2,t2)-u(x1,t2)ux(x1,t2)-u(x2,t1)ux(x2,t1)+u(x1,t1)ux(x1,t1)

We consider now several popular boundary conditionsMathworldPlanetmath for our equation:

  • Periodic Suppose that u is periodic in x with period p and satisfies the Kortwieg - de Vries equation. Then we integrals over t cancel against each other as do the endpoint terms, leaving us with

    x1x2𝑑xu2(x,t1)=x1x2𝑑xu2(x,t2).
  • Whole line Suppose that u satisfies the Kortwieg - de Vries equation and that u and its derivatives tend towards zero as x±. Then, taking the limit as x1- and x2+, we conclude that

    -+𝑑xu2(x,t1)=-+𝑑xu2(x,t2).
  • Finite interval We may also consider our differential equation on a finite interval and impose suitable boundary conditions at the endpoints. In this case, the integrals with respect to t do not cancel or automatically vanish as a consequence of the boundary conditions, so they will not, in general, give conserved quantities. Nevertheless, it is at least still possible to compute their value solely from the boundary data without having to know how the solution behaves in the interior.

Starting with the two vector fields given above, it is possible to generate more such vector fields and more conserved quantities. Since the Lie bracket of two vector fields with zero divergence also has zero divergence, we may

Title Korteweg - de Vries equation
Canonical name KortewegDeVriesEquation
Date of creation 2013-03-22 13:43:20
Last modified on 2013-03-22 13:43:20
Owner rspuzio (6075)
Last modified by rspuzio (6075)
Numerical id 19
Author rspuzio (6075)
Entry type Definition
Classification msc 35Q53
Synonym KdV equation