Partial Differential Equations

⚡ Key Takeaways

A partial differential equation (PDE) is an equation relating an unknown multivariable function to its partial derivatives. Key result: PDEs split into elliptic, parabolic, and hyperbolic types, each governing fundamentally different physical behavior (equilibrium, diffusion, waves). Applies to: heat flow, vibrations, electrostatics, fluid dynamics, quantum mechanics, finance. Foundations laid by: d'Alembert, Euler, Fourier, and Cauchy (18th–19th centuries).

💡 The Core Insight

A PDE is not solved by a single formula but by a function constrained by local relations among its rates of change. The type of the equation — elliptic, parabolic, or hyperbolic — is determined by the second-order terms alone and dictates which auxiliary conditions (boundary or initial data) make the problem well-posed. Geometry of the equation's characteristics tells you how information propagates.

DefinitionPartial Differential Equation

A PDE of order $k$ for an unknown function $u:\Omega\subseteq\mathbb{R}^n\to\mathbb{R}$ is a relation

1

\quad F\!\left(x,\,u,\,Du,\,D^2u,\,\dots,\,D^k u\right)=0,

where $D^j u$ denotes the collection of all partial derivatives of order $j$ , and $k$ is the highest order of derivative appearing.

$x=(x_1,\dots,x_n)\in\Omega$ : independent variables.
$u$ : the unknown (dependent) function.
$Du=(u_{x_1},\dots,u_{x_n})$ : the gradient.
$D^2u$ : the Hessian matrix of second partials $u_{x_ix_j}$ .
$F$ : a given function of all these arguments.

DefinitionLinearity and Classification by Form

A second-order PDE in two variables has the form

2

\quad A\,u_{xx}+2B\,u_{xy}+C\,u_{yy}+D\,u_x+E\,u_y+Fu=G,

where $A,B,C,D,E,F,G$ may depend on $x,y$ . It is:

linear if all coefficients depend only on $x,y$ ;
semilinear if $A,B,C$ depend only on $x,y$ but lower-order terms may depend on $u,u_x,u_y$ ;
quasilinear if $A,B,C$ may depend on $u,Du$ but the equation is linear in the highest derivatives;
fully nonlinear otherwise.

TheoremClassification by Discriminant

The second-order linear PDE [2] is classified at a point by the discriminant $\Delta = B^2 - AC$ :

3

\quad \begin{cases} \Delta < 0 & \text{elliptic (e.g. Laplace } u_{xx}+u_{yy}=0\text{)}\\[2pt] \Delta = 0 & \text{parabolic (e.g. heat } u_t = u_{xx}\text{)}\\[2pt] \Delta > 0 & \text{hyperbolic (e.g. wave } u_{tt}=u_{xx}\text{)} \end{cases}

The type may vary from point to point (e.g. the Tricomi equation $u_{xx}+x\,u_{yy}=0$ ).

ProofReduction to Canonical Form via Characteristics

We show why the discriminant governs the equation's character by seeking characteristic curves $\phi(x,y)=\text{const}$ along which the second-order part degenerates. 1. Seek a change of variables $\xi=\phi(x,y),\ \eta=\psi(x,y)$ that simplifies [2]. Under this change, the coefficient of $u_{\xi\xi}$ becomes

4

\quad \bar A = A\,\phi_x^2 + 2B\,\phi_x\phi_y + C\,\phi_y^2.

2. Setting $\bar A=0$ (so $\xi$ is a characteristic coordinate) and writing $m=-\phi_x/\phi_y=dy/dx$ along the curve, divide by $\phi_y^2$ :

5

\quad A\left(\frac{dy}{dx}\right)^2 - 2B\,\frac{dy}{dx} + C = 0.

3. Solving this quadratic for the slope of characteristics gives

6

\quad \frac{dy}{dx} = \frac{B \pm \sqrt{B^2 - AC}}{A}.

4. The number of real characteristic directions is governed by the radicand $B^2-AC=\Delta$ :

$\Delta>0$ : two real families → hyperbolic; canonical form $u_{\xi\eta}=\dots$ (wave-like).
$\Delta=0$ : one real family → parabolic; canonical form $u_{\eta\eta}=\dots$ (heat-like).
$\Delta<0$ : no real characteristics → elliptic; canonical form $u_{\xi\xi}+u_{\eta\eta}=\dots$ (Laplace-like).

5. Substituting characteristic coordinates eliminates cross terms and reduces each type to its standard model, completing the classification. $\blacksquare$

∎

JF

Joseph Fourier

1768–1830

Fourier, studying heat conduction for his 1822 treatise Théorie analytique de la chaleur, derived the heat equation and introduced the revolutionary idea that any function on an interval could be represented as a sum of sines and cosines. His method of separating variables and superposing trigonometric solutions remains the prototype for solving linear PDEs on bounded domains. Initially controversial — Lagrange doubted the convergence of his series — Fourier's work founded harmonic analysis and reshaped both mathematics and physics.

🧠 Intuition

Think of a PDE as a set of local rules a function must obey at every point, like a giant interlocking constraint network.

The Laplace equation $u_{xx}+u_{yy}=0$ says each value equals the average of its neighbors — like a stretched soap film or a steady temperature with no hot spots. Such functions are infinitely smooth and have no local maxima inside the region.
The heat equation $u_t = u_{xx}$ says the temperature rises where the profile is concave up (a valley fills in) and falls where concave down. It smears information out instantly and irreversibly — you cannot run it backward.
The wave equation $u_{tt}=u_{xx}$ says acceleration equals curvature: a taut string snaps back toward its average shape, producing oscillation. Signals travel at finite speed along characteristics and keep their shape.

The discriminant is just asking: do disturbances travel along definite paths (hyperbolic), spread everywhere at once (elliptic/parabolic), or sit in equilibrium (elliptic)?

Figure

Characteristic curves and the three PDE types in the (x,t) plane

Example1D Heat Equation by Separation of Variables

Given: A rod of length $L=\pi$ m with thermal diffusivity $\alpha=1\ \text{m}^2/\text{s}$ , insulated sides, ends held at $0^\circ$ C: $u(0,t)=u(\pi,t)=0$ . Initial temperature $u(x,0)=5\sin(2x)\ ^\circ$ C. Solve $u_t = u_{xx}$ .

Steps: 1. Separate: $u=X(x)T(t)$ gives $\dfrac{T'}{T}=\dfrac{X''}{X}=-\lambda$ . 2. Spatial problem $X''+\lambda X=0$ , $X(0)=X(\pi)=0$ gives eigenfunctions $X_n=\sin(nx)$ , $\lambda_n=n^2$ . 3. Temporal: $T_n'=-n^2 T_n \Rightarrow T_n=e^{-n^2 t}$ . 4. General solution: $u=\sum_n b_n \sin(nx)\,e^{-n^2 t}$ . 5. Match initial data $5\sin(2x)$ : only the $n=2$ term survives, $b_2=5$ , all others zero.

7

\quad \boxed{u(x,t)=5\,\sin(2x)\,e^{-4t}\ ^\circ\text{C}}

Check: at $t=0$ it equals $5\sin 2x$ ✓; decays with rate $n^2=4\,\text{s}^{-1}$ . After $t=\ln 2/4\approx0.173$ s the amplitude halves.

Exampled'Alembert Solution of the Wave Equation

Given: Infinite string, wave speed $c=2\ \text{m/s}$ , governed by $u_{tt}=4\,u_{xx}$ . Initial shape $u(x,0)=e^{-x^2}$ m, initially at rest $u_t(x,0)=0$ .

Steps: 1. General solution travels as two waves: $u(x,t)=F(x-ct)+G(x+ct)$ . 2. Rest condition $u_t(x,0)=-cF'+cG'=0 \Rightarrow F=G$ (up to constant). 3. Shape condition $F+G=e^{-x^2}\Rightarrow F=G=\tfrac12 e^{-x^2}$ .

8

\quad \boxed{u(x,t)=\tfrac12 e^{-(x-2t)^2}+\tfrac12 e^{-(x+2t)^2}\ \text{m}}

Interpretation: the initial bump splits into two half-height bumps moving left and right at $2$ m/s. At $t=1$ s, the peaks sit at $x=\pm 2$ m, each of height $0.5$ m.

⚠️ Common mistake

"The general solution of a PDE is a formula with constants, like an ODE." False. Because PDEs constrain functions of several variables, their general solutions involve arbitrary functions, not arbitrary constants. For example, the general solution of $u_x=0$ in two variables is $u(x,y)=f(y)$ for any function $f$ . Pinning down a unique solution requires boundary/initial data along curves or surfaces, not just numbers.

⚠️ Common mistake

"More data always makes a problem better-posed." False. Specifying too much (or the wrong kind of) data can make a problem ill-posed. The backward heat equation, or the Cauchy problem for Laplace's equation, can have solutions that fail to depend continuously on the data — tiny perturbations cause arbitrarily large changes.

CounterexampleHadamard's Ill-Posed Cauchy Problem

Consider Laplace's equation $u_{xx}+u_{yy}=0$ with Cauchy data

u(x,0)=0,\qquad u_y(x,0)=\tfrac1n\sin(nx).

The solution is $u(x,y)=\tfrac{1}{n^2}\sin(nx)\sinh(ny)$ . As $n\to\infty$ the data $\to 0$ uniformly, yet for any $y>0$ the solution grows like $\tfrac{1}{n^2}e^{ny}\to\infty$ . Hence the elliptic Cauchy problem is ill-posed — small data changes produce unbounded solution changes. Elliptic equations want boundary data, not Cauchy data.

connection

PDEs sit at a crossroads of analysis, geometry, and physics.

They generalize the Ordinary Differential Equations to many variables; separation of variables turns a PDE into ODEs.
Solving linear PDEs on intervals invokes Fourier Series and Sturm-Liouville Theory, whose eigenfunctions diagonalize the spatial operator.
The Laplace Equation and Harmonic Functions connect to complex analysis: real and imaginary parts of holomorphic functions are harmonic.
Distributional and energy methods rest on Functional Analysis and Sobolev Spaces, the modern setting for weak solutions.
Conservation laws and characteristics tie hyperbolic PDEs to Vector Calculus and the Divergence Theorem.

🔬 Where this appears

⏳ History

- 1746 — Jean d'Alembert derives the 1D wave equation for a vibrating string and gives its traveling-wave solution

F(x-ct)+G(x+ct)

.

- 1750s — Leonhard Euler extends wave studies; debate erupts over which "functions" are admissible as initial data (the "vibrating string controversy").

- 1752 — Euler writes down the continuity and Euler equations for ideal fluid flow.

- 1782 — Pierre-Simon Laplace studies the potential equation

\nabla^2 u=0

in celestial mechanics.

- 1813 — Siméon Poisson generalizes to

\nabla^2 u = -\rho

, modeling gravitational and electrostatic potentials with sources.

- 1822 — Joseph Fourier publishes the heat equation and Fourier-series method in Théorie analytique de la chaleur.

- 1845–1850 — Navier and Stokes complete the Navier–Stokes equations for viscous fluids.

- 1860s — James Clerk Maxwell unifies electricity and magnetism into PDEs predicting electromagnetic waves.

- 1902 — Jacques Hadamard formalizes well-posedness, exposing the ill-posed elliptic Cauchy problem and reframing what "solving" a PDE means.

open-problem

Navier–Stokes existence and smoothness is a Clay Millennium Prize Problem: it is unknown whether smooth, globally-defined solutions to the 3D incompressible Navier–Stokes equations always exist for smooth initial data, or whether they can form singularities ("blow up") in finite time. More broadly, the theory of nonlinear PDEs lacks a unified framework: regularity, uniqueness, and existence must often be established equation by equation. Even classification breaks down — many physically important equations are of mixed type or fully nonlinear, where characteristics, weak solutions, and entropy conditions (as in the theory of shock-forming conservation laws) become essential and delicate.

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also