Helmholtz equation

Last updated

In mathematics, the Helmholtz equation is the eigenvalue problem for the Laplace operator. It corresponds to the linear partial differential equation

Contents

where 2 is the Laplace operator, k2 is the eigenvalue, and f is the (eigen)function. When the equation is applied to waves, k is known as the wave number. The Helmholtz equation has a variety of applications in physics and other sciences, including the wave equation, the diffusion equation, and the Schrödinger equation for a free particle.

In optics, the Helmholtz equation is the wave equation for the electric field. [1]

The equation is named after Hermann von Helmholtz, who studied it in 1860. [2]

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 a time-independent form of the wave 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(r, t) is in fact separable:

Substituting this form into the wave equation and then simplifying, we obtain the following equation:

Notice that the expression on the left side depends only on r, whereas the right 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 the same constant value. This argument is key in the technique of solving linear partial differential equations by separation of variables. From this observation, we obtain two equations, one for A(r), the other for T(t):

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 ω = kc, where k is the wave number, and ω is the angular frequency (assuming a monochromatic field), the second equation becomes

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, whose exact form is determined by initial conditions, 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 solution to the spatial Helmholtz equation:

can be obtained for simple geometries 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 Émile Mathieu, leading to Mathieu's differential equation.

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).

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 vanishes 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

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

The general solution A then takes the form of a generalized Fourier series of terms involving products of Jn(km,nr) and the sine (or cosine) 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 j(kr) and y(kr) are the spherical Bessel functions, and Ym
(θ, φ)
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).

Writing r0 = (x, y, z) function A(r0) has asymptotics

where function f is called scattering amplitude and u0(r0) is the value of A at each boundary point r0.

Three-dimensional solutions given the function on a 2-dimensional plane

Given a 2-dimensional plane where A is known, the solution to the Helmholtz equation is given by: [3]

where

  • is the solution at the 2-dimensional plane,

As z approaches zero, all contributions from the integral vanish except for r=0. Thus up to a numerical factor, which can be verified to be 1 by transforming the integral to polar coordinates .

This solution is important in diffraction theory, e.g. in deriving Fresnel diffraction.

Paraxial approximation

In the paraxial approximation of the Helmholtz equation, [4] the complex amplitude A is expressed as

where u represents the complex-valued amplitude which modulates the sinusoidal plane wave represented by the exponential factor. Then under a suitable assumption, u approximately solves

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.

The assumption under which the paraxial approximation is valid is that the z derivative of the amplitude function u is a slowly varying function of z:

This condition is equivalent to saying that the angle θ between the wave vector k and the optical axis z is small: θ ≪ 1.

The paraxial form of the Helmholtz equation is found by substituting the above-stated expression for the complex amplitude into the general form of the Helmholtz equation as follows:

Expansion and cancellation yields the following:

Because of the paraxial inequality stated above, the 2u/∂z2 term is neglected in comparison with the k·∂u/∂z term. This yields the paraxial Helmholtz equation. Substituting u(r) = A(r) eikz then gives the paraxial equation for the original complex amplitude A:

The Fresnel diffraction integral is an exact solution to the paraxial Helmholtz equation. [5]

Inhomogeneous Helmholtz equation

Helmholtz source.png
Two sources of radiation in the plane, given mathematically by a function f, which is zero in the blue region
Helmholtz solution.png
The real part of the resulting field A, A is the solution to the inhomogeneous Helmholtz equation (∇2 + k2) A = −f.

The inhomogeneous Helmholtz equation is the equation

where ƒ : RnC is a 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) were 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

in spatial dimensions, for all angles (i.e. any value of ). Here where are the coordinates of the vector .

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

(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 f 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 H(1)
0
is a Hankel function, and

for n = 3. Note that we have chosen the boundary condition that the Green's function is an outgoing wave for |x| → ∞.

Finally, for general n,

where and . [6]

See also

Notes

  1. Blanche, Pierre-Alexandre (2014). Field guide to holography. SPIE field guides. Bellingham, Wash: SPIE Press. ISBN   978-0-8194-9957-8.
  2. Helmholtz Equation, from the Encyclopedia of Mathematics.
  3. Mehrabkhani, S., & Schneider, T. (2017). Is the Rayleigh-Sommerfeld diffraction always an exact reference for high speed diffraction algorithms?. Optics express, 25(24), 30229-30240.
  4. J. W. Goodman. Introduction to Fourier Optics (2nd ed.). pp. 61–62.
  5. Grella, R. (1982). "Fresnel propagation and diffraction and paraxial wave equation". Journal of Optics. 13 (6): 367–374. Bibcode:1982JOpt...13..367G. doi:10.1088/0150-536X/13/6/006.
  6. Björn Engquist; Hongkai Zhao (November 2018). "Approximate Separability of the Green's Function of the Helmholtz Equation in the High Frequency Limit". Communications on Pure and Applied Mathematics . 71 (11): 2220–2274. doi:10.1002/cpa.21755.

Related Research Articles

<span class="mw-page-title-main">Spherical coordinate system</span> 3-dimensional coordinate system

In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a given point in space is specified by three numbers, : the radial distance of the radial liner connecting the point to the fixed point of origin ; the polar angle θ of the radial line r; and the azimuthal angle φ of the radial line r.

<span class="mw-page-title-main">Wave equation</span> Differential equation important in physics

The wave equation is a second-order linear partial differential equation for the description of waves or standing wave fields such as mechanical waves or electromagnetic waves. It arises in fields like acoustics, electromagnetism, and fluid dynamics.

<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

<span class="mw-page-title-main">Gaussian beam</span> Monochrome light beam whose amplitude envelope is a Gaussian function

In optics, a Gaussian beam is a beam of electromagnetic radiation with high monochromaticity whose amplitude envelope in the transverse plane is given by a Gaussian function; this also implies a Gaussian intensity (irradiance) profile. This fundamental (or TEM00) transverse Gaussian mode describes the intended output of most (but not all) lasers, as such a beam can be focused into the most concentrated spot. When such a beam is refocused by a lens, the transverse phase dependence is altered; this results in a different Gaussian beam. The electric and magnetic field amplitude profiles along any such circular Gaussian beam (for a given wavelength and polarization) are determined by a single parameter: the so-called waist w0. At any position z relative to the waist (focus) along a beam having a specified w0, the field amplitudes and phases are thereby determined as detailed below.

<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).

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 continuum mechanics, the infinitesimal strain theory is a mathematical approach to the description of the deformation of a solid body in which the displacements of the material particles are assumed to be much smaller than any relevant dimension of the body; so that its geometry and the constitutive properties of the material at each point of space can be assumed to be unchanged by the deformation.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

Fourier optics is the study of classical optics using Fourier transforms (FTs), in which the waveform being considered is regarded as made up of a combination, or superposition, of plane waves. It has some parallels to the Huygens–Fresnel principle, in which the wavefront is regarded as being made up of a combination of spherical wavefronts whose sum is the wavefront being studied. A key difference is that Fourier optics considers the plane waves to be natural modes of the propagation medium, as opposed to Huygens–Fresnel, where the spherical waves originate in the physical medium.

Geometrical optics, or ray optics, is a model of optics that describes light propagation in terms of rays. The ray in geometrical optics is an abstraction useful for approximating the paths along which light propagates under certain circumstances.

This is a list of some vector calculus formulae for working with common curvilinear coordinate systems.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

In quantum physics, the scattering amplitude is the probability amplitude of the outgoing spherical wave relative to the incoming plane wave in a stationary-state scattering process. At large distances from the centrally symmetric scattering center, the plane wave is described by the wavefunction

<span class="mw-page-title-main">Stokes flow</span> Type of fluid flow

Stokes flow, also named creeping flow or creeping motion, is a type of fluid flow where advective inertial forces are small compared with viscous forces. The Reynolds number is low, i.e. . This is a typical situation in flows where the fluid velocities are very slow, the viscosities are very large, or the length-scales of the flow are very small. Creeping flow was first studied to understand lubrication. In nature, this type of flow occurs in the swimming of microorganisms and sperm. In technology, it occurs in paint, MEMS devices, and in the flow of viscous polymers generally.

<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.

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 fluid dynamics, the Oseen equations describe the flow of a viscous and incompressible fluid at small Reynolds numbers, as formulated by Carl Wilhelm Oseen in 1910. Oseen flow is an improved description of these flows, as compared to Stokes flow, with the (partial) inclusion of convective acceleration.

<span class="mw-page-title-main">Lagrangian mechanics</span> Formulation of classical mechanics

In physics, Lagrangian mechanics is a formulation of classical mechanics founded on the stationary-action principle. It was introduced by the Italian-French mathematician and astronomer Joseph-Louis Lagrange in his presentation to the Turin Academy of Science in 1760 culminating in his 1788 grand opus, Mécanique analytique.

<span class="mw-page-title-main">Potential flow around a circular cylinder</span> Classical solution for inviscid, incompressible flow around a cyclinder

In mathematics, potential flow around a circular cylinder is a classical solution for the flow of an inviscid, incompressible fluid around a cylinder that is transverse to the flow. Far from the cylinder, the flow is unidirectional and uniform. The flow has no vorticity and thus the velocity field is irrotational and can be modeled as a potential flow. Unlike a real fluid, this solution indicates a net zero drag on the body, a result known as d'Alembert's paradox.

Partial-wave analysis, in the context of quantum mechanics, refers to a technique for solving scattering problems by decomposing each wave into its constituent angular-momentum components and solving using boundary conditions.

References