Exam focus: Write the characteristic ODEs, carry the initial data along them, and check whether characteristics cross — that crossing is what makes the solution fail to exist globally.
Lec-22 Introduction to First Order PDE
NPTEL · Numerical Methods of ODE and PDE
Lec-37 Method of characteristics for Hyperbolic PDEs - I
NPTEL · Numerical Methods of ODE and PDE
Characteristics, and what happens when they cross.
Part B2023 Juneshow ▾hide ▴
Let u(x, y) be the solution of the Cauchy problem u⋅ux+uy=0 for x∈R,y>0, with u(x, 0) = x for x∈R. Which of the following is the value of u(2, 3)?
A.2
B.3
C.1/2✓
D.1/3
Solution
Characteristics: dx/dy = u with u constant along them, so x=x0+u⋅y and u=x0. Hence u = x/(1 + y), and u(2, 3) = 1/2.
Let u:R2→R solve ∂xu+2∂yu=0 on R2 with u(x, y) = sin x on the line y = 3x + 1, and let v:R2→R solve ∂xv+2∂yv=0 with v(x, 0) = sin x. Let S = [0,1] × [0,1]. Which of the following statements are true?
A.u changes sign in the interior of S.
B.u(x, y) = v(x, y) along a line in S.✓
C.v changes sign in the interior of S.✓
D.v vanishes along a line in S.✓
Solution
Both are constant on lines y − 2x = c. v = sin(x − y/2), which vanishes on x = y/2 and changes sign across it inside S. For u, the data line y = 3x + 1 meets the characteristic through (x, y) at x0=y−2x−1, so u = sin(y − 2x − 1); on S the argument lies in [−3,0]⊂(−π,0], so u ≤ 0 and does not change sign. u = v where y − 2x − 1 = x − y/2, i.e. on the line y = 2x + 2/3, which crosses S.
A.u(x, y) = F(x − y) for an arbitrary differentiable F✓
B.u(x, y) = F(x + y) for an arbitrary differentiable F
C.u(x, y) = F(xy) for an arbitrary differentiable F
D.u(x, y) = c, constants only
Solution
The characteristics of dx/1 = dy/1 are the lines x − y = const, and u is constant along them, so u = F(x − y). Check: u_x + u_y = F′ − F′ = 0. F(x + y) instead solves u_x − u_y = 0.
Consider the PDE x·u_x + y·u_y = 2u on x > 0. Which of the following are true?
A.The projected characteristic curves are the rays y = cx.✓
B.u(x,y)=x2 is a solution.✓
C.u(x, y) = xy is a solution.✓
D.The general solution is u=x2⋅F(y/x) for an arbitrary differentiable F.✓
Solution
Lagrange's auxiliary system dx/x = dy/y = du/2u gives first integrals y/x and u/x2, so u=x2F(y/x). Checks: u=x2 is F≡1(x⋅2x=2x2);u= xy is F(t) = t (xy + yx = 2xy). By Euler's relation these are exactly the functions positively homogeneous of degree 2 — all four statements hold.
Exam focus: B² − 4AC decides the type pointwise; the sign can change with the point (Tricomi). Reduce to canonical form by solving the characteristic ODE dy/dx = (B ± √(B²−4AC))/2A.
Lec-23 Introduction to Second Order PDE
NPTEL · Numerical Methods of ODE and PDE
The B² − 4AC classification; remember it is pointwise.
Part B2023 Decembershow ▾hide ▴
The partial differential equation x2uxx−2xy uxᵧ −3y2uyy+ux−uy=0 is
Consider the second-order PDE u_xx + 4u_xy + 4u_yy = 0. Which of the following are true?
A.The equation is parabolic everywhere.✓
B.The equation is hyperbolic everywhere.
C.The characteristic curves satisfy dy/dx = 2.✓
D.The equation has two distinct families of real characteristics.
Solution
With A = 1, B = 4, C = 4 the discriminant is B2−4AC = 16 − 16 = 0, so the equation is parabolic everywhere — (1) holds, (2) fails. The characteristics satisfy A(dy/dx)2−B(dy/dx) + C = 0, i.e. (dy/dx)2−4(dy/dx) + 4 = 0, giving the repeated root dy/dx = 2 and confirming (3). Because the root is repeated there is only one family of real characteristics, so (4) fails — that is precisely what distinguishes the parabolic case from the hyperbolic one.
Consider the PDE y·u_xx + u_yy = 0 (Tricomi-type). Which of the following are true?
A.The equation is elliptic in the half-plane y > 0.✓
B.The equation is hyperbolic in the half-plane y < 0.✓
C.The equation is parabolic exactly on the line y = 0.✓
D.The type of the equation is the same at every point of the plane.
Solution
With A = y, B = 0, C = 1 the discriminant is B2−4AC = −4y: negative (elliptic) for y > 0, positive (hyperbolic) for y < 0, zero on y = 0. The failure of (4) is the whole point of Tricomi-type equations — the type changes across a line, which is what makes transonic-flow models hard.
3. Laplace, heat and wave equations: separation of variables
Exam focus: Know d'Alembert cold, the separation-of-variables series for the three classical equations, and the qualitative differences: smoothing (heat), finite speed (wave), mean value (Laplace).
Lec-24 Finite Difference Approximations to Parabolic PDEs
NPTEL · Numerical Methods of ODE and PDE
The heat equation and its smoothing behaviour.
Lec-31 Finite Difference Approximations to Elliptic PDEs - I
NPTEL · Numerical Methods of ODE and PDE
Laplace's equation and the maximum principle.
Lec-35 Finite Difference Approximations to Hyperbolic PDEs - I
B.u(x0,t0) depends only on the initial data in [x0− ct0,x0+ ct0].✓
C.u(x, t) = 0 whenever |x| > 1 + ct.✓
D.Initial disturbances influence all of space instantaneously.
Solution
(1) is d'Alembert's formula and (2) its immediate consequence — the domain of dependence. (3) is finite propagation speed: outside the cone reached from [−1, 1] at speed c, both terms vanish. (4) belongs to the heat equation, whose kernel is positive everywhere for t > 0 — the standard contrast with the wave equation.
Let u solve the heat equation u_t = u_xx on 0 < x < L, 0 < t ≤ T, with continuous data. By the maximum principle, the maximum of u over the closed rectangle is attained
A.on the parabolic boundary: the bottom edge t = 0 or the sides x = 0, x = L✓
B.only at interior points
C.only on the top edge t = T
D.at the centre of the rectangle
Solution
The parabolic maximum principle: heat flows downhill, so new interior maxima cannot appear after time 0 — extremes are inherited from the initial data or pushed in through the lateral boundary. The top edge t = T is excluded from the 'boundary' precisely because it is determined by the past.