Partial Differential Equations: T. Muthukumar Tmk@iitk - Ac.in April 18, 2014
Partial Differential Equations: T. Muthukumar Tmk@iitk - Ac.in April 18, 2014
T. Muthukumar
[email protected]
Notations vii
1 PDE: An Introduction 1
1.1 Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Well-Posedness of PDE . . . . . . . . . . . . . . . . . . . . . . 9
1.3 Three Basic PDE: History . . . . . . . . . . . . . . . . . . . . 11
1.4 Continuity Equation . . . . . . . . . . . . . . . . . . . . . . . 13
3 Classification by Characteristics 61
3.1 Semilinear . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
iii
CONTENTS iv
4 Wave Equation 81
4.1 One Dimension . . . . . . . . . . . . . . . . . . . . . . . . . . 81
4.1.1 The Vibrating String: Derivation . . . . . . . . . . . . 81
4.1.2 Travelling Waves . . . . . . . . . . . . . . . . . . . . . 83
4.1.3 Domain of Dependence and Influence . . . . . . . . . . 88
4.1.4 Standing Waves: Separation of Variable . . . . . . . . 91
4.2 Higher Dimensions . . . . . . . . . . . . . . . . . . . . . . . . 95
4.2.1 Spherical Means . . . . . . . . . . . . . . . . . . . . . . 96
4.2.2 Odd Dimension . . . . . . . . . . . . . . . . . . . . . . 100
4.2.3 Inhomogeneous Wave equation . . . . . . . . . . . . . . 101
4.3 Eigenvalue Problem of Laplacian . . . . . . . . . . . . . . . . 103
4.3.1 In Rectangle . . . . . . . . . . . . . . . . . . . . . . . . 104
4.3.2 In Disk . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
Appendices 171
Bibliography 189
Index 191
CONTENTS vi
Notations
Symbols
Pn ∂2
∆ i=1 ∂x2i
∂
∇ , ∂ , . . . , ∂x∂n
∂x1 ∂x2
Function Spaces
Cc∞ (X) is the class of all infinitely differentiable functions on X with com-
pact support
General Conventions
vii
NOTATIONS viii
Chapter 1
PDE: An Introduction
(i) Modelling the problem or deriving the mathematical equation (in our
case it would be formulating PDE). The derivation process is usually a
result of conservation laws or balancing forces.
1.1 Definitions
Recall that the ordinary differential equations (ODE) dealt with functions
of one variable, u : Ω ⊂ R → R. The subset Ω could have the interval form
(a, b). The derivative of u at x ∈ Ω is defined as
u(x + h) − u(x)
u′ (x) := lim ,
h→0 h
1
CHAPTER 1. PDE: AN INTRODUCTION 2
provided the limit exists. The derivative gives the slope of the tangent line
at x ∈ Ω. How to generalise this notion of derivative to a function u :
Ω ⊂ Rn → R? These concepts are introduced in a course on multi-variable
calculus. However, we shall jump directly to concepts necessary for us to
begin this course.
Let Ω be an open subset of Rn and let u : Ω → R be a given function.
We denote the directional derivative of u at x ∈ Ω, along a vector ξ ∈ Rn , as
∂u u(x + hξ) − u(x)
(x) = lim ,
∂ξ h→0 h
provided the limit exists. The directional derivative of u at x ∈ Ω, along
the standard basis vectors ei = (0, 0, . . . , 1, 0, . . . , 0) is called the i-th partial
derivative of u at x and is given as
∂u u(x + hei ) − u(x)
ux i = (x) = lim .
∂xi h→0 h
Similarly, one can consider higher order derivatives, as well. We now in-
troduce Schwartz’s multi-index notation for derivative, which will be used
to denote a PDE in a concise form. A multi-index α ∈ Zn+ is a n-tuple
α = (α1 , . . . , αn ) of non-negative integers and let |α| = α1 + . . . + αn . If α
and β are two multi-indices, then α ≤ β means αi ≤ βi for all 1 ≤ i ≤ n.
Also, α ± β = (α1 ± β1 , . . . , αn ± βn ), α! = α1 ! . . . αn ! and xα = xα1 1 . . . xαnn
for any x ∈ Rn . A k-degree polynomial in n variables can be represented as
X
aα x α .
|α|≤k
(i) Note that u(x, t) = c is a solution of the PDE, for any constant c ∈ R.
Thus, we have a family of solutions depending on c.
i.e., the highest order derivative coefficient contains derivative only upto the
previous order. Finally, we say F is fully nonlinear if it depends nonlinearly
on the highest order derivatives.
Note that the classification is heirarchical, i.e., we have the inclusion
Thus, common sense tells that we classify a PDE based on the smallest class
it sits in the heirarchy.
CHAPTER 1. PDE: AN INTRODUCTION 6
(v) ux + sin(uy ) = u.
(xv) ux + u3 = 1.
(xvi) uxxyy + ex ux = y.
(vii) uxx + 4uxy + 4uyy = 0 for v = y − 2x and w = x. (should get uww = 0).
Any PDE not meeting the above criteria is said to be ill-posed. If too many
initial/boundary conditions are specified, then the PDE will have no solu-
tion. If too few initial/boundary conditions are specified, then the PDE will
not have unique solution. However, if the right amount of initial/boundary
conditions are specified, but at wrong places, then the solution may exist and
be unique except that that it may not depend continuously on the initial or
boundary data.
CHAPTER 1. PDE: AN INTRODUCTION 10
has the trivial solution u(x, t) = 0. We consider the same problem with a
small change in data, i.e.,
utt = uxx in R × (0, ∞)
u(x, 0) =0
ut (x, 0) = ε sin xε
has the trivial solution u(x, t) = 0. We consider the same problem with a
small change in data, i.e.,
utt = −uxx in R × (0, ∞)
u(x, 0) =0
ut (x, 0) = ε sin xε
and
with u(x, 0) = 1.
The fundamental question is, given a PDE, find these extra conditions
that make a PDE well-posedness.
∂ 2 u(x, t) ∂ 2 u(x, t)
= .
∂t2 ∂x2
This was then generalised to two and three dimensions by Euler (1759) and
D. Bernoulli (1762), i.e.,
∂ 2 u(x, t)
= ∆u(x, t),
∂t2
P3 ∂2
where ∆ = i=1 ∂x2i .
CHAPTER 1. PDE: AN INTRODUCTION 12
1
a mathematical formulation of conservation laws
CHAPTER 1. PDE: AN INTRODUCTION 13
∂u(x, t)
= ∆u(x, t),
∂t
P ∂2
where ∆ = 3i=1 ∂x 2 . The heat flow model was based on Newton’s law of
i
cooling.
Thus, by the beginning of 19th century, the three most important PDE’s
were identified.
Thus,
Z Z ! Z
1 h→0
ρ(x, t + h) dx − ρ(x, t + h) dx → ρ(x, t)v(x, t) · ν(x) dσ.
h Ωt+h Ωt ∂Ωt
In this chapter, we try to find the general solutions and particular of first
order PDE. A general first PDE has the form
F (∇u(x), u(x), x) = 0.
17
CHAPTER 2. FIRST ORDER PDE 18
Example 2.3. Find the first order PDE, by eliminating the arbitrary constants
a and b, satisfied by u
(i) u(x, y) = (x + a)(y + b)
(ii) u(x, y) = ax + by
Proof. (i) Differentiating the given equation w.r.t x and y, we get
ux = y + b, uy = x + a,
respectively. Eliminating a and b, we get
ux uy = u.
CHAPTER 2. FIRST ORDER PDE 19
ux = a, uy = b
xux + yuy = u.
Exercise 7. Find the first order PDE, by eliminating the arbitrary function
f , satisfied by u(x, y) = ex f (2x − y). (Answer: ux + 2uy − u = 0).
Exercise 8. Find the first order PDE, by eliminating the arbitrary function
f , satisfied by u(x, y) = e−4x f (2x − 3y). (Answer: 3ux + 2uy + 12u = 0).
This is a simple first order linear equation called the transport equation.
Solving
We wish to solve the transport equation ut + bux = 0 which describes the
motion of an object moving with constant speed b, as seen by a fixed observer
A.
Let us imagine that another observer B (say, on a skateboard) moving
with speed b observes the same object, moving in the direction of the object.
CHAPTER 2. FIRST ORDER PDE 20
For B the wave would appear stationary while for A, the fixed observer,
the wave would appear to travel with speed b. What is the equation of the
motion of the object for the moving observer B? To understand this we need
to identify the coordinate system for B relative to A. Let us fix a point x
at time t = 0. After time t, the point x remains as x for the fixed observer
A, while for the moving observer B, the point x is now x − bt. Therefore
the coordinate system for B is (ξ, s) where ξ = x − bt and s = t. Let v(ξ, s)
describe the motion of the object from B’s view. Then the PDE describing
the motion of object, as seen by B, is vs (ξ, s) = 0. Therefore, v(ξ, s) = f (ξ),
for some arbitrary function f (sufficiently differentiable), is the solution from
B’s perspective. To solve the problem from A’s perspective, we observe the
relations
ut = vξ ξt + vs st = −bvξ + vs and
ux = v ξ ξ x + v s s x = v ξ .
Therefore, ut + bux = −bvξ + vs + bvξ = vs and, hence, u(x, t) = v(ξ, s) =
f (ξ) = f (x − bt).
The equation of a line passing through (x, t) and parallel to (b, 1) is (x, t)+
s(b, 1), for all s ∈ (−t, ∞), i.e., (x + sb, t + s). Thus, for a fixed (x, t) ∈
Rn × (0, ∞), we set v(s) := u(x + sb, t + s) for all s ∈ (−t, ∞). Consequently,
Since (x, t) was arbitrary in Rn × (0, ∞), we have u(x, t) = g(x − tb) for
all x ∈ Rn and t ≥ 0. Thus, g(x − tb) is a classical solution to (2.2.1) if
g ∈ C 1 (Rn ). If g 6∈ C 1 (Rn ), we shall call g(x − tb) to be a weak solution of
(2.2.1).
∂u
(x, t) + b · ∇u(x, t) = f (x, t) in Rn × (0, ∞) (2.2.2)
∂t
As before, we set v(s) := u(x + sb, t + s), for all s ∈ R, and for any given
point (x, t) ∈ Rn × (0, ∞). Thus,
dv(s) ∂u
= b · ∇u(x + sb, t + s) + (x + sb, t + s) = f (x + sb, t + s).
ds ∂t
CHAPTER 2. FIRST ORDER PDE 22
Consider,
and set
dxj (s)
= Fpj (p(s), z(s), x(s)).
ds
Thus,
dx(s)
= ∇p F (p(s), z(s), x(s)). (2.4.1)
ds
Now substituting this in the first order equation, we get
dpi (s)
= −Fz (p(s), z(s), x(s))pi (s) − Fxi (p(s), z(s), x(s)).
ds
Thus,
dp(s)
= −Fz (p(s), z(s), x(s))p(s) − ∇x F (p(s), z(s), x(s)). (2.4.2)
ds
Similarly, we differentiate z(s) w.r.t s,
Xn
dz(s) dxj (s)
= uxj (x(s))
ds j=i
ds
n
X
= uxj (x(s))Fpj (p(s), z(s), x(s))
j=i
Thus,
dz(s)
= p(s) · ∇p F (p(s), z(s), x(s)). (2.4.3)
ds
We have 2n + 1 first order ODE called the characteristic equations of (2.5.1).
The steps derived above can be summarised in the following theorem:
Theorem 2.4.1. Let u ∈ C 2 (Ω) solve (2.5.1) and x(s) solve (2.4.1), where
p(s) = ∇u(x(s)) and z(s) = u(x(s)). Then p(s) and z(s) solve (2.4.2) and
(2.4.3), respectively, for all x(s) ∈ Ω.
CHAPTER 2. FIRST ORDER PDE 26
We end this section with remark that for linear, semi-linear and quasi-
linear PDE one can do away with (2.4.2), the ODE corresponding to p,
because for these problems (2.4.3) and (2.4.1) form a determined system.
However, for a fully nonlinear PDE one needs to solve all the 3 ODE’s to
compute the characteristic curve. The method of characteristics may be
generalised to higher order hyperbolic PDE’s.
(a) In first order linear problem, the ODE reduces to, as follows: Let
dz(s)
= b(x(s)) · p(s) = b(x(s)) · ∇u(x(s)) = −c(x(s))z(s).
ds
we have
dx(s)
= ∇p F = b(x(s)).
ds
Also,
dz(s)
= b(x(s), z(s)) · p(s) = b(x(s), z(s)) · ∇u(x(s)) = −c(x(s), z(s)).
ds
CHAPTER 2. FIRST ORDER PDE 27
dz
and = c(x(s), y(s), u(x(s), y(s))).
ds
The three ODE’s obtained are called characteristic equations. The union
of these characteristic (integral) curves give us the integral surface.
We shall now illustrate the method of characteristics for various examples
of first order PDE.
Example 2.5 (Linear Transport Equation). We have already solved the linear
transport equation by elementary method. We solve the same now using
method of characteristics. Consider the linear transport equation in two
variable,
ut + bux = 0, x ∈ R and t ∈ (0, ∞),
where the constant b ∈ R is given. Thus, the given vector field V (x, t) =
(b, 1, 0). The characteristic equations are
dx dt dz
= b, = 1, and = 0.
ds ds ds
Solving the 3 ODE’s, we get
x(s) = bs + c1 , t(s) = s + c2 , and z(s) = c3 .
Eliminating the parameter s, we get the curves (lines) x − bt = a constant
and z = a constant. Therefore, u(x, t) = f (x − bt) is the general solution,
for an arbitrary function f .
In the inhomogeneous equation case, where the ODE corresponding to z
is not zero, we intend to seek a function Φ(x, u) such that the solution u is
defined implicitly by Φ(x, u) = a constant. Suppose there is such a function
Φ then, by setting z := u,
Φxi + Φz uxi = 0 ∀i = 1, 2, . . . , n.
Assuming Φz 6= 0, we get
−Φxi
ux i = ∀i = 1, 2, . . . , n.
Φz
If u solves b(x, u) · ∇u(x) = c(x, u) then (b(w), c(w)) · ∇w Φ = 0 is a homoge-
neous first order equation of Φ in n + 1 variables. In this case,
dw(s) dΦ(w(s))
= (b, c) and = 0.
ds ds
CHAPTER 2. FIRST ORDER PDE 29
for some arbitrary smooth functions f and g. Compare this answer with
Example 2.3(ii).
Example 2.8. Let us compute the general solution (in terms of arbitrary
functions) of the first order PDE yu(x, y)ux (x, y) + xu(x, y)uy (x, y) = xy.
The characteristic equations are
dx dy dz
= yz, = xz and = xy.
ds ds ds
Hence,
dx dy
0 = x −y
ds ds
2
d(x ) d(y 2 )
= −
ds ds
d(x2 − y 2 )
= .
ds
Thus, x2 − y 2 = c1 and, similarly, x2 − z 2 = c2 . Hence, the general solution
is F (x2 − y 2 , x2 − z 2 ) = 0 for some arbitrary function F . Explicitly, for some
f and g,
Example 2.9. Let us compute the general solution (in terms of arbitrary
functions) of the first order PDE 2yux + uuy = 2yu2 . The characteristic
equations (ODE’s) are
dx dy dz
= 2y(s) = z(s) and = 2y(s)z 2 (s).
ds ds ds
Solving in the parametric form is quite cumbersome, because we will have a
second order nonlinear ODE of y, y ′′ (s) = 2y(y ′ )2 . However, for u 6= 0, we get
dz 2
dy
= 2y(s)z(s) solving which we get ln |z| = y 2 + c1 or z = c2 ey . Similarly,
dy y2 2
dx
= z
2y
= c22y
e
solving which we get c2 x + e−y = c3 or, equivalently,
2 2
xze−y + e−y = c3 . Thus, the general solution is
2 2
F (e−y (1 + xz), e−y z) = 0.
or 2 2
u(x, y) = f e−y (1 + xu) ey
CHAPTER 2. FIRST ORDER PDE 31
−1
u(x, y) =
x + f (y − x2 )
dx dy dz
= y(s) = −x(s) and = 2xyz.
ds ds ds
To avoid cumbersome ODE, let us begin by assuming y 6= 0, then dividing
the entire equation by y, we get
dx dy dz
=1 = −x(s)/y(s) and = 2xz.
ds ds ds
Solving which we get x(s) = s + c1 , y(s) = −s2 − 2c1 s + 2c2 and |z(s)| =
2
c3 es +2c1 s . Eliminating s between x and y, we get the characteristic curves
2
to be y 2 + x2 = a constant and z = c4 ex . Thus, the general solution is
2
F (y 2 + x2 , e−x z) = 0. Explicitly,
2
u(x, y) = f (y 2 + x2 )ex
Example 2.12. Let us compute the general solution (in terms of arbitrary
functions) of the first order PDE ux1 + ex1 ux2 + ex3 ux3 = (2x1 + ex1 )eu . The
characteristic equations (ODE’s) are
dx1 dx2 dx3 dz
=1 = e x1 = ex3 and = (2x1 + ex1 )ez(s) .
ds ds ds ds
Solving which we get x1 (s) = s + c1 , x2 (s) = ec1 es + c2 , e−x3 (s) = −s + c3 and
The idea is that for each x ∈ Ω, the graph of v is tangent to the graph of
u(·; a) for a = φ(x).
Example 2.18. The complete integral of the nonlinear PDE u2 (1+ |∇u|2 ) = 1
is u(x; a) = ±(1 − |x − a|2 )1/2 with |x − a| < 1. Now, solving Da u =
x−a
± (1−|x−a|2 )1/2 = 0 for a, we get a = φ(x) := x. Thus, the envelope is v(x) =
±1.
Example 2.22. Let us compute a complete integral of the first order PDE
ux uy = u(x, y). The equation is of the form F (p, z, x) = p1 p2 − z. The
characteristic equations are (using (2.4.1))
dx(s) dy(s)
, = (p2 (s), p1 (s)),
ds ds
(using (2.4.2))
dp(s)
= p(s)
ds
and (using (2.4.3))
dz(s)
= (p1 (s), p2 (s)) · (p2 (s), p1 (s)) = 2p1 (s)p2 (s)
ds
u(x, y) = (y − a)(x − b) + c3
dp(s)
= p(s)
ds
Thus, on integrating, we get p1 (s) = c1 es and p2 (s) = c2 es . Therefore,
√
p1 /p2 = a. Using this equation with p1 p2 = z, we get p1 = ± az and
CHAPTER 2. FIRST ORDER PDE 36
p
p2 = ± z/a. Now,
dz(s) dx(s) dy(s)
= (p1 (s), p2 (s)) · ,
ds ds ds
√ dx(s) p dy(s)
= ± az ± z/a
ds ds
1 dz(s) √ dx(s) √ dy(s)
√ = ± a + 1/ a
z ds ds ds
√ √ √
2 z = ± ax + y/ a + c3 .
Thus,
2
1 √ √
u(x, y) = b + ax + y/ a
2
Note that previous two examples compute two different complete integral
for same equation. However, in both examples, no choice of a and b will give
the zero solution u ≡ 0. Thus, u ≡ 0 is called singular solution.
Example 2.24. Let us find the complete integral, general solution and singular
solution of the fully non-linear PDE u2x + u2y = 1 + 2u. Since F (p, z, x) =
p21 + p22 − 1 − 2z, the ODEs are
dx(s) dy(s)
, = (2p1 (s), 2p2 (s)),
ds ds
dp(s)
= (2p1 , 2p2 )
ds
Thus, on dividing and integrating, we get p1 /p2 = a. Using the PDE, we get
(1 + a2 )p22 = 1 + 2z. Thus,
r r
1 + 2z 1 + 2z
p2 = ± p1 = ±a .
1 + a2 1 + a2
CHAPTER 2. FIRST ORDER PDE 37
Now,
dz(s) dx(s) dy(s)
= (p1 (s), p2 (s)) · ,
ds ds ds
r
1 + 2z dx(s) dy(s)
= ± a +
1 + a2 ds ds
1 dz(s) 1 dx(s) dy(s)
√ = ±√ a +
1 + 2z ds 1 + a2 ds ds
√ ax + y
1 + 2z = ± √ ± b.
1 + a2
Thus, 2
1 ax + y 1
u(x, y) = √ +b −
2 1 + a2 2
is a complete integral. Note that no choice of a and b will give the constant
solution u = −1/2. Thus, u ≡ −1/2 is called singular solution.
Exercise 9. Find the general solution (in terms of arbitrary functions) of the
given first order PDE
(i) xux + yuy = xe−u with x > 0. (Answer: u(x, y) = f (y/x) for some
arbitrary f ).
(iv) x(y 2 − u2 )ux − y(u2 + x2 )uy = (x2 + y 2 )u. (Answer: u(x, y) = xy f (x2 +
y 2 + u2 )).
2.6.1 Quasilinear
Consider the general first order qausilinear PDE with n independent variable
subset of the hyperplane whose final coordinate is zero. Let y = φ(x) and
v(y) := u(x). Then
Xn
∂φk
ux i = v yk .
k=1
∂xi
Using the above equation, the first order linear PDE becomes
n
X ∂φk
bi vy = c̃(v, y), (2.6.2)
i,k=1
∂xi k
(a(γ1 (r), γ2 (r)), b(γ1 (r), γ2 (r))) · (−γ2′ (r), γ1′ (r)) = (2, 3) · (−m, 1) = 3 − 2m.
Thus, the line is not a non-characteristic for m = 3/2, i.e., all lines with slope
3/2 is not a non-characteristic.
dx dy du
(r, s) = a(x, y, u); (r, s) = b(x, y, u) and (r, s) = c(x, y, u)
ds ds ds
such that on Γ, x(r, s) = γ1 (r), y(r, s) = γ2 (r) and z(r, s) = φ(r). The
non-characteristic property of Γ implies that the following Jacobian on Γ:
∂x ∂x
∂r
∂y
∂s
∂y = (a(γ1 , γ2 , φ), b(γ1 , γ2 , φ)) · (−γ2′ , γ1′ ) 6= 0.
∂r ∂s
By implicit function theorem (cf. Theorem C.0.2), one can solve for r and
s in terms of x and y, locally, in the neighbourhood of Γ. Set v(x, y) =
u(r(x, y), s(x, y)). We will show that v is a unique solution to (2.6.3). Note
that, by implicit function theorem, the value of v on Γ is
Moreover,
(ii) For the PDE given above, check for the characteristic property of the
following curves
(a) u(x, x) = x4 on y = x
(b) u(x, (3x − 1)/2) = x2 on y = (3x − 1)/2
(c) u(x, (3x − 1)/2) = e−4x .
Hence,
(ii) (a) Parametrise the curve y = x as Γ(r) : r 7→ (r, r). Thus γ1 (r) =
γ2 (r) = r. Since the coefficients of the PDE are a(r) = 2 and
b(r) = 3, we have
Hence Γ is non-characteristic.
(b) Parametrise the curve y = (3x − 1)/2 as Γ(r) : r 7→ (r, (3r − 1)/2).
Hence γ1 (r) = r and γ2 (r) = (3r − 1)/2 and
and
u(x, y) = (3x − 2y)4 e8(x−y) .
Thus, we have a unique solution u.
(b) Using the given condition, we have F (1, x2 e4x ) = 0. Either f (1) =
x2 e4x or f (x2 e4x ) = 1. The first case is not valid (multi-valued
function). The second case corresponds to z = e−4x which will not
satisfy the Cauchy data. Hence there is no solution u solving the
given PDE with the given data.
CHAPTER 2. FIRST ORDER PDE 43
(c) Once again using the given condition, we have F (1, x2 e4x ) = 0.
Either f (1) = x2 e4x or f (x2 e4x ) = 1. The first case is not valid
(multi-valued function). The second case corresponds to z = e−4x
which will satisfy the Cauchy data. Since there many choices of f
that satisfies f (x2 e4x ) = 1, we have infinite number of solutions (or
choices for) u that solves the PDE.
Example 2.27. (i) Find the general solution (in terms of arbitrary func-
tions) of the first order PDE ux (x, y) + uy (x, y) = 1.
(ii) For the PDE given above, check for the characteristic property of the
following curves
dx dy dz
= 1, = 1 and = 1.
ds ds ds
Hence,
x(s) = s + c1 y(s) = s + c2 and z(s) = s + c3 .
Thus, y − x = c4 and z − x = c5 or x − y = c4 and z − y = c5 . Hence,
for some f and g,
(ii) (a) Parametrise the curve x-axis as Γ(r) : r 7→ (r, 0). Thus γ1 (r) = r
and γ2 (r) = 0. Since the coefficients of the PDE are a(r) = 1 and
b(r) = 1, we have
Hence Γ is non-characteristic.
(b) Parametrise the curve y = x as Γ(r) : r 7→ (r, r). Hence γ1 (r) =
r = γ2 (r) and
(a) Now, u(x, 0) = φ(x) implies f (x) = φ(−x) + x or g(x) = φ(x), and
We know that the general solution of the transport equation is u(x, t) = u(x−
bt, 0). Thus, u(x, t) = φ(x − bt) is the unique solution of (2.6.4). We derive
the particular solution of the Cauchy problem (2.6.4) using parametrisation
CHAPTER 2. FIRST ORDER PDE 45
of the data curve Γ. The example also shows how the information on data
curve Γ is reduced as initial condition for the characteristic ODE’s. Note that
in the example below the data curve Γ is parametrised using the variable r
and the characteristic curves is parametrised using the variable s.
Example 2.28. We shall compute the solution of the Cauchy problem (2.6.4).
We first check for non-characteristic property of Γ. Note that Γ ≡ {(x, 0)},
the x-axis of xt-plane, is the (boundary) curve on which the value of u is given.
Thus, (Γ, φ) = {(x, 0, φ(x))} is the known curve on the solution surface of u.
We parametrize the curve Γ with r-variable, i.e., Γ = {γ1 (r), γ2 (r)} =
{(r, 0)}. Γ is non-characteristic, because (b, 1) · (0, 1) = 1 6= 0. The charac-
teristic equations are:
x(r, 0) = c1 (r) = r
We solve for r, s in terms of x, t and set u(x, t) = z(r(x, t), s(x, t)).
Example 2.29. Let Ω := {(x, y) ∈ R2 | x > 0, y > 0}. Let Γ := {(x, 0) | x >
0}. Consider the linear PDE
xuy (x, y) − yux (x, y) = u(x, y) in Ω
u(x, 0) = φ(x) on Γ.
The parametrisation of the initial curve is Γ(r) := (r, 0) for r > 0. Therefore,
(a(γ1 (r), γ2 (r)), b(γ1 (r), γ2 (r))) · (−γ2′ (r), γ1′ (r)) = (0, r) · (0, 1) = r 6= 0.
Note that
d2 x(r, s) d2 y(r, s)
= −x(r, s) and = −y(r, s).
ds ds
Then, x(r, s) = c1 (r) cos s + c2 (r) sin s and y(r, s) = c3 (r) cos s + c4 (r) sin s.
Using the initial condition, we get c1 (r) = r and c3 (r) = 0. Also,
dx(r, s)
0 = −y(r, 0) = |s=0 = −c1 (r) sin 0 + c2 (r) cos 0 = c2 (r).
ds
and, similarly, c4 (r) = r. Also, z(r, s) = c5 (r)es where c5 (r) = φ(r). Thus,
we have (x(r, s), y(r, s)) = (r cos s, r sin s), where r > 0 and 0 ≤ s ≤ π/2.
Hence, r = (x2 + y 2 )1/2 and s = arctan(y/x). Therefore, for any given (x, y),
we have p
u(x, y) = z(r, s) = φ(r)es = φ( x2 + y 2 )earctan(y/x) .
Example 2.30. Let Ω := {(x, y) ∈ R2 | y > 0}. Let Γ := {(x, 0) | x ∈ R}.
Consider the semi-linear PDE
ux (x, y) + uy (x, y) = u2 (x, y) in Ω
u(x, 0) = φ(x) on Γ.
The parametrisation of the initial curve is Γ(r) := (r, 0) for all r ∈ R. There-
fore,
(a(γ1 (r), γ2 (r)), b(γ1 (r), γ2 (r))) · (−γ2′ (r), γ1′ (r)) = (1, 1) · (0, 1) = 1 6= 0.
CHAPTER 2. FIRST ORDER PDE 47
2.6.2 Nonlinear
Definition 2.6.5. We say Γ ⊂ Ω ⊂ Rn is non-characteristic for the nonlin-
ear Cauchy problem
F (∇u(x), u(x), x) = 0 x ∈ Ω
(2.6.5)
u = g on Γ
if Γ is nowhere tangent to the Monge cone, i.e., there exists function v such
that F (v(r), g(r), γ1 (r), . . . , γn (r)) = 0 and vi = gxi , for i = 1, 2, . . . , n − 1
on Γ, and satisfies
n
X
Fpi (v, g, γ1 , . . . , γn ) · ν(x) 6= 0
i=1
Example 2.32. For a nonlinear PDE the concept P of characteristics also depend
on initial values. Consider the nonlinear PDE ni=1 u2xi = 1 with u(x′ , 0) =
φ(x′ ). Note that u = xi is a solution with φ(x′ ) = xi . However, any choice
of φ such that |∇x′ φ(x′ )|2 > 1 has no solution.
Suppose v is a solution of F (x, u, ∇u) = 0 in a neighbourhood of some
point x0 ∈ Γ.
Example 2.33. (i) Find the general solution (in terms of arbitrary func-
tions) of the first order PDE xux (x, y) + 2xuuy (x, y) = u(x, y).
(ii) For the PDE given above, check if the following curves in xy-plane are
non-characteristic and discuss the particular solutions of the PDE
(a, b) · (−γ2′ (r), γ1′ (r)) = (r, 4r2 ) · (−4r, 1) = −4r2 + 4r2 = 0.
(a, b)·(−γ2′ (r), γ1′ (r)) = (r, 4r3 )·(−9r2 , 1) = −9r3 +4r3 = −5r3 6= 0
(a, b) · (−γ2′ (r), γ1′ (r)) = (r, 2r3 ) · (−3r2 , 1) = −3r3 + 2r3 = −r3 6= 0
Example 2.34. We now give an example of a fully non linear PDE. Let Ω :=
{(x, y) ∈ R2 | x > 0}. Let Γ := {(0, y) | y ∈ R}. Consider the fully
non-linear PDE
ux uy = u(x, y) in Ω
u(0, y) = y 2 on Γ.
The parametrisation of the initial curve is Γ(r) := (0, r) for all r ∈ R. We
first look for the functions α and β such that
and
p2 (r, 0) = β(r) = 2r.
Thus, on integrating, we get p1 (r, s) = (r/2)es and p2 (r, s) = 2res , for all
s ∈ R. Using p, we solve for x to get x(r, s) = 2r(es − 1) and y(r, s) =
(r/2)(es + 1). Solving for z, we get
r2 2s
z(r, s) = (e + 1).
2
Solving r and s in terms of x and y, we get
4y − x x + 4y
r= and es = .
4 4y − x
(x+4y)2
Hence u(x, y) = z(r(x, y), s(x, y)) = 16
.
Example 2.35. Consider the fully non-linear PDE
ux uy = u(x, y) in R2
u(x, 1 + x) = x2 on Γ.
and
p2 (r, 0) = β(r) = r.
Thus, on integrating, we get p1 (r, s) = res and p2 (r, s) = res , for all s ∈ R.
Using p, we solve for x to get x(r, s) = res and y(r, s) = res + 1. Solving for
z, we get z(r, s) = r2 e2s . Note that there is no unique way of solving r and s
in terms of x and y. In fact, we have three possible represenation of u, viz.,
u = x2 , u = (y − 1)2 and u = x(y − 1). This is because the Jacobian is zero.
Of these three possibilities, only u = x(y − 1) satisfies the equation.
Example 2.36. For any given λ ∈ R, consider the fully non-linear PDE
2
ux + u2y = 1 in R2
u(x, x) = λx on Γ.
We can view α(r) = cos θ and β(r) = sin θ where θ is such that cos θ +sin θ =
λ. Since F (p, z, x) = p21 + p22 − 1, we have
Fp1 γ1′ − Fp2 γ2′ (r) = 2(p1 − p2 ) = 2(cos θ − sin θ) = 0, ∀θ = π/4 + kπ,
Thus, r
ay a2 y 2
p2 = − p1 = ± a + 2 .
z z
Now,
dz(s) dx(s) dy(s)
= (p1 (s), p2 (s)) · ,
ds ds ds
√ dx(s) 1 dz(s) dy(s)
∓ a = −p z + ay
ds z 2 + ay 2 ds ds
√ p
∓ ax ∓ b = − z 2 + ay 2
√
u2 (x, y) = (b + x a)2 − ay 2 .
√
Using the initial data, we get 9y 2 = (b + 2y a)2 − ay 2 which is satisfied with
b = 0 and a = 3. Thus, u(x, y)2 = 3(x2 − y 2 ).
Example 2.39. Consider the equation y(u2x − u2y ) + uuy = 0 with u(r2 , 0) =
2r. Note that F (p, z, x, y)√ = y(p21 − p22 ) + zp2 = 0. As in above example,
we get u2 (x, √
y) = (b + x a)2 − ay 2 . Using the initial conditions, we get
4r2 = (b + r2 a)2 . Differentiating this w.r.t to r and solving for r, we get
2 1 2
r =√ √ −b .
a a
CHAPTER 2. FIRST ORDER PDE 55
√
Substituting this in the equation of r, we get b = 1/ a. This gives
2
2 1 √ 1
u = √ +x a − ay 2 = (1 + xa)2 − ay 2 .
a a
We choose the positive root above to keep a > 0 so that all roots above made
sense. Therefore,
!2
p x y2
u2 (x, y) = x2 − y 2 1 + p −p .
x2 − y 2 x2 − y 2
Exercise 10. Find the general solution of the following PDE. Check if the
given data curve is non-characteristic or not. Also find the solution(s) (if it
exists) given the value of u on the prescribed curves.
2
(ix) e2y ux + xuy = xu2 . Also, when u(x, 0) = ex on y = 0.
(xxvi) uux + uy = 1 with u(2r2 , 2r) = 0 for r > 0. (Answer: No solution for
y 2 > 4x).
Exercise 11. Solve the equation xux + 2yuy = 0 with u(1, y) = ey . Does a
solutions exist with data on u(0, y) = g(y) or u(x, 0) = h(x)? What happens
to characteristic curves at (0, 0)?
2
Exercise 12. Solve the equation yux + xuy = 0 with u(0, y) = e−y . In which
region of the plane is the solution uniquely determined?
Exercise 13. Solve the equation ux + yuy = 0 with u(x, 0) = 1. Also, solve
the equation with u(x, 0) = x. If there is no solution, give reasons for non-
existence.
Exercise 14. (i) Find a complete integral of x(u2x + u2y ) = uux with u given
on the curve
(a) u(2y, y) = 5y
(b) u(0, r2 ) = 2r
(ii) Find a complete integral of 4uux − u3y = 0 with u given on the curve
u(0, r) = 4r.
Therefore,
x(r, s) = φ(r)s + r, and t(r, s) = s.
Solving r and s, in terms of x, t and z, we get s = t and r = x−zt. Therefore,
u(x, t) = φ(x − tu) is the solution in the implicit form.
Example 2.40. If the data φ(x) = c, some constant, then u(x, t) = c and the
characteristic curves are t = x/c.
Example 2.41. In the Burgers’ equation, suppose we choose the data φ(x) =
x, then
x
u(x, t) = .
1+t
Note that u has a singularity at (0, −1). This can be observed in the following
way: u takes the constant value c along the line t = x/c − 1 and all these
curves intersect at (0, −1) which means u is multiply defined at (0, −1) or,
rather undefined at (0, −1).
CHAPTER 2. FIRST ORDER PDE 59
Therefore,
1 x≤t
1−x
u(x, t) = t≤x≤1
1−t
0 x ≥ 1.
Note that for t ≤ 1 the solution behaves well, but for t ≥ 1, the characteristics
start crossing each other on the line t = x and u(x, t) takes both 0 and 1 on
the line t = x, for t ≥ 1. This situation is called the shock.
Example 2.43. In the Burgers’ equation, suppose we choose φ to be the
function
−1 x < −1
φ(x) = x −1 ≤ x ≤ 1
1 1 < x.
Then the characteristic curves are
−t + c c < −1
x = c(t + 1) −1 ≤ c ≤ 1
t+c 1 < c.
Therefore,
−1 x + t < −1
x
u(x, t) = −(t + 1) ≤ x ≤ (t + 1)
t+1
1 1 < x − t.
CHAPTER 2. FIRST ORDER PDE 60
Classification by
Characteristics
A general second order PDE is of the form F (D2 u(x), Du(x), u(x), x) = 0,
for each x ∈ Ω ⊂ Rn and u : Ω → R is the unknown. A Cauchy problem
poses the following: Given the knowledge of u on a smooth hypersurface
Γ ⊂ Ω can one find the solution u of the PDE? The knowledge of u on Γ is
said to be the Cauchy data.
What should be the minimum required Cauchy data for the Cauchy prob-
lem to be solved? Viewing the Cauchy problem as an initial value problem
corresponding to ODE, we know that a unique solution exists to the second
order ODE
′′
y (x) + P (x)y ′ (x) + Q(x)y(x) = 0 x ∈ I
y(x0 ) = y0
y ′ (x0 ) = y0′ .
61
CHAPTER 3. CLASSIFICATION BY CHARACTERISTICS 62
3.1 Semilinear
Consider the general second order qausilinear PDE with n independent vari-
able
F (x, u, Du, D2 u) := A(x) · D2 u − D(∇u, u, x), (3.1.1)
where A = Aij is an n × n matrix with entries Aij (x, u, ∇u), D2 u is the
2
Hessian matrix. The dot product in LHS is in Rn . Since we demand the
solution to be in C 2 , the mixed derivatives are equal and we can assume,
without loss generality that, A is symmetric. In fact if A is not symmetric,
we can replace A with As := 12 (A0 + At0 ), which is symmetric since A · D2 u =
As · D2 u.
Let Γ ⊂ Ω be an hypersurface defined implicitly as S(x) = 0 and ∇S 6= 0.
Given u and ∇u · ν on Γ, can we calculate all other derivatives of u on Γ?
To answer this, we map Γ to a hyperplane Γ0 by the map φ : Γ → Γ0 with
its components φ = (φ1 , φ2 , . . . , φn ) such that φn (x) = S(x) and Jφ (x) 6= 0,
non-zero Jacobian of φ, for all x ∈ Ω. For sufficiently regular hypersurfaces
such a mapping exists. Note that Γ0 is a subset of the hyperplane whose
final coordinate is zero. Let y = φ(x) and v(y) := u(x). Then
n
X n
X n
∂φk ∂φk ∂φl X ∂ 2 φk
ux i = v yk ux i x j = v yk yl + v yk .
k=1
∂xi k,l=1
∂xi ∂xj k=1 ∂xi ∂xj
Using the second equation the second order linear PDE becomes
n
X X n
∂φk ∂φl ∂ 2 φk
Aij vy y = D̃(∇v, v, y) − Aij vyk , (3.1.2)
i,j,k,l=1
∂xi ∂xj k l i,k=1
∂xi ∂xj
1
vyk yl (y) = lim [vyk (y1 , . . . , yl + h, . . . , yn−1 , 0) − vyk (y1 , . . . , yl , . . . , yn−1 , 0)].
h→0 h
Note that the elliptic case corresponds to the situation that every hy-
persurface S(x) = 0 with ∇S 6= 0 is non-characteristic corresponding the
elliptic operator, i.e., there are no real characteristics curves. Thus, one can
equivalently define a linear second order PDE to be elliptic at x if
n
X
Aij (x)ξi ξj 6= 0 ∀ξ ∈ Rn \ {0}.
i,j=1
Note that the geometry hidden in the above definition is very similar
to that we encountered in first order equation. Since ν = (−γ̇2 , γ̇1 ) is the
normal to Γ at each point, the above definition says that the curve is non-
characteristic if
2
X
Aij νi νj = Aγ̇22 − 2B γ̇1 γ̇2 + C γ̇12 6= 0
i,j=1
Therefore, the characteristic curves of (3.1.3) are given by the graphs whose
equation is √
dy B ± B 2 − AC
= .
dx A
Thus, we have three situations arising depending on the sign of the dis-
criminant, B 2 − AC. This classifies the given second order PDE based on
the sign of its discriminant d = B 2 − AC.
3.2 Quasilinear
All the arguments of semilinear PDE can be carried over to a quasilinear
PDE A(x, u(x), Du(x)). For each specific point x0 ∈ Ω, u(x0 ) = u0 and A0 =
A(x0 , u0 , ∇u(x0 )). The solutions to characteristic equation for a quasilinear
equation depends on the solution considered. Set
U := {(x, z, p) | x ∈ Ω, z ∈ R, p ∈ Rn }.
is said to be elliptic if the matrix Aij (x, z, p) is positive definite for each
(x, z, p) ∈ U . Further,
n
X
2
0 < α(x, z, p)|ξ| ≤ Aij (x, z, p)ξi ξj ≤ β(x, z, p)|ξ|2 ∀ξ ∈ Rn .
i,j=1
The bounds α(x, z, p) and β(x, z, p) are minimum and maximum eigen-
values, respectively. If β/α is uniformly bounded in U then PDE is uniformly
elliptic. The interesting thing about uniformly elliptic equation is that they
behave very similar to linear elliptic equations.
Example 3.1. Consider the minimal surface equation
!
∇u
∇· p = f (x)
1 + |∇u|2
where the second order coefficients are
2 −1/2 pi pj
Aij (x, z, p) = (1 + |p| ) δij −
1 + |p|2
and
1 1
α(x, z, p) = 2 3/2
β(x, z, p) = .
(1 + |p| ) (1 + |p|2 )1/2
Thus, the equation is not uniformly elliptic.
The minimal surface equation and the capillary equation are not uni-
formly elliptic.
3.3 Examples
Example 3.2 (Wave Equation). For a given c ∈ R, uyy − c2 uxx = 0 is hyper-
bolic. Since A = −c2 , B = 0 and C = 1, we have d = B 2 −AC = c2 > 0. The
eigen values of the coefficient matrix are 1, −c2 which have opposite sign.
Example 3.3 (Heat Equation). For a given c ∈ R, uy − cuxx = 0 is parabolic.
Since A = −c, B = 0 and C = 0, thus d = B 2 − AC = 0. The eigen values
of the coefficient matrix are 0, −c has a zero eigenvalue.
Example 3.4 (Laplace equation). uxx +uyy = 0 is elliptic. Since A = 1, B = 0
and C = 1, thus d = B 2 − AC = −1 < 0. The eigen values of the coefficient
matrix are 1, 1 which have same sign.
CHAPTER 3. CLASSIFICATION BY CHARACTERISTICS 68
Solution. (i) We have already seen the equation is hyperbolic and hence
it should have two characteristic curves. The characteristic curves are
given by the equation
√ √
dy B ± B 2 − AC ± c2 ∓1
= = 2
= .
dx A −c c
Thus, cy ± x = a constant is the equation of the two characteristic
curves. Note that the characteristic curves y = ∓x/c + y0 are boundary
of two cones in R2 with vertex at (0, y0 ).
(ii) We have already seen the equation is parabolic and hence it should
have one characteristic curve. The characteristic curve are given by the
equation √
dy B ± B 2 − AC
= = 0.
dx A
Thus, y = a constant is the equation of the characteristic curve. i.e.,
any horizontal line in R2 is a charateristic curve.
(iii) We have already seen the equation is elliptic and hence has no real
characteristics.
(iv) The equation is of mixed type. In the region x > 0, the characteristic
curves are y ∓ 2x3/2 /3 = a constant.
Exercise 15. Classify the following second order PDE in terms of the number
of characteristics:
Exercise 16. Classify the following second order PDE, in terms of the number
of characteristics, and find their characteristics, when it exists:
A similar procedure also makes a second order PDE into a system of first
order PDEs. Thus, we expect that our classification of second order PDE to
induce a classification for a system of first order PDE.
A general system of m first order linear PDE in n variables will be of the
form
X n
Aj (x, u)uxj = f (x, u) in Ω, (3.4.1)
j=1
on Γ.
Definition 3.4.2. The system of first order PDE (3.4.1) is called hyperbolic
if the m × m matrix
Xn
ξj Aj (x, u)
j=1
Definition 3.4.3. The system of first order PDE (3.4.1) is called elliptic if
n
!
X
det Aj (x, u)ξj = 0
j=1
only for ξ = 0.
Example 3.10 (Beltrami Equations). Consider the system of first order equa-
tions
W (x, y)ux (x, y) − b(x, y)vx (x, y) − c(x, y)vy (x, y) = 0
W (x, y)uy (x, y) + a(x, y)vx (x, y) + b(x, y)vy (x, y) = 0
6 0 and the matrix
where W, a, b, c are given such that W =
a b
b c
is positive definite. Set
W −b 0 −c
A1 = A2 = .
0 a W b
Then the system can be rewritten as
ux uy 0
A1 + A2 = .
vx vy 0
Then
n
!
X W ξ1 −bξ1 − cξ2
det Aj (x, u)ξj = = W (aξ12 + 2bξ1 ξ2 + cξ22 ) 6= 0
W ξ2 aξ1 + bξ2
j=1
only for ξ = 0.
Definition 3.6.1. For any PDE of the form (3.6.1) we define its discrimi-
nant as B 2 − AC.
ux = uw w x + uz z x ,
uy = uw w y + uz z y ,
uxx = uww wx2 + 2uwz wx zx + uzz zx2 + uw wxx + uz zxx
uyy = uww wy2 + 2uwz wy zy + uzz zy2 + uw wyy + uz zyy
uxy = uww wx wy + uwz (wx zy + wy zx ) + uzz zx zy + uw wxy + uz zxy
b2 − ac = (B 2 − AC)J 2 .
Therefore, we need to find w and z such that along the slopes of the charac-
teristic curves, √
dy B ± B 2 − AC −wx
= = .
dx A wy
This means that, using the parametrisation of the characteristic curves,
˙ = 0. Similarly for z. Thus, w and z
wx γ̇1 (s) + wy γ̇2 (s) = 0 and w(s)
are chosen such that they are constant on the characteristic curves.
The characteristic curves are found by solving
√
dy B ± B 2 − AC
=
dx A
and the coordinates are then chosen such that along the characteristic curve
w(x, y) √
2
and z(x, y) = a constant. Note that wx zy − wy zx =
= a constant
2
wy zy A B − AC 6= 0.
Example 3.11. Let us reduce the PDE uxx − c2 uyy = 0 to its canonical form.
Note that A = 1, B = 0, C = −c2 and B 2 − AC = c2 and the equation is
hyperbolic. The characteristic curves are given by the equation
√
dy B ± B 2 − AC
= = ±c.
dx A
Solving we get y ∓ cx = a constant. Thus, w = y + cx and z = y − cx. Now
writing
0 = 4c2 uwz
= uwz .
Example 3.12. Let us reduce the PDE uxx − x2 yuyy = 0 given in the region
{(x, y) | x ∈ R, x 6= 0, y > 0} to its canonical form. Note that A = 1, B = 0,
CHAPTER 3. CLASSIFICATION BY CHARACTERISTICS 77
Example 3.18. Let us reduce the PDE uxx + uxy + uyy = 0 to its canonical
form. Note that A = 1, B = 1/2, C = 1 and B 2 − AC = −3/4 < 0. The
PDE is elliptic. Solving the characteristic equation
√
dy 1 3
= ±i
dx 2 2
√ √
we get 2y = x ± i 3x + c. Let w = 2y − x and z = 3x.
Exercise 17. Rewrite the PDE in their canonical forms and solve them.
√
(a) uxx + 2 3uxy + uyy = 0
Wave Equation
81
CHAPTER 4. WAVE EQUATION 82
the density of the string and T > 0 denote the coefficient of tension of the
string. Let u(x, t) denote the vertical displacement of the string at the point
x and time t.
We shall imagine the string of length L as system of N objects, for N
sufficiently large. Think of N objects sitting on the string L at equidistant
(uniformly distributed). The position of the n-th object on the string is
given by xn = nL/N . One can think of the vibrating string as the harmonic
oscillator of N objects governed by the tension on the string (which behaves
like the spring). Let yn (t) = u(xn , t) denote the displacement of the object xn
at time t. The distance between any two successive objects is h = xn+1 −xn =
L/N . Then mass of each of the N object is mass of the string divided by N .
Since mass of the string is ρ×L, mass of each of the object xn , n = 1, 2, . . . , N ,
is ρh. Thus, by Newton’s second law, ρhyn′′ (t) is same as the force acting on
the n-th object. The force on xn is coming both from left (xn−1 ) and right
(xn+1 ) side. The force from left and right is given as T (yn−1 − yn )/h and
T (yn+1 − yn )/h, respectively. Therefore,
T
ρhyn′′ (t) = {yn+1 (t) + yn−1 (t) − 2yn (t)}
h
T
= {u(xn + h, t) + u(xn − h, t) − 2u(xn , t)}
h
′′ T u(xn + h, t) + u(xn − h, t) − 2u(xn , t)
yn (t) =
ρ h2
Note that assuming u is twice differentiable w.r.t the x variable, the term
on RHS is same as
T 1 u(xn + h, t) − u(xn , t) u(xn − h, t) − u(xn , t)
+
ρh h h
which converges to the second partial derivative of u w.r.t x as h → 0. The
h → 0 is the limit case of the N objects we started with. Therefore the
vibrating string system is governed by the equation
∂ 2u T ∂ 2u
=
∂t2 ρ ∂x2
where T is the tension and ρ is the density of the string. Equivalently,
∂ 2u 2
2∂ u
= c (4.1.1)
∂t2 ∂x2
where c2 = T /ρ, c > 0, on x ∈ (0, L) and t > 0.
CHAPTER 4. WAVE EQUATION 83
Exercise 18. Show that the wave equation (4.1.1), derived above can be
written as
uzz = uww in (w, z) ∈ (0, L) × (0, ∞).
under a new coordinate system (w, z). One may, in fact choose coordinate
such that the string is fixed between (0, π).
Proof. Set w = x/a and z = t/b, where a and b will be chosen appropriately.
Then, wx = 1/a and zt = 1/b. Therefore, ux = uw /a, ut = uz /b, a2 uxx = uww
and b2 utt = uzz . Choosing a = 1 and b = 1/c. Choosing a = L/π and
b = L/cπ makes the domain (0, π).
ux = u w w x + u z z x = u w + u z
ut = uw wt + uz zt = c(uw − uz )
uxx = uww + 2uzw + uzz
utt = c2 (uww − 2uzw + uzz )
Now that we have derived the general form of the solution of wave equa-
tion, we return to understand the physical system of a vibrating infinite
string. The initial shape (position at initial time t = 0) of the string is given
as u(x, 0) = g(x), where the graph of g on R2 describes the shape of the
string. Since we need one more data to identify the arbitrary functions, we
also prescribe the initial velocity of the string, ut (x, 0) = h(x).
Another interesting property that follows from the general solution is
that for any four points A, B, C and D that form a rectangle bounded by
characteristic curves in R × R+ then u(A) + u(C) = u(B) + u(D) because
u(A) = F (α) + G(β), u(C) = F (γ) + G(δ), u(B) = F (α) + G(δ) and u(D) =
F (γ) + G(β).
Theorem 4.1.1. Given g ∈ C 2 (R) and h ∈ C 1 (R), there is a unique C 2
solution u of the Cauchy initial value problem (IVP) of the wave equation,
utt (x, t) − c2 uxx (x, t) = 0 in R × (0, ∞)
u(x, 0) = g(x) in R (4.1.2)
ut (x, 0) = h(x) in R,
g ′ (x) + h(x)/c. Similarly, 2G′ (x) = g ′ (x) − h(x)/c. Integrating2 both these
equations, we get
Z
1 1 x
F (x) = g(x) + h(y) dy + c1
2 c 0
2
assuming they are integrable and the integral of their derivatives is itself
CHAPTER 4. WAVE EQUATION 85
and Z x
1 1
G(x) = g(x) − h(y) dy + c2 .
2 c 0
Since F (x) + G(x) = g(x), we get c1 + c2 = 0. Therefore, the solution to the
wave equation is given by (4.1.3).
Aliter. Let us derive the d’Alembert’s formula in an alternate way. Note
that the wave equation can be factored as
∂ ∂ ∂ ∂
+c −c u = utt − c2 uxx = 0.
∂t ∂x ∂t ∂x
∂ ∂
We set v(x, t) = ∂t − c ∂x u(x, t) and hence
vt (x, t) + cvx (x, t) = 0 in R × (0, ∞).
Notice that the above first order PDE obtained is in the form of homogeneous
linear transport equation (cf. (2.2.1)), which we have already solved. Hence,
for some smooth function φ,
v(x, t) = φ(x − ct)
and φ(x) := v(x, 0). Using v in the original equation, we get the inhomoge-
neous transport equation,
ut (x, t) − cux (x, t) = φ(x − ct).
Recall the formula for inhomogenoeus transport equation (cf. (2.2.2))
Z t
u(x, t) = g(x − at) + φ(x − a(t − s), s) ds.
0
Since u(x, 0) = g(x) and a = −c, in our case the solution reduces to,
Z t
u(x, t) = g(x + ct) + φ(x + c(t − s) − cs) ds
0
Z t
= g(x + ct) + φ(x + ct − 2cs) ds
0
Z
−1 x−ct
= g(x + ct) + φ(y) dy
2c x+ct
Z
1 x+ct
= g(x + ct) + φ(y) dy.
2c x−ct
CHAPTER 4. WAVE EQUATION 86
But φ(x) = v(x, 0) = ut (x, 0) − cux (x, 0) = h(x) − cg ′ (x) and substituting
this in the formula for u, we get
Z
1 x+ct
u(x, t) = g(x + ct) + (h(y) − cg ′ (y)) dy
2c x−ct
1
= g(x + ct) + (g(x − ct) − g(x + ct))
Z x+ct 2
1
+ h(y) dy
2c x−ct
Z
1 1 x+ct
= (g(x − ct) + g(x + ct)) + h(y) dy
2 2c x−ct
and
Ω2 = {(x, t) | x < ct}.
For (x, t) ∈ Ω1 , the solution
Z x+ct
1 1
u1 (x, t) = (g(x − ct) + g(x + ct)) + h(y) dy.
2 2c x−ct
and, hence,
By setting (
u1 (x, t) (x, t) ∈ Ω1
u(x, t) =
u2 (x, t) (x, t) ∈ Ω2
and the fact that alll derivatives of u are continuous along ct = x line, due
to the compatibility condition implies that u is a solution (4.1.4).
Corollary 4.1.3 (Dirichlet Condition). Given g ∈ C 2 [0, L], h ∈ C 1 [0, L] and
φ, ψ ∈ C 2 (0, ∞), there is a unique C 2 solution u of the homogeneous Cauchy
initial value problem (IVP) of the wave equation,
2
utt (x, t) − c uxx (x, t) = 0
in (0, L) × (0, ∞)
u(x, 0) = g(x) in [0, L]
ut (x, 0) = h(x) in [0, L] (4.1.5)
u(0, t) = φ(t) in (0, ∞)
u(L, t) = ψ(t) in (0, ∞),
where φ, ψ, g, h satisfies the compatibility condition
g(0) = φ(0), g ′′ (0) = φ′′ (0), h(0) = φ′ (0)
and
g(L) = ψ(0), g ′′ (L) = ψ ′′ (0), h(L) = ψ ′ (0).
Proof. We first consider case φ = ψ ≡ 0, homogeneous Dirichlet conditions.
We extend g as an odd function on [−L, L] by setting g(−x) = −g(x) for
x ∈ [0, L]. Therefore, g is an odd 2L-periodic function in R. Similarly,
we extend h and we have g ∈ C 2 (R) and h ∈ C 1 (R). Then, we have a
unique solution u ∈ C 2 solving the Cauchy problem in R × (0, ∞). But
note that v(x, t) = −u(−x, t) is also a solution to the same Cauchy problem.
Thus, u(x, t) = −u(−x, t) for all (x, t) ∈ R × R+ and u(0, t) = 0. Similarly,
w(x, t) = −u(2L − x, t) is also a solution and, hence, u(x, t) = −u(2L − x, t)
which implies u(L, t) = 0.
Consider the lines ct = x and ct = −x + cL, then we will obtain u in the
four regions as u1 , u2 , u3 and u4 . Then follow the idea similar to proof in
above theorem.
values between this interval. The interval [x − ct, x + ct] is called the domain
of dependence. Thus, the region of R × (0, ∞) on which the value of u(x, t)
depends forms a triangle with base [x − ct, x + ct] and vertex at (x, t). The
domain of dependence of (x, t) is marked in x-axis by the characteristic curves
passing through (x, t).
Given a point p on the x-axis what values of u on (x, t) will depend on
the value of g(p) and h(p). This region turns out to be a cone with vertex
at p and is called the domain of influence. The domain of influence is the
region bounded by the characteristic curves passing through p.
If the initial data g and h are supported in the interval Bx0 (R) then the
solution u at (x, t) is supported in the region Bx0 (R + ct). Consequently, if g
and h have compact support then the solution u has compact support in R
for all time t > 0. This phenomenon is called the finite speed of propagation.
The unit outward normal at each point of the boundary ∂T (x, t) is given by
ν = (ν1 , ν2 ) defined as
(0, −1)
(y, s) ∈ T0
1
ν(y, s) = √1+c2 (1, c) (y, s) ∈ T+
√ 1
1+c2
(−1, c) (y, s) ∈ T− .
Note that the second and third integral are just the direction derivatives of
u along the direction (1, c) and (−1, c) in the line T+ and T− , respectively.
Therefore, we have
Z Z x+ct
f ((y, s) dy ds = − h(y) dy + cu(x, t) − cu(x + ct, 0)
T (x,t) x−ct
+c[u(x, t) − u(x − ct, 0)]
1
u(x, t) = [g(x + ct) + g(x − ct)]
2 " #
Z x+ct Z t Z x+c(t−s)
1
+ h(y) dy + f (y, s) dy ds .
2c x−ct 0 x−c(t−s)
CHAPTER 4. WAVE EQUATION 91
with the initial condition G(x) := U (x, 0) = (g(x), h(x)). The solution
U (x, t) is given as (cf. Appendix G)
Z t
U (x, t) = S(t)G(x) + S(t − s)F (s) ds
0
u(0, t) = u(L, t) = 0.
We are also given the initial position u(x, 0) = g(x) (at time t = 0) and initial
velocity of the string at time t = 0, ut (x, 0) = h(x). Given g, h : [0, L] → R
such that g(0) = g(L) = 0 and h(0) = h(L), we need to solve the initial value
problem (4.1.5) with φ = ψ ≡ 0.
CHAPTER 4. WAVE EQUATION 92
Let us seek for solutions u(x, t) whose variables can be separated. Let
u(x, t) = v(x)w(t). Differentiating and substituting in the wave equation, we
get
v(x)w′′ (t) = c2 v ′′ (x)w(t)
Hence
w′′ (t) v ′′ (x)
= .
c2 w(t) v(x)
Since RHS is a function of x and LHS is a function t, they must equal a
constant, say λ. Thus,
v ′′ (x) w′′ (t)
= 2 = λ.
v(x) c w(t)
Using the boundary condition u(0, t) = u(L, t) = 0, we get
v(0)w(t) = v(L)w(t) = 0.
such that α = (c1 + c2 )/2 and β = (c1 − c2 )/2. Using the boundary condition
v(0) = 0, we get c1 = 0 and hence
√
v(x) = c2 sinh( λx).
√
Now using v(L) = 0, we have c2 sinh λL = 0. Thus, c2 = 0 and v(x) = 0.
We have seen this cannot be a solution.
CHAPTER 4. WAVE EQUATION 93
√
Finally, if λ < 0, then set ω = −λ. We need to solve the simple
harmonic oscillator problem
′′
v (x) + ω 2 v(x) = 0 x ∈ (0, L)
v(0) = v(L) = 0.
for some constant bk and λk = −(kπ/L)2 . It now remains to solve w for each
of these λk . For each k ∈ N, we solve for wk in the ODE
and Z t
z(x, t) = wfs (x, t − s) ds,
0
where fs = f (·, s), are solutions to (4.2.2) and (4.2.4).
Proof. Since wg solves (4.2.3) with h = g, we have
∂
v = wtg (x, t) = (wg ) = 0
∂t
on Rn × (0, ∞). Further, v(x, 0) = wtg (x, 0) = g(x) in Rn and vt (x, 0) =
wttg (x, 0) = c2 ∆wg (x, 0) = 0. Thus, v solves (4.2.2). Now, let us denote
CHAPTER 4. WAVE EQUATION 96
w(x, t) = wfs (x, t−s), for 0 < s < t, which satisfies (4.2.3) with h(·) = f (c,̇s).
Then,
Z t Z t
fs
fs
zt (x, t) = w (x, 0) + wt (x, t − s) ds = wtfs (x, t − s) ds
0 0
and
Z t
ztt (x, t) = wtfs (x, 0) + wttfs (x, t − s) ds
0
Z t
2
= f (x, t) + c ∆ wfs (x, t − s) ds
0
= f (x, t) + c2 ∆z.
as long as, the domain of dependence [r − ct, r + ct] is in (0, ∞). Otherwise,
we extend G and H to entire R as
G(x; r) r>0
G̃(x; r) = g(x) r=0
G(x; −r) r < 0
and
H(x; r) r>0
H̃(x; r) = h(x) r=0
H(x; −r) r < 0.
Invoking Lemma 4.2.3, we get M̃ (x; r, t) as
Z r+ct
1 1
(r − ct)G̃(r − ct) + (r + ct)G̃(r + ct) + y H̃(y) dy
2r 2rc r−ct
Since
u(x, t) = lim+ M̃ (x; r, t),
r→0
d
u(x, t) = (tG(ct)) + tH(ct) = ctG′ (ct) + G(ct) + tH(ct).
dt
Thus, rG̃, rH̃ ∈ C(R). Since G, H ∈ C 1 (R), rG̃, rH̃ ∈ C 1 (R). Further, since
G′′ and H ′′ are bounded as r → 0+, rG̃, rH̃ ∈ C 2 (R).
The domain of dependence of (x, t0 ) for the three dimensional wave equa-
tion is the boundary of the three dimensional sphere with radius ct.
The Hadamard’s method of descent is the technique of finding a solu-
tion of the two dimensional wave equation using the three dimensional wave
equation.
where p
ρ= (ξ − x)2 + (η − y)2 .
CHAPTER 4. WAVE EQUATION 100
Exercise 19. If n is odd, show that the correct form of M that satisfies the
one dimensional wave equation is
n−3
1 ∂ 2
(rn−2 M (r, t)).
r ∂r
Proof. The Poisson’s formula corresponding to the three dimension case gives
the formula for
Z
fs 1
w (x, t − s) = 2
f (y, s) dσy .
4πc (t − s) Sc(t−s) (x)
CHAPTER 4. WAVE EQUATION 102
and
Z t Z tZ
fs 1 f (y, s)
u(x, t) = w (x, t − s) ds = 2
dσy ds
0 4πc 0 Sc(t−s) (x) t − s
Z ct Z
1 f (y, t − τ /c)
= 2
dσy dτ [setting τ = c(t − s)]
4πc 0 Sτ (x) τ
Z
1 f (y, t − r/c)
= 2
dy
4πc Bct (x) r
where r = |x − y|. Similarly, one can derive the formulae for one and two
dimensions.
Note that in the three dimensional case the integrand is not taken at time
t, but at an earlier time t − r/c. Thus, the integrand in this case is called
retarded potential.
Example 4.1. Consider the wave equation
utt (x, t) − c2 uxx (x, t) = sin 3x in (0, π) × (0, ∞)
u(0, t) = u(π, t) = 0 in (0, ∞)
u(x, 0) = ut (x, 0) = 0 in (0, π).
We look for the solution of the homogeneous wave equation
wtt (x, t) − c2 wxx (x, t) = 0 in (0, π) × (0, ∞)
w(0, t) = w(π, t) = 0 in (0, ∞)
w(x, 0) = 0 in (0, π)
wt (x, 0) = sin 3x in (0, π).
By separation of variable technique, we know that the general solution of w
is ∞
X
w(x, t) = [ak cos(kct) + bk sin(kct)] sin(kx)
k=1
and ∞
X
w(x, 0) = ak sin(kx) = 0.
k=1
Thus, ak = 0, for all k ∈ N. Also,
∞
X
wt (x, 0) = bk ck sin(kx) = sin 3x.
k=1
CHAPTER 4. WAVE EQUATION 103
Note that, for all λ ∈ R, zero is a trivial solution of the Laplacian. Thus,
we are interested in non-zero λ’s for which the Laplacian has non-trivial
solutions. Such an λ is called the eigenvalue and corresponding solution uλ
is called the eigen function.
Note that if uλ is an eigen function corresponding to λ, then αuλ , for all
α ∈ R, is also an eigen function corresponding to λ. Let W be the real vector
space of all u : Ω → R continuous (smooth, as required) functions such that
u(x, y) = 0 on ∂Ω. For each eigenvalue λ of the Laplacian, we define the
subspace of W as
Though the above theorem assures the existence of eigenvalues for Lapla-
cian, it is usually difficult to compute them for a given Ω. In this course, we
shall compute the eigenvalues when Ω is a 2D-rectangle and a 2D-disk.
4.3.1 In Rectangle
Let the rectangle be Ω = {(x, y) ∈ R2 | 0 < x < a, 0 < y < b}. We wish to
solve the Dirichlet EVP in the rectangle Ω
−∆u(x, y) = λu(x, y) (x, y) ∈ Ω
u(x, y) = 0 (x, y) ∈ ∂Ω.
Hence
v ′′ (x) w′′ (y)
− =λ+ .
v(x) w(y)
Since LHS is function of x and RHS is function y and are equal they must
be some constant, say µ. We need to solve the EVP’s
4.3.2 In Disk
Let the disk of radius a be Ω = {(x, y) ∈ R2 | x2 + y 2 < a2 }. We wish to
solve the Dirichlet EVP in the disk Ω
−1 ∂ 1 ∂2u
r ∂r r ∂u
∂r
− r2 ∂θ2 = λu(r, θ) (r, θ) ∈ Ω
u(θ) = u(θ + 2π) θ∈R
u(a, θ) = 0 θ ∈ R.
Note that this none other than the Bessel’s equation. We already know that
for each k ∈ N ∪ {0}, we have the Bessel’s function Jk as a solution to the
Bessel’s
√ equation. Recall
√ the boundary condition on v, v(a) = 0. Thus,
y( λa) = 0. Hence λa should be a zero of the Bessel’s function.√
For each k ∈ N∪{0}, let zkl be the l-th zero of Jk , l ∈ N. Hence λa = zkl
2
and so λkl = zkl /a2 and y(x) = Jk (x). Therefore, v(r) = Jk (zkl r/a). For each
k ∈ N ∪ {0} and l ∈ N, we have
Heat Equation
5.1 Derivation
Let a homogeneous material occupy a region Ω ⊂ Rn with C 1 boundary. Let
k denote the thermal conductivity (dimensionless quantity) and c be the heat
capacity of the material. Let u(x, t) be a function plotting the temperature
of the material at x ∈ Ω in time t. The thermal energy stored at x ∈ Ω, at
time t, is cu(x, t). If v(x, t) denotes the velocity of (x, t), by Fourier law, the
thermal energy changes following the gradients of temperature, i.e.,
cu(x, t)v(x, t) = −k∇u.
The thermal energy is the quantity that is conserved (conservation law) and
satisfies the continuity equation (1.4.1). Thus, we have
k
ut − ∆u = 0.
c
If the material occupying the region Ω is non-homogeneous, anisotropic,
the temperature gradient may generate heat in preferred directions, which
themselves may depend on x ∈ Ω. Thus, the conductivity of such a material
at x ∈ Ω, at time t, is given by a n × n matrix K(x, t) = (kij (x, t)). Thus, in
this case, the heat equation becomes,
1
ut − div K∇u = 0.
c
The heat equation is an example of a second order equation in divergence
form. The heat equation gives the temperature distribution u(x, t) of the
107
CHAPTER 5. HEAT EQUATION 108
where g ∈ C(Ω). We end this section with a remark that under the transfor-
mation t 7→ −t, in contrast to the wave equation, the heat equation changes
to a background equation. This signifies that the heat equation describes
irreversible process, i.e., it is not possible to find the distribution of temper-
ature at an earlier time t < t0 , if the distribution is given at t0 .
(ii) (Neumann condition) ∇x u(x, t) · ν(x) = h(x, t), where nu(x) is the unit
outward normal of (x, t) ∈ ∂Ω × (0, T );
where σ(x) is the specific heat at x, ρ(x) is density of bar at x and κ(x) is
the thermal conductivity of the bar at x. If the bar is homogeneous, i.e, its
CHAPTER 5. HEAT EQUATION 109
and
v ′′ (x) = λv(x).
But we already know how to solve the eigenvalue problem involving v. For
each k ∈ N, we have the pair (λk , vk ) as solutions
to the EVP involving v,
2 2 kπx
where λk = −(kπ) /L and vk (x) = sin L some constants bk . For each
k ∈ N, we solve for wk to get
ln wk (t) = λk c2 t + ln α
2
where α is integration constant. Thus, wk (t) = αe−(kcπ/L) t . Hence,
kπx −(kcπ/L)2 t
uk (x, t) = vk (x)wk (t) = βk sin e ,
L
for some constants βk , is a solution to the heat equation. By superposition
principle, the general solution is
X∞ X∞
kπx −(kcπ/L)2 t
u(x, t) = uk (x, t) = βk sin e .
k=1 k=1
L
We now use the initial temperature on the circle to find the particular solu-
tion. We are given u(θ, 0) = g(θ). Thus,
∞
a0 X
g(θ) = u(θ, 0) = + [ak cos(kθ) + bk sin(kθ)]
2 k=1
Z π
1
bk = g(θ) sin(kθ) dθ.
π −π
Exercise 20. Solve the heat equation for 2D Rectangle and 2D Disk
As a first step, for each s ∈ (0, ∞), consider w(x, t; s) as the solution of the
homogeneous problem (auxiliary)
s
wt (x, t) − c2 ∆ws (x, t) = 0 in Ω × (s, ∞)
ws (x, t) = 0 in ∂Ω × (s, ∞)
ws (x, s) = f (x, s) on Ω × {s}.
ut\ − ∆u = fˆ
c t) = fˆ(ξ, t)
ubt (ξ, t) − ∆u(ξ,
X n
ût (ξ, t) − i2 ξj2 û(ξ, t) = fˆ(ξ, t)
j=1
where we have used the initial condition û(ξ, 0) = ĝ(ξ) of the ODE. Therefore,
by inverse Fourier formula,
Z Z t
−n/2 −|ξ|2 t iξ·x iξ·x−|ξ|2 (t−s) ˆ
u(x, t) = (2π) ĝ(ξ)e e + e f (ξ, s) ds dξ
Rn 0
Z Z
−n −iξ·y 2
= (2π) g(y)e dy e−|ξ| t eiξ·x dξ
Rn Rn
Z Z t Z
−n −iξ·y 2
+(2π) f (y, s)e dy eiξ·x−|ξ| (t−s) ds dξ
n n
Z R 0 Z R
−n iξ·(x−y)−|ξ|2 t
= (2π) g(y) e dξ dy
Rn Rn
Z Z t Z
−n iξ·(x−y)−|ξ|2 (t−s)
+(2π) f (y, s) e dξ ds dy
Rn 0 Rn
Z Z Z t
= g(y)K(x, y, t) dy + f (y, s)K(x, y, t − s) ds dy
Rn Rn 0
where Z
−n 2
K(x, y, t) = (2π) eiξ·(x−y)−|ξ| t dξ.
Rn
Note that
2
√
(x − y) √ (x − y) |x − y|2
iξ · (x − y) − |ξ| t = − ξ t − i √ · ξ t−i √ −
2 t 2 t 4t
CHAPTER 5. HEAT EQUATION 115
√ √
and, set η = ξ t − i (x−y)
√
2 t
. Therefore, dη = tdξ. Using this substituion
in K and simplifying, we get
2 /4t
K(x, y, t) = (4πt)−n/2 e−|x−y|
called the heat kernel or the fundamental solution of heat equation. The
function K can be motivated in another way. Note that if u(x, t) is a solu-
tion of the heat equation, then (u ◦ Tλ )(x, t) is also a solution of the heat
equation, where Tλ (x, t) = (λx, λ2 t) is a linear transformation for any λ 6= 0.
This scaling or dilation is called the parabolic scaling. Thus, we look for a
solution u(x, t) = v(t)w(r2 /t), where r = |x|. Substituting this separation
2
of variable in the heat equation, we derive v(t) = t−n/2 and w(t) = e−r /4t .
This motivates us to define the fundamental solution as
(
1 n/2 −r 2 /4t
− 4πt e x ∈ Rn , t > 0
K(x, t) :=
0 x ∈ Rn , t < 0.
CHAPTER 5. HEAT EQUATION 116
Chapter 6
The Laplacian
117
CHAPTER 6. THE LAPLACIAN 118
Therefore, ∇x Ru = Ot ∇y u and
(∆ ◦ R)u(x) = ∇x Ru · ∇x (Ru) = Ot ∇y u · Ot ∇y u = OOt ∇y u · ∇y u = ∆y u.
But ∆y u = (∆u)(Ox) = (R ◦ ∆)u(x).
A radial function is constant on every sphere about the origin. Since
Laplacian commutes with rotations, it should map the class of all radial
functions to itself.
In cartesian coordiantes, the n-dimensional Laplacian is given as
Xn
∂2
∆ := 2
.
i=1
∂x i
1 ∂ 2u 2
2 ∂ u
2
2 ∂ u ∂ 2u 1 ∂u
= sin θ + cos θ − 2 cos θ sin θ − .
r2 ∂θ2 ∂x2 ∂y 2 ∂x∂y r ∂r
Therefore,
∂ 2u 1 ∂ 2u ∂ 2u ∂ 2 u 1 ∂u
+ = + θ − .
∂r2 r2 ∂θ2 ∂x2 ∂y 2 r ∂r
and hence
∂ 2u 1 ∂ 2 u 1 ∂u
∆u = + + .
∂r2 r2 ∂θ2 r ∂r
where r ∈ [0, ∞), φ ∈ [0, π] (zenith angle or inclination) and θ ∈ [0, 2π)
(azimuth angle).
d2
Note that in one dimension, n = 1, ∆ = dx 2.
CHAPTER 6. THE LAPLACIAN 120
d2 v(r) (n − 1) dv(r)
∆u(x) = + .
dr2 r dr
Proof. Note that
p
∂r ∂|x| ∂( x21 + . . . + x2n )
= =
∂xi ∂xi ∂xi
1 2
= (x1 + . . . + x2n )−1/2 (2xi )
2
xi
= .
r
Thus,
Xn Xn
∂ ∂u(x) ∂ dv(r) xi
∆u(x) = =
i=1
∂xi ∂xi i=1
∂xi dr r
Xn
∂ 1 dv(r) n dv(r)
= xi +
i=1
∂xi r dr r dr
n
X x2 d dv(r) 1
i n dv(r)
= +
i=1
r dr dr r r dr
Xn
x2i 1 d2 v(r) 1 dv(r) n dv(r)
= 2
− 2 +
i=1
r r dr r dr r dr
r2 1 d2 v(r) 1 dv(r) n dv(r)
= 2
− 2 +
r r dr r dr r dr
2
d v(r) 1 dv(r) n dv(r)
= − +
dr2 r dr r dr
2
d v(r) (n − 1) dv(r)
= + .
dr2 r dr
Hence the result proved.
More generally, the Laplacian in Rn may be written in polar coordinates
as
∂2 n−1 ∂ 1
∆ := 2
+ + 2 ∆Sn−1
∂r r ∂r r
CHAPTER 6. THE LAPLACIAN 121
and uy (x, 0) = 0. Thus, the Cauchy problem has no solution unless g(x) = 0.
If g ≡ 0 then the solution is
2 3/2 2 3/2 2 3/2
u(x, y) = F x + y −F x− y +f x− y
3 3 3
(ii) (Neumann condition) ∇u·ν = g, where ν(x) is the unit outward normal
of x ∈ ∂Ω;
CHAPTER 6. THE LAPLACIAN 123
The second equality is called the compatibility condition. Thus, for a Neu-
mann problem the given data f, g must necessarily satisfy the compatibility
condition. Otherwise, the Neumann problem does not make any sense.
The aim of this chapter is solve
−∆u(x) = f (x) in Ω
one of the above inhomogeneous boudary condition on ∂Ω,
ux = v y and uy = −vx .
Conversely, for any holomorphic function f , its real part and imaginary part
are conjugate harmonic functions. This observation gives us more examples
of harmonic functions, for instance, since all complex polynomials f (z) = z m
are holomorphic we have (using the polar coordinates) u(r, θ) = rm cos mθ
and v(r, θ) = rm sin mθ are harmonic functions in R2 for all m ∈ N. Simi-
larly, since f (z) = log z = ln r + iθ is holomorphic in certain region, we have
u(r, θ) = ln r and v(r, θ) = θ are harmonic in R2 \ (0, 0) and R2 \ {θ = 0},
respectively.
Exercise 23. Show that there are infinitely many linearly independent har-
monic functions in the vector space C 2 (R2 ).
Two Dimensions
k
X
Pk (x, y) := ai xi y k−i
i=0
CHAPTER 6. THE LAPLACIAN 126
k
X
Qk (θ) = ai (cos θ)i (sin θ)k−i .
i=0
Three Dimensions
Consider a general homogeneous polynomial
X
Pk (x) := aα x α
|α|=k
(k+2)(k+1)
of degree k in R3 (three variables). Note that Pk contains k+2 k
= 2
coefficients. Then ∆Pk (x) is a homogeneous polynomial of degree k − 2 and,
hence, contains k(k−1)
2
coefficients. If Pk ∈ Hk (R3 ), i.e. ∆Pk (x) = 0, then all
the k(k−1)
2
coefficients should vanish. Thus, we have k(k−1) 2
equations relating
(k+2)(k+1)
the 2
coefficients of Pk and, hence, Hk (R3 ) is of dimension
(k + 2)(k + 1) − k(k − 1)
= 2k + 1.
2
The basis of the 2k +1 dimensional space Hk (R3 ) is given in terms of the Leg-
endre functions which we shall describe now. In spherical coordinates, x =
r sin φ cos θ, y = r sin φ sin θ and z = r cos φ. Thus, Pk (r, φ, θ) = rk R(φ)Q(θ)
where
X
R(φ)Q(θ) = aα (sin φ)α1 +α2 (cos φ)α3 (cos θ)α1 (sin θ)α2 .
|α|=k
dR d2 R dR
R′ (φ) = − sin φ and R′′ (φ) = sin2 φ 2 − cos φ
dw dw dw
In the new variable w, we get the Legendre equation
2 ′′ ′ m2
(1 − w )R (w) − 2wR (w) = − k(k + 1) R(w) w ∈ [−1, 1].
1 − w2
For each k ∈ N ∪ {0}, this has the Legendre polynomials, Rk,m (cos φ), as its
solutions. Therefore, in general,
However, we are interested only those Rk,m which gives a polynomial of degree
k in R3 . Thus, for m = 0, 1, . . . , k,
dk+m
Rk,m (w) = (1 − w2 )m/2 k+m
(1 − w2 )k .
dw
Note that, for each fixed k and all 1 ≤ m ≤ k, the collection
{Rk,0 (cos φ), cos mθRk,m (cos φ), sin mθRk,m (cos φ)} ⊂ Hk (R3 )
is also a zero set. The function Rk,0 (cos φ) and its constant multiples are
called zonal harmonics.
If 0 < m < k, then the spherical harmonics is of the form
dk+m
(α cos mθ + β sin mθ) sinm φ (1 − w2 )k .
dwk+m
If the first term is zero then tan mθ = −α/β. This corresponds to great
circle through the north pole and south pole of S2 and the angle between
the planes containing two consecutive great circle is π/m. The second term
vanishes on φ = 0 and φ = π corresponding to the north and south pole,
respectively. The third term vanishes on k − m latitude circle. Thus, we have
orthogonally intersecting family of circles which form the zero set which are
called tesseral harmonics.
If m = k then the spherical harmonics is of the form
Higher Dimensions
Consider a general homogeneous polynomial
X
Pk (x) := aα x α
|α|=k
of degree k in Rn (n variables). Note that Pk contains n+k−1 k
coefficients.
Then ∆Pk (x) is a homogeneous polynomial of degree k − 2 and, hence, con-
n+k−3 n
tains k−2 coefficients. If Pk ∈ Hk (R ), i.e. ∆Pk (x) = 0, then all the
n+k−3
n+k−3
k−2
coefficients should vanish. Thus, we have k−2
equations relating
n+k−1 n
k
coefficients of Pk and, hence, Hk (R ) is of dimension
n+k−1 n+k−3
ℓ := − .
k k−2
CHAPTER 6. THE LAPLACIAN 130
2π n/2
Definition 6.4.2. Let Ω be an open subset of Rn and wn = Γ(n/2)
(cf. Ap-
pendix E) be the surface area of the unit sphere S1 (0) of Rn .
(a) A function u ∈ C(Ω) is said to satisfy the first mean value property
(I-MVP) in Ω if
Z
1
u(x) = u(y) dσy for any Br (x) ⊂ Ω.
ωn rn−1 Sr (x)
(b) A function u ∈ C(Ω) is said to satisfy the second mean value property
(II-MVP) in Ω if
Z
n
u(x) = u(y) dy for any Br (x) ⊂ Ω.
ωn rn Br (x)
CHAPTER 6. THE LAPLACIAN 131
Exercise 25. Show that the first MVP and second MVP are equivalent. That
is show that u satisfies (a) iff u satisfies (b).
Owing to the above exercise we shall, henceforth, refer to the I-MVP and
II-MVP as just mean value property (MVP).
We shall now prove a result on the smoothness of a function satisfying
MVP.
Theorem 6.4.3. If u ∈ C(Ω) satisfies the MVP in Ω, then u ∈ C ∞ (Ω).
Proof. We first consider uε := ρε ∗ u, the convolution of u with mollifiers, as
introduced in Theorem F.0.10. where
Ωε := {x ∈ Ω | dist(x, ∂Ω) > ε}.
We shall now show that u = uε for all ε > 0, due to the MVP of u and the
radial nature of ρ. Let x ∈ Ωε . Consider
Z
uε (x) = ρε (x − y)u(y) dy
ZΩ
= ρε (x − y)u(y) dy (Since supp(ρε ) is in Bε (x))
Bε (x)
Z ε Z
= ρε (r) u(y) dσy dr (cf. Theorem E.0.8)
0 Sr (x)
Z ε
= u(x)ωn ρε (r)rn−1 dr (Using MVP of u)
Z ε0 Z
= u(x) ρε (r) dσy dr
0 Sr (0)
Z
= u(x) ρε (y) dy = u(x).
Bε (0)
Thus uε (x) = u(x) for all x ∈ Ωε and for all ε > 0. Since uε ∈ C ∞ (Ωε ) for
all ε > 0 (cf. Theorem F.0.10), we have u ∈ C ∞ (Ωε ) for all ε > 0.
CHAPTER 6. THE LAPLACIAN 132
Thus, v(r) = u(x) for all r > 0 and hence v is a constant function of r and
v ′ (s) = 0. But
Z
dv(s) 1
0= = ∆u(y) dy > 0
dr ωn rn−1 Bs (x)
S := {x ∈ Ω | u(x) = M }.
CHAPTER 6. THE LAPLACIAN 134
Hence equality will hold above only when u(y) = M for all y ∈ Br (x). Thus,
we have shown that for any x ∈ S, we have Br (x) ⊂ S. Therefore, S is open.
Since Ω is connected, the only open and closed subsets are ∅ or Ω. Since
S was assumed to be non-empty, we should have S = Ω. Thus, u ≡ M is
constant in Ω.
Corollary 6.4.8 (Weak maximum Principle). Let Ω be an open, bounded
subset of Rn . Let u ∈ C(Ω) be harmonic in Ω. Then
max u(y) = max u(y).
y∈Ω y∈∂Ω
The above inequality is true for all ε > 0. Thus, u(x) ≤ maxx∈∂Ω u(x), for all
x ∈ Ω. Therefore, maxΩ u ≤ maxx∈∂Ω u(x). and hence we have equality.
4
v ∈ C 2 (a, b) has a local maximum at x ∈ (a, b) then v ′ (x) = 0 and v ′′ (x) ≤ 0
CHAPTER 6. THE LAPLACIAN 135
Hence
C1
|Dα u(x)| ≤ kuk1,Br (x) for |α| = 1.
rn+1
Let now k ≥ 2 and α be a multi-index such that |α| = k. We assume the
induction hypothesis that the estimate to be proved is true for k − 1. Note
βu
that Dα u = ∂D ∂xi
for some i ∈ {1, 2, . . . , n} and |β| = k − 1. Moreover,
if u is harmonic then by differentiating the Laplace equation and using the
βu
equality of mixed derivatives, we have ∂D ∂xi
is harmoic for i = 1, 2, . . . , n.
Thus, following an earlier argument, we have
β n Z β
∂D u(x) nk ∂D u(y)
|Dα u(x)| = = dy
∂xi ωn rn Br/k (x) ∂xi
Z
nk n β
= D uν i dσ y
ωn r Sr/k (x)
n
nk β
≤ kD uk∞,Sr/k (x) .
r
It now only remains to estimate kDβ uk∞,Sr/k (x) . Let z ∈ Sr/k (x), then
B(k−1)r/k (z) ⊂ Br (x) ⊂ Ω. But, using induction hypothesis for k − 1, we
have
Ck−1 nk n+k
|Dα u(x)| ≤ kuk1,Br (x)
(k − 1)n+k−1 rn+k
C0 2(n+1)(k−1) nk (k − 1)k−1 k n+k
= kuk1,Br (x)
(k − 1)n+k−1 rn+k
n
C0 (2n+1 nk)k k 1
= kuk1,Br (x)
rn+k k−1 2n+1
n
C0 (2n+1 nk)k k 1
= n+k
kuk1,Br (x)
r 2(k − 1) 2
n
Ck k 1
≤ n+k
kuk1,Br (x) since ≤ 1.
r 2(k − 1) 2
CHAPTER 6. THE LAPLACIAN 137
Hence
Ck
|Dα u(x)| ≤ kuk1,Br (x) for |α| = k, ∀k ≥ 2.
rn+k
Proof. For any x ∈ Rn and r > 0, we have the estimate on the first derivative
as,
C1
|∇u(x)| ≤ kuk1,Br (x)
rn+1
2n+1 n
= kuk1,Br (x)
ωn rn+1
2n+1 n
≤ n+1
kuk∞,Rn ωn rn
ωn r
2n+1
= kuk∞,Rn → 0 as r → ∞.
r
Thus, ∇u ≡ 0 in Rn and hence u is constant.
In particular,
1
u(y) ≤ u(x) ≤ Cu(y) ∀x, y ∈ ω.
C
CHAPTER 6. THE LAPLACIAN 138
Thus, 1/2n u(y) ≤ u(x). Interchanging the role of x and y, we get 1/2n u(x) ≤
u(y). Thus, 1/2n u(y) ≤ u(x) ≤ 2n u(y) for all x, y ∈ ω such that |x − y| ≤ r.
Now, let x, y ∈ ω. Since ω is compact and connected in Ω, we can pick
points x = x1 , x2 , . . . , xm = y such that ∪m
i=1 Bi ⊃ ω, where Bi := Br/2 (xi )
and are sorted such that Bi ∩ Bi+1 6= ∅, for i = 1, 2, . . . , m − 1. Hence, note
that |xi+1 − xi | ≤ r. Therefore,
1 1 1 1
u(x) = u(x1 ) ≥ n
u(x2 ) ≥ 2n u(x3 ) ≥ . . . ≥ (m−1)n u(xm ) = (m−1)n u(y).
2 2 2 2
1
Thus, C can be chosen to be 2(m−1)n
.
(b) (Stability) |u1 (x) − u2 (x)| ≤ maxy∈∂Ω |g1 (y) − g2 (y)| for all x ∈ Ω.
Proof. The fact that there is atmost one solution to the Dirichlet problem
follows from the Theorem 6.5.1. Let w = u1 − u2 . Then w is harmonic.
(a) Note that w = g1 − g2 ≥ 0 on ∂Ω. Since g1 (x0 ) > g2 (x0 ) for some
x0 ∈ ∂Ω, then w(x) > 0 for all x ∈ ∂Ω. This proves the comparison
result.
We remark that the uniqueness result is not true for unbounded domains.
Example 6.3. Consider the problem (6.5.1) with g ≡ 0 in the domain Ω =
{x ∈ Rn | |x| > 1}. Obviously, u = 0 is a solution. But we also have a
non-trivial solution (
ln |x| n=2
u(x) =
|x|2−n − 1 n ≥ 3.
Example 6.4. Consider the problem (6.5.1) with g ≡ 0 in the domain Ω =
{x ∈ Rn | xn > 0}. Obviously, u = 0 is a solution. But we also have a
non-trivial solution u(x) = xn .
We have shown above that if a solution exists for (6.5.1) then it is unique
(cf. Theorem 6.5.1). So the question that remains to be answered is on
the existence of solution of (6.5.1), for any given domain Ω. In the modern
theory, there are three different methods to address this question of existence,
viz., Perron’s Method, Layer Potential (Integral Equations) and L2 methods.
CHAPTER 6. THE LAPLACIAN 140
1. w is superharmonic in Ω ∩ U
Theorem 6.6.3. The Dirichlet problem (6.5.1) is solvable for any arbitrary
bounded domain Ω and for any arbitrary g on ∂Ω iff all the points in ∂Ω are
regular.
Proof. One way is obvious. If (6.5.1) is solvable and x0 ∈ ∂Ω. Then, the
solution to
∆w = 0 in Ω
w = g on ∂Ω,
where g(x) = |x − x0 |, is a barrier function at x0 . Thus, any x0 ∈ ∂Ω
is a regular point. The converse is proved using the Perron’s method for
subharmonic functions.
Lemma 6.6.5. If Ω satisfies the exterior sphere condition then all boundary
points of Ω are regular.
Theorem 6.6.6. Any bounded domain with C 2 boundary satisfies the exte-
rior sphere condition.
Definition 6.6.7. A bounded domain Ω ⊂ Rn is said to satisfy the exterior
cone condition if for every point x0 ∈ ∂Ω there is a finite right circular cone
K with vertex at x0 such that K ∩ Ω = x0 .
Exercise 28. Any domain satisfying the exterior sphere condition also satisfies
the exterior cone condition.
Exercise 29. Every bounded Lipschitz domain satisfies the exterior cone con-
dition.
Lemma 6.6.8. If Ω satisfies the exterior cone condition then all boundary
points of Ω are regular.
/ Ω | |x − x0 | ≤ λi }.
diverges, where Ci := cap2 {x ∈
CHAPTER 6. THE LAPLACIAN 142
√ √
λx
If λ > 0, then v(x) = αe + βe− λx
. Equivalently,
√ √
v(x) = c1 cosh( λx) + c2 sinh( λx)
such that α = (c1 + c2 )/2 and β = (c1 − c2 )/2. Using the boundary condition
v(0) = 0, we get c1 = 0 and hence
√
v(x) = c2 sinh( λx).
√
Now using v(a) = 0, we have c2 sinh λa = 0. Thus, c2 = 0 and v(x) = 0.
We have seen this cannot √be a solution.
If λ < 0, then set ω = −λ. We need to solve
′′
v (x) + ω 2 v(x) = 0 x ∈ (0, a)
(6.7.1)
v(0) = v(a) = 0.
The constant δk are obtained by using the boundary condition u(x, b) = h(x)
which yields
∞
X
kπb kπx
h(x) = u(x, b) = δk sinh sin .
k=1
a a
Since h(0) = h(a) = 0, we know that h admits a Fourier Sine series. Thus
δk sinh kπb
a
is the k-th Fourier sine coefficient of h, i.e.,
−1 Z a
kπb 2 kπx
δk = sinh h(x) sin .
a a 0 a
Proof. Given the nature of the domain, we shall use the Laplace operator in
polar coordinates,
1 ∂ ∂ 1 ∂2
∆ := r + 2 2
r ∂r ∂r r ∂θ
where r is the magnitude component and θ is the direction component. Then
∂Ω is the circle of radius one. Then, solving for u(x, y) in the Dirichlet
problem is to equivalent to finding U (r, θ) : Ω → R such that
1∂ 1 ∂2U
r ∂r r ∂U
∂r
+ r2 ∂θ2 = 0 in Ω
U (r, θ + 2π) = U (r, θ) in Ω (6.7.2)
U (R, θ) = G(θ) on ∂Ω
where U (r, θ) = u(r cos θ, r sin θ), G : [0, 2π) → R is G(θ) = g(cos θ, sin θ).
Note that both U and G are 2π periodic w.r.t θ. We will look for solution
CHAPTER 6. THE LAPLACIAN 145
U (r, θ) whose variables can be separated, i.e., U (r, θ) = v(r)w(θ) with both
v and w non-zero. Substituting it in the polar form of Laplacian, we get
w d dv v d2 w
r + 2 2 =0
r dr dr r dθ
and hence
−r d dv 1 d2 w
r = .
v dr dr w dθ2
Since LHS is a function of r and RHS is a function of θ, they must equal a
constant, say λ. We need to solve the eigen value problem,
′′
w (θ) − λw(θ) = 0 θ∈R
w(θ + 2π) = w(θ) ∀θ.
To find the constants, we must use U (R, θ) = G(θ). If G ∈ C 1 [0, 2π], then G
admits Fourier series expansion. Therefore,
∞
a0 X k
G(θ) = + R ak cos(kθ) + Rk bk sin(kθ)
2 k=1
where Z π
1
ak = k G(θ) cos(kθ) dθ,
R π −π
Z π
1
bk = k G(θ) sin(kθ) dθ.
R π −π
Using this in the formula for U and the uniform convergence of Fourier series,
we get
Z " #
1 X r k
∞
1 π
U (r, θ) = G(η) + (cos kη cos kθ + sin kη sin kθ) dη
π −π 2 k=1 R
Z " #
1 X r k
∞
1 π
= G(η) + cos k(η − θ) dη.
π −π 2 k=1 R
CHAPTER 6. THE LAPLACIAN 147
k=1
R k=1
R 1 − Rr ei(η−θ)
2
R − rR cos(η − θ)
= −1
R2
+ r2 − 2rR cos(η − θ)
rR cos(η − θ) − r2
=
R2 + r2 − 2rR cos(η − θ)
in U (r, θ) we get
Z π
R2 − r 2 G(η)
U (r, θ) = dη.
2π −π R2 + r2 − 2rR cos(η − θ)
Note that the formula derived above for U (r, θ) can be rewritten in Carte-
sian coordinates and will have the form
Z
R2 − |x|2 g(y)
u(x) = 2
dy.
2πR SR (0) |x − y|
This can be easily seen, by setting y = R(x10 cos η +x20 sin η), we get dy = Rdη
and |x − y|2 = R2 + r2 − 2rR cos(η − θ). This is called the Poisson formula.
More generally, the unique solution to the Dirichlet problem on a ball of
radius R centred at x0 in Rn is given by Poisson formula
Z
R2 − |x − x0 |2 g(y)
u(x) = n
dy.
ωn R SR (x0 ) |x − y|
dw d2 w dw
w′ (φ) = − sin φ and w′′ (φ) = sin2 φ 2 − cos φ
dx dx dx
In the new variable x, we get the Legendre equation
We have already seen that this is a singular problem (while studying S-L
problems). For each k ∈ N ∪ {0}, we have the solution (wk , λk ) where
−∆u = f in Rn . (6.8.1)
v ′′ (r) (1 − n)
=
v ′ (r) r
′
ln v (r) = (1 − n) ln r + ln b
v ′ (r) = br(1−n)
Thus,
Z (
b
′ r
(2πr) for n = 2
1= v (r) dσ = 1−n n−1
Sr (0) br (r ωn ) for n ≥ 3.
The constant a can be chosen arbitrarly, but to keep things simple, we choose
a ≡ 0 for n ≥ 2. For convention sake, we shall add minus (“−”) sign (notice
the minus sign in (6.8.1)).
Definition 6.8.3. For any fixed x0 ∈ Rn We say K(x0 , x), defined as
(
1
− 2π ln |x − x0 | (n = 2)
K(x0 , x) := |x−x0 |2−n
ωn (n−2)
(n ≥ 3),
for all measurable subsets E of the measure space Rn . The Dirac measure
has the property that Z
dδx = 1
E
if x ∈ E and zero if x ∈
/ E. Also, for any integrable function f ,
Z
f (y) dδx = f (x).
Rn
−∆K(x0 , x) = δx0 in Rn .
CHAPTER 6. THE LAPLACIAN 153
Note that the above equation, as such, makes no sense because the RHS is a
set-function taking subsets of Rn as arguments, whereas K is a function on
Rn . To give meaning to above equation, one needs to view δx as a distribution
(introduced by L. Schwartz) and the equation should be interpreted in the
distributional derivative sense. The Dirac measure is the distributional limit
of the sequence of mollifiers, ρε , in the space of distributions.
In this section, we shall give a formula for the solution of the Poisson equation
(6.8.1) in Rn in terms of the fundamental solution.
Thus, for n = 2,
m2 1
|Im (x)| ≤ + | ln m| kD2 f k∞,Rn
2 2
and for n ≥ 3, we have
m2
|Im (x)| ≤ kD2 f k∞,Rn .
2(n − 2)
CHAPTER 6. THE LAPLACIAN 155
Thus, for n = 2,
|Jm (x)| ≤ m| ln m|k∇f k∞,Rn
and for n ≥ 3, we have
m
|Jm (x)| ≤ k∇f k∞,Rn .
(n − 2)
Hence, as m → 0, |Jm (x)| → 0. Now, to tackle the last term Km (x), we note
that a simple computation yields that ∇y K(y) = ωn−1 |y|n
y. Since we are in
the m radius sphere |y| = m. Also the unit vector ν outside of Sm (0), as a
boundary of Rn \ Bm (0), is given by −y/|y| = −y/m. Therefore,
1 1
∇y K(y) · ν = n+1
y·y = .
ωn m ωn mn−1
Z
Km (x) = − f (x − y)∇y K(y) · ν dσy
Sm (0)
Z
−1
= f (x − y) dσy
ωn mn−1 Sm (0)
Z
−1
= f (y) dσy
ωn mn−1 Sm (x)
Since f is continuous, for every ε > 0, there is a δ > 0 such that |f (x) −
f (y)| < ε whenever |x − y| < δ. When m → 0, we can choose m such that
m < δ and for this m, we see that Now, consider
Z
1
|Km (x) − (−f (x))| = f (x) − f (y) dσ y
ωn mn−1 Sm (x)
Z
1
= |f (x) − f (y)| dσy < ε.
ωn mn−1 Sm (x)
Thus, as m → 0, Km (x) → −f (x). Hence, u solves (6.8.1).
CHAPTER 6. THE LAPLACIAN 156
Remark 6.8.5. Notice that in the proof above, we have used the Green’s
identity eventhough our domain is not bounded (which is a hypothesis for
Green’s identity). This can be justified by taking a ball bigger than Bm (0)
and working in the annular region, and later letting the bigger ball approach
all of Rn .
Proof. We know that (cf. Theorem 6.8.4) u′ (x) := (K ∗ f )(x) solves (6.8.1),
the Poisson equation in Rn . Moreover, u′ is bounded for n ≥ 3, since K(x) →
0 as |x| → ∞ and f has compact support in Rn . Also, since u is given to be
a bounded solution of (6.8.1), v := u − u′ is a bounded harmonic function.
Hence, by Liouville’s theorem, v is constant. Therefore u = u′ + C, for some
constant C.
We turn our attention to studying Poisson equation in proper subsets of
R . Let Ω be an open bounded subset of Rn with C 1 boundary ∂Ω.
n
Proof. Let u and v be solutions of the Poisson equation with same boundary
conditions on ∂Ω. Then w := u − v is a harmonic function, ∆w = 0, with
homogeneous boundary condition on ∂Ω. By Green’s identity D.0.7, we have
Z Z
2
|∇w| dx = w(∇w · ν) dσ.
Ω ∂Ω
For the Drichlet, Neumann and Mixed case, the RHS is zero. For the Robin
condition the RHS is negative,
Z Z
w(∇w · ν) dσ = −c w2 dσ ≤ 0.
∂Ω ∂Ω
CHAPTER 6. THE LAPLACIAN 157
Lemma 6.8.8. Let f be bounded and locally Hölder continuous6 with expo-
nent γ ≤ 1 in Ω. Then u := K ∗ f ∈ C 2 (Ω), −∆u = f in Ω.
Theorem 6.8.9 (Existence). Let Ω be a bounded domain with all boundary
points being regular w.r.t Laplacian. The classical Dirichlet problem (6.8.2) is
solvable (hence uniquely) for any bounded, locally Hölder continuous function
f in Ω and continuous function g on ∂Ω.
Proof. Recall that K is a fundamental solution of −∆. Set w(x) := f ∗ K in
Rn then −∆w = f . Set v = u − w. Then (6.8.2) is solvable iff
−∆v = 0 in Ω
v = g − w on ∂Ω
any u ∈ C 2 (Ω) and vx (y) = K(y − x), the fundamental solution on Rn \ {x},
on the domain Ωm , we get
Z
u(y)∆y vx (y) dy
Ωm
Z Z
∂vx ∂u(y)
− vx (y)∆y u(y) dy = u(y) (y) − vx (y) dσy
Ωm ∂Ωm ∂ν ∂ν
Z Z
∂vx ∂u(y)
− vx (y)∆y u(y) dy = u(y) (y) − vx (y) dσy
Ωm ∂Ωm ∂ν ∂ν
Z Z Z Z
− + = +
Ω Bm (x) ∂Ω Sm (x)
Z
vx (y)∆y u(y) dy
Bm (x)
Z
∂vx
− u(y) (y) dσy
Sm (x) ∂ν
Z Z
∂u(y) ∂vx ∂u(y)
+ vx (y) dσy = u(y) (y) − vx (y) dσy
Sm (x) ∂ν ∂Ω ∂ν ∂ν
Z
+ vx (y)∆y u(y) dy
Z Ω
∂K ∂u(y)
Im (x) + Km (x) + Jm (x) = u(y) (y − x) − K(y − x) dσy
∂Ω ∂ν ∂ν
Z
+ K(y − x)∆y u(y) dy
Ω
Thus, (
m2 1
2 2
+ | ln m| kD2 uk∞,Ω for n = 2
|Im (x)| ≤ m2
2(n−2)
kD2 uk∞,Ω for n ≥ 3.
Hence, as m → 0, |Im (x)| → 0. Next, consider the term Km (x). Note that
−1
∇y K(y−x) = ωn |y−x|n (y−x). Since we are in the m radius sphere |y−x| = m.
CHAPTER 6. THE LAPLACIAN 159
Since u is continuous, for every ε > 0, there is a δ > 0 such that |u(x)−u(y)| <
ε whenever |x − y| < δ. When m → 0, we can choose m such that m < δ
and for this m, we see that Now, consider
Z
1
|Km (x) − (−u(x))| = u(x) − u(y) dσ y
ωn mn−1 Sm (x)
Z
1
= |u(x) − u(y)| dσy < ε.
ωn mn−1 Sm (x)
Thus, for n = 2,
(
m| ln m|k∇y uk∞,Ω for n = 2
|Jm (x)| ≤ m
|Jm (x)| ≤ (n−2) k∇y uk∞,Ω for n ≥ 3.
Now applying the second identity of Corollary D.0.7 for any u ∈ C 2 (Ω)
and v(y) = ψx (y), we get
Z Z
∂ψx ∂u
u − ψx dσy = (u∆y ψx − ψx ∆y u) dy.
∂Ω ∂ν ∂ν Ω
in (6.8.3), we get
Z Z
u(x) = (ψx (y) − K(y − x)) ∆y u dy + u∇ (ψx (y) − K(y − x)) · ν dσy .
Ω ∂Ω
The identity above motivates the definition of what is called the Green’s
function.
Proof. Let us fix x, y ∈ Ω. For a fixed m > 0, set Ωm = Ω \ (Bm (x) ∪ Bm (y))
and applying Green’s identity for v(·) := G(x, ·) and w(·) := G(y, ·), we get
Z Z
∂w(z) ∂v(z)
v(z) − w(z) dσz = v(z)∆z w(z) dz
∂Ωm ∂ν ∂ν Ωm
Z
− w(z)∆z v(z) dz
Ωm
Z
∂w(z) ∂v(z)
v(z) − w(z) dσz = 0
∂Ωm ∂ν ∂ν
Z Z
∂w ∂v ∂v ∂w
v −w dσz = w −v dσz
Sm (x) ∂ν ∂ν Sm (y) ∂ν ∂ν
Jm (x) + Km (x) = Jm (y) + Km (y).
CHAPTER 6. THE LAPLACIAN 162
Z
|Jm (x)| ≤ |v(z)∇z w(z) · ν| dσz
Sm (x)
Z
≤ k∇wk∞,Ω |v(z)| dσz
Sm (x)
Z
= k∇wk∞,Ω |ψx (z) − K(z − x)| dσz .
Sm (x)
Thus, for n = 2,
|Jm (x)| ≤ (2πmkψx k∞,Ω + m| ln m|) k∇wk∞,Ω
and for n ≥ 3, we have
n−1 m
|Jm (x)| ≤ ωn m kψx k∞,Ω + k∇wk∞,Ω .
(n − 2)
Hence, as m → 0, |Jm (x)| → 0. Now, consider the term Km (x),
Z
∂v(z)
Km (x) = − w(z) dσz
Sm (x) ∂ν
Z Z
∂K ∂ψx (z)
= w(z) (z − x) dσz − w(z) dσz .
Sm (x) ∂ν Sm (x) ∂ν
The second term goes to zero by taking the sup-norm outside the integral.
−1
To tackle the first term, we note that ∇z K(z − x) = ωn |z−x| n (z − x). Since
and
∂Rn+ = {x = (x1 , x2 , . . . , xn } ∈ Rn | xn = 0}.
To compute the Green’s function, we shall use the method of reflection.
The reflection technique ensures that the points on the boundary (along
which the reflection is done) remains unchanged to respect the imposed
Dirichlet condition.
−1 y − x⋆ y−x
∇y G(x, y) = ⋆ n
−
ωn |y − x | |y − x|n
−1
∇y G(x, y) = (x − x⋆ ).
ωn |y − x|n
CHAPTER 6. THE LAPLACIAN 164
2xn
∇y G(x, y) · ν = .
ωn |y − x|n
2xn
P (x, y) :=
ωn |y − x|n
It now remains to show that the u as defined above is, indeed, a solution of
(6.8.2) for Rn+ .
Exercise 30. Let f ∈ C(Rn+ ) be given. Let g ∈ C(Rn−1 ) be bounded. Then
u as given in (6.8.6) is in C 2 (Rn+ ) and solves (6.8.2).
Definition 6.8.15. For any x ∈ Rn \ {a}, we define its reflection along the
2 (x−a)
circle Sr (a) as x⋆ = r|x−a|2 + a.
The idea behind reflection is clear for the unit disk, i.e., when a = 0 and
r = 1, as x⋆ = |x|x2 . The above definition is just the shift of origin to a and
dilating the unit disk by r.
CHAPTER 6. THE LAPLACIAN 165
Therefore, |x−a|
r
|y−x⋆ | = |y−x| for all y ∈ Sr (a). For each fixed x ∈ Br (a), we
need to find a harmonic function ψx in Br (a) solving (6.8.4). Since K(· − x)
is harmonic in Br (a) \ {x}, we use the method of reflection to shift the
singularity of K at x to the complement of Br (a). Thus, we define
|x − a| ⋆
ψx (y) = K (y − x ) x 6= a.
r
2−n
For n ≥ 3, K |x−a| r
(y − x ⋆
) = |x−a|
r 2−n
K(y − x⋆ ). Thus, for n ≥ 3, ψx solves
(6.8.4), for x 6= a. For n = 2,
|x − a| ⋆ −1 |x − a|
K (y − x ) = ln + K(y − x⋆ ).
r 2π r
Hence ψx solves (6.8.4) for n = 2. Note that we are yet to identify a harmonic
function ψa corresponding to x = a. We do this by setting ψa to be the
constant function (
1
− 2π ln r (n = 2)
ψa (y) := r 2−n
ωn (n−2)
(n ≥ 3).
Thus, ψa is harmonic and solves (6.8.4) for x = a. Therefore, we define the
Green’s function to be
|x − a| ⋆
G(x, y) := K (y − x ) − K(y − x) ∀x, y ∈ Br (a), x 6= a and x 6= y
r
CHAPTER 6. THE LAPLACIAN 166
and (
1 r
− 2π ln |y−a|
(n = 2)
G(a, y) :=
1 2−n
ωn (n−2)
(r − |y − a|2−n ) (n ≥ 3).
We shall now compute the normal derivative of G. Recall that
−1
∇y K(y − x) = (y − x)
ωn |y − x|n
−|x−a|2−n
and one can compute ∇y K |x−a|
r
(y − x ⋆
) = r2−n ωn |y−x⋆ |n
(y−x⋆ ). Therefore,
−1 |x − a|2−n (y − x⋆ ) y−x
∇y G(x, y) = − .
ωn r2−n |y − x⋆ |n |y − x|n
If y ∈ Sr (a), we have
−1 |x − a|2 ⋆
∇y G(x, y) = (y − x ) − (y − x)
ωn |y − x|n r2
−1 |x − a|2
= − 1 (y − a)
ωn |y − x|n r2
Since the outward unit normal at any point y ∈ Sr (a) is 1r (y − a), we have
X n
−1 |x − a|2 1
∇y G(x, y) · ν = n 2
−1 (yi − ai )2
ωn |y − x| r r
2
i=1
−r |x − a|
= n
−1 .
ωn |y − x| r2
Definition 6.8.16. For all x ∈ Br (a) and y ∈ Sr (a), the map
r2 − |x − a|2
P (x, y) :=
rωn |y − x|n
is called the Poisson kernel for Br (a).
Now substituing for G in (6.8.5), we get the Poisson formula for u,
Z Z
r2 − |x − a|2 g(y)
u(x) = − G(x, y)f (y) dy + n
dσy . (6.8.7)
Br (a) rωn Sr (a) |y − x|
It now remains to show that the u as defined above is, indeed, a solution of
(6.8.2) for Br (a).
Exercise 31. Let f ∈ C(Br (a)) be given. Let g ∈ C(Sr (a)) be bounded.
Then u as given in (6.8.7) is in C 2 (Br (a)) and solves (6.8.2).
CHAPTER 6. THE LAPLACIAN 167
J(u) ≤ J(v) ∀v ∈ V.
J(u) ≤ J(v) ∀v ∈ V
J(u) ≤ J(u + tφ) (for any φ ∈ C 2 (Ω) such that φ = 0 on ∂Ω)
1
0 ≤ (J(u + tφ) − J(u))
t Z Z
1 1 2 2
0 ≤ t |∇φ| + 2t∇φ · ∇u dx − t f φ dx
t 2 Ω Ω
where ∂u
∂ν
:= ∇u · ν and ν = (ν1 , . . . , νn ) is the outward pointing unit normal
vector field of ∂Ω. Thus, the boundary imposed is called the Neumann
boundary condition. The solution of a Neumann problem is not necessarily
unique. If u is any solution of (6.9.1), then u + c for any constant c is also a
solution of (6.9.1). More generally, for any v such that v is constant on the
connected components of Ω, u + v is a solution of (6.9.1).
CHAPTER 6. THE LAPLACIAN 170
Appendices
171
Appendix A
Note that the gamma function is defined as an improper integral and its
existence has to be justified. Observe that for a fixed x > 0,
|e−t tx−1 | = e−t tx−1 ≤ tx−1 ∀t > 0,
R1
since for t > 0, |e−t | ≤ 1. Now, since 0 tx−1 dt exists, we have by comparison
R1
test the existence of the integral 0 e−t tx−1 dt. Now, for t → ∞, e−t tx−1 → 0
and hence the is a constant C > 0 such that
tx−1 e−t ≤ C/t2 ∀t ≥ 1.
R∞
Since 1 1/tR2 dt exists, we again have using comparison test the existence of
∞
the integral 1 e−t tx−1 dt. In fact, the gamma function can be defined for
any complex number z ∈ C such that Re(z) > 0.
It is worth noting that the gamma function Γ generalises the notion of
factorial of positive integers. This would be the first property we shall prove.
Exercise 32. Show that Γ(x + 1) = xΓ(x). In particular, for any positive
√
integer n, Γ(n + 1) = n!. Also, show that Γ(1) = 1 and Γ(1/2) = π.
Further, for any positive integer n,
√
Γ(n + 1/2) = (n − 1/2)(n − 3/2) . . . (1/2) π.
(Hint: Use integration by parts and change of variable).
173
APPENDIX A. THE GAMMA FUNCTION 174
Theorem A.0.1. Let f be positive and continuous on (0, ∞) and let log f
be convex on (0, ∞). Also, let f satisfy the recursive equation
f (x + 1) = xf (x) ∀x > 0
In particular, the above computation is true for the point p0 . Since r′ (t0 )
is the slope of the tangent at t0 to the curve C, we see that the vector ∇S(p0 )
is perpendicular to the tangent vector at p0 . Since this argument is true for
any curve that passes through p0 . We have that ∇S(p0 ) is normal vector to
the tangent plane at p0 . If, in particular, the equation of the surface is given
as S(x, y, z) = u(x, y) − z, for some u : R2 → R, then
175
APPENDIX B. NORMAL VECTOR OF A SURFACE 176
Appendix C
177
APPENDIX C. IMPLICIT FUNCTION THEOREM 178
Appendix D
Divergence Theorem
179
APPENDIX D. DIVERGENCE THEOREM 180
2π n/2
Γ(n/2)
and the volume of the ball B1 (0) in Rn is ωn /n. Consequently, for any x ∈ Rn
and the r > 0, the surface area of Sr (x) is rn−1 ωn and the volume of Br (x)
is rn ωn /n.
181
APPENDIX E. SURFACE AREA AND VOLUME OF DISK IN RN 182
Therefore,
Z Z n
Y
−π|x|2 2
In := e dx = e−πxi dx
Rn R i=1n
YZ
n
2
= e−πt dt
i=1 R
Z n
−πt2
= e dt = (I1 )n
R
Z Z 2(n/2)
−π|x|2 −πt2
n/2
e dx = e dt = (I1 )2 = (I2 )n/2
Rn R
Z n/2
−π|y|2
= e dy
R2
Z 2π Z ∞ n/2
−π|y|2
= e dy
0 0
Z 2π Z ∞ n/2
−πr 2
= e r dr dθ (since jacobian is r)
0 0
Z ∞ n/2
−πr 2
= 2π e r dr
0
Z ∞ n/2
−πs
= π e ds (by setting r2 = s)
0
Z ∞ n/2
−q
= e dq (by setting πs = q)
0
= (Γ(1))n/2 = 1.
Let ωn denote the surface area of the unit sphere S1 (0) in Rn , i.e.,
Z
ωn = dσ,
S1 (0)
APPENDIX E. SURFACE AREA AND VOLUME OF DISK IN RN 183
where t = s/r. Thus, the surface area of Sr (x) is rn−1 ωn . Similarly, volume
of a disk Br (x) is rn ωn /n.
APPENDIX E. SURFACE AREA AND VOLUME OF DISK IN RN 184
Appendix F
is in CRc∞ (Rn ) with supp(ρ) = B1 (0), ball with centre 0 and radius 1, where
−1
c−1 = |x|≤1 exp( 1−|x| 2 ) dx.
Thus, one can introduce a sequence of functions in Cc∞ (Rn ), called mol-
lifiers. For ε > 0, we set
( 2
cε−n exp( ε2−ε
−|x|2
) if |x| < ε
ρε (x) = (F.0.1)
0 if |x| ≥ ε,
185
APPENDIX F. MOLLIFIERS AND CONVOLUTION 186
Similarly, one can show that, for any tuple α, Dα uε (x) = (Dα ρε ∗ u)(x).
Thus, uε ∈ C ∞ (Ωε ).
Appendix G
Duhamel’s Principle
Notice that x0 e−at is a solution of the homogeneous ODE. Thus, the solution
x(t) can be given as
Z t
x(t) = S(t)x0 + S(t − s)f (s) ds
0
where S(t) is a solution operator of the linear equation, given as S(t) = e−at .
Consider the second order inhomogeneous ODE
′′
x (t) + a2 x(t) = f (t) in (0, ∞)
x(0) = x0 (G.0.2)
′
x (0) = x1 .
187
APPENDIX G. DUHAMEL’S PRINCIPLE 188
x′ (t) = ay(t).
Then
f (t)
y ′ (t) = − ax(t)
a
and the second order ODE can be rewritten as a system of first order ODE
with the initial condition X0 := X(0) = (x0 , x1 /a). We introduce the matrix
P (At)n
exponential eAt = ∞ n=1 n! . Then, multiplying the integration factor e
At
Notice that X0 e−At is a solution of the homogeneous ODE. Thus, the solution
X(t) can be given as
Z t
X(t) = S(t)X0 + S(t − s)F (s) ds
0
where S(t) is a solution operator of the linear equation, given as S(t) = e−At .
Bibliography
189
BIBLIOGRAPHY 190
Index
integral curve, 22
integral surface, 23
191