Skip to content

The Wave Equation

1.1 Derivation from Maxwell’s Equations

In vacuum, with no sources (ρ=0\rho = 0, J=0\mathbf{J} = \mathbf{0}), Maxwell’s equations give:

×(×E)=(E)2E=2E\nabla \times (\nabla \times \mathbf{E}) = \nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E} = -\nabla^2 \mathbf{E}

From Faraday’s law: ×E=B/t\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial tSo:

×(Bt)=t(×B)=μ0ε02Et2\nabla \times \left(-\frac{\partial \mathbf{B}}{\partial t}\right) = -\frac{\partial}{\partial t}(\nabla \times \mathbf{B}) = -\mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}

Hence:

2E=μ0ε02Et2\nabla^2 \mathbf{E} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}

This is the electromagnetic wave equation with wave speed c=1/μ0ε0c = 1/\sqrt{\mu_0 \varepsilon_0}.

The same equation holds for B\mathbf{B}:

2B=μ0ε02Bt2\nabla^2 \mathbf{B} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{B}}{\partial t^2}

1.2 General Solutions

The one-dimensional wave equation:

2ux2=1v22ut2\frac{\partial^2 u}{\partial x^2} = \frac{1}{v^2}\frac{\partial^2 u}{\partial t^2}

Has the general solution (d’Alembert’s solution):

u(x,t)=f(xvt)+g(x+vt)u(x, t) = f(x - vt) + g(x + vt)

Where ff is a wave travelling in the +x+x direction and gg in the x-x direction.

Proof. Substitute u=f(xvt)u = f(x - vt). Let ξ=xvt\xi = x - vt. Then u/x=f(ξ)\partial u/\partial x = f'(\xi) 2u/x2=f(ξ)\partial^2 u/\partial x^2 = f''(\xi), u/t=vf(ξ)\partial u/\partial t = -vf'(\xi) 2u/t2=v2f(ξ)\partial^2 u/\partial t^2 = v^2 f''(\xi). The wave equation gives f=(v2/v2)ff'' = (v^2/v^2)f''Which Is identically satisfied. The same holds for g(x+vt)g(x + vt). By linearity, the sum is also a solution. \blacksquare

Worked Example: Verifying a wave solution

Problem. A string under tension has the wave equation 2y/t2=1002y/x2\partial^2 y/\partial t^2 = 100\,\partial^2 y/\partial x^2. Verify that y(x,t)=0.03sin(5x50t)y(x,t) = 0.03\sin(5x - 50t) is a solution. Find the wave speed, wavelength, frequency, and propagation direction.

Solution. The wave speed is v=100=10v = \sqrt{100} = 10 m/s. Comparing with the standard form y=Asin(kxωt)y = A\sin(kx - \omega t): A=0.03A = 0.03 m, k=5k = 5 rad/m, ω=50\omega = 50 rad/s.

Verification: 2y/x2=25×0.03sin(5x50t)\partial^2 y/\partial x^2 = -25 \times 0.03\sin(5x - 50t) and 2y/t2=2500×0.03sin(5x50t)\partial^2 y/\partial t^2 = -2500 \times 0.03\sin(5x - 50t)So 2y/t2=100×2y/x2\partial^2 y/\partial t^2 = 100 \times \partial^2 y/\partial x^2. \checkmark

Wavelength: λ=2π/k=2π/51.26\lambda = 2\pi/k = 2\pi/5 \approx 1.26 m. Frequency: f=ω/(2π)=50/(2π)7.96f = \omega/(2\pi) = 50/(2\pi) \approx 7.96 Hz. Direction: the argument 5x50t5x - 50t indicates propagation in the +x+x direction.

1.3 Complex Representation

It is convenient to write monochromatic waves as:

E(r,t)=Re[E~ei(krωt)]\mathbf{E}(\mathbf{r}, t) = \mathrm{Re}\left[\tilde{\mathbf{E}}\, e^{i(\mathbf{k}\cdot\mathbf{r} - \omega t)}\right]

Where E~\tilde{\mathbf{E}} is the complex amplitude, k\mathbf{k} is the wave vector, and ω\omega is The angular frequency. The dispersion relation is ω=ck=ck\omega = ck = c|\mathbf{k}|.

The wave vector satisfies k=2π/λ|\mathbf{k}| = 2\pi/\lambda and ω=2πν\omega = 2\pi\nu.

When computing intensities, the complex representation simplifies the algebra. For a plane wave with Complex amplitude E~\tilde{E}:

I=12cε0E~2I = \frac{1}{2}c\varepsilon_0 |\tilde{E}|^2

Worked Example: Complex amplitude and intensity

Problem. A plane wave propagating in the +z+z direction has complex amplitude E~=(10+5i)x^\tilde{\mathbf{E}} = (10 + 5i)\hat{\mathbf{x}} V/m in vacuum. Find the real electric field, the Intensity, and the phase of the wave relative to cos(kzωt)\cos(kz - \omega t).

Solution. E~=102+52=12511.18|\tilde{E}| = \sqrt{10^2 + 5^2} = \sqrt{125} \approx 11.18 V/m. Phase: ϕ=arctan(5/10)=26.57°\phi = \arctan(5/10) = 26.57°.

Real field: E(z,t)=11.18cos(kzωt+26.57°)x^\mathbf{E}(z,t) = 11.18\cos(kz - \omega t + 26.57°)\,\hat{\mathbf{x}} V/m.

Intensity: I=12cε0E~2=12(3×108)(8.854×1012)(125)=0.166I = \frac{1}{2}c\varepsilon_0 |\tilde{E}|^2 = \frac{1}{2}(3 \times 10^8)(8.854 \times 10^{-12})(125) = 0.166 W/m2^2.

1.4 Standing Waves

When two waves of equal amplitude and frequency travel in opposite directions, their superposition Produces a standing wave. Consider:

u1=Asin(kxωt),u2=Asin(kx+ωt)u_1 = A\sin(kx - \omega t), \quad u_2 = A\sin(kx + \omega t)

Using the identity sinα+sinβ=2sinα+β2cosαβ2\sin\alpha + \sin\beta = 2\sin\frac{\alpha+\beta}{2}\cos\frac{\alpha-\beta}{2}:

u(x,t)=2Asin(kx)cos(ωt)u(x,t) = 2A\sin(kx)\cos(\omega t)

This is a standing wave with the following properties:

  • Nodes (points of zero displacement): kx=mπkx = m\piI.e., x=mλ/2x = m\lambda/2 for m=0,1,2,m = 0, 1, 2, \ldots
  • Antinodes (points of maximum displacement): kx=(m+1/2)πkx = (m + 1/2)\piI.e., x=(2m+1)λ/4x = (2m+1)\lambda/4.
  • All points between two nodes oscillate in phase (or in antiphase with adjacent segments).
  • The standing wave does not transport energy in either direction.

Standing waves on a string of length LL fixed at both ends. The boundary conditions u(0,t)=0u(0,t) = 0 And u(L,t)=0u(L,t) = 0 require sin(kL)=0\sin(kL) = 0So:

knL=nπ    λn=2Ln,fn=nv2L,n=1,2,3,k_n L = n\pi \implies \lambda_n = \frac{2L}{n}, \quad f_n = \frac{nv}{2L}, \quad n = 1, 2, 3, \ldots

The allowed frequencies are integer multiples of the fundamental frequency f1=v/(2L)f_1 = v/(2L). The general solution is a superposition of all normal modes:

u(x,t)=n=1(Ansinknxcosωnt+Bnsinknxsinωnt)u(x,t) = \sum_{n=1}^{\infty} \left(A_n \sin k_n x \cos \omega_n t + B_n \sin k_n x \sin \omega_n t\right)

Standing waves on a string fixed at one end (x=0x = 0) and free at the other (x=Lx = L). The free end requires u/xx=L=0\partial u/\partial x|_{x=L} = 0Giving cos(kL)=0\cos(kL) = 0So:

knL=(n+1/2)π    fn=(2n+1)v4L,n=0,1,2,k_n L = (n + 1/2)\pi \implies f_n = \frac{(2n+1)v}{4L}, \quad n = 0, 1, 2, \ldots

Only odd harmonics are present. A pipe open at one end and closed at the other behaves analogously For sound waves.

Worked Example: Guitar string harmonics

Problem. A guitar string of length L=0.65L = 0.65 m is fixed at both ends. The fundamental frequency Is f1=330f_1 = 330 Hz (the note E4). Find the wave speed, the frequency of the third harmonic, and the Positions of the nodes for the third harmonic.

Solution. Wave speed: v=2Lf1=2×0.65×330=429v = 2Lf_1 = 2 \times 0.65 \times 330 = 429 m/s.

Third harmonic: f3=3f1=990f_3 = 3f_1 = 990 Hz. Wavelength: λ3=2L/3=0.433\lambda_3 = 2L/3 = 0.433 m.

Nodes at x=mλ3/2=m(0.433/2)=0.217mx = m\lambda_3/2 = m(0.433/2) = 0.217m for m=0,1,2,3m = 0, 1, 2, 3. Positions: x=0,0.217,0.433,0.650x = 0, 0.217, 0.433, 0.650 m (the endpoints and two interior nodes).

1.5 Energy Transport by Waves

For a transverse wave on a string of linear mass density μ\mu under tension TTWith v=T/μv = \sqrt{T/\mu}:

Kinetic energy of an element dxdx:

dK=12μdx(ut)2dK = \frac{1}{2}\mu\,dx\left(\frac{\partial u}{\partial t}\right)^2

Potential energy (from stretching the string):

dU=12Tdx(ux)2dU = \frac{1}{2}T\,dx\left(\frac{\partial u}{\partial x}\right)^2

Proof that kinetic and potential energies are equal. For a rightward-travelling wave u=f(xvt)u = f(x - vt):

ut=vf(xvt)=vux\frac{\partial u}{\partial t} = -vf'(x - vt) = -v\frac{\partial u}{\partial x}

Since v2=T/μv^2 = T/\mu (i.e., T=μv2T = \mu v^2):

dU=12μv2(ux)2=12μ(ut)2=dKdU = \frac{1}{2}\mu v^2\left(\frac{\partial u}{\partial x}\right)^2 = \frac{1}{2}\mu\left(\frac{\partial u}{\partial t}\right)^2 = dK \quad \blacksquare

The total energy density (energy per unit length) is:

dEdx=dK+dU=μ(ut)2\frac{dE}{dx} = dK + dU = \mu\left(\frac{\partial u}{\partial t}\right)^2

Energy flux (power): The rate of energy transport past a point is:

P=Tuxut=μv(ut)2P = -T\frac{\partial u}{\partial x}\frac{\partial u}{\partial t} = \mu v\left(\frac{\partial u}{\partial t}\right)^2

For a sinusoidal wave u=Asin(kxωt)u = A\sin(kx - \omega t)The time-averaged power is:

P=12μvω2A2\langle P \rangle = \frac{1}{2}\mu v \omega^2 A^2

The intensity (power per unit area, generalising to 3D) is:

I=12ρvω2A2I = \frac{1}{2}\rho v \omega^2 A^2

Where ρ\rho is the mass density of the medium.

Worked Example: Wave energy and power on a string

Problem. A steel wire of diameter 1.0 mm and density 78007800 kg/m3^3 is under 500 N of tension. A sinusoidal wave of amplitude 5.0 mm and frequency 200 Hz propagates along it. Find (a) the wave Speed, (b) the average power, and (c) the intensity (power per unit cross-sectional area).

Solution. Cross-sectional area: Awire=π(0.5×103)2=7.85×107A_{\mathrm{wire} = \pi(0.5 \times 10^{-3})^2 = 7.85 \times 10^{-7}} m2^2. Linear density: μ=ρAwire=7800×7.85×107=6.12×103\mu = \rho A_{\mathrm{wire} = 7800 \times 7.85 \times 10^{-7} = 6.12 \times 10^{-3}} kg/m.

(a) v=T/μ=500/(6.12×103)=8.17×104=286v = \sqrt{T/\mu} = \sqrt{500/(6.12 \times 10^{-3})} = \sqrt{8.17 \times 10^4} = 286 m/s.

(b) ω=2π×200=1257\omega = 2\pi \times 200 = 1257 rad/s. P=12μvω2A2=12(6.12×103)(286)(1257)2(5.0×103)2\langle P \rangle = \frac{1}{2}\mu v \omega^2 A^2 = \frac{1}{2}(6.12 \times 10^{-3})(286)(1257)^2(5.0 \times 10^{-3})^2 =12(6.12×103)(286)(1.58×106)(2.5×105)=34.7= \frac{1}{2}(6.12 \times 10^{-3})(286)(1.58 \times 10^6)(2.5 \times 10^{-5}) = 34.7 W.

(c) Intensity: I=P/Awire=34.7/(7.85×107)=4.42×107I = \langle P \rangle/A_{\mathrm{wire} = 34.7/(7.85 \times 10^{-7}) = 4.42 \times 10^7} W/m2^2.

1.6 Wave Packets and Group Velocity

A real wave is never perfectly monochromatic. A wave packet is a superposition of plane waves With a narrow range of frequencies and wave vectors:

ψ(x,t)=A(k)ei(kxω(k)t)dk\psi(x,t) = \int_{-\infty}^{\infty} A(k)\, e^{i(kx - \omega(k)t)}\,dk

Where A(k)A(k) is the spectral amplitude, peaked around k0k_0 with width Δk\Delta k.

Expanding the dispersion relation around k0k_0:

ω(k)ω0+vg(kk0)+12α(kk0)2\omega(k) \approx \omega_0 + v_g(k - k_0) + \frac{1}{2}\alpha(k - k_0)^2

Where vg=dω/dkk0v_g = \left.d\omega/dk\right|_{k_0} is the group velocity and α=d2ω/dk2\alpha = d^2\omega/dk^2 is the group velocity dispersion (GVD).

Substituting and carrying out the Gaussian integral (for a Gaussian envelope A(k)A(k)):

ψ(x,t)2exp((xvgt)22σx2(t))|\psi(x,t)|^2 \propto \exp\left(-\frac{(x - v_g t)^2}{2\sigma_x^2(t)}\right)

Where σx(t)=σx(0)1+(αt/2σx2(0))2\sigma_x(t) = \sigma_x(0)\sqrt{1 + (\alpha t / 2\sigma_x^2(0))^2}.

The envelope moves at vgv_gWhile individual wave crests move at the phase velocity vp=ω/kv_p = \omega/k. The packet broadens over time due to GVD.

  • In a non-dispersive medium (ωk\omega \propto kSo vg=vpv_g = v_p): the packet propagates without distortion.
  • In a normally dispersive medium (d2ω/dk2>0d^2\omega/dk^2 \gt 0): vg<vpv_g \lt v_p and the packet broadens.

Relation to the wave equation. The 1D wave equation 2u/t2=v22u/x2\partial^2 u/\partial t^2 = v^2\partial^2 u/\partial x^2 Has dispersion relation ω=±vk\omega = \pm vkGiving vg=vp=vv_g = v_p = v — it is non-dispersive.

Worked Example: Group velocity in a dispersive medium

Problem. A dispersive medium has the dispersion relation ω=αk2\omega = \alpha k^2 with α=5.0×103\alpha = 5.0 \times 10^3 m2^2/s. For a wave packet centred at k0=200k_0 = 200 rad/m, Find (a) the phase velocity, (b) the group velocity, and (c) the time for the packet width to Double, given an initial width σ0=0.01\sigma_0 = 0.01 m.

Solution.

(a) Phase velocity: vp=ω/k=αk0=5.0×103×200=1.0×106v_p = \omega/k = \alpha k_0 = 5.0 \times 10^3 \times 200 = 1.0 \times 10^6 m/s.

(b) Group velocity: vg=dω/dk=2αk0=2×5.0×103×200=2.0×106v_g = d\omega/dk = 2\alpha k_0 = 2 \times 5.0 \times 10^3 \times 200 = 2.0 \times 10^6 m/s.

(c) GVD: αGVD=d2ω/dk2=2α=1.0×104\alpha_{\mathrm{GVD} = d^2\omega/dk^2 = 2\alpha = 1.0 \times 10^4} m2^2/s. The packet width doubles when 1+(αGVDt/2σ02)2=41 + (\alpha_{\mathrm{GVD}\, t / 2\sigma_0^2)^2 = 4}: αGVDt/2σ02=3\alpha_{\mathrm{GVD}\, t / 2\sigma_0^2 = \sqrt{3}} t=2σ023/αGVD=2(104)(1.732)/(104)=3.46×108t = 2\sigma_0^2\sqrt{3}/\alpha_{\mathrm{GVD} = 2(10^{-4})(1.732)/(10^4) = 3.46 \times 10^{-8}} s =34.6= 34.6 ns.

:::caution Common Pitfall The group velocity vg=dω/dkv_g = d\omega/dk describes the motion of the wave packet envelope and equals The energy transport velocity in lossless media. However, it is only equal to the signal velocity In regions of weak, normal dispersion. Near absorption resonances, vgv_g can exceed cc or become Negative — this does not violate causality, since the true signal velocity (front velocity) never Exceeds cc.

:::