markdown
The Heat, Wave, and Laplace Equationsmd 5fc5f7e
lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md
Download PDF

The Heat, Wave, and Laplace Equations

1Introduction

This lecture compares three representative partial differential equations according to the data they require, the regions on which information depends, conservation and dissipation, and the role of boundaries. The preceding classification identified the heat equation as parabolic, the wave equation as hyperbolic, and the Laplace equation as elliptic. We now explain how these distinctions affect standard constant-coefficient models.

data/lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md

2Standard Setting

Let κ>0 and c>0. We consider sufficiently regular classical solutions of constant-coefficient equations without forcing, heat sources, or damping. The heat and wave equations are written in one spatial dimension (x,t), whereas the Laplace equation is posed on ΩRn.

ModelTypeEquationTime order
HeatParabolicut=κuxx1
WaveHyperbolicutt=c2uxx2
LaplaceEllipticΔu=0No time variable

3An Equation Alone Does Not Define a Problem

A PDE problem combines an equation, a domain, and data. On the whole line, the heat and wave equations have no spatial boundary. On a bounded interval 0<x<L, both require endpoint conditions such as Dirichlet, Neumann, or Robin conditions.

  • Because the heat equation is first order in time, one prescribes one initial temperature, u(x,0)=f(x).
  • Because the wave equation is second order in time, one prescribes both displacement u(x,0)=f(x) and velocity ut(x,0)=g(x).
  • The Laplace equation has no time evolution and is posed as a boundary-value problem, for example with u|Ω=h.

Classical solutions of bounded-domain evolution problems also require compatibility where the initial surface meets the spatial boundary. A pure Neumann problem for the Laplace equation has a compatibility condition and is unique only up to an additive constant.

The preceding lecture on initial- and boundary-value problems develops these aspects of the problem specification.

data/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md

4Heat: Infinite Propagation and Smoothing

On R, let fCc(R), where Cc(R) denotes the continuous functions that vanish outside a bounded interval. The solution is expressed using the heat kernel

Gt(x)=14πκtexp(-x24κt),t>0,

as

u(x,t)=RGt(x-y)f(y)dy.

For t>0 this is smooth, and it converges uniformly to f as t0. Since Gt(x-y)>0, if f[PARSE ERROR: Undefined("Command(\"ge\")")]0 and f[PARSE ERROR: Undefined("Command(\"not\")")]0, then u(x,t)>0 at every x for every t>0. Thus a disturbance influences the entire space immediately in this model; this property is called infinite propagation speed.

The lecture on Fourier transforms and PDEs derives the heat kernel.

data/lecture/math/partial-differential-equations/fourier-transforms-and-pdes.lecture.n.md

On 0<x<L with homogeneous Dirichlet conditions,

ddt120Lu2dx=-κ0Lux2dx[PARSE ERROR: Undefined("Command(\"le\")")]0,

so the squared L2 norm is nonincreasing. Under homogeneous Neumann conditions,

ddt0Ludx=κ[ux]0L=0,

so the mean temperature is conserved. The same integration by parts shows that the squared L2 norm is nonincreasing under homogeneous Neumann conditions because its boundary term also vanishes. Boundary conditions therefore affect both conserved quantities and long-time behavior.

A mode is a spatial eigenfunction compatible with the boundary conditions whose time amplitude evolves independently under the linear PDE. For homogeneous Dirichlet conditions, sin(nπx/L) is such a spatial mode, and

un(x,t)=e-κ(nπ/L)2tsin(nπx/L).

Direct substitution verifies both the equation and the boundary conditions. Higher-frequency modes decay faster, which explains the smoothing mechanism. The next lecture develops the selection and superposition of these modes by separation of variables.

5Wave: Finite Propagation and Energy

For xR, t[PARSE ERROR: Undefined("Command(\"ge\")")]0, fC2(R), and gC1(R), d'Alembert's formula is

u(x,t)=f(x-ct)+f(x+ct)2+12cx-ctx+ctg(s)ds.

The value at (x,t) depends only on initial data in [x-ct,x+ct]. This domain of dependence expresses finite propagation with speed c. The factorization

t2-c2x2=(t-cx)(t+cx)

and the characteristic coordinates x-ct and x+ct underlie the formula. On a bounded interval, reflected waves must instead be treated according to the boundary conditions.

Define

E(t)=120L(ut2+c2ux2)dx.

Without forcing,

E(t)=c2[utux]0L.

For fixed endpoints u=0, one has ut=0 at the endpoints; for free endpoints ux=0. In either case the boundary flux vanishes and E(t)=E(0). Boundary input, forcing, or damping generally changes the energy.

data/lecture/math/partial-differential-equations/introduction-to-energy-methods.lecture.n.md

6Laplace: Global Boundary Dependence without Time Evolution

A C2 solution of Δu=0 is harmonic and satisfies the mean-value property. Unlike the preceding evolution equations, the Laplace equation describes a static global relation between boundary and interior values.

For the unit disk, let h be continuous and 2π-periodic, and write eiα=(cosα,sinα). With boundary data u(eiφ)=h(φ), the Poisson integral formula is

u(reiθ)=12π02π1-r21-2rcos(θ-φ)+r2h(φ)dφ,0[PARSE ERROR: Undefined("Command(\"le\")")]r<1.

It shows that one interior value depends on data from the entire boundary. This is static global dependence, not finite or infinite propagation in time.

If Ω is bounded and connected and uC2(Ω)C(Ω_) is harmonic, the maximum principle places the maximum and minimum on Ω. It implies uniqueness for harmonic functions with the same Dirichlet data. Existence requires a separate argument, and a pure Neumann problem remains indeterminate up to an additive constant.

data/lecture/math/partial-differential-equations/maximum-principle-basics.lecture.n.md

7Comparison on Common Axes

AxisHeatWaveLaplace
Basic dataOne initial valueTwo initial valuesBoundary values such as Dirichlet data
Bounded domainAlso requires lateral boundary conditionsAlso requires lateral boundary conditionsRequires spatial boundary conditions
Information dependenceReaches all space at every positive timeDomain of dependence with speed cNo temporal causality; globally coupled to the boundary
Integral quantitySquared norm is nonincreasing for homogeneous Dirichlet or Neumann data; mean is conserved for homogeneous Neumann dataEnergy is conserved in a closed systemStatic equilibrium rather than temporal conservation
Pointwise controlMaximum principle on the initial surface and lateral boundaryNo analogous maximum principleMaximum principle controlled by the spatial boundary
Typical long-time behaviorRelaxes to zero, the initial mean, or a stationary distribution according to boundary dataUndamped fixed-end solutions generally continue to oscillateNo time evolution

For the heat equation, the parabolic boundary consists of the initial surface and the spatial boundary, not the upper surface at the final time. A later lecture proves the maximum principles for the heat and Laplace equations.

The long-time row describes typical behavior for sufficiently regular solutions in the one-dimensional, constant-coefficient, source-free setting of this lecture, with homogeneous or fixed boundary data that are independent of time.

8Dimensional Check

If u is temperature, t is time, and x is length, then κ has dimensions length squared per time. In the wave equation, c has dimensions length per time and is the propagation speed. Accordingly, the diffusion length is κt, whereas the wave travel distance is ct.

9Scope of Validity

The formulas and identities above depend on constant coefficients, absence of sources, and sufficient regularity. For nonlinear diffusion, damping, forcing, inhomogeneous media, or weak data, conservation laws, propagation speeds, and maximum principles require separate analysis. This lecture does not complete the solution methods; it identifies the roles developed by separation of variables, Fourier analysis, maximum principles, and energy methods.

10Related Lectures

data/lecture/math/partial-differential-equations/separation-of-variables-and-fourier-series.lecture.n.md data/lecture/math/vector-calculus/laplacian-and-physical-applications.lecture.n.md
raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる