Poisson's equation is an

_{f}'' = free charge volume

_{i}'' (a _{''i''}. Poisson's equation can be utilized to solve this problem with a technique called Poisson surface reconstruction.
The goal of this technique is to reconstruct an _{i}'' and whose gradient at the points ''p_{i}'' equals the normal vectors n_{''i''}. The set of (''p_{i}'', n_{''i''}) is thus modeled as a continuous _{i}'' onto the nodes of the particular staggered grid cell containing ''p_{i}''. Kazhdan and coauthors give a more accurate method of discretization using an adaptive finite difference grid, i.e. the cells of the grid are smaller (the grid is more finely divided) where there are more data points. They suggest implementing this technique with an adaptive

Poisson Equation

at EqWorld: The World of Mathematical Equations

Poisson's equation

on

elliptic partial differential equation
Second-order linear partial differential equations (PDEs) are classified as either elliptic, hyperbolic, or parabolic. Any second-order linear PDE in two variables can be written in the form
:Au_ + 2Bu_ + Cu_ + Du_x + Eu_y + Fu +G= 0,\,
where ...

of broad utility in theoretical physics
Theoretical physics is a branch of physics that employs mathematical models and abstractions of physical objects and systems to rationalize, explain and predict natural phenomena. This is in contrast to experimental physics, which uses experim ...

. For example, the solution to Poisson's equation is the potential field caused by a given electric charge or mass density distribution; with the potential field known, one can then calculate electrostatic or gravitational (force) field. It is a generalization of Laplace's equation, which is also frequently seen in physics. The equation is named after French mathematician and physicist Siméon Denis Poisson
Baron Siméon Denis Poisson FRS FRSE (; 21 June 1781 – 25 April 1840) was a French mathematician and physicist who worked on statistics, complex analysis, partial differential equations, the calculus of variations, analytical mechanics, electri ...

.
Statement of the equation

Poisson's equation is $$\backslash Delta\backslash varphi\; =\; f$$ where $\backslash Delta$ is theLaplace operator
In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols \nabla\cdot\nabla, \nabla^2 (where \nabla is the ...

, and $f$ and $\backslash varphi$ are real or complex
Complex commonly refers to:
* Complexity, the behaviour of a system whose components interact in multiple ways so possible interactions are difficult to describe
** Complex system, a system composed of many components which may interact with each ...

-valued functions on a manifold
In mathematics, a manifold is a topological space that locally resembles Euclidean space near each point. More precisely, an n-dimensional manifold, or ''n-manifold'' for short, is a topological space with the property that each point has a n ...

. Usually, $f$ is given and $\backslash varphi$ is sought. When the manifold is Euclidean space
Euclidean space is the fundamental space of geometry, intended to represent physical space. Originally, that is, in Euclid's ''Elements'', it was the three-dimensional space of Euclidean geometry, but in modern mathematics there are Euclidean ...

, the Laplace operator is often denoted as and so Poisson's equation is frequently written as $$\backslash nabla^2\; \backslash varphi\; =\; f.$$
In three-dimensional Cartesian coordinate
A Cartesian coordinate system (, ) in a plane is a coordinate system that specifies each point uniquely by a pair of numerical coordinates, which are the signed distances to the point from two fixed perpendicular oriented lines, measured in ...

s, it takes the form
$$\backslash left(\; \backslash frac\; +\; \backslash frac\; +\; \backslash frac\; \backslash right)\backslash varphi(x,y,z)\; =\; f(x,y,z).$$
When $f\; =\; 0$ identically we obtain Laplace's equation.
Poisson's equation may be solved using a Green's function
In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.
This means that if \operatorname is the linear different ...

:
$$\backslash varphi(\backslash mathbf)\; =\; -\; \backslash iiint\; \backslash frac\backslash ,\; \backslash mathrm^3\backslash !\; r\text{'},$$
where the integral is over all of space. 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, such as the relaxation method
In numerical mathematics, relaxation methods are iterative methods for solving systems of equations, including nonlinear systems.
Relaxation methods were developed for solving large sparse linear systems, which arose as finite-difference disc ...

, an iterative algorithm.
Newtonian gravity

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, $$\backslash nabla\backslash cdot\backslash mathbf\; =\; -4\backslash pi\; G\backslash rho\; ~.$$ Since the gravitational field is conservative (andirrotational
In vector calculus, a conservative vector field is a vector field that is the gradient of some function. A conservative vector field has the property that its line integral is path independent; the choice of any path between two points does not ...

), it can be expressed in terms of a scalar potential
In mathematical physics, scalar potential, simply stated, describes the situation where the difference in the potential energies of an object in two different positions depends only on the positions, not upon the path taken by the object in tra ...

''Φ'',
$$\backslash mathbf\; =\; -\backslash nabla\; \backslash phi\; ~.$$
Substituting into Gauss's law
$$\backslash nabla\backslash cdot(-\backslash nabla\; \backslash phi)\; =\; -\; 4\backslash pi\; G\; \backslash rho$$
yields Poisson's equation for gravity,
$$\backslash nabla^2\; \backslash phi\; =\; 4\backslash pi\; G\; \backslash rho.$$
If the mass density is zero, Poisson's equation reduces to Laplace's equation. The corresponding Green's function can be used to calculate the potential at distance from a central point mass (i.e., the fundamental solution
In mathematics, a fundamental solution for a linear partial differential operator is a formulation in the language of distribution theory of the older idea of a Green's function (although unlike Green's functions, fundamental solutions do not a ...

). In three dimensions the potential is
$$\backslash phi(r)\; =\; \backslash dfrac\; .$$
which is equivalent to Newton's law of universal gravitation
Newton's law of universal gravitation is usually stated as that every particle attracts every other particle in the universe with a force that is proportional to the product of their masses and inversely proportional to the square of the distanc ...

.
Electrostatics

One of the cornerstones ofelectrostatics
Electrostatics is a branch of physics that studies electric charges at rest ( static electricity).
Since classical times, it has been known that some materials, such as amber, attract lightweight particles after rubbing. The Greek word for amb ...

is setting up and solving problems described by the Poisson equation. Solving the Poisson equation amounts to finding the electric potential
The electric potential (also called the ''electric field potential'', potential drop, the electrostatic potential) is defined as the amount of work energy needed to move a unit of electric charge from a reference point to the specific point in ...

for a given charge distribution $\backslash rho\_f$.
The mathematical details behind Poisson's equation in electrostatics are as follows ( SI units are used rather than Gaussian units
Gaussian units constitute a metric system of physical units. This system is the most common of the several electromagnetic unit systems based on cgs (centimetre–gram–second) units. It is also called the Gaussian unit system, Gaussian-cgs uni ...

, which are also frequently used in electromagnetism
In physics, electromagnetism is an interaction that occurs between particles with electric charge. It is the second-strongest of the four fundamental interactions, after the strong force, and it is the dominant force in the interactions of ...

).
Starting with Gauss's law for electricity (also one of Maxwell's equations
Maxwell's equations, or Maxwell–Heaviside equations, are a set of coupled partial differential equations that, together with the Lorentz force law, form the foundation of classical electromagnetism, classical optics, and electric circuits.
...

) in differential form, one has
$$\backslash mathbf\; \backslash cdot\; \backslash mathbf\; =\; \backslash rho\_f$$
where $\backslash mathbf\; \backslash cdot$ is the divergence operator, D = electric displacement field
In physics, the electric displacement field (denoted by D) or electric induction is a vector field that appears in Maxwell's equations. It accounts for the effects of free and bound charge within materials. "D" stands for "displacement", as in ...

, and ''ρdensity
Density (volumetric mass density or specific mass) is the substance's mass per unit of volume. The symbol most often used for density is ''ρ'' (the lower case Greek letter rho), although the Latin letter ''D'' can also be used. Mathematicall ...

(describing charges brought from outside).
Assuming the medium is linear, isotropic, and homogeneous (see polarization density
In classical electromagnetism, polarization density (or electric polarization, or simply polarization) is the vector field that expresses the density of permanent or induced electric dipole moments in a dielectric material. When a dielectric is ...

), we have the constitutive equation
In physics and engineering, a constitutive equation or constitutive relation is a relation between two physical quantities (especially kinetic quantities as related to kinematic quantities) that is specific to a material or substance, and approx ...

,
$$\backslash mathbf\; =\; \backslash varepsilon\; \backslash mathbf$$
where is the permittivity
In electromagnetism, the absolute permittivity, often simply called permittivity and denoted by the Greek letter ''ε'' (epsilon), is a measure of the electric polarizability of a dielectric. A material with high permittivity polarizes more in r ...

of the medium and E is the electric field
An electric field (sometimes E-field) is the physical field that surrounds electrically charged particles and exerts force on all other charged particles in the field, either attracting or repelling them. It also refers to the physical field f ...

.
Substituting this into Gauss's law and assuming is spatially constant in the region of interest yields
$$\backslash mathbf\; \backslash cdot\; \backslash mathbf\; =\; \backslash frac\; ~.$$
where $\backslash rho$ is a total volume charge density. In electrostatics, we assume that there is no magnetic field (the argument that follows also holds in the presence of a constant magnetic field). Then, we have that
$$\backslash nabla\; \backslash times\; \backslash mathbf\; =\; 0,$$
where is the curl operator. This equation means that we can write the electric field as the gradient of a scalar function (called the electric potential
The electric potential (also called the ''electric field potential'', potential drop, the electrostatic potential) is defined as the amount of work energy needed to move a unit of electric charge from a reference point to the specific point in ...

), since the curl of any gradient is zero. Thus we can write,
$$\backslash mathbf\; =\; -\backslash nabla\; \backslash varphi,$$
where the minus sign is introduced so that is identified as the electric potential energy per unit charge.
The derivation of Poisson's equation under these circumstances is straightforward. Substituting the potential gradient for the electric field,
$$\backslash nabla\; \backslash cdot\; \backslash mathbf\; =\; \backslash nabla\; \backslash cdot\; (\; -\; \backslash nabla\; \backslash varphi\; )\; =\; -\; ^2\; \backslash varphi\; =\; \backslash frac,$$
directly produces Poisson's equation for electrostatics, which is
$$\backslash nabla^2\; \backslash varphi\; =\; -\backslash frac.$$
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 from a central point charge (i.e., the Fundamental Solution) is:
$$\backslash varphi(r)\; =\; \backslash frac\; .$$
which is Coulomb's law of electrostatics. (For historic reasons, and unlike gravity's model above, the $4\; \backslash pi$ 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
In classical electromagnetism, magnetic vector potential (often called A) is the vector quantity defined so that its curl is equal to the magnetic field: \nabla \times \mathbf = \mathbf. Together with the electric potential ''φ'', the magnetic ve ...

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 symmetricGaussian
Carl Friedrich Gauss (1777–1855) is the eponym of all of the topics listed below.
There are over 100 topics all named after this German mathematician and scientist, all in the fields of mathematics, physics, and astronomy. The English eponymo ...

charge density
$$\backslash rho\_f(r)\; =\; \backslash frac\backslash ,e^,$$
where is the total charge, then the solution of Poisson's equation,
$$^2\; \backslash varphi\; =\; -\; ,$$
is given by
$$\backslash varphi(r)\; =\; \backslash frac\; \backslash frac\; \backslash ,\; \backslash operatorname\backslash left(\backslash frac\backslash right)$$
where is the error function
In mathematics, the error function (also called the Gauss error function), often denoted by , is a complex function of a complex variable defined as:
:\operatorname z = \frac\int_0^z e^\,\mathrm dt.
This integral is a special (non-elementar ...

.
This solution can be checked explicitly by evaluating .
Note that, for much greater than , the erf function approaches unity and the potential approaches the point charge
A point particle (ideal particle or point-like particle, often spelled pointlike particle) is an idealization of particles heavily used in physics. Its defining feature is that it lacks spatial extension; being dimensionless, it does not take u ...

potential
$$\backslash varphi\; \backslash approx\; \backslash frac\; \backslash frac\; ,$$
as one would expect. Furthermore, the error function approaches 1 extremely quickly as its argument increases; in practice for the relative error is smaller than one part in a thousand.
Surface reconstruction

Surface reconstruction is an inverse problem. The goal is to digitally reconstruct a smooth surface based on a large number of points ''ppoint cloud
Point or points may refer to:
Places
* Point, Lewis, a peninsula in the Outer Hebrides, Scotland
* Point, Texas, a city in Rains County, Texas, United States
* Point, the NE tip and a ferry terminal of Lismore, Inner Hebrides, Scotland
* Poin ...

) where each point also carries an estimate of the local surface normal
In geometry, a normal is an object such as a line, ray, or vector that is perpendicular to a given object. For example, the normal line to a plane curve at a given point is the (infinite) line perpendicular to the tangent line to the curve a ...

nimplicit function
In mathematics, an implicit equation is a relation of the form R(x_1, \dots, x_n) = 0, where is a function of several variables (often a polynomial). For example, the implicit equation of the unit circle is x^2 + y^2 - 1 = 0.
An implicit fun ...

''f'' whose value is zero at the points ''pvector
Vector most often refers to:
*Euclidean vector, a quantity with a magnitude and a direction
*Vector (epidemiology), an agent that carries and transmits an infectious pathogen into another living organism
Vector may also refer to:
Mathematic ...

field V. The implicit function ''f'' is found by integrating the vector field V. Since not every vector field is the gradient
In vector calculus, the gradient of a scalar-valued differentiable function of several variables is the vector field (or vector-valued function) \nabla f whose value at a point p is the "direction and rate of fastest increase". If the gra ...

of a function, the problem may or may not have a solution: the necessary and sufficient condition for a smooth vector field V to be the gradient of a function ''f'' is that the curl
cURL (pronounced like "curl", UK: , US: ) is a computer software project providing a library (libcurl) and command-line tool (curl) for transferring data using various network protocols. The name stands for "Client URL".
History
cURL was f ...

of V must be identically zero. In case this condition is difficult to impose, it is still possible to perform a least-squares
The method of least squares is a standard approach in regression analysis to approximate the solution of overdetermined systems (sets of equations in which there are more equations than unknowns) by minimizing the sum of the squares of the res ...

fit to minimize the difference between V and the gradient of ''f''.
In order to effectively apply Poisson's equation to the problem of surface reconstruction, it is necessary to find a good discretization of the vector field V. The basic approach is to bound the data with a finite difference grid. For a function valued at the nodes of such a grid, its gradient can be represented as valued on staggered grids, i.e. on grids whose nodes lie in between the nodes of the original grid. It is convenient to define three staggered grids, each shifted in one and only one direction corresponding to the components of the normal data. On each staggered grid we perform rilinear interpolationon the set of points. The interpolation weights are then used to distribute the magnitude of the associated component of ''noctree
An octree is a tree data structure in which each internal node has exactly eight children. Octrees are most often used to partition a three-dimensional space by recursively subdividing it into eight octants. Octrees are the three-dimensional ana ...

.
Fluid dynamics

For the incompressible Navier–Stokes equations, given by: $$\backslash begin\; +\; (\backslash cdot\backslash nabla)\; \&=\; -\backslash nabla\; p\; +\; \backslash nu\backslash Delta\; +\; \backslash \backslash \; \backslash nabla\backslash cdot\; \&=\; 0\; \backslash end$$ The equation for the pressure field $p$ is an example of a nonlinear Poisson equation: $$\backslash begin\; \backslash Delta\; p\; \&=\; -\backslash rho\; \backslash nabla\backslash cdot(\backslash cdot\; \backslash nabla\; )\; \backslash \backslash \; \&=\; -\backslash rho\backslash ,\; \backslash mathrm\backslash big((\backslash nabla\; )\; (\backslash nabla\; )\backslash big).\; \backslash end$$ Notice that the above trace is not sign-definite.See also

*Discrete Poisson equation
In mathematics, the discrete Poisson equation is the finite difference analog of the Poisson equation. In it, the discrete Laplace operator takes the place of the Laplace operator. The discrete Poisson equation is frequently used in numerical a ...

* Poisson–Boltzmann equation
The Poisson–Boltzmann equation is a useful equation in many settings, whether it be to understand physiology, physiological interfaces, polymer science, electron interactions in a semiconductor, or more. It aims to describe the distribution of th ...

* Helmholtz equation
* Uniqueness theorem for Poisson's equation
* Weak formulation Weak formulations are important tools for the analysis of mathematical equations that permit the transfer of concepts of linear algebra to solve problems in other fields such as partial differential equations. In a weak formulation, equations or con ...

References

Further reading

* * *External links

* {{springer, title=Poisson equation, id=p/p073290Poisson Equation

at EqWorld: The World of Mathematical Equations

Poisson's equation

on

PlanetMath
PlanetMath is a free, collaborative, mathematics online encyclopedia. The emphasis is on rigour, openness, pedagogy, real-time content, interlinked content, and also community of about 24,000 people with various maths interests. Intended to be ...

.
Potential theory
Partial differential equations
Electrostatics
Mathematical physics
Equations of physics
Electromagnetism