# Analysis and Optimal Condition of the Rear-Sound-Aided Control Source in Active Noise Control

Authors: Karel Kreuter; Yasuhide Kobayashi
Journal: Advances in Acoustics and Vibration
 (2011)

---

## Abstract

An active noise control scenario of simple ducts is considered. The previously suggested technique of using an single loudspeaker and its rear sound to cancel the upstream sound is further examined and compared to the bidirectional solution in order to give theoretical proof of its advantage. Firstly, a model with a new approach for taking damping effects into account is derived based on the electrical transmission line theory. By comparison with the old model, the new approach is validated, and occurring differences are discussed. Moreover, a numerical application with the consideration of damping is implemented for confirmation. The influence of the rear sound strength on the feedback-path system is investigated, and the optimal condition is determined. Finally, it is proven that the proposed source has an advantage of an extended phase lag and a time delay in the feedback-path system by both frequency-response analysis and numerical calculation of the time response.


---

## Body

## 1. Introduction

In the recent past the concept of active noise control (ANC) has arrest more attention. Due to the falling prices for digital signal processors (DSPs), ANC applications become more popular everyday. The Swinbanks’ (unidirectional) source [1] has frequently been applied to control noise scenarios to show its advantages over a standard bidirectional case with one control source. The implementation has shown to be of lower cost and at the same time better control performance. This has been confirmed by experiments in both adaptive and robust control setups [2, 3]. A welcome side effect was a lower controller gain and driving signals for the used loudspeakers [3, 4], which enables the usage of lower power loudspeakers. The advantage of the Swinbanks’ source has been theoretically proven for a practical setting on the duct length and loudspeaker locations [5, 6].The intention of the Swinbanks’ source is cancelling of the upstream sound, which is emitted by the downstream loudspeaker. In order to accomplish this, the upstream source is driven with the same amplitude as the downstream loudspeaker but opposite phase and proper time delay. Thus, the upstream sound at the mount of the upstream source is eliminated. The advantage in the control performance can be explained by the extension of the time period by which the sound emitted by the control source travels to the reference microphone. It has been proven that the Swinbanks’ source introduces an additional time delay corresponding with twice the distance between control source and duct end [5, 6]. Previous research pointed out on empirical basis that the closed-loop performance is improved when the actuator (control source) and the sensor (reference microphone) are of further distance [7, 8].In addition to the Swinbanks’ source, the rear-sound-aided source has been proposed to achieve similar advantage of Swinbanks’ source by using one loudspeaker [9]: the source is composed of a loudspeaker and a subduct which connects the rear side of the loudspeaker to an upstream junction in the main duct in order to attenuate the upstream sound by interfering the rear sound as shown in Figure 1. Experiments demonstrated that the proposed source has an advantage. It achieves a lower amplitude of error microphone signal by a smaller amplitude of control input over the conventional bidirectional source.Figure 1
Experimental apparatus used in [9].However, the following problems remain for the proposed source.(i)
The first principle model derived in [9] has not been well consistent with the experimental result of the frequency response.(ii)
No advantage of the proposed source has been shown theoretically: the additional time delay in the feedback-path transfer function has not been analyzed.In this paper, we concentrate on the second problem, because the novel contribution of showing an advantage in a theoretical way is also necessary to solve for the first problem in the future.In order to show the advantage of the proposed source, feasibility of the additional phase lag in the frequency response and the additional time delay in the time response are both investigated in this paper based on the frequency response analysis by the first principle model and the time-response simulation by the numerical calculation. In addition, the consistency between the frequency response analysis and the time response simulation will be shown to validate an assumption involved in the simulation model due to the junction of the proposed source, which is not necessary for more simpler structures of the conventional sources. To show consistency, the frequency response is numerically calculated by injecting a sinusoidal sound, where the damped wave equation is used to consider energy dissipation. Herewith, the steady-state response does not diverge for resonance frequencies.This paper is outlined as follows. Firstly, the first principle model [9] is derived based on the damped wave equation which later enables comparison of the frequency response of the first principle model with the numerical simulation. In Section 3, the damped wave equation will be directly implemented in a numerical simulation by using finite differences. We will demonstrate that the frequency response calculated by numerical implementation is consistent with the one given by the first principle model. Thus, the model will be available to examine both the additional phase lag and time delay in the feedback-path transfer function. Finally, the optimal condition achieving the longest time delay on the feedback-path transfer function will be given based on the first principle model in Section 4. The condition will be used in numerical simulations to directly visualize the additional time delay.


## 2. First Principle Model

In this section, the damped wave equation is used to derive frequency response functions in order to obtain a finite amplitude response at resonance frequencies. The damped wave equation is the wave equation with an additional term of the first partial derivative of pressure in time (see, e.g., [10]). The differences to the existing result of [9] will be explained, since the existing result is based on the (undamped) wave equation. Note that the derivation process of the damped wave equation will be described in this section, since the intermediate equations ((1) and (2) below) are necessary to exploit the existing result of the transmission line theory.The damped wave equation can be derived by two physical laws, the conservation of mass (1) and momentum (2). The acoustic velocity is given by u, p represents the pressure. Note that an additional term du is included in the conservation of momentum to take friction effects into account, where d is a damping coefficient,(1)1c02∂p∂t+ρ0∂u∂x=0,(2)ρ0∂u∂t+∂p∂x+du=0.
If the damping is neglected, the two equations are well known in acoustic literature (see, e.g., [11]) and are used to derive the wave equation without consideration of damping. However, to the best of our knowledge, (2) with d≠0 has not been used to derive a physical model in the field of active noise control. Usually, the damped wave equation (3) is solved without considering damping. The resulting undamped solution is used to approximate damping effect by modifying the wave number k:=ω/c0 (real number) to k:=(ω/c0)+jα (complex number) [12], which technique has also been used in [9]. The damped wave equation can be obtained by the derivation of (1) with respect to t and (2) with respect to x and combining both(3)1c02∂2p∂t2+dρ0c02∂p∂t=∂2p∂x2.
On the other hand, more general partial differential equations and the solution with damping consideration are common in the transmission line theory. It is well known for instantaneous voltage v(x,t) and current i(x,t) at a location x on a lossy transmission line that the following relations hold:(4)∂v∂x=-Ri-Li∂i∂t,∂i∂x=-Gv-C∂v∂t,
where R, L, C, and G are resistance, inductance, capacitance, and conductance per unit length of the line. One can verify that (4) equals (1) and (2) when using the following relations:(5)v(x,t)=p(x,t),i(x,t)=u(x,t),R:=d,L:=ρ0,G:=0,C:=1ρ0c02.
For harmonic signals(6)v(x,t):=Re{V(x)ejωt},i(x,t):=Re{I(x)ejωt},
the partial differential equation of transmission line theory can be rewritten as(7)dVdx=-ZI,dIdx=-YV,
where Z:=R+jωL and Y:=G+jωC are impedance and admittance per unit length. The general solution of (7) is known as(8)V=Ae-γx+Beγx,I=YZ(Ae-γx+Beγx),
where A and B are complex-valued integral constants. γ is given by(9)γ:=ZY:=α+jβ,α:=12R2+ω2L2G2+ω2C2-(ω2LC-RG),β:=12R2+ω2L2G2+ω2C2+(ω2LC-RG).Lemma 1.
In transmission line, suppose thatv(x,t) and i(x,t) are in steady state with the relations given in (6). Then, V(x) and I(x) at arbitrary position x is related with the values at origin, V(0) and I(0), by
(10)[V(x)I(x)]=[cosh⁡⁡(γx)-Z0sinh⁡⁡(γx)-1Z0sinh⁡⁡(γx)cosh⁡⁡(γx)][V(0)I(0)],
where Z0:=Z/Y is the characteristic impedance of the transmission line.Proof.
By substitutingx=0 into (8), A and B are written in V(0) and I(0). Then, the elimination of the constants A and B leads to (10).Likewise, a similar lemma holds for the sound propagation in a straight duct.Lemma 2.
Suppose thatp(x,t) and u(x,t) are in steady state with the relation given as
(11)p(x,t):=Re{p(x)ejωt},u(x,t):=Re{u(x)ejωt}.
Then,p(x) and u(x) at arbitrary position x is related with the values at origin, p(0) and u(0), by
(12)[p(x)u(x)]=[cosh⁡⁡(γx)-Z0sinh⁡⁡(γx)-1Z0sinh⁡⁡(γx)cosh⁡⁡(γx)][p(0)u(0)],(13)Z0∶=ρ0c01+djωρ0,γ:=α+jβ,(14)α∶=12c0ωω2+d2ρ02-ω2,β∶=12c0ωω2+d2ρ02+ω2,
where Z0 is called the acoustic characteristic impedance of the duct.Proof.
It is trivial from Lemma1 and the relations (5).The next lemma is used to derive first principle model so that series connections of ducts are simplified in the derivation process, which is an extended result to the undamped case in [9].Lemma 3.
Define a matrixT(l) with arbitrary complex numbers γ and Z0 by
(15)T(l)∶=[cosh⁡⁡(γl)-Z0sinh⁡⁡(γl)-1Z0sinh⁡⁡(γl)cosh⁡⁡(γl)].
Then, for any positive real number x and y, the following holds:
(16)T(x)T(y)=T(x+y).Proof.
It can be verified by multiplication and application of the theorems forsinh and cosh.The relationship of pressurep and acoustic velocity u between different points in space will serve as the basis for deriving a unifying transfer function which includes the proposed control source. We consider the duct system in Figure 2 which is the same as in [9], where SPK1 is a noise source, SPK2 and SPK3 are used as a control source, controller is connected between y and u to drive the control source based on the information of the reference microphone output y, and the error microphone output z is used to evaluate control performance. Suppose that SPK1, Ref.Mic., SPK3, SPK2, and Err.Mic. are located at x=0,ly,lv,lu, and lz, respectively, in the main duct of length L.Figure 2
Theoretical apparatus.There are three scenarios for the control source.(i)
Case A—Bidirectional source: SPK3 is not used; that is,v(t)=0 or H=0 and no subduct installed (LS=0). This can be described as bidirectional case since SPK2 emits waves in both directions.(ii)
Case B—Swinbanks’ source: SPK3 is used to cancel the upstream sound produced by SPK2 without subduct (LS=0). It is, therefore, driven with the opposite signal, but delayed v(t)=-u(t-(lu-lv)/c0). This can be represented by the transfer function H(s)=-es((lu-lv)/c0).(iii)
Case C—rear sound aided source (proposed one): The rear sound of SPK2 is redirected by a subduct and attached in a way that the rear sound will interfere with the front sound at the junction of ducts. Note here that the assumption is made, that the rear sound travelling in the subduct can be simulated by a subduct of lengthLS=lu-lv with an additional loudspeaker (SPK3) to cover the rear sound. The amplitude of SPK3 can be adjusted by the variable H; that is, u(t)=H·v(t), where H is treated as negative real number, whose magnitude is less than or equal to 1, in this paper. (H is set to −1 when the strength of the rear sound is as the same to the front one.) The realizability of this setting is under consideration, however, one possibility of this is to cover the inside wall of the subduct by glass wool to weaken the rear sound.In Figure2, (p•,u•) represents the pressure and the particle velocity at each position in the duct. For the modelling process, the input signals w, u, and v are assumed to be u0, uS and ũS, respectively, neglecting the loudspeaker dynamics. The outputs z and y are assumed to be pz and py, respectively, neglecting the microphone dynamics. Thus, u0:=w, uS:=u, ũS:=v, pz:=z, py:=y.According to Lemma2, the derivation process of the first principle model is eventually similar to the undamped case [9], where the derivation of the transfer function has been explained in detail. Hence, only the resultant transfer functions are shown as follows: (17)G0u=Z0(sinh⁡⁡(γ(L-lu))cosh⁡⁡(γLs)+sinh⁡⁡(γ(L-lv))⋅H)sinh⁡⁡(γ(L-lv))sinh⁡⁡(γLs)cosh⁡⁡(γlv)+cosh⁡⁡(γL)cosh⁡⁡(γLs),Gyu=G0u⋅cosh⁡⁡(γly),Gzu=-Z0(sinh⁡⁡(γ(lz-lu))+sinh⁡⁡(γ(lz-lv))cosh⁡⁡(γLs)⋅H)+G0u(sinh⁡⁡(γ(lz-lv))tanh⁡(γLs)cosh⁡⁡(γlv)+cosh⁡⁡(γlz)),G0w=Z0(sinh⁡⁡(γ(L-lv))tanh⁡(γLs)sinh⁡⁡(γlv)+sinh⁡⁡(γL))sinh⁡⁡(γ(L-lv))tanh⁡(γLs)cosh⁡⁡(γlv)+cosh⁡⁡(γL),Gyw=cosh⁡⁡(γly)G0w-Z0sinh⁡⁡(γly),Gzw=(sinh⁡⁡(γ(lz-lv))tanh⁡(γLs)cosh⁡⁡(γlv)+cosh⁡⁡(γlz))G0w-Z0(sinh⁡⁡(γ(lz-lv))tanh⁡(γls)sinh⁡⁡(γlv)+sinh⁡⁡(γlz)). It is obvious that for undamped case the derived functions and the ones given by [9] ought to be the same. This can be confirmed by setting the damping coefficient d=0 and using relationships between hyperbolic and trigonometric functions. Indeed, consistency has been confirmed for the undamped case. The major difference for the damped case is that the damping and acoustic impedance Z0 considered in [9] are constant for all frequencies, while those in the derived functions above are frequency dependent: the imaginary part of the complex wave number in [9], α of k=(ω/c0)-jα, does not depend on the angular frequency ω, while the one in the corresponding parameter -jγ=β-jα depends on ω. Specifically, the difference mainly arises in low frequency range, since α given by (14) might be regarded as a constant number in high-frequency range, whose limit is given by limω→∞α=d/(2ρ0c0). Figure 3 shows the frequency response of Case A (only Gzw is shown in the figure due to the limitation of the space) for the both considerations of damping, where the same parameters in [9] are used as follows: (18)L=3.61,lz=3.53,lu=2.32,ly=0.03,c0=340,ρ0=1.21.
The damping coefficient is set as d=72.6 which is obtained from the relationship d/(2ρ0c0)=30/340≈0.09 (0.09 is given in [9]) so that both considerations of damping result similar damping effect in high-frequency range. Although overall bode characteristics did not change, some minor differences appear due to the different approach of taking damping into account in low frequencies as it is shown in Figure 3(b).Comparison ofGzw (Case A) with [9].
(a)
Overview(b)
Zoom on low frequenciesWe here emphasize that the derived model in this paper is more relevant than that of [9], since the frequency response of the model can be directly compared with the one calculated by numerical implementation in the next section.


## 3. Numerical Implementation

The solution of the damped wave equation (3) using the transmission line theory shall be confirmed by comparing to an alternative approach, that is, numerical simulation by implementing a finite difference method, in which forward and backward difference are used to calculate the frequency response of the system.In this section, the numerical algorithm for Case A will first be derived. Based on this, the numerical algorithm for Case C will be explained.


### 3.1. Case A

The effect of a loudspeaker in duct can be considered in the wave equation by using a spacial delta function [8]. Thus, the damped wave equation with a loudspeaker can be described as follows:(19)1c02∂2p∂t2+dρ0c02∂p∂t=∂2p∂x2+ρ0u̇δ(x-lu).
For numerical analysis, time t and space x are discretized as follows:(20)t=ti∶=i⋅Δt,∀i∈{0,1,…,N},x=xj∶=j⋅Δx,∀j∈{0,1,…,M},N∶=TΔt,M∶=LΔx,ju∶=luΔx,
where T is a simulation period. After applying forward and backward difference on the damped wave equation (19), one obtains(21)p(ti+1,xj)-2p(ti,xj)+p(ti-1,xj)Δt2+dρ0p(ti,xj)-p(ti-1,xj)Δt=c02p(ti,xj+1)-2p(ti,xj)+p(ti,xj-1)Δx2+c02ρ0u̇δjjuΔx,
where δjju denotes the Kronecker delta to approximate the original property of integration being 1 by summation as follows:(22)1=∫0Lδ(x-lu)dx≈∑j=0M-1(δjjuΔx)Δx,δij={1,(i=j),0(i≠j).
Rearranging the recurrence relation for future time value and introducing auxiliary variables leads to(23)p(ti+1,xj)=p(ti,xj)⋅[2-s-2r]+p(ti-1,xj)⋅[s-1]+p(ti,xj+1)⋅[r]+p(ti,xj-1)⋅[r]+rΔxρ0u̇δjju,(24)s=Δt⋅dρ0,r=(c0⋅ΔtΔx)2.
Since the duct has two ends, boundary conditions have to be set accordingly. As it has been done in the analytical process, the pressure at the open end is set to zero; that is, p(ti,xM)=0 for all i. For the closed end at x=0, the pressure gradient has to be zero. Thus, the relation(25)∂p(t,x)∂x|x=0≈p(ti,x0)-p(ti,x-1)Δx=0⟺p(ti,x-1)=p(ti,x0)
holds with a pressure at a fictitious location x=-Δx=:x-1. This relation can be used in (23) for j=0 leading to(26)p(ti+1,x0)=p(ti,x0)⋅[2-s-r]+p(ti-1,x0)⋅[s-1]+p(ti,x1)⋅[r].By using the relations above, frequency response of the feedback-path systemGyu can be numerically calculated by injecting the source input as a sinusoidal signal u(ti)=-Ucosωti:  (1)
initial condition:p(t0,xj)=p(t1,xj)=0 for all j,(2)
fori=1,2,…,N-1 repeat the following:
(a)
at the open end:p(ti,x0)=0,(b)
at the closed end: (26),(c)
forj=1,2,…,M-1 repeat (23). (only at the source position j=ju, the last term rΔxρ0u̇(ti)δjju appears as rΔxρ0Usinti.)For the stability of the numerical algorithm above, parameters in Table1 are set. The steady state pressure variation at the microphone position p(ti,xjy) is evaluated to get gain and phase from the source input u to the microphone output y, where jy is defined as jy:=ly/Δx. The other frequency responses for Gzu, Gzw, and Gyw can be calculated in a similar way. For each transfer function 2000 logarithmically distributed frequencies have been simulated.Table 1
Numerical approximation settings.
ParameterNotationValueUnitSpatial step widthΔx0.01mSimulation timeT1sTime stepsN40000


### 3.2. Case C

As Case C consists of the main duct and an additional subduct, the previous derivation has to be adapted accordingly. First, similar to the main duct, an additionalz axis is introduced so that SPK3 is located at z=0 and the subduct is connected at z=LS as shown in Figure 2, by which the instantaneous pressure in the subduct is denoted as pS(t,z). Then, t and z are discretized in a similar manner as Case A(27)z=zk∶=k⋅Δx∀k∈{0,1,…,MS},MS∶=LSΔx,
where common spacing Δx is used. The input for SPK3 is specified as v(ti)=Hu(ti). The boundary condition for the closed end of the attached subduct is added like before. Furthermore, the pressure at the junction is assumed as the average of all surrounding pressure points(28)p(ti+1,xjv)=pS(ti+1,zMS)=p(ti+1,xjv-1)+p(ti+1,xjv+1)+pS(ti+1,zMS-1)3,
where jv:=lv/Δx. The relation above might now be only an assumption; however, the validity will be shown later. The numerical algorithm in detail is omitted from the paper, since it is parallel to Case A based on the additional setting for Case C. So far, two solutions of the frequency response based on the damped wave equation (3) have been obtained. In the remaining of this section, the results will be compared. Figure 4 shows frequency response of the feedback-path system for Cases A and C.Comparison of numerical and analytical solution, transfer functionGyu, H=-1.
(a)
Case A(b)
Case CApparently, both solutions match almost perfect as it can be seen in above figures. This testifies the application of transmission line theory and numerical simulation via finite difference method with the applied assumptions. However, the desired additional phase lag of Case C compared to Case A did not occur. In both cases, the total phase lag is about 2520 degree. A possible cause for the missing phase lag could be a different setting of the rear-sound strengthH, which leads to different bode diagrams for Gyu as pointed in [9].


## 3.1. Case A

The effect of a loudspeaker in duct can be considered in the wave equation by using a spacial delta function [8]. Thus, the damped wave equation with a loudspeaker can be described as follows:(19)1c02∂2p∂t2+dρ0c02∂p∂t=∂2p∂x2+ρ0u̇δ(x-lu).
For numerical analysis, time t and space x are discretized as follows:(20)t=ti∶=i⋅Δt,∀i∈{0,1,…,N},x=xj∶=j⋅Δx,∀j∈{0,1,…,M},N∶=TΔt,M∶=LΔx,ju∶=luΔx,
where T is a simulation period. After applying forward and backward difference on the damped wave equation (19), one obtains(21)p(ti+1,xj)-2p(ti,xj)+p(ti-1,xj)Δt2+dρ0p(ti,xj)-p(ti-1,xj)Δt=c02p(ti,xj+1)-2p(ti,xj)+p(ti,xj-1)Δx2+c02ρ0u̇δjjuΔx,
where δjju denotes the Kronecker delta to approximate the original property of integration being 1 by summation as follows:(22)1=∫0Lδ(x-lu)dx≈∑j=0M-1(δjjuΔx)Δx,δij={1,(i=j),0(i≠j).
Rearranging the recurrence relation for future time value and introducing auxiliary variables leads to(23)p(ti+1,xj)=p(ti,xj)⋅[2-s-2r]+p(ti-1,xj)⋅[s-1]+p(ti,xj+1)⋅[r]+p(ti,xj-1)⋅[r]+rΔxρ0u̇δjju,(24)s=Δt⋅dρ0,r=(c0⋅ΔtΔx)2.
Since the duct has two ends, boundary conditions have to be set accordingly. As it has been done in the analytical process, the pressure at the open end is set to zero; that is, p(ti,xM)=0 for all i. For the closed end at x=0, the pressure gradient has to be zero. Thus, the relation(25)∂p(t,x)∂x|x=0≈p(ti,x0)-p(ti,x-1)Δx=0⟺p(ti,x-1)=p(ti,x0)
holds with a pressure at a fictitious location x=-Δx=:x-1. This relation can be used in (23) for j=0 leading to(26)p(ti+1,x0)=p(ti,x0)⋅[2-s-r]+p(ti-1,x0)⋅[s-1]+p(ti,x1)⋅[r].By using the relations above, frequency response of the feedback-path systemGyu can be numerically calculated by injecting the source input as a sinusoidal signal u(ti)=-Ucosωti:  (1)
initial condition:p(t0,xj)=p(t1,xj)=0 for all j,(2)
fori=1,2,…,N-1 repeat the following:
(a)
at the open end:p(ti,x0)=0,(b)
at the closed end: (26),(c)
forj=1,2,…,M-1 repeat (23). (only at the source position j=ju, the last term rΔxρ0u̇(ti)δjju appears as rΔxρ0Usinti.)For the stability of the numerical algorithm above, parameters in Table1 are set. The steady state pressure variation at the microphone position p(ti,xjy) is evaluated to get gain and phase from the source input u to the microphone output y, where jy is defined as jy:=ly/Δx. The other frequency responses for Gzu, Gzw, and Gyw can be calculated in a similar way. For each transfer function 2000 logarithmically distributed frequencies have been simulated.Table 1
Numerical approximation settings.
ParameterNotationValueUnitSpatial step widthΔx0.01mSimulation timeT1sTime stepsN40000


## 3.2. Case C

As Case C consists of the main duct and an additional subduct, the previous derivation has to be adapted accordingly. First, similar to the main duct, an additionalz axis is introduced so that SPK3 is located at z=0 and the subduct is connected at z=LS as shown in Figure 2, by which the instantaneous pressure in the subduct is denoted as pS(t,z). Then, t and z are discretized in a similar manner as Case A(27)z=zk∶=k⋅Δx∀k∈{0,1,…,MS},MS∶=LSΔx,
where common spacing Δx is used. The input for SPK3 is specified as v(ti)=Hu(ti). The boundary condition for the closed end of the attached subduct is added like before. Furthermore, the pressure at the junction is assumed as the average of all surrounding pressure points(28)p(ti+1,xjv)=pS(ti+1,zMS)=p(ti+1,xjv-1)+p(ti+1,xjv+1)+pS(ti+1,zMS-1)3,
where jv:=lv/Δx. The relation above might now be only an assumption; however, the validity will be shown later. The numerical algorithm in detail is omitted from the paper, since it is parallel to Case A based on the additional setting for Case C. So far, two solutions of the frequency response based on the damped wave equation (3) have been obtained. In the remaining of this section, the results will be compared. Figure 4 shows frequency response of the feedback-path system for Cases A and C.Comparison of numerical and analytical solution, transfer functionGyu, H=-1.
(a)
Case A(b)
Case CApparently, both solutions match almost perfect as it can be seen in above figures. This testifies the application of transmission line theory and numerical simulation via finite difference method with the applied assumptions. However, the desired additional phase lag of Case C compared to Case A did not occur. In both cases, the total phase lag is about 2520 degree. A possible cause for the missing phase lag could be a different setting of the rear-sound strengthH, which leads to different bode diagrams for Gyu as pointed in [9].


## 4. Optimal Condition for the Longest Time Delay in Feedback-Path System

In order to examine the effect of a different setting of the rear-sound strengthH, the analytic solution is modified and compared to itself with several settings for H, as it can be seen in Figure 5.Figure 5
Comparison of transfer functionGyu of Case C for different rear-sound strengths H.Evidently, a lower value ofH will lead to the desired effect and, thus, shows an advantage over the bidirectional scenario. The reason for that behaviour can be found in the transfer function for G0u, since it is part of Gyu. For the explanation, the undamped equation will be used as it simplifies the next steps and did not show much influence on phase characteristics. The derivation will be shown for the transfer function of Case A, the same method can later be applied to Case C. Firstly, G0u will be simplified for Case A by setting Ls=0, H=0 and rewriting as exponential functions. After making the equation causal by multiplying with the longest time delay e-jkL to both denominator and numerator, and replacing jk by s/c0, the following is obtained:(29)G0u=ρ0c0e-s(lu/c0)-e-s(2L-lu)/c01+e-s(2L/c0).
This transfer function corresponds with a causal system depicted in Figure 6, where each signal block represents a time delay.Figure 6
Signal flow chart in frequency domain.It can be observed that the shortest time delay is given bye-s(lu/c0). As mentioned in the introduction, the value of the shortest time delay is critical for the control performance. Similar for Case C, the same procedure is applied and some variables for better overview are introduced. Since the denominator is not of relevance and long, it has been neglected here and the modified transfer function G0u can be derived as(30)G0u=4ρ0c0((e-s(3lu-2lv)/c0-e-s(2L-lu)/c0)+[(1/2)+H](e-s(lu/c0)-e-s(2L-2lv+lu)/c0)).
From this equation, it becomes clear that half rear-sound amplitude will eliminate the shortest time delay given by e-s(lu/c0) and leads to e-s(3lu-2lv)/c0 as the shortest time delay. Compared to the previous result for Case A (29), an additional time delay of e-s(2lu-2lv)/c0 could be achieved. This corresponds to twice the length of the subduct and can be seen in Figure 5. The time delay can also be verified by multiplying the transfer function Gyu of Case A by an additional, fictitious time delay. In Figure 7, it can be seen that this modified Case A and Case C with a rear-sound strength of −0.5 produce indeed a very similar phase lag.Figure 7
Comparison of phase characteristics for Case A, Case A with additional time delay, and Case C (H=-0.5).The numerical implementation reveals the additional time delay in feedback-path system by displaying the propagation of the incident waves. In Figures8 and 9, the wave propagation for Case A and Case C is shown. The y-axis depicts the acoustic pressure distribution for different points in time, the x-axis represents the position in the main duct, where the closed end is located at x=0. One complete period with a frequency of 1000 Hz is excited by the control loudspeaker at lu=2.32, and the propagation through the duct is observed. Note that the first upwards travelling wave is extinguished by the rear sound in Case C and previously derived time delay is obtained.Figure 8
Numeric time response, Case A.Figure 9
Numeric time response, Case C.


## 5. Conclusion

In this paper, as an additional approach in active noise control next to the Swinbanks' source, the proposed rear-sound-aided source has been analytically examined. The advantage of additional phase lag and time delay in the feedback-path transfer function of the proposed control source has been analyzed. To accomplish this, the damped wave equation has been used and based on this: a first principle model has been derived by the application of transmission line theory. The model has been confirmed by verifying that(i)
the previous model is included as the special case of an undamped scenario,(ii)
the frequency response is consistent with the one obtained by a numerical simulation based on a finite difference implementation of the damped wave equation.It has been shown, that the different consideration of damping induced some minor changes in low-frequency characteristics. Characteristics in middle and high-frequency range were consistent to the previous model, which implies validity of the previous results in [9]: The proposed control source has shown additional phase lag by assuming a weaker rear sound. The advantage has also been examined from a theoretical point of view in order to determine under which conditions an additional time delay is achieved. It has been demonstrated that the shortest time delay in the feedback-path transfer function is extended by 2LS/c0, when H=-0.5, where LS resembles the length of the subduct and c0 the speed of sound. This time delay has been also confirmed by comparison with the phase characteristic for the bidirectional source with additional time delay. In addition, the time delay has been directly confirmed in the time response simulation of the pressure variation by an incident wave.The presented results have analytically proved that the proposed control source induces additional time delay in the feedback-path transfer function as compared to the conventional bidirectional source by a suitable choice ofH. This ends up showing that the original setting of H=-1 is not suitable for the proposed source to maximize the additional time delay, which might be inconsistent with the result in [9], since the similar phase characteristic of Swinbanks’ source has been achieved by H=-1 in the experiment. To tackle this problem, the time response simulation in this paper will be used for comparison in detail with a measured time response by experiments. Further effect of the rear-sound amplitude on control performance will be examined.



---

*Source: 697108-2011-08-29.xml*