Method of Characteristics for First-Order PDEs

⚡ Key Takeaways

The Method of Characteristics converts a first-order partial differential equation (PDE) into a coupled system of ordinary differential equations (ODEs) along special curves called characteristics, along which the PDE collapses to a simple total-derivative relation. Key result: the solution surface of $a u_x + b u_y = c$ is foliated by curves on which information propagates, reducing the problem to integrating ODEs. Applies to quasilinear and fully nonlinear first-order PDEs in any number of variables. Origin: developed by Cauchy, Monge, and Charpit in the late 18th–19th centuries.

💡 The Core Insight

A first-order PDE prescribes the rate of change of $u$ in one special direction at each point. If you move along a curve tangent to that direction, the PDE tells you exactly how $u$ changes—an ODE. By stitching together these one-dimensional solutions, you reconstruct the entire solution surface. The PDE is really a family of ODEs in disguise.

DefinitionQuasilinear First-Order PDE

A quasilinear first-order PDE in two independent variables has the form

1

\quad a(x,y,u)\,u_x + b(x,y,u)\,u_y = c(x,y,u)

Derivation hint: This is the most general first-order PDE that is linear in the highest (first) derivatives $u_x, u_y$ ; coefficients may depend nonlinearly on $u$ itself. Here $u=u(x,y)$ is the unknown, $u_x = \partial u/\partial x$ , $u_y = \partial u/\partial y$ , and $a,b,c$ are given functions. The solution is a surface $z = u(x,y)$ in $(x,y,z)$ -space.

DefinitionCharacteristic Curves and the Characteristic System

The characteristic curves are the projections onto the $(x,y)$ -plane of the integral curves of the characteristic system of ODEs:

2

\quad \frac{dx}{dt} = a(x,y,u),\qquad \frac{dy}{dt} = b(x,y,u),\qquad \frac{du}{dt} = c(x,y,u)

Derivation hint: Parametrize a curve $(x(t),y(t))$ in the domain; the chain rule gives $du/dt = u_x \dot x + u_y \dot y$ . Choosing $\dot x = a$ , $\dot y = b$ makes the right-hand side equal $a u_x + b u_y = c$ . The full three-dimensional curves $(x(t),y(t),u(t))$ lying on the solution surface are called characteristic strips (in the fully nonlinear case) or simply characteristics.

DefinitionCauchy Problem and Initial Curve

The Cauchy problem supplies initial data on a curve $\Gamma$ in the domain, parametrized by $s$ :

3

\quad x = x_0(s),\quad y = y_0(s),\quad u = u_0(s) \quad\text{at } t=0

Derivation hint: We prescribe the value of $u$ along a one-dimensional curve $\Gamma$ ; characteristics emanate from each point of $\Gamma$ , sweeping out a two-dimensional surface as $(s,t)$ vary. The solution surface is generated by launching one characteristic from each point of $\Gamma$ .

TheoremExistence and Uniqueness via the Transversality Condition

Let $a,b,c \in C^1$ and let the initial curve $\Gamma$ be $C^1$ . If the transversality (non-characteristic) condition

4

\quad J = \det\begin{pmatrix} x_0'(s) & a(x_0,y_0,u_0)\\[2pt] y_0'(s) & b(x_0,y_0,u_0)\end{pmatrix} = x_0'(s)\,b - y_0'(s)\,a \neq 0

Derivation hint: $J$ is the Jacobian $\partial(x,y)/\partial(s,t)$ at $t=0$ ; it is nonzero precisely when $\Gamma$ is nowhere tangent to a characteristic direction, so the map $(s,t)\mapsto(x,y)$ is locally invertible by the inverse function theorem. holds along $\Gamma$ , then there exists a unique $C^1$ solution $u(x,y)$ in a neighborhood of $\Gamma$ satisfying the PDE and the initial data.

⚠️ Common mistake

If $J=0$ on $\Gamma$ —i.e., the initial curve is a characteristic—the Cauchy problem has either no solution or infinitely many. Data must be consistent with the ODE $du/dt=c$ for any solution to exist at all. Always check the transversality condition first.

ProofReduction of the PDE to ODEs

Step 1 — Geometry of the solution surface. Suppose $z=u(x,y)$ solves [1]. Write the surface implicitly as $F(x,y,z)=u(x,y)-z=0$ . Its normal vector is

5

\quad \nabla F = (u_x,\,u_y,\,-1).

Derivation hint: Take gradients of $F = u(x,y)-z$ ; the $z$ -component is $\partial(-z)/\partial z = -1$ .

Step 2 — The PDE as orthogonality. Equation [1] states $a u_x + b u_y - c = 0$ , which is exactly

6

\quad (a,\,b,\,c)\cdot(u_x,\,u_y,\,-1) = 0.

Derivation hint: Expand the dot product; the third term is $c\cdot(-1)=-c$ , and setting the sum to zero recovers [1]. So the vector field $(a,b,c)$ is tangent to the solution surface at every point.

Step 3 — Integral curves lie on the surface. A curve whose tangent is always $(a,b,c)$ —i.e., a solution of system [2]—stays within the surface once it starts there, because its velocity is perpendicular to the surface normal:

7

\quad \frac{d}{dt}F(x(t),y(t),u(t)) = \nabla F\cdot(\dot x,\dot y,\dot u) = (u_x,u_y,-1)\cdot(a,b,c)=0.

Derivation hint: Apply the chain rule to $F$ along the curve, then substitute [2] for the velocity and use [6]. Hence $F$ is constant ( $=0$ ) along characteristics: they remain on the solution surface.

Step 4 — Building the surface. Launch a characteristic from each point $(x_0(s),y_0(s),u_0(s))$ of $\Gamma$ . This yields $x=x(s,t)$ , $y=y(s,t)$ , $u=u(s,t)$ . When $J\neq 0$ (condition [4]), invert to get $s=s(x,y)$ , $t=t(x,y)$ and write $u(x,y)=u(s(x,y),t(x,y))$ . This function is $C^1$ and satisfies [1]. $\blacksquare$

∎

AC

Augustin-Louis Cauchy

1789–1857

Cauchy formalized the existence theory for first-order PDEs and laid the groundwork for what became the Cauchy problem. Building on Gaspard Monge's geometric picture of PDEs as families of surfaces and Paul Charpit's treatment of the fully nonlinear case (the Lagrange–Charpit method), Cauchy's rigorous approach—using the propagation of data along characteristic curves—turned a collection of geometric tricks into a systematic theory that underlies modern hyperbolic PDE analysis.

🧠 Intuition

Imagine dye dropped into a flowing river. The fluid velocity field is $(a,b)$ , and the dye concentration $u$ rides along with the current. A characteristic curve is simply a path traced by a floating particle. The equation $du/dt = c$ tells you how the concentration changes as you drift downstream—if $c=0$ , the dye value is carried unchanged along the flow (pure transport). The solution everywhere is just "look back upstream along your characteristic to where you started, and read off the initial value (adjusted by integrating $c$ )." Information cannot cross characteristics; each one is an isolated channel of communication.

Figure

Characteristics emanating from an initial curve

\Gamma

; the solution at point

P

is determined by tracing back along its characteristic to

\Gamma

.

ExampleLinear Transport Equation

Given: Solve $u_t + 2\,u_x = 0$ with initial condition $u(x,0)=\sin x$ . (Here $y=t$ , $a=1$ in the $t$ -direction, $b=2$ , $c=0$ .)

Step 1 — Characteristic ODEs. Treating $t$ as the parameter:

8

\quad \frac{dt}{d\tau}=1,\quad \frac{dx}{d\tau}=2,\quad \frac{du}{d\tau}=0.

Derivation hint: Match $u_t + 2u_x=0$ to $a u_t + b u_x = c$ with $a=1,b=2,c=0$ , then read off [2].

Step 2 — Integrate. From $dt/d\tau=1$ take $t=\tau$ . Then $x = 2t + x_0$ and $u=u_0$ (constant along characteristics).

Step 3 — Apply data. At $t=0$ , $x=x_0$ and $u_0=\sin x_0$ . Invert: $x_0 = x-2t$ .

Step 4 — Solution.

9

\quad u(x,t) = \sin(x-2t).

Derivation hint: Substitute $x_0 = x-2t$ into $u=u_0=\sin x_0$ ; the profile $\sin$ translates rightward at speed 2.

\boxed{\,u(x,t)=\sin(x-2t)\,}\qquad\text{(a wave moving right at speed } 2\text{ units/time)}

ExampleQuasilinear Equation with a Source

Given: Solve $u_x + u_y = u$ with $u(x,0)=x$ along $\Gamma:\ (x_0(s),y_0(s),u_0(s))=(s,0,s)$ .

Step 1 — Check transversality. $J = x_0' b - y_0' a = (1)(1)-(0)(1)=1\neq 0$ . ✓ Non-characteristic.

Step 2 — Characteristic system.

10

\quad \frac{dx}{dt}=1,\quad \frac{dy}{dt}=1,\quad \frac{du}{dt}=u.

Derivation hint: Read $a=1,b=1,c=u$ directly from $u_x+u_y=u$ and insert into [2].

Step 3 — Integrate with initial data.

11

\quad x = s+t,\qquad y = t,\qquad u = s\,e^{t}.

Derivation hint: $x(0)=s,\,y(0)=0,\,u(0)=s$ ; the equation $du/dt=u$ gives exponential growth $u=u(0)e^t$ .

Step 4 — Invert. $t=y$ and $s = x-t = x-y$ .

Step 5 — Solution.

12

\quad u(x,y) = (x-y)\,e^{\,y}.

Derivation hint: Substitute $s=x-y$ , $t=y$ into $u=se^t$ .

\boxed{\,u(x,y)=(x-y)e^{y}\,}

Check: $u_x = e^y$ , $u_y = -e^y+(x-y)e^y$ ; sum $= (x-y)e^y = u$ . ✓

⚠️ Common mistake

"Characteristics are always straight lines." False. They are straight only when $a$ and $b$ are constant (or, in the quasilinear case, when $u$ stays constant along them so the slope $b/a$ does not change). For genuinely nonlinear equations like the inviscid Burgers' equation $u_t + u u_x = 0$ , characteristic slopes depend on $u$ , so they generically bend, converge, and intersect.

CounterexampleShock Formation in Burgers' Equation

Consider $u_t + u\,u_x = 0$ with $u(x,0)=-x$ on $[-1,1]$ . Characteristics satisfy $dx/dt = u = u_0 = -x_0$ , so

13

\quad x = x_0 - x_0\,t = x_0(1-t),\qquad u = -x_0.

Derivation hint: $u$ is constant along characteristics, so its initial value $-x_0$ persists; the speed equals that constant. At $t=1$ all characteristics collapse to $x=0$ : the map $(x_0,t)\mapsto x$ becomes non-invertible ( $J=1-t=0$ ). The classical solution ceases to exist at $t=1$ even though the initial data was smooth. Beyond this, a weak solution with a shock wave is required, governed by the Rankine–Hugoniot condition. This shows the method's local nature: smooth solutions need not extend globally.

connection

The Method of Characteristics is the gateway to several neighboring areas. It generalizes the Linear Transport Equation to nonlinear settings and underpins the theory of Conservation Laws, where shock formation and weak solutions take center stage. The fully nonlinear version connects to Hamilton-Jacobi Equations via the Lagrange–Charpit equations, which mirror Hamilton's canonical equations in Classical Mechanics. Geometrically it ties to First-Order ODE Systems (the characteristic system itself) and to the Inverse Function Theorem (the transversality condition). In numerical analysis, it informs upwind schemes and the CFL Condition for stability.

🔬 Where this appears

⏳ History

- 1770s — Joseph-Louis Lagrange develops methods for linear and quasilinear first-order PDEs, introducing the idea of auxiliary ODEs (Lagrange's method).

- 1784 — Paul Charpit de Villecourt submits a memoir treating the general (fully nonlinear) first-order PDE; though never published, it founds the Lagrange–Charpit method.

- 1770s–1810s — Gaspard Monge gives the geometric interpretation of PDEs via families of surfaces, characteristic curves, and the "Monge cone."

- 1819 — Pierre-Simon Laplace and contemporaries refine the connection to envelope and contact geometry.

- 1830s — Augustin-Louis Cauchy establishes rigorous existence theory for the Cauchy problem along non-characteristic data.

- 1860s — William Rowan Hamilton's earlier optics work (1820s–30s) is recognized as the characteristic theory of the eikonal/Hamilton–Jacobi equation.

- 1900s — Sophus Lie embeds the theory within contact transformations and Lie group symmetry.

- 1950s — Peter Lax, Olga Oleinik, and others develop weak solutions, entropy conditions, and shock theory for nonlinear conservation laws, extending characteristics past their breakdown.

- 1980s — Crandall & Lions introduce viscosity solutions, giving a rigorous global theory for Hamilton–Jacobi equations where classical characteristics fail.

open-problem

The Method of Characteristics is fundamentally local: it guarantees a smooth solution only in a neighborhood of the initial curve. For nonlinear equations, characteristics may cross in finite time, producing shocks, and the classical solution no longer exists. The modern resolution—weak solutions with entropy/viscosity conditions—restores global existence and uniqueness, but identifying the physically correct weak solution (entropy condition selection) and extending the theory to systems of conservation laws in multiple space dimensions remains deep and partially open. Global well-posedness for multidimensional systems of hyperbolic conservation laws (e.g., compressible Euler with large data) is a major unsolved problem. Likewise, characteristics for fully nonlinear higher-dimensional PDEs intertwine with the still-active theory of viscosity solutions and weak KAM theory.

Method of Characteristics for First-Order PDEs

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also