ACOUSTICS Wave Equation
ACOUSTICS Wave Equation
Department of Physics
1 / 64
Chapter Two Roadmap
2 / 64
Chapter Objectives (Merged)
3 / 64
Microscopic Picture of Sound Propagation (Air Element)
4 / 64
Model: A Small Air “Box” (3D Picture)
∆y
∆x
5 / 64
Pressure Gradient & Net Force (1D Sketch)
air box
6 / 64
Step 1: Equation of Motion (Vector Form)
f⃗ = −(∇p) V
d⃗v d⃗v
−(∇p) V = ρ0 V ⇒ ρ0 = −∇p
dt dt
Meaning
A spatial pressure gradient accelerates fluid particles.
7 / 64
1D Form of the Motion Equation
Physical Meaning
Pressure gradient ⇒ acceleration. If pressure rises to the right, particles accelerate to the left.
8 / 64
Equation of Motion: Visual Snapshot
p(x, t) snapshot
x
du
∂p/∂x > 0 ⇒ dt
<0
9 / 64
Why Adiabatic? (Reason)
10 / 64
Adiabatic Law: Differential & Linear Form
Start with
PV γ = constant
Differentiate:
dP dV
d(PV γ ) = 0 ⇒ +γ =0
P V
Thus (for small changes around P0 , V0 ):
dV
dP = −γP0
V0
Let p be the incremental pressure and τ the incremental volume:
τ
p = −γP0
V0
11 / 64
Linearization with Time Derivatives (Useful Form)
dp P0 dτ
= −γ
dt V0 dt
Meaning
Pressure increase is proportional to the rate of fractional compression.
12 / 64
Bulk Modulus Interpretation
Key Insight
Large B (stiff medium) ⇒ large pressure change for the same compression.
13 / 64
Continuity: Geometry of Volume Change
air slice
ξ(x, t) ξ(x + ∆x, t)
x
∆x
14 / 64
Continuity Equation (Derivation)
Approximate
∂ξ τ ∂ξ
τ≈ ∆x ∆y ∆z ⇒ ≈
∂x V0 ∂x
Take time derivative and use u = ∂ξ/∂t:
∂ τ ∂ ∂ξ ∂u
= =
∂t V0 ∂x ∂t ∂x
15 / 64
Combine Gas Law + Continuity (Pressure–Velocity Link)
16 / 64
Insert Equation of Motion & Derive the Wave Equation
Differentiate in time:
∂2p
∂ ∂u
= −γP0
∂t 2 ∂x ∂t
From equation of motion:
∂u ∂p ∂u 1 ∂p
ρ0 =− ⇒ =−
∂t ∂x ∂t ρ0 ∂x
Substitute:
∂2p 1 ∂2p
2
= γP0
∂t ρ0 ∂x 2
Rearrange:
∂2p 1 ∂2p γP0
2
= 2 2 c2 =
∂x c ∂t ρ0
17 / 64
Interpretation of Wave Speed
s
γP0
c=
ρ0
Meaning
Sound travels faster in stiffer media and slower in denser media (all else equal).
18 / 64
3D Pressure Wave Equation (Extension)
In three dimensions:
1 ∂2p
∇2 p =
c 2 ∂t 2
where ∇2 = ∂xx + ∂yy + ∂zz .
z
plane-wave slice
19 / 64
General Solution (1D)
For
∂2p 1 ∂2p
2
= 2 2
∂x c ∂t
the general solution is
x x
p(x, t) = f1 t − + f2 t +
c c
x
f1 t − c : wave traveling in +x
x
f2 t + c : wave traveling in −x
20 / 64
Traveling-Wave Picture
t = t0
t = t0 + ∆t
moves right
21 / 64
Steady-State Harmonic Form
22 / 64
Harmonic Solution (Plane Progressive Wave)
p̃
ũ =
ρ0 c
23 / 64
Characteristic (Specific) Acoustic Impedance
Important
In a progressive plane wave, p̃ and ũ are in phase, so Z0 is purely resistive.
24 / 64
Characteristic Impedance: Phasor Picture
Im
p̃
Re
ũ
25 / 64
Rigidly Closed Tube: Setup
x
l
26 / 64
Boundary Condition at the Rigid End
u(l, t) = 0 ⇒ ũ(l) = 0
27 / 64
Standing-Wave Form (Pressure and Velocity)
28 / 64
Rigid Wall Standing Wave: Snapshot Diagram
p rigid wall
t=0
x
pressure nodes occur every λ/2
29 / 64
Resonance Condition (Modes in a Closed-Closed Tube)
For a tube closed at both ends, velocity nodes occur at both ends, giving:
nλ
l= (n = 1, 2, 3, . . . )
2
Fundamental: n = 1 ⇒ λ1 = 2l
Harmonics: n = 2, 3, . . .
30 / 64
Mode Shape (Velocity) — Simple TikZ Sketch
n = 1 (fundamental)
31 / 64
Specific Acoustic Impedance in a Tube
Purely reactive (imaginary): energy is stored and returned (no net transport).
Magnitude can become very large near velocity nodes.
32 / 64
Reactance Variation Along a Rigidly Closed Tube
6
Xs →
(reactive
+∞ part)
4
2
x (along the tube)
λ λ 3λ 2λ
2 2
−2
−4
→ −∞
−6
Caption: The reactive part of the tube’s specific acoustic impedance varies approximately as Xs ∝ − cot(kx). At positions
where sin(kx) = 0 (every λ/2), the reactance tends to ±∞ (vertical asymptotes). Between successive half-wavelengths
the sign of Xs alternates, indicating alternating compliant (spring-like) and inertive (mass-like) behaviour.
33 / 64
Low-Frequency Approximation (Interpretation)
If k(l − x) is small:
1 θ
cot θ ≈ − (θ ≪ 1)
θ 3
So a short air volume behaves like an acoustic compliance (spring-like storage).
Meaning
At low frequencies, air mainly stores potential energy (compression), like a spring.
34 / 64
Freely Traveling Spherical Wave (Concept)
35 / 64
Chapter Two Summary (Merged)
36 / 64
Quick Practice
3
1 If ρ0 = 1.2 kg/m and c = 343 m/s, compute Z0 .
2 In a rigidly closed tube, where do velocity nodes occur?
3 Explain why sound in air is closer to adiabatic than isothermal.
4 A point source doubles its distance from you. What happens to pressure amplitude (approx.)?
37 / 64
Tips
38 / 64
Q1(a)(i): Propagation of Sound
(a)(i) With a clear labelled illustration, explain how sound propagates through a medium.
Wave propagation (c)
39 / 64
Q1(a)(i): Solutions
3 Cause of motion: pressure gradient produces force on a small air element ⇒ acceleration.
4 Energy transport: disturbance travels with speed
s
B
c= (elasticity + inertia).
ρ0
Key sentence
Sound propagates because the medium has bulk modulus (restoring force) and density (inertia).
40 / 64
Q1(b)(i): From Displacement to Pressure
(b)(i) The displacement is S(x, t) = a cos(kx − ωt) in a medium of bulk modulus B. Derive the pressure
fluctuation p ′ (x, t).
Expected idea
Pressure fluctuation is proportional to strain (fractional compression):
∂S ∂S
strain ≈ ⇒ p ′ ∝ −B .
∂x ∂x
41 / 64
Q1(b)(i): Solution Plan
1 Strain (1D): for a longitudinal wave, the fractional change in length/volume is approximately
∂S
ϵ= .
∂x
2 Bulk modulus definition:
∆P ∆V
B = −V ⇒ p ′ = −B .
∆V V
For 1D motion, ∆V /V ≈ ϵ = ∂S/∂x.
3 Substitute S(x, t) = a cos(kx − ωt):
∂S
= −ka sin(kx − ωt).
∂x
4 Hence
∂S
p ′ = −B = Bka sin(kx − ωt).
∂x
Note on signs
Some texts write p ′ = +B ∂S/∂x depending on how compression is defined. State your convention and keep it
consistent.
42 / 64
Q1(b)(ii): Max Displacement from Max Pressure
(b)(ii) Max pressure variation pmax = 30 Pa. Find max displacement amplitude a in air: ρ0 = 1.22 kg m−3 ,
f = 1 kHz, c = 330 m s−1 .
43 / 64
Q1(b)(ii): Worked Solution (Units Shown)
pmax 30
umax = = ≈ 7.46 × 10−2 m s−1
ρ0 c (1.22)(330)
ω = 2πf = 2π(1000) ≈ 6.283 × 103 rad s−1
umax 7.46 × 10−2
a= ≈ ≈ 1.19 × 10−5 m
ω 6.283 × 103
a ≈ 1.2 × 10−5 m = 12 µm
44 / 64
Q1(c)(i): Why Adiabatic in Air?
PV γ = constant.
45 / 64
Q2(a): Derive EOM and Continuity
(a)(i) Derive the equation of motion. (a)(ii) Derive the equation of continuity.
Targets
∂u ∂p ′ ∂ρ′ ∂u
ρ0 =− , + ρ0 = 0.
∂t ∂x ∂t ∂x
46 / 64
Q2(a): Solutions / Solution Plan (Clear Steps)
Equation of motion (Newton’s 2nd law):
1 Choose an air element of volume V = A∆x.
2 Pressure forces: p(x)A on left and p(x + ∆x)A on right.
3 Net force:
∂p ′ ∂p ′
Fx = p(x)A − p(x + ∆x)A ≈ − A∆x = − V.
∂x ∂x
4 Mass m = ρ0 V , acceleration ∂u/∂t.
5 Conclude:
∂u ∂p ′
ρ0 =− .
∂t ∂x
Continuity:
1 Start from 1D continuity:
∂ρ ∂(ρu)
+ = 0.
∂t ∂x
2 Linearize with ρ = ρ0 + ρ′ and neglect products ρ′ u:
∂ρ′ ∂u
+ ρ0 = 0.
∂t ∂x
47 / 64
Q2(b): Adiabatic Law to Wave Equation
Final destination
∂ 2 p′ 1 ∂ 2 p′ γP0
2
= 2 , c2 = .
∂x c ∂t 2 ρ0
48 / 64
Q2(b): Solution Plan
1 Start with PV γ = C . Differentiate:
2 Divide by PV γ :
dP dV
+γ = 0.
P V
3 For small oscillations about (P0 , V0 ) and V = V0 − τ :
τ ∂p ′ γP0 ∂τ
p ′ = −γP0 ⇒ =− .
V0 ∂t V0 ∂t
4 Relate volume strain to displacement:
τ ∂ξ ∂ξ
≈ , u= .
V0 ∂x ∂t
5 Combine with:
∂u ∂p ′ ∂ρ′ ∂u
ρ0 =− and + ρ0 =0
∂t ∂x ∂t ∂x
′
to eliminate u (and ρ ) and obtain the wave equation.
49 / 64
Q11(a): General Solution of the 1D Wave Equation
Meaning
f1 : wave travelling in +x direction; f2 : wave travelling in −x direction.
50 / 64
Q11(a):
1 Start with
∂2p 1 ∂2p
2
= 2 2.
∂x c ∂t
2 Factorize the operator:
∂ 1 ∂ ∂ 1 ∂
− + p = 0.
∂x c ∂t ∂x c ∂t
3 Hence either factor acting on p gives a first-order PDE. Integrating yields:
x x
p = f1 t − + f2 t + .
c c
51 / 64
Q7(c): Derive Zs = p/u for a Plane Progressive Wave
(c) Derive Zs = p/u for a freely travelling plane progressive wave and show Zs = ρ0 c.
52 / 64
Q7(c): Solution Plan (Phasor Substitution Shown)
53 / 64
Q8: Rigidly Closed Tube (Structured)
sin k(l − x)
u(x) = u0 .
sin(kl)
(a) State boundary conditions. (b) Derive specific acoustic impedance Zs at l ′ = l − x. (c) For l ′ < λ/20, show
′
1
Zs ≈ jωC a
, where Ca = ρ l c 2 . (d) Locate first three pressure/velocity nodes for given f , l. (e) Sketch pattern for
0
n = 2.
54 / 64
Q8: Tube Diagram
x
length l
55 / 64
Q8(a,b): Solutions
u(l) = 0.
u(0) = u0 .
p(l ′ )
Zs (l ′ ) = = −jρ0 c cot kl ′ .
u(l ′ )
(Derivation route: write incident + reflected plane waves, apply u = 0 at rigid termination, then take ratio p/u.)
56 / 64
Q8(c): Small-Length Approximation (Explain the Physics)
Given:
Zs (l ′ ) = −jρ0 c cot kl ′ .
Physical meaning
For very short air columns, the tube behaves like a spring (compliance) that stores potential energy.
57 / 64
Q12(a): Spherical Progressive Wave Impedance
(a) For a spherical free-progressive wave, write p(r , t) and ur (r , t) and derive:
jkr
Zs = ρ0 c .
1 + jkr
58 / 64
Q12(a): Solution Plan (Show the Missing Step Everyone Forgets)
p̂ jkr
= ρ0 c .
ûr 1 + jkr
59 / 64
p2
Q15(b): Derive ⟨I ⟩ = ρ0 c
p2
(b) Derive I = for plane progressive waves.
ρ0 c
Core identities
p̂
I (t) = p ′ (t) u(t), Z0 = = ρ0 c.
û
60 / 64
Q15(b): Solution Plan (RMS vs Peak Made Clear)
1 Instantaneous intensity:
I (t) = p ′ (t) u(t).
2 For a progressive plane wave, p ′ and u are in phase and
p ′ (t)
u(t) = .
ρ0 c
3 Hence:
[p ′ (t)]2
I (t) = .
ρ0 c
4 Time-average over one cycle gives:
2
prms
⟨I ⟩ = .
ρ0 c
√
If you use peak amplitude p̂, then prms = p̂/ 2.
61 / 64
Q20(a,b): Full First-Principles Derivation (Chapter Two Master)
(a) Derive continuity and equation of motion using a small box of air. (b) Combine with adiabatic law to obtain
the full pressure wave equation (show every differentiation step).
Exam note
Write the derivation in a strict logical order (definitions → laws → elimination → wave equation).
62 / 64
Q20(a,b): Perfect Answer Structure (Step-by-Step)
Write in this order (and box each final equation):
1 Define perturbations:
p = P0 + p ′ , ρ = ρ0 + ρ′ , |p ′ | ≪ P0 , |ρ′ | ≪ ρ0 .
2 Equation of motion (Newton):
∂u ∂p ′
ρ0 =− .
∂t ∂x
3 Continuity (linearized):
∂ρ′ ∂u
+ ρ0 =0 .
∂t ∂x
4 Adiabatic state equation:
γP0
p ′ = c 2 ρ′ , c2 = .
ρ0
5 Differentiate continuity w.r.t. t, substitute EOM, then eliminate ρ′ using p ′ = c 2 ρ′ .
6 Arrive at the pressure wave equation:
∂ 2 p′ 1 ∂ 2 p′
2
= 2 .
∂x c ∂t 2
63 / 64
Quick Drill (Revision Check)
64 / 64