DTU — Acoustics Laboratory · 1989 M.Sc. Thesis Archive · §3.4

Sound Propagation in Ducts
with Absorbent Linings Figure 3.10(b) — Modal attenuation for a bulk-reacting liner

Abstract
This page re-computes Figure 3.10(b) from the thesis: the axial attenuation of the (0,0) and (0,2) modes in a rectangular duct lined on two parallel sides with a bulk-reacting absorber. Scott's transcendental equation is solved for each frequency using Newton–Raphson for the (0,0) mode and an arctangent-fixed-point iteration for the (0,2) mode, after the method of Christie (1971). Lining properties are taken from the Delany–Bazley formulae as modified by Mechel. The composite least-attenuation curve is the envelope min{Re(Γ)} of the two modes. Move the sliders to explore the parameter space.

Controls

Result

Computing…
attenuation · dB per unit duct width FIG. 3.10(b) — RECONSTRUCTED
(0,0) mode
(0,2) mode
least-attenuated envelope
Figure 3.10(b). Produced using the 'Extended' algorithm (re-implemented in JavaScript). Γ is the attenuation; the dashed curve is the plane-wave (0,0) mode, the solid curve the (0,2) mode. Current parameters shown below.
σ
8500 N·s/m⁴
ℓ
0.250 m (duct half-width)
d
0.050 m (lining thickness)
b = d/ℓ
0.200
resolution
20 Hz
c
344.25 m/s (air, 20 °C)
ρ
1.205 kg/m³

Method

For one mode of propagation the acoustic potentials Φ (air) and Θ (lining) share a common axial propagation constant Γ. Applying continuity of normal particle velocity and acoustic pressure at the lining–air interface, with a rigid outer wall, leads to Scott's transcendental equation, rewritten by Christie in the dimensionless form:

f(w) = w·tan w + M·√(w² + p²)·tan(b·√(w² + p²)) = 0 (2.2.9)

where w = √((Γℓ)² + (kℓ)²), p = √((hℓ)² − (kℓ)²), b = d/ℓ, and M = ρ/ρ′ is the ratio of air density to complex lining density. The lining wave number h and characteristic impedance Za are obtained from the Delany–Bazley/Mechel formulae with normalised flow resistivity A = σ/(ρ f).

Once w is found, the duct propagation constant follows from Γ = √((w/ℓ)² − k²), and the attenuation is 20·log₁₀(e)·Re{Γ} ≈ 8.686·Re{Γ} dB per metre. Here it is normalised to the duct half-width ℓ.

The iteration is seeded from the Morse (local-reaction) starting points w(0,0) ≈ 0 and w(0,2) ≈ π, then continued across the frequency axis using linear extrapolation of the previous solution as the next initial value — exactly as in the 1988 BASIC code of the original thesis.