Problems and Solutions For Partial Differential Equations
Problems and Solutions For Partial Differential Equations
Problems and Solutions For Partial Differential Equations
for
Partial Differential Equations
by
Willi-Hans Steeb
International School for Scientific Computing
at
University of Johannesburg, South Africa
Yorick Hardy
Department of Mathematical Sciences
at
University of South Africa, South Africa
Contents
1 Linear Partial Differential Equations
12
56
Bibliography
68
Index
69
vi
Chapter 1
Problem 1.
equation
is given by
1 (x x0 + 2t)2
u(x, t) = exp
.
4t
4t
Problem 2.
mf
.
xm
(1)
(2)
Dxm (f f ) = 0,
(3)
Problem 3.
for m
odd
(1)
(2)
2mE
~2
with
G = 0.
Problem 5. Consider an electron of mass m confined to the x y plane
and a constant magnetic field B parallel to the z-axis, i.e.
0
B = 0 .
B
The Hamilton operator for this two-dimensional electron is given by
p + eA)2
1
= (
H
=
((
px + eAx )2 + (
py + eAy )2 )
2m
2m
where A is the vector potential with
B=A
and
,
py = i~ .
x
y
(i) Show that B can be obtained from
0
A = xB
0
px = i~
yB
A = 0 .
0
(ii) Use the second choice for A to find the Hamilton operator H.
(iii) Show that
px ] = 0.
[H,
(iv) Let k = px /~. Make the ansatz for the wave function
(x, y) = eikx (y)
= E reduces to
and show that the eigenvalue equation H
mc2
~2 d 2
2
+
(y y0 ) (y) = E(y)
2m dy 2
2
where
eB
~k
,
y0 :=
.
m
eB
(v) Show that the eigenvalues are given by
1
En = n +
~c ,
n = 0, 1, 2, . . . .
2
c :=
Problem 6.
1
i~
=
+ V (x) .
t
2m
:= ,
v := =
.
Problem 7.
where
c(x, t)
.
x
The diffusion coefficient D(x) depends as follows on x
j(x, t) = D(x)
D(x) = D0 (1 + gx)
g 0.
i = + V (x)
t
2
where
=
2
2
2
+
+
.
x21
x22
x23
Let
:= ||2 = ,
(1)
v := =
= =
x1
x2
x3
(x, t)
~2 2
=
(x, t) + V (x, t)(x, t).
t
2m
where the functions and S are real. Find the partial differential equations
are and S.
= E of a particle
Problem 10. Consider the Schrodinger equation H
on the torus. A torus surface can be parametrized by the azimuthal angle
and its polar angle
x(, ) = (R + a cos()) cos()
y(, ) = (R + a cos()) sin()
z(, ) = a sin()
where R and a are the outer and inner radius of the torus, respectively such
that the ratio a/R lies between zero and one.
(i) Find H.
(ii) Apply the separation ansatz
(, ) = exp(im)()
where m is an integer.
Problem 11.
t 0,
< x <
where
0
G(x, t|x , t ) = p
(x x0 )2
exp
4D(t t0 )
4D(t t0 )
1
.
Here G(x, t|x0 , t0 ) is called the fundamental solution of the diffusion equation
obtained for the initial data (x x0 ) at t = t0 , where denotes the Dirac
delta function.
(i) Let u(x, 0) = f (x) = exp(x2 /(2)). Find u(x, t).
(ii) Let u(x, 0) = f (x) = exp(|x|/). Find u(x, t).
Problem 12. Let f : R3 R3 be a differentiable function. Consider the
first order vector partial differential equation
f = kf
B=
1 E
,
c2 t
E=
B
,
t
divE = 0,
divB = 0
with B = 0 H.
(i) Assume that E2 = E3 = 0 and B1 = B3 = 0. Simplify Maxwells
equations.
(ii) Now assume that E1 and B2 only depends on x3 and t with
E1 (x3 , t) = f (t) sin(k3 x3 ),
where k3 is the third component of the wave vector k. Find the system of
ordinary differential equations for f (t) and g(t) and solve it and thus find
the dispersion relation. Note that
E3 /x2 E2 /x3
E = E1 /x3 E3 /x1 .
E2 /x1 E1 /x2
Problem 14.
tial equation
1
(cv)
k
Problem 15.
+
+
f = 0.
x1
x2
x3
,
x
Dt =
+ u(x, t) .
t
x
Cex
ex
on
[0, 1] [0, 1]
u(x, 1) = 2,
u(0, y) = 1,
,
x2 j,k
(x)2
2u
y 2
u(1, y) = 2.
j,k
and then solve the linear equation. Consider the cases x = y = 1/3 and
x = y = 1/4.
Problem 19.
+
+ k 2 = 0.
x2
y 2
v(x, y) = 2xy,
e
(x, y) = (u(x,
y), v(x, y)).
admits the solution u(t, ) = sin() sin(t). Does this solution satisfies the
boundary condition u(t, = ) = 0 and the initial condition u(t = 0, ) =
0?
Problem 22.
= ~
H
2
3
X
k
1 2
+ ((q2 q1 d)2 + (q3 q2 d)2 )
2
mj qj
2
j=1
where d is the distance between two adjacent atoms. Apply the linear
transformation
(q1 , q2 , q3 ) = (q2 q1 ) d
(q1 , q2 , q3 ) = (q3 q2 ) d
1
(m1 q1 + m2 q2 + m3 q3 )
X(q1 , q2 , q3 ) =
M
where M := m1 +m2 +m3 and show that the Hamilton operator decouples.
Problem 23.
2
2
= 1
H
+
+
6g(x1 x2 )(x2 x3 )
2m x21
x22
x23
2
p
+ V (x, t) (x, t) = i~
2m
t
= i~
2m x2
t
(x, t)
p2
(p, t) = i~
.
2m
t
Consider the momentum space approach with the solution
(p, t) = 0 (p) exp(ip2 t/2m~)
where (p, 0) = 0 (p) is the initial momentum distribution. We define
x
:= i~(/p). Find h
xi(t).
Problem 27.
Show that
2
2
1
x + (y 1)2
x + (y + 1)2
1
T (x, y, t) = exp
exp
4t
4t
t
t
satisfies this partial differential equation.
Problem 30.
= 2 +
(xu)
t
x
x
(xu) = 0.
2+
x
x
(ii) Find time dependent solution of the inital value problem with
u(t = 0, x) = N0 (x 0 ) =
2
1
e(x0 ) /(20 )
20
with 0 0 the variance and 0 the mean. With = 0 we have the delta
function at 0.
Problem 31. Consider the Schrodinger equation
N
2
X
X
11
describing a one-dimensional Bose gas with the -fnction repulsive interaction. Show that for N = 2 the eigenvalue problem can be solved with a
exponential ansatz.
Problem 32.
n n cos( nt)
n=0
2 P (t, )
= P (t, )
t2
P (0, )
= 0,
t
n n = (),
n=0
P (t, 0) = (0).
e
P (z, )
exp(z 2 ln(/x))
(x)
dx.
x
Chapter 2
Nonlinear Partial
Differential Equations
Problem 2.
v
,
t
x
y
M=
+u
.
y
x
x
c = constant
(1)
(2)
13
Problem 3. Show that the nonlinear nondispersive part of the Kortewegde Vries equation
u
u
+ ( + u)
=0
(1)
t
x
possesses shock wave solutions that are intrinsically implicit
u(x, t) = f (x ( + u(x, t))t)
(2)
(1)
where a1 , a2 and a3 are non-zero constants. It contains dispersive, dissipative and nonlinear terms.
(i) Find a solution of the form
u(x, t) =
b1
(1 + exp(b2 (x + b3 t + b4 )))2
(2)
(1)
1
u(x, t),
6
(t) = 5t.
(2)
v
2v
+ av(1 6v).
=
x2
(3)
(4)
j=0
Show that
u(x, y, t) =
1 + exp
!!2
(x y/ 2) (5/ 6)t
2
2
+
x2
2
2
2
2
= 0.
2
y
x y xy
(1)
(i) Show that this equation has an implicit solution for of the form
x = X(, y, t) = f (, t)y + h(, t)
(2)
L = (ut c2 u2x ) +
cos
u .
(1)
2
m
Here c is a limiting velocity while m, , and are real parameters. ut and
ux are the partial derivatives of with respect to t and x, respectively. In
the terminology of quantum field theory, m is the mass associated with the
normal modes of the linearized theory, while /m2 is a dimensionless, coupling constant that measures the strength of the interaction between these
normal modes. In classical theory m is proportional to the characteristic
frequency of these normal modes.
15
(i) Let
x
t
,
t ,
u mu
m
m
and set c = 1. Show that then the Lagrangian density becomes
x
L=
m4
((u2t u2x ) + 2 cos u)
2
(2)
(3)
m2 2
(u + u2x 2 cos u) + .
2 t
(4)
(ii) By choosing
m4
(5)
show that the minimum energy of the theory is made zero and (4) can be
written as
m4 2
(u + u2x + 2(1 cos u)).
(6)
H=
2 t
=
u
(t = 0) = b(x).
t
(2)
Equation (1) can describe processes which develop in time, since it is of the
hyperbolic type if 1 + (u/x)2 (u/t)2 > 0. Show that the hyperbolic
condition for (1) implies for the initial conditions that
1 + a02 (x) b2 (x) > 0.
Problem 10.
(3)
(1)
(2)
u(j)
.
x
(3)
(i) Let
u(0) (x, t) = ln(x ct).
(4)
(5)
(6)
(7)
4 2 c
a + 2a2 tan2 (a(x ct)).
3
6
(8)
+u
+
(u) = 0
t
x x
:=
2
2
2
+ 2 + 2.
2
x
y
z
(1)
(2)
Here u is the normalized deviation of the ion density from the average.
Exact solitonlike solutions exist in one, two and three space dimensions.
They depend on the independent variables through the combination
x ct,
respectively, where
u
u c
u=0
2
(4)
17
and
2
:=
,
x2
1
:=
1
:= 2
r r
r
r
2
(5)
for the three cases. (i) Show that for the one-dimensional case we obtain
the soliton solution
u(x, t) = 3c sech2 (c1/2 (x ct x0 )/2).
(6)
(ii) Show that the flat soliton (6) is unstable with respect to nonaligned
perturbations.
Problem 12.
w 2 w
+
+ 2|w|2 w = 0,
t
x2
= 1
(1)
= + h(),
(2)
:= x t
(3)
(1)
6
c1
+
2
2
5 x /Dx + y /Dy + z 2 /Dz + c2 (t + t0 )
24
x2 /Dx + y 2 /Dy + z 2 /Dz
(x2 /Dx + y 2 /Dy + z 2 /Dz + c2 (t t0 ))2
(2)
satisfies (1).
Problem 14. The Burgers equation determines the motion of a pressureless fluid subjected to dissipation
u
u
2u
+u
=
.
t
x
x2
(1)
(2)
2 v
.
v x
(3)
(ii) Solving (1) with respect to u we can introduce the iteration formula
u(j+1) =
(j)
(j)
ux
(u(j)
xx ut ),
j = 0, 1, 2, . . .
(4)
where ut u/t etc.. Show that for u(0) = v n we find the sequence
u(1) = (n 1)
1 v
v x
(5)
2 v
fixed point
(6)
v x
Remark. Thus the Hopf-Cole transformation is an attractor (fixed point)
of the iteration (4).
u(2) =
Problem 15.
u
=
u
.
t
x
x
2 !
x
1
x
1
for t > 0, < 1
x
x1
x1
1
u(x, t) =
0
for t > 0, > 1
x1
(1)
(2)
where
x1 := (6t)1/3
(3)
u + u = u2
(1)
x2
t
19
2
1
5
1 tanh
x vt
.
2 6
6
(2)
Problem 17. Show that the coupled system of partial differential equations
2n 2n
2
2E
E
= nE,
=
(|E|2 )
+
(1)
i
2
2
2
t
x
t
x
x2
admits the solution
is
s2
2
E(x, t) = E0 (x st) exp
x+i
t
(2a)
2
4
n(x, t) =
where
22
cosh2 (x st x0 )
(2b)
p
2(1 s)2
E0 (x st) =
cosh((x st x0 ))
(2c)
u1
u2
u3
+
+
=0
x
y
z
flow velocity field, and are the constant kinematic viscosity and fluid
density. Here p is the fluid pressure. We have
u1
u1
1
u1 x + u2 u
y + u3 z
u2
u2
2
u u := u1 u
x + u2 y + u3 z
u3
u3
3
u1 u
x + u2 y + u3 z
Find the time evolution of
u3
v := curlu u =
y
u1
z
u2
x
u2
z
u3
x
u1
y
Show that
u
u
+u
=0
t
x
(1)
=x
(2)
with
t(x, t) = t,
Problem 20.
x
(x, t) = u(x, t),
u
(
x(x, t), t(x, t)) = x.
(3)
(1)
(1)
1 u
,
u x
t(x, t) = t,
x
(x, t) = x
(2)
=
t
x
Problem 22.
u
u
2 +
.
x
(3)
(1)
u
(
x(x, t), t(x, t)) = u(x, t).
(2)
21
(1)
= 3
t
t
(2)
u(s, t)ds .
(3)
u
(
x(x, t), t(x, t)) = u(x, t) (1)
the equation
3u
u
= u3 3
t
x
can be recast into the differential equation
!
2
u
3u
1 u
1 u
2u
3 2
u
2 = 0.
t x
3
u
x
2 x
(2)
(3)
x, t)
1 v(
v(
x, t) x
(4)
1
u(s, t)
3 u(s, t)
2 u 1 u
=
ds =
u(s, t)
ds = u 2
.
2
t
t
s3
x 2 x
u (s, t)
(5)
Problem 25.
u
u
(
x, t) =
,
x
x
(x, t) = u(x, t),
t(x, t) = t
(1)
(2)
=u
2
t
x
Problem 26.
A(
x)
x
.
(3)
u
=0
t
(1)
u0
v
= 0.
t
(3)
u(x,t)
f (z)dz
u0
u
t
(4a)
!
Z u(x,t)
v
u
= exp
f (z)dz
(4b)
x1
x
1
u0
!
!
2
Z u(x,t)
Z u(x,t)
2v
u
2u
=
exp
exp
.
f
(z)dz
f
(z)dz
f
(u)
x21
x21
x1
u0
u0
(4c)
Problem 27.
It may be regarded as describing solitary waves in a two-dimensional Josephson junction. The Lamb substitution is given by
u(x, y, t) = 4 tan1 [M (x, y, t)] .
Find the equation for M .
(1)
23
(1)
2
2
2
2
2
2
x0
x1
xn
(2)
(3)
where w(x) and (w) are new unknown functions. Show that (1) takes the
form
d
d2
w +
w w = F ()
(4)
dw
dw2
where
2
2
2
w
w
w
w w :=
.
(5)
x0
x1
xn
(ii) Show that (4) is equivalent to the following equation
!
!
2
n
Pn
d
d
P
d
d2
w
+
+ 2 (w w )+
F () = 0 (6)
dw
Pn
dw
dw2
dw Pn
where Pn (w) is an arbitrary polynomial of degree n in w, and = 1, 0, 1.
Moreover Pn dPn /dw.
(iii) Assume that satisfies
!
d2 d Pn
+
= F ().
(7)
dw2
dw Pn
Show that (6) takes the form
d
dw
Pn
w
Pn
!
+
d2
(w w ) = 0.
dw2
(8)
Pn
,
Pn
w w =
(9)
is also a solution of (8), and in this way we obtain a solution of (1) provided
satisfies (7).
2u 2u
u + u3 = 0.
t2
x2
(i) Show that (1) can be derived from the Lagrangian density
L(ut , ux , u) =
1 2
(u u2x ) g(u)
2 t
(1)
(2)
1 2
1
(u 1)2 (u4 2u2 + 1).
4
4
(3)
1
X2 =
2
u
t
2
+
u
x
2 !
1
(u2 1)2 .
4
(6)
L
L
+
=0
t ut
x ux
u
(7)
u + u3 = 0.
t2
x2
(1)
25
(2)
where
2 (1 v 2 )1
(3)
and x0 is a constant. These solutions are of kink (K) and antikink (K)
type traveling at constant velocity v. These solutions do not tend to zero at
infinity, but they do connect two minimum states of the potential 14 (u2 1)2 .
(ii) The energy and momentum densities T1 and T2 for the kink and antikink
are obtained by substituting (2) into (5) of problem 1. Show that
T1 (x, t) =
2
cosh4 ((/ 2)(x vt x0 ))
2
(4)
2v
(5)
cosh4 ((/ 2)(x vt x0 )).
2
(iii) Show that by integrating over x we find the energy and the momentum
Z +
Z +
4
4 v
E :=
T1 dx =
,
P :=
T2 dx = .
(6)
3
3
2
2
T2 (x, t) =
(iv) Show that the densities T1 and T2 are localized in space, in contrast
with uK and uK .
(v) We associate the mass
M 2 := E 2 P 2
(7)
xT1 dx
Xc := Z
.
+
T1 dx
(8)
Show that the kink and antikink solutions of (6) take the values
M2 =
8
,
9
Xc = vt + x0 .
(9)
y(x, y) = y,
u
(
x(x, y), y(x, y)) = x.
(1)
=K 1+
+
(2)
2
2
x y
xy
x
y
is hodograph invariant.
Problem 32. Show that the nonlinear equation of the Born-Infeld type
for a scalar field u is obtained by varying the Lagrangian
L(ux , uy , uz , ut ) = 1 (1 + u2x + u2y + u2z u2t )1/2
(1)
x
y
z
t
x2
y
z
t
u
x
2
2u
x2
u
y
2
2u
y 2
u
z
2
2u
z 2
u
t
2
2u 2u
2u
2
2
t
xy xy
u u 2 u
u u 2 u
u u 2 u
u u 2 u
u u 2 u
2
+2
+2
+2
= 0.
x x xz
y z yz
t x tx
t y ty
t z tz
(2)
1 2
m4
2 2
L(ut , ux , u) = (ut c ux ) +
cos
u .
(1)
2
m
Here c is a limiting velocity while m, , and are real parameters. ut and
ux are the partial derivatives of u with respect to t and x, respectively. In
the terminology of quantum field theory, m is the mass associated with the
normal modes of the linearized theory, while /m2 is a dimensionless, coupling constant that measures the strength of the interaction between these
normal modes. In classical theory m is proportional to the characteristic
frequency of these normal modes. Let
x
x
,
m
t
,
m
mu
u
(2)
and set c = 1. (i) Show that then the Lagrangian density becomes
L(ut , ux , u) =
m4
((u2t u2x ) + 2 cos u)
2
(3)
27
H=
(4)
(5)
the minimum energy of the theory is made zero and (4) can be written as
H=
m4 2
(u + u2x + 2(1 cos u)).
2 t
(6)
Problem 34. Show that the following theorem holds. The conservation
law
(2)
(3)
1
T (u/x, u/t, u)
(4a)
F (u/x, u/t, u)
.
T (/x; /t; u)
(4b)
T 0 (u/x0 , u/t0 , u) =
F 0 (u/x0 , /t0 , u) =
t0 (x, t) = t,
1
= 0 0,
x
T x
F0
= 0 0 + 0.
t
T x
t
(5)
2w
2w
2
+ 2 + 2 (u u) = 0
2
x
y
x
(1b)
x
,
t
y(x, y, t) =
y
,
t
1
t(x, y, t) = ,
t
(2a)
ix2
iy 2
+
4t
4t
u
(
x(x, y, t), y(x, y, t), t(x, y, t))
(2b)
1
w(
x(x, y, t), y(x, y, t), t(x, y, t)).
(2c)
t2
Show that u
and w
satisfy the same equation with the subscripts x, y, t
replaced by x
, y, t.
w(x, y, t) =
Problem 36. An simplified analog of the Boltzmann equation is constructed as follows. It is one-dimensional and the velocities of the molecules
are allowed to take two discrete values, c, only. Thus the distribution
function in the Boltzmann equation, f (x, v, t), is replaced by two functions
u+ (x, t) and u (x, t) denoting the density of particles with velocity +c or
c, repsectively, at point x and time t. The gas is not confined, but x varies
over all points of the real line (, ). It is further assumed that there
are only two types of interaction, viz., two + particles go over into two
particles, and vice versa, the probalitity for both processes to occur within
one unit of time being the same number . Then the Boltzmann equation,
which in absence of external forces reads
Z
f (x, v, t)
+vf (x, v, t) = d3 v1 d|vv1 |(|vv1 |, )[f 0 f10 f f1 ] (1)
t
where f 0 is the final distribution, i.e., the distribution after a collision.
Under the assumption described above (1) translates into a system of two
equations
u+
u+
+c
= (u2 u2+ ),
t
x
u
u
c
= (u2+ u2 )
t
x
(2)
where c and are positive constants. This model is called the Carleman
model. The Carleman model is rather unphysical. However with its aid one
can prove almost all those results which one would like to obtain for the
Boltzmann equation itself as, for instance, the existence of solutions for a
wide class of initial conditions or a rigorous treatment of the hydrodynamic
limit. (i) Show that as for the Boltzmann equation, the H theorem holds
for the Caleman model: The quantity
Z
(u+ (x, t) ln u+ (x, t) + u (x, t) ln u (x, t))dx
(3)
never decreases in time. (ii) Show that there exists the following generalizations of the H theorem. Let f be concave function, defined on the half-line
(0, ) which is once continuously differentiable. Let
Z
Sf := (f (u+ ) + f (u ))dx
(4)
29
Show that
d
Sf (u) 0.
(5)
dt
Thus, not only does entropy never decrease, but the same is true for all
quasientropies.
(iii) Show that as a consequence, all Renyi entropies
Z
S := (1 )1 ln (u
(6)
+ (x, t) + u (x, t))dx
never decrease. In information-theretical language, this means that all sensible measures of the lack of information are nondecreasing. In other words,
information is lost, or chaos is approached, in the strongest possible way.
Problem 37. For the functions u, v : [0, ) R R we consider the
Cauchy problem
u u
+
= v 2 u2 ,
t
x
u(0, x) = u0 (x),
v
v
= u2 v 2
t
x
(1a)
v(0, x) = v0 (x).
(1b)
(2)
D S
+
= 2DS
t
x
(3)
and S 2 D2 0.
(4)
(ii) Find explict solutions assuming that S and D are conjugate harmonic
functions.
Problem 38. The Carleman model as an approximation of the Boltzmann equation is given by the non-linear equations
u+
u+
+c
= (u2 u2+ ),
t
x
u
u
c
= (u2+ u2 ).
t
x
(i) Show that if u (x, 0) 0, then u (x, t) 0 for all t > 0.
(1a)
(1b)
Z
(u+ + u )dx = 0.
(2)
The two properties have to hold, of course, in any model that is considered
to be of some physical relevance.
(iii) Let
Z
H
(u+ ln u+ + u ln u )dx
(3)
(4)
(u+ ln u+ + u ln u ) + c (u+ ln u+ u ln u ) 0.
t
x
(5)
(f (u+ ) + f (u ))dx
(6)
= (f02 f+ f ).
t
x
(1a)
(1b)
(1c)
31
u := f+ f
(2)
+
(u) = 0,
t
x
Problem 40.
(u)
+
(f+ + f ) = 0.
t
x
(3)
u
+
f (u) = 0
t
x
(1)
(u) +
q(u) = 0
t
x
(2)
holds identically for any smooth vector field u(x, t) which staisfies (1).
(i) Show that (2) follows from (1) if
m
X
fj
q
=
,
u
u
u
j
k
k
j=1
k = 1, . . . , m.
(3)
(ii) Show that for m = 1, every convex function (n) is an entropy for (1)
with entropy flux
Z u
q(u) =
0 ()df ().
(4)
0
Problem 41.
y x
(1)
(2)
Z
c
Z
f (x)dx +
f (x) := lim
0
f (x)dx .
(3)
c+
The Cauchy principal value can be found by applying the residue theorem.
Let
f (x, t)
u(x, t) := i ln
(4)
f (x, t)
where
f (x, t) :=
N
Y
(x xj (t))
(5a)
j=1
j = 1, 2, . . . , N,
xn 6= xm
for n 6= m.
(5b)
u
= (ln(f f ))
(6)
H
t
t
follows from (2), (4) and (5).
(ii) Show that (1) is transformed into the form
N
X
j=1
N
dxj
1
dxj X
1
1
+
=
QN
j=1 (x xj (t))
QN
j=1 (x xj (t))
QN
j=1 (x xj (t))
QN
j=1 (x xj (t))
(7)
(iii) Show that (7) is equivalent to the equation
1
(f f ) = (f 2 f 2 ) = =(f 2 ).
t
2i
(8)
(iv) Show that by multiplying both sides of (7) by x xn (t) and then
putting x = xn , we obtain
QN
dxn
1
j=1 (xn (t) xj (t))
= QN
dt
2i j=1(j6=n) (xn (t) xj (t))
Problem 42.
n = 1, 2, . . . , N.
(9)
(1)
!
.
33
(2)
into (1) and require that w(z) satisfies an ordinary differential equation.
This is the so-called direct method and z is the so-called reduced variable.
Problem 43.
3
X
=0
M + F ()
= 0.
x
(1)
: R4 C4 ,
1
2
=
3
4
(2)
M is a positive constant,
:= ( 0 , ) 1 + 2 3 4
1
2
3
4
(3)
3 :=
1
0
0
1
(5)
(6)
u M u + 0 u + F (
uu)u = 0.
(7)
i
k
xk
k=1
(ii) Let
0
v(r)
1
.
u(x) :=
cos
iw(r)
i
sin e
(8)
(9a)
dv
= w(F (v 2 w2 ) (M + )).
dr
(9b)
(10)
u = v 2 w2 .
u
(11)
Problem 44. Consider a one-dimensional system to describe the electronbeam plasma system
ne
+
(ne (ue + V )) = 0
t
x
(1a)
ue
e
ue
1 pe
+ (ue + V )
=
E
t
x
me
me ne x
(1b)
E
E
+ (ue + V )
= 4eni (ue + V ).
(1c)
t
x
where ne is the density of the beam-electron fluids, me the beam-electron
mass, ue the bulk fluid velocity, p = Kb Te ne the particle fluid pressure, e the
charge on an electron and E = . We study the simplest case of uniform
plasma in the absence of an external electromagnetic field, and assume
that the positive ions are taken to form a fixed, neutralizing background
of uniform density ni = N0 = const throughout the present analysis. The
electrons move with a beam drift velocity V corresponding to the ions.
Assume that electrostatic perturbation is sinusoidal
n(t)
1
ne (x, t) = N0 1 +
(sin(kx) + cos(kx)) ,
n(t) N0
(2a)
N0
2
ue (x, t) = u(t)(cos(kx) sin(2kx)),
Define
X(t) := n(t),
Y (t) :=
1
ku(t),
2
Z(t) :=
e k
E(t),
me p 2
q :=
Kb Te k 2
,
m e N0 2
p :=
(3)
4e2
.
me
35
Show that when (2) and (3) are substituted into (1a)-(1c), we obtain
dX
= kV X XY + 2N0 Y
dt
dY
= qX + Y 2 pZ
dt
dZ
= N0 Y Y Z + kV Z.
dt
Problem 45.
(1)
subject to periodic boundary conditions in the interval [0, L], with initial
conditions u(x, 0) = u0 (x). We only consider solutions with zero spatial
average. We recall that for L 2 all initial conditions evolve into u(x, t) =
0. We expand the solution for u in the Fourier series
u(x, t) =
an (t) exp(ikn x)
(2)
n=
(3)
Here a
denotes the complex conjugate of a. Since we choose solutions with
zero average we have a0 = 0. (i) Show that inserting the series expansion
(2) into (1) we obtain the following system for the time evolution of the
Fourier amplitudes
X
dan
1
+ (kn4 kn2 ikn3 )an + ikn
(am anm + a
m an+m ) = 0.
dt
2
m=0
(4)
(ii) Show that keeping only the first five modes we obtain the system
da1
+ (1 ik 3 )a1 + ik(
a1 a2 + a
2 a3 + a
3 a4 + a
4 a5 ) = 0
dt
da2
+ (2 8ik 3 )a2 + ik(a21 + 2
a1 a3 + 2
a2 a4 + 2
a3 a5 ) = 0
dt
da3
+ (3 27ik 3 )a3 + 3ik(a1 a2 + a
1 a4 + a
2 a5 ) = 0
dt
da4
+ (4 64ik 3 )a4 + 2ik(a22 + 2a1 a3 + 2
a1 a5 ) = 0
dt
(5a)
(5b)
(5c)
(5d)
k := 2 ,
L
n := kn4 kn2 .
(5e)
(6)
Problem 46.
4 4
=
+ 4 (4)
t
x
x z
z x
"
2 2 #
2R
=R
+4
+
+
+
t
x
x z
z x
Ti
z
x
z
(1)
where is the Prandtl number, R the Rayleigh number, the stream
function, and a function measuring the difference between the profile of
temperature and a profile linearly decreasing with height and time, namely
= T T0 +
2
F t + Rz.
Ti
(2)
T0 T1 1
.
T1 r
(3)
This quantity is a constant once the fluid and the experimental setting are
chosen. The dimensionless temperatures appearing in the last term of the
second of equations (1) are now expressed by
T0 = Rc
(r + 1)
,
T1 =
Rc
.
(4)
=0
z
Show that the equation
dX
= X + Y,
dt
at z = 0, H
at z = 0, H
at x = kL,
kZ
dY
= rX Y XZ + W
dt
(5)
Z
1+
r
Z
1+
r
,
dZ
= bZ + XY
dt
37
(6)
can be derived from (1) and the boundary condition (5) via a Fourier expansion of the function (x, z, t) and (x, z, t).
Problem 47.
(1)
(2)
one obtains two real equations which, after suitable linear combination and
division by c21 , can be written as
!
2 !
2
2
2
2 R
2
+R
= (1 + )
+ + ( + )R R (3a)
t
t
x2
x
= (1 + 2 )
r2
(1 + (1 )R2 )R2 . (3b)
R 2 + 2 R 2
2 t
t
x
x
Here we have introduced c2 := c1 and := 1/c1 . (ii) Make an expansion
in of the form
R := R0 + 2 R2 + ,
Problem 48.
equation,
:= 1 (1 + 2 1 + )
(4)
w 2 w
+
+ Qw|w|2 = 0
(1)
t
x2
where Q is a constant. (i) Show that a discretization with the periodic
boundary conditions wj+N wj is given by
i
k = 1, 2, . . .
(2)
where
(a)
(1)
wj
:= wj
and
(2)
(b) wj
:=
1
(wj+1 + wj1 ).
2
(3)
(ii) Show that both schemes are of second-order accuracy. (iii) Show that
in case (2a) there are first integrals, the L2 norm,
I :=
N
1
X
j=0
|wj |2
(4)
N
1
X
j=0
1
|wj+1 wj |2
4
Q|wj | .
h2
2
(5)
The Poisson brackets are the standard ones. Thus when N = 2 the system
is integrable. This system has been used as a model for a nonlinear dimer.
(iv) Show that the Hamilton function of scheme (2b) is given by (h = 1)
H = i
N
1
X
j=0
1
4
ln(1 + Qwj wj )
Q
2
(6)
(7)
{qm , qn } = {pm , pn } = 0.
(8)
and
Problem 49.
where
:=
(1)
+ 2+ 2 2
2
x
y
z
t
(2)
1
0
0
1
,
1 :=
0
1
1
0
,
2 :=
0
i
i
0
,
3 :=
1 0
0 1
(4)
+ i1
+ i3
+ 2 .
x
y
z
t
(5)
Let a
denote complex conjugates. Let 0 2, 0 2, < <
and < < be arbitrary parameters. We set
U := exp(i1 exp(i2 e 1 )).
(6)
39
Assume that and are solutions of (1) and (2), respectively. Let
u := i.
(7)
is then given by
The B
acklund transformation B
i = B(,
, )
(8)
where B(,
, ) is the B
acklund transformation operator. The functions
and are related by
1
1
au = sin
u
U.
(9)
2
2
The B
acklund transformation works as follows. Let (respectively ) be a
solution of (1) (respectively 3)) then solve (9) for (respectively ). The
solution then solves (3) (respectively (1)).
(i) Show for any U such that
= U 1
U
(10)
all solutions of (9) must be of the form (i.e. plane travelling wave with
speed v less than one)
u(x, y, z, t) = f (),
:= kx + ly + mz t
(11)
v=
k2
2
.
+ l2 + m2
(ii) Thus show that we have a Backlund transformation between the ordinary differential equations
d2
= sin ,
d 2
d2
= sinh
d 2
defined by
du
= 2ei sin
d
1
u
.
2
u
u 2 u
2u
4
+ 6u
+
3 2 = 0.
(1)
2
x
t
x x
y
K(x, x; y, t)
x
(2)
where
Z
K(x, z; y, t) + F (x, z; y, t) +
(3)
2F
F
2F
+
=0
x2
z 2
y
(4)
(5)
(6)
and
Problem 51. Show that Hirotas operators Dxn (f g) and Dxm (f g) given
in (1) and (2) can be written as
Dxn (f
Dxm Dtn (f
g) =
n
X
(1)(nj) n! j f nj g
g) =
,
j!(n j)! xj xnj
j=0
n
m X
X
(1)(m+nji) m!
j=0 i=0
j!(m j)!
(3)
n!
i+j f n+mij g
.
i!(n i)! ti xj tni xmj
(4)
Problem 52.
(1)
2 ln f (x, t)
.
x2
(2)
41
+f 4 4
xt x t
x
x x3
2f
x2
2
= 0.
(3)
Problem 53.
(4)
u
3u 5u
u
u 2 u
+
15u
+
= 0.
+ 45u2
+ 15
t
x
x x2
x3
x5
(1)
2 ln f (x, t)
x2
(2)
Problem 54.
(3)
(1)
2 ln f (x, t)
x2
(2)
Problem 55.
(3)
2
3
H3 = w,
= .
Problem 56.
form
6=
u
=u
+ C
t
x
x
=
(u)
t x
Problem 57.
2
(+1)/2 (x, t).
+1
6= 1.
u
u 3 u
6u
+
= 0.
t
x x3
For a steady-state pulse solution we make the ansatz
c
b
2
(x ct)
u(x, t) = sech
2
2
where
1
.
cosh(y)
Find the condition on b and c such that this ansatz is a solution of the
Korteweg-de Vries equation.
sech(y) :=
Problem 58. Consider the Korteweg-de Vries equation and its solution
given in the previous problem. Show that
Z p
|u(x, t)|dx = .
43
Hint. We have
Z
Z
sech(s)ds
Problem 59.
1
ds = 2 arctan(es ).
cosh(s)
Show that
u(x, t) = 12
u
+u
+
=0
t
x
x
p
u
p
+ p
+u
=0
t
x
x
where u(x, t) is the velocity field, (x, t) is the density field and p(x, t) is
the pressure field. Here t is the time, x is the space coordinate and is the
ratio of specific heats. Find the linearized equation around u
e, pe, e.
Problem 61. Consider the classical it Heisenberg ferromagnetic equation
S
2S
=S
t
x2
where S = (S1 , S2 , S3 )T , S12 +S22 +S32 = 1 and denotes the vector product
The natural boundary conditions are S(x, t) (0, 0, 1) as |x| .
(i) Find partial differential equation under the stereographic projection
S1 =
2u
,
Q
S2 =
where Q = 1 + u2 + v 2 .
(ii) Perform a Painleve test.
2v
,
Q
S3 =
1 + u2 + v 2
Q
t
x2
Both systems of differential equations arise as consistency conditions of a
system of linear partial differential equations
= U ,
x
=V
t
+ [U, V ] = 0.
t
x
Two systems of nonlinear partial differential equations that are integrable
if there is an invertible 2 2 matrix g which depends on x and t such that
U1 = gU2 g 1 +
g 1
g ,
x
V1 = gV2 g 1 +
g 1
g .
t
v (Ju) = 0
t
x
where u = (u1 , u2 , u3 )T , v = (v1 , v2 , v3 ), u2 = v2 = 1, J = diag(j1 , j2 , j3 )
is a 3 3 diagonal matrix and denotes the vector product. Consider the
linear mapping M : R3 so(3)
0
u3 u2
M (u) = u3
0
u1
u2 u1
0
where so(3) is the simple Lie algebra of the 3 3 skew-symmetric matrices.
Rewrite the system of partial differential equations using M (u).
Problem 63.
45
0
u3 u2
M (u) = u3
0
u1
u2 u1
0
where so(3) is the simple Lie algebra of the 3 3 skew-symmetric matrices.
Rewrite the system using M (u).
Problem 64.
Let
u(x, t) =
t
.
x
u
+
= ,
x 2
=
+u
t
x2
x
Problem 66.
2
+u
,
2
x
x t
K :=
u
+ .
x
2
1 u 1 u
D +
2 x
2
,
x
D1 f (x) :=
f (s)ds.
(i) Show that applying the recursion operator R to the right-hand side of
the Burgers equation results in the partial differential equation
3u 3 2u 3 2 2u
u
+ u 2+ u
=
.
t
x3
2 x
4 x
(ii) Show that this partial differential equation can also be derived from the
linear partial differential equation
3
=
t
x3
and the transformation
(x, t) = exp
Note that
D
Z
1
2
u
x
u(s, t)ds .
= u.
Problem 67.
tions
Problem 68.
f
+ f 2 = 0,
x
g
+ f g = 0,
x
f
+ fg = 0
t
g
+ g 2 = 0.
t
u 3 u
u
6u
+
= 0.
t
x x3
Show that
2
log (x, t)
x2
is a solution of the Kortweg-de Vries equation, where
cj ck
3
3
(x, t) = det jk +
exp((j + k )t (j + k )x)
.
j + k
j,k=1,...,N
u(x, t) = 2
47
u
x
2
= 0.
Let v(x, t) = exp(u(x, t)). Find the partial differential equation for v.
Problem 70.
tion
1 2
+ ||2 = 0.
+
t
2 x2
Show that
(x, t) = 2sech(s) exp(i(s, t))
with
s := 2(x (t)),
(t) = 2t + 0 ,
(s, t) :=
s + (t)
(t) = 4(2 + 2 )t + 0 .
2u 2u
u + u(u2 + v 2 ) = 0,
t2 x2
(i) Show that
2v 2v
v + v(u2 + v 2 ) = 0.
t2 x2
2
v(x, t) = 0
u 2 u
+
+ 2|u|2 u = 0,
t
x2
< x <
w
= Mw
x2
with x1 = x, x2 = it and
L=
3
X
z S
=1
M =i
3
X
z S
,,=1
3
X
S
+ 2z1 z2 z3
z1 S
x
=1
provide a Lax pair. Here 1 , 2 , 3 are the Pauli spin matrices and the spectral parameters (z1 , z2 , z3 ) constitue an algebraic coordinate of an elliptic
curve defined by
z2 z2 =
Problem 74.
1
(J J ),
4
, = 1, 2, 3.
49
v(t, x)
Rs
.
w(t) + 0 dsv(s, t)
u 2 u
+ 2c|u|2 = 0
+
t
x2
where x R and c = 1. Show that it admits the solution
1 x 2
A
ln(t)
u(x, t) = exp it
+ 2cA2
+
.
4 t
t
t
t
i
= 0.
x t
x x3
y 2
Consider the one-soliton solution
u(x, y, t) = 2k12 sech2 (k1 x + k2 y t)
where sech(x) 1/ cosh(x) 2/(ex + ex ). Find the dispersion relation
(k1 , k2 ).
Problem 77. Find the partial differential equation given by the condition
u
u/z1
det
= 0.
u/z2 2 u/z1 z2
Find a solution of the partial differential equation.
Problem 78. A one-dimensional Schrodinger equation with cubic nonlinearity is given by
i~
~2 2
=
g,
t
2m x2
:=
where g > 0.
(i) Show that the partial differential equation admits the (soliton) solution
mv
~
m2 v 2
2u
s (x, t) = exp i
(x ut)
, 2 =
1
.
~
gm cosh((x vt))
~2
v
u
u(x0 , y) 0
M
+ 4 3 + 6u
+3
dx + = 0.
t
x
x
x
y
Problem 80.
u
x
2
g = g11 (x1 , x2 )dx1 dx1 +g12 (x1 , x2 )dx1 dx2 +g21 (x1 , x2 )dx2 dx1 +g22 (x1 , x2 )dx2 dx2
with g12 = g21 and the gjk are smooth functions of x1 , x2 . Let det(g)
g11 g22 g12 g21 6= 0 and
g11
g22
g12
2
R=
det g11 /x1 g22 /x1 g12 /x1 .
(det(g))2
g11 /x2 g22 /x2 g12 /x2
51
w 2 w
+ 2|w|2 w = 0.
+
t
x2
j1 0 0
S
2S
=S
+ S (JS), J = 0 j2 0
t
x2
0 0 j3
where S = (S1 , S2 , S3 )T and S12 +S22 +S32 = 1. Express the partial differental
equation using p(x, t) and p(x, t) given by
p
p
S1 = 1 p2 cos(q), S2 = 1 p2 sin(q), S3 = p.
Problem 84.
= sin(u).
x21
x22
c2 t2
Let
u(x1 , x2 , t) = 4 arctan(v(x1 , x2 , t)).
Find the partial differential equation for v. Separate this partial differential equation into a linear part and nonlinear part. Solve these partial
differential equations to find solutions for the two-dimensional sine-Gordon
equation. Note that
sin(4) 4 sin() cos() 8 sin3 () cos()
sin(arctan()) =
,
1 + 2
cos(arctan()) =
and therefore
sin(4 arctan(v)) =
1
1 + 2
4v(1 v 2 )
.
(1 + v 2 )2
Furthermore
2
2v
arctan(v) =
x21
(1 + v 2 )2
v
x1
2
+
1 2v
.
1 + v 2 x21
2u
u
u u u
+a
+b
+
=0
xt
x
t
x t
where a, b are constants. Show that the equation can be linearized with the
transformation
u(x, t) = bx at + ln(v(x, t).
Problem 88.
L
2
L
L
+
= 0.
t vt
x vx
x2 vxx
v
53
Problem 89.
L
1
L = (vx )3 (vxx )2 .
vt
2
c
(x ct) .
2
p
dx |u(x, t)| = .
Problem 90.
Show that
u+ (x, t) = 4 arctan(exp(d1 (x vt X) cosh())
u (x, t) = 4 arctan(exp(d1 (x vt X) cosh())
are solutions of the one-dimensional sine-Gordon equation, where v =
c tanh(), d = c/. Discuss.
Problem 91. Let (s, t) be the curvature and (s, t) be the torsion with
s and t being the arclength and time, respectively. Consider the complex
valued function
Z s
w(s, t) = (s, t) exp(i
ds0 (s0 , t)).
0
2
r
r 2 r = b
t
s
s
where b is the binormal unit vector, then w(s, t) satisfies the nonlinear
Schr
odinger equation
i
2
1
w + 2 w + |w|2 w = 0.
t
s
2
u
u
+A
=0
t
x
where u = (u1 , u2 , u3 )T and the 3 3 matrix is given by
u1
u2
u3
u2 u1
0
u3
0
u1
where and are constants. Along a characteristic curve C : x = x(s), t =
t(s) for the system of partial differential equation one has
det(A I3 ) = 0,
x0 (s)
.
t0 (s)
= Z2 V W
t
x
Z
1
1
= Z2 + V W
t
2
2
(i) Show that in the limit 0 the Navier-Stokes equation are invariant
under the scaling transformation ( > 0)
r r,
v h v,
t 1h t.
(ii) Show that for finite one finds invariance of the Navier-Stokes equation
if 1+h .
Problem 95.
k 2 /2
,
(cosh((kx t)/2))2
k3 =
2
ln(F (x, t)).
x2
55
Chapter 3
Problem 1.
(1)
x
2 (x1 , x2 , ) = x1 + x2
u
(
x1 (x1 , x2 ), x
2 (x1 , x2 ), ) =
u(x1 , x2 )
1 u(x1 , x2 )
t0 = c3 t,
u(x, t) = c2 u0 (x0 , t0 )
(2a)
(2b)
57
, (u) :=
.
x0 x1
xn
x0
x1
xn
Show that the equation is invariant under the Poincare algebra
,
,...,
x0 x1
xn
x0
+ x1
, x0
+ x2
, . . . , x0
+ xn
x1
x0
x2
x0
xn
x0
xj
xk
,
xk
xk
j, k = 1, 2, . . . , n
j 6= k .
,
x
Z1 =
Z2 =
,
t
Z3 =
+t ,
Z5 = u
x
t
x
x
u
Z6 = u + t ,
Z7 = x
+t +u
.
t
u
x
t
u
Z4 = x
Problem 5.
(2)
(1)
+ (v ) + v = 0
t
(1a)
v
1
+ (v )v + (p + H2 ) (H )H = 0
t
2
H
+ (v )H + H v (H )v = 0
t
H=0
p
+ (v )
=0
(1b)
(1c)
(1d)
(1e)
v1
= ia1 v2 v3 ,
x
v3
= ia3 v1 v2 .
t
(1)
similarity variable
(2)
and
v1 (x, t) = t
1+
2
v2 (x, t) = t
1+
2
1 ialnx
e
3 (s) exp(i3 (s))
x
where i 0, 0 i < 2 and = 1.
v3 (x, t) =
(3)
Problem 8.
The motion of an inviscid, incompressible-ideal fluid is
goverened by the system of equations
u
+ u u = p
t
(1)
u=0
(2)
59
(i) Show that the Lie symmetry group of the Euler equations in three
dimensions is generated by the vector fields
a
=a
+ a0
a00 x ,
x
u
p
vc = c
s1 = x
+ c0
c00 z ,
z
w
p
+y
+z
+t ,
x
y
z
t
r2 = y
x
+v
u ,
x
y
u
v
rx = z
vb = b
s2 = t
+ b0
b00 y
y
v
p
v0 =
u
v
w
2p
t
u
v
w
p
ry = z
x
+w
u
x
z
u
w
y
+w
v
,
y
z
v
w
vq = q
(4)
1 u
=
.
t
x u2 x
(1)
Within the jet bundle formalism we consider instead of (1) the submanifold
u2x
uxx
+
2
=0
u2
u3
(2)
uxxx
6ux uxx
6u3x
+
= 0, . . .
u2
u3
u4
(3)
(4)
F (x, t, u, ux , uxx ) ut
and its differential consequences
Fx (x, t, u, ux , ) utx
with the contact forms
= du ux dx ut dt,
The non-linear partial differential equation (1) admits the Lie point symmetries (infinitesimal generator)
Xv = ux
,
u
Tv = ut
,
Vv = (xux + u) .
(5)
u
u
The subscript v denotes the vertical vector fields. The non-linear partial
differential equation (1) also admits Lie-B
acklund vector fields. The first in
the hierarchy is
uxxx
9ux uxx
12u3x
.
(6)
U=
+ 5
u3
u4
u
u
Sv = (xux 2tut )
+ ct
+u
.
t
x
u
61
First solve the first order autonomous system of differential equations (initial value problem) which belongs to the vector field V , i.e.
dt0
= t0 ,
d
dx0
= ct0 ,
d
du0
= u0 .
d
From this solution of the initial value problem find the transformation between the prime and unprime system. Then using differentiation and the
chain rule show that the prime and unprime system have the same form.
Problem 11.
w
+ w = F (w)
t
t/4
w(t, x vt)
|w(t, x)|2 dx
Rn
is a conserved quantity.
(v) Show that the Hamilton density
Z
1
1
|w(t, x)|2 + G(|w(t, x)|2 )dx
H(w) =
2
2
n
R
is a conserved quantity.
Problem 12. For a barotropic fluid of index the Navier-Stokes equation
read, in one space dimension
c
c
( 1) v
+v
+
c
=0
t
x
2
x
continuity equation
v
v
2
c
2v
+v
+
c
=0
t
x ( 1) x x2
Eulers equation
(1a)
(1b)
= v,
x
v
1
=
(v 2 + c2 ).
t
x 2
(2)
c(x, t) = g(s)
x
t
(3)
(4)
Show that the Navier-Stokes equation yields the following system of ordinary differential equations
df
1 dg
s g+s f
+ g(2f 1) = 0
(5a)
ds
2 ds
2
d2 f
df
+
2
ds
ds
dg
4
+ s(1 2f ) + 2f (1 f ) = 2g s + g .
s
ds
(5b)
(iii) Show that the continuity equation (1a) admits of a first integral, expressing the law of mass conservation, namely
s2 g(2f 1) = C
(6)
(1)
63
ex2 u
.
x2
u
(2)
Construct a similarity ansatz from the symmetry vector field
+ c2
+ c3 x1
+u
(3)
Z = c1
x1
x2
x1
u
Z1 =
,
x1
Z2 =
,
x2
Z3 = x1
+u ,
x1
u
Z4 = ex2
w(s)
,
(3t)2/3
s=
x
(3t)1/3
where s the similarity variable and find the ordinary differential equation.
Problem 15. Let u be the velocity field and p the pressure. Show that
the Navier-Stokes equations
u
+ (u )u u + p = 0,
t
divu = 0
,
t
,
x1
,
x2
x3
+
, t
+
, t
+
x1
u1
x2
u2
x3
u3
x2
+ u1
u2
,
x1
x2
x1
u2
u1
x2
x3
+ u2
u3
,
x3
x2
u3
u2
x3
x1
+ u3
u1
x1
x3
u1
u3
2t + x1
+ x2
+ x3
u1
u2
u3
2p .
t
x1
x2
x3
u1
u2
u3
p
t
S
1 |Q|2
=
Problem 17.
1
= F 1 + G2 ,
1
= A1 + B2
2
= H1 F 2 ,
2
= C1 A2 .
Problem 18.
y
+ 2u1
+ 3u2
.
t
y
u1
u2
65
Problem 19.
xu
+t
.
t 2 x
u
Problem 20. Find the Lie symmetry vector fields for the Monge-Ampere
equations
2 2
2u 2u
u
= K
2
2
t x
xt
where K = +1, K = 0, K = 1.
Problem 21. Let u = (u1 , u2 , u3 )T , v = (v1 , v2 , v3 )T and u21 +u22 +u23 = 1,
v12 + v22 + v32 = 1. Let be the vector product. The O(3) chiral field
equations are the system of partial differential equations
u u
+
+ u Rv = 0,
t
x
v v
+ v Ru = 0
t
x
~2
+ U =
2m
U = 4Gm2 ||2
where m is the mass, G the gravitational constant and the energy eigenvalues. The normalization consition is
Z
|(x1 , x2 , x3 )|2 dx1 dx2 dx3 = 1 .
R3
is invariant under
x0 (x, t)
ct0 (x, t)
=
cosh()
sinh()
sinh()
cosh()
x
ct
67
68
Bibliography
Bibliography
Books
Courant R. and Hilbert D.
Methods of Mathematical Physics, Volume 1, Wiley-VCH Verlag (2004)
Courant R. and Hilbert D.
Methods of Mathematical Physics, Volume 2, Wiley-VCH Verlag (2008)
DiBenedetto E.
Partial Differential Equations, second edition, Springer (2010)
John F.
Partial Differential Equations, Volume I, Springer 1982
Sneddon I.
Elements of Partial Differential Equations, McGraw-Hill (1957)
Steeb W.-H.
The Nonlinear Workbook, fourth edition, World Scientific Publishing, Singapore, 2008
Steeb W.-H.
Problems and Solutions in Theoretical and Mathematical Physics, Second
Edition, Volume I: Introductory Level
World Scientific Publishing, Singapore (2003)
Steeb W.-H.
Problems and Solutions in Theoretical and Mathematical Physics, Second
Edition, Volume II: Adanced Level
World Scientific Publishing, Singapore (2003)
Papers
Ouroushev D., Martinov N. and Grigorov A., An approach for solving
the two-dimensional sine-Gordon equation, J. Phys. A: Math. Gen. 24,
L527-528 (1991)
Index
Antikink, 25
Boltzmann equation, 28, 31
Born-Infeld equation, 57
Born-Infeld type equation, 26
Boussinesq-Oberbech equations, 36
Broadwell model, 30
Burgers equation, 20, 45
Carleman model, 28, 29
Cauchy principal value, 32
Central difference scheme, 7
Chiral field equations, 44
Cole-Hopf transformation, 20
Conservation law, 24, 27
Conserved density, 24
Conserved flux, 24
dAlembert equation, 23
Diffusion equation, 5
Direct method, 33
Dispersion relation, 6, 49
Drift diffusion equation, 1, 4
Electron-beam plasma system, 34
Energy center, 25
Entropy density, 30
Entropy flux, 30
Euler-Lagrange equation, 24
Fisher equation, 13
Fitzhugh-Nagum equation, 32
Fokker-Planck equation, 10
Galilean transformation, 65
Ginzburg-Landau equation, 37
H thoerm, 28
Harry-Dym equation, 57
Heisenberg ferromagnetic equation,
43
Hilbert transform, 32
Hodograph invariant, 25
Hodograph transfromation, 25
Hopf-Cole transformation, 18
Inverse scattering transform, 40
Kadomtsev Petviashvilli equation, 49
Kadomtsev-Petviashvili equation, 41
Kink, 25
KLamb substitution, 22
Korteweg-de Vries equation, 13, 15
Korteweg-de Vries-Burgers equation,
13
Landau-Lifshitz equation, 44
Langmuir oscillations, 19
Laplace equation, 7
Laplace transformation, 11
Lax pair, 12
Lie-Backlund vector fields, 60
Magneto-Hydro-Dynamics equations,
57
Maxwells equation, 6
Monge-Ampere equation, 26
Navier-Stokes equation, 19, 61
Nonlinear Dirac equation, 33
Nonlinear one-dimensional Schrodinger
equation, 37
Nonlinear Schrodinger equation, 17
69
70
Pauli matrices, 33
Pauli spin matrices, 38
periodic boundary conditions, 37
Polar form, 9
Polytropic gas dynamics, 42
Product ansatz, 4
Raman scattering equations, 58
Recursion operator, 45
Reduced variable, 33
Renyi entropies, 29
Residue theorem, 32
Riemann invariants, 42
Sawada-Kotera equation, 41
Shock wave solutions, 13
Sine-Gordon equation, 12, 14, 26, 38
Sine-Hilbert equation, 31
Standing waves, 33
Stationary states, 33
torus, 5
Traveling wave solution, 12
Two-dimensional nonlinear Schrodinger
equation, 27
Velocity potential, 62
Zakharov-Kuznetsov equation, 16
Index