# Acoustic Response of a Sinusoidally Perturbed Hard-Walled Duct

Authors: Silvio Del Giudice; Giancarlo Bernasconi
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Acoustic wave propagation in hard-walled ducts is of interest in many fields including vehicle design, musical instruments acoustics, and architectural and environmental noise-control. For the case of small sinusoidal perturbation of the cross-section, it is possible to derive simple though approximate analytical formulas of its plane wave acoustic reflection and transmission spectral response that resembles the optical situation of uniform Bragg gratings. The proof is given here, starting from the “horn equation” and then exploiting the coupled-modes theory. Examples of the results obtained with these analytical formulas are shown for some sinusoidally perturbed ducts and compared to results obtained through a numerical method, revealing a very good agreement.


---

## Body

## 1. Introduction

The propagation of waves in periodic media has received much attention in the past in different fields of physics: a comprehensive review can be found in Elachi [1], with references on the propagation of acoustic waves in ducts with sinusoidally perturbed walls [2–4]. An interesting feature of wave propagation in periodic media discussed in the review is the existence of stopbands and passbands related to the medium periodicities.In the last decades much research effort was dedicated to the theoretical and experimental study of elastic wave propagation in periodic waveguides. For instance, Fokkema [5] dealt with periodic boundaries of elastic media; other authors studied waves propagating along periodically corrugated plates [6–10] and along ducts [11–15].The purpose of this paper is not to advance the research work accomplished so far, but rather to provide an approximate simplification of the established theory, when proper hypotheses are satisfied. Attention is indeed limited to acoustic propagation in hard-walled ducts whose cross-section undergoes a small sinusoidal perturbation with respect to a reference mean value. The work holds for any filling fluid, typically air, provided that the hard-wall hypothesis is verified.Like in Munday et al. [15], the analysis is restricted to one-dimensional (plane wave) propagation, where the waveguide geometry is defined simply by the cross-section along the axial coordinate. The starting point of the theoretical analysis is the Webster horn equation, as performed by Nagarkar and Finch [16], who studied sinusoidal horns, and by Griffiths and Steinke [17], who reviewed the theory of one-dimensional wave propagation in locally periodic media consisting of an arbitrary number of identical cells and showed the acoustic solution for some particular geometries. Lau and Campos [18] also solved the acoustic wave equation for one-dimensional propagation along a duct with a small wall sinusoidal perturbation: the exact solutions were obtained as power series expansions around the middle of the duct.Recently, Hawwa [19] analyzed sound waves in a circular cylindrical duct having a geometric periodicity at its wall, by solving numerically the wave equations.In optics, or more general in electromagnetics, a periodically perturbed medium is called a Bragg grating, or simply a multilayer medium: the transmission and reflection of a uniform grating can be expressed with simple closed-form formulas (Kogelnik [20]). The acoustic analog of the uniform Bragg grating is a duct whose cross-section sinusoidally varies but, to the authors’ knowledge, a simple formula for the acoustic response of such a waveguide is not available in literature.This paper exploits thus the Bragg gratings theory to solve the acoustic “horn equation,” obtaining simple formulas for the reflection and transmission spectral response of the waveguide as a function of acoustic and geometric parameters, under the hypothesis of a small sinusoidal cross-section variation compared to the mean reference value.Simple and closed-form solutions are advantageous in modeling/inversion procedures, as, for example, in bore reconstruction, design of noise-control devices (employed, e.g., for jet engines or HVAC systems), and even monitoring of transportation pipelines.The following sections describe the scenario, the mathematical derivation, and the result for some example cases.The solution provided can take into account also wave attenuation, typically hard-wall losses (boundary layer friction), whose terms are added a posteriori.


## 2. Theory

The scenario considered here is a circular hard-walled duct with a small periodic perturbation of the cross-section along the axial coordinate (Figure1). Two similar though analytically different geometric cases are considered: one is when a sinusoidal function describes the variation of the cross-section; the other when a sinusoidal function describes the variation of the radius of the circular cross-section.Figure 1
Duct with a periodic perturbation of the cross-section.As in Lau and Campos [18], the acoustic wavelength is supposed to be larger than the transverse dimensions of the duct, so that only the fundamental longitudinal mode exists. Moreover, the changes in cross-section are supposed not to be too rapid with respect to the transverse dimension, so that the wavefronts remain approximately plane.Under these hypotheses the governing geometrical parameters reduce to the cross-section area (the shape can be neglected), and the starting physical law is the one-dimensional “horn equation” [17]:
(1)∂2Ψ∂t2=v2(∂2Ψ∂x2+1SdSdxdΨdx).
The field variable Ψ is the pressure over ambient, S is the cross-section area, v is the phase velocity, x is the waveguide axial coordinate, and t is the time.Let us consider harmonic waves:(2)Ψ=ψ(x)e-iωt.The substitution of (2) in the wave equation (1) yields
(3)∂2ψ∂x2+1SdSdxdψdx+β2ψ=0,
where β=ω/v is the acoustic wavenumber.From this point, the mathematical steps and the approximations introduced follow the computation of the transfer function for a uniform optical fiber Bragg grating (Erdogan [21] and Kogelnik [20]), with adaptation to the acoustic waveguide case.The pressure field is expressed as a linear combination of the fundamental modes propagating in the opposite directions(4)ψ(x)=A(x)e-iβx+B(x)eiβx,
where, according to the adopted time convention, A and B are the amplitudes of the waves propagating in the -x and +x directions, respectively.The first derivative of (4) is
(5)ψ′=A′e-iβx-iβAe-iβx+B′eiβx+iβBeiβx,
with prime denoting derivative with respect to x.In the hypothesis of weak coupling between the two modes, the second derivatives ofA and B are neglected, because A′′≪β2A and B′′≪β2B; therefore, the second derivative of (4) becomes
(6)ψ′′≅-2iβA′e-iβx-β2Ae-iβx+2iβB′eiβx-β2Beiβx.Substituting the expressions for derivatives, and multiplying byeiβx, the wave equation (3) yields
(7)S′S(A′-iβA+B′e2iβx+iβBe2iβx)-2iβA′+2iβB′e2iβx=0.
The solution of (7) requires a system of two differential equations for A and B. Therefore A′ and B′ are alternately isolated from (7), neglecting in each resulting equation the dependency on the derivative of the other coefficient. This is justified by the fact that  S and so A and B are slowly variable functions of  x [21].The resulting system is(8)A′=iβS′/SS′/S-2iβA+-iβS′/SS′/S-2iβBe2iβx,B′=-iβS′/SS′/S+2iβB+iβS′/SS′/S+2iβAe-2iβx,
that can be expressed in this way:
(9)A′=-ic11A-ic12Be2iβx,B′=-ic22B-ic21Ae-2iβx.
Under the hypothesis of slow section perturbations, and far from the null frequency, it is |S′/S|≪|2iβ|, and the coefficients cij in (9) become
(10)c11=-βS′/SS′/S-2iβ≅S′2iS,c12=βS′/SS′/S-2iβ≅-S′2iS,c21=-βS′/SS′/S+2iβ≅-S′2iS,c22=βS′/SS′/S+2iβ≅S′2iS.
Now, the two aforementioned geometric cases are considered; first, when the cross-section is a cosine function of  x; that is,
(11)S(x)=scos(γx)+s0,
and second, when the radius of the circular cross-section is a cosine function of x; that is,
(12)S(x)=π(scos(γx)+s0)2.γ is the perturbation wavenumber,  s is the perturbation amplitude, and  s0 is the mean value.The ratio between the derivative of the cross-section and the cross-section itselfS′/S, for s≪s0, is approximately
(13)S′(x)S(x)≅-4κsin(γx)=2iκ(e+iγx-e-iγx),
where κ=sγ/4s0for the cross-section sinusoidal perturbation and κ=sγ/2s0for the radius sinusoidal perturbation.The system (9) becomes
(14)A′=-iκ(e+iγx-e-iγx)A+iκ(e+iγx-e-iγx)Be2iβx,B′=-iκ(e+iγx-e-iγx)B+iκ(e+iγx-e-iγx)Ae-2iβx.
The terms that contain a rapidly oscillating dependence with x can be neglected: they correspond to the complex exponentials with high phase constant (in magnitude) compared to the others, that is, e±iγx, e2iβx+iγx, and e-2iβx-iγx [21].The result is(15)A′=-iκBei(2β-γ)x,B′=iκAe-i(2β-γ)x.
By performing in the system (15) the following substitutions:
(16)a=Ae-i(β-γ/2)x,b=Bei(β-γ/2)x,
and calling σ=β-γ/2, one obtains
(17)a′=-iσa-iκbb′=iσb+iκa.
The system (17) is a standard system of two coupled first-order ordinary differential equations with constant coefficients, for which closed-form solutions can be found, when appropriate boundary conditions are specified.SinceA corresponds to the backward propagating wave and B to the forward propagating one, the boundary conditions, for a forward propagating wave impinging the waveguide at  x=0, are
(18)A(L)=0B(0)=1⟹a(L)=0b(0)=1,
where L is the length of the sinusoidally perturbed duct.The solution of system (17) with boundary conditions (18) is obtained by means of linear algebra:
(19)a(x)=iκ(e-δLeδx-eδLe-δx)eδL(iσ-δ)-e-δL(iσ+δ),b(x)=eδLe-δx(iσ-δ)-e-δLeδx(iσ+δ)eδL(iσ-δ)-e-δL(iσ+δ).The reflection and transmission spectral responsesR(ω) and T(ω) correspond, respectively, to a(0) and b(L). Hence, finally
(20)R(ω)=iSsinhδLcoshδL-iPsinhδL,T(ω)=1coshδL-iPsinhδL,
where
(21)δ=κ2-σ2,S=κδ,P=σδ.According to their definition, the spectral responsesR and T refer to the acoustic pressure amplitude: the corresponding power spectral responses can be found by taking the square of their magnitude. In absence of attenuation, the sum of the power spectral responses is unitary, as expected.Equations (20) are the acoustic analog of the optical formulas for the uniform Bragg grating in Kogelnik [20], and the parameters κ and σ are expressed here as a function of acoustic wave parameters and waveguide geometric parameters.The reflection spectral response has a maximum forσ=0, corresponding to the frequency fmax=γv/4π. For the same frequency the transmission response has a minimum.Finally, even if the considered one-dimensional wave equation does not contain a loss term, one may add the absorption phenomenon a posteriori, typically due to the duct walls, by redefining the acoustic wavenumber(22)β=ωv+iα,
where α(ω) is the absorption coefficient computed with mean cross-section parameters and v(ω) is the corresponding phase velocity.


## 3. Examples

Some examples of the acoustic behavior for a sinusoidally perturbed duct, filled with air at standard conditions (20°C, 1 atm), are presented here. They reveal the link between the passbands/stopbands response and the waveguide geometrical parameters, and they permit inferring the validity limit of the approximate formula.Wave attenuation and dispersion are included too, using (22), and α(ω) and v(ω) are computed according to the wide-tube approximation (well described by Tijdeman [22]). Closed-form solutions are compared with full waveform numerical methods.We consider a sinusoidal duct with lengthL, parameterized by the radius of the circular cross-section r:
(23)r(x)=scos(γx)+s0,0<x<L.
The average radius s0 is set to 2 cm, and the radius perturbation frequency γ/2π is set to 50m-1.The other parameters are set to the following values:Case (i):L=2m, s=0.02mm;Case (ii):L=20m, s=0.2mm;Case (iii):L=2m, s=2mm.Figure2 shows the reflection and transmission responses, in magnitude, computed with (20), for the different cases, in a frequency interval centered at the frequency fmax. The wavelength corresponding to fmax is not much higher than the duct transverse dimension, as required.Reflection ((a), (c), (e)) and transmission ((b), (d), (f)) spectral responses, in magnitude, of a sinusoidally perturbed duct computed with (20) (solid line) and numerical results (crosses). ((a)-(b)) Case (i); ((c)-(d)) Case (ii); ((e)-(f)) Case (iii).
(a)
(b)
(c)
(d)
(e)
(f)The crosses in Figure2 are computed with an “exact” one-dimensional simulator, based on the computation method described by Munday et al. [15], but with transmission coefficients included: the agreement between the approximate and the exact solution is very good.The parameters values have been chosen to progressively increase thes/s0 ratio, since a hypothesis for the approximation is s≪s0: the effect of increasing this ratio is to widen the stopband. Moreover, Case (ii) investigates the effect of extending the duct length L, which is the other geometric parameter of the periodic structure: a higher length results in a sharper transition of the stopband.The difference between approximate and theoretical values is barely noticeable even whens/s0 ratio is as large as 1/10, which may be considered the limit of validity of this approximation. Furthermore, the analytic formula is able to correctly simulate the attenuation phenomenon (this is apparent in the transmission response) which was not theoretically justified, but simply added a posteriori.The phases are not shown here, but they have been verified as well.


## 4. Conclusion

The acoustic spectral response of a sinusoidally perturbed hard-wall duct has been derived and given in a simple formula, by following the optical analog of Bragg gratings. The formula is the same as in optics, with the electromagnetic parameters replaced by their equivalent acoustic parameters and periodic duct geometry.Results are valid for small cross-section perturbations and in this case successful comparisons with a numerical method are shown, even in case of wave attenuation.The availability of simple analytical formulas permits a direct analysis of the link between the acoustic response and the duct geometrical parameters and the design of efficient modeling/inversion procedures in the fields of bore/pipe reconstruction, noise control, and so forth.Finally it can be noticed that, even if the derivation strictly requires a sinusoidal perturbation, any cross-section deformation can be decomposed in sinusoidal functions and therefore the results can be applied, provided that the underlying hypotheses are satisfied, to a much broader range of scenarios.



---

*Source: 267291-2013-11-12.xml*