Jump to content

Helmholtz equation

From Wikipedia, the free encyclopedia

This is an old revision of this page, as edited by Army1987 (talk | contribs) at 20:56, 27 January 2009 (Paraxial approximation: formatting). The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

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

where ∇2 is the Laplacian, k is a constant, and the unknown function A = A(r) is defined on n-dimensional Euclidean space Rn (typically n = 1, 2, or 3, when the solution to this equation makes physical sense[discuss]).

Motivation and uses

The Helmholtz equation often arises in the study of physical problems involving partial differential equations (PDEs) in both space and time. The Helmholtz equation, which represents the time-independent form of the original equation, results from applying the technique of separation of variables to reduce the complexity of the analysis.

For example, consider the wave equation

Separation of variables begins by assuming that the wave function u(rt) 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.)

Rearranging the first equation, we obtain the Helmholtz equation:

Likewise, after making the substitution

the second equation becomes

where k is the wave vector and ω is the angular frequency.

Harmonic solutions

It is relatively easy to show that solutions to the Helmholtz equation will take the form:

which corresponds to the time-harmonic solution

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.

Solving the Helmholtz equation using separation of variables

The general solution to the spatial Helmholtz equation

can be obtained using separation of variables.

Vibrating membrane

The two-dimensional analogue of the vibrating string is the vibrating membrane, with the edges clamped to be motionless. The Helmholtz equation was solved for many basic shapes in the 19th century: the rectangular membrane by Siméon Denis Poisson in 1829, the equilateral triangle by Gabriel Lamé in 1852, and the circular membrane by Alfred Clebsch in 1862. The elliptical drumhead was studied by Emile Mathieu, leading to Mathieu's differential equation. The solvable shapes all correspond to shapes whose dynamical billiard table is integrable, that is, not chaotic. When the motion on a correspondingly-shaped billiard table is chaotic, then no closed form solutions to the Helmholtz equation are known. The study of such systems is known as quantum chaos, as the Helmholtz equation and similar equations occur in quantum mechanics.

If the edges of a shape are straight line segments, then a solution is integrable or knowable in closed-form only if it is expressible as a finite linear combination of plane waves that satisfy the boundary conditions (zero at the boundary, i.e., membrane clamped).

An interesting situation happens with a shape where about half of the solutions are integrable, but the remainder are not. A simple shape where this happens is with the regular hexagon. If the wavepacket describing a quantum billiard ball is made up of only the closed-form solutions, its motion will not be chaotic, but if any amount of non-closed-form solutions are included, the quantum billiard motion becomes chaotic. Another simple shape where this happens is with an "L" shape made by reflecting a square down, then to the right.

If the domain is a circle of radius a, then it is appropriate to introduce polar coordinates r and θ. The Helmholtz equation takes the form

We may impose the boundary condition that A vanish if r=a; thus

The method of separation of variables leads to trial solutions of the form

where Θ must be periodic of period 2π. This leads to

and

It follows from the periodicity condition that

and that n must be an integer. The radial component R has the form

where the Bessel function Jn(ρ) satisfies Bessel's equation

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 frequencies are given by

The general solution A then takes the form of a doubly infinite sum of terms involving products of

These solutions are the modes of vibration of a circular drumhead.

Three-dimensional solutions

In spherical coordinates, the solution is:

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 has asymptotics

when

where function ƒ is called scattering amplitude and is the value of A at each boundary point .

Paraxial approximation

The paraxial approximation of the Helmholtz equation is:

where is the transverse part of the Laplacian.

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.

Expansion and cancellation yields the following:

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.

Inhomogeneous Helmholtz equation

The inhomogeneous Helmholtz equation is the equation

where is a given function with compact support, and This equation is very similar to the screened Poisson equation, and would be identical if the plus sign 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 , where the vertical bars denote the Euclidean norm.

With this condition, the solution to the inhomogeneous Helmholtz equation is the convolution

(notice this integral is actually over a finite region, since has compact support). Here, 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 of the space. One has

for n = 1,

for n = 2, where is a Hankel function, and

for n = 3.

References

  • M. Abramowitz and I. Stegun eds., Handbook of Mathematical functions with Formulas, Graphs and Mathematical Tables, National Bureau of Standards. Washington, D. C., 1964.
  • Riley, K.F., Hobson, M.P., and Bence, S.J. (2002). Mathematical methods for physics and engineering, Cambridge University Press, ch. 19. ISBN 0-521-89067-5.
  • McQuarrie, Donald A. (2003). Mathematical Methods for Scientists and Engineers, University Science Books: Sausalito, California, Ch. 16. ISBN 1-891389-24-6.
  • Bahaa E. A. Saleh and Malvin Carl Teich (1991). Fundamentals of Photonics. New York: John Wiley & Sons. ISBN 0-471-83965-5. Chapter 3, "Beam Optics," pp. 80–107.
  • A. Sommerfeld, Partial Differential Equations in Physics, Academic Press, New York, New York, 1949.
  • Howe, M. S. (1998). Acoustics of fluid-structure interactions. Cambridge; New York: Cambridge University Press. ISBN 0-521-63320-6.