ON SOUND RADIATION FROM AN OPEN-ENDED NON-UNIFORMLY LINED CYLINDRICAL NOZZLE

Size: px
Start display at page:

Download "ON SOUND RADIATION FROM AN OPEN-ENDED NON-UNIFORMLY LINED CYLINDRICAL NOZZLE"

Transcription

1 24th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES ON SOUND RADIATION FROM AN OPEN-ENDED NON-UNIFORMLY LINED CYLINDRICAL NOZZLE L. M. B. C. Campos, J. M. G. S. Oliveira Secção e Mecânica Aeroespacial, Instituto Superior Técnico, Lisboa Coex, Portugal Keywors: uct, impeance, noise attenuation, soun raiation Abstract Acoustic liners are wiely use in jet engine inlet an exhaust ucts, as a passive means of noise reuction. One possible way of liner optimization is the use of a non-uniform acoustic impeance istribution. It has been shown that such a liner can lea to an attenuation in the acoustic amplitue of some of the moes present in the uct. However, the soun attenuation perceive by an observer in the far-fiel is arguably the most important effect to be achieve. These issues may be illustrate by consiering the raiation of soun from a cylinrical uct with impeance varying circumferentially, axially or in both irections an the acoustic pressure at the far-fiel in each situation. The raiation of soun from an open pipe is represente by a pressure istribution on a isk, viz. the exit plane. The raiation of soun in free space, i.e. without obstacles, is specifie by a Kirchhof integral. The source istribution on the uct exit plane, which allows the evaluation of raiation integrals, is specifie by the raial, axial, an circumferential moes, in the cylinrical nozzle. The evaluation of the raiation integrals show that (i) the total acoustic fiel consists of a spherical wave multiplying a sum of raial moes n = 1,..., an azimuthal moes of o orer only m = 1,3,5,...; (ii) each moe consists of a monopole term an a ipole term an epen on the frequency an the raial wavenumbers, etermine by the bounary conition at the uct wall. When the impeance istribution varies circumferentially or axially, the evaluation of the wavenumbers involves the etermination of the roots of an infinite eterminant, while when the impeance aries both axially an circumferentially, the roots of a oubly infinite eterminant have to be calculate. In the case of external noise of an aircraft, the observer is on the groun, at a istance much greater than the uct iameter, an the raiation integrals for an observer in the far-fiel can be simplifie, since the ipole term is weak. The evaluation of the acoustic pressure for an observer in the far-fiel, shows that it epens on the raial wavenumbers in the nozzle, which are specifie by the wall bounary conitions, an thus epen on the acoustic impeance istribution. This allows comparison of har-walle nozzles, with liners with constant impeance an non-uniform liners, the latter with impeance istribution varying circumferentially, axially or in both irections. The far-fiel acoustic pressure is specifie by a irectivity factor that is calculate for several values of the parameters for each of the cases mentione above. 1 Introuction In a cylinrical nozzle noise can be absorbe in its interior by vortical flow [1 1] an at the walls by acoustic liners, which may have uniform [11 16] or non-uniform [17 24] impeances. The soun fiel receive by an observer in the farfiel is etermine by raiation out of the open 1

2 L. M. B. C. CAMPOS, J. M. G. S. OLIVEIRA en of the nozzle [25 32] or inlet of a fan [33,34]. In the case of nozzle there is a refraction effect which can be represente either by a vortex sheet [35] or an irregular shear layer [1,2] issuing from the lip. The aim of the present paper is to relate the acoustic fiel raiate through the nozzle exit to the far-fiel to the acoustic moes in a nozzle with non-uniform wall impeance. The raiation from the isk in the nozzle exit plane, to an observer in the far-fiel, consists, to leaing orer, of monopole an ipole terms. The acoustic pressure istribution in the nozzle exit plane, is expresse in terms of uct moes, allowing the evaluation of raiation integrals. The latter involve the acoustic eigenfunction, corresponing to the eigenvalues for propagating or cut-on moes, in the case of rigi walls or walls with uniform impeance, an also walls with impeance varying circumferentially, axially or in both irections. A numerical application is mae, incluing the calculation of raial eigenvalues an plotting o acoustic eigenfunctions for uniform liner uct in comparison with weakly an strongly circumferentially non-uniform liners; the benefits of the non-uniform lining are assesse, not incluing the effects of transmission across the irregular an turbulent shear layer issuing from the jet nozzle lip. 2 Soun raiation from an open pipe The raiation of soun from an open pipe (Fig. 1) is represente by a pressure istribution on a isk (2.1), viz. the exit plane; the raiation integrals are simplifie (2.3) for an observer in the far-fiel (2.2). 2.1 Source istribution on a circular isk The raiation of soun in free space, i.e. without obstacles, is specifie by the Kirchhof integral: p(x) = 1 4π D 1 x y e iω(t x y /c) q(y) 3 y, (1) for a spatial source istribution of strength q at position y in the omain D, with frequency ω, y r α z Fig. 1 Soun raiation from the isk on the exit plane of a cylinrical pipe to an observer in the far-fiel an raiation to the observer at x, in a homogeneous meium at rest, for which the soun spee is c. In the case of a isk of raius a on the XOY - plane with centre at the origin, the position of the source is written in polar coorinates (R,α): α 2π, R a; D θ ρ y = R(e x cosα + e y sinα), (2) an the position of the observer in spherical coorinates (r, θ, ϕ): x = r(e x cosθ + e z sinθ), (3) where ϕ = because the XOZ-plane can be taken through the observer. The raiation integral (1) becomes in this case p(r,θ,t) = e iωt 4π 2π R α R D(R,α) ei(ω/c)d(r,α) q(r,α) (4) where q(r,α) is the source istribution, an D(R,α) = x y = R 2 +r 2 2Rr sinθcosα 1/2, (5) is the istance between observer an source. 2.2 Reception by an observer in the far-fiel In the case of external noise of an aircraft, the observer is on the groun, at a istance much P x 2

3 On soun raiation from an open-ene non-uniformly line cylinrical nozzle greater than the nozzle iameter, R 2 a 2 << r 2 : D(R,α) = r Rsinθcosα + O(R 2 /r 2 ), (6a) an the istance may be simplifie by (6a), an the inverse istance by: 1/D(R,α) = 1/r + R/r 2 sinθcosα + O(R 2 /r 2 ). (6b) Substitution of (6a, 6b) in the raiation integral (4) specifies p(r,θ,t) = e iωt 4π R r 2π α R q(r,α) [1 + Rr ] sinθcosα e i(ω/c)(r Rsinθcosα) (7) the acoustic fiel receive by the observer in the far-fiel. 2.3 Decomposition into monopole an ipole terms The acoustic pressure receive in the far-fiel (7) may be written p(r,θ,t) = eiω(r/c t) (I 1 + I 2 ), (8) 4πr as a spherical wave raiate from the origin (or isk or nozzle centre) to the observer, multiplie by monopole (9a) an ipole (9b) terms I 1 = I 2 = 2π 2π α Rq(R,α)Re i(ωr/c)sinθcosα) (9a) α R q(r,α) Re i(ωr/c)sinθcosα) (R/r)sinθcosα where the latter is relate to the former by I 2 = 1 r (9b) (iω/c) I 1, (1) since ifferentiation with regar to iω/c is equivalent to multiplication by Rsinθcosα. Substitution of (9a) an (1) in (8) yiels: { } p(r,θ,t) = eiω(r/c t) 1 + i 4πr (ωr/c) 2π α Rq(R,α)Re i(ωr/c)sinθcosα, (11) as the acoustic pressure receive by an observer in the far-fiel, from a source istribution q(r, α) on a isk. The latter is specifie next in terms of the acoustic moes of a nozzle. 3 Moal structure in a cylinrical nozzle The source istribution on the nozzle exit plane (3.2), which allows the evaluation of raiation integrals (3.3), is specifie by the raial, axial, an circumferential moes, in the cylinrical nozzle. 3.1 Raial, axial an circumferential moes The moes in a cylinrical nozzle are specifie by the solution of the classical wave equation in cylinrical coorinates: { 1 R R R R R 2 α z } c 2 t 2 Q =, (12) in the absence of mean flow. Since the coefficients of the wave equation epen only on raius R, it is convenient to use a Fourier ecomposition in (z,α,t), viz.: Q(R,α,z,t) = e iωt + + e imα k e ikz P m (R,k), m= (13) a wave of frequency ω, with longituinal wavenumber k an circumferential wavenumber m. Substitution of (13) in (12) shows that: [R 2 2 R 2 + R R + κ2 R 2 m 2 ] P m (R,k) =, (14) the raial epenence is specifie, for a cylinrical nozzle, by a Bessel function of orer m: P m (R,k) = J m (κr), (15a) 3

4 L. M. B. C. CAMPOS, J. M. G. S. OLIVEIRA with raial wavenumber κ 2 = (ω/c) 2 k 2, (15b) specifie by a bounary conition at the nozzle wall R = a, which specifies the raial moes κ nm with n = 1,...,, which may be istinct for each circumferential orer m. 3.2 Amplitues of soun generation in a nozzle Thus the acoustic fiel in the nozzle consists of a superposition (13, 15a) of: Q(R,α,z,t) = e iωt + m= e imα J m (κ mn R)e ikmnz, (16) n=1 where m is the azimuthal an n the raial orer, the raial wavenumbers κ mn are specifie by the bounary conition at the nozzle wall R = a, the axial wavenumbers are relate by (15b), viz.: k mn = (ω/c) 2 (κ mn ) 2, (17) an the amplitues of each moe are specifie by the source istribution in the nozzle. The pressure istribution q(r,α)e iωt in the nozzle exit plane z = is thus specifie by: q(r,α) = Q(R,α,,) = + e imα J m (κ mn R), (18) m= n=1 which may be substitute in the raiation integral (11) to specify the soun fiel receive by the observer in the far-fiel. Substitution of (17) into (11) shows that the soun fiel receive by the observer in the far-fiel p(r,θ,t) = eiω(r/c t) 4πr + m= n=1 p mn (r,θ), (19a) is specifie by a superposition of raiation integrals for each moe: { } p mn (r,θ) = 1 i (ωr/c) 2π α RRe i(ωr/c)sinθcosα e imα J m (κ mn R), which are evaluate next. (19b) 3.3 Evaluation of raiation integrals for each moe The acoustic pressure in the far-fiel for the n-th raial an m-th azimuthal moe is given by an integration in R: { } p mn (r,θ) = 1 + i (ωr/c) RRJ m (κ mn R)I m ((ωr/c)sinθ), where the integration in α appears in ζ = (ωr/c)sinθ : I m (ζ) = 2π The integral may be re-written: e imα+iζcosα α. (2a) (2b) π ( I m (ζ) = e imα e iζcosα + ( 1) m e iζcosα) α π [ = cos(mα) e iζcosα + ( 1) m e iζcosα] α, which is expressible in terms of Bessel functions [36]: by: J m (ζ) = i m π π e izcosα cos(mα)α, (21) I m (ζ) = πi m [J m (ζ) + ( 1) m J m ( ζ)]. (22) The series expansion for the Bessel functions J m (ζ) = (ζ/2) m { ( ζ 2 /4) s /(s!(s + m)!) }, s= (23) 4

5 On soun raiation from an open-ene non-uniformly line cylinrical nozzle implies that so that J m ( ζ) = ( ) m J m (ζ), (24) I m (ζ) = 2πi m J m (ζ). (25) 4 Acoustic effects of non-uniform impeance The evaluation of the acoustic pressure for an observer in the far-fiel (4.1), shows that it epens on the raial wavenumbers in the nozzle, which are specifie by the wall bounary conitions (4.2). This allows comparison of har-walle nozzles, with liners with constant impeance an non-uniform liners, the latter with impeance istribution varying circumferentially, axially or in both irections. 4.1 Selection of cut-off an cut-on moes When substituting equations (25) in (19b): { } p mn (r,θ) = 1 + i (ωr/c) J m (κ mn R)J m ((ωr/c)sinθ)rr, (26) the property [37] of Bessel functions: implies that J m (ζ) = ( ) m J m (ζ), (27) p mn (r,θ) = p m,n (r,θ), (28) so that the total pressure fiel (19a) simplifies to: p(r,θ,t) = eiω(r/c t) 2πr n=1 + m= The secon term of (26) involves: (ωr/c) J m((ωr/c)sinθ) = R r p mn (r,θ). (29) (ωr/c) J m((ωr/c)sinθ) = R r sinθj m((ωr/c)sinθ), (3) where prime enotes erivative of the Bessel function with regar to its argument: J m(ζ) = J m 1 (ζ) m ζ J m(ζ). (31) Substitution of (31) in (3) yiels (ωr/c) J m((ωr/c)sinθ) = ( ) R ωr r sinθj m 1 c sinθ mc ωr J m ( ) ωr c sinθ, (32) which specifies the ipole term in the acoustic pressure receive in the far-fiel: i a r p mn (r,θ) = a 2 1 { J m (κ mn as) J m (sωsinθ) [ m ] } Ω J m(sωsinθ) sinθj m 1 (sωsinθ) ss, (33) where a imensionless raial istance (34a) an a imensionless frequency (34b) were introuce: s = R/a Ω = ωa/c. (34a) (34b) Thus: (i) the total acoustic fiel (29) consists of a spherical wave multiplying a sum of raial moes n = 1,..., an azimuthal moes of orer m =,1,..., ; (ii) each moe (33) consists of a monopole term (first in the curly brackets) an a ipole term (in square brackets); (iii) the parameters are the imensionless frequency (34b) calculate from the nozzle raius a, the ratio of the latter to the istance of the observer a/r in the ipole term (which is weak because a 2 << r 2 for observers in the far-fiel), an κ mn a which are the raial wavenumbers, etermine by the bounary conition at the uct wall. 4.2 Rigi, impeance an non-uniform walls The simplest case (I) is a nozzle with rigi walls, for which the normal velocity at the wall is zero, 5

6 L. M. B. C. CAMPOS, J. M. G. S. OLIVEIRA implying from the momentum equation that the normal erivative of the pressure is zero: 1 p = iωv n (R = a) = ρ r, (35) r=a where ρ is the mass ensity. Thus: = J m(κ mn a), (36) so that the raial wavenumbers κ mn are etermine by the zeros j mn of the erivative of the Bessel function J m : J m( j mn ) = : κ mn = j mn /a. (37) Since these zeros are real, the raial wavenumbers κ mn are real, an the corresponing (17) axial wavenumbers: k mn a = Ω 2 ( j mn ) 2, (38) are: (i) either real, for propagating or cut-on moes, if j mn Ω; (ii) or imaginary, for evanescent or cut-off moes, if j mn > Ω. Since the zeros of the erivative of the Bessel function J m form an unboune sequence j m1, j m2,..., there is a finite number of cut-on moes, larger for larger imensionless frequency. The cut-off moes make a negligible contribution to raiation to the far-fiel, so the sum in the total acoustic pressure (29) is restricte to cut-on moes. The amplitue of the latter cut-on moes has been foun to be weakly epenent on moe orer for turbomachinery noise, an thus may be taken as a constant factor, an omitte together with the spherical wave term e iω(r/c t) /2πr an a 2, which are common factors, regarless of the wall conition, an thus o not affect the comparison between rigi line walls. For a rigi wall, the farfiel acoustic pressure is thus specifie by a irectivity factor: P(θ) = m= j mn <Ω 1 J m ( j mn s)j m (sωsinθ)s s, n=1 (39) where only the monopole term was consiere, since it ominates the ipole term. In the case (II) of a wall with uniform impeance Z, the raial wavenumbers are specifie by the roots of: iz J m(κ mn a) = J m (κ mn a), (4) where Z is the specific impeance, i. e. impeance ivie by that of a plane wave: the Z = Z /ρc. (41) In this case the raial wavenumbers κ mn are generally complex, an the axial wavenumbers (17) also: k mn a = Ω 2 (κ mn a) 2. (42) Although the istinction between cut-on an cutoff moes is not so clear in the case of an impeance wall, the sum for the total acoustic fiel is taken for the cut-on moes of a rigi wall, as in (39), but for the complex raial wavenumbers, so that (39) remains vali, an can be evaluate (cfr. [36]). The acoustic pressure receive in the far-fiel is: p mn (θ) = 2ik mn a [ (κ mn a) 2 (Ωsinθ) 2] 1 [ Ω sinθjmn (κ mn a)j m(ωsinθ) κ mn aj m (Ωsinθ)J (κ mn a) ] (43) We procee to consier the calculation of raial eigenvalues k mn in the case of liners of circumferentially or axially non-uniform impeance or both. 4.3 Axially an /or circumferentially nonuniform liners This expression (43) can also be use in the case of non-uniform wall impeance Z(θ, z), except that raial wavenumbers are no longer etermine as the roots of (4). In the case of a circumferentially non-uniform istribution, which is represente by the first two terms of a Fourier series: Z A (θ) = Z (1 + 2εcosθ), (44) the raial wavenumbers are specifie exactly by the roots of an infinite eterminant, [ et i Ω ] κa J m(κa)δ mm Z m mj m(κa) =, (45) 6

7 On soun raiation from an open-ene non-uniformly line cylinrical nozzle where all the impeance Fourier coefficients Z m = except Z an Z ±1 = εz. In the case of an axially non-uniform wall impeance over a length L of uct, represente by the first two terms of a Fourier series: Z B (z) = Z [1 + 2δcos(2πz/L)], (46) the raial wavenumbers are specifie exactly by the roots of an infinite eterminant, [ ] et i 1 + 2πl κa J m(κa)δ ll Z l l J m(κa) =, (47) where all the impeance Fourier coefficients Z l = except Z an Z ±1 = δz. In the case of impeance istribution varying both axially an raially, represente by the first two terms of the ouble Fourier series: Z C (θ,z) = Z (1 + εcosθ)[1 + δcos(2πz/l)], (48) the raial wavenumbers are specifie exactly by the roots of a ouble infinite eterminant, i.e. an infinite eterminant whose terms are infinite eterminants in a combination of (45) an (47). The irectivity of the acoustic pressure in the far-fiel is specifie by (43) also in the case of non-uniform impeance, except that the raial wavenumbers are given by the roots of (45) for circumferentially non-uniform impeance (44), the roots of (47) for axially nonuniform impeance (46), an roots of a combination of (45) an (47) for impeance varying both axially an circumferentially (48). 5 Noise reuction for far-fiel observer The eigenvalues an eigenfunctions are calculate for circumferentially non-uniform liners, to assess the noise reuction benefit relative to uniform liners. 5.1 Eigenvalues for the imensionless raial wavenumbers Arguably the best measure of the effectiveness of an acoustic liner is the effect on the reuction of far-fiel noise. As a numerical example consier a nozzle of raius a = 1m, an a wave frequency f = 1kHz or ω = 2π f = s 1, corresponing, for a soun spee c = 34ms 1, to a imensionless frequency or Helmholtz number Ω = ωr/c = The conition j mn < Ω for the roots j mn of the erivatives of the Bessel functions J ( j mn ) =, specifies the propagating or cut-on moes. The raial wavenumbers can be etermine using (4) for a uniform specific impeance impeance: Z = 2.5 i.4, (49) an (45), for circumferentially non-uniform impeances (44), with relative amplitue of the harmonic ε =.1 +.1i, ε =.3i. (5a) (5b) The former is esignate comparatively the weakly (5a) an the latter the strongly (5b) non uniform liner, although in both cases the impeance is the same (49) on the mean, an the variation from the mean is small in absolute terms in both cases. 5.2 Soun attenuation ue to non-uniform liner In orer to assess the soun attenuation ue to the non-uniform liner, the eigenfunctions for the acoustic pressure (39) with non-uniform p mn (θ;z,ε) an uniform p mn (θ;z,) liner are compare in a logarithmic scale: p (θ,z,ε) mn (θ;z,ε) log 1 p mn (θ,z,), (51) in ecibels; the effect on energy woul be the ouble of that given by (51), an is plotte as a function of azimuthal angle θ in Figs. 2 an 3. In Fig. 2 the axisymmetric moe m = is consiere for increasing raial orers n = 1,...,6. The funamental raial moe n = 1 (top left) is not much affecte by the non-uniform impeance, but the effect becomes visible for the secon 7

8 L. M. B. C. CAMPOS, J. M. G. S. OLIVEIRA a.4 b c e.4 f Fig. 2 Comparison on a ecimal logarithmic scale (51) of the ratio of amplitues of the acoustic pressure (43), for a non-uniform (44) an a uniform (49) liner with the same mean impeance, in the cases of weak (- -) an strong ( ) non-uniformity, as a function of azimuthal angle θ, for axisymmetric moes m = an (a) the funamental raial moe n = 1 an the harmonics (b) n = 2, (c) n = 3, () n = 4, (e) n = 5, an (f) n = 6. n = 2 (top right) an thir n = 3 (mile left) raial moes. Significant amplitue variation occur for the fourth n = 4 (mile right), fifth n = 5 (bottom left), an sixth n = 6 (bottom right) raial moes; these amplitue variations are mostly reuctions, although some increases occur. Broaly similar conclusions can be rawn from Fig. 3, which consiers the funamental raial moe n = 1 for several azimuthal orers m =,2,...,14. The axisymmetric moe m = 8

9 On soun raiation from an open-ene non-uniformly line cylinrical nozzle a.1 b c e.1 f g.1 h Fig. 3 As Fig. 2, for the funamental axial moe n = 1 an (a) axisymmetric m =, an nonaxisymmetric moes of increasing even azimuthal orer (b) m = 2, (c) m = 4, () m = 6, (e) m = 8, (f) m = 1, (g) m = 12, an (h) m = 14. 9

10 L. M. B. C. CAMPOS, J. M. G. S. OLIVEIRA (top left) is not much affecte by the non-uniform liner. The azimuthal moes m = 2 (top right) an m = 4 (secon row, left) show visible reuctions. The higher orer azimuthal moes m = 6 (thir row, left), m = 1 (thir row, right), m = 12 (bottom left), an m = 14 (bottom right) show a clear amplitue reuction ue to the strongly nonuniform liner over most irections, although an increase is seen for some irections. It is clear that a non-uniform liner, with relative impeance variations in the range 2 3%, can provie an overall soun reuction. 6 Discussion The effect of non-uniform impeance on soun raiation to the far-fiel has been assesse on the basis of the acoustic pressure on the isk corresponing to the nozzle exit plane, without accounting for ege iffraction effects. The latter can be quite important [35], in particular in the presence of a mean flow, when a turbulent an irregular shear layer [1,2] is issue from the nozzle lip. The latter causes spectral an irectional broaening of soun, which further reuces the noise levels. Although the latter effects have not been moelle here, it is clear that an attenuation of the input soun fiel incient on the shear layer from the interior of the jet will result in a further attenuate soun fiel transmitte to an observer in the far-fiel outsie the shear layer. The attenuations shown for the basic in-uct soun fiel are clear even though: (i) the non-uniform liner (44) has the same average impeance as the uniform liner; (ii) only one liner impeance harmonic was use, involving a single parameter ε, which was not optimize. By consiering multi-harmonic impeance istribution ] Z(θ) = Z [1 + L l=1 ε l cos(lθ), (52) an optimizing the parameter ε l, greater noise reuction coul be obtaine. In the present example a relative impeance variation ε of 2 3% was shown to be sufficient to affect significantly the soun fiels. References [1] L. M. B. C. Campos. The spectral broaening of soun in turbulent shear layers. Part 1: The transmission of soun through turbulent shear layers. Journal of Flui Mechanics, 89: , [2] L. M. B. C. Campos. The spectral broaening of soun in turbulent shear layers. Part 2: The spectral broaening of soun an aircraft noise. Journal of Flui Mechanics, 89: , [3] H. E. Plumblee an P. E. Doak. Duct noise raiation through a jet flow. Journal of Soun an Vibration, 65(4): , [4] P. A. Nelson an C. L. Morfey. Aeroynamic soun prouction in low spee flow ucts. Journal of Soun an Vibration, 79(2): , [5] M. S. Howe. The amping of soun by turbulent wall shear layers. Journal of the Acoustical Society of America, 98: , [6] C. K. W. Tam an L. Auriault. The wave moes in ucte swirling flows. Journal of Flui Mechanics, 371:1 2, [7] L. M. B. C. Campos an P. G. T. A. Serrão. On the acoustics of an exponential bounary layer. Philosophical Transactions of the Royal Society of Lonon, Series A, 356: , [8] L. M. B. C. Campos, J. M. G. S. Oliveira, an M. H. Kobayashi. On soun propagation in a linear flow. Journal of Soun an Vibration, 219(5):739 77, [9] L. M. B. C. Campos an M. H. Kobayashi. On the reflection an transmmission of soun in a thick shear layer. Journal of Flui Mechanics, 42:1 24, 2. [1] M. Willatzen. Soun propagation in a moving flui confine by cilinrical walls exact series solutions for raially epenent flow profiles. Acustica, 87: , 21. [11] A. D. Rawlins. Raiation of soun from an unflange rigi cylinrical uct with an acoustically absorbing internal surface. Proceeings of the Royal Society of Lonon, Series A, 361:65 91, [12] K. Ogimoto an G. W. Johnston. Moal raiation impeances for semi-infinite unflange circular ucts incluing flow effects. Journal of 1

11 On soun raiation from an open-ene non-uniformly line cylinrical nozzle Soun an Vibration, 62(4):598 65, [13] W. Koch an W. Möhring. Eigensolutions for liners in uniform mean flow ucts. AIAA Journal, 21(2):2 213, [14] W. Koch. Raiation of soun from a twoimensional acoustically line uct. Journal of Soun an Vibration, 55(2): , [15] S. W. Rienstra. Contributions to the theory of soun propagation in ucts with bulk-reacting lining. Journal of the Acoustical Society of America, 77(5): , May [16] D. A. Bies, C. H. Hansen, an G E Briges. Soun propagation in rectangular an circular cross-section ucts with flow an bulkreacting liner. Journal of Soun an Vibration, 146(1):47 8, [17] M. S. Howe. The attenuation of soun in a ranomly line uct. Journal of Soun an Vibration, 87(1):83 13, [18] W. R. Watson. An acoustic evaluation of circumferentially segmente uct liners. AIAA Journal, 22(9): , September [19] C. R. Fuller. Propagation an raiation of soun from flange circular ucts with circumferentially varying wall amittances, I: semiinfinite ucts. Journal of Soun an Vibration, 93(3):321 34, [2] C. R. Fuller. Propagation an raiation of soun from flange circular ucts with circumferentially varying wall amittances, II: finite ucts with sources. Journal of Soun an Vibration, 93(3): , [21] P. G. Vaiya. The propagation of soun in ucts line with circumferentially non-uniform amittance of the form η + η q exp(iqθ). Journal of Soun an Vibration, 1(4): , [22] B. Regan an J. Eaton. Moelling the influence of acoustic liner non-uniformities on uct moes. Journal of Soun an Vibration, 219(5): , [23] L. M. B. C. Campos an J. M. G. S. Oliveira. On the acoustic moes in a cylinrical nozzle with an arbitrary impeance istribution. In Proceeings of the 23 r ICAS Congress, 22. [24] L. M. B. C. Campos an J. M. G. S. Oliveira. On the optimization of non-uniform acoustic liners on annular nozzles. To be publishe, 23. [25] S. D. Savkar. Raiation of cylinrical uct acoustic moes with flow mismatch. Journal of Soun an Vibration, 42: , [26] W. Koch. Attenuation of soun in multi-element acoustically line rectangular ucts in the absence of mean flow. Journal of Soun an Vibration, 52(4): , [27] S. W. Rienstra. Acoustic raiation from a semiinfinite uct in a uniform subsonic mean flow. Journal of Soun an Vibration, 94(2): , [28] R. Martinez. Diffracting open-ene pipe treate as a lifting surface. AIAA Journal, 26(4):396 44, [29] A. Snakovska, H. Iczak, an B. Bogusz. Moal analysis of the acoustic fiel raiate from an unflange cylinrical uct theory an measurement. Acustica, 82:21 26, [3] A. Snakowska an H. Iczak. Preiction of multitone soun raiation from a circular uct. Acustica, 83: , [31] P. Joseph an C. L. Morfey. Multimoe raiation from an unflange, semi-infinite circular uct. Journal of the Acoustical Society of America, 15(5):259 26, May [32] S. T. Hocter. Exact an approximate irectivity patterns of the soun raiate from cylinrical uct. Journal of Soun an Vibration, 227(2):397 47, [33] S. J. Horowitz, R. K. Sigman, an B. T. Zinn. An iterative metho for preicting turbofan inlet acoustics. AIAA Journal, 2(12): , December [34] H. D. Meyer. Effect of inlet reflections on fan noise raiation. AIAA Journal, 34(9): , September [35] R. M. Munt. Acoustic raiation from a semiinfinite circular uct in an uniform subsonic mean flow. Journal of Flui Mechanics, 83:69 64, [36] M. Abramowitz an I. A. Stegun. Hanbook of Mathematical Functions. Dover, [37] E. T. Whittaker an G. N. Watson. A Course of Moern Analysis. Cambrige University Press,

Study on aero-acoustic structural interactions in fan-ducted system

Study on aero-acoustic structural interactions in fan-ducted system Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Duct spinning mode s particle velocity imaging with in-duct circular microphone array

Duct spinning mode s particle velocity imaging with in-duct circular microphone array AIAA Aviation 16-2 June 214, Atlanta, GA 2th AIAA/CEAS Aeroacoustics Conference AIAA 214-3112 Duct spinning moe s particle velocity imaging with in-uct circular microphone array Qingkai Wei Peking University,

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

On the use of leaky modes in open waveguides for the sound propagation modeling in street canyons

On the use of leaky modes in open waveguides for the sound propagation modeling in street canyons On the use of leaky moes in open waveguies for the soun propagation moeling in street canyons Arien Pelat, a Simon Félix, an Vincent Pagneux LAUM, CNRS, Université u Maine, avenue Olivier Messiaen, 7285

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

RETROGRADE WAVES IN THE COCHLEA

RETROGRADE WAVES IN THE COCHLEA August 7, 28 18:2 WSPC - Proceeings Trim Size: 9.75in x 6.5in retro wave 1 RETROGRADE WAVES IN THE COCHLEA S. T. NEELY Boys Town National Research Hospital, Omaha, Nebraska 68131, USA E-mail: neely@boystown.org

More information

Ultra-thin Acoustic Metasurface-Based Schroeder Diffuser

Ultra-thin Acoustic Metasurface-Based Schroeder Diffuser Ultra-thin Acoustic Metasurface-Base Schroeer Diffuser Yifan Zhu, Xuong Fan, Bin Liang *, Jianchun Cheng *, an Yun Jing * Key Laboratory of Moern Acoustics, MOE, Institute of Acoustics, Department of Physics,

More information

P. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom

P. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

Improvements of a parametric model for fan broadband and tonal noise

Improvements of a parametric model for fan broadband and tonal noise Improvements of a parametric model for fan broadband and tonal noise A. Moreau and L. Enghardt DLR - German Aerospace Center, Mueller-Breslau-Str. 8, 10623 Berlin, Germany antoine.moreau@dlr.de 4083 Engine

More information

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2)

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2) 1.4.3 Elementary solutions to Laplace s equation in the spherical coorinates (Axially symmetric cases) (Griffiths 3.3.) In the spherical coorinates (r, θ, φ), the Laplace s equation takes the following

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations

Optimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,

More information

Computational Aeroacoustics

Computational Aeroacoustics Computational Aeroacoustics Duct Acoustics Gwénaël Gabard Institute of Sound and Vibration Research University of Southampton, UK gabard@soton.ac.uk ISVR, University of Southampton, UK ERCOFTAC Computational

More information

1 Boas, p. 643, problem (b)

1 Boas, p. 643, problem (b) Physics 6C Solutions to Homework Set #6 Fall Boas, p. 643, problem 3.5-3b Fin the steay-state temperature istribution in a soli cyliner of height H an raius a if the top an curve surfaces are hel at an

More information

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010 NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

The Three-dimensional Schödinger Equation

The Three-dimensional Schödinger Equation The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

PHY 114 Summer 2009 Final Exam Solutions

PHY 114 Summer 2009 Final Exam Solutions PHY 4 Summer 009 Final Exam Solutions Conceptual Question : A spherical rubber balloon has a charge uniformly istribute over its surface As the balloon is inflate, how oes the electric fiel E vary (a)

More information

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL

EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Research Article Numerical Analysis of Inhomogeneous Dielectric Waveguide Using Periodic Fourier Transform

Research Article Numerical Analysis of Inhomogeneous Dielectric Waveguide Using Periodic Fourier Transform Microwave Science an Technology Volume 2007, Article ID 85181, 5 pages oi:10.1155/2007/85181 Research Article Numerical Analysis of Inhomogeneous Dielectric Waveguie Using Perioic Fourier Transform M.

More information

AIAA The Webster Equation Revisited

AIAA The Webster Equation Revisited I -5 The Webster Equation Revisite S.W. Rienstra Einhoven University of Technology Einhoven, NL 8th I/CES eroacoustics Conference 7-9 June Breckenrige, Colorao For permission to copy or to republish, contact

More information

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Wave Propagation in Grounded Dielectric Slabs with Double Negative Metamaterials

Wave Propagation in Grounded Dielectric Slabs with Double Negative Metamaterials 6 Progress In Electromagnetics Research Symposium 6, Cambrige, US, March 6-9 Wave Propagation in Groune Dielectric Slabs with Double Negative Metamaterials W. Shu an J. M. Song Iowa State University, US

More information

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR 5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Sound diffraction by the splitter in a turbofan rotor-stator gap swirling flow

Sound diffraction by the splitter in a turbofan rotor-stator gap swirling flow Nationaal Lucht- en Ruimtevaartlaboratorium National Aerospace Laboratory NLR Sound diffraction by the splitter in a turbofan rotor-stator gap swirling flow R.J. Nijboer This report is based on a presentation

More information

Partial Differential Equations

Partial Differential Equations Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.

More information

AIEEE Physics Model Question Paper

AIEEE Physics Model Question Paper IEEE Physics Moel Question Paper ote: Question o. 11 to 1 an 1 to consist of Eight (8) marks each for each correct response an remaining questions consist of Four (4) marks. ¼ marks will be eucte for inicating

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

Math 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions

Math 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

V = Flow velocity, ft/sec

V = Flow velocity, ft/sec 1 Drag Coefficient Preiction Chapter 1 The ieal force acting on a surface positione perpenicular to the airflow is equal to a ynamic pressure, enote by q, times the area of that surface. Dynamic pressure

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical

More information

S10.G.1. Fluid Flow Around the Brownian Particle

S10.G.1. Fluid Flow Around the Brownian Particle Rea Reichl s introuction. Tables & proofs for vector calculus formulas can be foun in the stanar textbooks G.Arfken s Mathematical Methos for Physicists an J.D.Jackson s Classical Electroynamics. S0.G..

More information

6 Wave equation in spherical polar coordinates

6 Wave equation in spherical polar coordinates 6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Section 2.7 Derivatives of powers of functions

Section 2.7 Derivatives of powers of functions Section 2.7 Derivatives of powers of functions (3/19/08) Overview: In this section we iscuss the Chain Rule formula for the erivatives of composite functions that are forme by taking powers of other functions.

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

A Model of Electron-Positron Pair Formation

A Model of Electron-Positron Pair Formation Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.

More information

The continuity equation

The continuity equation Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal

More information

Effect of Rotation on Thermosolutal Convection. in a Rivlin-Ericksen Fluid Permeated with. Suspended Particles in Porous Medium

Effect of Rotation on Thermosolutal Convection. in a Rivlin-Ericksen Fluid Permeated with. Suspended Particles in Porous Medium Av. Theor. Appl. Mech., Vol. 3,, no. 4, 77-88 Effect of Rotation on Thermosolutal Convection in a Rivlin-Ericksen Flui Permeate with Suspene Particles in Porous Meium A. K. Aggarwal Department of Mathematics

More information

Electromagnetism Answers to Problem Set 3 Spring Jackson Prob. 2.1: Charge above a grounded plane. (z d)

Electromagnetism Answers to Problem Set 3 Spring Jackson Prob. 2.1: Charge above a grounded plane. (z d) Electromagnetism 76 Answers to Problem Set 3 Spring 6. Jackson Prob..: Charge above a groune plane (a) Surface charge ensity E z (ρ, z) = φ z = q [ (z ) [ρ + (z ) (z + ) 3/ [ρ + (z + ) 3/ Evaluate E z

More information

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie

More information

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor DOI: 10.24352/UB.OVGU-2017-110 TECHNISCHE MECHANIK, 37, 2-5, (2017), 347-357 submitte: June 15, 2017 Simple Electromagnetic Motor Moel for Torsional Analysis of Variable Spee Drives with an Inuction Motor

More information

MAE 210A FINAL EXAM SOLUTIONS

MAE 210A FINAL EXAM SOLUTIONS 1 MAE 21A FINAL EXAM OLUTION PROBLEM 1: Dimensional analysis of the foling of paper (2 points) (a) We wish to simplify the relation between the fol length l f an the other variables: The imensional matrix

More information

Single Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design

Single Arm, Centrifugal, Water Turbine for Low Head and Low Flow Application: Part 1- Theory and Design Energy an Power 2018, 8(2): 51-55 DOI: 10.5923/j.ep.20180802.03 Single Arm, Centrifugal, Water Turbine for Low ea an Low Flow Application: Part 1- Theory an Design Kiplangat C. Kononen 1, Augustine B.

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

PH 132 Exam 1 Spring Student Name. Student Number. Lab/Recitation Section Number (11,,36)

PH 132 Exam 1 Spring Student Name. Student Number. Lab/Recitation Section Number (11,,36) PH 13 Exam 1 Spring 010 Stuent Name Stuent Number ab/ecitation Section Number (11,,36) Instructions: 1. Fill out all of the information requeste above. Write your name on each page.. Clearly inicate your

More information

VIBRATIONS OF A CIRCULAR MEMBRANE

VIBRATIONS OF A CIRCULAR MEMBRANE VIBRATIONS OF A CIRCULAR MEMBRANE RAM EKSTROM. Solving the wave equation on the isk The ynamics of vibrations of a two-imensional isk D are given by the wave equation..) c 2 u = u tt, together with the

More information

Approaches for Predicting Collection Efficiency of Fibrous Filters

Approaches for Predicting Collection Efficiency of Fibrous Filters Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Dusty Plasma Void Dynamics in Unmoving and Moving Flows

Dusty Plasma Void Dynamics in Unmoving and Moving Flows 7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological

More information

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9 Oregon State University PH 213 Spring Term 2018 Prep 1 Suggeste finish ate: Monay, April 9 The formats (type, length, scope) of these Prep problems have been purposely create to closely parallel those

More information

Microwave Reflection from the Region of Electron Cyclotron Resonance Heating in the L-2M Stellarator )

Microwave Reflection from the Region of Electron Cyclotron Resonance Heating in the L-2M Stellarator ) Microwave Reflection from the Region of Electron Cyclotron Resonance Heating in the L-2M Stellarator German M. BATANOV, Valentin D. BORZOSEKOV, Nikolay K. KHARCHEV, Leoni V. KOLIK, Eugeny M. KONCHEKOV,

More information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012 arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,

More information

Planar sheath and presheath

Planar sheath and presheath 5/11/1 Flui-Poisson System Planar sheath an presheath 1 Planar sheath an presheath A plasma between plane parallel walls evelops a positive potential which equalizes the rate of loss of electrons an ions.

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

Physics 504, Lecture 7 Feb 14, Energy Flow, Density and Attenuation. 1.1 Energy flow and density. Last Latexed: February 11, 2011 at 10:54 1

Physics 504, Lecture 7 Feb 14, Energy Flow, Density and Attenuation. 1.1 Energy flow and density. Last Latexed: February 11, 2011 at 10:54 1 Last Latexe: February, at :54 Physics 54, Lecture 7 Feb 4, Energy Flow, Density an ttenuation We have seen that there are iscrete moes for electromagnetic waves with E, B e ikz it corresponing, for real

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

Cut-on, cut-off transition of sound in slowly varying flow ducts

Cut-on, cut-off transition of sound in slowly varying flow ducts Cut-on, cut-off transition of sound in slowly varying flow ducts Sjoerd W. Rienstra 19-3-1 Department of Mathematics and Computing Science, Eindhoven University of Technology, The Netherlands, s.w.rienstra@tue.nl

More information

The Principle of Least Action and Designing Fiber Optics

The Principle of Least Action and Designing Fiber Optics University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic

More information

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ Department of Physics University of arylan College Park, arylan PHYSICS 4 Prof. S. J. Gates Fall 5 Final Exam Dec. 6, 5 This is a OPEN book examination. Rea the entire examination before you begin to work.

More information

Students need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling.

Students need encouragement. So if a student gets an answer right, tell them it was a lucky guess. That way, they develop a good, lucky feeling. Chapter 8 Analytic Functions Stuents nee encouragement. So if a stuent gets an answer right, tell them it was a lucky guess. That way, they evelop a goo, lucky feeling. 1 8.1 Complex Derivatives -Jack

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Mathematical Methods wks 5,6: PDEs

Mathematical Methods wks 5,6: PDEs Mathematical Methos wks 5,6: PDEs John Magorrian, magog@thphys.ox.ac.uk These are work-in-progress notes for the secon-year course on mathematical methos. The most up-to-ate version is available from http://www-thphys.physics.ox.ac.uk/people/johnmagorrian/mm.

More information

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper

More information

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

Research Article Relativistic Propagation of Linearly/Circularly Polarized Laser Radiation in Plasmas

Research Article Relativistic Propagation of Linearly/Circularly Polarized Laser Radiation in Plasmas ISRN Optics Volume 213, Article ID 642617, 8 pages http://x.oi.org/1.1155/213/642617 Research Article Relativistic Propagation of Linearly/Circularly Polarize Laser Raiation in Plasmas Sonu Sen, 1 Meenu

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II) Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution

More information

ABSTRACT. Analysis of the Initial Spray from Canonical Fire Suppression Nozzles. Professor André Marshall, Department of Fire Protection Engineering

ABSTRACT. Analysis of the Initial Spray from Canonical Fire Suppression Nozzles. Professor André Marshall, Department of Fire Protection Engineering ABSTRACT Title of Document: Analysis of the Initial Spray from Canonical Fire Suppression Nozzles Ning Ren, Master of Science, 007 Directe By: Professor Anré Marshall, Department of Fire Protection Engineering

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

Sound Transmission in an Extended Tube Resonator

Sound Transmission in an Extended Tube Resonator 2016 Published in 4th International Symposium on Innovative Technologies in Engineering and Science 3-5 November 2016 (ISITES2016 Alanya/Antalya - Turkey) Sound Transmission in an Extended Tube Resonator

More information

Longitudinal Waves in a Rotating Solid Cylinder Immersed in an Inviscid Fluid Using Chebyshev Polynomial

Longitudinal Waves in a Rotating Solid Cylinder Immersed in an Inviscid Fluid Using Chebyshev Polynomial Longituinal Waves in a Rotating Soli Cyliner Immerse in an Invisci Flui Using Chebyshev Polynomial R. Selvamani *1, P.Ponnusamy 1 Department o Mathematics, Karunya University, Coimbatore, TamilNau, Inia

More information

NUMERICAL STUDY OF THERMAL RADIATIONS AND THERMAL STRATIFICATION MECHANISMS IN MHD CASSON FLUID FLOW. and Sardar Muhammad BILAL c

NUMERICAL STUDY OF THERMAL RADIATIONS AND THERMAL STRATIFICATION MECHANISMS IN MHD CASSON FLUID FLOW. and Sardar Muhammad BILAL c NUMERICAL STUDY OF THERMAL RADIATIONS AND THERMAL STRATIFICATION MECHANISMS IN MHD CASSON FLUID FLOW Khalil Ur REHMAN b c * Noor Ul SABA b Iffat ZEHRA c Muhamma Yousaf MALIK ab an Sarar Muhamma BILAL c

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads

Vibration Analysis of Railway Tracks Forced by Distributed Moving Loads IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,

More information

SOLUTION & ANSWER FOR KCET-2009 VERSION A1 [PHYSICS]

SOLUTION & ANSWER FOR KCET-2009 VERSION A1 [PHYSICS] SOLUTION & ANSWER FOR KCET-009 VERSION A [PHYSICS]. The number of significant figures in the numbers.8000 ---- 5 an 7.8000 5 significant igits 8000.50 7 significant igits. β-ecay means emission of electron

More information

NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS

NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This

More information