Evaluation of the Green's function for the mixed potential integral equation (MPIE) method in the time domain for layered media

In the use of the time-domain integral equation (TDIE) method for the analysis of layered media, it is important to have the time-domain layered medium Green's function computed for many source-to-field distances /spl rho/ and time instants t. In this paper, a numerical method is used that computes the mixed potential Green's functions G/sub v/(/spl rho/,t) and G/sub A/(/spl rho/,t) for a multilayered medium for many /spl rho/'s and t's simultaneously. The method is applicable to multilayered media and for lossless or lossy dispersive media. Salient features of the method are: 1) the use of complex /spl omega/ so that the surface wave poles are lifted off the real k/sub /spl rho// axis such that pole extractions are not required; 2) the use of half-space extraction so that the integrand for the Sommerfeld integral decays exponentially along the k/sub /spl rho// axis to obtain fast convergence of the integral; and 3) the use of the fast Hankel transform so that the Green's function is calculated for many values of /spl rho/ simultaneously. For a four-layer medium, we illustrate the numerical results by a three-dimensional plot of /spl rho/G/sub v/(/spl rho/,t) versus /spl rho/ and t and demonstrate the space-time evolution of these Green's functions. For a maximum frequency range of 8 GHz, the method requires only a few CPU minutes to compute a table of 100 (points in /spl rho/) /spl times/ 168 (points in t) uniformly spaced values of G/sub v/(/spl rho/,t) on an 867-MHz Pentium PC.

[1]  K. Michalski,et al.  Multilayered media Green's functions in integral equation formulations , 1997 .

[2]  E. Michielssen,et al.  Analysis of transient electromagnetic scattering from closed surfaces using a combined field integral equation , 2000 .

[3]  Leung Tsang,et al.  Surface electric fields and impedance matrix elements of stratified media , 2000 .

[4]  Ramachandra Achar,et al.  Simulation of high-speed interconnects , 2001, Proc. IEEE.

[5]  L. Tsang,et al.  Closed-form spatial electric field Green's functions of microstrip structures using the fast Hankel transform and the matrix pencil method , 2000 .

[6]  Jen-Tsai Kuo,et al.  Fast full-wave analysis of planar microstrip circuit elements in stratified media , 1998 .

[7]  J. Mosig,et al.  An Overview of Integral Equation Techniques Applied to Microstrip Antenna Design , 1985, 1985 15th European Microwave Conference.

[9]  J. Mosig.  Integral equation Technique , 1989 .

[11]  Andrew J. Poggio,et al.  CHAPTER 4 – Integral Equation Solutions of Three-dimensional Scattering Problems , 1973 .

[12]  L. Tsang,et al.  A sparse-matrix/canonical grid method for analyzing densely packed interconnects , 2000 .

[13]  Chi Hou Chan,et al.  Complex image method for sources in bounded regions of multilayer structures , 1994 .

[14]  G. E. Howard,et al.  A closed-form spatial Green's function for the thick microstrip substrate , 1991 .

[15]  James R. Wait,et al.  A pulsed dipole in the earth , 1976 .

[16]  Leung Tsang,et al.  Numerical evaluation of the transient acoustic waveform due to a point source in a fluid-filled borehole , 1979 .

[17]  Niels B. Christensen,et al.  Optimized fast Hankel transform filters , 1990 .

[18]  Leung Tsang,et al.  Surface electric fields and spatial derivatives of Green's function of layered media based on half-space extraction , 2000 .