ELECTROMAGNETIC FIELD IN THE PRESENCE OF A THREE-LAYERED SPHERICAL REGION

Size: px
Start display at page:

Download "ELECTROMAGNETIC FIELD IN THE PRESENCE OF A THREE-LAYERED SPHERICAL REGION"

Transcription

1 Progre In Electromagnetic Reearch, PIER 45, , 004 ELECTROMAGNETIC FIELD IN THE PRESENCE OF A THREE-LAYERED SPHERICAL REGION K. Li and S.-O. Park School of Engineering Information and Communication Univerity 58-4 Hwaam-dong, Yuung-gu, Daejeon, , South Korea H.-Q. Zhang Shaanxi Atronomical Obervatory Chinee Academy of Science Xi an, Shaanxi, , P. R. China Abtract In thi paper, the region of interet conit of the pherical dielectric earth, coated with a dielectric layer under the air. The imple explicit formula have been derived for the electromagnetic field of a vertical electric dipole and vertical magnetic dipole in the preence of three-layered pherical region, repectively. Next, baing on the above reult, the formula for the ix component of the field in the air generated by a horizontal electric dipole are derived by uing reciprocity. The computation how that the trapped urface wave of electric type can be excited efficiently and the trapped urface wave of magnetic type can not be excited when the thickne l of the dielectric layer i omewhat and atifie the condition 0 < k1 k 0 l< π. When the thickne l of the dielectric layer i omewhat and atifie the condition π k < 1 k 0 l<π, the trapped urface wave of magnetic type can be excited. Thee formula and the computation can be applied to the urface communication at the lower frequencie. 1 Introduction The Field in the Air Generated by a Vertical Electric Dipole over the Surface of the Dielectric-Coated Spherical Electrically Earth Alo with China Reearch Intitute of Radiowave Propagation, Xinxiang, Henan, , P. R. China

2 104 Li, Park, and Zhang 3 The Field in the Air Generated by a Magnetic Dipole over the Surface of the Dielectric-Coated Spherical Electrically Earth 4 The Field in the Air Generated by a Horizontal Electric Dipole over the Surface of the Dielectric- Coated Spherical Electrically Earth 5 Concluion Acknowledgment Reference 1. INTRODUCTION The old problem of the electromagnetic field generated by a vertical electric dipole and horizontal electric dipole in a three-layered planar region ha been tudied and debated by everal invetigator becaue of it ueful application in everal cae [1 5]. It i neceary to reanalyze the old problem again. Only everal year ago, the complete, imple and accurate formula were derived for the electromagnetic field generated by a vertical electric dipole in a three-layered planar region. In [6], it ha been drawn a concluion that the trapped urface wave can be induced by a vertical electric dipole in a three-layered planar region. The wave number of the trapped urface wave i between that of the wave number k 0 in the air and the wave number k 1 in the dielectric layer, which i different with that preented in [4]. The amplitude of the trapped wave attenuate a r 1/ in the r direction near the urface and attenuate exponentially a e gz in the z direction. Thi development, naturally, rekindled the interet in the pherical problem. Recently, Pan and Zhang [8] and Li et al. [9] had derived the accurate and imple formula for the electromagnetic field of a vertical electric dipole and horizontal electric dipole over the urface of a pherical conducting earth coated with a dielectric layer, repectively. When applied to practical ituation uch a the wave propagation over the dielectric-coated earth, the earth hould be regarded a an electrically dielectric phere. It i neceary to tudy the problem furthermore. In the preent tudy, the imple analytical formula for the electromagnetic field in the air are outlined when a dipole ource (vertical electric dipole, vertical magnetic dipole, or horizontal electric dipole) i located near the urface of the dielectric-coated pherical electrically earth. Firt, along uch reearch line of [8] and [9], we will preent the formula for the field of vertical electric dipole and

3 EM field in three-layered pherical region 105 vertical magnetic dipole over the urface of a dielectric-coated pherical electrically earth. Next, baing on thee reult, we will derive the formula for the ix component of the electromagnetic field generated by a horizontal electric dipole over the dielectric-coated pherical electrically earth. It hould be noted that the condition of the trapped urface wave of electric type and thoe of the trapped urface wave of magnetic type are preented clearly in thi paper. Thee formula and the computation can be applied to the urface communication at the lower frequencie. Figure 1. Geometry of a vertical electric dipole in the preence of a three-layered pherical region.. THE FIELD IN THE AIR GENERATED BY A VERTICAL ELECTRIC DIPOLE OVER THE SURFACE OF THE DIELECTRIC-COATED SPHERICAL ELECTRICALLY EARTH Firt, we will analyze the field of a vertical electric dipole in the preence of three-layered pherical region. In the preent tudy, the geometry of the problem conit of the air (region 0, r > a), the dielectric layer with it thickne l (region 1, a l < r < a), and the earth (region, r<a l), a illutrated in Figure 1. The wave number in the three region are k 0 = ω µ 0 ε 0, k 1 = ω µ 0 (ε 0 ε r1 + iσ 1 /ω),

4 106 Li, Park, and Zhang k = ω µ 0 (ε 0 ε r + iσ /ω) (1) where ε r i the relative permittivity, σ i the conductivity. Auming that the dipole i repreented by the current denity ẑidlδ(x)δ(y)δ(z b), where b = z + a, and z > 0 denote the height of the dipole above the urface of the dielectric-coated earth, and ue i made of the time dependence e ( iωt), the nonzero component of the field E r, E θ, and H ϕ, which the field i defined a the field of electric type, can be expreed a follow, ( ) E r (i) = r + k i (U i r) () E (i) θ = 1 r θ r (U ir) (3) H ϕ (i) = iωε i r θ (U ir). (4) Here the potential function U i i the olution of in the calar Helmholtz equation ( + ki )U i = 0 (5) where i = 0, 1,. In pherical coordinate, becaue of ymmetry, the expanded form (5) i written a ( 1 r r U ) i 1 + r r r in θ ( θ in θ U i θ ) + k i U i =0. (6) Uing the eparation of variable method, we aume that a olution for U i can be written in the form of U i = 1 r Z i(r)φ(θ) (7) then [ 1 d in θ Φ(θ) ] + ν(ν + 1)Φ(θ) = 0 (8) in θ dθ θ d [ ] Z i (r) dr + ki ν(ν +1) 1 ki Z i (r) = 0. (9) r The equation (8) i cloely related to the well-known Legendre differential equation whoe olution i taken a the Legendre function of firt kind P ν (co(π θ)) [8]. In Region 0, the olution of (9) ha been addreed in [8]. Z 0 (r) =AW (t y) (10)

5 EM field in three-layered pherical region 107 where W (t) i the Airy function of the econd kind, and t and y are defined a [ a 3 ] /3 [ ] ν(ν +1) t = ν(ν +1) a k0 (11) [ ] ν(ν +1) 1/3 ( ) 1/3 y = a 3 z k 0 z. (1) k 0 a In (10), r = a + z and z a, z>0denote the height above the urface of the dielectric-coated earth. We aume that both the dipole and the point of obervation cloe to the pherical interface between the air and the dielectric layer. In the whole text we take z a and z r a. In Region 1, a l r a, taking into account l a, we take r a. Then (9) i implified a d Z 1 (r) dr +(k1 k0)z 1 (r) = 0 (13) The olution of (13) i Z 1 (r) =Be i k1 k 0 [r (a l)] + Ce i k1 k 0 [r (a l)]. (14) According to the impedance boundary condition at r = a l, weget E (1) θ η 1 H ϕ (1) = g (15) r=a l At lower frequencie, the urface impedance g can be approximated a, g k ( ) 1 k1 1 (16) k k Then, e = C B = k1 k 0 k 1 g. (17) k 1 k 0 + k 1 g At the boundary r = a, with the condition E (0) θ z=a = E (1) θ z=a and H ϕ (0) z=a = H ϕ (1) z=a,weget ( )1 AW 3 k0 = ib k1 k 0 a k 0 (ei k1 k 0 l e e i k1 k 0 l ) (18) A k0w = B k1(e i k1 k 0 l + e e i k1 k 0 l ). (19)

6 108Li, Park, and Zhang Then, the differential equation i obtained a follow where q = ik 0 k 1 k 0 k 1 W (t) qw (t) = 0 (0) ( k0 a ) e e i k1 k 0 l 1+ e e i. (1) k1 k 0 l In the pherical conducting cae g 0, then e = 1, the above formula reduce to q = k 0 k1 k ( ) 1 0 k0 a 3 tan( k1 k 0 l). () k 1 The expreion of () coincide with that given in the literature [8]. If t are the root of the differential equation (37). Then the potential function U 0 in the air i expreed a U 0 = 1 A F (z r )P ν (co(π θ)) (3) r where the height-gain function F (z r )i F (z r )= W (t y). (4) W (t ) With the imilar method in [8], the coefficient A in (3) will be determined. ( ) 1/3 A = k0idl F(z ) e iνπ ka ωε 0 t q. (5) When ν i very large and θ i not cloe to 0 and π [8]. 1 P ν (co(π θ)) = πk 0 a in θ exp [i(k 0aθ + t x)] exp [i(ν + 1 ] 4 )π (6) where x =( k 0a )1/3 θ. Finally, the accurate and complete formula for the component E r, E θ and ηh ϕ of the electromagnetic field generated by a vertical

7 EM field in three-layered pherical region 109 electric dipole at (a + z, 0, 0) over the dielectric-coated pherical earth are expreed a follow. E r(a + z r,θ,ϕ) E θ (a + z r,θ,ϕ) = iidl η η 0 H ϕ (a + z r,θ,ϕ) λa F (z ) F (z r ) t q e itx ei(kaθ+ π 4 ) πx if (z ) F(z) k 0 z z=z r [ ( ) ] /3 e itx θ in θ k 0 (t q ) 1+ t (7) ka F (z ) F (z r ) [ ] (t q ) 1+ t ( e itx ka )/3 where r = a + z r. To the aphalt- and cement-coated pherical earth or ice-coated pherical eawater at the lower frequencie, it i eaily to compute the component of the electromagnetic field veru the propagating ditance at different thicknee of the dielectric layer. If the thickne of the dielectric layer i not zero, the parameter t can be computed directly by uing the method preented in [8]. Auming that the radiu of the earth i taken a a = 6370 km, the dielectric layer a l < r < a(region 1) i characterized by the relative permittivity ε r1 = 15, and conductivity σ 1 = 10 5 S/m, the earth in the region r < a l (Region ) i characterized by the relative permittivity ε r = 100, and conductivity σ = 4 S/m, t can be calculated eaily. Fig. and 3 how both the real part and the imaginary part of t 1 and t t 5, repectively. It hould be noted that q i alway imaginary and will increae with the thickne l from zero when the thickne l of the dielectric layer atify the condition 0 < k1 k 0 l < π. From Fig., it can be een that the real part of t 1 increae ignificantly when l>40 m and the imaginary part of t 1 approache zero when the thickne of the dielectric layer l i between 40 m and 50 m. It mean that the firt propagation mode will attenuate very lowly along the urface of the pherical ground when the thickne of the dielectric layer. Therefore, the firt term may be conidered a the trapped urface wave [8]. In thi paper, we will define the trapped urface wave generated by a vertical electric dipole a the trapped wave of electric type. Taking into account the conductivity σ 1 of the dielectric layer, the imaginary part of t 1 increae ignificantly when l>50 m.

8 110 Li, Park, and Zhang Re t l (m) (a) The real part of t Re t l (m) (b) The imaginary part of t 1 Figure. The value of t 1 veru the thickne of the dielectric layer l.

9 EM field in three-layered pherical region 111 In contrat with the effect on E r and H ϕ by the thickne l, we will compute E r and H ϕ at different thicknee l. Auming a = 6370 km, f = 100 khz, ε r1 = 15, ε r = 100, σ 1 = 10 5 S/m, σ = 4 S/m, and z r = z =0, E r in V/m and H ϕ in A/m due to unit vertical electric dipole at l =0m,l = 50 m, and l = 100 m are plotted in Fig. 4 and 5, repectively. From Fig. 5, it can be een that the trapped urface wave of electric magnetic can be excited efficiently when the thickne l of the dielectric layer i between 40 m and 50 m and atifie the condition 0 < k1 k 0 l<π. Compared with the computation preented in [8], it i concluded that the electromagnetic field i affected largely by the conductivity of the earth and the electromagnetic field attenuate rapidly for the dielectric-coated earth with high conductivity. In the cae of the perfect conducting earth, the attenuation of the electromagnetic field over the urface of the dielectric-coated pherical conducting earth will be larget. 3. THE FIELD IN THE AIR GENERATED BY A MAGNETIC DIPOLE OVER THE SURFACE OF THE DIELECTRIC-COATED SPHERICAL ELECTRICALLY EARTH If the excitation ource i replaced by a vertical magnetic dipole with it moment M = Ida, and da i the area of the loop, the nonzero component of the electromagnetic field in the air are E ϕ, H r, and H θ (TE mode), which the field i defined a the field of electric type. E ϕ (i), H r (i), and H (i) θ can be expreed a follow, E ϕ (i) = iωµ 0 r θ (V ir) (8) ( ) H r (i) = r + k (V i r) (9) H (i) θ = 1 r θ r (V ir). (30) where the potential function V i i the olution of in the calar Helmholtz equation ( + ki )V i =0, i =0, 1. (31) It i noted that the urface impedance of magnetic type i E (1) ϕ η 1 H (1) θ = g (3) r=a l

10 11 Li, Park, and Zhang Re t 5 Re t - Re t Re t 4 Re t Re t l (m) (a) The real part of t t Re t 5 Re t - Re t Re t 4 Re t Re t l (m) (b) The imaginary part of t t 5 Figure 3. The value of t t 5 veru the thickne of the dielectric layer l.

11 EM field in three-layered pherical region 113 1e-7 l = 0 m l = 50 m l = 100 m The abolute value 1e-8 1e The propagating ditance (km) Figure 4. E r in V/m due to unit vertical electric dipole veru the propagating ditance ρ at l = 0 m, l = 50 m, and l = 100 m, repectively: f = 100 khz, a = 6370 km, ε r1 = 15, ε r = 100, σ 1 =10 5 S/m, σ = 4 S/m, and z r = z =0. 1e-9 l = 0 m l = 50 m l = 100 m The abolute value 1e-10 1e-11 1e The propagating ditance (km) Figure 5. H ϕ in A/m due to unit vertical electric dipole veru the propagating ditance ρ at l = 0 m, l = 50 m, and l = 100 m, repectively: f = 100 khz, a = 6370 km, ε r1 = 15, ε r = 100, σ 1 =10 5 S/m, σ = 4 S/m, and z r = z =0.

12 114 Li, Park, and Zhang Follow the ame tep in Section, we will obtain the formula of the component E ϕ, H r, and H θ in the air generated by a vertical magnetic dipole over the dielectric earth. E ϕ η 0 H r η 0 H θ where G (z )G (z r ) t (q h ) = ωµ 0Ida e i(kaθ+ π 4 ) G (z )G (z r ) πx λa θ in θ t (q h ) i G (z ) G(z) k 0 z z=z r t (q h ) q h = i k 1 k 0 k 0 ( k0 a m = ) 1 3 e it x e it x e it x. (33) 1 m e i k1 k 0 l 1+ m e i. (34) k1 k 0 l k 1 k 0 g + k 1 k 1 k 0 g k 1. (35) In the pherical perfectly conducting cae g 0, then m = 1, the above formula reduce to ( ) 1 q h k0 a 3 k1 = k 0 ]. (36) k 0 tan [k 1 k0 l t are the root of the following differential equation W (t ) q h W (t ) = 0 (37) Here, W (t ) i the Airy function of the econd kind. The height-gain function G (z r )i G (z r )= W (t y) (38) W (t ) Auming that all the parameter of the dielectric layer a l< r<a(region 1) and the earth in the region r<a l (Region ) are ame with thoe in Section. Both the real part and the imaginary part of t 1 t 5 are calculated at f = 100 khz and plotted in Fig. 6, repectively. Fig. 6 how both the real part and the imaginary part of

13 EM field in three-layered pherical region Re t 5 Re t 4 Re t 1 - Re t 5 3 Re t 3 Re t Re t l (m) (a) The real part of t 1 t Im t 5 6 Im t 4 Im t 1 - Im t Im t 3 Im t Im t l (m) (b) The imaginary part of t 1 t 5 Figure 6. The value of t 1 t 5 layer l. veru the thickne of the dielectric

14 116 Li, Park, and Zhang t 1 t 5. It hould be noted that q i alway negative and will increae from negative infinite to zero with the thickne l when the thickne l of the dielectric layer atify the condition 0 < k1 k 0 l< π. From Fig. 6, it i een that the urface wave cannot be excited by a vertical magnetic dipole when the thickne l atifie the condition 0 < k1 k 0 l<π. When the thickne l of the dielectric layer atifie the condition π k < 1 k 0 l <π, q i alway poitive and will increae from zero to infinite with the thickne l. When the thickne of the layer i larger than a certain value, the real part of t 1 increae ignificantly and the imaginary part of t 1 approache to zero. From the dicuion in Section, it i een that the urface wave can be excited by a vertical magnetic dipole when the thickne l of the dielectric layer i larger than a certain value and atifie the condition π k < 1 k 0 l<π. In thi paper, the trapped urface wave generated by a vertical electric dipole i defined a the trapped urface wave of magnetic type. In fact, the thickne l of the electric layer uch a ice or permafrot i le than 00 m, the above condition cannot be atified. In contrat with the effect on E ϕ and H r by the thickne l, we will compute E ϕ and H r at different thicknee l. Auming a = 6370 km, f = 100 khz, ε r1 = 15, ε r = 100, σ 1 = 10 5 S/m, σ = 4 S/m, and z r = z =0, E ϕ in V/m and H r in A/m due to unit vertical electric dipole at l =0m,l = 50 m, and l = 100 m are plotted in Fig. 7 and 8, repectively. 4. THE FIELD IN THE AIR GENERATED BY A HORIZONTAL ELECTRIC DIPOLE OVER THE SURFACE OF THE DIELECTRIC-COATED SPHERICAL ELECTRICALLY EARTH From the above expreion of the electromagnetic field of the vertical electric dipole and thoe of the vertical magnetic dipole in the air, uing reciprocity theorem, the accurate and complete formula of the vertical component Er he and Hr he for the electromagnetic field generated by a horizontal electric dipole upon the urface of the dielectric-coated pherical conducting earth in the air can be developed eaily. E he r = iηidhe ei(k 0aθ+ π 4 ) πx co ϕ λa θ in θ if (z r ) F(z) k 0 z z=z (t q )[1+ t ( k 0 a )/3 ] eit x (39)

15 EM field in three-layered pherical region 117 1e-10 1e-11 1e-1 l = 0 m l = 50 m l = 100 m The abolute value 1e-13 1e-14 1e-15 1e-16 1e-17 1e The propagating ditance (km) Figure 7. E ϕ in V/m due to unit vertical electric dipole veru the propagating ditance ρ at l = 0 m, l = 50 m, and l = 100 m, repectively: f = 100 khz, a = 6370 km, ε r1 = 15, ε r = 100, σ 1 =10 5 S/m, σ = 4 S/m, and z r = z =0. 1e-1 1e-13 1e-14 l = 0 m l = 50 m l = 100 m The abolute value 1e-15 1e-16 1e-17 1e-18 1e-19 1e-0 1e The propagating ditance (km) Figure 8. H r in A/m due to unit vertical electric dipole veru the propagating ditance ρ at l = 0 m, l = 50 m, and l = 100 m, repectively: f = 100 khz, a = 6370 km, ε r1 = 15, ε r = 100, σ 1 =10 5 S/m, σ = 4 S/m, and z r = z =0.

16 118Li, Park, and Zhang Hr he = iidhe G (z )G (z r ) λa t (q h ) e it x (40) In pherical coordinate, from Maxwell equation, it i eaily to expre the ret four component Eθ he, ϕ, Hθ he and Hϕ he in term of Er he and H he r ( r + k ( ) r + k ei(k 0aθ+ π 4 ) θ in θ πx in ϕ ) r + k (rhϕ he Er he ) = iωε 0 θ + 1 in θ ( ) r + k (reθ he ) = Er he r θ + iωµ 0 in θ ( ) (rhθ he ) = iωε 0 Er he in θ ϕ From (41), we get (re he ϕ ) = 1 in θ H he r ϕ H he r r ϕ + Hr he r θ E he r r ϕ iωµ 0 H he r θ (41) (4) (43) (44) ν(ν +1) (H he Er he ϕ )=iωε 0 r θ + 1 in θ H he r r ϕ (45) Taking into account the relation ν(ν +1)/r k a, the component H he ϕ i formulated a follow. Hϕ he = iidhe λa ei(k 0aθ+ π 4 ) θ in θ if (z r ) F(z) k 0 z t q 1 + k 0 a in θ πx co ϕ z=z e itx G (z) k 0 z z=zr G (z ) t (q h ) e it x (46) Similarly, the ret three component Eθ he, Hhe θ and Eϕ he are written a Hθ he = Idhe λa k 0 a in θ ei(k0aθ+ π 4 ) πx in ϕ θ in θ i F (z r ) F(z) k 0 z z=z (t q )[1 + t ( k 0 a )/3 ] eit x

17 EM field in three-layered pherical region 119 G (z) k 0 z z=zr G (z ) t (q h ) e it x (47) Eθ he = ηidhe λa i ei(k 0aθ+ π 4 ) θ in θ F (z) k 0 z 1 + k 0 a in θ Eϕ he = iηidhe λa ei(k 0aθ+ π 4 ) θ in θ i k 0 a in θ + πx co ϕ z=zr F(z) k 0 z t q G (z r ) G (z ) t (q h ) z=z πx in ϕ F (z) k 0 z G (z r ) G (z ) t (q h ) e itx ] e it x z=zr F(z) k 0 z z=z (t q )[1 + t ( k 0 a )/3 ] eit x ] e it x (48) (49) where upercript he deignate the horizontal electric dipole ource. From the above derivation, it can be een that the electromagnetic field generated by a horizontal electric dipole include the electromagnetic field of electric type and the field of magnetic type. Obviouly, for the electromagnetic field of a horizontal electric dipole, the trapped urface wave of electric type can be excited efficiently when the thickne l of the dielectric layer i omewhat and atifie the condition 0 < k1 k 0 l<π and the trapped urface wave of magnetic type can be excited when the thickne l of the dielectric layer i omewhat and atifie the condition π k < 1 k 0 l<π. 5. CONCLUSION In the above analyi, the imple explicit formula have been derived for the electromagnetic field over the urface of a dielectric-coated pherical electrically earth when a dipole ource (a vertical electric dipole, vertical magnetic dipole or horizontal electric dipole) i in the air. The computation how that both the trapped urface wave of electric type and that of magnetic type can be exited efficiently under the different condition. Alo, the electromagnetic field i affected largely by the conductivity of the earth.

18 10 Li, Park, and Zhang ACKNOWLEDGMENT Thi work wa upported by the National Reearch Laboratory (NRL) of Minitry of Science and Technology, in Korea under contract M REFERENCES 1. King, R. W. P., The electromagnetic field of a horizontal electric dipole in the preence of a three-layered region, Journal of Applied Phy., Vol. 69, No. 1, , King, R. W. P., The electromagnetic field of a horizontal electric dipole in the preence of a three-layered region: upplement, Journal of Applied Phy., Vol. 74, No. 8, , King, R. W. P. and S. S. Sandler, The electromagnetic field of a vertical electric dipole in the preence of a three-layered region, Radio Sci., Vol. 9, No. 1, , Wait, J. R., Comment on The electromagnetic field of a vertical electric dipole in the preence of a three-layered region by R. W. P. King and S. S. Sandler, Radio Sci., Vol. 33, No., 51 53, King, R. W. P. and S. S. Sandler, Reply, Radio Sci., Vol. 33, No., 55 56, Zhang, H. Q. and W. Y. Pan, Electromagnetic field of a vertical electric dipole on a perfect conductor coated with a dielectric layer, Radio Science, Vol. 37, No. 4, /000RS00348, Fock, V. A., Electromagnetic Diffraction and Propagation Problem, 191 1, Pergamon Pre, Ltd., Oxford, Pan, W. Y. and H. Q. Zhang, Electromagnetic field of a vertical electric dipole on the pherical conductor covered with a dielectric layer, Radio Sci., Vol. 38, No. 3, /00RS00689, Li, K., S. O. Park, and H. Q. Zhang, Electromagnetic field generated by a vertical magnetic dipole upon the urface of a dielectric-coated pherical conducting earth, Journal of Electromagnetic Wave and Application, Vol. 1, No. 1, , Dec Gradhteyn, I. S. and I. M. Ryzhik, Table of Integral, Serie, and Product, Academic Pre, New York, 1980.

19 EM field in three-layered pherical region 11 Kai Li received the degree of B.Sc. degree in Phyic, M.Sc. degree in Radio Phyic and Ph.D. degree Atrophyic from Fuyang Normal Univerity, Anhui, China, in 1990, Xidian Univerity, Xi an, Shaanxi, China, in 1994 and Shaanxi Atronomical Obervatory, the Chinee Academy of Science, Shaanxi, China, in 1998, repectively. From Aug to Dec. 000, he had been on the faculty of China Reearch Intitute Radiowave Propagation (CRIRP). From Jan. 001 to Dec. 00, he wa a potdoctoral fellow at Information and Communication Univerity, Taejon, South Korea. Since Jan. 003, he ha been a reearch fellow in the School of Electrical and Electronic Engineering, Nanyang Technology Univerity (NTU), Singapore. Hi current reearch interet include electromagnetic theory and radio wave propagation. Dr. Li i a enior member of Chinee Intitute of Electronic (CIE) and a member of Chinee Intitute of Space Science (CISS). Seong-Ook Park wa born in KyungPook, Korea, on December, He received the B.S. degree in electrical engineering in 1987 from KyungPook National Univerity, Korea. In 1989, he received the M.S. degree in electrical engineering from Korea Advanced Intitute of Science and Technology, Seoul. He i currently working toward the Ph.D. degree from Arizona State Univerity. From March 1989 to Augut 1993, he wa an reearch engineer with Korea Telecom, Korea, working with microwave ytem and network. Later, he joined the Telecommunication Reearch Center at Arizona State Univerity. Hi reearch interet include analytical and numerical technique in the area of microwave integrated circuit. Mr. Park i a member of Phi Kappa Phi Scholatic Honor Societie. Hong-Qi Zhang received the degree of B.Sc. degree in Radio Phyic and Ph.D. in degree Atrophyic from Lanzhou Univerity at Lanzhou, China, in 1986 and Shaanxi Atronomical Obervatory, the Chinee Academy of Science, Shaanxi, China, in 001, repectively. From 1986, he ha been on the faculty of China Reearch Intitute Radiowave Propagation (CRIRP). At 1995 he advanced to hi preent rank of enior engineer at CRIRP. Hi reearch interet include wave propagation, thundertorm location technology, and tranient electromagnetic field. Dr. Zhang won the Chinee Academy of Science award for hi excellent Ph.D. thei in 001.

ELECTROMAGNETIC FIELD IN THE AIR GENERATED BY A HORIZONTAL ELECTRIC DIPOLE LOCATED IN THE SPHERICAL ELECTRICALLY EARTH COATED WITH A DIELECTRIC LAYER

ELECTROMAGNETIC FIELD IN THE AIR GENERATED BY A HORIZONTAL ELECTRIC DIPOLE LOCATED IN THE SPHERICAL ELECTRICALLY EARTH COATED WITH A DIELECTRIC LAYER J. of Electromagn. Wave and Appl., Vol. 1, No. 1, 199{11, ELECTROMAGNETIC FIELD IN THE AIR GENERATED BY A HORIZONTAL ELECTRIC DIPOLE LOCATED IN THE SPHERICAL ELECTRICALLY EARTH COATED WITH A DIELECTRIC

More information

ELECTROMAGNETIC FIELD FROM A HORIZONTAL ELECTRIC DIPOLE IN THE SPHERICAL ELECTRICALLY EARTH COATED WITH N-LAYERED DIELECTRICS

ELECTROMAGNETIC FIELD FROM A HORIZONTAL ELECTRIC DIPOLE IN THE SPHERICAL ELECTRICALLY EARTH COATED WITH N-LAYERED DIELECTRICS Progre In Electromagnetic Reearch, PIER 54, 1 44, 005 ELECTROMAGNETIC FIELD FROM A HORIZONTAL ELECTRIC DIPOLE IN THE SPHERICAL ELECTRICALLY EARTH COATED WITH N-LAYERED DIELECTRICS K. Li The Electromagnetic

More information

ELECTROMAGNETIC FIELD FOR A HORIZONTAL ELECTRIC DIPOLE BURIED INSIDE A DIELECTRIC LAYER COATED HIGH LOSSY HALF SPACE

ELECTROMAGNETIC FIELD FOR A HORIZONTAL ELECTRIC DIPOLE BURIED INSIDE A DIELECTRIC LAYER COATED HIGH LOSSY HALF SPACE Progress In Electromagnetics Research, PIER 5, 63 86, 5 ELECTROMAGNETIC FIELD FOR A HORIZONTAL ELECTRIC DIPOLE BURIED INSIDE A DIELECTRIC LAYER COATED HIGH LOSSY HALF SPACE H.-Q. Zhang National Time Service

More information

Lecture 23 Date:

Lecture 23 Date: Lecture 3 Date: 4.4.16 Plane Wave in Free Space and Good Conductor Power and Poynting Vector Wave Propagation in Loy Dielectric Wave propagating in z-direction and having only x-component i given by: E

More information

Bogoliubov Transformation in Classical Mechanics

Bogoliubov Transformation in Classical Mechanics Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How

More information

PHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased PHYS B - HW #6 Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed. Problem. from Griffith Show that the following, A µo ɛ o A V + A ρ ɛ o Eq..4 A

More information

EELE 3332 Electromagnetic II Chapter 10

EELE 3332 Electromagnetic II Chapter 10 EELE 333 Electromagnetic II Chapter 10 Electromagnetic Wave Propagation Ilamic Univerity of Gaza Electrical Engineering Department Dr. Talal Skaik 01 1 Electromagnetic wave propagation A changing magnetic

More information

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement

More information

Tuning of High-Power Antenna Resonances by Appropriately Reactive Sources

Tuning of High-Power Antenna Resonances by Appropriately Reactive Sources Senor and Simulation Note Note 50 Augut 005 Tuning of High-Power Antenna Reonance by Appropriately Reactive Source Carl E. Baum Univerity of New Mexico Department of Electrical and Computer Engineering

More information

A novel protocol for linearization of the Poisson-Boltzmann equation

A novel protocol for linearization of the Poisson-Boltzmann equation Ann. Univ. Sofia, Fac. Chem. Pharm. 16 (14) 59-64 [arxiv 141.118] A novel protocol for linearization of the Poion-Boltzmann equation Roumen Tekov Department of Phyical Chemitry, Univerity of Sofia, 1164

More information

FI 3221 ELECTROMAGNETIC INTERACTIONS IN MATTER

FI 3221 ELECTROMAGNETIC INTERACTIONS IN MATTER 6/0/06 FI 3 ELECTROMAGNETIC INTERACTION IN MATTER Alexander A. Ikandar Phyic of Magnetim and Photonic CATTERING OF LIGHT Rayleigh cattering cattering quantitie Mie cattering Alexander A. Ikandar Electromagnetic

More information

A note on the bounds of the error of Gauss Turán-type quadratures

A note on the bounds of the error of Gauss Turán-type quadratures Journal of Computational and Applied Mathematic 2 27 276 282 www.elevier.com/locate/cam A note on the bound of the error of Gau Turán-type quadrature Gradimir V. Milovanović a, Miodrag M. Spalević b, a

More information

Effects of vector attenuation on AVO of offshore reflections

Effects of vector attenuation on AVO of offshore reflections GEOPHYSICS, VOL. 64, NO. 3 MAY-JUNE 1999); P. 815 819, 9 FIGS., 1 TABLE. Effect of vector attenuation on AVO of offhore reflection J. M. Carcione ABSTRACT Wave tranmitted at the ocean bottom have the characteritic

More information

in a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim

in a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College

More information

BACKSCATTER FROM A SPHERICAL INCLUSION WITH COMPLIANT INTERPHASE CHARACTERISTICS. M. Kitahara Tokai University Shimizu, Shizuoka 424, Japan

BACKSCATTER FROM A SPHERICAL INCLUSION WITH COMPLIANT INTERPHASE CHARACTERISTICS. M. Kitahara Tokai University Shimizu, Shizuoka 424, Japan BACKSCATTER FROM A SPHERICAL INCLUSION WITH COMPLIANT INTERPHASE CHARACTERISTICS M. Kitahara Tokai Univerity Shimizu, Shizuoka 424, Japan K. Nakagawa Total Sytem Intitute Shinjuku, Tokyo 162, Japan J.

More information

AN EFFICIENT INTEGRAL TRANSFORM TECHNIQUE OF A SINGULAR WIRE ANTENNA KERNEL. S.-O. Park

AN EFFICIENT INTEGRAL TRANSFORM TECHNIQUE OF A SINGULAR WIRE ANTENNA KERNEL. S.-O. Park AN EFFICIENT INTEGRAL TRANSFORM TECHNIQUE OF A SINGULAR WIRE ANTENNA KERNEL S.-O. Park Department of Electronics Engineering Information and Communications University 58-4 Hwaam-dong, Yusung-gu Taejon,

More information

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.

More information

MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE

MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE Proceeding 3rd Annual Conference IEEE/EMBS Oct.5-8, 1, Itanbul, TURKEY MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE S.S. Şeker, B. Yeldiren Boğaziçi Univerity,

More information

Lecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)

Lecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas) Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.

More information

Control Systems Analysis and Design by the Root-Locus Method

Control Systems Analysis and Design by the Root-Locus Method 6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If

More information

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Nonlinear Single-Particle Dynamics in High Energy Accelerators Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction

More information

0 of the same magnitude. If we don t use an OA and ignore any damping, the CTF is

0 of the same magnitude. If we don t use an OA and ignore any damping, the CTF is 1 4. Image Simulation Influence of C Spherical aberration break the ymmetry that would otherwie exit between overfocu and underfocu. One reult i that the fringe in the FT of the CTF are generally farther

More information

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,

More information

Online supplementary information

Online supplementary information Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne

More information

III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES

III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation

More information

Advanced Smart Cruise Control with Safety Distance Considered Road Friction Coefficient

Advanced Smart Cruise Control with Safety Distance Considered Road Friction Coefficient International Journal of Computer Theory and Engineering, Vol. 8, No. 3, June 06 Advanced Smart Cruie Control with Safety Ditance Conidered Road Friction Coefficient Doui Hong, Chanho Park, Yongho Yoo,

More information

Analysis of Step Response, Impulse and Ramp Response in the Continuous Stirred Tank Reactor System

Analysis of Step Response, Impulse and Ramp Response in the Continuous Stirred Tank Reactor System ISSN: 454-50 Volume 0 - Iue 05 May 07 PP. 7-78 Analyi of Step Repone, Impule and Ramp Repone in the ontinuou Stirred Tank Reactor Sytem * Zohreh Khohraftar, Pirouz Derakhhi, (Department of hemitry, Science

More information

Physics 2212 G Quiz #2 Solutions Spring 2018

Physics 2212 G Quiz #2 Solutions Spring 2018 Phyic 2212 G Quiz #2 Solution Spring 2018 I. (16 point) A hollow inulating phere ha uniform volume charge denity ρ, inner radiu R, and outer radiu 3R. Find the magnitude of the electric field at a ditance

More information

Representation Formulas of Curves in a Two- and Three-Dimensional Lightlike Cone

Representation Formulas of Curves in a Two- and Three-Dimensional Lightlike Cone Reult. Math. 59 (011), 437 451 c 011 Springer Bael AG 14-6383/11/030437-15 publihed online April, 011 DOI 10.1007/0005-011-0108-y Reult in Mathematic Repreentation Formula of Curve in a Two- and Three-Dimenional

More information

Research Article Reliability of Foundation Pile Based on Settlement and a Parameter Sensitivity Analysis

Research Article Reliability of Foundation Pile Based on Settlement and a Parameter Sensitivity Analysis Mathematical Problem in Engineering Volume 2016, Article ID 1659549, 7 page http://dxdoiorg/101155/2016/1659549 Reearch Article Reliability of Foundation Pile Baed on Settlement and a Parameter Senitivity

More information

Reflection and transmission of obliquely incident graphene plasmons by discontinuities in

Reflection and transmission of obliquely incident graphene plasmons by discontinuities in Home Search Collection Journal About Contact u My IOPcience Reflection and tranmiion of obliquely incident graphene plamon by dicontinuitie in urface conductivity: obervation of the Brewter-like effect

More information

Introduction to Laplace Transform Techniques in Circuit Analysis

Introduction to Laplace Transform Techniques in Circuit Analysis Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found

More information

Computers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order

Computers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order Computer and Mathematic with Application 64 (2012) 2262 2274 Content lit available at SciVere ScienceDirect Computer and Mathematic with Application journal homepage: wwweleviercom/locate/camwa Sharp algebraic

More information

Progress In Electromagnetics Research M, Vol. 18, 31 42, 2011

Progress In Electromagnetics Research M, Vol. 18, 31 42, 2011 Progre In Electromagnetic Reearch M, Vol. 18, 31 42, 2011 UNIFORM PLANE WAVE REFLECTION FROM PEC PLANE EMBEDDED IN A NONLINEAR MEDIUM M. M. Hunain and M. J. Mughal Faculty of Electronic Engineering GIK

More information

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization Finite Element Analyi of a Fiber Bragg Grating Accelerometer for Performance Optimization N. Baumallick*, P. Biwa, K. Dagupta and S. Bandyopadhyay Fiber Optic Laboratory, Central Gla and Ceramic Reearch

More information

Research Article Existence for Nonoscillatory Solutions of Higher-Order Nonlinear Differential Equations

Research Article Existence for Nonoscillatory Solutions of Higher-Order Nonlinear Differential Equations International Scholarly Reearch Network ISRN Mathematical Analyi Volume 20, Article ID 85203, 9 page doi:0.502/20/85203 Reearch Article Exitence for Nonocillatory Solution of Higher-Order Nonlinear Differential

More information

Automatic Control Systems. Part III: Root Locus Technique

Automatic Control Systems. Part III: Root Locus Technique www.pdhcenter.com PDH Coure E40 www.pdhonline.org Automatic Control Sytem Part III: Root Locu Technique By Shih-Min Hu, Ph.D., P.E. Page of 30 www.pdhcenter.com PDH Coure E40 www.pdhonline.org VI. Root

More information

Coupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries

Coupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries Coupling of Three-Phae Sequence Circuit Due to Line and Load Aymmetrie DEGO BELLAN Department of Electronic nformation and Bioengineering Politecnico di Milano Piazza Leonardo da inci 01 Milano TALY diego.ellan@polimi.it

More information

Dyadic Green s Function

Dyadic Green s Function EECS 730 Winter 2009 c K. Sarabandi Dyadic Green Function A mentioned earlier the application of dyadic analyi facilitate imple manipulation of field ector calculation. The ource of electromagnetic field

More information

Riemann s Functional Equation is Not a Valid Function and Its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr.

Riemann s Functional Equation is Not a Valid Function and Its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr. Riemann Functional Equation i Not a Valid Function and It Implication on the Riemann Hypothei By Armando M. Evangelita Jr. armando78973@gmail.com On Augut 28, 28 ABSTRACT Riemann functional equation wa

More information

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281 72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition

More information

Hylleraas wavefunction for He. dv 2. ,! r 2. )dv 1. in the trial function. A simple trial function that does include r 12. is ) f (r 12.

Hylleraas wavefunction for He. dv 2. ,! r 2. )dv 1. in the trial function. A simple trial function that does include r 12. is ) f (r 12. Hylleraa wavefunction for He The reaon why the Hartree method cannot reproduce the exact olution i due to the inability of the Hartree wave-function to account for electron correlation. We know that the

More information

BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL

BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)

More information

Lecture 10 Filtering: Applied Concepts

Lecture 10 Filtering: Applied Concepts Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering

More information

Lecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs)

Lecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs) Lecture 4 Topic 3: General linear model (GLM), the fundamental of the analyi of variance (ANOVA), and completely randomized deign (CRD) The general linear model One population: An obervation i explained

More information

Unbounded solutions of second order discrete BVPs on infinite intervals

Unbounded solutions of second order discrete BVPs on infinite intervals Available online at www.tjna.com J. Nonlinear Sci. Appl. 9 206), 357 369 Reearch Article Unbounded olution of econd order dicrete BVP on infinite interval Hairong Lian a,, Jingwu Li a, Ravi P Agarwal b

More information

Dimension Effect on Dynamic Stress Equilibrium in SHPB Tests

Dimension Effect on Dynamic Stress Equilibrium in SHPB Tests International Journal of Material Phyic. ISSN 97-39X Volume 5, Numer 1 (1), pp. 15- International Reearch Pulication Houe http://www.irphoue.com Dimenion Effect on Dynamic Stre Equilirium in SHPB Tet Department

More information

FUNDAMENTALS OF POWER SYSTEMS

FUNDAMENTALS OF POWER SYSTEMS 1 FUNDAMENTALS OF POWER SYSTEMS 1 Chapter FUNDAMENTALS OF POWER SYSTEMS INTRODUCTION The three baic element of electrical engineering are reitor, inductor and capacitor. The reitor conume ohmic or diipative

More information

ELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER

ELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER Progress In Electromagnetics Research Letters, Vol. 31, 55 64, 2012 ELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER Y.-J. Zhi

More information

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang

More information

How a charge conserving alternative to Maxwell s displacement current entails a Darwin-like approximation to the solutions of Maxwell s equations

How a charge conserving alternative to Maxwell s displacement current entails a Darwin-like approximation to the solutions of Maxwell s equations How a charge conerving alternative to Maxwell diplacement current entail a Darwin-like approximation to the olution of Maxwell equation 12 ab Alan M Wolky 1 Argonne National Laboratory 9700 South Ca Ave

More information

Riemann s Functional Equation is Not Valid and its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr.

Riemann s Functional Equation is Not Valid and its Implication on the Riemann Hypothesis. Armando M. Evangelista Jr. Riemann Functional Equation i Not Valid and it Implication on the Riemann Hypothei By Armando M. Evangelita Jr. On November 4, 28 ABSTRACT Riemann functional equation wa formulated by Riemann that uppoedly

More information

Math 273 Solutions to Review Problems for Exam 1

Math 273 Solutions to Review Problems for Exam 1 Math 7 Solution to Review Problem for Exam True or Fale? Circle ONE anwer for each Hint: For effective tudy, explain why if true and give a counterexample if fale (a) T or F : If a b and b c, then a c

More information

Manprit Kaur and Arun Kumar

Manprit Kaur and Arun Kumar CUBIC X-SPLINE INTERPOLATORY FUNCTIONS Manprit Kaur and Arun Kumar manpreet2410@gmail.com, arun04@rediffmail.com Department of Mathematic and Computer Science, R. D. Univerity, Jabalpur, INDIA. Abtract:

More information

EE Control Systems LECTURE 14

EE Control Systems LECTURE 14 Updated: Tueday, March 3, 999 EE 434 - Control Sytem LECTURE 4 Copyright FL Lewi 999 All right reerved ROOT LOCUS DESIGN TECHNIQUE Suppoe the cloed-loop tranfer function depend on a deign parameter k We

More information

DEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM

DEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM DEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM Brad M. Brown Matt N. Straer R. Paneer Selvam Univerity of Arkana Department of Civil Engineering 4190 Bell Engineering Center Fayetteville,

More information

Standard Guide for Conducting Ruggedness Tests 1

Standard Guide for Conducting Ruggedness Tests 1 Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM

More information

One Class of Splitting Iterative Schemes

One Class of Splitting Iterative Schemes One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi

More information

CHARGING OF DUST IN A NEGATIVE ION PLASMA (APS DPP-06 Poster LP )

CHARGING OF DUST IN A NEGATIVE ION PLASMA (APS DPP-06 Poster LP ) 1 CHARGING OF DUST IN A NEGATIVE ION PLASMA (APS DPP-06 Poter LP1.00128) Robert Merlino, Ro Fiher, Su Hyun Kim And Nathan Quarderer Department of Phyic and Atronomy, The Univerity of Iowa Abtract. We invetigate

More information

Modeling of seed magnetic island formation

Modeling of seed magnetic island formation EUROFUSION WPMST1-PR(16) 15656 IG Miron et al. Modeling of eed magnetic iland formation Preprint of Paper to be ubmitted for publication in 43rd European Phyical Society Conference on Plama Phyic (EPS)

More information

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014 Phyic 7 Graduate Quantum Mechanic Solution to inal Eam all 0 Each quetion i worth 5 point with point for each part marked eparately Some poibly ueful formula appear at the end of the tet In four dimenion

More information

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,

More information

Stratified Analysis of Probabilities of Causation

Stratified Analysis of Probabilities of Causation Stratified Analyi of Probabilitie of Cauation Manabu Kuroki Sytem Innovation Dept. Oaka Univerity Toyonaka, Oaka, Japan mkuroki@igmath.e.oaka-u.ac.jp Zhihong Cai Biotatitic Dept. Kyoto Univerity Sakyo-ku,

More information

Chapter 4. The Laplace Transform Method

Chapter 4. The Laplace Transform Method Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination

More information

Euler-Bernoulli Beams

Euler-Bernoulli Beams Euler-Bernoulli Beam The Euler-Bernoulli beam theory wa etablihed around 750 with contribution from Leonard Euler and Daniel Bernoulli. Bernoulli provided an expreion for the train energy in beam bending,

More information

Lecture 15 - Current. A Puzzle... Advanced Section: Image Charge for Spheres. Image Charge for a Grounded Spherical Shell

Lecture 15 - Current. A Puzzle... Advanced Section: Image Charge for Spheres. Image Charge for a Grounded Spherical Shell Lecture 15 - Current Puzzle... Suppoe an infinite grounded conducting plane lie at z = 0. charge q i located at a height h above the conducting plane. Show in three different way that the potential below

More information

Annex-A: RTTOV9 Cloud validation

Annex-A: RTTOV9 Cloud validation RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh

More information

SOME RESULTS ON INFINITE POWER TOWERS

SOME RESULTS ON INFINITE POWER TOWERS NNTDM 16 2010) 3, 18-24 SOME RESULTS ON INFINITE POWER TOWERS Mladen Vailev - Miana 5, V. Hugo Str., Sofia 1124, Bulgaria E-mail:miana@abv.bg Abtract To my friend Kratyu Gumnerov In the paper the infinite

More information

ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS. Volker Ziegler Technische Universität Graz, Austria

ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS. Volker Ziegler Technische Universität Graz, Austria GLASNIK MATEMATIČKI Vol. 1(61)(006), 9 30 ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS Volker Ziegler Techniche Univerität Graz, Autria Abtract. We conider the parameterized Thue

More information

Interaction of Pile-Soil-Pile in Battered Pile Groups under Statically Lateral Load

Interaction of Pile-Soil-Pile in Battered Pile Groups under Statically Lateral Load Interaction of Pile-Soil-Pile in Battered Pile Group under Statically Lateral Load H. Ghaemadeh 1*, M. Alibeikloo 2 1- Aitant Profeor, K. N. Tooi Univerity of Technology 2- M.Sc. Student, K. N. Tooi Univerity

More information

A Buckling Problem for Graphene Sheets. J. Galagher 1, Y. Milman 2, S. Ryan 3, D. Golovaty 3, P. Wilber 3, and A. Buldum 4

A Buckling Problem for Graphene Sheets. J. Galagher 1, Y. Milman 2, S. Ryan 3, D. Golovaty 3, P. Wilber 3, and A. Buldum 4 A Buckling Problem for Graphene Sheet J. Galagher 1, Y. Milman 2, S. Ryan 3, D. Golovaty 3, P. Wilber 3, and A. Buldum 4 1 Department of Phyic, Rocheter Intitute of Technology, Rocheter, NY 14623, USA

More information

THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER

THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Proceeding of IMAC XXXI Conference & Expoition on Structural Dynamic February -4 Garden Grove CA USA THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Yung-Sheng Hu Neil S Ferguon

More information

Estimating floor acceleration in nonlinear multi-story moment-resisting frames

Estimating floor acceleration in nonlinear multi-story moment-resisting frames Etimating floor acceleration in nonlinear multi-tory moment-reiting frame R. Karami Mohammadi Aitant Profeor, Civil Engineering Department, K.N.Tooi Univerity M. Mohammadi M.Sc. Student, Civil Engineering

More information

OMNÈS REPRESENTATIONS WITH INELASTIC EFFECTS FOR HADRONIC FORM FACTORS

OMNÈS REPRESENTATIONS WITH INELASTIC EFFECTS FOR HADRONIC FORM FACTORS FIELD THEORY, GRAVITATION AND PARTICLE PHYSICS OMNÈS REPRESENTATIONS WITH INELASTIC EFFECTS FOR HADRONIC FORM FACTORS IRINEL CAPRINI National Intitute of Phyic and Nuclear Engineering POB MG 6, Bucharet,

More information

Supplementary information. Dendritic optical antennas: scattering properties and fluorescence enhancement

Supplementary information. Dendritic optical antennas: scattering properties and fluorescence enhancement Supplementary information Dendritic optical antenna: cattering propertie and fluorecence enhancement Ke Guo 1, Aleandro Antoncecchi 1, Xuezhi Zheng 2, Mai Sallam 2,3, Ezzeldin A. Soliman 3, Guy A. E. andenboch

More information

696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of

696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q  ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of Vol 12 No 7, July 2003 cfl 2003 Chin. Phy. Soc. 1009-1963/2003/12(07)/0695-05 Chinee Phyic and IOP Publihing Ltd Lie ymmetrie and conerved quantitie of controllable nonholonomic dynamical ytem Fu Jing-Li(ΛΠ±)

More information

1. Basic introduction to electromagnetic field. wave properties and particulate properties.

1. Basic introduction to electromagnetic field. wave properties and particulate properties. Lecture Baic Radiometric Quantitie. The Beer-Bouguer-Lambert law. Concept of extinction cattering plu aborption and emiion. Schwarzchild equation. Objective:. Baic introduction to electromagnetic field:

More information

Stability. ME 344/144L Prof. R.G. Longoria Dynamic Systems and Controls/Lab. Department of Mechanical Engineering The University of Texas at Austin

Stability. ME 344/144L Prof. R.G. Longoria Dynamic Systems and Controls/Lab. Department of Mechanical Engineering The University of Texas at Austin Stability The tability of a ytem refer to it ability or tendency to eek a condition of tatic equilibrium after it ha been diturbed. If given a mall perturbation from the equilibrium, it i table if it return.

More information

Pulsed Magnet Crimping

Pulsed Magnet Crimping Puled Magnet Crimping Fred Niell 4/5/00 1 Magnetic Crimping Magnetoforming i a metal fabrication technique that ha been in ue for everal decade. A large capacitor bank i ued to tore energy that i ued to

More information

Math Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK

Math Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI

More information

On the Isentropic Forchheimer s Sound Waves Propagation in a Cylindrical Tube Filled with a Porous Media

On the Isentropic Forchheimer s Sound Waves Propagation in a Cylindrical Tube Filled with a Porous Media 5th WSEAS Int. Conf. on FLUID MECHANICS (FLUIDS') Acapulco, Mexico, January 5-7, On the Ientropic Forchheimer Sound Wave Propagation in a Cylindrical Tube Filled with a Porou Media H. M. Dwairi Civil Engineering

More information

A BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH

A BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Mathematical and Computational Application Vol. 11 No. pp. 181-191 006. Aociation for Scientific Reearch A BATCH-ARRIVA QEE WITH MTIPE SERVERS AND FZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Jau-Chuan

More information

Simulation and Analysis of Linear Permanent Magnet Vernier Motors for Direct Drive Systems

Simulation and Analysis of Linear Permanent Magnet Vernier Motors for Direct Drive Systems Available online at www.ijpe-online.com vol. 3, no. 8, December 07, pp. 304-3 DOI: 0.3940/ijpe.7.08.p.3043 Simulation and Analyi of Linear Permanent Magnet Vernier Motor for Direct Drive Sytem Mingjie

More information

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou

More information

Theoretical study of the dual harmonic system and its application on the CSNS/RCS

Theoretical study of the dual harmonic system and its application on the CSNS/RCS Theoretical tudy of the dual harmonic ytem and it application on the CSNS/RCS Yao-Shuo Yuan, Na Wang, Shou-Yan Xu, Yue Yuan, and Sheng Wang Dongguan branch, Intitute of High Energy Phyic, CAS, Guangdong

More information

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr 0.1 Related Rate In many phyical ituation we have a relationhip between multiple quantitie, and we know the rate at which one of the quantitie i changing. Oftentime we can ue thi relationhip a a convenient

More information

New bounds for Morse clusters

New bounds for Morse clusters New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität

More information

Notes on Phase Space Fall 2007, Physics 233B, Hitoshi Murayama

Notes on Phase Space Fall 2007, Physics 233B, Hitoshi Murayama Note on Phae Space Fall 007, Phyic 33B, Hitohi Murayama Two-Body Phae Space The two-body phae i the bai of computing higher body phae pace. We compute it in the ret frame of the two-body ytem, P p + p

More information

Robustness analysis for the boundary control of the string equation

Robustness analysis for the boundary control of the string equation Routne analyi for the oundary control of the tring equation Martin GUGAT Mario SIGALOTTI and Mariu TUCSNAK I INTRODUCTION AND MAIN RESULTS In thi paper we conider the infinite dimenional ytem determined

More information

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning

More information

Simulation Study on the Shock Properties of the Double-Degree-of-Freedom Cushioning Packaging System

Simulation Study on the Shock Properties of the Double-Degree-of-Freedom Cushioning Packaging System Proceeding of the 7th IAPRI World Conference on Packaging Simulation Study on the Shock Propertie of the Double-Degree-of-Freedom Cuhioning Packaging Sytem Xia Zhu, Qiaoqiao Yan, Xiaoling Yao, Junbin Chen,

More information

An Analytical Solution of the Radiative Transfer Equation for Inhomogeneous Finite Medium with Fresnel Boundary Conditions

An Analytical Solution of the Radiative Transfer Equation for Inhomogeneous Finite Medium with Fresnel Boundary Conditions rab ournal of Nuclear Science pplication, 46(3), (4-5) 3 n nalytical Solution of the Radiative Tranfer Equation for Inhomogeneou Finite Medium with Frenel Boundary Condition. Elghazaly Reactor & Neutron

More information

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004 18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem

More information

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -

More information

c = ω k = 1 v = ω k = 1 ²µ

c = ω k = 1 v = ω k = 1 ²µ 3. Electromagtic Wave 3.5. Plama wave A plama i an ionized ga coniting of charged particle (e.g., electron and ion). Variou wave can be excited eaily in a plama. Wave phenomena have been an important ubject

More information

NULL HELIX AND k-type NULL SLANT HELICES IN E 4 1

NULL HELIX AND k-type NULL SLANT HELICES IN E 4 1 REVISTA DE LA UNIÓN MATEMÁTICA ARGENTINA Vol. 57, No. 1, 2016, Page 71 83 Publihed online: March 3, 2016 NULL HELIX AND k-type NULL SLANT HELICES IN E 4 1 JINHUA QIAN AND YOUNG HO KIM Abtract. We tudy

More information

The type 3 nonuniform FFT and its applications

The type 3 nonuniform FFT and its applications Journal of Computational Phyic 206 (2005) 1 5 Short Note The type 3 nonuniform FFT and it application June-Yub Lee a, *, Lelie Greengard b a Department of Mathematic, Ewha Woman Univerity, 11-1 Daehyundong,

More information

Study of a Freely Falling Ellipse with a Variety of Aspect Ratios and Initial Angles

Study of a Freely Falling Ellipse with a Variety of Aspect Ratios and Initial Angles Study of a Freely Falling Ellipe with a Variety of Apect Ratio and Initial Angle Dedy Zulhidayat Noor*, Ming-Jyh Chern*, Tzyy-Leng Horng** *Department of Mechanical Engineering, National Taiwan Univerity

More information

ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS

ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS VOLKER ZIEGLER Abtract We conider the parameterized Thue equation X X 3 Y (ab + (a + bx Y abxy 3 + a b Y = ±1, where a, b 1 Z uch that

More information

Factor Sensitivity Analysis with Neural Network Simulation based on Perturbation System

Factor Sensitivity Analysis with Neural Network Simulation based on Perturbation System 1402 JOURNAL OF COMPUTERS, VOL. 6, NO. 7, JULY 2011 Factor Senitivity Analyi with Neural Network Simulation baed on Perturbation Sytem Runbo Bai College of Water-Conervancy and Civil Engineering, Shandong

More information