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

Acoustic modes are independent spatial patterns in which sound can propagate through a duct. For a uniform duct, the acoustic pressure satisfies the Helmholtz equation:

where is the acoustic wavenumber. The pressure may be expanded as

where:

• is the transverse mode shape.
• is the axial wavenumber.
• and are the downstream- and upstream-travelling modal amplitudes.
• The rigid-wall boundary condition is at the duct wall.

The transverse eigenfunctions are orthogonal, allowing the total sound field to be separated into individual modes.

For small acoustic disturbances, write pressure and density as and . The linearized governing equations are:

• Continuity equation:
  

• Momentum equation:
  

• Isentropic relation:
  

Differentiate continuity with respect to time and substitute the divergence of the momentum equation:

Using gives

For harmonic motion, , this becomes the Helmholtz equation:

Boundary conditions at the duct wall determine the permitted acoustic modes.

For a duct mode with transverse eigenvalue , the axial wavenumber in a stationary medium is

• A mode is cut-on when . In this case, is real and the mode propagates while carrying acoustic energy along the duct.
• A mode is cut-off when . Then is imaginary, and the mode decays exponentially with axial distance.
• At cut-off, and .

The cut-off angular frequency is

or

The value of depends on duct geometry, dimensions, wall boundary conditions, and mode order. The plane mode has and is therefore cut-on at all nonzero frequencies in an ideal uniform duct.

In cylindrical coordinates , a circular-duct pressure mode can be written as

where:

• is the circumferential mode order.
• is the radial mode order.
• is the Bessel function of the first kind.
• is the transverse eigenvalue.
• is the axial wavenumber.

For a rigid wall of radius , the normal acoustic velocity must vanish. Therefore,

which gives

Each pair defines a distinct spinning and radial pressure pattern. Modes with rotate circumferentially, while the mode is the plane mode.

Uniform mean flow convects acoustic disturbances and makes upstream and downstream propagation different. For axial mean velocity , the modal dispersion relation is

Its main consequences are:

• Downstream-propagating waves are convected by the flow and have a different axial wavelength from upstream waves.
• Upstream waves encounter the flow and may have shorter wavelengths.
• For a plane wave, the approximate axial wavenumbers are
  
• The cut-on condition is modified. For a uniform subsonic flow, a mode is cut-on when
  
• Mean flow changes acoustic energy flux, impedance, refraction, and liner attenuation.

Velocity and temperature gradients in a real nacelle add refraction and mode-coupling effects beyond the uniform-flow model.

Hard-wall duct:

• The wall is assumed perfectly rigid.
• Normal acoustic velocity is zero:
  
• From linear momentum, this corresponds to
  
• In the ideal model, the wall absorbs no acoustic energy.

Acoustically lined duct:

• The wall has a finite acoustic impedance .
• Pressure and normal velocity satisfy
  
• Equivalently, for harmonic motion in a stationary fluid,
  
  with the sign depending on the normal and time conventions.
• The impedance is generally complex: its resistance represents energy absorption, while its reactance represents stored acoustic energy.

A liner therefore produces complex modal wavenumbers and attenuates selected frequency ranges, whereas a hard wall mainly guides the sound.

Fan-generated sound propagating toward the inlet travels through an annular or nearly circular duct before radiating from the intake lip. Important effects include:

• Modal propagation: The source excites circumferential and radial duct modes, each with its own cut-off frequency and axial wavenumber.
• Mean-flow convection: Intake flow modifies upstream modal wavelengths, energy flux, and cut-on behavior.
• Area variation: Changes in inlet area cause reflection, transmission, and coupling between modes.
• Acoustic liners: Treated surfaces absorb sound over selected frequency bands.
• Boundary-layer effects: Nonuniform flow near the wall alters liner response and modal attenuation.
• Intake-lip scattering: Duct modes are diffracted at the lip and converted into free-field radiation.
• Geometry and incidence: Intake shape and flight condition affect refraction and directivity.

The observed far-field inlet noise is therefore determined jointly by the fan source, internal nacelle transmission, and radiation at the inlet opening.

Inlet side:

• Sound propagates upstream against the intake flow.
• The inlet lip strongly influences diffraction and forward-arc directivity.
• The duct is often relatively cool, so sound speed and density variations are moderate.
• Inlet liners are commonly used to attenuate fan tones and broadband noise.

Exhaust side:

• Sound generally propagates downstream with the mean flow.
• The annular bypass duct, outlet guide vanes, bifurcations, and nozzle can scatter or couple modes.
• Temperature and velocity gradients may be stronger and can significantly refract sound.
• Radiation from the bypass or core nozzle interacts with jet noise and shear layers.

In both paths, modal cut-on, duct area changes, impedance treatments, and termination geometry determine transmission. However, the opposite propagation direction relative to flow and the different thermal and geometric environments lead to different attenuation and directivity.

At an open nacelle termination, a guided duct mode encounters a sudden change from the internal duct to the external domain. It is partly reflected and partly converted into outward-propagating waves through diffraction at the lip.

Far-field directivity depends on:

• Modal order: Different circumferential and radial modes produce different radiation lobes.
• Cut-on ratio: A mode close to cut-off tends to radiate at large angles, while a strongly cut-on mode generally radiates closer to the duct axis.
• Frequency and aperture size: The Helmholtz number determines how directional the opening is.
• Termination geometry: Lip thickness, curvature, scarfing, and nozzle shape affect diffraction.
• Mean flow: Flow changes phase matching and convects or refracts the radiated field.
• Modal phase and amplitude: Interference among simultaneously radiated modes modifies the total pattern.

Radiation efficiency can be represented through a modal radiation impedance, whose real part measures radiated acoustic power and whose imaginary part represents near-field energy storage.

The transfer element method divides a complex duct or nacelle into simpler axial elements. Each element relates acoustic state variables at its inlet and outlet through a transfer matrix:

For a plane-wave model, a common state vector is

where is acoustic pressure and is acoustic volume velocity. For multimodal analysis, the state contains vectors of modal pressure and velocity amplitudes.

The complete system is assembled by cascading element matrices:

The method can represent:

• Uniform and varying-area duct sections.
• Liners and impedance boundaries.
• Area discontinuities and junctions.
• Mean-flow effects.
• Source and termination conditions.

Once the global relation is formed, reflection, transmission, insertion loss, and radiated sound can be calculated.

For harmonic plane waves in a uniform duct of area ,

With the convention, the corresponding volume velocity is

Eliminating and between the states at and gives

where the characteristic impedance based on volume velocity is

Thus,

For a lossless reciprocal element, . Reversing the harmonic convention changes the signs of the imaginary terms but not the physical result.

If a duct consists of consecutive elements, each element satisfies

Successive substitution gives

Important assembly points are:

• Matrix order follows the physical sequence of elements.
• Pressure and volume velocity continuity are imposed at ordinary interfaces.
• Junction matrices are introduced for abrupt area changes, branches, liners, or other discontinuities.
• A prescribed source may be represented by an incident-wave amplitude, pressure, velocity, impedance, or an additional source vector.
• An outlet impedance supplies the condition
  
• A radiation impedance may be used at an open nacelle termination.

Solving the resulting linear equations gives inlet reflection, outlet transmission, internal pressure, and acoustic power. For many modes or long lossy ducts, scattering matrices may be numerically more stable than direct transfer-matrix multiplication.

Transfer matrix:

• Relates the complete acoustic state at one axial station to another.
• Is easily cascaded by ordinary matrix multiplication.
• Is convenient for compact one-dimensional or low-order models.
• Can become ill-conditioned when strongly evanescent modes cause exponentially large and small matrix entries.

Scattering matrix:

• Relates outgoing wave amplitudes to incoming wave amplitudes:
  
• Directly provides reflection and transmission coefficients.
• Is generally more stable for long systems, lossy liners, and many evanescent modes.
• Requires an appropriate matrix composition rule when sections are cascaded.

Thus, transfer matrices are conceptually simple and efficient for modest systems, whereas scattering matrices are preferred when numerical stability is critical.

Let the duct area be , acoustic pressure be , particle velocity be , and volume velocity be . Under the quasi-one-dimensional, lossless, no-mean-flow approximation, the harmonic momentum equation is

The harmonic continuity equation is

Therefore, the state equation is

Eliminating gives Webster's horn equation:

An element transfer matrix is obtained by integrating the first-order state equation across the element. If is approximated as constant over a short element, the uniform-duct matrix may be used with a representative area. More accurate elements integrate the varying coefficients directly or use analytical horn solutions.

A practical construction procedure is:

1. Define the geometry: Specify and divide the duct into short axial elements.
2. Choose the state vector: Commonly use .
3. Form the state equation:
   
   where
   
4. Approximate each element: Evaluate the area at the midpoint, interpolate it, or use an exact conical or exponential-horn solution.
5. Compute the element matrix: For a short element of length ,
   
6. Cascade the matrices:
   
7. Refine the mesh: Reduce until transmission and reflection results converge.

The elements should be short relative to both the acoustic wavelength and the geometric variation scale. Multimodal models additionally require overlap integrals to represent coupling between local modes.

A varying duct area changes the local acoustic impedance and modal eigenfunctions. This causes two main effects:

• Reflection: For plane waves, the local characteristic impedance based on volume velocity is
  
  A rapid change in creates an impedance mismatch, so part of an incident wave is reflected. Gradual area changes generally reduce reflection.

• Mode coupling: In a multimodal duct, the transverse mode shapes depend on the local cross-section. As the geometry varies, a mode entering the section may project onto several local modes. Energy can therefore transfer between radial or circumferential mode orders.

Coupling becomes stronger when:

• Geometry changes rapidly compared with the wavelength.
• Modes have similar axial wavenumbers.
• The cross-section becomes asymmetric.
• Struts, bifurcations, or liners disturb circumferential uniformity.

A sufficiently slow and symmetric variation may be treated by an adiabatic approximation in which each mode evolves with little conversion to other modes.

The principal fan noise mechanisms are:

• Rotor-alone tonal noise: Periodic blade loading and thickness effects generate tones at the blade-passing frequency and its harmonics.
• Rotor-stator interaction noise: Rotor wakes and potential disturbances interact with outlet guide vanes, producing strong discrete tones.
• Inflow-distortion interaction: Atmospheric turbulence, boundary layers, or distorted intake flow create unsteady blade loading.
• Broadband turbulence-interaction noise: Random turbulence interacting with blades and vanes produces a continuous spectrum.
• Trailing-edge noise: Turbulent boundary layers scatter at blade trailing edges.
• Tip-clearance noise: Leakage vortices and their interaction with the casing and blade passages generate tonal and broadband components.
• Shock-associated noise: At transonic or supersonic blade-tip speeds, shocks and shock nonuniformities can produce multiple pure tones, often called buzz-saw noise.

The resulting source field is distributed over the fan and is most effectively represented in terms of frequency, circumferential order, radial mode content, amplitude, and phase.

If a rotor with blades rotates at revolutions per second, the blade-passing frequency is

The th harmonic occurs at

For interaction between a rotor having blades and a stator having vanes, the Tyler-Sofrin relation gives possible circumferential mode orders as

where is the rotor harmonic index and is any integer representing the spatial harmonic of the vane row.

The mode order determines the number and rotational sense of circumferential pressure lobes. However, generation of a circumferential order does not guarantee far-field radiation. The mode must also have at least one cut-on radial order at the tone frequency. Cut-off modes decay inside the nacelle, while cut-on modes can propagate to the inlet or exhaust and radiate.

A fan source may be represented at a reference plane by modal amplitudes rather than by a single acoustic pressure. The source vector can be written as

where each complex amplitude specifies the magnitude and phase of a circumferential-radial mode.

A complete model includes:

• Tone frequencies, such as the blade-passing frequency and its harmonics.
• Circumferential orders selected by rotor-stator interaction.
• Radial mode amplitudes and phases.
• Separate upstream- and downstream-going source components.
• Broadband cross-spectral matrices when sources are random and partially coherent.

The source may enter a transfer relation as an additive vector:

where represents source forcing. Alternatively, incident modal wave amplitudes may be prescribed on each side of an actuator-disk model. The source representation is then combined with nacelle propagation, reflection, liner attenuation, and radiation models.

An integrated prediction procedure consists of the following stages:

1. Specify operating conditions: Define fan speed, blade and vane counts, mean-flow profiles, temperature, density, and sound speed.
2. Model the source: Determine tonal and broadband modal source amplitudes, including rotor-stator orders
   
3. Determine local modes: Calculate transverse eigenfunctions, cut-off frequencies, axial wavenumbers, and modal impedances for each nacelle section.
4. Discretize the nacelle: Divide the inlet and exhaust paths into uniform, varying-area, lined, and discontinuity elements.
5. Construct element operators: Include propagation, mean flow, liner impedance, mode coupling, area changes, and losses.
6. Assemble the system: Cascade transfer matrices or combine scattering matrices while applying continuity conditions.
7. Apply terminations: Use modal radiation impedances at the inlet lip and exhaust nozzle, together with source-plane conditions.
8. Solve for modal amplitudes: Obtain reflected and transmitted waves throughout the nacelle.
9. Calculate radiation: Convert exit-plane modes into far-field pressure and directivity using radiation or diffraction models.
10. Evaluate metrics: Compute sound pressure level, acoustic power, insertion loss, tone levels, and directivity.
11. Validate and refine: Check modal truncation, element-size convergence, power balance, and agreement with experimental data.

This approach separates source generation, internal transmission, and external radiation while retaining their interaction through reflection and impedance boundary conditions.