Laplace's Equation
Laplace's Equation
Laplace's equation
In mathematics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace
who first studied its properties. This is often written as:
where ∆ = ∇2 is the Laplace operator and φ is a scalar function. In general, ∆ = ∇2 is the Laplace–Beltrami or
Laplace–de Rham operator.
Laplace's equation and Poisson's equation are the simplest examples of elliptic partial differential equations. The
general theory of solutions to Laplace's equation is known as potential theory. The solutions of Laplace's equation
are the harmonic functions, which are important in many fields of science, notably the fields of electromagnetism,
astronomy, and fluid dynamics, because they can be used to accurately describe the behavior of electric,
gravitational, and fluid potentials. In the study of heat conduction, the Laplace equation is the steady-state heat
equation.
Definition
In three dimensions, the problem is to find twice-differentiable real-valued functions f, of real variables x, y, and z,
such that
In Cartesian coordinates
In cylindrical coordinates,
In spherical coordinates,
In Curvilinear coordinates,
or
where ∇ • is the divergence operator (also symbolized "div") which maps vectors to scalars, and ∇ is the gradient
operator (also symbolized "grad") which maps vectors to vectors. (hence, the Laplacian Δf ≝ div grad f, maps the
vector function f to a scalar magnitude; specifically it maps the grad (the partial derivatives ) of f to a scalar
(function).)
Laplace's equation 2
If the right-hand side is specified as a given function, h(x, y, z), i.e., if the whole equation is written as
Boundary conditions
The Dirichlet problem for Laplace's
equation consists of finding a solution
φ on some domain D such that φ on
the boundary of D is equal to some
given function. Since the Laplace
operator appears in the heat equation,
one physical interpretation of this
problem is as follows: fix the
temperature on the boundary of the
domain according to the given
specification of the boundary
condition. Allow heat to flow until a
stationary state is reached in which the
Laplace's Equation on an annulus (inner radius r=2 and outer radius R=4) with Dirichlet
temperature at each point on the
Boundary Conditions: u(r=2)=0 and u(R=4)=4sin(5*θ)
domain doesn't change anymore. The
temperature distribution in the interior
will then be given by the solution to the corresponding Dirichlet problem.
The Neumann boundary conditions for Laplace's equation specify not the function φ itself on the boundary of D, but
its normal derivative. Physically, this corresponds to the construction of a potential for a vector field whose effect is
known at the boundary of D alone.
Solutions of Laplace's equation are called harmonic functions; they are all analytic within the domain where the
equation is satisfied. If any two functions are solutions to Laplace's equation (or any linear homogeneous differential
equation), their sum (or any linear combination) is also a solution. This property, called the principle of
superposition, is very useful, e.g., solutions to complex problems can be constructed by summing simple solutions.
Analytic functions
The real and imaginary parts of a complex analytic function both satisfy the Laplace equation. That is, if z = x + iy,
and if
then the necessary condition that f(z) be analytic is that the Cauchy-Riemann equations be satisfied:
Laplace's equation 3
Therefore u satisfies the Laplace equation. A similar calculation shows that v also satisfies the Laplace equation.
Conversely, given a harmonic function, it is the real part of an analytic function, f(z) (at least locally). If a trial form
is
The Laplace equation for φ implies that the integrability condition for ψ is satisfied:
and thus ψ may be defined by a line integral. The integrability condition and Stokes' theorem implies that the value
of the line integral connecting two points is independent of the path. The resulting pair of solutions of the Laplace
equation are called conjugate harmonic functions. This construction is only valid locally, or provided that the path
does not loop around a singularity. For example, if r and θ are polar coordinates and
However, the angle θ is single-valued only in a region that does not enclose the origin.
The close connection between the Laplace equation and analytic functions implies that any solution of the Laplace
equation has derivatives of all orders, and can be expanded in a power series, at least inside a circle that does not
enclose a singularity. This is in sharp contrast to solutions of the wave equation, which generally have less regularity.
There is an intimate connection between power series and Fourier series. If we expand a function f in a power series
inside a circle of radius R, this means that
with suitably defined coefficients whose real and imaginary parts are given by
Therefore
which is a Fourier series for f. These trigonometric functions can themselves be expanded, using multiple angle
formulae.
Laplace's equation 4
Fluid flow
Let the quantities u and v be the horizontal and vertical components of the velocity field of a steady incompressible,
irrotational flow in two dimensions. The condition that the flow be incompressible is that
then the incompressibility condition is the integrability condition for this differential: the resulting function is called
the stream function because it is constant along flow lines. The first derivatives of ψ are given by
and the irrotationality condition implies that ψ satisfies the Laplace equation. The harmonic function φ that is
conjugate to ψ is called the velocity potential. The Cauchy-Riemann equations imply that
Thus every analytic function corresponds to a steady incompressible, irrotational fluid flow in the plane. The real
part is the velocity potential, and the imaginary part is the stream function.
Electrostatics
According to Maxwell's equations, an electric field (u,v) in two space dimensions that is independent of time satisfies
and: where ρ is the charge density. The first Maxwell equation is the integrability condition for
the differential
which is the Poisson equation.It is important to note that the Laplace equation can be used in three-dimensional
problems in electrostatics and fluid flow just as in two dimensions.
Fundamental solution
A fundamental solution of Laplace's equation satisfies
where the Dirac delta function δ denotes a unit source concentrated at the point (x′, y′, z′). No function has this
property, but it can be thought of as a limit of functions whose integrals over space are unity, and whose support (the
region where the function is non-zero) shrinks to a point (see weak solution). It is common to take a different sign
convention for this equation than one typically does when defining fundamental solutions. This choice of sign is
often convenient to work with because −Δ is a positive operator. The definition of the fundamental solution thus
implies that, if the Laplacian of u is integrated over any volume that encloses the source point, then
Laplace's equation 5
The Laplace equation is unchanged under a rotation of coordinates, and hence we can expect that a fundamental
solution may be obtained among solutions that only depend upon the distance r from the source point. If we choose
the volume to be a ball of radius a around the source point, then Gauss' divergence theorem implies that
It follows that
on a sphere of radius r that is centered around the source point, and hence
Note that, with the opposite sign convention (used in Physics), this is the potential generated by a point particle, for
an inverse-square law force, arising in the solution of Poisson equation. A similar argument shows that in two
dimensions
where log(r) denotes the natural logarithm. Note that, with the opposite sign convention, this is the potential
generated by a pointlike sink (see point particle), which is the solution of the Euler equations in two-dimensional
incompressible flow.
Green's function
A Green's function is a fundamental solution that also satisfies a suitable condition on the boundary S of a volume V.
For instance,
may satisfy
and u assumes the boundary values g on S, then we may apply Green's identity, (a consequence of the divergence
theorem) which states that
The notations un and Gn denote normal derivatives on S. In view of the conditions satisfied by u and G, this result
simplifies to
Thus the Green's function describes the influence at (x′, y′, z′) of the data f and g. For the case of the interior of a
sphere of radius a, the Green's function may be obtained by means of a reflection (Sommerfeld, 1949): the source
point P at distance ρ from the center of the sphere is reflected along its radial line to a point P' that is at a distance
Laplace's equation 6
Note that if P is inside the sphere, then P' will be outside the sphere. The Green's function is then given by
where R denotes the distance to the source point P and R' denotes the distance to the reflected point P. A
consequence of this expression for the Green's function is the 'Poisson integral formula. Let ρ, θ, and φ be
spherical coordinates for the source point P. Here θ denotes the angle with the vertical axis, which is contrary to the
usual American mathematical notation, but agrees with standard European and physical practice. Then the solution
of the Laplace equation inside the sphere is given by
where
A simple consequence of this formula is that if u is a harmonic function, then the value of u at the center of the
sphere is the mean value of its values on the sphere. This mean value property immediately implies that a
non-constant harmonic function cannot assume its maximum value at an interior point.
Electrostatics
In free space the Laplace equation of any electrostatic potential must equal zero since ρ (charge density) is zero in
free space.
Taking the gradient of the electric potential we get the electrostatic field
Taking the divergence of the electrostatic field, we obtain Poisson's equation, that relates charge density and electric
potential
In the particular case of the empty space (ρ = 0) Poisson's equation reduces to Laplace's equation for the electric
potential.
Using a uniqueness theorem and showing that a potential satisfies Laplace's equation (second derivative of V should
be zero i.e. in free space) and the potential has the correct values at the boundaries, the potential is then uniquely
defined.
A potential that doesn't satisfy Laplace's equation together with the boundary condition is an invalid electrostatic
potential.
References
• Evans, L. C. (1998). Partial Differential Equations. Providence: American Mathematical Society.
ISBN 0-8218-0772-2.
• Petrovsky, I. G. (1967). Partial Differential Equations. Philadelphia: W. B. Saunders.
• Polyanin, A. D. (2002). Handbook of Linear Partial Differential Equations for Engineers and Scientists. Boca
Raton: Chapman & Hall/CRC Press. ISBN 1-58488-299-9.
• Sommerfeld, A. (1949). Partial Differential Equations in Physics. New York: Academic Press.
Laplace's equation 7
External links
• Hazewinkel, Michiel, ed. (2001), "Laplace equation" [1], Encyclopedia of Mathematics, Springer,
ISBN 978-1-55608-010-4
• Laplace Equation (particular solutions and boundary value problems) [2] at EqWorld: The World of Mathematical
Equations.
• Laplace Differential Equation [3] on PlanetMath
• Example initial-boundary value problems [4] using Laplace's equation from exampleproblems.com.
• Weisstein, Eric W., "Laplace’s Equation [5]", MathWorld.
• Module for Laplace’s Equation by John H. Mathews [6]
• Find out how boundary value problems governed by Laplace's equation may be solved numerically by boundary
element method [7]
References
[1] http:/ / www. encyclopediaofmath. org/ index. php?title=p/ l057470
[2] http:/ / eqworld. ipmnet. ru/ en/ solutions/ lpde/ lpde301. pdf
[3] http:/ / planetmath. org/ encyclopedia/ LaplaceDifferentialEquation. html
[4] http:/ / www. exampleproblems. com/ wiki/ index. php/ PDE:Laplaces_Equation
[5] http:/ / mathworld. wolfram. com/ LaplacesEquation. html
[6] http:/ / math. fullerton. edu/ mathews/ c2003/ DirichletProblemMod. html
[7] http:/ / www. ntu. edu. sg/ home/ mwtang/ bemsite. htm
Article Sources and Contributors 8
License
Creative Commons Attribution-Share Alike 3.0
//creativecommons.org/licenses/by-sa/3.0/