Definition and Basic Form

Standard Equation

Laplace's equation is a second-order linear partial differential equation (PDE) defined as:

∇²φ = 0

Here, ∇² denotes the Laplacian operator, φ = φ(x, y, z) is a twice-differentiable scalar function.

In Cartesian Coordinates

Expanded in 3D Cartesian coordinates:

∂²φ/∂x² + ∂²φ/∂y² + ∂²φ/∂z² = 0

Dimensions

Laplace's equation applies in any dimension n, with the Laplacian defined as sum of second partial derivatives over all spatial variables.

Derivation and Physical Interpretation

Origin from Potential Theory

Describes equilibrium states where potential fields have no local maxima or minima inside domain.

Steady-State Heat Equation

Obtained by setting time derivative in heat equation to zero:

∂u/∂t = α∇²u → steady state: ∇²u = 0

Electrostatics

Electric potential in charge-free region satisfies Laplace's equation, derived from Maxwell's equations.

Mathematical Properties

Linearity

Equation is linear: if φ₁ and φ₂ satisfy, then any linear combination also satisfies.

Elliptic Type

Classified as elliptic PDE, implying smoothness of solutions inside domain.

Maximum Principle

Solution attains maxima and minima on the boundary; no internal extrema unless constant.

Uniqueness

Boundary value problems for Laplace's equation have unique solutions under suitable conditions.

Boundary Conditions

Dirichlet Boundary Condition

Function value φ specified on boundary ∂Ω.

Neumann Boundary Condition

Normal derivative ∂φ/∂n specified on boundary.

Mixed (Robin) Condition

Combination of Dirichlet and Neumann conditions.

Well-Posedness

Correct boundary conditions ensure existence, uniqueness, and stability of solutions.

Analytical Solution Methods

Separation of Variables

Decomposes PDE into ODEs solvable by eigenfunction expansions.

Integral Transform Methods

Fourier and Laplace transforms convert PDE to algebraic equations.

Green's Functions

Express solution in terms of integral over boundary data using fundamental solution.

Conformal Mapping

Transforms complex domains to simpler ones preserving Laplace's equation form.

Harmonic Functions

Definition

Functions φ satisfying Laplace's equation are harmonic: ∇²φ = 0.

Mean Value Property

Value at a point equals average over any surrounding sphere.

Analyticity

Harmonic functions are infinitely differentiable and real-analytic inside domain.

Subharmonic and Superharmonic Functions

Generalizations used in potential theory and PDE analysis.

Applications in Physics and Engineering

Electrostatics

Potential distribution in charge-free regions; design of capacitors, shielding.

Fluid Dynamics

Velocity potential for incompressible, irrotational flow satisfies Laplace's equation.

Heat Conduction

Steady-state temperature distributions in solids without internal heat sources.

Gravitational Fields

Potential fields in vacuum; planetary and astrophysical modeling.

Image Processing

Smoothing and interpolation using harmonic functions.

Numerical Methods

Finite Difference Method (FDM)

Discretizes domain; approximates derivatives by difference quotients.

Finite Element Method (FEM)

Divides domain into elements; uses variational formulation for approximate solutions.

Boundary Element Method (BEM)

Reduces problem dimension by formulating integral equations on boundary.

Iterative Solvers

Gauss-Seidel, Successive Over-Relaxation accelerate convergence for linear systems.

Numerical MethodAdvantagesDisadvantages
Finite Difference MethodSimple implementation, good for regular gridsLess flexible for complex geometries
Finite Element MethodHandles complex domains, flexible meshingMore complex implementation, higher computational cost
Boundary Element MethodReduced dimensionality, efficient for infinite domainsRequires Green's function, less general applicability

Classical Examples

Laplace's Equation in a Rectangle

Solution by separation of variables with Dirichlet boundary conditions on edges.

Potential Around a Cylinder

Analytic solution using cylindrical coordinates and boundary conditions on surface.

Steady-State Temperature in a Circular Disk

Radial symmetry reduces PDE to ODE; Bessel functions may appear in solution.

Example: 2D Laplace Equation in rectangle (0 < x < a, 0 < y < b)∇²φ = ∂²φ/∂x² + ∂²φ/∂y² = 0Boundary conditions: φ(0,y) = 0, φ(a,y) = 0 φ(x,0) = 0, φ(x,b) = f(x)Solution by separation: φ(x,y) = Σ A_n sinh(nπy/a) sin(nπx/a) where A_n determined by Fourier sine series of f(x) 

Laplacian Operator and Generalizations

Definition

Laplacian ∇² = div(grad) = Σ ∂²/∂x_i² over all spatial variables x_i.

Coordinate Systems

Expression varies in cylindrical, spherical, and other curvilinear coordinates.

Generalized Laplacians

Includes Laplace-Beltrami operator on manifolds, fractional Laplacians in nonlocal PDEs.

Extensions and Related Equations

Poisson's Equation

Generalization: ∇²φ = f(x), with source term f(x).

Heat Equation

Time-dependent PDE: ∂u/∂t = α∇²u, Laplace's equation is steady-state form.

Wave Equation

Second-order hyperbolic PDE involving Laplacian in spatial variables.

Helmholtz Equation

Stationary wave equation: ∇²φ + k²φ = 0, reduces to Laplace's if k=0.

Common Problems and Challenges

Ill-Posed Problems

Incorrect boundary data can lead to non-existence or non-uniqueness of solutions.

Singularities

Points where boundary conditions or domain geometry induce solution singularities.

Numerical Stability

Discretization errors can accumulate; requires careful mesh design and solver choice.

References

  • Evans, L.C., Partial Differential Equations, American Mathematical Society, vol. 19, 2010, pp. 279-291.
  • Farlow, S.J., Partial Differential Equations for Scientists and Engineers, Dover Publications, 1993, pp. 45-67.
  • Strauss, W.A., Partial Differential Equations: An Introduction, Wiley, vol. 10, 2007, pp. 120-140.
  • Arfken, G.B., Weber, H.J., Mathematical Methods for Physicists, Academic Press, 7th ed., 2013, pp. 654-670.
  • John, F., Partial Differential Equations, Springer-Verlag, 1982, pp. 112-138.