Trefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies.
|
|
- Barry Edmund Turner
- 5 years ago
- Views:
Transcription
1 Trefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies. S. Yu. Reutsiy Magnetohydrodynamic Laboratory, P. O. Box 136, Mosovsi av.,199, 61037, Kharov, Uraine. e mail: reutsiy@synet.harov.com Abstract A new numerical method for scattering from inhomogeneous bodies is presented. The cases of E and H polarizated incident wave scattered by an infinite 2D cylinder are considered. The scattered field is looed for in two different domains. The first one is a bounded region inside the scattering body with an inhomogeneous permittivity ε(x, y). The second one is an unbounded homogeneous region outside the scatterer. An approximate solution for the scattered field inside the scatterer is looed for by applying the QTSM technique. The method of discrete sources is used to approximate the scattered field in the unbounded region outside the scattering body. A comparison of the numerical and analytic solutions is performed. 1 Introduction Recently various numerical methods have been developed to compute electromagnetic scattering by bodies of irregular form. Boundary methods: the extended boundary condition method (EBCM) [1], [2] or T matrix method [3], the methods based on the generalized multipole technique (GMT) [4], [5] and the method of discrete sources (MDS) [6], [7] are the fastest and most powerful tools in this field. In this paper a new numerical method of the Trefftz type is presented for this goal. We consider the 2D problem of scattering of a plane, electromagnetic wave by a penetrable continuously inhomogeneous dielectric cylinder. The axis of the cylinder is along OZ. We denote the section of the cylinder in the XY plane as Ω. The region surrounding the cylinder is free space. It is assumed throughout that all fields are taen to be harmonic and that the time dependence, taen to be e iωt, has been factored out. Suppose also that the magnetic permeability constant µ is a fixed constant everywhere: µ = µ 0. Under these assumptions Maxwell s system is: roth = i 0 εe, rote = i 0 H. (1) 1
2 It is written in dimensionless form. Here E and H are vectors of electric and magnetic fields correspondingly, ε is the relative permittivity, the wave number 0 = 2πl s /λ, λ is the wave length of an incident field and l s is the scaling length ( e.g. a typical size of a scatterer). The relative permittivity is a smooth enough complex valued function inside Ω. And it is taen to be 1 outside the scatterer (free space). The both cases of polarization of the incident wave E polarization: E (inc) = {0, 0, E (inc) z }, H (inc) = {H (inc) x, H (inc) y, 0} and H polarization: E (inc) = {E (inc) x, E y (inc), 0}, H (inc) = {0, 0, H z (inc) } are considered. Under the assumptions listed before the scattered field {E (s), H (s) } has the same polarization as the incident one. Let us use the following notations: u (i) E z (s) (x, y), E polarization (x, y) = (x, y) Ω H z (s) (x, y), H polarization u (e) (x, y) = E z (s) (x, y), E polarization H z (s) (x, y), H polarization (x, y) R 2 \ Ω Using Maxwell s system (1) one can get the systems of the scalar equations: in the E polarization case and [ ] u (e) (x, y) = 0, (x, y) R 2 \ Ω (2) [ ] + 0ε 2 (i) (x, y) u (i) (x, y) = 0, (x, y) Ω (3) [ ( 1 x ε (i) (x, y) [ ] u (e) (x, y) = 0, (x, y) R 2 \ Ω (4) x ) + ( y 1 ε (i) (x, y) y in the case of the H polarized incident wave. The boundary conditions are: E polarization : u (i) u (e) = u (inc), H polarization : u (i) u (e) = u (inc), ) ] u (i) (x, y) = 0, (x, y) Ω (5) u (i) u(e) 1 u (i) ε (i) u(e) = u(inc) = u(inc) on Ω. (6) on Ω. (7) Here denotes the derivative in the direction of the normal vector n = (n x, n y ) of Ω and we denote u (inc) =E z (inc), or H z (inc). We also suppose that the exterior field u (e) is an outgoing cylindrical wave at a large distance from the scatterer: u (e) 1 r exp( i 0 r), r, r = x 2 + y 2. (8) 2
3 2 Numerical algorithm The method presented belongs to the group of boundary methods when an approximate solution is looed for in the form of a linear combination K u (x, y) = q Ψ (x, y), where the trial functions Ψ (x, y) satisfy exactly the corresponding partial differential equation (PDE) but do not necessary satisfy the boundary conditions which are imposed on the solution. As far as the homogeneous exterior region R 2 \Ω is considered one can use a version of the method of discrete sources (MDS) a very effective technique for scattering problems in homogeneous medium developed in the last few years and which is gaining more and more interest because of its many practical applications. The description of the recent development of MDS and other references can be found in [7]. An approximate solution is looed for in the form of a linear combination of Green s functions which satisfy the radiation condition at infinity: K e u (e) (x, y) = q (e) ) K e H(2) 0 ( 0 (x x ) 2 + (y y ) 2 = q (e) Ψ(e) (x, y). Here H (2) 0 denotes the Hanel function of the second ind and zero order. The source points (x, y ) are placed inside the scatterer Ω and q (e) To write the similar approximate solution are the free parameters. K i u (i) (x, y) = q (i) Ψ(i) for inhomogeneous region Ω a new numerical technique developed (x, y) (9) in [8], [9], [10] is used. First, it assumes applying the embedding procedure to the irregular solution domain. In this wor it is assumed that the solution domain Ω may be embedded in the square Ω 0 = [ 1, 1] [ 1, 1]. Of course, if this is not the case originally, then appropriate translation and scaling operations may be performed to mae it so. Let us consider each polarization separately. 2.1 E polarization According to the embedding idea the initial differential operator L (i) = ε(i) (x, y) is replaced with L (i) = ε(i) (x, y). Here ε (i) (x, y) is a truncated series ε (i) (x, y) = +N n,m= N G nm e iπ(nx+my), (x, y) Ω 0, (10) which approximates ε (i) (x, y) for (x, y) Ω. A numerical technique called the C expansion procedure proposed by Smelov is used to get this approximation. With more details this technique is described in [8], [9], [10]. 3
4 The trial functions Ψ (i) (x, y) are taen as solutions of the PDE: ] L (i) Ψ (i) [ = + 0 ε 2 (i) (x, y) Ψ (i) = I M,l(x ζ ), (11) where I M,l (x ζ) is a δ-shaped source function which essentially differs from zero only inside some neighborhood of the source point ζ =(ξ, η), i.e. this is an approximate Dirac s δ function. However, this is an infinitely differentiable function in the form of truncated series over the trigonometric system exp[iπ(nx + my)]: I M,l (x ζ) = M n,m= M C n,m (ζ) exp[iπ(nx + my)], (12) where the regularization coefficients are: C n,m (ζ) = 1 4 r n(m, l) r m (M, l) exp[ iπ(nξ + mη)], (13) r 0 (M, l) = 1, r n (M, l) = σ l n(m), σ n (M) = sin nπ (M + 1) / nπ (M + 1). Here σ n (M) are so called the Lanczos sigma factors. Note that the regularization coefficients r n (M, l) coincide with σ l n(m) only in the case of the trigonometric expansion. The general method of the regularization procedure is presented in [11]. An approximate solution of (11) is looed for in the same form of a truncated series: Ψ (i) M (x, y) = n,m= M U () n,me iπ(nx+my). (14) Integrating over Ω 0 with the weight functions e iπ(nx+my), n, m = 0, ±1,..., ±M one gets: π 2 (n 2 + m 2 )U () nm i n, j m N These equations can be rewritten in the form: with Introducing the (2M + 1) 2 vectors: M i,j=1 G n i,m j U () i,j A i,j n,mu () i,j = C n,m () A i,j n,m = π 2 (n 2 + m 2 )δ i,n δ j,m + 2 0G n i,m j. = C n,m (ζ )=C () n,m. (15) U = {U () () () () M, M,..., U M,+M,..., U +M, M,..., U +M,+M }, C = {C () M, M,..., C() M,+M,..., C() +M, M,..., C() +M,+M } 4
5 they can be presented as a set of (2M + 1) 2 (2M + 1) 2 linear systems with the same matrix and different vectors of the right hand side: ÂU = C, = 1,..., K i. (16) It is important to note that the coefficients A i,j n,m of the matrix are obtained in an analytic way without any numerical integration. This is due to the extension of the initial irregular domain Ω to the square Ω 0 and replacing the initial complex permittivity ε (i) (x, y) with the truncated series ε (i) (x, y). Solving these systems one gets the trial functions for expansion (9). Note that using these trial functions means that in fact the Helmholtz equation is solved with the nonzero right hand side: K i [ + 0 ε 2 (i) (x, y)]u (i) = q (i) I M,l(ζ ). (17) When the source points ζ are removed from Ω inside Ω 0, then the right hand side of (17) is a small value for (x, y) Ω. So, we introduce an additional error in an approximate solution when these basis functions are used. But as the numerical experiments have shown, this error can be of the same level as the one due to the approximate satisfaction of the boundary condition. 2.2 H polarization Here we begin with an approximation of the function ν(x, y) = 1/ε (i) (x, y). It is assumed that it can be approximated by the expansion similar to (10): ν(x, y) = +N n,m= N G nm e iπ(nx+my), (x, y) Ω 0. (18) The trial functions Ψ (i) (x, y) are taen now as solutions of the PDE (cf. (11)): [ ( ν(x, y) ) + ( ν(x, y) ) ] Ψ (i) x x y y = I M,l(x ζ ). (19) An approximate solution of (19) is looed for in the same form as the one in the E polarization case (see (14)). By integrating over Ω 0 with the weight functions e iπ(nx+my), n, m = 0, ±1,..., ±M one gets: π 2 i n, j m N (in + jm)g n i,m j U () i,j + 2 0U () n,m = C () n,m. (20) It can be rewritten in the form of a linear system lie (16) with respect to the unnowns U () i,j. The coefficients of the matrix now are: A i,j n,m = π 2 (in + jm)g n i,m j + 2 0δ i,n δ j,m. 5
6 Note that, as above, when the approximation (18) is nown, then the coefficients can calculated in an analytic way without any numerical integration. Having the two linear combinations K e u (e) (x, y) = q (e) Ψ(e) (x, y), u(i) (x, y) = K i q (i) Ψ(i) (x, y) (21) we get the free parameters q (e), q(i) as a solution of an over determined system of algebraic equations. The system is obtained from the boundary conditions (6) or (7) using the collocation procedure at the collocation points (x j, y j ), j = 1,..., N c distributed uniformly on Ω. The resulting linear system with 2 N c equations and K e + K i unnowns q (e), q(i) procedure. is solved by the least squares 3 Numerical examples To test the algorithm described above the axially symmetric problems are considered. Let Ω be a dis of radius a centered at the point (0, 0) in the XY plane and let (r, ψ) be the polar coordinates: x = r cos ψ, y = r sin ψ. (22) All the calculations were performed with a = 0.5. Following [12] the source points are placed as far as possible from the boundary Ω. So, they are uniformly distributed on the circles of the radiuses r (i) = 0.95 and r (e) = 0.1 for approximation of the inner and outer solution correspondingly. Without loss of generality, the incident field is taen in the form of a plane wave propagating along the X axis: u inc (x) = exp(i 0 x) = exp(i 0 r cos ψ) = J 0 ( 0 r) + 2 i m J m ( 0 r) cos mψ. (23) Here J m (... ) denotes the Bessel function of the first ind and m th order. 1) Constant permittivity. The relative permittivity inside Ω is taen to be constant ε (i) = const. In this case both systems (2), (3) and (4), (5), written in polar coordinates, are: 2 u (e) r r u (e) r + 1 r 2 2 u (e) ψ u (e) = 0, r > a, (24) 2 u (i) r r u (i) r The boundary conditions, e.g. for H polarization, are: + 1 r 2 2 u (i) ψ ε (i) u (i) = 0, r < a. (25) 6
7 u (i) (a, ψ) u (e) (a, ψ) = J 0 ( 0 a) + 2 i m J m ( 0 a) cos mψ, (26) 1 u (i) u(e) ε (i) (a, ψ) r r (a, ψ) = 0J 0( 0 a) i m J m( 0 a) cos mψ (27) for all 0 ψ 2π. Here we denote J m(z) = dj m (z)/dz. This problem has an analytic solution in the form of expansions over the orthogonal system of the functions cos mψ: where H (2) m u (e) (r, ψ) = u (i) (r, ψ) = m=0 m=0 A m H (2) m ( 0 r) cos mψ, r > a, (28) ) B m J m ( ε (i) 0 r cos mψ, r < a, (29) denotes the Hanel function of the second ind and m th order. Using orthogonality of the cos mψ basis one gets a 2 2 linear system for each pair A m, B m. We denote the resulting solution as u (e) anlt and u(i) anlt and use it to compare with the numerical solution u(e) num and u (i) num. Note that in this case the trial functions can be found in the analytical way. Indeed, as it follows from (18) and so, eq. (20) can be easily resolved Suppose 1/ε (i), n = m = 0; G n,m = 0, otherwise. U () n,m = ε (i) C () n,m ε (i) 2 0 π2 (n 2 + m 2 ). ε (i) 2 0 π 2 ( n 2 + m 2) 0, n, m = 1, 2,... To estimate the accuracy of the calculations we use the maximal relative error. For the internal solution it is defined as: e in = max n=1,...,n t { (i) u num(x (i) n, y n (i) ) u (i) anlt (x(i) n, y n (i) ) u (i) num(x (i) n, y n (i) ) + u (i) anlt (x(i) n, y n (i) ) The similar value e out is used to compare u (e) num and u (e) anlt. We use N t }. (30) = 100 test points x (i) n, y n (i) which are distributed uniformly inside Ω to compute e in. To compute e out we choose N t = 100 test points uniformly located on the circles of radius r m = a(m + 1), m = 1,..., 10 outside Ω. These errors are placed in Table 1 for different numbers of harmonics M and degrees of freedom (DOF) 7
8 K i = K e (we tae them equal). We also compute the errors e bcd and e bcn in approximation of the boundary conditions: e bcd = 1 N c N c j=1 [ u (i) j u (e) j ] 2, u (inc) j ebcn = [ 1 N c 1 u N c ε (i) j=1 (i) j u (e) j (inc) u j ] 2, (31) where summation is performed using all the collocation points. Table 1: H polarization. The constant permittivity ε (i) = 1.2, the wave number 0 = 1. The maximal relative errors e in, e out inside and outside the scatterer. DOF M = 10 M = 20 e in e out e bcd e bcn e in e out e bc e bcn This example demonstrates two different inds of error in the approximate solution. The error of the first ind is caused by approximate satisfaction of the boundary condition on Ω. It is inherent in all boundary methods and it depends on the DOF which are used to approximate the boundary data. For a small number of harmonics M, increasing of the DOF decreases the boundary error e bc. However, this does not improve the accuracy of approximation of the scattered field. This can be explained by using of the source functions I(x ξ) instead of the Dirac functions δ(x ξ). In this case we do not have enough harmonics to obtain the appropriate condensed ones. When M increases, I(x ξ) approaches the delta function. As a result, the errors in computation of the scattered field have the same level as the ones in the boundary conditions. 2) Variable permittivity, E polarization. The relative permittivity inside Ω is taen as ε (i) (x, y) = 1 + β 2 (x 2 + y 2 ) 1 + β 2 r 2, β = const. (32) Here eqs. (2), (3) tae the form: 2 u (i) r r u (i) r + 1 r 2 2 u (i) ψ (1 + β 2 r 2 )u (i) = 0, r < a, (33) 2 u (e) r r u (e) r + 1 r 2 2 u (e) ψ u (e) = 0, r > a. (34) 8
9 The boundary conditions are: u (i) (a, ψ) u (e) (a, ψ) = J 0 ( 0 a) + 2 i m J m ( 0 a) cos mψ, u (i) u(e) (a, ψ) r r (a, ψ) = 0J 0( 0 a) i m J m( 0 a) cos mψ for all 0 ψ 2π. Similar to the example considered above this problem has an analytic solution in the form of expansion over the same system cos mψ: u (e) (r, ψ) = A m H m (2) ( 0 r) cos mψ, r > a, (35) u (i) (r, ψ) = m=0 B m V m (r 0, β) cos mψ, r < a, (36) m=0 where V m (r 0, β) is a solution of the equation [r 2 d2 dr 2 + r d ] dr m2 + 0(1 2 + β 2 r 2 )r 2 V m (r 0, β) = 0 bounded at the point r = 0. It can easily be verified that in this case V m (r/ 0 0, β) satisfies the equation [r 2 d2 dr 2 + r d ] dr m2 + (1 + β2 0 2 r 2 )r 2 V m (r/ 0 0, β) = 0. So, with a fixed m it depends on the only parameter α = β/ 0. We denote a solution of the equation bounded at the point ρ = 0 as Z m (ρ α). So we get [ρ 2 d2 dρ 2 + ρ d ] dρ m2 + (1 + α 2 ρ 2 )ρ 2 Z m = 0 (37) V m (r 0, β) = Z m ( 0 r β 0 ). The function Z m (ρ α) is looed for in the form of a power series: Z m (ρ α) = ρ m n=0 Substituting in (37) we get the recursive formula: b n (m, α)ρ 2n. (38) b 0 b 1 = 4(m + 1), b n = b n 1 + α 2 b n 2 4n(m + n) (39) with the free parameter b 0. Note that the value α = 0 corresponds to the scattering by a homogeneous circular cylinder. In this case the functions V m should be replaced by the Bessel s functions of the first ind J m. This can be provided by taing b 0 = 1 m!2 m (40) 9
10 because in this case Z m (ρ 0) = J m (ρ). (41) Numerical experiments show the rapid converges of (38) at least for ρ, α satisfying 0 ρ 0.5, 0 α 5. The summation is truncated as b n (m, α) becomes less than We use finite sums with the number of terms equal to 20 in evaluation (35), (36). This provides a truncating error less than in all of the calculations performed. The data in Table 2 were obtained for the wave number 0 = 2. The coefficient β = 1. In all of the calculations performed we use N = 6 (see (10)). This provides the maximal absolute error e a 10 6 in the approximation of ε (i) (x, y) inside Ω. Table 2: E polarization. Scattering by an inhomogeneous cylinder with the relative permittivity ε (i) (x, y) = 1 + x 2 + y 2 ; the wave number 0 = 2; the number of harmonics M = 20. error DOF= 10 DOF= 15 DOF= 20 e in e out ) Variable permittivity. H polarization. The relative permittivity inside Ω is given in (32). In this case eqs. (4), (5), written in polar coordinates, are: 2 u (e) r r u (e) r + 1 r 2 2 u (e) ψ u (e) = 0, r > a, (42) 2 u (i) r β2 r 2 u (i) r (1 + β 2 r 2 ) r u (i) r 2 ψ 2 + 0(1 2 + β 2 r 2 )u (i) = 0, r < a (43) with the boundary conditions: for all 0 ψ 2π. u (i) (a, ψ) u (e) (a, ψ) = J 0 ( 0 a) + 2 i m J m ( 0 a) cos mψ, 1 u (i) u(e) 1 + β 2 a 2 (a, ψ) r r (a, ψ) = 0J 0( 0 a) i m J m( 0 a) cos mψ Looing for the approximate solution inside the scatterer in the form (36) one gets the equation for V m : r 2 ( 1 + β 2 r 2) d 2 V m dr 2 + r ( 1 β 2 r 2) dv m dr + ( 1 + β 2 r 2) [ 2 0r 2 ( 1 + β 2 r 2) m 2] V m = 0 10
11 Denoting z = βr, c = 2 0/β 2 one gets: z 2 ( 1 + z 2) d 2 V m dz 2 + z ( 1 z 2) dv m dz + ( 1 + z 2) [ cz 2 ( 1 + z 2) m 2] V m = 0. (44) An approximate solution of (44) bounded at z = 0 is sought in the form of the series: V m (z, c) = b n (m, c) z m+2n. (45) The following recurrence formulae are obtained using the Mathematica 3 pacage: n=0 b n = cb n 3 + 2cb n 2 + [ 4 ( n 2 3n + 2 ) + c + (4n 6) m ] b n 1, 4n (m + n) b 0 = 1 2 m m!, b 1 = b 0 (2m c) 4(m + 1), b 2 = 2b 0c + b 1 (c + 2m). 8 (m + 2) As the numerical experiments shown the series (45) converge at least for z, c satisfying z 1, c 5. The summation is truncated as b n (m, c) becomes less than The unnown coefficients A m, B m are obtained as a solution of the sequence of 2 2 linear systems following from orthogonality of the cos mψ basis. In Table 3 the analytic solution is used to examine the one obtained by the numerical algorithm presented in the previous section. The wave number is: 0 = 1 and ν(x, y) = (x 2 + y 2 ). The values e in and e out are the same as the ones in Table 1. The data presented correspond to the same number of harmonics M = 20 in approximate solutions and the two different approximations of ν(x, y): with N = 1 and N = 5. Table 3: H polarization. Scattering by an inhomogeneous cylinder with ν(x, y) = 1+1.5(x 2 +y 2 ); the wave number 0 = 1; the number of harmonics M = 20. DOF N = 1 N = 5 e in e out e in e out One can see that for N = 1 corresponding to 9 terms expansion (18) increasing of DOF from 10 to 30 does not increase the exactness of the approximate solution. This means that the main 11
12 error in the solution is caused by the error in approximation of ν(x, y). This is the third ind of error of approximation in the method presented. When N increases, the differential operator of (19) approaches the initial one in (5) and the total error begin to decrease with the growth of DOF. 4 Conclusion The method presented can be regarded as a generalization of MDS onto the problems of scattering by electrical inhomogeneous bodies. When an accurate approximation of the complex permittivity ε (i) ( or of the inverse value 1/ε (i) ) is nown in the form of truncated Fourier series, then it provides solution of the scattering problem with a high precision. The method presented can be extended to 3D scattering. At least it seems to be quite simple in the scalar acoustic case. These extensions and 3D electromagnetic scattering are topics of the further investigations. References [1] T. Wriedt and A. Doicu. Formulations of the extended boundary condition method for three dimensional scattering using the method of discrete sources. J. Modern Opt., 45: , [2] T. Wriedt and U. Comberg. Comparison of computational scattering methods. J. Quant. Spectrosc. Radiat. Transfer, 60: , [3] M. I. Mishcheno, L. D. Travis, and D. W. Macowsi. T matrix computations of light scattering by nonspherical particles: a review. J. Quant. Spectrosc. Radiat. Transfer, 55: , [4] Ch. Hafner. The Generalized Multipole Technique for computational electromagnetics. Artech House Boos, Boston, [5] N. Kuster. Multiple Multipole Method for simulating EM problems involving biological bodies. IEEE Transactions on Biomedical Engineering, 40: , [6] Y. Leviatan and A. Boag. Analysis of electromagnetic scattering using a Current Model Method. Comput. Phys. Comm., 68: , [7] A. Doicu, Y. Eremin, and T. Wriedt. Acoustic and Electromagnetic Scattering Analysis using Discrete Sources. Academic Press, San Diego,
13 [8] S. Reutsiy and B. Tirozzi. Trefftz Spectral Method for elliptic equations of general type. Comp. Ass. Mech. and Eng. Sci., 8(4): , [9] S. Reutsiy and B. Tirozzi. A new boundary method for electromagnetic scattering from inhomogeneous bodies. J. Quant. Spectrosc. Radiat. Transfer, 72: , [10] S. Reutsiy. A boundary method of Trefftz type with approximate trial functions. Engineering Analysis with Boundary Elements, 37: , [11] S. Reutsiy and B. Tirozzi. Quasi Trefftz spectral method for separable linear elliptic equations. Comp. Math. with Applications, 37:47 74, [12] A. Bogomolny. Fundamental solutions method for elliptic boundary value problems. SIAM J. Numer. Anal., 22: ,
The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationNumerical computation of the Green s function for two-dimensional finite-size photonic crystals of infinite length
Numerical computation of the Green s function for two-dimensional finite-size photonic crystals of infinite length F.Seydou 1,Omar M.Ramahi 2,Ramani Duraiswami 3 and T.Seppänen 1 1 University of Oulu,
More informationOn Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material Surfaces
Vol. 114 2008) ACTA PHYSICA POLONICA A No. 6 A Optical and Acoustical Methods in Science and Technology On Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationElectromagnetic Scattering from an Anisotropic Uniaxial-coated Conducting Sphere
Progress In Electromagnetics Research Symposium 25, Hangzhou, China, August 22-26 43 Electromagnetic Scattering from an Anisotropic Uniaxial-coated Conducting Sphere You-Lin Geng 1,2, Xin-Bao Wu 3, and
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationEigenfunction Expansions of Source-excited Electromagnetic Fields on Open Cylindrical Guiding Structures in Unbounded Gyrotropic Media
128 Eigenfunction Expansions of Source-excited Electromagnetic Fields on Open Cylindrical Guiding Structures in Unbounded Gyrotropic Media A. V. Kudrin 1, E. Yu. Petrov 1, and T. M. Zaboronkova 2 1 University
More informationComparative Analysis of Techniques for Source Radiation in Cylindrical EBG with and without Periodic Discontinuities
1398 Progress In Electromagnetics Research Symposium Abstracts, St Petersburg, Russia, 22 25 May 2017 Comparative Analysis of Techniques for Source Radiation in Cylindrical EBG with and without Periodic
More informationApplication of Spatial Bandwidth Concepts to MAS Pole Location for Dielectric Cylinders
arquette University e-publications@arquette Electrical and Computer Engineering Faculty Research and Publications Engineering, College of --20 Application of Spatial Bandwidth Concepts to AS Pole Location
More informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationSEVERAL PRODUCTS OF DISTRIBUTIONS ON MANIFOLDS
Novi Sad J. Math. Vol. 39, No., 2009, 3-46 SEVERAL PRODUCTS OF DISTRIBUTIONS ON MANIFOLDS C. K. Li Abstract. The problem of defining products of distributions on manifolds, particularly un the ones of
More informationBEHAVIOR OF THE REGULARIZED SAMPLING INVERSE SCATTERING METHOD AT INTERNAL RESONANCE FREQUENCIES
Progress In Electromagnetics Research, PIER 38, 29 45, 2002 BEHAVIOR OF THE REGULARIZED SAMPLING INVERSE SCATTERING METHOD AT INTERNAL RESONANCE FREQUENCIES N. Shelton and K. F. Warnick Department of Electrical
More informationRECTANGULAR CONDUCTING WAVEGUIDE FILLED WITH UNIAXIAL ANISOTROPIC MEDIA: A MODAL ANALYSIS AND DYADIC GREEN S FUNCTION
Progress In Electromagnetics Research, PIER 5, 111 19, 000 RECTANGULAR CONDUCTING WAVEGUIDE FILLED WITH UNIAXIAL ANISOTROPIC MEDIA: A MODAL ANALYSIS AND DYADIC GREEN S FUNCTION S. Liu, L. W. Li, M. S.
More informationAn adaptive fast multipole boundary element method for the Helmholtz equation
An adaptive fast multipole boundary element method for the Helmholtz equation Vincenzo Mallardo 1, Claudio Alessandri 1, Ferri M.H. Aliabadi 2 1 Department of Architecture, University of Ferrara, Italy
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationA half submerged metal sphere (UIC comprehensive
Problem 1. exam) A half submerged metal sphere (UIC comprehensive A very light neutral hollow metal spherical shell of mass m and radius a is slightly submerged by a distance b a below the surface of a
More informationPhysics 504, Lecture 9 Feb. 21, 2011
Last Latexed: February 17, 011 at 15:8 1 1 Ionosphere, Redux Physics 504, Lecture 9 Feb. 1, 011 Let us return to the resonant cavity formed by the surface of the Earth a spherical shell of radius r = R
More informationGeneral review: - a) Dot Product
General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationScattering cross-section (µm 2 )
Supplementary Figures Scattering cross-section (µm 2 ).16.14.12.1.8.6.4.2 Total scattering Electric dipole, a E (1,1) Magnetic dipole, a M (1,1) Magnetic quardupole, a M (2,1). 44 48 52 56 Wavelength (nm)
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationElectromagnetic Scattering from a PEC Wedge Capped with Cylindrical Layers with Dielectric and Conductive Properties
0. OZTURK, ET AL., ELECTROMAGNETIC SCATTERING FROM A PEC WEDGE CAPPED WIT CYLINDRICAL LAYERS... Electromagnetic Scattering from a PEC Wedge Capped with Cylindrical Layers with Dielectric and Conductive
More informationClassical Scattering
Classical Scattering Daniele Colosi Mathematical Physics Seminar Daniele Colosi (IMATE) Classical Scattering 27.03.09 1 / 38 Contents 1 Generalities 2 Classical particle scattering Scattering cross sections
More informationME equations. Cylindrical symmetry. Bessel functions 1 kind Bessel functions 2 kind Modifies Bessel functions 1 kind Modifies Bessel functions 2 kind
Δϕ=0 ME equations ( 2 ) Δ + k E = 0 Quasi static approximation Dynamic approximation Cylindrical symmetry Metallic nano wires Nano holes in metals Bessel functions 1 kind Bessel functions 2 kind Modifies
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
More informationANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR THE THREE DIMENSIONAL TIME-HARMONIC MAXWELL PROBLEM
MATHEMATICS OF COMPUTATION Volume 77, Number 261, January 2008, Pages 1 10 S 0025-5718(07)02037-6 Article electronically published on September 18, 2007 ANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR
More informationScattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.
Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June
More informationECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 17
ECE 634 Intermediate EM Waves Fall 16 Prof. David R. Jacson Dept. of ECE Notes 17 1 General Plane Waves General form of plane wave: E( xz,, ) = Eψ ( xz,, ) where ψ ( xz,, ) = e j( xx+ + zz) The wavenumber
More informationAnalysis of eddy currents in a gradient coil
Analysis of eddy currents in a gradient coil J.M.B. Kroot Eindhoven University of Technology P.O.Box 53; 56 MB Eindhoven, The Netherlands Abstract To model the z-coil of an MRI-scanner, a set of circular
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More information. (70.1) r r. / r. Substituting, we have the following equation for f:
7 Spherical waves Let us consider a sound wave in which the distribution of densit velocit etc, depends only on the distance from some point, ie, is spherically symmetrical Such a wave is called a spherical
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationPHY752, Fall 2016, Assigned Problems
PHY752, Fall 26, Assigned Problems For clarification or to point out a typo (or worse! please send email to curtright@miami.edu [] Find the URL for the course webpage and email it to curtright@miami.edu
More informationAnalytical Study of Formulation for Electromagnetic Wave Scattering Behavior on a Cylindrically Shaped Dielectric Material
Research Journal of Applied Sciences Engineering and Technology 2(4): 307-313, 2010 ISSN: 2040-7467 Maxwell Scientific Organization, 2010 Submitted Date: November 18, 2009 Accepted Date: December 23, 2009
More informationAPPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION AND REFLECTION BY A CONDUCTING SHEET
In: International Journal of Theoretical Physics, Group Theory... ISSN: 1525-4674 Volume 14, Issue 3 pp. 1 12 2011 Nova Science Publishers, Inc. APPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION
More informationPropagation of Radio Frequency Waves Through Density Filaments
PSFC/JA-15-13 Propagation of Radio Frequency Waves Through Density Filaments A. K. Ram and K. Hizanidis a May 015 a National Technical University of Athens (part of HELLAS) School of Electrical and Computer
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationColumbia University Department of Physics QUALIFYING EXAMINATION
Columbia University Department of Physics QUALIFYING EXAMINATION Monday, January 9, 217 3:PM to 5:PM Classical Physics Section 2. Electricity, Magnetism & Electrodynamics Two hours are permitted for the
More informationAn Inverse Problem for Gibbs Fields with Hard Core Potential
An Inverse Problem for Gibbs Fields with Hard Core Potential Leonid Koralov Department of Mathematics University of Maryland College Park, MD 20742-4015 koralov@math.umd.edu Abstract It is well known that
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationAXIALLY SLOTTED ANTENNA ON A CIRCULAR OR ELLIPTIC CYLINDER COATED WITH METAMATERIALS
Progress In Electromagnetics Research, PIER 1, 329 341, 2 AXIALLY SLOTTED ANTENNA ON A CIRCULAR OR ELLIPTIC CYLINDER COATED WITH METAMATERIALS A-K. Hamid Department of Electrical/Electronics and Computer
More informationElectromagnetic (EM) Waves
Electromagnetic (EM) Waves Short review on calculus vector Outline A. Various formulations of the Maxwell equation: 1. In a vacuum 2. In a vacuum without source charge 3. In a medium 4. In a dielectric
More informationUniversity of Illinois at Chicago Department of Physics. Electricity & Magnetism Qualifying Examination
University of Illinois at Chicago Department of Physics Electricity & Magnetism Qualifying Examination January 7, 28 9. am 12: pm Full credit can be achieved from completely correct answers to 4 questions.
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationA method for creating materials with a desired refraction coefficient
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. A method
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More informationEXISTENCE OF GUIDED MODES ON PERIODIC SLABS
SUBMITTED FOR: PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS June 16 19, 2004, Pomona, CA, USA pp. 1 8 EXISTENCE OF GUIDED MODES ON PERIODIC SLABS Stephen
More informationAntennas and Propagation. Chapter 2: Basic Electromagnetic Analysis
Antennas and Propagation : Basic Electromagnetic Analysis Outline Vector Potentials, Wave Equation Far-field Radiation Duality/Reciprocity Transmission Lines Antennas and Propagation Slide 2 Antenna Theory
More informationPEMC PARABOLOIDAL REFLECTOR IN CHIRAL MEDIUM SUPPORTING POSITIVE PHASE VELOC- ITY AND NEGATIVE PHASE VELOCITY SIMULTANE- OUSLY
Progress In Electromagnetics Research Letters, Vol. 10, 77 86, 2009 PEMC PARABOLOIDAL REFLECTOR IN CHIRAL MEDIUM SUPPORTING POSITIVE PHASE VELOC- ITY AND NEGATIVE PHASE VELOCITY SIMULTANE- OUSLY T. Rahim
More informationSystem of Two Dielectric Cylinders Involving Charge Sources: I. Calculation of the Electric Field
Technical Physics, Vol. 5, No., 5, pp. 39 4. Translated from Zhurnal Tekhnicheskoœ Fiziki, Vol. 75, No., 5, pp.. Original Russian Text Copyright 5 by Emets. THEORETICAL AND MATHEMATICAL PHYSICS System
More informationModeling Focused Beam Propagation in scattering media. Janaka Ranasinghesagara, Ph.D.
Modeling Focused Beam Propagation in scattering media Janaka Ranasinghesagara, Ph.D. Teaching Objectives The need for computational models of focused beam propagation in scattering media Introduction to
More informationMATRIX REPRESENTATIONS FOR MULTIPLICATIVE NESTED SUMS. 1. Introduction. The harmonic sums, defined by [BK99, eq. 4, p. 1] sign (i 1 ) n 1 (N) :=
MATRIX REPRESENTATIONS FOR MULTIPLICATIVE NESTED SUMS LIN JIU AND DIANE YAHUI SHI* Abstract We study the multiplicative nested sums which are generalizations of the harmonic sums and provide a calculation
More informationWave Theory II (7) Physical Optics Approximation
Wave Theory II (7) Physical Optics Approximation Jun-ichi Takada (takada@ide.titech.ac.jp) In this lecture, the physical optics approximation (), which is classified as a semi-analytical technique, is
More informationMaxwell s equations in Carnot groups
Maxwell s equations in Carnot groups B. Franchi (U. Bologna) INDAM Meeting on Geometric Control and sub-riemannian Geometry Cortona, May 21-25, 2012 in honor of Andrey Agrachev s 60th birthday Researches
More informationUndulator Radiation Inside a Dielectric Waveguide
Undulator Radiation Inside a Dielectric Waveguide A.S. Kotanjyan Department of Physics, Yerevan State University Yerevan, Armenia Content Motivation On features of the radiation from an electron moving
More informationElectromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient
168 Int. J. Computing Science and Mathematics, Vol. Electromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient A.G. Ramm Mathematics epartment, Kansas
More informationl=0 The expansion coefficients can be determined, for example, by finding the potential on the z-axis and expanding that result in z.
Electrodynamics I Exam - Part A - Closed Book KSU 15/11/6 Name Electrodynamic Score = 14 / 14 points Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try
More informationwhere d is the vibration direction of the displacement and c is the wave velocity. For a fixed time t,
3 Plane waves 3.1 Plane waves in unbounded solid Consider a plane wave propagating in the direction with the unit vector p. The displacement of the plane wave is assumed to have the form ( u i (x, t) =
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationNumerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material
Chapter 30 Numerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material Hiroshi Maeda Additional information is available at the end of the chapter http://dx.doi.org/10.5772/50555
More informationScalar electromagnetic integral equations
Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how
More informationLight Scattering Group
Light Scattering Inversion Light Scattering Group A method of inverting the Mie light scattering equation of spherical homogeneous particles of real and complex argument is being investigated The aims
More informationElectromagnetic Theory I
Electromagnetic Theory I Final Examination 18 December 2009, 12:30-2:30 pm Instructions: Answer the following 10 questions, each of which is worth 10 points. Explain your reasoning in each case. Use SI
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationLecture 3 Fiber Optical Communication Lecture 3, Slide 1
Lecture 3 Optical fibers as waveguides Maxwell s equations The wave equation Fiber modes Phase velocity, group velocity Dispersion Fiber Optical Communication Lecture 3, Slide 1 Maxwell s equations in
More informationA cylinder in a magnetic field (Jackson)
Problem 1. A cylinder in a magnetic field (Jackson) A very long hollow cylinder of inner radius a and outer radius b of permeability µ is placed in an initially uniform magnetic field B o at right angles
More informationPhysics 505 Fall 2005 Homework Assignment #8 Solutions
Physics 55 Fall 5 Homework Assignment #8 Solutions Textbook problems: Ch. 5: 5., 5.4, 5.7, 5.9 5. A circular current loop of radius a carrying a current I lies in the x-y plane with its center at the origin.
More informationCONTROLLED CHANGES IN SPECTRA OF OPEN QUASI-OPTICAL RESONATORS
Progress In Electromagnetics Research B, Vol. 16, 85 105, 2009 CONTROLLED CHANGES IN SPECTRA OF OPEN QUASI-OPTICAL RESONATORS L. G. Velycho and Y. K. Sireno Institute of Radiophysics and Electronics National
More informationarxiv: v1 [cond-mat.mes-hall] 13 Feb 2019
Numerical Analysis of the Helmholtz Green s Function for Scalar Wave Propagation Through a Nano-hole on a Plasmonic Layer Désiré Miessein, 1, 2 Norman J. M. Horing, 2 Godfrey Gumbs, 3, 4 and Harry Lenzing
More informationUniversidad de Cantabria Facultad de Ciencias
Universidad de Cantabria Facultad de Ciencias ON LIGHT SCATTERING BY NANOPARTICLES WITH CONVENTIONAL AND NON-CONVENTIONAL OPTICAL PROPERTIES PH.D. THESIS Braulio García-Cámara Santander, July 2010 2 Theoretical
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More information1. Reminder: E-Dynamics in homogenous media and at interfaces
0. Introduction 1. Reminder: E-Dynamics in homogenous media and at interfaces 2. Photonic Crystals 2.1 Introduction 2.2 1D Photonic Crystals 2.3 2D and 3D Photonic Crystals 2.4 Numerical Methods 2.4.1
More informationMAT389 Fall 2016, Problem Set 4
MAT389 Fall 2016, Problem Set 4 Harmonic conjugates 4.1 Check that each of the functions u(x, y) below is harmonic at every (x, y) R 2, and find the unique harmonic conjugate, v(x, y), satisfying v(0,
More informationIN THE reconstruction of material shapes and properties, we
1704 IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 46, NO. 11, NOVEMBER 1998 Nonlinear Inversion in TE Scattering Bert Jan Kooij and Peter M. van den Berg Abstract A method for reconstructing
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationGeneralized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields
PHYSICAL REVIEW E 71, 5661 5 Generalized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields D. R. Lytle II Department of Electrical and Computer Engineering,
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationOptimization of Resonances in Photonic Crystal Slabs
Optimization of Resonances in Photonic Crystal Slabs Robert P. Lipton a, Stephen P. Shipman a, and Stephanos Venakides b a Louisiana State University, Baton Rouge, LA, USA b Duke University, Durham, NC,
More informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationMultiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions
Multiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions March 6, 2013 Contents 1 Wea second variation 2 1.1 Formulas for variation........................
More informationMemoirs on Differential Equations and Mathematical Physics
Memoirs on Differential Equations and Mathematical Physics Volume 34, 005, 1 76 A. Dzhishariani APPROXIMATE SOLUTION OF ONE CLASS OF SINGULAR INTEGRAL EQUATIONS BY MEANS OF PROJECTIVE AND PROJECTIVE-ITERATIVE
More informationVolume and surface integral equations for electromagnetic scattering by a dielectric body
Volume and surface integral equations for electromagnetic scattering by a dielectric body M. Costabel, E. Darrigrand, and E. H. Koné IRMAR, Université de Rennes 1,Campus de Beaulieu, 35042 Rennes, FRANCE
More informationNOTE. Application of Contour Dynamics to Systems with Cylindrical Boundaries
JOURNAL OF COMPUTATIONAL PHYSICS 145, 462 468 (1998) ARTICLE NO. CP986024 NOTE Application of Contour Dynamics to Systems with Cylindrical Boundaries 1. INTRODUCTION Contour dynamics (CD) is a widely used
More informationOn the Spectrum of Volume Integral Operators in Acoustic Scattering
11 On the Spectrum of Volume Integral Operators in Acoustic Scattering M. Costabel IRMAR, Université de Rennes 1, France; martin.costabel@univ-rennes1.fr 11.1 Volume Integral Equations in Acoustic Scattering
More informationRadio Propagation Channels Exercise 2 with solutions. Polarization / Wave Vector
/8 Polarization / Wave Vector Assume the following three magnetic fields of homogeneous, plane waves H (t) H A cos (ωt kz) e x H A sin (ωt kz) e y () H 2 (t) H A cos (ωt kz) e x + H A sin (ωt kz) e y (2)
More informationThe concept of perfect electromagnetic conductor (PEMC) has been defined by medium conditions of the form [1, 2]
Progress In Electromagnetics Research B, Vol. 5, 169 183, 2008 REFLECTION AND TRANSMISSION OF WAVES AT THE INTERFACE OF PERFECT ELECTROMAGNETIC CONDUCTOR PEMC I. V. Lindell and A. H. Sihvola Electromagnetics
More informationSound 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 informationOperation of Evanescent Wave Intensity Using Metamaterials of Negative Permittivity and Permeability
PIERS ONLINE, VOL. 3, NO. 8, 2007 3 Operation of Evanescent Wave Intensity Using Metamaterials of Negative Permittivity and Permeability Yu. N. Barabanenkov, M. Yu. Barabanenkov 2, and S. A. Nikitov Institute
More informationCLASSICAL ELECTRICITY
CLASSICAL ELECTRICITY AND MAGNETISM by WOLFGANG K. H. PANOFSKY Stanford University and MELBA PHILLIPS Washington University SECOND EDITION ADDISON-WESLEY PUBLISHING COMPANY Reading, Massachusetts Menlo
More informationGraduate Diploma in Engineering Circuits and waves
9210-112 Graduate Diploma in Engineering Circuits and waves You should have the following for this examination one answer book non-programmable calculator pen, pencil, ruler No additional data is attached
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationPDEs in Spherical and Circular Coordinates
Introduction to Partial Differential Equations part of EM, Scalar and Vector Fields module (PHY2064) This lecture Laplacian in spherical & circular polar coordinates Laplace s PDE in electrostatics Schrödinger
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationApproximation Methods in QM
Chapter 3 Approximation Methods in QM Contents 3.1 Time independent PT (nondegenerate)............... 5 3. Degenerate perturbation theory (PT)................. 59 3.3 Time dependent PT and Fermi s golden
More informationA Spectral Time-Domain Method for Computational Electrodynamics
A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation
More information