logo

Helmholtz Equation 📂Partial Differential Equations

Helmholtz Equation

Definition

The following partial differential equation is called the Helmholtz equation.

2u(x)+k2u(x)=Δu(x)+k2u(x)=(Δ+k2)u(x)=0,xRn \nabla^{2}u(x) + k^{2} u(x) = \Delta u(x) + k^{2} u(x) = (\Delta + k^{2} )u(x) = 0,\quad x \in \mathbb{R}^{n}

Here, 2=Δ\nabla ^{2} = \Delta is the Laplacian.

Explanation

It can also be expressed in the form of Δu=λu-\Delta u = \lambda u. Hence it is sometimes called the eigenvalue equation for the Laplace operator.

It can be derived from the wave equation, thus it is also referred to as the reduced wave equation1.

Unlike the wave equation, which includes derivatives with respect to both time and space, the Helmholtz equation lacks the time term, making it a partial differential equation that depends only on the spatial variables.

Derivation

The wave equation is as follows:

Δu(x,t)1c22u(x,t)t2=0 \Delta u(x,t) - \dfrac{1}{c^{2}}\dfrac{\partial^{2} u(x,t)}{\partial t^{2}} = 0

Here, cc represents the velocity of the wave.


Method 1

The solution to the wave equation, that is, the wave function, is as follows:

u(x,t)=u(x)u(t)=eikxeiωt=ei(kxωt) u (x,t) = u(x)u(t) = e^{ikx} e^{-i\omega t} = e^{i(kx - \omega t)}

Here, x,tx , t represent space and time, respectively, and k,ωk, \omega represent the wave number and angular frequency. When the speed of the wave is cc, the following relationship holds:

k=ωc k = \dfrac{\omega}{c}

Therefore, if we figure out utt(x,t)u_{tt}(x,t), it is as follows:

utt(x,t)=2t2ei(kxωt)=(iω)2ei(kxωt)=ω2ei(kxωt) u_{tt}(x,t) = \dfrac{\partial ^{2}}{\partial t^{2}}e^{i(kx - \omega t)} = (-i \omega)^{2}e^{i(kx - \omega t)} = -\omega^{2}e^{i(kx - \omega t)}

Substituting this into the wave equation yields the Helmholtz equation.

Δu1c22ut2= 0    Δei(kxωt)+ω2c2ei(kxωt)= 0    (Δeikx+k2eikx)eiωt= 0    Δeikx+k2eikx= 0    Δu(x)+k2u(x)= 0 \begin{align*} && \Delta u - \dfrac{1}{c^{2}}\dfrac{\partial^{2} u}{\partial t^{2}} =&\ 0 \\[1em] \implies && \Delta e^{i(kx - \omega t)} + \dfrac{\omega^{2}}{c^{2}} e^{i(kx - \omega t)} =&\ 0 \\[1em] \implies && \left( \Delta e^{ikx} + k^{2} e^{ikx} \right) e^{-i\omega t}=&\ 0 \\[1em] \implies && \Delta e^{ikx} + k^{2} e^{ikx}=&\ 0 \\[1em] \implies && \Delta u(x) + k^{2} u(x) =&\ 0 \end{align*}

Method 2

Taking the Fourier transform of the wave equation with respect to tt gives:

Δu(x,t)1c2utt(x,t)= 0    Δu^(x,ω)1c2utt^(x,ω)= 0 \begin{align*} && \Delta u(x,t) - \dfrac{1}{c^{2}}u_{tt}(x,t) =&\ 0 \\ \implies && \widehat{\Delta u}(x,\omega) - \dfrac{1}{c^{2}} \widehat{u_{tt}}(x,\omega) =&\ 0 \end{align*}

Here, using the property of Fourier transform u^(ω)=ω2u^(ω)\widehat{u^{\prime \prime}}(\omega) = - \omega^{2} \widehat{u}(\omega) for the second term yields:

Δu^(x,ω)+ω2c2u^(x,ω)= 0    Δu^(x,ω)+k2u^(x,ω)= 0    Δu(x,t)+k2u(x,t)= 0 \begin{align*} && \widehat{\Delta u}(x,\omega) + \dfrac{\omega^{2}}{c^{2}} \widehat{u}(x,\omega) =&\ 0 \\[1em] \implies && \widehat{\Delta u}(x,\omega) + k^{2} \widehat{u}(x,\omega) =&\ 0 \\[1em] \implies && \Delta u(x, t) + k^{2} u(x, t) =&\ 0 \end{align*}

Assuming that u(x,t)u(x,t) is separable in terms of variables,

Δu(x,t)+k2u(x,t)= 0    Δu(x)u(t)+k2u(x)u(t)= 0    Δu(x)+k2u(x)= 0 \begin{align*} && \Delta u(x, t) + k^{2} u(x, t) =&\ 0 \\[1em] \implies && \Delta u(x) u(t) + k^{2} u(x) u(t) =&\ 0 \\[1em] \implies && \Delta u(x) + k^{2} u(x) =&\ 0 \end{align*}

Method 3

Assuming that it is separable in variables as u(x,t)=u(x)v(t)u(x, t) = u(x)v(t), let’s organize the equation as follows:

2u(x,t)x2= 1c22u(x,t)t2    d2udx2v= 1c2d2vdt2u    1ud2udx2= 1c21vd2vdt2 \begin{align*} && \dfrac{\partial^{2} u(x, t)}{\partial x^{2}} =&\ \dfrac{1}{c^{2}}\dfrac{\partial^{2} u(x, t)}{\partial t^{2}} \\[1em] \implies && \dfrac{d^{2} u}{d x^{2}} v =&\ \dfrac{1}{c^{2}}\dfrac{d^{2} v }{d t^{2}} u \\[1em] \implies && \dfrac{1}{u}\dfrac{d^{2} u}{d x^{2}} =&\ \dfrac{1}{c^{2}} \dfrac{1}{v}\dfrac{d^{2} v }{d t^{2}} \end{align*}

Then, since the left-hand side is independent of tt and the right-hand side is independent of xx, it can be concluded that both sides are constants relative to xx and tt. Let’s call that constant k2-k^{2}. Then, we obtain:

1ud2udx2= k2    d2udx2= k2u    d2udx2+k2u= 0 \begin{align*} && \dfrac{1}{u}\dfrac{d^{2} u}{d x^{2}} =&\ -k^{2} \\[1em] \implies && \dfrac{d^{2} u}{d x^{2}} =&\ -k^{2}u \\[1em] \implies && \dfrac{d^{2} u}{d x^{2}} + k^{2}u =&\ 0 \end{align*}


  1. David Colton and Rainer Kress, Inverse Acoustic and Electromagnetic Scattering Theory (4th Edition, 2019), p15 ↩︎