Unit 5: SOUND GENERATION, PROPAGATION, AND RADIATION IN/FROM AN AEROENGINE NACELLE

ASE417 — Aeroacoustics 3 min read

I. Orientation — Governing Framework of Aeroengine Duct Acoustics

Aeroengine nacelle acoustics describes how pressure disturbances generated mainly by the fan propagate through annular inlet and bypass ducts, interact with mean flow and acoustic liners, and radiate into the far field. For small disturbances, the governing equations follow from linearized conservation of mass, momentum, and energy.

Acoustic variables: The total pressure, density, and velocity are decomposed into steady mean and small fluctuating parts:
TEXT
```
  p_total = p̄ + p′,    ρ_total = ρ̄ + ρ′,    u_total = Ū + u′
```
Here, (p̄,\rhō,Ū) are mean quantities and (p′,\rho′,u′) are acoustic perturbations.
Linear-acoustic assumption: Products such as (p′u′) are neglected in the governing field equations because perturbation amplitudes are small compared with mean quantities.
Wave equation: In a stationary, uniform, inviscid medium, acoustic pressure satisfies:
TEXT
```
  ∇²p′ − (1/c₀²)(∂²p′/∂t²) = 0
```
where (p′) is acoustic pressure, (c₀) is sound speed, (t) is time, and (\nabla^2) is the spatial Laplacian.
Harmonic convention: A single-frequency field is commonly represented as (p′=\Re{\hat p e^{i\omega t}}), where (\hat p) is complex pressure amplitude, (i^2=-1), and (\omega=2\pi f) is angular frequency for frequency (f).
Important nacelle effects:
- Geometry: Area variation and annular cross-sections alter modal propagation.
- Mean flow: Inlet or exhaust flow changes wavelength, cut-on conditions, and radiation direction.
- Liners: Impedance-treated walls absorb selected frequencies.
- Boundaries: Fan faces, duct openings, bifurcations, and nozzles reflect or scatter waves.

II. Duct Acoustic Fields — Modes, Cut-On, and Energy Transport

A duct confines sound laterally, so its acoustic field separates into cross-sectional eigenfunctions and axial travelling or decaying waves.

A. Basic theory of sound propagation in ducts

The central principle is that duct sound propagates as discrete modes rather than as an unrestricted plane wave.

Modal representation: For a uniform duct aligned with coordinate (x), harmonic pressure can be expanded as:
TEXT
```
  p̂(x,r,θ) = Σ A_mn Φ_mn(r,θ)e^(−ik_x,mn x)
```
where (A{mn}) is modal amplitude, (\Phi{mn}) is a cross-sectional eigenfunction, (m) is circumferential order, (n) is radial order, and (k_{x,mn}) is axial wavenumber.
Rigid-wall condition: Zero normal velocity at a solid wall gives:
TEXT
```
  ∂p̂/∂n_w = 0
```
where (n_w) denotes the outward wall-normal direction. This condition determines the allowable transverse eigenvalues.
Dispersion relation without flow:
TEXT
```
  k_x,mn² = k₀² − μ_mn²,    k₀ = ω/c₀
```
Here, (k₀) is the free-space acoustic wavenumber and (\mu_{mn}) is the transverse eigenvalue fixed by duct shape and boundary conditions.
Cut-on and cut-off:
1. Cut-on mode: If (k₀>\mu{mn}), then (k{x,mn}) is real and the mode transports acoustic energy axially.
2. Cut-off mode: If (k₀<\mu{mn}), then (k{x,mn}) is imaginary and amplitude decays exponentially with distance.
Plane-wave mode: The lowest rigid-duct mode has (\mu{00}=0), uniform pressure over the cross-section, and (k{x,00}=k₀). Higher modes become cut-on as frequency increases.
Circular and annular ducts: Circular-duct eigenfunctions involve Bessel functions; annular-duct eigenfunctions combine Bessel functions of the first and second kinds to satisfy conditions at both hub and casing radii.
Mean-flow dispersion: For uniform axial flow of Mach number (M),
TEXT
```
  (ω − Mc₀k_x)² = c₀²(k_x² + μ²)
```
where (k_x) is axial wavenumber and (\mu) is transverse eigenvalue. Upstream and downstream waves therefore have different wavelengths.
Acoustic power: Time-averaged axial power is obtained from acoustic intensity:
TEXT
```
  W = (1/2) Re ∫_S p̂ û_x* dS
```
where (W) is acoustic power, (S) is duct area, (\hat u_x) is axial velocity amplitude, and the asterisk denotes complex conjugation.

III. Nacelle Acoustic Environment — Internal Transmission and External Radiation

An aeroengine nacelle is an annular, non-uniform, lined duct containing mean flow, rotating machinery, structural obstructions, and open boundaries.

A. Sound propagation in an aeroengine nacelle

Nacelle propagation is governed by the evolution and scattering of duct modes between the fan, inlet lip, bypass nozzle, and exterior field.

Propagation paths:
1. Forward radiation: Fan-generated sound travels through the inlet and radiates from the nacelle lip.
2. Rearward radiation: Sound travels through the bypass duct and exits through the fan nozzle; turbine and jet sources may contribute separately.
Annular modal structure: Hub and casing boundaries produce circumferential-radial modes ((m,n)). Large (|m|) or (n) generally gives a higher cut-on frequency and a more complex cross-sectional pressure pattern.
Mode scattering: Area changes, pylons, bifurcations, liners, and non-axisymmetric flow redistribute energy between modes. In an exactly axisymmetric, uniform duct, circumferential order (m) remains uncoupled.
Liner boundary condition: A locally reacting liner is represented by acoustic impedance:
TEXT
```
  Z(ω) = p̂/û_n = R + iX
```
where (\hat u_n) is wall-normal velocity, (R) is resistance responsible for dissipation, and (X) is reactance associated with stored acoustic energy.
Flow–liner interaction: Grazing flow alters effective impedance and can produce boundary-layer refraction. Liner performance therefore depends on frequency, sound level, flow Mach number, and modal incidence angle.
Inlet radiation: At the open inlet, internal modes couple to free-space waves. A cut-on mode radiates preferentially at an angle related approximately to its axial and transverse wavenumber components; cut-off modes contribute mainly near the lip.
Flight effects: Forward motion and inlet flow cause convection and refraction. They modify apparent directivity and create different propagation conditions on the upstream and downstream sides.
Installation effects: The wing, fuselage, and pylon can reflect or shield nacelle noise. Consequently, isolated-engine directivity may differ substantially from aircraft-level radiation.

IV. Transfer-Based Modelling — Relating Acoustic States Across Components

The transfer element method represents each duct segment or acoustic component by an algebraic relation between state variables at its inlet and outlet.

A. Fundamental idea of the transfer element method

The method reduces a distributed wave problem to a sequence of compact elements whose matrices can be assembled to predict system transmission.

State vector: For plane-wave modelling, a convenient state is:
TEXT
```
  q(x) = [p̂(x), Û(x)]ᵀ
```
where (\hat p) is pressure amplitude, (\hat U=S\hat u_x) is acoustic volume velocity, (S) is cross-sectional area, and superscript (T) denotes transpose.
Element relation:
TEXT
```
  q₂ = T_e q₁
```
where (q_1) and (q_2) are states at the element boundaries and (T_e) is the element transfer matrix.
Uniform lossless duct element:
TEXT
```
  T_e = [ cos(k₀L)       −iZ_c sin(k₀L) ]
        [ −i sin(k₀L)/Z_c   cos(k₀L)    ]
```
Here, (L) is element length and (Z_c=\rho_0c_0/S) is characteristic impedance based on volume velocity.
Assembly: For consecutive elements,
TEXT
```
  T_total = T_N T_(N−1) ... T₂ T₁
```
where (T_j) is the matrix of element (j) and (N) is the number of elements. Matrix order follows the physical propagation sequence.
Boundary conditions: Source impedance, termination impedance, prescribed pressure, or radiation impedance closes the system. Reflection and transmission coefficients then follow from the solved boundary states.
Strengths and limitations: The method is computationally efficient for quasi-one-dimensional systems, but a two-variable state cannot capture strong higher-mode coupling; multimodal transfer or scattering matrices are then required.

V. Non-Uniform Duct Elements — Discretization of Area Variation

A varying duct is modelled by preserving pressure and volume-flow continuity while allowing local area and characteristic impedance to change.

A. Construction of transfer element for a varying cross-section duct

The transfer element is constructed from the governing horn equation or from a product of short locally uniform segments.

Webster horn equation: For slowly varying area (S(x)),
TEXT
```
  (1/S)d/dx[S(dp̂/dx)] + k₀²p̂ = 0
```
where (x) is axial position and (S(x)) is local area. The approximation assumes predominantly plane, locally axial propagation.
First-order state equations:
TEXT
```
  dp̂/dx = −iωρ₀Û/S
  dÛ/dx = −iωS p̂/(ρ₀c₀²)
```
where (\rho_0) is mean density. These equations directly generate a spatial transfer operator.
Segmented construction:
1. Divide the duct into short lengths (\Delta x_j).
2. Evaluate area (S_j) at each segment midpoint.
3. Form the local uniform-duct matrix using (L=\Delta xj) and (Z{c,j}=\rho_0c_0/S_j).
4. Multiply the matrices in propagation order.
Interface conditions: Pressure and volume velocity remain continuous at an ideal abrupt area junction:
TEXT
```
  p̂_left = p̂_right,    Û_left = Û_right
```
Particle velocity itself changes because (\hat u_x=\hat U/S).
Accuracy control: Segment length should be small relative to both wavelength and geometric variation scale. Rapid area changes, separated flow, or significant transverse modes require multidimensional finite-element or mode-matching treatment.

VI. Fan Acoustic Sources — Tonal and Broadband Representation

Fan source modelling specifies the modal pressure or velocity field injected at a source plane so that nacelle propagation and radiation can be calculated separately.

A. Fan noise source modelling

A useful fan model distinguishes deterministic blade-passing tones from stochastic broadband noise and assigns each component a frequency, mode content, amplitude, and phase.

Blade-passing frequency:
TEXT
```
  f_BPF = BN/60
```
where (B) is rotor blade count and (N) is rotational speed in revolutions per minute. Harmonics occur at (n f_{BPF}), where (n) is a positive integer.
Rotor–stator interaction modes: Circumferential orders commonly satisfy the Tyler–Sofrin relation:
TEXT
```
  m = nB − sV
```
where (m) is circumferential mode order, (V) is stator-vane count, and (s) is any integer indexing spatial aliases. Only cut-on members efficiently propagate.
Source mechanisms:
- Tonal loading noise: Periodic blade forces and rotor–stator wake interaction generate coherent BPF harmonics.
- Broadband noise: Turbulence ingestion, boundary-layer interaction, tip flow, and wake turbulence generate continuous spectra.
- Thickness noise: Periodic displacement of fluid by rotating blades contributes deterministic radiation.
Modal source specification: At a chosen fan plane, the source may be represented by modal amplitudes (A{mn}^{+}) and (A{mn}^{-}), where superscripts indicate downstream- and upstream-travelling waves.
Broadband statistics: Random modal sources are described using a cross-spectral density matrix:
TEXT
```
  C_ab(f) = E[A_a(f)A_b*(f)]
```
where (C_{ab}) measures frequency-dependent coherence between modes (a) and (b), and (E[\cdot]) denotes ensemble averaging.
Coupling to nacelle models: Source modal vectors are propagated through transfer or scattering matrices, combined with liner losses and reflections, and matched to an external radiation model to obtain acoustic power and directivity.
Model limitations: Predictions depend strongly on assumed blade loading, turbulence spectra, phase relationships, and source location. High-fidelity computational aeroacoustics or calibrated rig data is needed when nonlinear flow, shocks, or strong rotor–stator coupling dominates.