volume | PIER Journals

2006-10-26

A Rapidly-Convergent, Mixed-Partial Derivative Boundary Condition Green's Function for an Anisotropic Half-Space: Perfect Conductor Case

John Jarem

Dr. John Jarem

Electrical and Computer Engineering Department

University of Alabama in Huntsville

USA

Homepage

Progress In Electromagnetics Research, Vol. 67, 39-112, 2007

doi:10.2528/PIER06072804

Abstract

The problem of determining the Green's function of an electric line source located in a permeable, anisotropic (Ëœμ_xx, Ëœμ_xy, Ëœμ_yx, Ëœμ_yy and Ëœε_zz nonzero interacting parameters) half-space above a Perfect Magnetic Conductor (PMC) ground plane (called the TM^z case herein) for the case where image theory cannot be applied to find the Green's function of the PMC ground plane system has been studied. Monzon [2, 3] studied the Green's function TE^z problem dual to the present one for two cases; (1) when the system was unbounded, anisotropic space where Ëœε_xx, Ëœε_xy, Ëœε_yx, Ëœε_yy and Ëœμ_zz were the nonzero interacting parameters; and (2) when the scattering system was an anisotropic half-space located above a Perfect Electric Conductor (PEC) ground plane and where Ëœε_xx, Ëœε_yy and Ëœμ_zz were the nonzero interacting parameters and Ëœε_xy = Ëœε_yx = 0. Monzon [2] referred to the latter ground plane case as the case where "usual image" theory could be used to find the Green's function of the system. The Green's function for the TM^z-PMC case studied herein was derived by introducing and using a novel, linear coordinate transformation, namely Ëœx = (Ïƒ_P /Ï„)Ëœx + Ïƒ_MËœy, Ëœy = Ëœy, (Eqs. (6c,d), herein). This transformation, a modification of that used by [2, 3], reduced Maxwell's anisotropic equations of the system to a non-homogeneous, Helmholtz wave equation from which the Green's function, G, meeting boundary conditions, could be determined. The coordinate transformation introduced was useful for the present PMC ground plane problem because it left the position of the PMC ground plane and all lines parallel to it, unchanged in position from the original coordinate system, thus facilitating imposition of EM boundary conditions at the PMC ground plane. In transformed (or primed) coordinates, for the TMz-PMC and TE^z-PEC ground plane problems, respectively, the boundary conditions for the Green's functions G_TM â‰¡ E_z and G_TE â‰¡ H_z were shown to be [Î±*(âˆ‚G/âˆ‚Ëœx) + (âˆ‚G/âˆ‚Ëœy) |_Ëœy=0 = 0, where G = G_TM or G_TE, Î± = Î±_TM or Î±_TE, where Î±_TM and Î±_TE are complex constants (dually related to each other by Ëœμ â†” Ëœε), and where Ëœy' = 0 is the position of the ground plane. An interesting result of the analysis was that the constant Î±TE (as far as the author knows) coincidently turned out to be the same as the first of two important constants, namely S₁(dl_v/dl) and S2(dl_v/dl), which were used by Monzon [2, 3] to formulate integral equations (based on Green's second theorem) from which EM scattering from anisotropic objects could be studied. Spatial Fourier transform (k-space) techniques were used to determine the Green's function of the Helmholtz wave equation expressed in transformed coordinates which satisfied the mixed-partial derivative boundary condition of the system. The Green's function G was expressed as a sum of a "free space" Green's function gf (proportional to a Hankel function H⁽²⁾₀ and assumed excited by the line source in unbounded space) and a homogeneous Green's function g whose spectral amplitude was chosen such that, when g was added to g_f , the sum G = g_f + g, satisfied boundary conditions. The k-space, Sommerfeld integrals making up g turned out to converge slowly, and so, using contour integration in the complex plane (Appendix B), g was expressed in a rapidly convergent form, and thus one leading to its practical numerical evaluation. Extensive numerical testing of how well the Green's functions G = g_f + g satisfied boundary conditions and how well the homogeneous Green's function g satisfied Green's second theorem was performed. Excellent verification of the numerical and analytical procedures were found and displayed in the error tables of Tables 1-4 of the paper. Plots illustrating the novel coordinate transformation introduced were presented and several numerical plots of the Green's function developed herein were presented. The application of the present work to find the Green's function for the case when a negative-index, anisotropic metamaterial is adjacent to a conducting ground plane is discussed. Application of the Green's function theory developed herein to multi-layer anisotropic systems is discussed.

Citation

Copy Export

John Jarem, "A Rapidly-Convergent, Mixed-Partial Derivative Boundary Condition Green's Function for an Anisotropic Half-Space: Perfect Conductor Case," Progress In Electromagnetics Research, Vol. 67, 39-112, 2007.
doi:10.2528/PIER06072804

References

1. Monzon, J. C. and N. J. Damaskos, "Two-dimensional scattering by a homogeneous anisotrpoic rod," IEEE Transactions on Antennas and Propagation, Vol. 34, No. 10, 1243-1249, 1986.
doi:10.1109/TAP.1986.1143739

2. Monzon, J. C., "On a surface integral representation for homogeneous anisotropic regions: two-dimensional case," IEEE Transactions on Antennas and Propagation, Vol. 36, No. 10, 1401-1406, 1988.
doi:10.1109/8.8627

3. Monzon, J. C., "The equivalence principle for two-dimensional anisotropies," IEEE Transactions on Antennas and Propagation, Vol. 39, No. 12, 1781-1784, 1991.
doi:10.1109/8.121601

4. Shul'ga, S. N., "Scattering of electromagnetic wave by an inhomogeneous anisotropic inclusion immersed in anisotropic halfspace," Journal of Communications Technology and Electronics, Vol. 43, No. 10, 1086-1090, 1998.

5. Shul'ga, S. N., "Two-dimensional problem of scattering of a wave beam by the anisotropic half-space with an anisotropic inclusion," Optics and Spectroscopy, Vol. 87, No. 3, 465-471, 1999.

6. Zhuk, N. P., S. N. Shul'ga, and A. G. Yarovoi, "Two-dimensional scattering of electromagnetic waves from a permeable inclusion in an anisotropic medium," Technical Physics, Vol. 43, No. 1, 75-79, 1998.
doi:10.1134/1.1258940

7. Zhuck, N. P. and A. G. Yarovoy, "Two-dimensional scattering from a inhomogeneous dielectric cylinder embedded in a stratified medium: case of TM polarization," IEEE Transactions on Antennas and Propagation, Vol. 42, No. 1, 16-21, 1994.
doi:10.1109/8.272296

8. Goryushko, D. N.A. V. Malyuskin, S. N. Shulga, and A. A. Shmat'ko, "Wavetrains scattering by anisotropic composite absorbing layer with metal-dielectric substrate,'' 2002 12th Int. Conference Microwave & Telecommunications Technology," (CriMiCo'2002), 9-13, 2002.

9. Malyuskin, A. V.D. N. Goryushko, A. A. Shmat'ko, and S. N. Shulga, "Scattering of a wave beam by an inhomogeneous anisotropic chiral layer," MMET'02Pr oceedings, 566-568, 2002.

10. Malyuskin, A. V.M. P. Perepechai, and S. N. Shulga, "Effective electromagnetic parameters of strongly fluctuating statistically layered bianisotropic medium," MMET'2000 Proceedings, 361-363, 2000.

11. Simsek, E., J. Liu, and Q. H. Liu, "A spectral integral method (SIM) for layered media," IEEE Transactions on Antennas and Propagation, Vol. 54, No. 6, 1742-1749, 2006.
doi:10.1109/TAP.2006.875500

12. Monzon, J. C., "On the application of the Sommerfeld representation in a two-dimensional rotationally invariant anisotropic medium," IEEE Transactions on Antennas and Propagation, Vol. 38, No. 7, 1028-1034, 1990.
doi:10.1109/8.55600

13. Monzon, J. C., "Two-dimensional Green's functions for a rotationally invariant anisotropic medium," IEEE Transactions on Antennas and Propagation, Vol. 38, No. 5, 616-624, 1990.
doi:10.1109/8.53489

14. Monzon, J. C., "Two-dimensional surface integral representations for a rotationally invariant anisotropic medium: applications to scattering," IEEE Transactions on Antennas and Propagation, Vol. 39, No. 1, 108-111, 1991.
doi:10.1109/8.64444

15. Monzon, J. C., "Three-dimensional field expansion in the most general rotationally symmetric anisotropic material: application to scattering by a sphere," IEEE Transactions on Antennas and Propagation, Vol. 37, No. 6, 728-735, 1989.
doi:10.1109/8.29359

16. Monzon, J. C., "Radiation and scattering in homogeneous general biisotropic regions," IEEE Transactions on Antennas and Propagation, Vol. 38, No. 2, 227-235, 1990.
doi:10.1109/8.45125

17. Ren, W., "Contributions to the electromagnetic wave theory of bounded homogeneous anisotropic media," Physical Review E, Vol. 47, No. 1, 664-673, 1993.
doi:10.1103/PhysRevE.47.664

18. Lindell, I. V., A. H. Sihvola, S. A. Tretyakov, and A. J. Viitanen, Electromagnetic Waves in Chiral and Bi-isotropic Media, Artech House, 1994.

19. Jarem, J. M. and P. P. Banerjee, Computational Methods for Electromagnetic and Optical Systems, Marcel Dekker Inc., 2000.

20. Shulga, S. N. and O. V. Bagatskaya, "Analysis of a waveguide Tjunction with a 2D scatterer in the interaction region via Green's theorem approach," International Conference on Antenna Theory and Techniques, 785-788, 2003.

21. Felson, L. B. and N. Marcuvitz, Radiation and Scattering of Waves, Prentice-Hall Inc., 1973.

22. Pendr, J. B. and D. R. Smith, "The quest for the superlens," Scientific American, Vol. 295, No. 1, 60-67, 2006.

23. Itoh, T. and A. A. Oliner (eds.), "Special Issue on Metamaterial, IETMAB," IEEE Trans. on Microwave Theory and Techniques, Vol. 53, No. 4, 1413-1417, 2005.
doi:10.1109/TMTT.2005.845126

24. Engheta, N. and R. W. Ziolokwski, "A positive future for doublenegative metamaterials," IEEE Trans. on Microwave Theory and Techniques, Vol. 53, No. 4, 1535-1556, 2005.
doi:10.1109/TMTT.2005.845188

25. Abramowitz, M. and I. A. Stegun, Handbook of Mathematical Functions, Dover Publications Inc., 1972.

26. Ishimaru, A., Electromagnetic Wave Propagation, Radiation, and Scattering, Prentice Hall, 1991.

27. Visual Numerics, Inc., "Bessel functions," (Bessel function subroutines BSJS/DBSJS and BSYS/DBSYS), No. 7, 103-106, 2006.

28. Jarem, J. M., "Rigorous coupled wave thoery of anisotropic, azimuthally-inhomogenoues cylindrical systems," Progress In Electromagnetics Research, Vol. 19, 109-127, 1998.
doi:10.2528/PIER97103100