In mathematics, Poisson's equation is a partial differential equation of elliptic type with broad utility in mechanical engineering and theoretical physics. It arises, for instance, to describe the potential field caused by a given charge or mass density distribution; with the potential field known, one can then calculate the associated gravitational or electrostatic field. It is a generalization of Laplace's equation, which is also frequently seen in physics. The equation is named after the French mathematician, geometer, and physicist Siméon Denis Poisson.
Statement of the equation
Poisson's equation is
where is the Laplace operator, and and are real or complex-valued functions on a manifold. Usually, is given and is sought. When the manifold is Euclidean space, the Laplace operator is often denoted as ∇2 and so Poisson's equation is frequently written as
In three-dimensional Cartesian coordinates, it takes the form
When identically we obtain Laplace's equation.
Poisson's equation may be solved using a Green's function; a general exposition of the Green's function for Poisson's equation is given in the article on the screened Poisson equation. There are various methods for numerical solution. The relaxation method, an iterative algorithm, is one example.
In the case of a gravitational field g due to an attracting massive object of density ρ, Gauss's law for gravity in differential form can be used to obtain the corresponding Poisson equation for gravity,
Since the gravitational field is conservative (and irrotational), it can be expressed in terms of a scalar potential Φ,
Substituting into Gauss's law
yields Poisson's equation for gravity,
If the mass density is zero, Poisson's equation reduces to Laplace's equation. Using Green's Function, the potential at distance r from a central point mass m (i.e., the fundamental solution) is
which is equivalent to Newton's law of universal gravitation.
One of the cornerstones of electrostatics is setting up and solving problems described by the Poisson equation. Solving the Poisson equation amounts to finding the electric potential φ for a given charge distribution .
Substituting this into Gauss's law and assuming ε is spatially constant in the region of interest yields
In the absence of a changing magnetic field, B, Faraday's law of induction gives
where ∇× is the curl operator and t is the time.
The derivation of Poisson's equation under these circumstances is straightforward. Substituting the potential gradient for the electric field,
directly produces Poisson's equation for electrostatics, which is
Solving Poisson's equation for the potential requires knowing the charge density distribution. If the charge density is zero, then Laplace's equation results. If the charge density follows a Boltzmann distribution, then the Poisson-Boltzmann equation results. The Poisson–Boltzmann equation plays a role in the development of the Debye–Hückel theory of dilute electrolyte solutions.
Using Green's Function, the potential at distance r from a central point charge Q (ie: the Fundamental Solution) is:
which is Coulomb's law of electrostatics. (For historic reasons, and unlike gravity's model above, the factor appears here and not in Gauss's law.)
The above discussion assumes that the magnetic field is not varying in time. The same Poisson equation arises even if it does vary in time, as long as the Coulomb gauge is used. In this more general context, computing φ is no longer sufficient to calculate E, since E also depends on the magnetic vector potential A, which must be independently computed. See Maxwell's equation in potential formulation for more on φ and A in Maxwell's equations and how Poisson's equation is obtained in this case.
Potential of a Gaussian charge density
If there is a static spherically symmetric Gaussian charge density
where Q is the total charge, then the solution φ(r) of Poisson's equation,
is given by
where erf(x) is the error function.
This solution can be checked explicitly by evaluating ∇2φ.
Note that, for r much greater than σ, the erf function approaches unity and the potential φ(r) approaches the point charge potential
as one would expect. Furthermore the erf function approaches 1 extremely quickly as its argument increases; in practice for r > 3σ the relative error is smaller than one part in a thousand.
This technique reconstructs the implicit function f whose value is zero at the points pi and whose gradient at the points pi equals the normal vectors ni. The set of (pi, ni) is thus a sampling of a continuous vector ﬁeld V. The implicit function f is found by integrating the vector ﬁeld V. Since not every vector ﬁeld is the gradient of a function, the problem may or may not have a solution: the necessary and sufﬁcient condition for a smooth vector ﬁeld V to be the gradient of a function f is that the curl of V must be identically zero. In case this condition is difﬁcult to impose, it is still possible to perform a least-squares fit to minimize the difference between V and the gradient of f.
- Poisson Equation at EqWorld: The World of Mathematical Equations
- Evans, Lawrence C. (1998), Partial Differential Equations, Providence (RI): American Mathematical Society, ISBN 0-8218-0772-2
- Mathews, Jon; Walker, Robert L. (1970), Mathematical methods of physics (2nd ed.), New York: W. A. Benjamin, ISBN 0-8053-7002-1
- Polyanin, Andrei D. (2002), Handbook of Linear Partial Differential Equations for Engineers and Scientists, Boca Raton (FL): Chapman & Hall/CRC Press, ISBN 1-58488-299-9