Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 30 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
30
Dung lượng
1,96 MB
Nội dung
Electromagnetic Wave Scattering from Material Objects Using HybridMethods 21 4.5 5 5.5 6 6.5 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f [GHz] TE ξ =0 ° ξ=45 ° ξ=90 ° 4.5 5 5.5 6 6.5 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f [GHz] ξ =0 ° ξ=45 ° ξ=90 ° TM Fig. 23. Power reflection coefficients of the fundamental space harmonics versus frequency for one-hundred-layered square lattice periodic arrays of metallic posts embedded in dielectric cylinders from Fig 22. Parameters of the structure: h = 19.5mm, d = 25mm, r 0 = 0.01h, r 1 = 0.09h, R = 0.19h, ψ 1,2 = 10 ◦ , ε r = 20 structures are identical but rotated by 90 ◦ with respect to each other. For the same plane wave illuminating both configurations, they produce stop bands which only slightly overlap. When half of the structure (i.e. 10 last or first arrays) are being rotated with respect to the other half one obtains the effect of stop band shifting. The stop bands, which are almost identical in width, can be shifted from one bandwidth to another. The case of 90 ◦ rotation of stacks is presented in Fig. 24(b). When only every other periodic array are being rotated the produced stop band is widening andin the case of 90 ◦ rotation it embraces both stop bands as can be seen in Fig. 24(c). 3.3.4 Tunneling effect An interesting effect of wave tunneling can be obtained in the structure under investigation. This effect, along with the "growing evanescent envelope" for field distributions, was previously observed in metamaterial medium (negative value of real permittivity and permeability) and a structure composed of a pair of only-epsilon-negative and only-mu-negative layers Alu & Engheta (2003). This effect was also discussed in Alu & Engheta (2005) for periodically layered stacks of frequency selective surfaces (FSS). It was shown in Alu & Engheta (2005) that a complete electromagnetic wave tunneling may be achieved through a pair ofdifferent stacked FSSs which are characterized by dual behaviors, even though each stack is completely alone opaque (operates in its stop band). Similar effect can be obtained for the structure composed of a pair of identical stacks of periodic arrays of cylindrical posts rotated by 90 ◦ with respect to each other. This effect can also be controlled by introducing a gap d between the stacks (see Fig. 25). The calculation of a total scattering matrix for a pair of such stacks boils down to cascading the scattering matrix of a stack calculated for TE wave excitation with the scattering matrix calculated for TM wave excitation. The tunneling effect has been obtained for the periodic structure described in Fig. 25. The stop bands are formed in the same frequency range for both TE and TM waves. Therefore, we obtain a pair of stacks with dual behavior both of which operate in their stop bands. In the equivalent circuit analogy one stack is represented by a periodical line loaded with capacitors, while the other one is loaded with inductances. 47 Electromagnetic Wave Scattering from Material Objects Using Hybrid Methods 22 Will-be-set-by-IN-TECH 5 5.5 6 6.5 7 7.5 8 8.5 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 f[GHz] 5 5.5 6 6.5 7 7.5 8 8.5 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] h l plane wave H k E R 0 R 0 plane wave H k E plane wave H k E l h plane wave H k E a) b) c) Fig. 24. Power reflection coefficients of the fundamental space harmonics versus frequency for normal incidence of TE wave on a periodic structures; Parameters of the structures h = 20mm, l = h, r = 0.06h , R = 0.35h, ε r1 = 3, ε r2 = 2.5, number of sections 20. Fig. 25. Schematic 3-D representation of a periodic structure under investigation Fig. 26 illustrates the power reflection coefficients for the normal incidence of a TE polarized plane wave on stacks of periodic structuresof dielectric cylinders with double dielectric inclusions. The characteristics for scattering from structure 1, 2 and3 are illustrated. The results show that stacks 1 and 2 are completely opaque in presented frequency ranges. However, when half of these configurations are rotated by 90 ◦ with respect to the other half, forming the structure 3, the tunneling effect can be observed. The obtained configurations enable the signal from a very narrow frequency range to tunnel through the structure. This tunneling effect can be controlled by adjusting the distance d between stacks (see Fig. 25). Fig. 27 shows the characteristics of the power reflection coefficients for the normal incidence of a TE polarized plane wave on structure 3 for different values of distance d.Itcanbeclearly seen that this value is directly connected to the frequency of the tunneled wave. 48 BehaviourofElectromagneticWavesinDifferentMediaandStructuresElectromagnetic Wave Scattering from Material Objects Using HybridMethods 23 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] structure 1 structure 2 structure 3 plane wave H k E plane wave H k E plane wave H k E single unit cell Fig. 26. Power reflection coefficients of the fundamental space harmonics versus frequency for normal incidence of TE wave on a periodic structures; Parameters of the structures: h = 20mm, l = 20mm, d = l, R = 0.48h, ε r = 1.5, inclusion - two dielectric cylinders r = 0.16h, ε rc = 2, displacement from the center .24h number of sections 20; 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] 5.4 5.6 5.8 6 6.2 6.4 6.6 6.8 7 7.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R 0 f[GHz] d=1.2l d=1.4l d=1.6l d=1.8l d=2l d=2.2l Fig. 27. Power reflection coefficients of the fundamental space harmonics versus frequency for normal incidence of TE wave on a periodic structures described in Fig. 26 for different values of distance d between stacks; dashed line (red) - structure 1, dash-dot line (green) - structure 2; solid line (blue) - structure 3 49 Electromagnetic Wave Scattering from Material Objects Using Hybrid Methods 24 Will-be-set-by-IN-TECH 4. Conclusion In the chapter a hybrid method ofelectromagnetic wave scattering from structures containing complex cylindrical or spherical objects is presented. Depending of the investigated post geometry different numerical techniques were utilized such as mode-matching technique, method of moments and finite difference method defined in the frequency domain. The proposed approach enables to determine the scattering parameters of open and closed structures containing the configuration of cylindrical objects of arbitrary cross-section and axially symmetrical posts. The proposed technique rests on defining the collateral cylindrical or spherical object containing the investigated element and then utilizing the analytical iterative model for determining scattering parameters of arbitrary configuration of objects. The obtained solution can be combined with the arbitrary external excitation which allows analyzing the variety of open and closed microwave structures. The convergence of the method have been analyzed during the numerical studies. Additionally, in order to verify the correctness of the developed method the research of a number of open and closed microwave structures such as beam shaping configurations, resonators, filters and periodic structures have been conducted. The obtained numerical results have been verified by comparing them with the ones obtained form alternative numerical methods or own measurements. A good agreement between obtained results was achieved. 5. Acknowledgement This work was supported inpart by the Polish Ministry of Science and Higher Education under Contract N515 501740, decision No 5017/B/T02/2011/40 andinpart from sources for science in the years 2010-2012 under COST Action IC0803, decision No 618/N-COST/09/2010/0. 6. References Aiello, G., Alfonzetti, S. & Dilettoso, E. (2003). Finite-element solution of eddy-current problems in unbounded domains by means of the hybrid FEM-DBCI method, MAGN 39(3): 1409 – 1412. Alessandri, F., Giordano, M., Guglielmi, M., Martirano, G. & Vitulli, F. (2003). A new multiple-tuned six-port Riblet-type directional coupler in rectangular waveguide, IEEE Trans. Microw. Theory Tech. 51(5): 1441 – 1448. Alu, A. & Engheta, N. (2003). Pairing an epsilon-negative slab with a mu-negative slab: Resonance, tunneling and transparency, IEEE Trans. Antennas Propag. 51(10): 2558–2571. Alu, A. & Engheta, N. (2005). Evanescent growth and tunneling through stacks of frequency-selective surfaces, IEEE An tennas Wireless Propag. Letters 4: 417–420. Arndt, F., Catina, V. & Brandt, J. (2004). Efficient hybrid MM/MoM technique for the CAD of circular combline filters with resonators of more general shape, IEEE Microwave Symposium Digest 3: 1407 – 1410. Aza, G., Encinar, J. A., Zapata, J. & Lambea, M. (1998). Full-wave analysis of cavity-backed and probe-fed microstrip patch arrays by a hybrid mode-matching generalized scattering matrix and finite-element method, IEEE Trans. Antennas Propag. 46(2): 234 – 242. Dittloff, J., Arndt, F. & Grauerholz, D. (1988). Optimum design of waveguide E-plane stub-loaded phase shifters, IEEE Trans. Microw. Theory Tech. 36(3): 582–587. 50 BehaviourofElectromagneticWavesinDifferentMediaandStructuresElectromagnetic Wave Scattering from Material Objects Using HybridMethods 25 Elsherbeni, A. Z., Hamid, M. & Tian, G. (1993). Iterative scattering of a Gaussian beam by an array of circular conducting and dielectric cylinders, J. of Electromag. Waves App. 7(10): 1323–1342. Esteban, H., Cogollos, S., Boria, V., Blas, A. S. & Ferrando, M. (2002). A new hybrid mode-matching/numerical method for the analysis of arbitrarily shaped inductive obstacles and discontinuities in rectangular waveguides, IEEE Trans. Microw. Theory Tech. 50(4): 1219–1224. Gimeno, B., Cruz, J. L., Navarro, E. A. & Such, V. (1994). A polarizer rotator system for three-dimensional oblique incidence, IEEE Trans. Antennas Propag. 42(7): 912–919. Hamid, A K., Ciric, I. & Hamid, M. (1991). Iterative solution of the scattering by an arbitrary configuration of conducting or dielectric spheres, IEE Proc. H Microw., Antennas Propag. 138(6): 565 – 572. Kildal, P S., Kishk, A. A. & Tengs, A. (1996). Reduction of forward scattering from cylindrical objects using hard surfaces, IEEE Trans. Antennas Propag. 44: 1509–1520. Kushta, T. & Yasumoto, K. (2000). Electromagnetic scattering from periodic arrays of two circular cylinders per unit cell, Prog. Electromag. Research 29: 69–85. Kusiek, A. & Mazur, J. (2009). Analysis of scattering from arbitrary configuration of cylindrical objects using hybrid finite-difference mode-matching method, Prog. Electromag. Research 97: 105 – 127. Kusiek, A. & Mazur, J. (2010). Hybrid finite-difference/mode-matching method for analysis of scattering from arbitrary configuration of rotationally-symmetrical posts, Prog. Electromag. Research 110: 23 – 42. Kusiek, A. & Mazur, J. (2011). Application of hybrid finite-difference mode-matching method to analysis ofstructures loaded with axially-symmetrical posts, Microw. Opt. Tech. Lett. 53(1): 189–194. Lech, R. & Mazur, J. (2007). Analysis of circular cavity with cylindrical objects, IEEE Trans. Microw. Theory Tech. 55(10): 2115 – 2123. Lech, R., Polewski, M. & Mazur, J. (2006). Electromagnetic wave scattering by a set of axially symmetric ω-medium cylinder, Microw. Opt. Tech. Lett. 48(9): 1880 – 1884. Mrozowski, M. (1994). A hybrid PEE-FDTD algorithm for accelerated time domain analysis ofelectromagnetic waves, IEEE Microw. Guided Wave Lett. 4: 323 – 325. Mrozowski, M., Okoniewski, M. & Stuchly, M. (1996). A hybrid PEE-FDTD method for efficient field modeling in cyllindrical coordinates, Elect. Lett. 32: 194 – 195. Polewski, M., Lech, R. & Mazur, J. (2004). Rigorous modal analysis ofstructures containing inhomogeneous dielectric cylinders, IE EE Trans. Microw. Theory Tech. 52(5): 1508–1516. Polewski, M., Lech, R. & Mazur, J. (2006). Radiation phenomena from pseudochiral cylinders at plane wave incidence, Microw. Opt. Tech. Lett. 48(9): 1880 – 1884. Polewski, M. & Mazur, J. (2002). Scattering by an array of conducting lossy dielectric, ferrite and pseudochiral cylinders, Prog. Electromag. Research pp. 283–310. Quick Wave 3D (QWED) (n.d.). http://www.qwed.com.pl/. Rogier, H. (1998). A new hybrid FDTD-BIE approach to model electromagnetic scattering problems, IEEE Microw. Guided Wave Lett. 8(3): 138 – 140. Roy, T., Sarkar, T. K., DjordjeviC, A. R. & Salazar-Palma, M. (1996). A hybrid method solution of scattering by conducting cylinders (TM case), IEEE Trans. Microw. Theory Tech. 44(12): 2145 – 2151. 51 Electromagnetic Wave Scattering from Material Objects Using Hybrid Methods 26 Will-be-set-by-IN-TECH Rubio, J., Arroyo, J. & Zapata, J. (1999). Analysis of passive microwave circuits by using a hybrid 2-D and 3-D finite-element mode-matching method, IEEE Trans. Microw. Theory Tech. 47(9): 1746 – 1749. Sabbagh, M. E. & Zaki, K. (2001). Modeling of rectangular waveguide junctions containing cylindrical posts, Prog. Electromag. Research 33: 299–331. Sharkawy, M. A., Demir, V. & Elsherbeni, A. Z. (2006). Plane wave scattering from three dimensional multiple objects using the iterative multiregion technique based on the FDFD method, IEEE Trans. Antennas Pr opag. 54(2): 666–673. Shen, T., Zaki, K. A. & Wang, C. (2000). Tunable dielectric resonators with dielectric tuning disks, IEEE Trans. Microw. Theory Tech. 48: 2439 ˝ U2445. Toyama, H. & Yasumoto, K. (2005). Electromagnetic scattering from periodic arrays of composite circular cylinder with internal cylindrical scatterers, Prog. Electromag. Research 52: 355–358. Uher, J., Bornemann, J. & Rosenberg, U. (1993). Waveguide Components for Antenna Feed Systems: Theory and CAD, Artech House Antennas and Propagation Library, Norwood. Waterman, P. C. (1971). Symmetry, unitarity and geometry inelectromagnetic scattering, Physical Review D 3: 825–839. Xu, F. & Hong, W. (2004). Analysis of two dimensions sparse multicylinder scattering problem using DD-FDTD method, IEEE Trans. Antennas Propag. 52(10): 2612 – 2617. Yasumoto, K. (2000). Generalized method for electromagnetic scattering by two-dimensional periodic discrete composites using lattice sums, Proc. Int. Conf. Microwave and Millimeter Wave Technology (9): 29–34. Yasumoto, K., Toyama, H. & Kushta, T. (2004). Accurate analysis of two-dimensional electromagnetic scattering from multilayered periodic arrays of circular cylinders using lattice sums technique, IEEE Trans. Antennas Propag. 52: 2603–2611. Yasumoto, K. & Yoshitomi, K. (1999). Efficient calculation of lattice sums for freespace periodic green’s function, IEEE Trans. Antennas Propag. 47(6): 1050–1055. 52 BehaviourofElectromagneticWavesinDifferentMediaandStructures 4 The Eigen Theory ofElectromagneticWavesin Complex Media Shaohua Guo Zhejiang University of Science and Technology P. R. China 1. Introduction Since J. C. Maxwell presented the electromagnetic field equations in 1873, the existence ofelectromagneticwaves has been verified in various medium (Kong, 1986; Monk, 2003). But except for Helmholtz’s equation ofelectromagneticwavesin isotropic media, the laws of propagation ofelectromagneticwavesin anisotropic media are not clear to us yet. For example, how many electromagneticwaves are there in anisotropic media? How fast can these electromagneticwaves propagate? Where are propagation direction and polarization direction of the electromagnetic waves? What are the space patterns of these waves? Although many research works were made in trying to deduce the equations ofelectromagneticwavesin anisotropic media based on the Maxwell’s equation (Yakhno, 2005, 2006; Cohen, 2002; Haba, 2004), the explicit equations ofelectromagneticwavesin anisotropic media could not be obtained because the dielectric permittivity matrix and magnetic permeability matrix were all included in these equations, so that only local behaviourofelectromagnetic waves, for example, in a certain plane or along a certain direction, can be studied. On the other hand, it is a natural fact that electric and magnetic fields interact with each other in classical electromagnetics. Therefore, even if most of material studies deal with the properties due to dielectric polarisation, magneitc materials are also capable of producing quite interesting electro-magnetic effects (Lindellm et al., 1994). From the bi-anisotropic point of view, magnetic materials can be treated as a subclass of magnetoelectric materials. The linear constitutive relations linking the electric and magnetic fields to the electric and magnetic displacements contain four dyadics, three of which have direct magnetic contents. The magnetoelectric coupling has both theoretical and practical significance in solid state physics and materials science. Though first predicted by Pierre Curie, magnetoeletric coupling was originally through to be forbidden because it violates time-reversal symmetry, until Laudau and Lifshitz (Laudau & Lifshitz, 1960) pointed out that time reversal is not a symmetry operation in some magnetic crystal. Based on this argument, Dzyaloshinskii (Dzyaloshinskii, 1960) predicted that magnetoelectric effect should occur in antiferromagnetic crystal Cr 2 O 3 , which was verified experimentally by Astrov (Astrov, 1960). Since then the magnetoelectric coupling has been observed in single-phase materials where simultaneous electric and magnetic ordering coexists, andin two-phase composites where the participating phase are pizoelectric and piezomagnetic (Bracke & Van Vliet,1981; Van Run et al., 1974) . Agyei and Birman (Agyei & Birman, 1990) carried out a detailed BehaviourofElectromagneticWavesinDifferentMediaandStructures 54 analysis of the linear magnetoelectric effect, which showed that the effect should occur not only in some magnetic but also in some electric crystals. Pradhan (Pradhan, 1993) showed that an electric charge placed in a magnetoelectric medium becomes a source of induced magnetic field with non-zero divergence of volume integral. Magnetoelectric effect in two- phase composites has been analyzed by Harshe et al. ( Harshe et al., 1993), Nan (Nan, 1994) and Benveniste (Benveniste, 1995). Broadband transducers based on magnetoelectric effect have also been developed (Bracke & Van Vliet, 1981). Although the development mentioned above, no great progress in the theories ofelectromagneticwavesin bi-anisotropic media because of the difficulties in deal with the bi-coupling in electric field and magnetic one of the Maxell’s equation and the bi-anisotropic constitutive equation by classical electromagnetic theory. Recently there is a growing interest modeling and analysis of Maxwell’s equations (Lee & Madsen, 1990; Monk, 1992; Jin et al., 1999). However, most work is restricted to simple medium such as air in the free space. On the other hand, we notice that lossy and dispersive media are ubiquitous, for example human tissue, water, soil, snow, ice, plasma, optical fibers and radar-absorbing materials. Hence the study of how electromagnetic wave interacts with dispersive media becomes very important. Some concrete applications include geophysical probing and subsurface studied of the moon and other planets (Bui et al., 1991), High power and ultra-wide-band radar systems, in which it is necessary to model ultra-wide-band electromagnetic pulse propagation through plasmas (Dvorak & Dudley, 1995), ground penetrating radar detection of buried objects in soil media (liu & Fan, 1999). The Debye medium plays an important role inelectromagnetic wave interactions with biological and water-based substances (Gandhi & Furse, 1997). Until 1990, some paper on modeling of wave propagation in dispersive media started making their appearance in computational electromagnetics community. However, the published papers on modeling of dispersive media are exclusively restricted to the finite-difference time-domain methods and the finite element methods (Li & Chen, 2006; Lu et al., 2004). To our best knowledge, there exist only few works in the literature, which studied the theoretical model for the Maxwell’s equation in the complex anisotropic dispersive media, and no explicit equations ofelectromagneticwavesin anisotropic dispersive media can be obtained due to the limitations of classical electromagnetic theory. Chiral materials have been recently an interesting subject. In a chiral medium, an electric or magnetic excitation will produce simultaneously both electric and magnetic polarizations. On the other hand, the chiral medium is an object that cannot be brought into congruence with its mirror image by translation and rotation. Chirality is common in a variety of naturally occurring and man-made objects. From an operation point of view, chirality is introduced into the classical Maxwell equations by the Drude-Born-Fedorov relative constitutive relations in which the electric and magnetic fields are coupled via a new materials parameter (Lakhtakia, 1994; Lindell et al., 1994), the chirality parameter. These constitutive relations are chosen because they are symmetric under time reversality and duality transformations. In a homogeneous isotropic chiral medium the electromagnetic fields are composed of left-circularly polarized (LCP) and right- circularly polarized (RCP) components (Jaggard et al., 1979; Athanasiadis & Giotopoulos, 2003), which have different wave numbers and independent directions of propagation. Whenever an electromagnetic wave (LCP, RCP or a linear combination of them) is incident upon a chiral scatterer, then the scattered field is composed of both LCP and RCP components and therefore both LCP and RCP far-field patterns are derived. Hence, in the vector problem we need to specify two The Eigen Theory ofElectromagneticWavesin Complex Media 55 directions of propagation and two polarizations. In recent years, chiral materials have been increasingly studied and there is a growing literature covering both their applications and the theoretical investigation of their properties. It will be noticed that the works dealing with wave phenomena in chiral materials have been mainly concerned with the study of time- harmonic waves which lead to frequency domain studies (Lakhtakia et al., 1989; Athanasiadis et al., 2003). In this chapter, the idea of standard spaces is used to deal with the Maxwell’s electromagnetic equation (Guo, 2009, 2009, 2010, 2010, 2010). By this method, the classical Maxwell’s equation under the geometric presentation can be transformed into the eigen Maxwell’s equation under the physical presentation. The former is in the form of vector and the latter is in the form of scalar. Through inducing the modal constitutive equations of complex media, such as anisotropic media, bi-anisotropic media, lossy media, dissipative media, and chiral media, a set of modal equations ofelectromagneticwaves for all of those media are obtained, each of which shows the existence ofelectromagnetic sub-waves, meanwhile its propagation velocity, propagation direction, polarization direction and space pattern can be completely determined by the modal equations.This chapter will make introductions of the eigen theory to reader in details. Several novel theoretical results were discussed in the different parts of this chapter. 2. Standard spaces ofelectromagneticmediaIn anisotropic electromagnetic media, the dielectric permittivity and magnetic permeability are tensors instead of scalars. The constitutive relations are expressed as follows , =⋅ = ⋅D ε EBμ H (1) Rewriting Eq.(1) in form of scalar, we have , == ε μ ii jj ii jj DEB H (2) where the dielectric permittivity matrix ε and the magnetic permeability matrix μ are usually symmetric ones, and the elements of the matrixes have a close relationship with the selection of reference coordinate. Suppose that if the reference coordinates is selected along principal axis of electrically or magnetically anisotropic media, the elements at non-diagonal of these matrixes turn to be zero. Therefore, equations (1) and (2) are called the constitutive equations ofelectromagneticmedia under the geometric presentation. Now we intend to get rid of effects of geometric coordinate on the constitutive equations, and establish a set of coordinate-independent constitutive equations ofelectromagneticmedia under physical presentation. For this purpose, we solve the following problems of eigen-value of matrixes. () () , −− λγ II ε φ=0 μ ϕ=0 (3) where () 1,2,3= λ i i and () 1,2,3= γ i i are respectively eigen dielectric permittivity and eigen magnetic permeability, which are constants of coordinate-independent. () 1,2,3= φ i i and () 1,2,3= ϕ i i are respectively eigen electric vector and eigen magnetic vector, which show the electrically principal direction and magnetically principal direction of anisotropic media, and are all coordinate-dependent. We call these vectors as standard spaces. Thus, the matrix BehaviourofElectromagneticWavesinDifferentMediaandStructures 56 of dielectric permittivity and magnetic permeability can be spectrally decomposed as follows , ΤΤ == εΦΛΦ μ ΨΠΨ (4) where [ ] 123 ,,= λλλΛ diag and [ ] 123 ,,= γγγ Π diag are the matrix of eigen dielectric permittivity and eigen magnetic permeability, respectively. { } 123 ,,= Φ φφφ and { } 123 ,,= Ψ ϕϕϕ are respectively the modal matrix of electric mediaand magnetic media, which are both orthogonal and positive definite matrixes, and satisfy T = I ΦΦ , T = I ΨΨ . Projecting the electromagnetic physical qualities of the geometric presentation, such as the electric field intensity vector E , magnetic field intensity vector H , magnetic flux density vector B and electric displacement vector D into the standard spaces of the physical presentation, we get Τ =⋅ * D D Φ , Τ =⋅ * E E Φ (5) Τ =⋅ * BB Ψ , Τ =⋅ * HH Ψ (6) Rewriting Eqs.(5) and (6) in the form of scalar, we have * Τ =⋅ i = 1,2,3D φ ii D , * Τ =⋅ i = 1,2,3E φ ii E (7) * Τ =⋅ i = 1,2,3B ϕ ii B , * Τ =⋅ i = 1,2,3H ϕ ii H (8) These are the electromagnetic physical qualities under the physical presentation. Substituting Eq. (4) into Eq. (1) respectively, and using Eqs.(5) and (6) yield ** = i = 1,2,3 λ iii DE (9) ** = i = 1,2,3 γ iii BH (10) The above equations are just the modal constitutive equations in the form of scalar. 3. Eigen expression of Maxwell’s equation The classical Maxwell’s equations in passive region can be written as ×=∇HD ∇ t , ×=−∇ E B ∇ t (11) Now we rewrite the equations in the form of matrix as follows 11 22 33 0 0 0 −∂ ∂ ∂−∂ =∇ −∂ ∂ zy zx t yx HD HD HD (12) or [ ] { } { } Δ=∇ t HD (13) [...]... only one kind ofelectromagnetic wave in isotropic crystal, which is identical with the classical result; there are two kinds ofelectromagneticwavesin uniaxial crystal; three kinds ofelectromagneticwavesin biaxial crystal and three kinds of distorted 72 BehaviourofElectromagneticWavesinDifferentMediaandStructureselectromagneticwavesin monoclinic crystal Also for bi-anisotropic media, ... W1( 3) [φ1 ,φ2 , 3 ] , Wele = W1( ) [ϕ 1 ,ϕ 2 ,ϕ 3 ] 3 (87) 33 T T {1,1,1} , ϕ 1* = {1,1,1} , ξ1 = 1 33 Then the eigen-qualities and eigen-operators of an isotropic lossy media are respectively shown as follows where, φ1* = * E1 = 3 ( E1 + E2 + E3 ) 3 (88) 64 BehaviourofElectromagneticWavesinDifferentMediaandStructures * H1 = 3 ( H1 + H2 + H3 ) 3 1 − ∂2 + ∂2 + ∂2 x y z 3 ( 1* = ) (89)... a diagonal matrix Thus Eqs (33 ) and (34 ) can be uncoupled in the form of scalar T i* H i* + ε 0γ i∇ t2 H i* = 0 i = 1, 2, 3 (35 ) i* Ei* + ε 0γ i∇ t2 Ei* = 0 i = 1, 2, 3 (36 ) Eqs. (35 ) and (36 ) are the modal equations ofelectromagneticwavesin anisotropic magnetics 59 The Eigen Theory ofElectromagneticWavesin Complex Media 5 Electromagneticwavesin bi-anisotropic media 5.1 Bi-anisotropic constitutive... fully dynamical theory of piezoelectromagnetic waves Acta Mech., Vol 215, No 3, (December 2010), pp 33 5 -34 4, ISSN 0001-5970 [17] Haba, Z (2004) Green functions and propagation ofwavesin strongly inhomogeneous media Journal of Physis A: Mathematical and General, Vol 37 , No 9, (September 2004), pp 9295- 930 2, ISSN 030 5-4470 The Eigen Theory ofElectromagneticWavesin Complex Media 75 [18] Harshe, G.;... waves are E1 = Ae − k2 x ⋅ e ( i k1 x − ω t ) = A⋅e ( i k1 x − ω t ) (1 23) It is an attenuated sub -waves 8 Electromagneticwavesin chiral media 8.1 The constitutive equation of chiral media The constitutive equations of chiral media are the following D = ε ⋅ E - χ ⋅ ∇t H (124) B = χ ⋅ ∇t E + μ ⋅ H (125) 68 BehaviourofElectromagneticWavesinDifferentMediaandStructures where χ is the matrix of. .. (December 1974), pp 1710-1714, ISSN 0022-2461 76 BehaviourofElectromagneticWavesinDifferentMediaandStructures [35 ] Yakhno, V.; Yakhno, T & Kasap, M (2006) A novel approach for modeling and simulation ofelectromagneticwavesin anisotropic dielectrics Internatioal Journal of Solids and Structures, Vol 43, No 12, (December 2006), pp 6261-6276, ISSN 002076 83 ... two electromagneticwavesin Dzyaloshinskii’s bi-anisotropic media, and the electromagneticwavesin bi-anisotropic medium will go faster duo to the bi-coupling between electric field and magnetic one For isotropic dispersive medium, the electromagnetic wave is an attenuated sub -waves And for chiral medium, there exist different propagating states of electromagnetivc wavesindifferent frequency band,... , ς 3 ] Then we have Di* = λi Ei* − ς i∇ t H i* ( 130 ) Bi* = ς i ∇t Ei* + γ i H i* ( 131 ) Eqs.( 130 ) and ( 131 ) are just the modal constitutive equations for anisotropic chiral media 8.2 Eigen equations ofelectromagneticwavesin chiral media Substituting Eqs ( 130 ) and ( 131 ) into Eqs ( 23) and (24), respectively, we have {Δ } E * i * i {Δ } H * i = −∇t {ϕi } (ς i ∇t Ei* + γ i H i* ) i = 1, 2, 3 ( 132 )... (106) Eqs (105) and (106) are just the modal constitutive equations for the general dispersive media 7.2 Eigen equations ofelectromagneticwavesin dispersive media Substituting Eqs (105) and (106) into Eqs ( 23) and (24), respectively, we have {Δ } E * i * i ( 1 2 = −∇ t {ϕi } γ i H i* + γ i( ) H i* + γ i( ) H i* + ) (107) 66 BehaviourofElectromagneticWavesinDifferentMediaandStructures {Δ... matrix, and using Eq (4), we have T T D* = Λ E * + GH * (41) B* = GE * + Π H * (42) Rewriting the above in indicial notation, we get Di* = λi Ei* + gij H * j i = 1,2 ,3 j = 1,2 ,3 ( 43) Bi* = γ i H i* + gij E* j i = 1,2 ,3 j = 1,2 ,3 (44) Eqs ( 43) and (44) are just the modal constitutive equations for bi-anisotropic media 5.2 Eigen equations ofelectromagneticwavesin bi-anisotropic media Substituting Eqs ( 43) . Eqs. (33 ) and (34 ) can be uncoupled in the form of scalar ** 2* 0 0 1,2 ,3+ ∇ = = εγ ii iti HHi (35 ) ** 2* 0 0 1,2 ,3+ ∇= = εγ ii iti EEi (36 ) Eqs. (35 ) and (36 ) are the modal equations of electromagnetic. Behaviour of Electromagnetic Waves in Different Media and Structures 64 () * 11 23 3 3 =++HHHH (89) () *222 1 1 3 =−∂+∂+∂ xyz (90) So, the equation of electromagnetic wave in lossy. are being rotated the produced stop band is widening and in the case of 90 ◦ rotation it embraces both stop bands as can be seen in Fig. 24(c). 3. 3.4 Tunneling effect An interesting effect of wave