# Evaluation of the Trapped Surface Wave of a Vertical Electric Dipole Based on Undetermined Coefficient Method in the Presence ofN-Layered Region: A Graphical Approach

Authors: Hong Lei Xu; Ting Ting Gu; Yong Zhu; Xiao Wei; Liang Sheng Li; Hong Cheng Yin
Journal: International Journal of Antennas and Propagation
 (2019)

---

## Abstract

In previous studies, the trapped surface wave, which is defined by the residue sums, has been addressed in the evaluation of the Sommerfeld integrals describing electromagnetic field of a vertical dipole in the presence of three-layered or four-layered region. But unfortunately, the existing computational scheme cannot provide analytical solution of the field in the presence of theN-layered region when N > 4. The scope of this paper is to overcome the limitations in root finding algorithm implied by the previous approach and provide solution of poles in stratified media. A set of pole equations following with explicit expressions are derived based on the undetermined coefficient method, which enable a graphical approach to obtain initial values of real roots. Accordingly, the generated trapped surface wave components are computed when both the observation point and the electric dipole source are on or near the surface of a dielectric-coated conductor. Validity, efficiency, and accuracy of the proposed method are illustrated by numerical examples.


---

## Body

## 1. Introduction

It is known that the electromagnetic fields radiated by a vertical or horizontal electric dipole in stratified media are interesting and practically important in many cases [1–3]. In the past decades, this problem has led to many published papers and available achievements on analytical solutions [4–21], especially the three-layered or four-layered cases.The general integral representation of the electromagnetic field due to a dipole source has been addressed by King et al. [6] in the presence of the N-layered region, which was developed from the original expressions firstly formulated by Sommerfeld in 1909 [7] in the presence of half-spaces. Efforts have been made in a series of works [6–8] by investigators to derive analytical expressions for the Sommerfeld integrals, which leads to a better understanding of describing electromagnetic radiation from the dipole source than numerical solution, as well as time savings with respect to conventional techniques used to evaluate the integrals. In Chapter 15 of the monograph [6], the propagation of the electromagnetic pulses radiated by a horizontal electric dipole with delta-function excitation in the presence of the three-layered region was processed analytically, which demonstrated that the total field on or near the air-dielectric boundary is determined primarily by lateral wave, where the amplitude of the field along the boundary is 1/ρ2. Unfortunately, the integrals cannot be evaluated by means of the mentioned analytical procedure for electromagnetic field in the presence of the N-layered region.In the comments by Wait and Mahoud et al. in 1998 [9, 10], and studies by other pioneers, particularly including Collin [11, 12] and Zhang and Pan [13], the three-layered structure was reconsidered by the use of asymptotic methods, contour integration, and branch cuts, where it is pointed out that the trapped surface wave, which is determined by residue sums of the poles, can be excited efficiently by a dipole source with the amplitude of the field 1/ρ1/2 and should not be neglected over a dielectric-coated conductor [13]. To extend the study, the Sommerfeld integrals have been evaluated in the four-layered case by a similar method [14, 15], correspondingly. The details are summarized in a recent book by Li [16].These new developments on the analytical results for the electromagnetic field in three-layered and four-layered structures [9–16], where it was proved that the trapped surface wave defined by the residue sums as the dominant wave propagates along the surface of air-dielectric boundary at long propagation distance [13], aroused interest in the study on properties of the trapped surface wave in the N-layered structure. However, few literatures relate to a computational scheme to solve the poles, for which the surface impedance at the air-dielectric boundary is in expression of a recursive form due to multireflections [6, 17]. Specially, the discrete pole root λj∗ is hard to solve by analytical and numerical solutions due to multivalue properties of the recursive equation over four-layered media.In a recent study by Cross, the solution of poles relying on numerical root finding algorithms is addressed in the four-layered structure to approximate the scenario of a leaking water-pipe, buried in the shallow subsurface [18], where it is suggested that the analytical solution of poles requires efficiency and enhancement of accuracies in stratified media.Following the research line, we are attempting to derive the pole equation to release the difficulties in root finding in theN-layered structure. In the analysis, a set of pole equations with explicit expressions is developed for a computational scheme, so as to make it possible analytical evaluation of the trapped surface wave based on the undetermined coefficient method. The solution of obtained equations in three-layered to six-layered structures is carried out, respectively, as illustrative examples in lossless case by a graphical approach. The obtained equations offer advantages in terms of time savings with respect to standard pole equation in root finding procedures. In addition, computation and discussion are carried out by investigating the full-wave analytical expressions, as well as their interfering behaviour, which guarantees correctness by evaluating the residue sums' contributions to the fields, but also allows us to gain useful insight into the physics of the problem.


## 2. Formulation of the Problem

The 3D geometry under consideration and its 2D cylindrical coordinate system are shown in Figures1 and 2, respectively, where the vertical electric dipole in the ϕ=0 direction is located at 0,0,d. Region z>0 is the upper half-space occupied by air, and the lower half-space is composed of a success of n+1 horizontal layers, each with arbitrary thickness ljj=1,2,…,n and arbitrary wave number kj, in which(1)k0=ωμ0ε0,kj=ωμ0ε0εrj+iσjω,j=1,2,…,n+1.Figure 1
A multilayered model representing the simplest approximation of wave propagation.Figure 2
Geometrical profile of the horizontalN-layered model.In the monograph by King et al. [6], the general integral formulas are addressed for the electromagnetic field due to a vertical electric dipole in the presence of N-layered media. The integrated formulas in Region 0 (Air) are written in forms of(2)
B
0
ϕ
=
i
μ
0
4
π
∫
0
∞
e
i
γ
0
z
−
d
+
e
i
γ
0
z
+
d
−
Q
n
+
1
e
i
γ
0
z
+
d
γ
0
−
1
J
1
λ
ρ
λ
2
d
λ
,
(3)
E
0
ρ
=
i
ω
μ
0
4
π
k
0
2
∫
0
∞
±
e
i
γ
0
z
−
d
+
e
i
γ
0
z
+
d
−
Q
n
+
1
e
i
γ
0
z
+
d
J
1
λ
ρ
λ
2
d
λ
,
z
≥
d
,
0
<
z
<
d
,
(4)
E
0
z
=
−
ω
μ
0
4
π
k
0
2
∫
0
∞
e
i
γ
0
z
−
d
+
e
i
γ
0
z
+
d
−
Q
n
+
1
e
i
γ
0
z
+
d
γ
0
−
1
J
0
λ
ρ
λ
3
d
λ
.In aforementioned formulas,n represents for the number of intermediate dielectric layers and the reflection coefficient Qn is expressed by(5)Qn=−γ0−k02/ωμ0Zs00γ0+k02/ωμ0Zs00,in which Zs00 represents the surface impedance. If the procedure is continued to the boundary at z′=z=0 between Region 0 and Region 1, the desired surface impedance is obtained. It is(6)Zs0′0=ωμ0γ0k02tanhtanh−1γ1k02γ0k12tanh−iγ1l1+tanh−1γ2k12γ1k22tanh−iγ2l2+tanh−1γ3k22γ2k32tanh−iγ3l3+⋯+tanh−1γnkn−12γn−1kn2tanh−iγnln+tanh−1γn+1kn2γnkn+12…,(7)γi=ki2−λ2;i=0,1,2,…,n+1.The first and second terms in (2)–(4) stand for the direct wave and the ideal reflected wave, respectively, and have been previously solved in [6] by King et al. Considering that γ0 and γ1 are functions of λ, which is made of the relations between Bessel and Hankel function:(8)Jnλρ=12Hn1λρ+Hn2λρ,Hn1−λρ=Hn2λρ−1n+1.Thus, the third terms in (2)–(4) are rewritten as follows:(9)B0ϕ3=−ik024πω∫−∞∞Zs00H11λρeiγ0z+dλ2γ0+k02/ωμ0Zs0′0γ0dλ,(10)E0ρ3=−i4π∫−∞∞Zs00H11λρeiγ0z+dλ2γ0+k02/ωμ0Zs0′0dλ,(11)E0z3=14π∫−∞∞Zs00H01λρeiγ0z+dλ3γ0+k02/ωμ0Zs0′0γ0dλ.By similar treatment of asymptotic methods, contour integration and branch cuts applied in three-layered and four-layered structures, which is addressed in detail by Li [16], the terms from (9)–(11) can be considered as the combination of the branch cut integrals and contribution of the residues of poles, which are defined by the trapped-surface-wave-term and the lateral-wave-term, respectively. It is necessary to shift the contours around all the branches at λ1=k0, λ2=k1, …, λj=kn+1, respectively. Since the approximation of each branch cut integral refers to a lot of mathematical derivations, which is beyond the main scope of this paper, the lateral-wave-term is not addressed in illustrative examples for simplicity. In the next, attention is made on evaluation of the trapped-surface-wave-term by inviting analytical techniques to evaluate the discrete pole residues.


## 3. Evaluation of the Trapped Surface Wave in the Presence of theN-Layered Region

### 3.1. Residue of the Poles for the Trapped-Surface-Wave-Term

Following the integration procedure addressed in [16], the trapped-surface-wave-term is defined by the sum of pole residues. In the N-layered structure, the terms of trapped surface wave due to a vertical electric dipole are written as(12)B0ϕSur=k022ω∑jZs0′λj∗H11λj∗ρeiγ0λj∗z+dλj∗2qN′λj∗γ0λj∗,E0ρSur=12∑jZs0′λj∗H11λj∗ρeiγ0λj∗z+dλj∗2qN′λj∗,E0zSur=i2∑jZs0′λj∗H01λj∗ρeiγ0λj∗z+dλj∗3qN′λj∗γ0λj∗,in which the function qN′ is expressed by(13)qN′=∂∂λγ0+k02ωμ0Zs0′0n=N−2,and λj∗j=1,2,… are the discrete roots of the pole equation of the electromagnetic field in the presence of the N-layered region, which is in a variant of the surface impedance function Zs0′0 at the air-to-dielectric boundary where n=N−2 represents the number of intermediate dielectric layers of the N-layered structure, written as(14)qN=γ0+k02ωμ0Zs0′0n=N−2=0.Considering the pole equation by (14) with substitution of (6) in the form of a recursive expression, it is necessary to derive the explicit pole equation first. Suppose (14) can be rewritten with respect to the function GNn, so that(15)qN=γ0+γ1k02k12GNnλ=0,in which GNn is defined by(16)GNnλ=k12ωμ0γ1Zs0′0.The superscriptn of function GNn represents the recursive order of the function that is in a variant of λ. In the meanwhile, (15) can be expanded if the order n≥1 by taking into account of (6) with the function GNn−1. Specifically, it is expressed by a set of sequential equations, as follows:(17)qN=2=γ0+γ1k02k12G20λ,qN=3=γ0+γ1k02k12G31λ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22G30λ,qN=4=γ0+γ1k02k12G42λ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22G41λ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22⋅tanh−iγ2l2+tanh−1γ3k22γ2k32G40λ,…qN=γ0+γ1k02k12GNnλ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22GNn−1λ.It is seen from equation (17) that the recursive order n of function GNn is reduced if the following identity is satisfied:(18)GNj+1λ=tanh−iγmlm+tanh−1γm+1km2γmkm+12GNjλ,m=N−j−1.Combining equations from (17), it is inferred that GNjλj=0≡1 is applied in each equation when n=N−2. Consequently, the equations are adapted by exploiting both (18) and GNjλj=0≡1, as follows:(19)N=3:tanhiγ1l1+G30λ=γ2k12γ1k22+G30λγ2k12γ1k22tanhiγ1l1,whenG30λ=1,(20)N=4:tanhiγ2l2+G41λ=γ2k12γ1k22+G41λγ2k12γ1k22tanhiγ2l2,whenG40λ=1,(21)N=5:tanhiγ3l3+G52λ=γ2k12γ1k22+G52λγ2k12γ1k22tanhiγ3l3,whenG50λ=1,…In (19), the explicit expression of pole equation in the three-layered structure is easy to obtain by G30=1. Similarly, for general case with n > 3, such as in the four-layered and five-layered structures by (20) and (21), respectively, where the function GNn−1 with respect to the variable λ has been applied, the explicit equation is defined by suppressing GN0=1. Through mathematical derivations, the variant coefficients are derived accordingly in expression of recursive form, as follows: (22)N≥4:G30λ=−γ0k12γ1k02,N≥5:G41λ=tanhiγ1l1+G30λγ2k12/γ1k22+G30λγ2k12/γ1k22tanhiγ1l1,N≥6:G52λ=tanhiγ2l2+G41λγ3k22/γ2k32+G41λγ3k22/γ2k32tanhiγ2l2.…Therefore, it is concluded that the general expression of pole equation for electromagnetic field in the presence ofN-layered region is derived from (14), in expression of(23)tanhiγnln+GNn−1λ=γ2k12γ1k22+GNn−1λγ2k12γ1k22tanhiγnln,whenGN0=1,where(24)GNjλ=tanhiγjlj+GNj−1λγj+1kj2/γjkj+12+GNj−1λγj+1kj2/γjkj+12tanhiγjlj.


### 3.2. Solution of Poles: A Graphical Approach

In order to evaluate the pole residues analytically, the set of equations from (19)–(21) are developed into explicit expressions with substitution of equations from (22). For convenience, the pole equation in the presence of three-layered region obtained from (19) is rewritten as(25)N=3:tanhiγ1l1=γ2k12γ1k22+γ0k12γ1k02−γ0k12γ1k02γ2k12γ1k22tanhiγ1l1.Analogously, the pole equation is derived for electromagnetic field of a vertical electric dipole in the presence ofN-layered region by substituting (23) with (24), iteratively. Specifically, the pole equations for the fields in the four-layered to six-layered structures can be written as follows:(26)N=4:tanhiγ2l2=γ3k22γ2k32+γ1k22γ2k12⋅γ1k02γ0k12×tanhiγ1l1−A1λtanhiγ1l1−B1λ−γ3k22γ2k32γ1k22γ2k12⋅γ1k02γ0k12×tanhiγ1l1−A1λtanhiγ1l1−B1λ⋅tanhiγ2l2,(27)N=5:tanhiγ3l3=γ4k32γ3k42+γ2k32γ3k22⋅γ0k12γ1k02×tanhiγ2l2−A2λtanhiγ2l2−B2λ−γ4k32γ3k42γ2k32γ3k22⋅γ0k12γ1k02×tanhiγ2l2−A2λtanhiγ2l2−B2λ⋅tanhiγ3l3,(28)N=6:tanhiγ4l4=γ5k42γ4k52+γ3k42γ4k32⋅γ0k12γ1k02×tanhiγ3l3−A3λtanhiγ3l3−B3λ−γ5k42γ4k52γ3k42γ4k32⋅γ0k12γ1k02×tanhiγ3l3−A3λtanhiγ3l3−B3λ⋅tanhiγ4l4.…By a few mathematical efforts, the functionsAjλ and Bjλ from (26)–(28) are derived through iterative substitutions, which are listed in Appendix, correspondingly. It is noted that the pole equation can be expressed in analogous form for the electromagnetic field of a vertical dipole in the presence of N-layered media. Accordingly, substituting the equations from (26)–(28) with functions Ajλ, the explicit expression of equations are obtained. By transpositions and rearrangements, the expression of pole equation in the presence of N-layered region can be described as(29)1+γn+1kn2γnkn+12⋅γn−1kn2γnkn−12⋅∑m=1n−1∏j=1mCn∗j,mtanhiγjlj∑k=1n−1∏j=1mCn∗m−j,n−mtanhiγjljtanhiγnln=γn+1kn2γnkn+12+γn−1kn2γnkn−12⋅1−∑m=1n−1∏j=1mCn∗j,mtanhiγjlj∑m=1n−1∏j=1mCn∗m−j,n−mtanhiγjlj,in which the coefficients Cn∗j,m, j = 1, …, m; m = 1, …, n are a set of coefficients defined by positive or negative ratio of γjkm and γmkj, which are determined from (26)–(28). If the bottom half-space is considered as perfectly conducting layer as depicted in Figure 3, the expression of (29) are simplified in condition of kn+1⟶∞, which reduces to(30)tanhiγnln=γn−1kn2γnkn−12⋅1−∑m=1n−1∏j=1mCn∗j,mtanhiγjlj1−∑m=1n−1∏j=1mCn∗m−j,n−mtanhiγjlj.Figure 3
TheN-layered model with the bottom half-space as conducting.It is seen from (29) and (30) that the obtained equations for the fields are expressed in form of addition, subtraction, and multiplication of functions tanhiγjlj with coefficients Cn∗, which offer advantages in terms of time savings with respect to standard numerical root finding procedures by (14). Specifically, the obtained pole equations for the electromagnetic field over a conductor coated by two-layered to four-layered dielectrics are in terms of the following equations:(31)N=4:tanhiγ2l2=γ1k22γ2k12⋅γ1k02γ0k12tanhiγ1l1−A1λtanhiγ1l1−B1λ,(32)N=5:tanhiγ3l3=γ2k32γ3k22⋅γ0k12γ1k02tanhiγ2l2−A2λtanhiγ2l2−B2λ,(33)N=6:tanhiγ4l4=γ3k42γ4k32⋅γ0k12γ1k02tanhiγ3l3−A3λtanhiγ3l3−B3λ.…It is convenient to validate (31) in derivation followed, which is exactly the same with (4.24) presented in [16] for the electromagnetic field of a vertical electric dipole on the surface of a conductor coated by two dielectric layers.(34)qN=4Per=γ0−iγ1k02k12tanγ1l1+tanh−1γ2k12γ1k22tanγ2l2=γ0−iγ1k02k12tanγ1l1+γ2k12γ1k22tanγ2l2×1−γ2k12γ1k22tanγ1l1tanγ2l2−1=0.Through the aforementioned analysis, the computational scheme to evaluate the residue sums for the trapped surface wave is proposed by exploiting the derived equations from (25)–(28). A graphical approach is applied in the present study to obtain initial values of real roots. For instance of a five-layered structure interpreted in Figure 4, the real roots are obtained from intersections of functions fλ and gλ, where the intermediate dielectric layers are assumed to be lossless.(35)fN=5λ=tanγ3l3;gN=5λ=right side of equation27i.Figure 4
A five-layered model with the bottom half-space as perfectly conducting.


## 3.1. Residue of the Poles for the Trapped-Surface-Wave-Term

Following the integration procedure addressed in [16], the trapped-surface-wave-term is defined by the sum of pole residues. In the N-layered structure, the terms of trapped surface wave due to a vertical electric dipole are written as(12)B0ϕSur=k022ω∑jZs0′λj∗H11λj∗ρeiγ0λj∗z+dλj∗2qN′λj∗γ0λj∗,E0ρSur=12∑jZs0′λj∗H11λj∗ρeiγ0λj∗z+dλj∗2qN′λj∗,E0zSur=i2∑jZs0′λj∗H01λj∗ρeiγ0λj∗z+dλj∗3qN′λj∗γ0λj∗,in which the function qN′ is expressed by(13)qN′=∂∂λγ0+k02ωμ0Zs0′0n=N−2,and λj∗j=1,2,… are the discrete roots of the pole equation of the electromagnetic field in the presence of the N-layered region, which is in a variant of the surface impedance function Zs0′0 at the air-to-dielectric boundary where n=N−2 represents the number of intermediate dielectric layers of the N-layered structure, written as(14)qN=γ0+k02ωμ0Zs0′0n=N−2=0.Considering the pole equation by (14) with substitution of (6) in the form of a recursive expression, it is necessary to derive the explicit pole equation first. Suppose (14) can be rewritten with respect to the function GNn, so that(15)qN=γ0+γ1k02k12GNnλ=0,in which GNn is defined by(16)GNnλ=k12ωμ0γ1Zs0′0.The superscriptn of function GNn represents the recursive order of the function that is in a variant of λ. In the meanwhile, (15) can be expanded if the order n≥1 by taking into account of (6) with the function GNn−1. Specifically, it is expressed by a set of sequential equations, as follows:(17)qN=2=γ0+γ1k02k12G20λ,qN=3=γ0+γ1k02k12G31λ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22G30λ,qN=4=γ0+γ1k02k12G42λ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22G41λ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22⋅tanh−iγ2l2+tanh−1γ3k22γ2k32G40λ,…qN=γ0+γ1k02k12GNnλ=γ0+γ1k02k12tanh−iγ1l1+tanh−1γ2k12γ1k22GNn−1λ.It is seen from equation (17) that the recursive order n of function GNn is reduced if the following identity is satisfied:(18)GNj+1λ=tanh−iγmlm+tanh−1γm+1km2γmkm+12GNjλ,m=N−j−1.Combining equations from (17), it is inferred that GNjλj=0≡1 is applied in each equation when n=N−2. Consequently, the equations are adapted by exploiting both (18) and GNjλj=0≡1, as follows:(19)N=3:tanhiγ1l1+G30λ=γ2k12γ1k22+G30λγ2k12γ1k22tanhiγ1l1,whenG30λ=1,(20)N=4:tanhiγ2l2+G41λ=γ2k12γ1k22+G41λγ2k12γ1k22tanhiγ2l2,whenG40λ=1,(21)N=5:tanhiγ3l3+G52λ=γ2k12γ1k22+G52λγ2k12γ1k22tanhiγ3l3,whenG50λ=1,…In (19), the explicit expression of pole equation in the three-layered structure is easy to obtain by G30=1. Similarly, for general case with n > 3, such as in the four-layered and five-layered structures by (20) and (21), respectively, where the function GNn−1 with respect to the variable λ has been applied, the explicit equation is defined by suppressing GN0=1. Through mathematical derivations, the variant coefficients are derived accordingly in expression of recursive form, as follows: (22)N≥4:G30λ=−γ0k12γ1k02,N≥5:G41λ=tanhiγ1l1+G30λγ2k12/γ1k22+G30λγ2k12/γ1k22tanhiγ1l1,N≥6:G52λ=tanhiγ2l2+G41λγ3k22/γ2k32+G41λγ3k22/γ2k32tanhiγ2l2.…Therefore, it is concluded that the general expression of pole equation for electromagnetic field in the presence ofN-layered region is derived from (14), in expression of(23)tanhiγnln+GNn−1λ=γ2k12γ1k22+GNn−1λγ2k12γ1k22tanhiγnln,whenGN0=1,where(24)GNjλ=tanhiγjlj+GNj−1λγj+1kj2/γjkj+12+GNj−1λγj+1kj2/γjkj+12tanhiγjlj.


## 3.2. Solution of Poles: A Graphical Approach

In order to evaluate the pole residues analytically, the set of equations from (19)–(21) are developed into explicit expressions with substitution of equations from (22). For convenience, the pole equation in the presence of three-layered region obtained from (19) is rewritten as(25)N=3:tanhiγ1l1=γ2k12γ1k22+γ0k12γ1k02−γ0k12γ1k02γ2k12γ1k22tanhiγ1l1.Analogously, the pole equation is derived for electromagnetic field of a vertical electric dipole in the presence ofN-layered region by substituting (23) with (24), iteratively. Specifically, the pole equations for the fields in the four-layered to six-layered structures can be written as follows:(26)N=4:tanhiγ2l2=γ3k22γ2k32+γ1k22γ2k12⋅γ1k02γ0k12×tanhiγ1l1−A1λtanhiγ1l1−B1λ−γ3k22γ2k32γ1k22γ2k12⋅γ1k02γ0k12×tanhiγ1l1−A1λtanhiγ1l1−B1λ⋅tanhiγ2l2,(27)N=5:tanhiγ3l3=γ4k32γ3k42+γ2k32γ3k22⋅γ0k12γ1k02×tanhiγ2l2−A2λtanhiγ2l2−B2λ−γ4k32γ3k42γ2k32γ3k22⋅γ0k12γ1k02×tanhiγ2l2−A2λtanhiγ2l2−B2λ⋅tanhiγ3l3,(28)N=6:tanhiγ4l4=γ5k42γ4k52+γ3k42γ4k32⋅γ0k12γ1k02×tanhiγ3l3−A3λtanhiγ3l3−B3λ−γ5k42γ4k52γ3k42γ4k32⋅γ0k12γ1k02×tanhiγ3l3−A3λtanhiγ3l3−B3λ⋅tanhiγ4l4.…By a few mathematical efforts, the functionsAjλ and Bjλ from (26)–(28) are derived through iterative substitutions, which are listed in Appendix, correspondingly. It is noted that the pole equation can be expressed in analogous form for the electromagnetic field of a vertical dipole in the presence of N-layered media. Accordingly, substituting the equations from (26)–(28) with functions Ajλ, the explicit expression of equations are obtained. By transpositions and rearrangements, the expression of pole equation in the presence of N-layered region can be described as(29)1+γn+1kn2γnkn+12⋅γn−1kn2γnkn−12⋅∑m=1n−1∏j=1mCn∗j,mtanhiγjlj∑k=1n−1∏j=1mCn∗m−j,n−mtanhiγjljtanhiγnln=γn+1kn2γnkn+12+γn−1kn2γnkn−12⋅1−∑m=1n−1∏j=1mCn∗j,mtanhiγjlj∑m=1n−1∏j=1mCn∗m−j,n−mtanhiγjlj,in which the coefficients Cn∗j,m, j = 1, …, m; m = 1, …, n are a set of coefficients defined by positive or negative ratio of γjkm and γmkj, which are determined from (26)–(28). If the bottom half-space is considered as perfectly conducting layer as depicted in Figure 3, the expression of (29) are simplified in condition of kn+1⟶∞, which reduces to(30)tanhiγnln=γn−1kn2γnkn−12⋅1−∑m=1n−1∏j=1mCn∗j,mtanhiγjlj1−∑m=1n−1∏j=1mCn∗m−j,n−mtanhiγjlj.Figure 3
TheN-layered model with the bottom half-space as conducting.It is seen from (29) and (30) that the obtained equations for the fields are expressed in form of addition, subtraction, and multiplication of functions tanhiγjlj with coefficients Cn∗, which offer advantages in terms of time savings with respect to standard numerical root finding procedures by (14). Specifically, the obtained pole equations for the electromagnetic field over a conductor coated by two-layered to four-layered dielectrics are in terms of the following equations:(31)N=4:tanhiγ2l2=γ1k22γ2k12⋅γ1k02γ0k12tanhiγ1l1−A1λtanhiγ1l1−B1λ,(32)N=5:tanhiγ3l3=γ2k32γ3k22⋅γ0k12γ1k02tanhiγ2l2−A2λtanhiγ2l2−B2λ,(33)N=6:tanhiγ4l4=γ3k42γ4k32⋅γ0k12γ1k02tanhiγ3l3−A3λtanhiγ3l3−B3λ.…It is convenient to validate (31) in derivation followed, which is exactly the same with (4.24) presented in [16] for the electromagnetic field of a vertical electric dipole on the surface of a conductor coated by two dielectric layers.(34)qN=4Per=γ0−iγ1k02k12tanγ1l1+tanh−1γ2k12γ1k22tanγ2l2=γ0−iγ1k02k12tanγ1l1+γ2k12γ1k22tanγ2l2×1−γ2k12γ1k22tanγ1l1tanγ2l2−1=0.Through the aforementioned analysis, the computational scheme to evaluate the residue sums for the trapped surface wave is proposed by exploiting the derived equations from (25)–(28). A graphical approach is applied in the present study to obtain initial values of real roots. For instance of a five-layered structure interpreted in Figure 4, the real roots are obtained from intersections of functions fλ and gλ, where the intermediate dielectric layers are assumed to be lossless.(35)fN=5λ=tanγ3l3;gN=5λ=right side of equation27i.Figure 4
A five-layered model with the bottom half-space as perfectly conducting.


## 4. Computation and Discussion

### 4.1. Illustrative Examples by a Graphical Approach

Primary objective of this section is to show the graphical approach to evaluate the poles by equations from (25)–(28) representing the three-layered to six-layered structures, respectively. In Figure 5 and 6, the real roots are obtained in lossless case, where the bottom half-space is considered as perfectly conducting and the electric lengths of intermediate dielectrics are identical. The computation coefficients of relative permittivity in the intermediate dielectrics are chosen by εr1=2.65, εr2=4.0, εr3=6.0, and εr4=8.0, respectively, where each electric length of the intermediate dielectrics are chosen by kili=0.45π and kili=1.45π, respectively.Figure 5
Graphical approach for root finding in evaluations of poles by equations from (25)–(27) at the operating frequency f = 100 MHz for the electromagnetic field due to a vertical dipole over a perfect conductor coated by multilayered lossless dielectrics, with (a) N = 3, kili=0.45π; (b) N = 4, kili=1.45π; (c) N = 5, kili=0.45π; (d) N = 5, kili=1.45π; (e) N = 3, kili=1.45π; (f) N = 4, kili=1.45π; (g) N = 5, kili=0.45π; and (h) N = 5, kili=1.45π.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)Figure 6
Graphical approach for root finding in evaluations of poles by (33) for the electromagnetic field due to a vertical dipole over a perfect conductor coated by four-layered lossless dielectrics, with (a) kili=0.45π; (b) kili=0.45π; (c) kili=1.45π; and (d) kili=1.45π.
(a)
(b)
(c)
(d)As illustrated in Figure5, the roots of equations from (25) to (27) are obtained from intersections of functions tanγnln and the right side terms gNλ on the condition of kn+1⟶∞. It is seen from Figures 5(c)–5(h), the terms of right side of (27)/i representing five-layered structure have imaginary part in the region of k1≤γj∗≤k2, since the function γ1=k12−γj∗21/2 is imaginary in the region of γj>k1. In Figure 6, the roots of (28) are plotted on condition of k5⟶∞ for electromagnetic field of a vertical dipole in the six-layered structure. In the computations, the intermediate dielectrics are chosen as lossless and the bottom half-space is perfectly conducting, whereas in the general cases where the roots may not be always real, the resulted roots are considered as initial values for numerical iterations.In Figure7, by applying the roots of poles, the electromagnetic field of a vertical dipole is computed, correspondingly. The total field, the trapped surface wave, and DRL waves (direct wave, reflected waves, and lateral wave) are computed with the same parameters in Figures 5(b) and 5(f), respectively, at the operating frequency f = 100 MHz. It is seen from Figure 7 that the trapped surface wave propagates along the air-to-dielectric boundary as the dominate wave on the condition of k1<k2<k3, and k3⟶∞, in the presence of four-layered region. When k0<λ<k1, γjλ=iλ2−kj21/2 is always a positive imaginary number, and the terms including the factor will attenuate exponentially as in the z direction. In Figure 8, the curves of the trapped surface wave versus propagation distance are plotted in three-layered to five-layered structures, respectively, when both the observation point and dipole source are placed on the surface of air-dielectric boundary. The electric lengths of each layer of intermediate dielectrics are identical with kili=0.45π. It is seen that field strength for trapped surface wave disperses as increased by layers of coating dielectrics.Figure 7
The electric componentE0zρ,z of electromagnetic field due to a vertical dipole on the conductor coated by two dielectric layers.Figure 8
The electric component of trapped surface waveE0zSρ,z generated by a vertical dipole in the presence of three-layered to five-layered regions.One would ask whether the derived equation in expression of (23) is advantageous in terms of computation over previous numerical procedures. This aspect is illustrated in Table 1, which shows the computation complexity and nested recursion times taken by the proposed method, and the previous numerical solution, to calculate the poles of the trapped surface wave generated by a vertical electric dipole lying on a horizontal layered medium at z = 0.Table 1
Comparison of computational feasibility and complexity for root finding algorithms. (n: represents the order of computations).
Method
Analytical solution by graphical approach
Computation complexity of root finding algorithm
Nested recursion
Three-layered case
In lossless case
Easy
1
Four-layered case
In lossless case
Complicated
2
N
≥
5
Unavailable
Complicated
N − 2
N
≤
6 by proposed equations (25)–(28)
In lossless case
Easy
0
N
>
6 by proposed equations (23) and (24)
In lossless case
Complicated
0


### 4.2. Interfering Behaviour of the Electromagnetic Field above the Surface of a Dielectric-Coated Conductor

To investigate the propagation properties of the electromagnetic field of a vertical electric dipole over a dielectric-coated conductor, the spatial distributions inρ^−z^ plane are plotted in Figures 9 to 11 at the operating frequency f = 100 MHz, when both the observation point and radiating source are located at the surface of air-dielectric boundary. In Figure 9, the DR fields (including the direct wave and ideal reflected wave) are computed by the first two terms in equations from (2)–(4). It is noted that the DR terms in the multilayered structures are the same with those for the uniform half-space. For a three-layered structure, the total fields (including the DR waves, lateral wave, and trapped surface wave) are plotted in Figure 10 with the relative permittivity and electric length of the intermediate dielectric chosen as εr1=2.65, and k1l1=0.45π, respectively. It is seen that the DR wave, lateral wave, and the trapped surface wave are combined in the total field to produce an interference pattern.Figure 9
Spatial distributions inρ^−z^ plane of the DR terms of the magnetic field component B0ϕ in dB due to a vertical electric dipole in the presence of N-layered regions.Figure 10
Spatial distributions inρ^−z^ plane of the electromagnetic field in dB due to a vertical electric dipole above a dielectric-coated conductor. (a) E0zρ,z and (b) B0Φρ,z.
(a)
(b)Figure 11
Spatial distributions inρ^−z^ plane of the electric field Ezρ,z in dB due to a vertical electric dipole above the conductor coated by two dielectric layers at f=100MHz; εr1=2.65, εr2=4, and z=d=0m with (a) k1l1=k2l2=0.45π; (b) k1l1=k2l2=1.45π; (c) trapped surface wave; and (d) lateral wave.
(a)
(b)
(c)
(d)In order to investigate the properties of the electromagnetic field, the total field (including the DR waves, lateral wave, and trapped surface wave) is computed in Figures11(a) and 11(b), respectively. The lateral wave and trapped surface wave are computed in Figures 11(c) and 11(d), respectively. The computation coefficients of relative permittivity in the intermediate dielectrics are chosen as the same in Figures 5(c) and 5(d), while each electric length of the intermediate dielectrics are chosen by kili=0.45π and kili=1.45π, respectively. It is seen from Figures 11(c) and 11(d) that in the z^ direction and the trapped surface, wave attenuates rapidly away from the surface of air and dielectrics.


## 4.1. Illustrative Examples by a Graphical Approach

Primary objective of this section is to show the graphical approach to evaluate the poles by equations from (25)–(28) representing the three-layered to six-layered structures, respectively. In Figure 5 and 6, the real roots are obtained in lossless case, where the bottom half-space is considered as perfectly conducting and the electric lengths of intermediate dielectrics are identical. The computation coefficients of relative permittivity in the intermediate dielectrics are chosen by εr1=2.65, εr2=4.0, εr3=6.0, and εr4=8.0, respectively, where each electric length of the intermediate dielectrics are chosen by kili=0.45π and kili=1.45π, respectively.Figure 5
Graphical approach for root finding in evaluations of poles by equations from (25)–(27) at the operating frequency f = 100 MHz for the electromagnetic field due to a vertical dipole over a perfect conductor coated by multilayered lossless dielectrics, with (a) N = 3, kili=0.45π; (b) N = 4, kili=1.45π; (c) N = 5, kili=0.45π; (d) N = 5, kili=1.45π; (e) N = 3, kili=1.45π; (f) N = 4, kili=1.45π; (g) N = 5, kili=0.45π; and (h) N = 5, kili=1.45π.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)Figure 6
Graphical approach for root finding in evaluations of poles by (33) for the electromagnetic field due to a vertical dipole over a perfect conductor coated by four-layered lossless dielectrics, with (a) kili=0.45π; (b) kili=0.45π; (c) kili=1.45π; and (d) kili=1.45π.
(a)
(b)
(c)
(d)As illustrated in Figure5, the roots of equations from (25) to (27) are obtained from intersections of functions tanγnln and the right side terms gNλ on the condition of kn+1⟶∞. It is seen from Figures 5(c)–5(h), the terms of right side of (27)/i representing five-layered structure have imaginary part in the region of k1≤γj∗≤k2, since the function γ1=k12−γj∗21/2 is imaginary in the region of γj>k1. In Figure 6, the roots of (28) are plotted on condition of k5⟶∞ for electromagnetic field of a vertical dipole in the six-layered structure. In the computations, the intermediate dielectrics are chosen as lossless and the bottom half-space is perfectly conducting, whereas in the general cases where the roots may not be always real, the resulted roots are considered as initial values for numerical iterations.In Figure7, by applying the roots of poles, the electromagnetic field of a vertical dipole is computed, correspondingly. The total field, the trapped surface wave, and DRL waves (direct wave, reflected waves, and lateral wave) are computed with the same parameters in Figures 5(b) and 5(f), respectively, at the operating frequency f = 100 MHz. It is seen from Figure 7 that the trapped surface wave propagates along the air-to-dielectric boundary as the dominate wave on the condition of k1<k2<k3, and k3⟶∞, in the presence of four-layered region. When k0<λ<k1, γjλ=iλ2−kj21/2 is always a positive imaginary number, and the terms including the factor will attenuate exponentially as in the z direction. In Figure 8, the curves of the trapped surface wave versus propagation distance are plotted in three-layered to five-layered structures, respectively, when both the observation point and dipole source are placed on the surface of air-dielectric boundary. The electric lengths of each layer of intermediate dielectrics are identical with kili=0.45π. It is seen that field strength for trapped surface wave disperses as increased by layers of coating dielectrics.Figure 7
The electric componentE0zρ,z of electromagnetic field due to a vertical dipole on the conductor coated by two dielectric layers.Figure 8
The electric component of trapped surface waveE0zSρ,z generated by a vertical dipole in the presence of three-layered to five-layered regions.One would ask whether the derived equation in expression of (23) is advantageous in terms of computation over previous numerical procedures. This aspect is illustrated in Table 1, which shows the computation complexity and nested recursion times taken by the proposed method, and the previous numerical solution, to calculate the poles of the trapped surface wave generated by a vertical electric dipole lying on a horizontal layered medium at z = 0.Table 1
Comparison of computational feasibility and complexity for root finding algorithms. (n: represents the order of computations).
Method
Analytical solution by graphical approach
Computation complexity of root finding algorithm
Nested recursion
Three-layered case
In lossless case
Easy
1
Four-layered case
In lossless case
Complicated
2
N
≥
5
Unavailable
Complicated
N − 2
N
≤
6 by proposed equations (25)–(28)
In lossless case
Easy
0
N
>
6 by proposed equations (23) and (24)
In lossless case
Complicated
0


## 4.2. Interfering Behaviour of the Electromagnetic Field above the Surface of a Dielectric-Coated Conductor

To investigate the propagation properties of the electromagnetic field of a vertical electric dipole over a dielectric-coated conductor, the spatial distributions inρ^−z^ plane are plotted in Figures 9 to 11 at the operating frequency f = 100 MHz, when both the observation point and radiating source are located at the surface of air-dielectric boundary. In Figure 9, the DR fields (including the direct wave and ideal reflected wave) are computed by the first two terms in equations from (2)–(4). It is noted that the DR terms in the multilayered structures are the same with those for the uniform half-space. For a three-layered structure, the total fields (including the DR waves, lateral wave, and trapped surface wave) are plotted in Figure 10 with the relative permittivity and electric length of the intermediate dielectric chosen as εr1=2.65, and k1l1=0.45π, respectively. It is seen that the DR wave, lateral wave, and the trapped surface wave are combined in the total field to produce an interference pattern.Figure 9
Spatial distributions inρ^−z^ plane of the DR terms of the magnetic field component B0ϕ in dB due to a vertical electric dipole in the presence of N-layered regions.Figure 10
Spatial distributions inρ^−z^ plane of the electromagnetic field in dB due to a vertical electric dipole above a dielectric-coated conductor. (a) E0zρ,z and (b) B0Φρ,z.
(a)
(b)Figure 11
Spatial distributions inρ^−z^ plane of the electric field Ezρ,z in dB due to a vertical electric dipole above the conductor coated by two dielectric layers at f=100MHz; εr1=2.65, εr2=4, and z=d=0m with (a) k1l1=k2l2=0.45π; (b) k1l1=k2l2=1.45π; (c) trapped surface wave; and (d) lateral wave.
(a)
(b)
(c)
(d)In order to investigate the properties of the electromagnetic field, the total field (including the DR waves, lateral wave, and trapped surface wave) is computed in Figures11(a) and 11(b), respectively. The lateral wave and trapped surface wave are computed in Figures 11(c) and 11(d), respectively. The computation coefficients of relative permittivity in the intermediate dielectrics are chosen as the same in Figures 5(c) and 5(d), while each electric length of the intermediate dielectrics are chosen by kili=0.45π and kili=1.45π, respectively. It is seen from Figures 11(c) and 11(d) that in the z^ direction and the trapped surface, wave attenuates rapidly away from the surface of air and dielectrics.


## 5. Conclusion

In summary, a computational scheme is presented to evaluate the residue sums for the fields of a vertical electric dipole oriented perpendicular to a stratified medium. The pole equation is reconstructed with explicit expression in the presence of theN-layered region. Therefore, a set of explicit equations is derived to enable analytical solutions on roots finding algorithms. Computations and analysis of electromagnetic field and its spatial distributions are carried out on and above a planar multidielectric-coated perfect conductor. The lateral wave and the trapped surface wave are combined in the total field to produce an interference pattern. It has been shown how the proposed equations are accurate and have significantly less computation complexity than previous numerical procedures.



---

*Source: 1657587-2019-10-20.xml*