Helmholtz Equation
Helmholtz Equation
Helmholtz Equation
Two sources of radiation in the plane, given mathematically by a function which is zero in the blue region.
The real part of the resulting field A, A is the solution to the inhomogeneous Helmholtz equation
The Helmholtz equation, named for Hermann von Helmholtz, is the elliptic partial differential equation
Contents
[hide]
1 Motivation and uses o 1.1 Harmonic solutions 2 Solving the Helmholtz equation using separation of variables o 2.1 Vibrating membrane o 2.2 Three-dimensional solutions 3 Paraxial approximation
Separation of variables begins by assuming that the wave function u(r, t) is in fact separable:
Substituting this form into the wave equation, and then simplifying, we obtain the following equation:
Notice the expression on the left-hand side depends only on r, whereas the right-hand expression depends only on t. As a result, this equation is valid in the general case if and only if both sides of the equation are equal to a constant value. From this observation, we obtain two equations, one for A(r), the other for T(t):
and
where we have chosen, without loss of generality, the expression k2 for the value of the constant. (It is equally valid to use any constant k as the separation constant; k2 is chosen only for convenience in the resulting solutions.)
for arbitrary (complex-valued) constants C and D, which will depend on the initial conditions and boundary conditions, and subject to the dispersion relation
We now have Helmholtz's equation for the spatial variable r and a second-order ordinary differential equation in time. The solution in time will be a linear combination of sine and cosine functions, with angular frequency of , while the form of the solution in space will depend on the boundary conditions. Alternatively, integral transforms, such as the Laplace or Fourier transform, are often used to transform a hyperbolic PDE into a form of the Helmholtz equation. Because of its relationship to the wave equation, the Helmholtz equation arises in problems in such areas of physics as the study of electromagnetic radiation, seismology, and acoustics.
and
and that n must be an integer. The radial component R has the form
where the Bessel function Jn() satisfies Bessel's equation 2Jn'' + Jn' + (2 n2)Jn = 0, and =kr. The radial function Jn has infinitely many roots for each value of n, denoted by m,n. The boundary condition that A vanishes where r=a will be satisfied if the corresponding wavenumbers are given by
The general solution A then takes the form of a doubly infinite sum of terms involving products of
This solution arises from the spatial solution of the wave equation and diffusion equation. Here and are the spherical Bessel functions, and
are the spherical harmonics (Abramowitz and Stegun, 1964). Note that these forms are general solutions, and require boundary conditions to be specified to be used in any specific case. For infinite exterior domains, a radiation condition may also be required (Sommerfeld, 1949). For function A(r0) has asymptotics
where function is called scattering amplitude and u0(r0) is the value of A at each boundary point r0.
where
This equation has important applications in the science of optics, where it provides solutions that describe the propagation of electromagnetic waves (light) in the form of either paraboloidal waves or Gaussian beams. Most lasers emit beams that take this form. In the paraxial approximation, the complex magnitude of the electric field E becomes
where A represents the complex-valued amplitude of the electric field, which modulates the sinusoidal plane wave represented by the exponential factor. The paraxial approximation places certain upper limits on the variation of the amplitude function A with respect to longitudinal distance z. Specifically:
and
These conditions are equivalent to saying that the angle between the wave vector k and the optical axis z must be small enough so that
The paraxial form of the Helmholtz equation is found by substituting the above-stated complex magnitude of the electric field into the general form of the Helmholtz equation as follows.
Because of the paraxial inequalities stated above, the 2A/z2 factor is neglected in comparison with the A/z factor. The yields the Paraxial Helmholtz equation.
where : Rn C is a given function with compact support, and n = 1, 2, 3. This equation is very similar to the screened Poisson equation, and would be identical if the plus sign (in front of the k term) is switched to a minus sign. In order to solve this equation uniquely, one needs to specify a boundary condition at infinity, which is typically the Sommerfeld radiation condition
uniformly in
with
With this condition, the solution to the inhomogeneous Helmholtz equation is the convolution
(notice this integral is actually over a finite region, since f has compact support). Here, G is the Green's function of this equation, that is, the solution to the inhomogeneous Helmholtz equation with equaling the Dirac delta function, so G satisfies
The expression for the Green's function depends on the dimension n of the space. One has
for n = 1,
for n = 2, where
for n = 3. Note that we have chosen the boundary condition that the Green's function is an outgoing wave for .