The First Order Quasi-Linear Pdes: 2.1 Characteristic Curves and Integral Surfaces
The First Order Quasi-Linear Pdes: 2.1 Characteristic Curves and Integral Surfaces
The First Order Quasi-Linear Pdes: 2.1 Characteristic Curves and Integral Surfaces
The first order quasi-linear PDEs have the following general form:
where α ≥ 0. The distinguished feature of this class is that one can solve it through a
system of ODEs at least locally. Such a method is called method of characteristic. We will
first introduce the general theory in two variables in first two sections, then we generalize
it into general cases. We then discuss in details on the semi-linear case in two variables
with some examples. The transport equation will be presented in fourth section with a
direct approach.
Keeping this equation in mind, the vector field (a, b, c) (called characteristic direction)
therefore plays an essential role in solving (2.3). On some interval I ⊂ R, this vector field
determines the following parametric curves (called characteristic curves)
23
24 CHAPTER 2. THE FIRST ORDER QUASI-LINEAR PDES
From the definition of S, we know that, for each point on S, there is a characteristic
curve passing through. Therefore, S is the union of all characteristic curves. Now, it is not
hard to see from argument we made above that we did prove the following result.
2.2. CAUCHY PROBLEM 25
f (φ, ψ) = 0,
Theorem 2.2.1 Suppose that x0 (s), y0 (s) and u0 (s) are continuous differentiable func-
tions of s in a closed interval 0 ≤ s ≤ 1 and that a, b and c are functions of x, y and
u with continuous first order partial derivatives with respect to their arguments in some
domain D of (x, y, z)-space containing the initial curve
Remark 2.2.2 The curve C, on which the Cauchy data is prescribed, is the projection of
the initial curve γ onto the (x, y)-plane. The condition (2.8) is precisely equivalent to the
curve C is not characteristic in the sense of section 1.2.
Roughly speaking, a proof of this theorem can be described as follows: One starts with
a time t = 0 and solve the equation (2.5) with the initial data (2.7) which gives a unique
solution
x = X(s, t), y = Y (s, t), z = Z(s, t), s ∈ [0, 1], t ∈ [0, T ], (2.10)
for some positive T , such that
(X(s, 0), Y (s, 0), Z(s, 0)) = (x0 (s), y0 (s), u0 (s)).
26 CHAPTER 2. THE FIRST ORDER QUASI-LINEAR PDES
Now, with the condition (2.8), one can solve (s, t) in terms of x and y to obtain
s = φ(x, y), t = ψ(x, y).
Substituting these into z = Z(s, t), one has
u = Z(s, t) = Z(φ(x, y), ψ(x, y)) ≡ u(x, y),
the solution of the Cauchy problem. Of course, in general, the solutions exist only in a
neighborhood of the initial curve γ, and such a solution is call a local solution. The rigorous
proof of this theorem involves the implicit function theorem and we shall not pursuit it
here.
which satisfies the initial datum given in (2.13). Now, similar to (2.8), we propose the
condition %
a1 a2 · · · an (s) %%
∂f1 ∂f2
··· ∂fn %
∂s1 ∂s1 ∂s1 %
J = det .. .. ... .. % =' 0, (2.15)
. . . %%
∂f1 ∂f2
· · · ∂s∂fn−1
n %
∂sn−1 ∂sn−1 S
2.4. SOME EXAMPLES 27
which is equivalent to that the initial surface is not characteristic. Under this condition
(2.15), one can solve from (2.14) for s = (s1 , s2 , · · · , sn−1 ) and t in terms of x1 , · · · , xn , and
then upon to a substitution, the solution u = Z(x1 , · · · , xn ) to (2.11) with initial condition
(2.13) is obtained.
where v > 0 is a constant, and u0 is a regular function. First of all, one easily verifies that
the initial curve
γ : t = 0, x = τ, z = u0 (x)
is not characteristic, as the determinant defined in (2.8) is 1.
In view of Example 2.4.1, we solve the following
)
ẋ(s) = v,
(2.18)
ṫ(s) = 1,
where “ ˙ ” indicates the derivative with respect to the parameter s. By imposing the
initial conditions at s = 0, )
x(0) = x0 ,
t(0) = t0 ,
28 CHAPTER 2. THE FIRST ORDER QUASI-LINEAR PDES
Since the initial condition is assigned on t = 0, which implies (x0 , t0 ) is a point of the curve
and t0 = 0. So the characteristic curve may be rewritten, eliminating the parameter s, as
x(t) = vt + x0 . (2.19)
1
This is a line in (x, t)-plane intersecting {t = 0} at x0 with the slope or, equivalently
v
speed v.
Now let φ(t) = u(x(t), t) be a solution of (2.17). Along the characteristic curve (2.19)
we have:
i.e., the solution of (2.17) is constant along the characteristics (2.19). So you can resolve
the ordinary differential equation for φ:
φ(t) = φ(0),
i.e.
u(vt + x0 , t) = u(x(t), t) = u(x(0), 0) = u0 (x0 ).
To determine the value of solution u at a generic point in (x, t)-plane, one simply determines
the point x0 of intersection between the characteristic passing through (x, t) and the axis
{t = 0}, that is, from x = vt + x0 to obtain
x0 = x − vt.
u(x, t) = u0 (x − vt).
S : x = s1 , y = s2 , z = 0, w = φ(s1 , s2 ),
2.4. SOME EXAMPLES 29
we compute
x 2y 1
∂x ∂y ∂z s1 2s2 1
J = det ∂s1 ∂s1 ∂s1 = det 1 0 0 = −1 '= 0.
∂x ∂y ∂z 0 1 0
∂s2 ∂s2 ∂s2
Therefore, we expect the problem (2.20) has a unique solution near S. We solve the initial
value problem
dx dy
= x, = 2y,
dt dt
dz dw
= 1, = 3w,
dt
dt
(x, y, z, w)|t=0 = (s1 , s2 , 0, φ(s1 , s2 )),
to obtain
x = s1 et , y = s2 e2t , z = t, w = φ(s1 , s2 )e3t .
From the first three equations in the above, we have
t = z, s1 = xe−z , s2 = ye−2z .
where c is a constant.
For this problem, we note that the coefficients (x, t, c) is singular at (0, 0, c) for c '= 0,
where the equation does not make sense except for u ≡ 0. The initial curve is
γ : x = s, t = 1, z = f (s),
Therefore,
x
s= , τ = log t
t
30 CHAPTER 2. THE FIRST ORDER QUASI-LINEAR PDES
τ β
γ1 : x(0, τ ) = τ, t(0, τ ) = , z(0, τ ) = τ
α α
for s = 0, and α '= 0, β are arbitrary constants. One easily checks that both u(x, t) = αβ x,
and u(x, t) = βt are solutions for the equation with the new initial data. We found infinitely
many solutions in this case.
In view of example 2.4.1 and 2.4.2, one easily find the characteristic curves as the solution
of
dx
= u(x, t), x(0) = x0 .
dt
Also, we know that the solution u(x, t) is constant along each characteristic curves. There-
fore, the characteristic curves are the family of straight lines
So, the solution of u at each point (x, t) is thus determined by the characteristic line from
certain x0 on {t = 0} passing through (x, t) where u(x, t) = u0 (x0 ).
We remark that this process cannot go far in general. Suppose there are two point
a < b on x-axis such that u0 (a) > u0 (b). We know that the characteristic line issuing from
(a, 0) will overtake the characteristic line from (b, 0) at
b−a
t= .
u0 (a) − u0 (b)
In this case, we could not assign a value for u(x, t) at this intersection point, which corre-
sponds to two different values. Actually, the solution blows up in this case in its derivative.
More precisely, we carry out the following calculations. Let w = ux , which satisfies
wt + uwx = −w2 .
2.5. LINEAR TRANSPORT EQUATION 31
ut + b • Du = 0, x ∈ Rn , t > 0 (2.23)
arises from many important applications, such as particle mechanics, kinetic theory. Here
b = (b1 , b2 , · · · , bn ) is a constant vector. Of course, one can easily solve this PDE using the
characteristic method. Here, we will use a direct approach instead. Actually, the (2.23)
can be rewritten as
(Du, ut ) • (b, 1) = 0,
which means the function u(x, t) is constant in the direction of (b, 1) on the (x, t)-space.
Therefore, if we know the value of u at any point of the straight line through (x, t) with
the direction (b, 1) in (x, t)-space, we know the value of u(x, t).
Now let us consider the initial-value problem
)
ut + b • Du = 0, x ∈ Rn , t > 0,
(2.24)
u(x, 0) = g(x), x ∈ Rn ,
where g(x) is a given function. Fix a point (x, t), the line through (x, t) in the direction
ofd (b, 1) is given in the parametric form with the parameter s as
(x + sb, t + s), s ∈ R.
The line hits the initial plane {t = 0} when s = −t at the point (x − bt, 0). We know that
u is constant on the line, we thus obtain
We see from the above that if (2.24) has a smooth solution, then it is given by (2.25).
Conversely, if g is continuous differentiable, the (2.25) is the unique solution of (2.24). We
also remark that the solution in the form of (2.25) reads as the information traveling in a
constant velocity b, such a solution is thus called traveling wave solution.
32 CHAPTER 2. THE FIRST ORDER QUASI-LINEAR PDES
Remark 2.5.1 When g(x) is not C 1 , problem (2.24) does not admit a smooth solution.
However, even when g is not continuous, the formula (2.25) does provide a reasonable
candidate for a solution. We may informally declare that (2.25) is a weak solution of
(2.24). This makes sense even if g is discontinuous and so is u. This notion is very useful
in nonlinear PDEs.
which satisfies
Therefore, 0 0
u(x, t) − g(x − bt) = z(0) − z(−t) = ż(s) ds
−t
0 0
= f (x + sb, t + s) ds (2.27)
−t
0 t
= f (x + (s − t)b, s) ds,
0
2.6 Problems
Problem 1. Solve the following problems (x, y, z ∈ R, t > 0)
prove that u ≡ 0.
yux − xuy = 0
• (a) x2 ux − y 2 uy = 0, u → ex as y → ∞;
Problem 7. Show that u1 = ex , and u2 = e−y are solutions of the nonlinear equation
(ux + uy )2 − u2 = 0,
but that their sum (ex + e−y ) is not a solution of the equation.
Problem 10. Find the family of curves which represent the general solution of the PDE
Determine the particular member of the family which contains the line u = x and y = 0.