Automatic Incorporation of Surface Wave Poles in Discrete Complex Image Method

Abstract

Abstract-Discrete complex image method is introduced to get a closed-form dyadic Green's function by a sum of spherical waves. However, the simulation result by the traditional discrete complex image method is only valid in near-field for several wavelengths. In this paper, we analyze the form of spectral domain dyadic Green's function in the whole k ρ plane and the variety of valid range of simulation results by different sampling paths in two-level discrete complex image method. Consequently, for dyadic Green's function, surface wave pole contribution both in spectral domain and spatial domain is clarified. We introduce the automatic incorporation of surface wave poles in discrete complex image method without extracting surface wave poles. The contribution of surface wave poles in spectral domain and spatial domain dyadic Green's function is further confirmed in the new method. Besides, this method can represent dyadic Green's function by spherical waves in the layer where the source and field points are. So it satisfies the splitting requirement and consequently reduces the computational complexity dramatically especially for objects with large scale inẑ direction.

Modal Analysis of Surface Plasmon Resonance Sensor Coupled to Periodic Array of Core-Shell Metallic Nanoparticles

“…When solving problems involving integrals like Eq. (10), several recent approaches have been proposed [36,37], all consist of the evaluation of spectral functions using the Sommerfeld Identity with variants of the discrete complex image method (DCIM) as an acceleration tool.…”

Radiation of a Horizontal Dipole in the Presence of a Three-Layered Region and Microstrip Antenna

“…5 and 6, the total field, both the trapped surface waves of electric type (TM) and magnetic type (TE), and the lateral wave are calculated from (27) for E r 0Θ (r 0 , Θ, Φ) and (28) for E r 0Φ (r 0 , Θ, Φ), respectively. Graphs of the amplitudes of the components E r 0Θ (r 0 , Θ, Φ) and E r 0Φ (r 0 , Θ, Φ) for a unit horizontal dipole vs. the angular degrees of Θ are computed and shown in Figs.…”

Numerical simulation of a grounding grid buried in horizontal multilayered earth based on ‘T’ typical basic elements in the frequency domain

“…In addition, extracting the surface wave components numerically is complicated [14,15], and the surface wave terms cannot be obtained for analytical calculations of mutual impedance and the induction coefficients. Thus, we employ an improved two-level DCIM without extracting the surface wave poles to reduce the approximation complexity, and to extend the valid range for planar structures that contain both lossless and lossy media [16]. The improved two-level DCIM is used to study the multilayered earth model and thus the quasi-static term of Green's function can be extracted easily.…”