Skip to content

The Wave Equation | Physics - Wyatt's Notes

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 t So:

×(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}

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.

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.

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) = 0 So:

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) = 0 So:

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

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.

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

Mistake 1: Confusing phase velocity with group velocity Phase velocity vp=ω/kv_p = \omega/k is the speed of individual wave crests, while group velocity vg=dω/dkv_g = d\omega/dk is the speed of the wave packet envelope (and hence of information). In non-dispersive media they are equal, but in dispersive media they differ. A common error is assuming energy travels at the phase velocity, which is incorrect when dispersion is present.

Mistake 2: Forgetting that standing waves do not transport energy A standing wave is a superposition of two counter-propagating waves of equal amplitude. The net energy flux is zero because the energy flowing in one direction is exactly cancelled by energy flowing in the other. Students sometimes assume standing waves carry energy because individual travelling waves do. The energy in a standing wave oscillates between kinetic and potential forms but does not propagate.

Mistake 3: Misapplying superposition to non-linear media The superposition principle applies only to linear media. In non-linear optics (e.g., high-intensity lasers in certain crystals), waves interact and the total field is not directly the sum of individual fields. The wave equation itself is derived from Maxwell’s equations, which are linear in vacuum, but material responses can introduce non-linearity.

  • Maxwell’s Equations: The electromagnetic wave equation is derived from Maxwell’s equations in free space by taking curls of Faraday’s and Ampere-Maxwell’s laws.
  • Electrodynamics: Electrodynamics provides the foundation for understanding how changing electric and magnetic fields propagate as waves.
  • Interference: The superposition of wave solutions to the wave equation produces interference patterns used in Young’s double slit and Fabry-Perot interferometers.
  • Calculus
  • Linear Algebra
  • Vector Calculus
  • Quantum Computing