Chapter - 11

Key Words: Sound sources, Plane waves, Reflection, Impedance, Absorption, Acoustics

Reflection of Plane Waves

🎯 Learning Objectives

By the end of this topic, you should be able to:

  • Explain and apply Snell’s Law to acoustic wave reflection.
  • Derive and manipulate the reflection factor for oblique incidence.
  • Understand the role of specific acoustic impedance and boundary conditions.
  • Compute the absorption coefficient for locally and non-locally reacting surfaces.
  • Analyze limiting cases such as hard, soft, and matched boundaries.
  • Model plane-wave fields using complex notation and boundary interactions.
  • Apply these results to wall materials, absorbers, underwater surfaces, and architectural acoustics.

Introduction:

A plane acoustic wave interacting with a boundary such as a wall, floor, or interface exhibits reflection, transmission, and sometimes absorption depending on the boundary’s acoustic impedance. This interaction is fundamental in architectural acoustics, noise control, underwater acoustics, and transducer design.
When a plane wave strikes a large, smooth, flat surface, the reflection is specular, meaning that

The angle of incidence equals the angle of reflection - Snell’s Law.

For an incident plane wave approaching at an angle θ\theta, the amplitude and phase of the reflected wave depend on the boundary’s surface impedance.

Geometry of Incidence and Reflection

The angles of the incident and reflected waves satisfy:

θi=θr\theta_i = \theta_r
Signal domain illustration

Figure 1. Free-field propagation from elementary sources in direction of maximum sound radiation

The reflected wave may experience a change in amplitude and a phase shift due to interactions at the boundary.

Mathematical Description of Incident and Reflected Plane Waves:

Incident Wave

A monochromatic plane wave incident on a boundary is expressed as

pi(x,y,t)=p^ej(ωtkxcos(θ)kysin(θ))\underline{p}_i(x,y,t) = \hat{p} e^{j(\omega t - kx \cos(\theta)- ky \sin(\theta))}

where:

  • p^\hat{p} is the complex amplitude,
  • k=ω/ck = \omega/c is the wavenumber,
  • θ\theta is the angle of incidence.

Reflected Wave

The reflection factor for a plane wave at oblique incidence is given by

pr(x,y,t)=p^Rej(ωtkxcos(θ)kysin(θ))\underline{p}_r(x,y,t) = \hat{p} \underline{R} e^{j(\omega t - kx \cos(\theta)- ky \sin(\theta))}

with RR denoting the reflection factor, R=Rejγ R = |\underline{R}| e^{jγ}

It is related to the wall impedance, ZZ, by

R=prpi=Zcos(θ)Z0Zcos(θ)+Z0\underline{R} = \frac{\underline{p}_r}{\underline{p}_i} = \frac{\underline{Z} \cos(\theta)-Z_0}{\underline{Z} \cos(\theta)+Z_0}

where:

  • Z0=ρ0cZ_0 = \rho_0 c is the characteristic impedance of air,
  • Z\underline{Z} is the surface impedance,
  • θ\theta is the angle of incidence.

The wall impedance, ZZ, is defined as the ratio of sound pressure to the normal component of particle velocity, both determined at the wall.c The impedance may be independent of the angle of incidence This phenomenon is called a local reaction. Due to this local reaction adjacent sections of the same wall surface are independent from each other, so that no tangential waves are transmitted along the wall surface. This is a good approximation for heavy walls, for walls with low bending stiffness and for porous absorbers with high flow resistivity.

Example:

Signal domain illustration

Figure 1. Free-field propagation from elementary sources in direction of maximum sound radiation


Python Example: Plane Wave (Oblique Incidence)

2-D visualization of a plane wavefront

Press Run Code: Output will appear here.

Python Example: Reflection From an Impedance Boundary

Press Run Code: Output will appear here.

Absorption Coefficient

The absorption coefficient describes the fraction of incident acoustic energy not returned as reflection:

α=pi2pr2pi2=1R2\alpha = \frac{|\underline{p}_i|^2 - |\underline{p}_r|^2}{|\underline{p}_i|^2} = 1 - |\underline{R}|^2

This is crucial in architectural acoustics and noise control.

Specific Impedance and the Reflection Fact

The specific normalized impedance is:

ζ=1Z0(pvn)=1cosθ1+R1R\zeta = \frac{1}{Z_0} \left( \frac{\underline{p}}{\underline{v}_n} \right) = \frac{1}{\cos\theta} \frac{1 + \underline{R}}{1 - \underline{R}}

For locally reacting surfaces, absorption becomes:

α=4[ζ]cosθ(1+[ζ]cosθ)2+([ζ]cosθ)2\alpha = \frac{4 \Re[\underline{\zeta}] \cos\theta}{(1 + \Re[\underline{\zeta}] \cos\theta)^2 + (\Im[\underline{\zeta}] \cos\theta)^2}

This formulation captures complex impedance surfaces including membranes, porous layers, and micro-perforated materials.


Examples of Wall Impedances

We can define some extreme cases for perpendicular incidence. Real materials approximate these conditions quite well. A heavy concrete wall, for instance, represents a hard wall and the ocean surface represents a soft wall for underwater sound.
Some other examples are listed as follows

  • Matched wall: Z=Z0\underline{Z} = Z_0 ; R=0\underline{R} = 0 ; α=1\alpha = 1
  • Hard wall: Z=\underline{Z} = \infty ; R=1\underline{R} = 1 ; α=0\alpha = 0
  • Soft wall: Z=0\underline{Z} = 0 ; R=1\underline{R} = -1 ; α=0\alpha = 0

Evanescent and Inhomogeneous Waves Near Boundaries:

When the angle of incidence approaches grazing incidence θ90\theta \rightarrow 90^\circ, the component of the acoustic wavenumber normal to the boundary becomes:

kn=kcosθk_n = k \cos\theta

As cosθ0\cos\theta \to 0, the normal wavenumber tends to zero, and beyond this limit becomes purely imaginary.
This yields an evanescent wave, whose pressure field decays exponentially away from the boundary:

pe(y)eknyp_e(y) \propto e^{-|k_n| y}

Evanescent waves:

  • Do not transport energy away from the surface,
  • Can store reactive (non-radiative) energy,
  • Strongly influence near-boundary pressure enhancement,
  • Become significant in micro-perforated absorbers, membrane resonators, and soft surfaces.

In architectural acoustics, they contribute to surface impedance behavior and can modify the effective absorption at extreme incidence angles.

Non-Locally Reacting Surfaces and Surface Wave Propagation:

The assumption of local reaction (impedance independent of incidence angle) fails for:

  • Panels with finite bending stiffness,
  • Multi-layer constructions (e.g., drywall + air cavity),
  • Porous materials with low flow resistivity,
  • Viscoelastic or tensioned membranes.

In such cases, the boundary supports tangential (lateral) wave motion, and the impedance becomes angle-dependent:

Z=Z(ω,θ)\underline{Z} = \underline{Z}(\omega, \theta)

This leads to:

  • Structural–acoustic coupling,
  • Surface-guided waves,
  • Modal conversion between acoustic and structural modes,
  • Frequency-dependent scattering.

A classical example is a thin elastic plate, where bending waves propagate with dispersion:

ω2=Dρshkb4\omega^2 = \frac{D}{\rho_s h} k_b^4

Here, DD is flexural rigidity, ρs\rho_s surface density, and kbk_b bending wavenumber.
Such surfaces no longer satisfy simple reflection formulas and require solving coupled boundary-value problems.

Energy Conservation and Complex Reflection Coefficients:

For general boundaries, the reflection coefficient is complex:

R=Rejγ\underline{R} = |\underline{R}| e^{j\gamma}

It carries:

  • Magnitude → how much energy is returned,
  • Phase → time shift upon reflection.

The conservation law reads:

A+R+T=1A + R + T = 1

where:

  • AA is absorbed energy,
  • R=R2R = |\underline{R}|^2 is reflected energy,
  • TT is transmitted energy (zero for infinitely massive walls).

The phase γ\gamma contains important information:

  • Soft walls exhibit a phase inversion near π\pi
  • Hard walls approximate zero phase shift
  • Reactive surfaces produce frequency-dependent phase changes, affecting standing wave patterns and room response.

Pressure–Velocity Boundary Layer Effects:

Very close to the boundary (within a few millimeters), viscosity and thermal conduction create a boundary layer where:

  • Particle velocity is reduced due to shear,
  • Pressure remains nearly unchanged,
  • Local impedance deviates from its ideal value,
  • Additional losses occur, especially at high frequencies.

The effective acoustic impedance becomes:

Zeff=Zideal+Zviscothermal\underline{Z}_{\text{eff}} = Z_{\text{ideal}} + Z_{\text{viscothermal}}

These corrections are essential for:

  • Micro-perforated panels,
  • High-frequency absorbers,
  • Precision impedance tube measurements.

Surface Impedance of Porous Materials: Delany–Bazley–Miki Models:

Porous materials used in acoustic treatment possess complex, frequency-dependent impedance.
A classical empirical model is:

Z=Z0(1+C1(fσ)C2jC3(fσ)C4)\underline{Z} = Z_0 \left( 1 + C_1 \left( \frac{f}{\sigma} \right)^{-C_2} - j C_3 \left( \frac{f}{\sigma} \right)^{-C_4} \right)

where:

  • σ\sigma is flow resistivity,
  • C1,C2,C3,C4C_1, C_2, C_3, C_4 are empirical constants.

The impedance becomes highly angle-dependent, making the simple reflection factor insufficient.
These models are essential in room acoustics and absorber design.

Complex Pressure Fields and Standing Waves:

The superposition of incident and reflected waves yields:

p(x)=p^(ejkx+Re+jkx)p(x) = \hat{p}(e^{-jkx} + \underline{R} e^{+jkx})

This produces a standing wave ratio (SWR):

S=1+R1RS = \frac{1 + |\underline{R}|}{1 - |\underline{R}|}

Even small reflection magnitudes (R<0.1|\underline{R}| < 0.1) can significantly modify the pressure distribution inside rooms and ducts.
This explains:

  • Bass buildup in corners (high reflection + room modes),
  • Nulls at listening positions (destructive interference),
  • Sensitivity in impedance measurement setups.

Plane Waves vs. Spherical Waves at Finite Distances

The plane-wave reflection theory assumes far-field incidence.
For sources close to the wall:

  • Waves have spherical curvature,
  • Incidence is not uniform across the surface,
  • Near-field terms become significant.

The reflection coefficient becomes position-dependent:

R(r,θ)R(θ)\underline{R}(r,\theta) \neq \underline{R}(\theta)

This is important in:

  • Loudspeaker-wall interactions,
  • Source-to-boundary distance effects,
  • Near-field microphone calibration.

Impedance of Multi-Layer Surfaces (Transfer-Matrix Method):

Real-world walls often consist of multiple layers.
Their combined impedance is computed using the transfer-matrix method (TMM):

[p0v0]=T[pLvL]\begin{bmatrix}p_0 \\v_0\end{bmatrix}=\mathbf{T}\begin{bmatrix}p_L \\v_L\end{bmatrix}

where T\mathbf{T} is the product of layer matrices.
This method captures:

  • Resonances (e.g., mass–spring–mass systems),
  • Phase changes through layers,
  • Frequency-dependent transmission loss.

Such models are essential for wall design, studio isolation, and underwater sonar applications.

Reflection from Curved Surfaces:

For curved boundaries, Snell’s law applies locally, but curvature introduces:

  • Focusing or defocusing of reflected rays,
  • Spatially varying reflection coefficients,
  • Modified standing wave fields,
  • Potential caustics (regions of high intensity).

This becomes important in:

  • Whispering galleries,
  • Curved auditorium walls,
  • Reflective domes and parabolic surfaces.

Numerical Approaches to Plane-Wave Reflection:

In complex geometries, analytical solutions fail.
Numerical methods used include:

  • Boundary Element Method (BEM),
  • Finite Element Method (FEM),
  • Transfer-Matrix Method (TMM),
  • Ray-tracing with frequency-dependent reflection,
  • Parabolic Equation (PE) models (underwater acoustics).

These methods allow modeling:

  • Arbitrary surface geometry,
  • Frequency-dependent impedance,
  • Complex scattering scenarios.

Spherical Wave Above Impedance Plane:

A point source at SS radiates a spherical wave. The total field at RR contains a contribution from the direct sound travelling along the vector r0\vec{r}_0 and another component reflected from the plane.

Signal domain illustration

Figure 1. Spherical wave reflection above impedance plane

A plane wave incidence and a spherical wave incidence differ significantly as the impedance plane is hit at various angles and the total reflection must be integrated:

p=jωρ0Q^ejkr04πr0ω2ρ0Q^4πcθJ0(krsinθ)ejk(z+z0)cosθR(θ)sinθdθp=\frac{j \omega \rho_0 \hat{Q} e^{jkr_0}}{4 \pi r_0}-\frac{\omega^2 \rho_0 \hat{Q}}{4 \pi c}\int_{\theta}{J_0(kr\sin\theta)\, e^{-jk(z+z_0)\cos \theta}\,\underline{R}(\theta)\,\sin\theta\, d\theta}

The first term represents the direct sound and the latter the contribution of the reflection.

Note that the reflection term is an integral over various angles between the wavefront and the surface.

An approximation of the entire expression assuming a constant angle of incidence, θ0\theta_0, is

p=jωρ0Q^ejkr04πr0+jωρ0Q^4πr1R(θ0)p=\frac{j \omega \rho_0 \hat{Q} e^{jkr_0}}{4 \pi r_0}+\frac{j\omega \rho_0 \hat{Q}}{4 \pi r_1}\,\underline{R}(\theta_0)

This equation assumes a constant angle and, thus, a plane wave.
Now, the contribution of the reflection can be related to another point source, called the image source, which is apparently located below the surface radiating a spherical wave whose amplitude is reduced by R(θ0)\underline{R}(\theta_0).

Even though a spherical wave is present, the reflection is calculated for a constant angle θ0\theta_0. The sound field model is effectively switched to a plane-wave description when reflection occurs. This image source model is even accurate when R=1R = 1 or R=1R = -1. For other reflection factors, the approximation is sufficiently accurate when zz, z0z_0, and rλr \gg \lambda.

In a reflecting plane, plane waves or spherical waves can be divided into elementary sources (Huygens principle). Together these sources build up an interference field that produces a plane or spherical wave. This model is the basis for the simplicity of Eq. (which is exact when R=1|\underline{R}| = 1).

If, however, the plane has discontinuities in geometry or impedance, the Huygens superposition is disturbed and the reflected field exhibits scattering and diffraction. Scattering results from reflection by an object; diffraction results from interaction with the edges or boundaries of an object.


🧪 Interactive Examples







📝 Key Takeaways

  • Plane-wave reflection is governed by Snell’s law, ensuring that the angle of incidence equals the angle of reflection for smooth, large surfaces.
  • The reflection factor depends on the surface impedance and incidence angle and determines both the amplitude and phase of the reflected wave.
  • Surface impedance encapsulates the mechanical and acoustic behavior of the boundary and can be real (resistive), imaginary (reactive), or complex.
  • The absorption coefficient quantifies the fraction of acoustic energy not returned as reflection and depends strongly on the impedance and angle of incidence.
  • Locally reacting surfaces have angle-independent impedance, while non-locally reacting surfaces support tangential wave motion and require structural–acoustic models.
  • Evanescent waves arise near grazing incidence and influence near-boundary pressure fields without transporting energy away from the surface.
  • Hard, soft, and matched boundaries represent important limiting cases and guide practical design in architectural, underwater, and engineering acoustics.
  • Advanced models such as Delany–Bazley, transfer-matrix methods, and complex impedance formulations provide realistic predictions for layered or porous materials.

🧠 Quick Quiz

1) What determines the phase shift experienced by a reflected plane wave?

2) For a matched boundary where Z=Z0Z = Z_0, what is the reflection coefficient RR?

3) What is the defining characteristic of a locally reacting surface?

4) Which phenomenon occurs when the normal component of the wavenumber becomes purely imaginary?

5) The absorption coefficient α\alpha is given by 1R21 - |R|^2. What does this physically represent?