PIER | |

Progress In Electromagnetics Research | ISSN: 1070-4698, E-ISSN: 1559-8985 |

Home > Vol. 67 > pp. 39-112
## A RAPIDLY-CONVERGENT, MIXED-PARTIAL DERIVATIVE BOUNDARY CONDITION GREEN'S FUNCTION FOR AN ANISOTROPIC HALF-SPACE: PERFECT CONDUCTOR CASEBy J. M. Jarem
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 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 ^{z}TE problem dual to the present one for two cases; (1) when the system was unbounded, anisotropic space where Ëœε^{z}_{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-PMC case studied herein was derived by introducing and using a novel, linear coordinate transformation, namely ^{z}Ëœ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 were shown to be _{TM} â‰¡ E_{z} and G_{TE} â‰¡ H_{z}[Î±*(âˆ‚G/âˆ‚Ëœx) + (âˆ‚G/âˆ‚Ëœy) |, where _{Ëœy=0} = 0G = G or _{TM}G, _{TE}Î± = Î± or _{TM}Î±, where _{TE}Î± and _{TM}Î± are complex constants (dually related to each other by Ëœμ â†” Ëœε), and where _{TE}Ëœ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 and _{1}(dl_{v}/dl)S2(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 (_{v}/dl)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 ^{(2)}_{0}g whose spectral amplitude was chosen such that, when g was added to g , the sum _{f}G = g, satisfied boundary conditions. The _{f} + gk-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 satisfied boundary conditions and how well the homogeneous Green's function _{f} + gg 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.
2. Monzon, J. C., "On a surface integral representation for homogeneous anisotropic regions: two-dimensional case," 3. Monzon, J. C., "The equivalence principle for two-dimensional anisotropies," 4. Shul'ga, S. N., "Scattering of electromagnetic wave by an inhomogeneous anisotropic inclusion immersed in anisotropic halfspace," 5. Shul'ga, S. N., "Two-dimensional problem of scattering of a wave beam by the anisotropic half-space with an anisotropic inclusion," 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," 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," 8. Goryushko, D. N., A. V. Malyuskin, S. N. Shulga, and A. A. Shmat'ko, 9. Malyuskin, A. V., D. N. Goryushko, A. A. Shmat'ko, and S. N. Shulga, 10. Malyuskin, A. V., M. P. Perepechai, and S. N. Shulga, 11. Simsek, E., J. Liu, and Q. H. Liu, "A spectral integral method (SIM) for layered media," 12. Monzon, J. C., "On the application of the Sommerfeld representation in a two-dimensional rotationally invariant anisotropic medium," 13. Monzon, J. C., "Two-dimensional Green's functions for a rotationally invariant anisotropic medium," 14. Monzon, J. C., "Two-dimensional surface integral representations for a rotationally invariant anisotropic medium: applications to scattering," 15. Monzon, J. C., "Three-dimensional field expansion in the most general rotationally symmetric anisotropic material: application to scattering by a sphere," 16. Monzon, J. C., "Radiation and scattering in homogeneous general biisotropic regions," 17. Ren, W., "Contributions to the electromagnetic wave theory of bounded homogeneous anisotropic media," 18. Lindell, I. V., A. H. Sihvola, S. A. Tretyakov, and A. J. Viitanen, 19. Jarem, J. M. and P. P. Banerjee, 20. Shulga, S. N. and O. V. Bagatskaya, 21. Felson, L. B. and N. Marcuvitz, 22. Pendr, J. B. and D. R. Smith, "The quest for the superlens," 23. Itoh, T. and A. A. Oliner (eds.), "Special Issue on Metamaterial, IETMAB," 24. Engheta, N. and R. W. Ziolokwski, "A positive future for doublenegative metamaterials," 25. Abramowitz, M. and I. A. Stegun, 26. Ishimaru, A., 27. Visual Numerics, Inc., "Bessel functions," 28. Jarem, J. M., "Rigorous coupled wave thoery of anisotropic, azimuthally-inhomogenoues cylindrical systems," |

© Copyright 2014 EMW Publishing. All Rights Reserved