Helmholtz decomposition

From Wikipedia, the free encyclopedia
Jump to: navigation, search

In physics and mathematics, in the area of vector calculus, Helmholtz's theorem,[1][2] also known as the fundamental theorem of vector calculus,[3][4][5][6][7][8][9] states that any sufficiently smooth, rapidly decaying vector field in three dimensions can be resolved into the sum of an irrotational (curl-free) vector field and a solenoidal (divergence-free) vector field; this is known as the Helmholtz decomposition. It is named after Hermann von Helmholtz.[10]

This implies that any such vector field F can be considered to be generated by a pair of potentials: a scalar potential φ and a vector potential A.

Statement of the theorem[edit]

Let F be a vector field on a bounded domain V in R3, which is twice continuously differentiable, and let S be the surface that encloses the domain V. Then F can be decomposed into a curl-free component and a divergence-free component:[11]

\mathbf{F}=-\boldsymbol{\nabla}\Phi+\boldsymbol{\nabla}\times\mathbf{A},

where

\Phi\left(\mathbf{r}\right)=\frac{1}{4\pi}\int_{V}\frac{\boldsymbol{\nabla}'\cdot\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\frac{1}{4\pi}\oint_{S}\mathbf{\hat{n}}'\cdot\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'


\mathbf{A}\left(\mathbf{r}\right)=\frac{1}{4\pi}\int_{V}\frac{\boldsymbol{\nabla}'\times\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\frac{1}{4\pi}\oint_{S}\mathbf{\hat{n}}'\times\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'


If V is R3 itself (unbounded), and F vanishes faster than 1/r as r\rightarrow\infty, then the second component of both scalar and vector potential are zero. That is,[12]

\Phi\left(\mathbf{r}\right)=\frac{1}{4\pi}\int_{\text{all space}}\frac{\boldsymbol{\nabla}'\cdot\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'


\mathbf{A}\left(\mathbf{r}\right)=\frac{1}{4\pi}\int_{\text{all space}}\frac{\boldsymbol{\nabla}'\times\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'

Derivation[edit]

Suppose we have a vector function \mathbf{F}\left(\mathbf{r}\right) of which we know the curl, \boldsymbol{\nabla}\times\mathbf{F}, and the divergence, \boldsymbol{\nabla}\cdot\mathbf{F}, in the domain and the fields on the boundary. Writing the function using delta function in the form

\delta^3\left(\mathbf{r}-\mathbf{r}'\right)=-\frac{1}{4\pi}\nabla^{2}\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|},
\mathbf{F}\left(\mathbf{r}\right)=\int_{V}\mathbf{F}\left(\mathbf{r}'\right)\delta^3\left(\mathbf{r}-\mathbf{r}'\right)\mathrm{d}V'=\int_{V}\mathbf{F}\left(\mathbf{r}'\right)\left(-\frac{1}{4\pi}\nabla^{2}\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\right)\mathrm{d}V'=-\frac{1}{4\pi}\nabla^{2}\int_{V}\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'.

Using the identity

\nabla^{2}\mathbf{a}=\boldsymbol{\nabla}\left(\boldsymbol{\nabla}\cdot\mathbf{a}\right)-\boldsymbol{\nabla}\times\left(\boldsymbol{\nabla}\times\mathbf{a}\right).

We get,

\mathbf{F}\left(\mathbf{r}\right)=-\frac{1}{4\pi}\left[\boldsymbol{\nabla}\left(\boldsymbol{\nabla}\cdot\int_{V}\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'\right)-\boldsymbol{\nabla}\times\left(\boldsymbol{\nabla}\times\int_{V}\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'\right)\right]
=-\frac{1}{4\pi}\left[\boldsymbol{\nabla}\left(\int_{V}\mathbf{F}\left(\mathbf{r}'\right)\cdot\boldsymbol{\nabla}\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'\right)+\boldsymbol{\nabla}\times\left(\int_{V}\mathbf{F}\left(\mathbf{r}'\right)\times\boldsymbol{\nabla}\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'\right)\right].

Noting that, \boldsymbol{\nabla}\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}=-\boldsymbol{\nabla}'\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|},\, we can rewrite the last expression as

\mathbf{F}\left(\mathbf{r}\right)=-\frac{1}{4\pi}\left[-\boldsymbol{\nabla}\left(\int_{V}\mathbf{F}\left(\mathbf{r}'\right)\cdot\boldsymbol{\nabla}'\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'\right)-\boldsymbol{\nabla}\times\left(\int_{V}\mathbf{F}\left(\mathbf{r}'\right)\times\boldsymbol{\nabla}'\frac{1}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'\right)\right].


Then using the vectorial identities

\mathbf{a}\cdot\boldsymbol{\nabla}\psi=-\psi\left(\boldsymbol{\nabla}\cdot\mathbf{a}\right)+\boldsymbol{\nabla}\cdot\left(\psi\mathbf{a}\right)

and

\mathbf{a}\times\boldsymbol{\nabla}\psi=\psi\left(\boldsymbol{\nabla}\times\mathbf{a}\right)-\boldsymbol{\nabla}\times\left(\psi\mathbf{a}\right),

we get

\mathbf{F}\left(\mathbf{r}\right)=-\frac{1}{4\pi}\left[-\boldsymbol{\nabla}\left(
-\int_{V}\frac{\boldsymbol{\nabla}'\cdot\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
+\int_{V}\boldsymbol{\nabla}'\cdot\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
\right)-\boldsymbol{\nabla}\times\left(
\int_{V}\frac{\boldsymbol{\nabla}'\times\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\int_{V}\boldsymbol{\nabla}'\times\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
\right)\right].


Take advantage of the divergence theorem, the equation can be rewritten as

\mathbf{F}\left(\mathbf{r}\right)=-\frac{1}{4\pi}\left[-\boldsymbol{\nabla}\left(
-\int_{V}\frac{\boldsymbol{\nabla}'\cdot\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
+\oint_{S}\mathbf{\hat{n}}'\cdot\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'
\right)-\boldsymbol{\nabla}\times\left(
\int_{V}\frac{\boldsymbol{\nabla}'\times\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\oint_{S}\mathbf{\hat{n}}'\times\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'
\right)\right]
=
-\boldsymbol{\nabla}\left[\frac{1}{4\pi}\int_{V}\frac{\boldsymbol{\nabla}'\cdot\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\frac{1}{4\pi}\oint_{S}\mathbf{\hat{n}}'\cdot\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'\right]
+\boldsymbol{\nabla}\times\left[\frac{1}{4\pi}\int_{V}\frac{\boldsymbol{\nabla}'\times\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\frac{1}{4\pi}\oint_{S}\mathbf{\hat{n}}'\times\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'\right].

with outward surface normal  \mathbf{\hat{n}}' .

Define

\Phi\left(\mathbf{r}\right)\equiv\frac{1}{4\pi}\int_{V}\frac{\boldsymbol{\nabla}'\cdot\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\frac{1}{4\pi}\oint_{S}\mathbf{\hat{n}}'\cdot\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'


\mathbf{A}\left(\mathbf{r}\right)\equiv\frac{1}{4\pi}\int_{V}\frac{\boldsymbol{\nabla}'\times\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
-\frac{1}{4\pi}\oint_{S}\mathbf{\hat{n}}'\times\frac{\mathbf{F}\left(\mathbf{r}'\right)}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}S'


Hence

\mathbf{F}=-\boldsymbol{\nabla}\Phi+\boldsymbol{\nabla}\times\mathbf{A}


Another derivation from the Fourier transform[edit]

Write F as a Fourier transform :

\vec{\mathbf{F}}(\vec{r}) = \iiint \vec{\mathbf{G}}(\vec{\omega}) e^{\displaystyle i \, \vec{\omega} \cdot \vec{r}} d\vec{\omega}

The Fourier transform of a scalar field is a scalar field, and the Fourier transform of a vector field is a vector field of same dimension.

Now consider the following scalar and vector fields :

\begin{array}{lll} G_\Phi(\vec{\omega}) =   i\, \frac{\displaystyle \vec{\mathbf{G}}(\vec{\omega}) \cdot \vec{\omega}}{||\vec{\omega}||^2} & \quad\quad &
\vec{\mathbf{G}}_\mathbf{A}(\vec{\omega}) = i\, \vec{\omega} \times \left( \vec{\mathbf{G}}(\vec{\omega}) + i G_\Phi(\vec{\omega}) \, \vec{\omega} \right)  \\
 && \\
\Phi(\vec{r}) = \displaystyle \iiint G_\Phi(\vec{\omega}) e^{\displaystyle i \, \vec{\omega} \cdot \vec{r}} d\vec{\omega} & & \vec{\mathbf{A}}(\vec{r}) = \displaystyle \iiint  \vec{\mathbf{G}}_\mathbf{A}(\vec{\omega}) e^{\displaystyle i \, \vec{\omega} \cdot \vec{r}} d\vec{\omega} \end{array}

Hence

 \vec{\mathbf{G}}(\vec{\omega}) = - i \,\vec{\omega} \, G_\Phi(\vec{\omega})  + i \, \vec{\omega} \times \vec{\mathbf{G}}_\mathbf{A}(\vec{\omega})

\begin{array}{lll}\vec{\mathbf{F}}(\vec{r}) &=& \displaystyle - \iiint i \, \vec{\omega}\, G_\Phi(\vec{\omega})  \, e^{\displaystyle i \, \vec{\omega} \cdot \vec{r}} d\vec{\omega}
+  \iiint i \, \vec{\omega} \times \vec{\mathbf{G}}_\mathbf{A}(\vec{\omega}) e^{\displaystyle i \, \vec{\omega} \cdot \vec{r}} d\vec{\omega} \\
&=& - \boldsymbol{\nabla} \Phi(\vec{r}) +  \boldsymbol{\nabla} \times \vec{\mathbf{A}}(\vec{r})
\end{array}

Fields with prescribed divergence and curl[edit]

The term "Helmholtz Theorem" can also refer to the following. Let C be a solenoidal vector field and d a scalar field on R3 which are sufficiently smooth and which vanish faster than 1/r2 at infinity. Then there exists a vector field F such that

\nabla \cdot \mathbf{F} = d and  \nabla \times \mathbf{F} = \mathbf{C};

if additionally the vector field F vanishes as r → ∞, then F is unique.[12]

In other words, a vector field can be constructed with both a specified divergence and a specified curl, and if it also vanishes at infinity, it is uniquely specified by its divergence and curl. This theorem is of great importance in electrostatics, since Maxwell's equations for the electric and magnetic fields in the static case are of exactly this type.[12] The proof is by a construction generalizing the one given above: we set

\mathbf{F} = - \nabla\,(\mathcal{G} (d)) + \nabla \times (\mathcal{G}(\mathbf{C})),

where \mathcal{G} represents the Newtonian potential operator. (When acting on a vector field, such as ∇ × F, it is defined to act on each component.)

Differential forms[edit]

The Hodge decomposition is closely related to the Helmholtz decomposition, generalizing from vector fields on R3 to differential forms on a Riemannian manifold M. Most formulations of the Hodge decomposition require M to be compact.[13] Since this is not true of R3, the Hodge decomposition theorem is not strictly a generalization of the Helmholtz theorem. However, the compactness restriction in the usual formulation of the Hodge decomposition can be replaced by suitable decay assumptions at infinity on the differential forms involved, giving a proper generalization of the Helmholtz theorem.

Weak formulation[edit]

The Helmholtz decomposition can also be generalized by reducing the regularity assumptions (the need for the existence of strong derivatives). Suppose Ω is a bounded, simply-connected, Lipschitz domain. Every square-integrable vector field u ∈ (L2(Ω))3 has an orthogonal decomposition:

\mathbf{u}=\nabla\varphi+\nabla \times \mathbf{A}

where \varphi is in the Sobolev space H1(Ω) of square-integrable functions on Ω whose partial derivatives defined in the distribution sense are square integrable, and A ∈ H(curl,Ω), the Sobolev space of vector fields consisting of square integrable vector fields with square integrable curl.

For a slightly smoother vector field u ∈ H(curl,Ω), a similar decomposition holds:

\mathbf{u}=\nabla\varphi+\mathbf{v}

where \varphi ∈ H1(Ω) and v ∈ (H1(Ω))d.

Longitudinal and transverse fields[edit]

A terminology often used in physics refers to the curl-free component of a vector field as the longitudinal component and the divergence-free component as the transverse component.[14] This terminology comes from the following construction: Compute the three-dimensional Fourier transform of the vector field F. Then decompose this field, at each point k, into two components, one of which points longitudinally, i.e. parallel to k, the other of which points in the transverse direction, i.e. perpendicular to k. So far, we have

\mathbf{F}\left(\mathbf{k}\right) = \mathbf{F}_t\left(\mathbf{k}\right) + \mathbf{F}_l\left(\mathbf{k}\right)
\mathbf{k} \cdot \mathbf{F}_t\left(\mathbf{k}\right) = 0.
\mathbf{k}\times \mathbf{F}_l\left(\mathbf{k}\right) = \mathbf{0}.


Now we apply an inverse Fourier transform to each of these components. Using properties of Fourier transforms, we derive:

\mathbf{F}\left(\mathbf{r}\right) = \mathbf{F}_t\left(\mathbf{r}\right)+\mathbf{F}_l\left(\mathbf{r}\right)
\boldsymbol{\nabla} \cdot \mathbf{F}_t\left(\mathbf{r}\right) = 0
\boldsymbol{\nabla}\times \mathbf{F}_l\left(\mathbf{r}\right) = \mathbf{0}


Since \boldsymbol{\nabla}\times\left(\boldsymbol{\nabla}\Phi\right)=0 and \boldsymbol{\nabla}\cdot\left(\boldsymbol{\nabla}\times\mathbf{A}\right)=0,

we can get

\mathbf{F}_t=\boldsymbol{\nabla}\times\mathbf{A}=\frac{1}{4\pi}\boldsymbol{\nabla}\times\int_V\frac{\boldsymbol{\nabla}'\times\mathbf{F}}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'
\mathbf{F}_l=-\boldsymbol{\nabla}\Phi=-\frac{1}{4\pi}\boldsymbol{\nabla}\int_V\frac{\boldsymbol{\nabla}'\cdot\mathbf{F}}{\left|\mathbf{r}-\mathbf{r}'\right|}\mathrm{d}V'


so this is indeed the Helmholtz decomposition.[15]

See also[edit]

Notes[edit]

  1. ^ On Helmholtz's Theorem in Finite Regions. By Jean Bladel. Midwestern Universities Research Association, 1958.
  2. ^ Hermann von Helmholtz. Clarendon Press, 1906. By Leo Koenigsberger. p357
  3. ^ An Elementary Course in the Integral Calculus. By Daniel Alexander Murray. American Book Company, 1898. p8.
  4. ^ Vector Analysis. By Josiah Willard Gibbs. Yale University Press, 1901.
  5. ^ Electromagnetic theory, Volume 1. By Oliver Heaviside. "The Electrician" printing and publishing company, limited, 1893.
  6. ^ Elements of the differential calculus. By Wesley Stoker Barker Woolhouse. Weale, 1854.
  7. ^ An Elementary Treatise on the Integral Calculus: Founded on the Method of Rates Or Fluxions. By William Woolsey Johnson. John Wiley & Sons, 1881.
    See also: Method of Fluxions.
  8. ^ Vector Calculus: With Applications to Physics. By James Byrnie Shaw. D. Van Nostrand, 1922. p205.
    See also: Green's Theorem.
  9. ^ A Treatise on the Integral Calculus, Volume 2. By Joseph Edwards. Chelsea Publishing Company, 1922.
  10. ^ See:
    • H. Helmholtz (1858) "Über Integrale der hydrodynamischen Gleichungen, welcher der Wirbelbewegungen entsprechen" (On integrals of the hydrodynamic equations which correspond to vortex motions), Journal für die reine und angewandte Mathematik, 55: 25-55. On page 38, the components of the fluid's velocity (u, v, w) are expressed in terms of the gradient of a scalar potential P and the curl of a vector potential (L, M, N).
    • However, Helmholtz was largely anticipated by George Stokes in his paper: G. G. Stokes (presented: 1849 ; published: 1856) "On the dynamical theory of diffraction," Transactions of the Cambridge Philosophical Society, vol. 9, part I, pages 1-62; see pages 9-10.
  11. ^ "Helmholtz' Theorem". University of Vermont. 
  12. ^ a b c David J. Griffiths, Introduction to Electrodynamics, Prentice-Hall, 1999, p. 556.
  13. ^ Cantarella, Jason; DeTurck, Dennis; Gluck, Herman (2002). "Vector Calculus and the Topology of Domains in 3-Space". The American Mathematical Monthly 109 (5): 409–442. JSTOR 2695643. 
  14. ^ [0801.0335] Longitudinal and transverse components of a vector field
  15. ^ Online lecture notes by Robert Littlejohn

References[edit]

General references[edit]

  • George B. Arfken and Hans J. Weber, Mathematical Methods for Physicists, 4th edition, Academic Press: San Diego (1995) pp. 92–93
  • George B. Arfken and Hans J. Weber, Mathematical Methods for Physicists International Edition, 6th edition, Academic Press: San Diego (2005) pp. 95–101

References for the weak formulation[edit]

  • C. Amrouche, C. Bernardi, M. Dauge, and V. Girault. "Vector potentials in three dimensional non-smooth domains." Mathematical Methods in the Applied Sciences, 21, 823–864, 1998.
  • R. Dautray and J.-L. Lions. Spectral Theory and Applications, volume 3 of Mathematical Analysis and Numerical Methods for Science and Technology. Springer-Verlag, 1990.
  • V. Girault and P.A. Raviart. Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms. Springer Series in Computational Mathematics. Springer-Verlag, 1986.

External links[edit]