Poisson's equation

Last updated
Simeon Denis Poisson Simeon Poisson.jpg
Siméon Denis Poisson

Poisson's equation is an elliptic partial differential equation of broad utility in theoretical physics. 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 the corresponding 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 who published it in 1823. [1] [2]

Contents

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: 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, an iterative algorithm.

Applications in physics and engineering

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. Gauss's law for gravity is

Since the gravitational field is conservative (and irrotational), it can be expressed in terms of a scalar potential ϕ:

Substituting this into Gauss's law, yields Poisson's equation for gravity:

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 r from a central point mass m (i.e., the fundamental solution). In three dimensions the potential is which is equivalent to Newton's law of universal gravitation.

Electrostatics

Many problems in electrostatics are governed by the Poisson equation, which relates the electric potential φ to the free charge density , such as those found in conductors.

The mathematical details of Poisson's equation, commonly expressed in SI units (as opposed to Gaussian units), describe how the distribution of free charges generates the electrostatic potential in a given Region (mathematics).

Starting with Gauss's law for electricity (also one of Maxwell's equations) in differential form, one has where is the divergence operator, D is the electric displacement field, and ρf is the free-charge density (describing charges brought from outside).

Assuming the medium is linear, isotropic, and homogeneous (see polarization density), we have the constitutive equation where ε is the permittivity of the medium, and E is the electric field.

Substituting this into Gauss's law and assuming that ε is spatially constant in the region of interest yields In electrostatics, we assume that there is no magnetic field (the argument that follows also holds in the presence of a constant magnetic field). [3] Then, we have that 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), since the curl of any gradient is zero. Thus we can write where the minus sign is introduced so that φ is identified as the electric potential energy per unit charge. [4]

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

Specifying the 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 a Green's function, the potential at distance r from a central point charge Q (i.e., the fundamental solution) is which is Coulomb's law of electrostatics. (For historical 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 class of cases, 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 an appropriate 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. [5] This solution can be checked explicitly by evaluating 2φ.

Note that for r much greater than σ, approaches unity, [6] and the potential φ(r) approaches the point-charge potential, as one would expect. Furthermore, the error 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. [6]

Surface reconstruction

Surface reconstruction is an inverse problem. The goal is to digitally reconstruct a smooth surface based on a large number of points pi (a point cloud) where each point also carries an estimate of the local surface normal ni. [7] Poisson's equation can be utilized to solve this problem with a technique called Poisson surface reconstruction. [8]

The goal of this technique is to reconstruct an 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 modeled as a continuous vector field V. The implicit function f is found by integrating the vector field V. Since not every vector field is the gradient 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 of V must be identically zero. In case this condition is difficult to impose, it is still possible to perform a least-squares 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 trilinear interpolation on the set of points. The interpolation weights are then used to distribute the magnitude of the associated component of ni onto the nodes of the particular staggered grid cell containing pi. 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. [8] They suggest implementing this technique with an adaptive octree.

Fluid dynamics

For the incompressible Navier–Stokes equations, given by

The equation for the pressure field is an example of a nonlinear Poisson equation: Notice that the above trace is not sign-definite.

See also

Related Research Articles

<span class="mw-page-title-main">Laplace's equation</span> Second-order partial differential equation

In mathematics and physics, 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 or where is the Laplace operator, is the divergence operator, is the gradient operator, and is a twice-differentiable real-valued function. The Laplace operator therefore maps a scalar function to another scalar function.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

<span class="mw-page-title-main">Potential flow</span> Velocity field as the gradient of a scalar function

In fluid dynamics, potential flow or irrotational flow refers to a description of a fluid flow with no vorticity in it. Such a description typically arises in the limit of vanishing viscosity, i.e., for an inviscid fluid and with no vorticity present in the flow.

<span class="mw-page-title-main">Electric potential</span> Line integral of the electric field

Electric potential is defined as the amount of work/energy needed per unit of electric charge to move the charge from a reference point to a specific point in an electric field. More precisely, the electric potential is the energy per unit charge for a test charge that is so small that the disturbance of the field under consideration is negligible. The motion across the field is supposed to proceed with negligible acceleration, so as to avoid the test charge acquiring kinetic energy or producing radiation. By definition, the electric potential at the reference point is zero units. Typically, the reference point is earth or a point at infinity, although any point can be used.

<span class="mw-page-title-main">Gauss's law</span> Foundational law of electromagnetism relating electric field and charge distributions

In physics, Gauss's law, also known as Gauss's flux theorem, is one of Maxwell's equations. It is an application of the divergence theorem, and it relates the distribution of electric charge to the resulting electric field.

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 , (where is the nabla operator), or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).

In fluid dynamics, Stokes' law gives the frictional force – also called drag force – exerted on spherical objects moving at very small Reynolds numbers in a viscous fluid. It was derived by George Gabriel Stokes in 1851 by solving the Stokes flow limit for small Reynolds numbers of the Navier–Stokes equations.

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

<span class="mw-page-title-main">Green's function</span> Impulse response of an inhomogeneous linear differential operator

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.

<span class="mw-page-title-main">Electrostatics</span> Study of stationary or slow-moving electric charges

Electrostatics is a branch of physics that studies slow-moving or stationary electric charges.

In physics and mathematics, the Helmholtz decomposition theorem or the fundamental theorem of vector calculus states that certain differentiable vector fields can be resolved into the sum of an irrotational (curl-free) vector field and a solenoidal (divergence-free) vector field. In physics, often only the decomposition of sufficiently smooth, rapidly decaying vector fields in three dimensions is discussed. It is named after Hermann von Helmholtz.

<span class="mw-page-title-main">Charge density</span> Electric charge per unit length, area or volume

In electromagnetism, charge density is the amount of electric charge per unit length, surface area, or volume. Volume charge density is the quantity of charge per unit volume, measured in the SI system in coulombs per cubic meter (C⋅m−3), at any point in a volume. Surface charge density (σ) is the quantity of charge per unit area, measured in coulombs per square meter (C⋅m−2), at any point on a surface charge distribution on a two dimensional surface. Linear charge density (λ) is the quantity of charge per unit length, measured in coulombs per meter (C⋅m−1), at any point on a line charge distribution. Charge density can be either positive or negative, since electric charge can be either positive or negative.

The method of image charges is a basic problem-solving tool in electrostatics. The name originates from the replacement of certain elements in the original layout with fictitious charges, which replicates the boundary conditions of the problem.

<span class="mw-page-title-main">Inhomogeneous electromagnetic wave equation</span> Equation in physics

In electromagnetism and applications, an inhomogeneous electromagnetic wave equation, or nonhomogeneous electromagnetic wave equation, is one of a set of wave equations describing the propagation of electromagnetic waves generated by nonzero source charges and currents. The source terms in the wave equations make the partial differential equations inhomogeneous, if the source terms are zero the equations reduce to the homogeneous electromagnetic wave equations, which follow from Maxwell's equations.

In physics, the Green's function for the Laplacian in three variables is used to describe the response of a particular type of physical system to a point source. In particular, this Green's function arises in systems that can be described by Poisson's equation, a partial differential equation (PDE) of the form where is the Laplace operator in , is the source term of the system, and is the solution to the equation. Because is a linear differential operator, the solution to a general system of this type can be written as an integral over a distribution of source given by : where the Green's function for Laplacian in three variables describes the response of the system at the point to a point source located at : and the point source is given by , the Dirac delta function.

<span class="mw-page-title-main">Mathematical descriptions of the electromagnetic field</span> Formulations of electromagnetism

There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.

<span class="mw-page-title-main">Retarded potential</span> Type of potential in electrodynamics

In electrodynamics, the retarded potentials are the electromagnetic potentials for the electromagnetic field generated by time-varying electric current or charge distributions in the past. The fields propagate at the speed of light c, so the delay of the fields connecting cause and effect at earlier and later times is an important factor: the signal takes a finite time to propagate from a point in the charge or current distribution to another point in space, see figure below.

In mathematics, vector spherical harmonics (VSH) are an extension of the scalar spherical harmonics for use with vector fields. The components of the VSH are complex-valued functions expressed in the spherical coordinate basis vectors.

In physics, Gauss's law for gravity, also known as Gauss's flux theorem for gravity, is a law of physics that is equivalent to Newton's law of universal gravitation. It is named after Carl Friedrich Gauss. It states that the flux of the gravitational field over any closed surface is proportional to the mass enclosed. Gauss's law for gravity is often more convenient to work from than Newton's law.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

References

  1. Jackson, Julia A.; Mehl, James P.; Neuendorf, Klaus K. E., eds. (2005), Glossary of Geology, American Geological Institute, Springer, p. 503, ISBN   9780922152766
  2. Poisson (1823). "Mémoire sur la théorie du magnétisme en mouvement" [Memoir on the theory of magnetism in motion]. Mémoires de l'Académie Royale des Sciences de l'Institut de France (in French). 6: 441–570. From p. 463: "Donc, d'après ce qui précède, nous aurons enfin:selon que le point M sera situé en dehors, à la surface ou en dedans du volume que l'on considère." (Thus, according to what preceded, we will finally have: depending on whether the point M is located outside, on the surface of, or inside the volume that one is considering.) V is defined (p. 462) as where, in the case of electrostatics, the integral is performed over the volume of the charged body, the coordinates of points that are inside or on the volume of the charged body are denoted by , is a given function of and in electrostatics, would be a measure of charge density, and is defined as the length of a radius extending from the point M to a point that lies inside or on the charged body. The coordinates of the point M are denoted by and denotes the value of (the charge density) at M.
  3. Griffiths, D. J. (2017). Introduction to Electrodynamics (4th ed.). Cambridge University Press. pp. 77–78.
  4. Griffiths, D. J. (2017). Introduction to Electrodynamics (4th ed.). Cambridge University Press. pp. 83–84.
  5. Salem, M.; Aldabbagh, O. (2024). "Numerical Solution to Poisson's Equation for Estimating Electrostatic Properties Resulting from an Axially Symmetric Gaussian Charge Density Distribution". Mathematics. 12 (13): 1948. doi: 10.3390/math12131948 .
  6. 1 2 Oldham, K. B.; Myland, J. C.; Spanier, J. (2008). "The Error Function erf(x) and Its Complement erfc(x)". An Atlas of Functions. New York, NY: Springer. pp. 405–415. doi:10.1007/978-0-387-48807-3_41. ISBN   978-0-387-48806-6.
  7. Calakli, Fatih; Taubin, Gabriel (2011). "Smooth Signed Distance Surface Reconstruction" (PDF). Pacific Graphics. 30 (7).
  8. 1 2 Kazhdan, Michael; Bolitho, Matthew; Hoppe, Hugues (2006). "Poisson surface reconstruction". Proceedings of the fourth Eurographics symposium on Geometry processing (SGP '06). Eurographics Association, Aire-la-Ville, Switzerland. pp. 61–70. ISBN   3-905673-36-3.

Further reading