Good
Good
“The first process ... in the effectual study of sciences must be one of sim-
plification and reduction of the results of previous investigations to a form
in which the mind can grasp them.”
James Clerk Maxwell
Ax2 + Bxy + Cy 2 + Dx + Ey + F = 0,
is positive, zero, or negative. If this is true at all points, then the equation
is said to be hyperbolic, parabolic, or elliptic in a domain. In the case of
two independent variables, a transformation can always be found to reduce
the given equation to canonical form in a given domain. However, in the
case of several independent variables, it is not, in general, possible to find
such a transformation.
To transform equation (4.1.2) to a canonical form we make a change of
independent variables. Let the new variables be
Assuming that ξ and η are twice continuously differentiable and that the
Jacobian
ξx ξy
J = ,
(4.1.5)
ηx ηy
ux = uξ ξx + uη ηx , uy = uξ ξy + uη ηy ,
uxx = uξξ ξx + 2uξη ξx ηx + uηη ηx2 + uξ ξxx + uη ηxx ,
2
(4.1.6)
uxy = uξξ ξx ξy + uξη (ξx ηy + ξy ηx ) + uηη ηx ηy + uξ ξxy + uη ηxy ,
uyy = uξξ ξy2 + 2uξη ξy ηy + uηη ηy2 + uξ ξyy + uη ηyy .
where
4.2 Canonical Forms 93
The resulting equation (4.1.7) is in the same form as the original equation
(4.1.2) under the general transformation (4.1.4). The nature of the equation
remains invariant under such a transformation if the Jacobian does not
vanish. This can be seen from the fact that the sign of the discriminant
does not alter under the transformation, that is,
B ∗2 − 4A∗ C ∗ = J 2 B 2 − 4AC , (4.1.9)
which can be easily verified. It should be noted here that the equation can
be of a different type at different points of the domain, but for our purpose
we shall assume that the equation under consideration is of the single type
in a given domain.
The classification of equation (4.1.2) depends on the coefficients A (x, y),
B (x, y), and C (x, y) at a given point (x, y). We shall, therefore, rewrite
equation (4.1.2) as
These two equations are of the same type and hence we may write them in
the form
dζ = ζx dx + ζy dy = 0.
Thus,
dy ζx
=− , (4.2.3)
dx ζy
φ1 (x, y) = c1 , c1 = constant,
φ2 (x, y) = c2 , c2 = constant.
ξ = φ1 (x, y) , η = φ2 (x, y) ,
uξη = H1 , (4.2.7)
4.2 Canonical Forms 95
α = ξ + η, β = ξ − η, (4.2.8)
This form is called the second canonical form of the hyperbolic equation.
(B) Parabolic Type
In this case, we have B 2 − 4AC = 0, and equations (4.2.5) and (4.2.6)
coincide. Thus, there exists one real family of characteristics, and we obtain
only a single integral ξ = constant (or η = constant).
Since B 2 = 4AC and A∗ = 0, we find that
√ √ 2
A∗ = Aξx2 + Bξx ξy + Cξy2 = A ξx + C ξy = 0.
in which the coefficients assume the same form as the coefficients in equation
(4.1.11). With the use of (4.2.13), the equations A∗ = C ∗ = 0 become
2
Aαx + Bαx αy + Cαy2 − Aβx2 + Bβx βy + Cβy2
+i [2Aαx βx + B (αx βy + αy βx ) + 2Cαy βy ] = 0,
Aαx2 + Bαx αy + Cαy2 − Aβx2 + Bβx βy + Cβy2
−i [2Aαx βx + B (αx βy + αy βx ) + 2Cαy βy ] = 0,
or,
(A∗∗ − C ∗∗ ) + iB ∗∗ = 0, (A∗∗ − C ∗∗ ) − iB ∗∗ = 0.
A∗∗ = C ∗∗ and B ∗∗ = 0.
where H5 = (H4 /A∗∗ ). This is called the canonical form of the elliptic
equation.
We close this discussion of canonical forms by adding an important
comment. From mathematical and physical points of view, characteristics
or characteristic coordinates play a very important physical role in hyper-
bolic equations. However, they do not play a particularly physical role in
parabolic and elliptic equations, but their role is somewhat mathematical
4.2 Canonical Forms 97
y 2 uxx − x2 uyy = 0.
Here
A = y2 , B = 0, C = −x2 .
Thus,
1 2 1 2 1 2 1 2
ξ= y − x , η= y + x .
2 2 2 2
From the relations (4.1.6), we have
ux = uξ ξx + uη ηx = −xuξ + xuη ,
uy = uξ ξy + uη ηy = yuξ + yuη ,
uxx = uξξ ξx2 + 2uξη ξx ηx + uηη ηx2 + uξ ξxx + uη ηxx
= x2 uξξ − 2x2 uξη + x2 uηη − uξ + uη .
uyy = uξξ ξy2 + 2uξη ξy ηy + uηη ηy2 + uξ ξyy + uη ηyy
= y 2 uξξ + 2y 2 uξη + y 2 uηη + uξ + uη .
y 2 uηη = 0.
Thus,
uηη = 0 for y = 0.
4.3 Equations with Constant Coefficients 99
uxx + x2 uyy = 0,
2y − ix2 = c1 , 2y + ix2 = c2 .
Thus, if we write
ξ = 2y − ix2 , η = 2y + ix2 ,
and hence,
1 1
α= (ξ + η) = 2y, β= (ξ − η) = −x2 ,
2 2i
we obtain the canonical form
1
uαα + uββ = − uβ .
2β
It should be remarked here that a given partial differential equation may
be of a different type in a different domain. Thus, for example, Tricomi’s
equation
is elliptic for x > 0 and hyperbolic for x < 0, since B 2 − 4AC = −4x. For
a detailed treatment, see Hellwig (1964).
ξ = y − λ1 x, η = y − λ2 x, (4.3.3)
where
√
B+ B 2 − 4AC
λ1,2 = . (4.3.4)
2A
The linear second-order partial differential equation with constant coeffi-
cients may be written in the general form as
where D1 , E1 , and F1 are constants. Here, since the coefficients are con-
stants, the lower order terms are expressed explicitly.
When A = 0, equation (4.3.1) does not hold. In this case, the charac-
teristic equation may be put in the form
2
−B (dx/dy) + C (dx/dy) = 0,
Integration gives
x = c1 , x = (B/C) y + c2 ,
ξ = x, η = x − (B/C) y. (4.3.8)
uξη = 0. (4.3.10)
where φ and ψ are arbitrary functions, and λ1 and λ2 are given by (4.3.3).
(B) Parabolic Type
When B 2 − 4AC = 0, the equation is of parabolic type, in which case
only one real family of characteristics exists. From equation (4.3.4), we find
that
λ1 = λ2 = (B/2A) ,
y = (B/2A) x + c1 ,
B 2 − 4AC = 0,
uηη = 0. (4.3.14)
102 4 Classification of Second-Order Linear Equations
y = x + c1 , y = (x/4) + c2 .
ξ = y − x, η = y − (x/4) ,
α = ξ + η, β = ξ − η,
in the form
1 1 8
uαα − uββ = uα − uβ − .
3 3 9
Example 4.3.2. The equation
ξ = y + 2x, η = y,
dx2 − c2 dt2 = 0.
Therefore,
x + ct = ξ = constant, x − ct = η = constant.
Thus, the characteristics are straight lines, which are shown in Figure 4.3.1.
The characteristics form a natural set of coordinates for the hyperbolic
equation.
In terms of new coordinates ξ and η defined above, we obtain
−4c2 uξη = 0.
Since c = 0, we have
uξη = 0.
uη = ψ1 (η) .
y + tanh x = constant.
Hence,
ξ = y + tanh x, η = y − tanh x.
B 2 − 4AC = + i sech2 x.
Integrating gives
y + i tanh x = constant.
Thus,
ξ = y + i tanh x, η = y − i tanh x.
1 1
α= (ξ + η) = y, β= (ξ − η) = tanh x.
2 2i
In terms of these new variables, equation (4.3.24b) can be transformed into
the canonical form
2β
uαα + uββ = uβ , |β| < 1. (4.3.26)
1 − β2
Example 4.3.6. Consider the equation
uxx + (2 cosecy) uxy + cosec2 y uyy = 0. (4.3.27)
ξ = x + cos y and η = y.
uηη = 0, for y = 0.
where f (ξ) and g (ξ) are arbitrary functions. In terms of the independent
variables x and y, we have
y y
u (x, y) = y f +g .
x x
108 4 Classification of Second-Order Linear Equations
1 8
vξ = v− .
3 9
This can be easily integrated by separating the variables. Integrating with
respect to ξ, we have
8 1 (ξ/3)
v= + e F (η) .
3 3
Integrating with respect to η, we obtain
8 1
u (ξ, η) = η + g (η) eξ/3 + f (ξ) ,
3 3
where f (ξ) and g (η) are arbitrary functions. The general solution of the
given equation becomes
8 1 1 x 1 (y−x)
u (x, y) = y− + g y− e3 + f (y − x) .
3 4 3 4
Hence, we obtain
uξη = 0.
Integration yields
y = c1 and y = 3x + c2 .
ξ uξη + uη = 0.
1
u (x, y) = f (y − 3x) + g (y) . (4.4.4)
y
Equation (4.4.2) has coefficients A = 1, B = 2, C = 1, D = E = F =
G = 0. Hence, B 2 − 4AC = 0, the equation is parabolic. The characteristic
equation is
dy
= 1,
dx
and the characteristics are
ξ = y − x = c1 and η = y.
uηη = 0.
u (x, y) = y f (y − x) + g (y − x) . (4.4.5)
y = (1 − 2i) x + c1 , y = (1 + 2i) x + c2 ,
and hence,
ξ = y − (1 − 2i) x, η = y − (1 + 2i) x,
f ′′ (ξ) + 2 ξf ′ (ξ) = 0.
which give
u0 1
A= √ and B = u0 .
π 2
Thus, the final solution is
√x
1 4κt 2 1
u (x, t) = u0 √ e−α dα + .
π 0 2
4.6 Exercises
3
(xi) 3 uxx + 4 uxy − 4 uyy = 0.
114 4 Classification of Second-Order Linear Equations
(c) uxx + 5uxy + 4uyy + 7uy = sin x, (d) uxx + uyy + 2ux + 8uy + u = 0,
uxx + xuyy = 0,
4.6 Exercises 115
1
(ii) uαα + uββ + 3β = 0, x > 0.
Show that the characteristic curves for x < 0 are cubic parabolas.
8. Use the polar coordinates r and θ (x = r cos θ, y = r sin θ) to transform
the Laplace equation uxx + uyy = 0 into the polar form
1 1
∇2 u = urr + ur + 2 uθθ = 0.
r r
9. (a) Using the cylindrical polar coordinates x = r cos θ, y = r sin θ, z = z,
transform the three-dimensional Laplace equation uxx + uyy + uzz = 0
into the form
1 1
urr + ur + 2 uθθ + uzz = 0.
r r
(b) Use the spherical polar coordinates (r, θ, φ) so that x = r sin φ cos θ,
y = r sin φ sin θ, z = r cos φ to transform the three-dimensional Laplace
equation uxx + uyy + uzz = 0 into the form
2 1 1
urr + ur + 2 (sin φ uφ )φ + 2 2 uθθ = 0.
r r sin φ r sin φ
ut = κ (uxx + uyy ) ,
uxx + uyy = 0,
u (x, 0) = f (x) and uy (x, 0) = g (x) .
12. Classify each of the following equations and reduce it to canonical form:
(a) y uxx − x uyy = 0, x > 0, y > 0; (b) uxx + sech4 x uyy = 0,
(c) y 2 uxx + x2 uyy = 0, (d) uxx − sech4 x uyy = 0,
into the canonical form. Use the canonical form to find the general
solution.
14. Classify each of the following equations for u (x, t):
where p (x), c (x, t), a (x, t), and α (x) are given functions that take
only positive values in the (x, t) plane. Find the general solution of the
equation in (d).