Sarabjeet Singh (1), Colin H. Hansen (2) and Carl Q. Howard (2)

Size: px
Start display at page:

Download "Sarabjeet Singh (1), Colin H. Hansen (2) and Carl Q. Howard (2)"

Transcription

1 Acoustics 8 Geeong, Victoria, Austraia 4 to 6 November 8 Acoustics and ustainabiity: How shoud acoustics adapt to meet future demands? A detaied tutoria for evauating in-duct net acoustic power transmission in a circuar duct with an attached cyindrica Hemhotz resonator using method arabjeet ingh (1), Coin H. Hansen () and Car Q. Howard () (1) Vipac Engineers and cientists Ltd., PO Box 419, Kent Town, A 571, Austraia () choo of Mechanica Engineering, The University of Adeaide, A 5, Austraia ABTRACT The method has been extensivey expored by many researchers for anaysing acoustic duct systems. However, there does not exist a comprehensive source that eaborates the use of method for evauating the performance of acoustic fiters. The current approach describes a detaied step-by-step method for evauating induct net acoustic power transmission for a harmonic pane wave traveing in a circuar duct with an attached cyindrica Hemhotz resonator using the method. The net acoustic power transmission is evauated using two different methods: (1) the product of acoustic pressure and acoustic voume veocity at the duct exit, and () estimates of the in-duct sound fied using the two-microphone technique. INTRODUCTION The method (aso known as transmission matrix or four-poe parameter presentation) has been widey accepted as a too for anaysing compex systems due to its computationa efficiency and fexibiity. Many previous researchers (Munja 1987, Pierce 1989, nowdon 1971, and To and Doige 1978) have used the method for anaysing acoustic duct systems. Athough many research papers, in addition to standard acoustic text books (Beranek and Ver 199, Kinser et a. 199, and Pierce 1989), provide a brief overview of the method and its associated uses, they do not incorporate detaied anayses of a particuar system. This paper describes in detai a step-by-step method for approximating the downstream net acoustic power transmission for a harmonic pane wave traveing inside a circuar duct to which a cyindrica Hemhotz resonator (HR) is attached. It begins with a description of the of a uniform circuar duct which incudes a the steps and pertinent cacuations in order to estimate the in-duct net acoustic power transmission. Then, the buiding of the of a couped duct-hr system is described and the net acoustic power transmission inside the duct downstream of the HR is approximated. The net acoustic power transmission inside the duct was estimated by using the estimates of acoustic pressure and acoustic voume veocity at the terminating end of the duct. However, in addition to the above stated method, net acoustic power transmission inside the duct was aso estimated using the two-microphone method deveoped by Chung and Baser (Chung and Baser 198) and extended by Åbom (Åbom 1989). The purpose of estimating the acoustic power transmission using the atter method is to incorporate additiona detaied cacuations using the method. Furthermore, the atter method provides a means to verify the estimations of the former method. UNIFORM CIRCULAR DUCT Figure 1 shows a schematic of a circuar duct of (physica) ength, radius a, and cross-sectiona area. v p s piston-driven (source) end open end radiating into free space x x v p Figure 1. A schematic of a uniform circuar duct modeed as being driven by a constant ampitude piston source mounted at one (eft) end and open at the other (right). Acoustics 8 1

2 4-6 November 8, Geeong, Austraia Proceedings of ACOUTIC 8 Aso shown at the end of the duct is the end-correction factor,, which wi be discussed ater. The eft end of the duct, referred to as a source end, was modeed as being driven by a constant ampitude piston with a unit acoustic voume veocity, and the right end of the duct, referred to as a duct exit, was modeed as an open end radiating into free space. Pane Wave Assumption As the method works on a principe of pane wave propagation inside the duct, it is important to estimate the cut-on frequency, which is defined as the frequency beow which ony pane waves propagate inside the duct. The cut-on frequency is given by: 1.841c f c (1) πa where, c is the speed of sound in air (343 m/s at C), and a is the radius of the duct. For the current paper, the dimensions of the circuar duct were taken as: diameter.1555 m and ength 3 m. The cut-on frequency for the above dimensioned duct is 19 Hz. Therefore, the pane wave mode is the ony mode which propagates beow 19 Hz. Radiation Impedance For the purpose of modeing the end of the duct as open and radiating into free space, the cacuation of radiation impedance is required. It is approximated by assuming a circuar piston ocated in an infinite baffe, and is the ratio of the force exerted by the piston on the acoustic fied to the veocity of the piston. Radiation impedance is given by (Beranek and Ver 199, Bies and Hansen 3, Pierce 1989, Kinser et a. 198, Munja 1987): r R + jx () The rea part of radiation impedance, R, is termed as radiation reactance and represents the energy radiated away from the open end in the form of sound waves. The imaginary part, X, is termed as radiation admittance and represents the mass oading of the fuid (air) just outside the open end. It is we known that the theoretica expression for the impedance of an unfanged open duct, with pane waves propagating inside it, is given by (Beranek and Ver 199, Bies and Hansen 3, Pierce 1989, Kinser et a. 198, Munja 1987): ( ka) + j(. 6) ka 4 where, is the radiation impedance, ρ is the density of the fuid medium, c is the speed of sound, k is the wave number, and a is the radius of the duct. End-Correction The term end-correction refers to an additiona ength to the actua ength of a tube or an orifice that accounts for the entrained mass of the fuid that vibrates outside the tube s or orifice s opening. The expression for the end-correction of an unfanged open end of a duct radiating into free space can be found in standard acoustica text books (Beranek and Ver 199, Bies and Hansen 3, Kinser et a. 199, and Munja 1987). It is given by: (3) where, a is the radius of the duct. Transfer Matrix of a Circuar Duct The of a circuar duct of uniform crosssectiona area,, and ength,, is given by (Munja 1987): p cos r j qr sin Yr ( k ˆ ) jy ( ˆ ) r sin k pr 1 ( k ˆ ) cos( k ˆ ) q r 1 The subscript r impies that the duct is defined/denoted as an c eement r. The quantity Y r is the characteristic impedance, and k ˆ k( 1+ iη) is the compex wave number. Com- pex wave numbers are introduced in order to convenienty incorporate damping in the system, even though this formuation impies hysteretic damping when in fact the system is better described as viscousy damped. The quantity η is the oss factor, which is equivaent to twice the critica damping ratio for a viscous system. The quantities, p r, p r-1 and q r, q r-1 represent the acoustic pressures and acoustic mass veocities at the extreme ends of the duct, respectivey (input and output sides). These quantities can be r.m.s vaues, ampitudes or instantaneous, provided they are a the same type. Rewriting equation (5) by reating acoustic voume veocity and acoustic pressure instead of acoustic mass veocity and acoustic pressure, respectivey gives: v cos r pr j sin ( ˆ k) j sin( k ˆ ) ( k ˆ ) cos( k ˆ ) vr 1 pr 1 where, v r and v r-1 represent acoustic voume veocities at the input and output sides of the duct (eement r), respectivey. The above shown transition in equation (6) affected the definition of the eements of the compared to those defined by Munja (1987) (in equation (5)). A pivota reason for impementing this transition is to enabe a quick comparison between theoretica predictions and experimenta resuts. Athough no experimenta resuts are discussed in this paper, in most of the reevant experimenta investigations reported in the iterature, a oudspeaker backed by a sma air-tight cavity, as an acoustic source, is attached to one end of the duct. And, generay, sound pressure eve measurements inside the duct are normaised by the acoustic voume veocity of the oudspeaker. Therefore, from a practica point of view, a reation between acoustic voume veocity and acoustic pressure is much preferred over a reation between acoustic mass veocity and acoustic pressure. Athough acoustic mass veocity is directy proportiona to acoustic partice veocity which, in turn, is directy proportiona to acoustic voume veocity, the equations shown in this paper eiminate the need to derive a reationship between acoustic voume veocity and pressure. As indicated in figure 1, v and p are acoustic voume veocity and pressure, respectivey, at the input side of the duct. imiary, v and p have the same definition at the output side. Reating these state variabes gives the resutant transfer matrix of the duct, which is described beow: (5) (6). 6a (4) Acoustics 8

3 Proceedings of ACOUTIC 8 v cos j sin ( k ˆ eff ) j sin( k ˆ eff ) ( k ˆ ) cos( k ˆ ) eff eff v p A very important point which needs to be considered is the end-correction factor of the unfanged open end of the duct. Whie estimating the eements of the, the endcorrection factor,, must be added to the physica ength of the duct. Therefore, in equation (7), eff respresents the effective ength of the duct, which corresponds to the sum of physica duct ength,, and the end-correction factor,. In order to make the evauation of equation (7) easier for reevant future cacuations, the eements of the above described were assigned aphabetica characters as shown beow: A C ( ) B v D p where, A D cos kˆ ( ), B ˆ j sin k( + ) + ( ) c ρ. and C ˆ j sin k( + ) (7) (8) ( ), ρ c As described earier, the in-duct net acoustic power transmission is cacuated using two different methods, (1) the product of acoustic pressure and veocity at the duct exit, and () the two-micrphone method. These methods are described in the foowing text. METHOD 1: PRODUCT OF ACOUTIC PREURE AND VOLUME VELOCITY AT THE DUCT EXIT The equations shown beow describe the procedure to estimate the acoustic pressure and voume veocity at the duct exit. Estimation of these two variabes wi in turn faciitate the estimation of in-duct net acoustic power transmission. From equation (8), we get v A v + B p (9) p By using, equation (9) can be written as: v p v A + B p (1) oving equation (1) for p with respect to v gives p v A + B (11) Equation (11) gives the estimate of the ratio of acoustic pressure p at the duct exit to the input voume veocity v. is the radiation impedance of an unfanged open end of the duct and can be cacuated using equation (3). 4-6 November 8, Geeong, Austraia B v A B + v A (1) Equation (1) gives the measure of the ratio of acoustic voume veocity v at the duct exit to the input voume veocity v. If v is set to unity (1 m 3 /s), then p and v woud be p (13) A + B B A + B v A (14) Equations (13) and (14) represent the estimates of acoustic pressure and acoustic voume veocity, respectivey, at the duct exit; the duct being driven by a constant ampitude piston at the source end with a unit voume veocity. In-Duct Net Acoustic Power Transmission The net acoustic power fux associated with a harmonic pane wave, which takes into account the pertinent incident and refected panes waves, traveing inside a duct is given by (Fahy, 1995): 1 W * net Re p v (15) where, p and v represent the compex ampitudes of pressure and voume veocity, respectivey, and the * indicates the compex conjugate. This expression can be written as: 1 W net PV cos( φ) (16) where, P and V are the scaar ampitudes that are the modui of the compex ampitudes of the pressure and voume veocity, respectivey shown in equation (15), and φ is the phase ange between them. Equation (16) represents the net acoustic power transmission inside the duct as the product of acoustic pressure and voume veocity ampitudes. ubstituting the vaues of equations (13) and (14) in equation (16) provides the estimate of net acoustic power transmission in the duct. Assuming that the anaysis in the preceding sections used ampitude quantities, equation (16) cacuates the estimate of the net acoustic power transmission in the duct when the fuid present inside the duct is assumed to be driven by a constant ampitude harmonic piston source with a unit voume veocity. Figure shows a pot of the net acoustic power transmission in the uniform circuar duct of diameter equas.1555 m and ength equas 3 m. imiary, soving equation (9) for v with respect to v, which incorporates the use of equation (11), gives: Acoustics 8 3

4 4-6 November 8, Geeong, Austraia Proceedings of ACOUTIC 8 Acoustic Impedance at Arbitrary Location Inside a Duct From standard acoustic text books (Beranek and Ver 199, Bies and Hansen 3, Kinser et a. 199, and Pierce 1989), the specific acoustic impedance at any ocation inside a duct can be expressed as x jkx + jkx c a+ exp + a exp jkx + jkx a+ exp a exp ρ (17) where, a + and a - are the moda ampitudes of incident and refected acoustic waves, respectivey, the - and + signs represent the propagation of acoustic wave in +ve and ve x directions, respectivey, and x is the ocation aong the duct. Figure. Net acoustic power transmission as a function of harmonic frequency in a uniform circuar duct driven by a constant ampitude harmonic piston source (of unit acoustic voume veocity) mounted at one (eft) end and open at the other (right). Athough the above discussed method showed the equations and procedure required to estimate the in-duct net acoustic power transmission, the use of the method is better iustrated in the second anaysis method to foow. METHOD : TWO-MICROPHONE TECHNIQUE The two-microphone technique was deveoped by Chung and Baser (Chung and Baser 198) and requires an estimate of acoustic pressure at two different ocations inside a duct. The foowing paragraphs detai the procedure for estimating the acoustic pressure at different ocations inside a duct. Acoustic Pressure and Veocity at Arbitrary Locations Inside a Duct Figure 3 shows a schematic of a uniform duct of ength depicting the end-correction,, at the open end of the duct. v piston-driven (source) end v x open end radiating into free space v Using equation (17), the vaues of acoustic impedance at ocations x and x are given by: a+ + a (18) a+ a jk + jk c a+ exp + a exp jk + jk a+ exp a exp ρ (19) Eiminating constants a + and a - faciitates a reationship between and, which is given by: + j tan( k) () 1+ j tan( k) If the impedance at a duct termination is known, the vaue of the acoustic impedance at any point inside the duct can be cacuated using equation (). In order to distinguish between the source and radiation impedance, equation () was rewritten with a sight modification in terms of nomencature. Denoting the acoustic impedance at x by s (source impedance) and at x by r (radiation impedance), gives: r + j tan( k) s (1) 1+ j r tan( k) p p x x p Finay, the vaue of acoustic impedance at ocation x ooking into the duct from the source end (x ) can be given by (utiising equation (1)): x -x x x Figure 3. A schematic of a uniform duct depicting ocation x aong the duct where acoustic pressure and acoustic voume veocity need to be estimated. uppose that acoustic pressure and acoustic voume veocity at ocation x need to be estimated, and are denoted by p x and v x, respectivey. As evident from equations (11) and (1), the estimation of acoustic pressure and voume veocity at a particuar ocation requires the corresponding vaue of acoustic impedance at that particuar ocation. Therefore, estimation of acoustic pressure and voume veocity at ocation x requires estimating the acoustic impedance at ocation x, which is detaied beow. s j tan( kx) x () j s tan( kx) The aternate expression for cacuating the vaue of x, by using the vaue of r instead of s, is given by: r + j tan x 1+ j s tan ( k( x) ) ( k( x) ) (3) Equation (3) shows the acoustic impedance at ocation x ooking into the duct from the open end (x ). 4 Acoustics 8

5 Proceedings of ACOUTIC 8 End-Corrections Before using equations () and (3) for cacuating the vaue of x, the end-correction of an unfanged open end of a duct must be accounted. After incorporating the vaue of the end-correction in equations () and (3), x is expressed as: s j tan x j s tan r + j tan x 1+ j r tan ( k( x )) ( k( x )) ( k( x + )) ( k( x + )) (4) (5) Once the vaue of acoustic impedance at a particuar ocation in the duct is known, the acoustic pressure and veocity can be cacuated using equations (11) and (1), respectivey. The acoustic pressure and voume veocity at ocation x with respect to unit voume veocity, v, are given by: p x x (6) Ax + Bx x Bx x Ax + Bx v x x Ax where, ( kx ˆ cos ) and j sin( kx ˆ ) A x B x. (7) Knowedge of the acoustic pressure at two appropriatey spaced different ocations inside a duct estimates the in-duct net acoustic power transmission by using the decomposition of sound fied which is discussed beow. Decomposition of the ound Fied Decomposition is a technique for determining the ampitudes of the acoustic waves propagating each way inside a duct. It faciitates the separation of modes into incident and refected parts. Once the acoustic power associated with the incident and refected waves is determined, the net acoustic power transmission is simpy their difference. Åbom (Åbom 1989) extended the two-microphone technique deveoped by Chung and Baser (Chung and Baser 198) and documented the scheme for in-duct moda decomposition. However, in this paper we are not concerned with higher order modes and we are ony decomposing the sound fied into right-traveing and eft-traveing pane waves. Generay, for a uniform straight duct with rigid was and an arbitrary cross-sectiona shape carrying an axia mean fow, the sound fied in the duct can be expressed as: n ik x ik x [ aˆ n + + a n + exp + ˆ exp ] Ψ n pˆ (8) where, pˆ is the acoustic pressure, aˆ +, aˆ are the moda ampitudes of acoustic pressure associated with the incident (+ve x direction) and refected waves (-ve x direction), respectivey, k ˆ + n, k ˆ n are the axia wave numbers in the positive and 4-6 November 8, Geeong, Austraia negative directions, respectivey, n is the moda number, and Ψ is the eigenfunction for mode n. n For the case of a uniform circuar duct with rigid-was and carrying no axia fow, with ony pane waves propagating inside it, equation (8) can be written as: ikx ikx pˆ + aˆ + exp + aˆ exp (9) The terms aˆ + and aˆ can be cacuated by either estimating or experimentay measuring the acoustic pressure at two ocations in a duct and using the two-microphone technique. Expressing the resuts of such a measurement in a matrix formuation yieds (Åbom 1989): where, pˆ is a [ 1] pˆ Maˆ (3) coumn vector containing the estimates or measures of the acoustic pressures at two different oca- moda matrix containing the propagation tions, M is a [ ] terms, and â is a [ 1] coumn vector containing the unknown moda ampitudes. Writing vector pˆ in terms of transfer function between the two pressure measuring or estimating ocations yieds: pˆ pˆ r H (31) where, pˆ r is the acoustic pressure at one of the two pressure measuring or estimating ocations, and H is the [ 1] coumn vector containing the transfer function between the two ocations. Using equation (31), equation (3) can be written as: oving equation (3) for â gives: aˆ p ˆ rh Maˆ (3) 1 p ˆ M H (33) r aˆ pˆ r + ˆ iks iks (34) a exp exp H1 For this anaysis, pˆ r was taken as the acoustic pressure estimate at the ocation coser to the source end, H 1 represents the pressure transfer function between two pressure estimates, and s is the axia distance between the two pressure estimating ocations. oving equation (34), for the estimates of â+ and â, yieds: + iks ( exp ) pˆ r H aˆ iks iks exp exp iks ( exp ) pˆ r H aˆ 1 + iks iks exp exp (35) (36) The acoustic power associated with the incident and refected waves is expressed as: aˆ + W+ (37) Acoustics 8 5

6 4-6 November 8, Geeong, Austraia Proceedings of ACOUTIC 8 aˆ W (38) The in-duct net acoustic power transmission in x direction is estimated as the difference between W + and W -, and is expressed as: W net W+ W aˆ + aˆ In-Duct Net Acoustic Power Transmission (39) Assuming that the anaysis in the preceding sections used ampitude quantities, equation (39) cacuates the estimate of the net acoustic power transmission in the duct when the fuid present in the duct is assumed to be driven by a constant ampitude harmonic piston source with a unit voume veocity. The net acoustic power estimates from the two mehods described above exacty match with one another. COUPLED DUCT HR YTEM Mounting a resonator on a duct reduces the noise transmission at a particuar frequency which is a characteristic of the resonator s geometry. Comparison of the acoustic power transmission inside the duct without and with the resonator can be used to evauate its performance. In order to deveop the compete of the duct- HR system, it was discretised into three eements: (1) eement - section of the duct upstream of the HR, () eement 1 - section of the duct downstream of the HR, and (3) eement HR - the Hemhotz resonator. Figure 4 shows a schematic of a HR mounted on to a duct of ength at ocation x. The figure aso iustrates the three eements described above, aong with the neck-cavity and neck-duct interfaces. eement HR s p s, v s v HR p s1, v s1 v, p v 1, p 1 neck-cavity interface neck-duct interface eement eement 1 v p r A key issue reated to the theoretica anaysis of a duct-hr system is the incorporation of two end-correction factors for the neck of the HR; (1) the neck-cavity interface endcorrection factor, and () the neck-duct interface end correction factor, both of which must be added to the physica neck ength to get the effective neck ength. The expression for the neck-cavity interface end-correction factor cosey approximates the geometry of a cyindrica piston radiating in a tube. It is given by (eamat and Ji ): r δ e.8r (4) R where, r denotes the radius of the neck of a HR, and R denotes the radius of the cavity of a HR. Another equation which has been widey used in the iterature for the neck-cavity interface end-correction factor was presented by Ingard (Ingard 1953). His equation sighty differs from equation (4) and was for an orifice in an anechoicay terminated duct. However, equation (4) is considered to be better than the equation presented by Ingard and has been used in this current paper. On the other hand, the neck-duct interface end-correction factor has not been modeed anayticay due to the difficuty in interpreting the non-panar sound fied in the region of the duct adjacent to the neck opening. However, Ji (5) derived two-curve fitting expressions for the circuar neck-duct interface end-correction factor based on the ratio of neck and duct diameter using Boundary Eement Anaysis. The corresponding expressions are given by: r r δ i r a a ; r.4 a r r r ; >.4 a a (41) where, r denotes the radius of the neck of a HR, and a denotes the radius of a duct to which a HR is attached. Transfer Matrix of a HR Figure 4 shows a uniform duct with a HR attached as a side branch to one its was at ocation x, aong with the reevant notations required to deveop the compete of the duct-hr system. The compex ampitude of the voume veocity entering the duct, before the HR's ocation, is denoted by v. imiary, it is v 1 in the section after the resonator's ocation and v HR fowing into the HR. From the continuity condition, the voume veocity and pressure are ocay conserved at junction x. Therefore, ( x) v ( x) v ( x) v HR + 1 (4) ( x) p ( x) p ( x) p HR 1 (43) x x x Figure 4. A schematic of a couped duct-hr system depicting its division into three eements aong with neck-duct and neck-cavity interfaces, and reevant notations. End-Corrections of a HR s Neck Putting equations (4) and (43) in a matrix form yieds: vhr + v p 1 vhr v p1 p 1 (44) (45) 6 Acoustics 8

7 Proceedings of ACOUTIC 8 From equation (43), p ( x) p ( x) HR 1 HR and aso, HR vhr p. Here, HR is the acoustic impedance just outside the opening of the resonator, which can be referred to as the input point impedance. Therefore, equation (45) can be re-written as: 1 1 v HR p 1 (46) 1 Finay, the matrix 1 represents the HR 1 of the Hemhotz resonator. Input Point Impedance of a HR The impedance at the opening of the HR, i.e. HR, can be estimated by appying the equation to the neck and cavity of the HR. The neck section of the resonator is denoted by s 1 having cross-sectiona area A n and ength 1, and the cavity section is denoted by s having cross-sectiona area A c and ength. The ength 1 represents the effective ength of the neck, which incudes two end-corrections at each side of the neck, which were described in the previous section beow figure 4 (equations (4) and (41)). Reating the acoustic voume veocity and pressure at the opening of the neck v s1, p s1 (see figure 4) and at the end of the cavity v s, p s (see figure 4), the equation for the Hemhotz resonator is expressed as: or vs1 vs ps1 ps of the neck of the cavity vs1 As1 ps1 Cs1 Bs1 As Ds1 Cs Bsvs Ds ps where, ( ˆ A D cos k ), j sin( kˆ ) C s 1 B s s 1 s1 1 (47) (48) B s 1 1, and j sin( kˆ 1 ). imiary, A ( ˆ s Ds cos k ) sin( ˆ j k ), and j sin( kˆ ). C s, Let us denote the resutant eements of equation (48) as: vs1 As ps1 Cs Bsvs Ds ps (49) As the end of the resonator (cavity) is bocked (rigidy terminated), there wi be no fow at the end of the resonator, hence, v. ubstituting v in equation (49), the s s acoustic pressure and voume veocity at the opening of the neck becomes v B p (5) s1 s s 4-6 November 8, Geeong, Austraia p D p (51) s1 s s Dividing equation (51) by (5), gives the estimate of the input point impedance of the HR, as: ps1 Ds HR vs1 Bs Compete Transfer Matrix of a Duct-HR ystem (5) Referring to figure 4, for the convenience of presenting the compete of the duct-hr system, et us denote: (1) eement, the first section of the duct upstream of the HR, by subscript, and () eement 1, the second section of the duct downstream of the HR, by subscript 1. The compete method of the duct-hr system shown in figure 4 can be expressed as: v p of eement of eement HR of eement 1 A C B 1 D 1 A1 HR C 1 1 B1 v D1 p (53) (54) Let the resutant of matrix mutipication of equation (54) be denoted by: A C B v D p (55) In-Duct Net Acoustic Power Transmission Downstream of the HR The dimensions of the duct were kept same as were considered in the anaysis of the uniform circuar duct (diameter of.1555 m and ength of 3 m). The dimensions of the cyindrica Hemhotz resonator were: cavity diameter.131 m, cavity ength.7 m, neck diameter.55 m, and physica neck ength.93 m. The HR was mounted at a distance of.5 m from the source end of the duct. Method 1: Product of the acoustic pressure and voume veocity at the duct exit. The transfer matrices of three eements of the duct-hr system shown in figure 4, eement, eement 1, and eement HR, were written using equation (54) as per the dimensions of the duct-hr system described in the immediate preceding paragraph. The acoustic pressure and acoustic voume veocity at the duct exit (open end radiating into free space), x, were estimated in the simiar way as described in equations (11) and (1), and are given by: p v A + B B v A + B v A (56) (57) Acoustics 8 7

8 4-6 November 8, Geeong, Austraia Proceedings of ACOUTIC 8 where, p and v represent the acoustic pressure and voume veocity at the duct exit, v denotes the input voume veocity, is the radiation impedance at the unfanged open end of the duct (aso denoted by r ), and A, B, C and D are the eements of the resutant compete (equation (55)). Mutipying the estimates of acoustic pressure (equation (56)) and acoustic voume veocity (equation (57)) at the duct exit as per equation (16) resuts in the net acoustic power transmission inside the duct downstream of the HR. Method : Two-Microphone Decompostion of ound Fied The net acoustic power transmission inside the duct downstream of the HR was aso estimated by using the twomicrophone in-duct decomposition method and exacty the same resuts as for Method 1 were obtained. Figure 5 shows a pot of the net acoustic power transmission in the duct downstream of the HR. For comparison purposes, this figure aso incudes the pot of the net acoustic power transmission in the duct without the resonator. factors of the neck of the HR in addition to the actua dimensions of the duct-hr system. The neck-cavity interface endcorrection factor is we documented. However, there exists some uncertainty concerning the end-correction factor at the neck-duct interface, due to the difficuty in anayticay modeing the compex sound fied in the regions of the duct adjacent to the neck opening. However, empirica equations have been formuated by Ji (5). Even though Ji found that his empirica expressions agreed with his experimenta resuts, in another study (ingh et a. 6 and ingh 6), it was observed that the experimenta performance of a duct mounted HR did not match the corresponding theoretica estimations based on the empirica end-correction factor reported by Ji. ingh (6) deveoped and soved a three-dimensiona finite eement mode of the above described duct-hr system using the software package ANY. Unike the theoretica anaysis of a duct-hr system in which the end-correction factors have to be cacuated, ANY automaticay determines and incorporates the end-correction factors during its soution phase, and therefore, ANY resuts are considered more reiabe than the theoretica cacuations which differ sighty from the ANY resuts. In order to further identify the vaidity of the ANY resuts, experiments were conducted, and the experimentay measured performance of the HR mounted on to a duct was found to agree with the ANY resuts. The difficuty in accuratey estimating the neck duct interface end correction factor compromises the abiity of the transfer matrix method to provide an accurate performance estimation of the duct mounted HR. Figure 5. Net acoustic power transmission in a uniform circuar duct with a HR mounted on to it at a distance of.5 m from the soure end; the duct being driven by a constant ampitude harmonic piston source mounted at one (eft) end and open at the other (right). DICUION The effect of mounting the HR on to the duct, as evident from Figure 5, minimised the net acoustic power transmission by db at 4 Hz. Generay, a frequency at which the maximum reduction of in-duct net acoustic power transmission downstream of a HR occurs is considered to be the resonance frequency of the HR, which is 4 Hz in the current anaysis. However, when a HR is mounted on to a duct, a couped system is created whose resonance frequency is different to that of the HR as a stand-aone device. Therefore, the frequency at which the maximum reduction of in-duct net acoustic power transmission downstream of a HR occurs does not correspond to the resonance frequency of the HR (ingh et a. 6 and ingh 6). The acoustic performance obtained using the method for the circuar duct-cyindrica HR system shoud ony be considered as approximate. This is because of the duct-hr system reated imitations of the method, which are discussed beow. DUCT-HR YTEM RELATED LIMITATION OF THE TRANFER MATRIX METHOD The key issue in the deveopment of the duct-hr system is the incorporation of the two end-correction Aso, as of the magnitude of the damping that occurs in a practica duct-hr system cannot be accuratey modeed, the method with damping excuded can resut in unreaisticay high estimates of acoustic power reduction. However, in this paper, damping has been accounted for in the method by using the vaue of oss factor, η, equa to.4. CONCLUION A detaied method for evauating the acoustic performance of a cyindrica Hemhotz resonator mounted on to a uniform circuar duct using the method has been presented. The in-duct net acoustic power transmission was estimated using two different methods: (1) product of the acoustic pressure and acoustic voume veocity at the openended duct exit, and () in-duct decomposition of sound fied. The net acoustic power resuts from both methods exacty matched each other. The method for the described duct-hr system appication ony provides a rough estimate of the acoustic performance of the duct mounted HR in terms of the frequency corresponding to the maximum net acoustic power reduction and the amount of net acoustic power reduction. The imitations corresponding to the duct-hr system are aso described. REFERENCE Åbom, M. 1989, Moda decomposition in ducts based on transfer function measurements between microphone pairs, Journa of ound and Vibration, 135(1), Beranek, L.L. and Ver, I.L. 199, Noise and Vibration Contro Engineering, John Wiey and ons. Bies, D.A. and Hansen, C.H. 3, Engineering Noise Contro, pon Press. 8 Acoustics 8

9 Proceedings of ACOUTIC November 8, Geeong, Austraia Chung, J.Y. and Baser, D.A. 198, Transfer function method for measuring in-duct acoustic properties. I theory, II experiments, Journa of Acoustica ociety of America, 68(3), Fahy, F.J. 1995, ound Intensity, E&F pon, econd Edition. Ingard, U. 1953, On the theory and design of acoustic resonators, Journa of Acoustica ociety of America, 5(6), Ji,.L. 5, Acoustic ength correction of a cosed cyindrica side-branched tube, Journa of ound and Vibration, 83, Kinser, L.E., Frey, A.R., Coppens, A.B. and anders, J.V. 198, Fundamentas of Acoustics, John Wiey and ons. Munja, M.L. 1987, Acoustics of Ducts and Muffers, John Wiey and ons. Pierce, A.D. 1989, Acoustics, The Acoustica ociety of America. eamat, A. and Ji,.L., Circuar asymmetric Hemhotz resonators, Journa of the Acoustica ociety of America, 17 (5), ingh,., Hansen, C.H. and Howard, C.Q. 6, The eusive cost function for tuning adaptive Hemhotz resonators, Proceedings of the Austraian Acoustica ociety Conference 6, Christchurch, New eaand ingh,. 6, Tona noise attenuation in ducts by optimising adaptive Hemhotz resonators, Masters Thesis, The University of Adeaide, Austraia. nowdon, J.C. 1971, Mechanica four poe parameters and thier appications, Journa of ound and Vibration, 15(3), To, C.W.. and Doige, A.G. 1978, A transient testing technique for the determination of matrix parameters of acoustic systems. I theory and principes, II experimenta procedures and resuts, Journa of ound and Vibration, 6(), 7-3. Acoustics 8 9

Proceedings of Meetings on Acoustics

Proceedings of Meetings on Acoustics Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.

More information

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage Magnetic induction TEP Reated Topics Maxwe s equations, eectrica eddy fied, magnetic fied of cois, coi, magnetic fux, induced votage Principe A magnetic fied of variabe frequency and varying strength is

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

Chapter 4. Moving Observer Method. 4.1 Overview. 4.2 Theory

Chapter 4. Moving Observer Method. 4.1 Overview. 4.2 Theory Chapter 4 Moving Observer Method 4.1 Overview For a compete description of traffic stream modeing, one woud reuire fow, speed, and density. Obtaining these parameters simutaneousy is a difficut task if

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

(Refer Slide Time: 2:34) L C V

(Refer Slide Time: 2:34) L C V Microwave Integrated Circuits Professor Jayanta Mukherjee Department of Eectrica Engineering Indian Intitute of Technoogy Bombay Modue 1 Lecture No 2 Refection Coefficient, SWR, Smith Chart. Heo wecome

More information

P. T. Thawani Climate Control Operations, Ford Motor Company, S. Commerce Drive, Dearborn, Michigan 48120

P. T. Thawani Climate Control Operations, Ford Motor Company, S. Commerce Drive, Dearborn, Michigan 48120 Prediction of the vibro-acoustic transmission oss of panar hose-pipe systems M. L. Munja Department of Mechanica Engineering, Indian Institute of Science, Bangaore 560 012, India P. T. Thawani Cimate Contro

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE 3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses

More information

Forces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment

Forces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment Forces of Friction When an object is in motion on a surface or through a viscous medium, there wi be a resistance to the motion This is due to the interactions between the object and its environment This

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

Traffic data collection

Traffic data collection Chapter 32 Traffic data coection 32.1 Overview Unike many other discipines of the engineering, the situations that are interesting to a traffic engineer cannot be reproduced in a aboratory. Even if road

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, we as various

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Chapter 2 Theoretical Preliminaries of Acoustics

Chapter 2 Theoretical Preliminaries of Acoustics Chapter 2 Theoretica Preiminaries of Acoustics In this chapter, we review some of the fundamentas of acoustics and introduce the spherica harmonic expansion of a sound fied, which is the basis for the

More information

High Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings

High Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings Purdue University Purdue e-pubs Internationa Compressor Engineering Conference Schoo of Mechanica Engineering 2010 High Efficiency Deveopment of a Reciprocating Compressor by Carification of Loss Generation

More information

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE

CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE Progress In Eectromagnetics Research, PIER 30, 73 84, 001 CONCHOID OF NICOMEDES AND LIMACON OF PASCAL AS ELECTRODE OF STATIC FIELD AND AS WAVEGUIDE OF HIGH FREQUENCY WAVE W. Lin and Z. Yu University of

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations .615, MHD Theory of Fusion ystems Prof. Freidberg Lecture : The Moment Equations Botzmann-Maxwe Equations 1. Reca that the genera couped Botzmann-Maxwe equations can be written as f q + v + E + v B f =

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

SEA Subsystem Attribute Correction Based on Stiffness Multipliers

SEA Subsystem Attribute Correction Based on Stiffness Multipliers Internationa Conference on Eectrica, Mechanica and Industria Engineering (ICEMIE 016) SEA Subsystem Attribute Correction Based on Stiffness Mutipiers Jie Gao1,* and Yong Yang 1 Coege of hysics and Information

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

Elements of Kinetic Theory

Elements of Kinetic Theory Eements of Kinetic Theory Statistica mechanics Genera description computation of macroscopic quantities Equiibrium: Detaied Baance/ Equipartition Fuctuations Diffusion Mean free path Brownian motion Diffusion

More information

PHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I

PHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I 6 n terms of moment of inertia, equation (7.8) can be written as The vector form of the above equation is...(7.9 a)...(7.9 b) The anguar acceeration produced is aong the direction of appied externa torque.

More information

Effect of Oxygen Injection into Argon Induction Plasmas on Chemically Non-Equilibrium Conditions

Effect of Oxygen Injection into Argon Induction Plasmas on Chemically Non-Equilibrium Conditions Proceedings of 17th Internationa Symposium on Pasma Chemistry, Toronto, Canada, August 7-12, 25 Effect of Oxygen Injection into Argon Induction Pasmas on Chemicay Non-Equiibrium Conditions Nobuhiko Atsuchi

More information

Research on liquid sloshing performance in vane type tank under microgravity

Research on liquid sloshing performance in vane type tank under microgravity IOP Conference Series: Materias Science and Engineering PAPER OPEN ACCESS Research on iquid soshing performance in vane type tan under microgravity Reated content - Numerica simuation of fuid fow in the

More information

Numerical simulation of javelin best throwing angle based on biomechanical model

Numerical simulation of javelin best throwing angle based on biomechanical model ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Transmission Matrix Model of a Quarter-Wave-Tube with Gas Temperature Gradients

Transmission Matrix Model of a Quarter-Wave-Tube with Gas Temperature Gradients Transmission Matrix Model of a Quarter-Wave-Tube with Gas Temperature Gradients Carl Howard School of Mechanical Engineering, University of Adelaide, South Australia, Australia ABSTRACT A transmission

More information

Function Matching Design of Wide-Band Piezoelectric Ultrasonic Transducer

Function Matching Design of Wide-Band Piezoelectric Ultrasonic Transducer Function Matching Design of Wide-Band Piezoeectric Utrasonic Transducer Yingyuan Fan a, Hongqing An b Weifang Medica University, Weifang, 261053, China a yyfan@126.com, b hongqingan01@126.com Abstract

More information

Ultrasonic Measurements of Kinematic Viscosity for Analize of Engine Oil Parameters

Ultrasonic Measurements of Kinematic Viscosity for Analize of Engine Oil Parameters th European Conference on Non-Destructive Testing (ECNDT 04), October 6-0, 04, Prague, Czech Repubic More Info at Open Access Database www.ndt.net/?id=6344 Utrasonic Measurements of Kinematic Viscosity

More information

Candidate Number. General Certificate of Education Advanced Level Examination June 2010

Candidate Number. General Certificate of Education Advanced Level Examination June 2010 Centre Number Surname Candidate Number For Examiner s Use Other Names Candidate Signature Examiner s Initias Genera Certificate of Education Advanced Leve Examination June 2010 Question 1 2 Mark Physics

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation 22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................

More information

Transmission Matrix Model of a Quarter-Wave-Tube with Gas Temperature Gradients

Transmission Matrix Model of a Quarter-Wave-Tube with Gas Temperature Gradients Proceedings of Acoustics 2013 Victor Harbor Transmission Matrix Model of a Quarter-Wave-Tube with Gas Temperature Gradients Carl Howard School of Mechanical Engineering, University of Adelaide, South Australia,

More information

MA 201: Partial Differential Equations Lecture - 10

MA 201: Partial Differential Equations Lecture - 10 MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary

More information

14-6 The Equation of Continuity

14-6 The Equation of Continuity 14-6 The Equation of Continuity 14-6 The Equation of Continuity Motion of rea fuids is compicated and poory understood (e.g., turbuence) We discuss motion of an idea fuid 1. Steady fow: Laminar fow, the

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

Multiple Beam Interference

Multiple Beam Interference MutipeBeamInterference.nb James C. Wyant 1 Mutipe Beam Interference 1. Airy's Formua We wi first derive Airy's formua for the case of no absorption. ü 1.1 Basic refectance and transmittance Refected ight

More information

Distribution Systems Voltage Profile Improvement with Series FACTS Devices Using Line Flow-Based Equations

Distribution Systems Voltage Profile Improvement with Series FACTS Devices Using Line Flow-Based Equations 16th NATIONAL POWER SYSTEMS CONFERENCE, 15th-17th DECEMBER, 010 386 Distribution Systems otage Profie Improvement with Series FACTS Devices Using Line Fow-Based Equations K. enkateswararao, P. K. Agarwa

More information

Simplified analysis of EXAFS data and determination of bond lengths

Simplified analysis of EXAFS data and determination of bond lengths Indian Journa of Pure & Appied Physics Vo. 49, January 0, pp. 5-9 Simpified anaysis of EXAFS data and determination of bond engths A Mishra, N Parsai & B D Shrivastava * Schoo of Physics, Devi Ahiya University,

More information

Parallel-Axis Theorem

Parallel-Axis Theorem Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states

More information

APPENDIX C FLEXING OF LENGTH BARS

APPENDIX C FLEXING OF LENGTH BARS Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or

More information

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity

More information

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,

More information

TM Electromagnetic Scattering from 2D Multilayered Dielectric Bodies Numerical Solution

TM Electromagnetic Scattering from 2D Multilayered Dielectric Bodies Numerical Solution TM Eectromagnetic Scattering from D Mutiayered Dieectric Bodies Numerica Soution F. Seydou,, R. Duraiswami, N.A. Gumerov & T. Seppänen. Department of Eectrica and Information Engineering University of

More information

SECTION A. Question 1

SECTION A. Question 1 SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second

More information

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

International Journal of Advance Engineering and Research Development

International Journal of Advance Engineering and Research Development Scientific Journa of Impact Factor (SJIF): 4.4 Internationa Journa of Advance Engineering and Research Deveopment Voume 3, Issue 3, March -206 e-issn (O): 2348-4470 p-issn (P): 2348-6406 Study and comparison

More information

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL

More information

Self Inductance of a Solenoid with a Permanent-Magnet Core

Self Inductance of a Solenoid with a Permanent-Magnet Core 1 Probem Sef Inductance of a Soenoid with a Permanent-Magnet Core Kirk T. McDonad Joseph Henry Laboratories, Princeton University, Princeton, NJ 08544 (March 3, 2013; updated October 19, 2018) Deduce the

More information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur odue 2 naysis of Staticay ndeterminate Structures by the atri Force ethod Version 2 E T, Kharagpur esson 12 The Three-oment Equations- Version 2 E T, Kharagpur nstructiona Objectives fter reading this

More information

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c) A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,

More information

O9e Fringes of Equal Thickness

O9e Fringes of Equal Thickness Fakutät für Physik und Geowissenschaften Physikaisches Grundpraktikum O9e Fringes of Equa Thickness Tasks 1 Determine the radius of a convex ens y measuring Newton s rings using ight of a given waveength.

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

MA 201: Partial Differential Equations Lecture - 11

MA 201: Partial Differential Equations Lecture - 11 MA 201: Partia Differentia Equations Lecture - 11 Heat Equation Heat conduction in a thin rod The IBVP under consideration consists of: The governing equation: u t = αu xx, (1) where α is the therma diffusivity.

More information

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006 expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

TWO- AND THREE-DIMENSIONAL SIMULATION OF A RISING BUBBLE AND FALLING DROPLET USING LEVEL SET METHOD

TWO- AND THREE-DIMENSIONAL SIMULATION OF A RISING BUBBLE AND FALLING DROPLET USING LEVEL SET METHOD European Conference on Computationa Fuid Dynamics ECCOMAS CFD 2006 P. Wesseing, E. Oñate, J. Périaux (Eds) TU Deft, The Netherands, 2006 TWO- AND THREE-DIMENSIONAL SIMULATION OF A RISING BUBBLE AND FALLING

More information

Measurement of acceleration due to gravity (g) by a compound pendulum

Measurement of acceleration due to gravity (g) by a compound pendulum Measurement of acceeration due to gravity (g) by a compound penduum Aim: (i) To determine the acceeration due to gravity (g) by means of a compound penduum. (ii) To determine radius of gyration about an

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary

More information

Effect of effective length of the tube on transmission loss of reactive muffler

Effect of effective length of the tube on transmission loss of reactive muffler Effect of effective length of the tube on transmission loss of reactive muffler Gabriela Cristina Cândido da SILVA 1 ; Maria Alzira de Araújo NUNES 1 1 University of Brasilia, Brazil ABSTRACT Reactive

More information

Combining reaction kinetics to the multi-phase Gibbs energy calculation

Combining reaction kinetics to the multi-phase Gibbs energy calculation 7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation

More information

FRIEZE GROUPS IN R 2

FRIEZE GROUPS IN R 2 FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the

More information

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18 Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms

More information

Impedance correction for a branched duct in a thermoacoustic air-conditioner

Impedance correction for a branched duct in a thermoacoustic air-conditioner 20-22 November 2006, Christchurch, New ealand Impedance correction for a branched duct in a thermoacoustic air-conditioner Carl Howard School of Mechanical Engineering, The University of Adelaide, Adelaide,

More information

A nodal collocation approximation for the multidimensional P L equations. 3D applications.

A nodal collocation approximation for the multidimensional P L equations. 3D applications. XXI Congreso de Ecuaciones Diferenciaes y Apicaciones XI Congreso de Matemática Apicada Ciudad Rea, 1-5 septiembre 9 (pp. 1 8) A noda coocation approximation for the mutidimensiona P L equations. 3D appications.

More information

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

Physics 566: Quantum Optics Quantization of the Electromagnetic Field

Physics 566: Quantum Optics Quantization of the Electromagnetic Field Physics 566: Quantum Optics Quantization of the Eectromagnetic Fied Maxwe's Equations and Gauge invariance In ecture we earned how to quantize a one dimensiona scaar fied corresponding to vibrations on

More information

The influence of temperature of photovoltaic modules on performance of solar power plant

The influence of temperature of photovoltaic modules on performance of solar power plant IOSR Journa of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vo. 05, Issue 04 (Apri. 2015), V1 PP 09-15 www.iosrjen.org The infuence of temperature of photovotaic modues on performance

More information

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport

More information

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated

More information

1. Measurements and error calculus

1. Measurements and error calculus EV 1 Measurements and error cacuus 11 Introduction The goa of this aboratory course is to introduce the notions of carrying out an experiment, acquiring and writing up the data, and finay anayzing the

More information

Science & Technologies COMPARISON OF MASS TRANSFER COEFFICIENTS DETERMINED BY DIFFERENT METHODS IN DISTILLATION COLUMN WITH THREE TRAYS

Science & Technologies COMPARISON OF MASS TRANSFER COEFFICIENTS DETERMINED BY DIFFERENT METHODS IN DISTILLATION COLUMN WITH THREE TRAYS COMPARION OF MA TRANFER COEFFICIENT DETERMINED BY DIFFERENT METHOD IN DITIATION COUMN WITH THREE TRAY Dian Radev, Dorin Georgiev, Desisava Koeva, Mariana Karaivanova Facuty of Technica ciences, Prof. D-r.

More information

arxiv: v1 [math.ca] 6 Mar 2017

arxiv: v1 [math.ca] 6 Mar 2017 Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear

More information

Acoustic Nondestructive Testing and Measurement of Tension for Steel Reinforcing Members

Acoustic Nondestructive Testing and Measurement of Tension for Steel Reinforcing Members Acoustic Nondestructive Testing and Measurement of Tension for Stee Reinforcing Members Part 1 Theory by Michae K. McInerney PURPOSE: This Coasta and Hydrauics Engineering Technica Note (CHETN) describes

More information

Radiation Fields. Lecture 12

Radiation Fields. Lecture 12 Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence

More information

Aerodynamic Design of Turbine Blades Using an Adjoint Equation Method

Aerodynamic Design of Turbine Blades Using an Adjoint Equation Method 43rd AIAA Aerospace Sciences Meeting and Exhibit, January 10-13, 2005, Reno, Nevada Aerodynamic Design of Turbine ades Using an Adjoint Equation Method Hsiao-Yuan Wu, and Feng Liu Department of Mechanica

More information

Modal analysis of a multi-blade system undergoing rotational motion

Modal analysis of a multi-blade system undergoing rotational motion Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa

More information

Version 2.2 NE03 - Faraday's Law of Induction

Version 2.2 NE03 - Faraday's Law of Induction Definition Version. Laboratory Manua Department of Physics he University of Hong Kong Aims o demonstrate various properties of Faraday s Law such as: 1. Verify the aw.. Demonstrate the ighty damped osciation

More information

Acoustic-Structure Simulation of Exhaust Systems by Coupled FEM and Fast BEM

Acoustic-Structure Simulation of Exhaust Systems by Coupled FEM and Fast BEM Acoustic-Structure Simuation of Exhaust Systems by Couped FEM and Fast BEM Lothar Gau, Michae Junge Institute of Appied and Experimenta Mechanics, University of Stuttgart Pfaffenwadring 9, 755 Stuttgart,

More information

Electromagnetic Waves

Electromagnetic Waves Eectromagnetic Waves Dispacement Current- It is that current that comes into existence (in addition to conduction current) whenever the eectric fied and hence the eectric fux changes with time. It is equa

More information

A sta6s6cal view of entropy

A sta6s6cal view of entropy A sta6s6ca view of entropy 20-4 A Sta&s&ca View of Entropy The entropy of a system can be defined in terms of the possibe distribu&ons of its moecues. For iden&ca moecues, each possibe distribu&on of moecues

More information

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest

More information

ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS. Abstract

ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS. Abstract ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS Dr. Costas P. Antonopouos, University of Patras, Greece Assoc. Prof. Thanasis C. Triantafiou, University of Patras, Greece Abstract

More information

Faculty of Machine Building. Technical University of Cluj Napoca

Faculty of Machine Building. Technical University of Cluj Napoca Facuty of Machine Buiding Technica University of Cuj Napoca CONTRIBUTIONS TO THE CALCULATION AND ANALYSIS OF DYNAMIC ABSORBERS, WITH APPLICATIONS ON BALANCING MECHANICAL SYSTEMS PhD THESIS 11 Scientific

More information

Incorporation of surface tension to interface energy balance in crystal growth

Incorporation of surface tension to interface energy balance in crystal growth Cryst. Res. Techno. 42, No. 9, 914 919 (2007) / OI 10.1002/crat.200710927 Incorporation of surface tension to interface energy baance in crysta growth M. Yidiz and. ost* Crysta Growth aboratory, epartment

More information