ODE: Laplace's Equation

Laplace's Equation describes the steady-state distribution of a potential, such as temperature, electric potential, or fluid velocity potential, in a region. Mastering it is essential for engineers because it models equilibrium phenomena where time-dependent transients have decayed, leaving a stable pattern governed entirely by the shape of the domain and its boundaries. Solving it unlocks the ability to design electrostatic shields, predict stress points, and model ideal fluid flow around objects.

The Foundation: Defining Laplace's Equation and Harmonic Functions

In two dimensions, Laplace's equation is a second-order partial differential equation (PDE) given by: $${\nabla }^{2}u=\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}=0$$ Here, $u(x,y)$ represents the steady-state potential function. A function that satisfies Laplace's equation in a domain is called a harmonic function. The equation's physical meaning is profound: at any interior point, the value of $u$ is the average of its values in the immediate neighborhood. There are no local maxima or minima inside the domain; all extremes must occur on the boundary. This is the essence of a steady-state: no sources, sinks, or internal energy generation.

The core mathematical challenge is a boundary value problem. We don't have initial conditions as with wave or heat equations; instead, we prescribe values or derivatives of $u$ along the entire boundary of the domain. The most common type is the Dirichlet problem, where the solution $u$ is specified on the boundary. The Neumann problem specifies the normal derivative (the flux) on the boundary. Solving these problems means finding the unique harmonic function inside that matches these boundary conditions.

The Primary Solution Technique: Separation of Variables

The most powerful analytical method for solving Laplace's equation on simple domains is separation of variables. The technique assumes the solution can be written as a product of functions, each depending on only one coordinate: $u(x,y)=X(x)Y(y)$. Substituting into $\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}=0$ and dividing by $XY$ yields: $$\frac{X^{\prime \prime }(x)}{X(x)}=-\frac{Y^{\prime \prime }(y)}{Y(y)}$$ Since the left side depends only on $x$ and the right only on $y$, both must equal a constant, called the separation constant. This leads to two ordinary differential equations. The choice of this constant (positive, negative, or zero) and the subsequent solutions are dictated by the boundary conditions to avoid trivial solutions.

The geometry of the domain dictates the coordinate system. In rectangular domains, we use Cartesian coordinates $(x,y)$. The separation of variables typically leads to sinusoidal and exponential/hyperbolic function families. For example, in a rectangle where $u=0$ on three sides and a specified function on the fourth, the solution becomes a Fourier sine series in one direction, with coefficients modified by hyperbolic sines in the other direction to satisfy the zero conditions.

In circular domains (disks, annuli, or exterior regions), polar coordinates $(r,\theta )$ are natural. Here, Laplace's equation transforms to: $$\frac{{\partial }^{2}u}{\partial r^2}+\frac{1}{r}\frac{\partial u}{\partial r}+\frac{1}{r^2}\frac{{\partial }^{2}u}{\partial {\theta }^2}=0$$ Applying separation of variables, $u(r,\theta )=R(r)\Theta (\theta )$, leads to a crucial constraint: $\Theta (\theta )$ must be $2\pi$-periodic. This forces the separation constant to be $n^2$ for integer $n$. The radial equation becomes a Cauchy-Euler (equidimensional) equation, with solutions $r^n$ and $r^{-n}$ (or $\ln r$ for $n=0$). The full solution is a Fourier series in $\theta$ with radial coefficients: $$u(r,\theta )=A_0+B_0\ln r+\sum _{n=1}^{\infty }\left[(A_n{r}^n+B_n{r}^{-n})\cos (n\theta )+(C_n{r}^n+D_n{r}^{-n})\sin (n\theta )\right]$$ The constants $A_n,B_n,C_n,D_n$ are determined by the boundary conditions. For a problem inside a full disk of radius $a$, finiteness at $r=0$ requires $B_n=0$ and $B_0=0$.

Key Properties: The Maximum Principle and Uniqueness

Harmonic functions possess elegant theoretical properties that guarantee and characterize solutions. The maximum principle (and its companion, the minimum principle) states that a non-constant harmonic function attains its maximum and minimum values only on the boundary of the domain. This has immediate practical consequences: if you know the boundary temperatures of a steady-state plate, you instantly know that no interior point can be hotter than the hottest boundary point or colder than the coldest.

A direct corollary is the uniqueness of solutions to the Dirichlet problem. If two harmonic functions $u_1$ and $u_2$ satisfy the same Dirichlet boundary conditions, then their difference $v=u_1-u_2$ is harmonic and zero on the boundary. By the maximum and minimum principles, $v$ must be zero everywhere inside the domain, so $u_1=u_2$. This tells engineers that if they find any solution matching the boundary data, it is the solution. Uniqueness for Neumann problems holds only up to an additive constant.

Engineering Applications: Electrostatics and Fluid Flow

The power of Laplace's equation lies in its universal appearance across physics. In electrostatics, in a charge-free region, the electric potential $V$ satisfies ${\nabla }^{2}V=0$. Solving this with boundary conditions (e.g., voltage on conductors) allows you to map the electric field $\vec{E}=-\nabla V$. For instance, finding the potential between two coaxial cylinders at different potentials is a direct application in polar coordinates, crucial for designing capacitors.

In ideal, incompressible, irrotational fluid flow, the velocity potential $\varphi$ (where $\vec{v}=\nabla \varphi$) also satisfies Laplace's equation. Streamlines are contours of the harmonic conjugate function $\psi$. Solving boundary value problems lets you model flow around obstacles. The classic example is flow around a cylinder: the solution in polar coordinates combines a uniform flow term ($A_{1}r\cos \theta$) and a dipole term ($B_1{r}^{-1}\cos \theta$) to satisfy the boundary condition of no flow into the cylinder's surface. This model is foundational for aerodynamics and hydrodynamics.

Common Pitfalls

1. Misapplying boundary conditions to determine the separation constant. A common error is to incorrectly assume the sign of the separation constant $k$ before applying boundary conditions. Always let the physical constraints of the problem (e.g., boundedness, periodicity) dictate whether $k$ is positive, negative, or zero. In a rectangular domain with two homogeneous Dirichlet conditions on opposite sides, a positive $k^2$ leads to a sinusoidal $X(x)$, which can satisfy those conditions for specific eigenvalues.
2. Ignoring finiteness or boundedness at singular points. In polar coordinates, the origin $r=0$ is a singular point of the equation. For a problem defined on a full disk, you must reject solutions like $r^{-n}$ or $\ln r$ because they blow up at the origin, even if they are mathematically valid. Similarly, for an exterior problem ($r>a$), you often require the solution to remain bounded as $r\to \infty$.
3. Confusing the roles of different function families. In rectangular coordinates, the solution along an axis with homogeneous boundary conditions typically yields eigenfunctions (sines/cosines). The solution in the perpendicular direction involves the corresponding hyperbolic sines/cosines or exponentials, which are not periodic. Do not mistakenly try to impose periodicity in this direction.
4. Overlooking the constant term and the logarithmic term in polar solutions. The general solution in polar coordinates includes the terms $A_0$ and $B_0\ln r$. The $B_0\ln r$ term is often discarded for interior disk problems, but it is vital for problems with net flux (like a line charge in electrostatics) or in annular regions. Always include it in your general setup before applying boundary conditions to see if it must be zero.

Summary

• Laplace's equation ${\nabla }^{2}u=0$ models steady-state potential distributions. Its solutions are harmonic functions, which have the mean-value property and obey the maximum principle.
• The primary analytical solution method is separation of variables. The coordinate system must match the domain: Cartesian for rectangular domains and polar for circular domains, leading to solutions built from trigonometric, exponential, hyperbolic, and power functions.
• Solving Laplace's equation is a boundary value problem (Dirichlet or Neumann). The maximum principle ensures uniqueness for the Dirichlet problem, a critical guarantee for engineering applications.
• Key applications include calculating the electric potential in charge-free electrostatics and the velocity potential for ideal, irrotational fluid flow.
• Success requires careful attention to boundary conditions to select the correct separation constant and to reject mathematically valid but physically nonsensical solutions that are unbounded at key points.