Poisson's equation (Aakash Dave)



Poisson's equation

In mathematics, Poisson's equation is a partial differential equation with broad utility in electrostatics, mechanical engineering and theoretical physics. It is named after the French mathematician, geometer and physicist Siméon-Denis Poisson. The Poisson equation is
\Delta\varphi=f
where Δ is the Laplace operator, and f and φ are real or complex-valued functions on a manifold. When the manifold is Euclidean space, the Laplace operator is often denoted as {\nabla}^2 and so Poisson's equation is frequently written as
{\nabla}^2 \varphi = f.
In three-dimensional Cartesian coordinates, it takes the form
\left( \frac{\partial^2}{\partial x^2} + \frac{\partial^2}{\partial y^2} + \frac{\partial^2}{\partial z^2} \right)\varphi(x,y,z) = f(x,y,z).
For vanishing f, this equation becomes Laplace's equation
\Delta \varphi = 0. \!
The Poisson equation may be solved using a Green's function; a general exposition of the Green's function for the Poisson 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.

Electrostatics

One of the cornerstones of electrostatics is the posing and solving of problems that are described by the Poisson equation. Finding φ for some given f is an important practical problem, since this is the usual way to find the electric potential for a given charge distribution.
The derivation of Poisson's equation in electrostatics follows. SI units are used and Euclidean space is assumed.
Starting with Gauss' law for electricity (also part of Maxwell's equations) in a differential control volume, we have:
\mathbf{\nabla} \cdot \mathbf{D} = \rho_f
\mathbf{\nabla} \cdot is the divergence operator.
\mathbf{D} is the electric displacement field.
\rho_f \, is the free charge density (describing charges brought from outside).
Assuming the medium is linear, isotropic, and homogeneous (see polarization density), then:
\mathbf{D} = \varepsilon \mathbf{E}
\varepsilon is the permittivity of the medium.
\mathbf{E} is the electric field.
By substitution and division, we have:
\mathbf{\nabla} \cdot \mathbf{E} = \frac{\rho_f}{\varepsilon}
In the absence of a changing magnetic field, \mathbf{B}, Faraday's law of induction gives:
\nabla \times \mathbf{E} = -\dfrac{\partial \mathbf{B}} {\partial t} = 0
\nabla \times is the curl operator.
t is time.
Since the curl of the electric field is zero, it is defined by a scalar electric potential field, \varphi (see Helmholtz decomposition).
\mathbf{E} = -\nabla \varphi
Eliminating \mathbf{E} by substitution, we have a form of the Poisson equation:
\nabla \cdot \nabla \varphi = {\nabla}^2 \varphi = -\frac{\rho_f}{\varepsilon}.
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.
(Note: Although 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. However, in this more general context, computing \varphi is no longer sufficient to calculate \mathbf{E}, since the latter also depends on the magnetic vector potential, which must be independently computed.)

Potential of a Gaussian charge density
If there is a static spherically symmetric Gaussian charge density ρf(r):
\rho_f(r) = \frac{Q}{\sigma^3\sqrt{2\pi}^3}\,e^{-r^2/(2\sigma^2)},
where Q is the total charge, then the solution φ (r) of Poisson's equation,
{\nabla}^2 \varphi = - { \rho_f \over \varepsilon } ,
is given by
\varphi(r) = { 1 \over 4 \pi \varepsilon } \frac{Q}{r}\,\mbox{erf}\left(\frac{r}{\sqrt{2}\sigma}\right)
where erf(x) is the error function. This solution can be checked explicitly by evaluating {\nabla}^2 \varphi. Note that, for r much greater than σ, the erf function approaches unity and the potential φ (r) approaches the point charge { 1 \over 4 \pi \varepsilon_0 } {Q \over r} , 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.

 


 

Subscribe to: Posts (Atom)