# On the Radiation from a Short Current-Carrying Straight Wire Oriented Perpendicular to a Stratified Medium.

1. INTRODUCTIONThe study of the radiation characteristics of dipole antennas situated in proximity to stratified media has attracted the interest of several scientists in the past years [1-20]. This is because electric dipoles are extensively used in a variety of engineering applications, especially in the areas of close-to-the-surface radio communication, geophysical prospecting, radio remote sensing, and hypethermia [2-7,1115,22, 23]. Despite the number of contributions from researchers, to date accurate analytical expressions for the fields, valid regardless of the operating frequency, can be derived only under the assumption that both the source and the observation points lie on the surface of the material medium [10,23]. If this condition cannot be met, closed-form expressions for the fields may still be obtained, but at the price of introducing approximations that limit their applicability to electrically dense media [4-7,11] or specified frequency ranges [4-6,17-20]. Examples of contributions in this direction are the solutions valid for the quasi-static, indermediate, and far-field frequency ranges, which result from reducing the Sommerfeld integrals to forms amenable to asymptotic techniques such as the saddle-point method [13,17-20].

The present paper focuses on the problem of a short current-carrying straight wire oriented perpendicular to a stratified ground. The scope of the work is to overcome the limitations implied by the approach described in [10,23], and determine exact analytical expressions for the generated fields that are valid for any position of the source and observation points in the air space above the medium. To accomplish this task, first the direct field and the ideal reflected field are extracted from the integral representation for the magnetic vector potential. Next, the remaining part of the magnetic vector potential is written as the sum of branch-cut integrals and closed- contour integrals around the poles of the integrand. Finally, the hyperbolic branch cuts are extracted from the integrands of the branch-cut integrals and replaced with equivalent pole sets [23], so as to make it possible analytical integration. Each field component is given as the sum of two exponential terms (describing the direct and the ideal reflected fields) and three fast convergent sequences of residue sums, that represent the above surface ground wave, the lateral wave and the surface waves. The obtained formulas for the fields offer advantages in terms of time savings with respect to standard numerical integration procedures. Moreover, investigating the expressions for the different contributions to the fields, as well as their outcomes, allows to gain useful insight into the physics of the problem, and this feature makes the proposed solution advantageous over numerical simulation tools employed to solve electromagnetic boundary value problems.

2. FORMULATION OF THE PROBLEM

The geometrical configuration under study is sketched in Fig. 1. The emitter, modeled as a vertical electric dipole of moment 1 x [e.sup.j[omega]t] A x m, is located at height h above an N-layer conducting medium. The nth layer of the medium has dielectric permittivity [[epsilon].sub.n], electrical conductivity [[sigma].sub.n], and thickness [d.sub.n] = [z.sub.n] - [z.sub.n-1], being [z.sub.n] the position of the upper bound of the layer. The bottom or first layer is assumed to be a semi-infinite region. The magnetic permeability is everywhere equal to that of free space [[mu].sub.0] We introduce a cylindrical coordinate system ([rho], [phi], z) such that the air-medium interface coincides with the plane z=0 and the short wire antenna is placed on the z axis. The symmetry about the z axis implies that the electromagnetic field originates from a z-directed magnetic vector potential [17]. As a consequence, the field components produced in the air space may be expressed as

[E.sub.[rho]] = -j[omega]/[k.sup.2.sub.0] [[partial derivative].sup.2]A/[partial derivative][rho][partial derivative]z, (1)

[E.sub.z] = j[omega]/[k.sup.2.sub.0] ([[partial derivative].sup.2]/[partial derivative][[rho].sup.2] + 1/[rho] [partial derivative]/[partial derivative][rho]) A, (2)

[H.sub.[phi]] = -1/[[mu].sub.0] [partial derivative]A/[[partial derivative].sub.[rho]], (3)

where A is the non-null component of the vector potential, given by

[mathematical expression not reproducible], (4)

being

[mathematical expression not reproducible] (5)

and

R = [square root of ([[rho].sup.2] + [(z - h).sup.2])], R' = [square root of ([[rho].sup.2] + [(z + h).sup.2])], (6)

[u.sub.n] = - [square root of ([[lambda].sup.2] - [k.sup.2.sub.n])], Re[[u.sub.n]] > 0, (7)

[k.sup.2.sub.n] = [[omega].sup.2][[mu].sub.n][[epsilon].sub.n] - j[omega][[mu].sub.n][[sigma].sub.n]. (8)

In Eq. (4), the two exponential terms represent the primary field and the ideal field reflected by the medium (that is the field of a negative image), while [A.sup.s] is a correction term caused by the imperfect conductivity of the layered half-space. On the other hand, [Z.sub.n] is the intrinsic impedance of the nth layer (n=0 denotes free space), given by

[Z.sub.n] = [u.sub.n]/[[sigma].sub.n] + j[omega][[epsilon].sub.n], (9)

while [[??].sub.n] is the surface impedance at z=[z.sub.n], defined by the recursive formula

[[??].sub.n] = [Z.sub.n] [[??].sub.n-1] + [Z.sub.n] tanh ([u.sub.n] [d.sub.n])/[Z.sub.n] + [[??].sub.n-1] tanh ([u.sub.n] [d.sub.n]), n = 2, ..., N, (10)

with

[[??].sub.1] = [Z.sub.1]. (11)

The aim of this paper is to analytically evaluate Eq. (5) and, subsequently, the field components in Eqs. (1)-(3). To this end, it is first convenient to rewrite Eq. (5) as

[A.sup.s] ([rho], z, h) = [[mu].sub.0]/4[pi] [[integral].sup.[infinity].sub.- [infinity]] G([u.sub.0], [u.sub.1], [lambda], z, h) [H.sup.(1).sub.0] ([[lambda].sub.p])[lambda]d[lambda], (12)

with

[mathematical expression not reproducible], (13)

where the dependence on the square roots [u.sub.2], [u.sub.s], ..., [u.sub.N] has not been made explicit to signify that the function G is even in such quantities [17]. Next, proceeding as in [17,24] makes it possible to deform the integration contour so that it is wrapped around the singularities of the integrand located in the upper-half of the complex plane, as shown in Fig. 2(a). This leads to express [A.sup.s] as

[mathematical expression not reproducible], (14)

that is as the sum of two branch-cut integrals associated with -[k.sub.0] and - [k.sub.1] plus a number of closed contour integrals around the poles of G, that is [[lambda].sup.(-1).sub.1], [[lambda].sup.(-1).sub.2], [[lambda].sup.(-1).sub.3], .... Notice that the integral along the upper infinite semi-circle has been neglected in Eq. (14). In fact, this contribution cannot but be identically zero, because the Hankel function in the integrand decays exponentially with increasing [absolute value of ([lambda])] in the upper half-plane. The explicit form of [I.sup.(-1)] may be easily obtained through residue calculus. It yields

[I.sup.(-1)] = j[[mu].sub.0]/2 [summation over (m)] [A.sup.(-1).sub.m] (z, h) [H.sup.(1).sub.0] [[[lambda].sup.(-1).sub.m][rho]], (15)

where [A.sup.(-1).sub.m] is the residue of the function [lambda]G at [lambda]=[[lambda].sup.(-1).sub.m], that is

[mathematical expression not reproducible]. (16)

On the other hand, each of the integrals [I.sup.(0)] and [I.sup.(1)] may be evaluated by extracting the relevant branch cut from the integrand and replacing it with an equivalent set of poles. Branch-cut extraction procedure consists of decomposing G into its even and odd parts with respect to [u.sub.0] or [u.sub.1]. For instance, [I.sup.(0)] requires a decomposition of G with respect to [u.sub.0], as follows

G ([u.sub.0], [u.sub.1], [lambda], z, h) = [G.sup.(0).sub.e] ([lambda], z, h) + [G.sup.(0).sub.o] ([lambda], z, h) [u.sub.0], (17)

where [G.sup.(0).sub.[epsilon]] and [G.sup.(0).sub.o] are even function of [u.sub.0], given by

[G.sup.(0).sub.e] ([lambda], z, h) = i [G ([u.sub.0], [u.sub.1], [lambda], z, h) + G (-[u.sub.0], [u.sub.1], [lambda], z, h)], (18)

[G.sup.(0).sub.o] ([lambda], z, h) = 1/2[u.sub.0] [G ([u.sub.0], [u.sub.1], [lambda], z, h) -G (-[u.sub.0], [u.sub.1], [lambda], z, h)]. (19)

Since [G.sup.(0).sub.e] does not exhibit any branch cut running from -[k.sub.0], it will not contribute to [I.sup.(0)]. As a consequence, it yields

[mathematical expression not reproducible], (20)

and, analogously, for [I.sub.1] one obtains

[mathematical expression not reproducible], (21)

with

[G.sup.(1).sub.o] ([lambda], z, h) = 1/2[u.sub.1] [G ([u.sub.0], [u.sub.1], [lambda], z, h) -G ([u.sub.0], -[u.sub.1], [lambda], z, h)]. (22)

Pole sets equivalent to the extracted branch cuts may be determined by iteratively applying the square root algorithm to the computation of the square roots [u.sub.0] and [u.sub.1], as described in [10, 23,24]. At the Zth iteration, the square root algorithm generates the following partial-fraction representation for [u.sub.n] [23]

[mathematical expression not reproducible], (23)

being [M.sub.l]=[2.sup.l-1] - 1, and

[[lambda].sup.(n).sub.m] = -[square root of ([k.sup.2.sub.n] - [t.sup.2.sub.m])], (24)

[mathematical expression not reproducible], (25)

with

[t.sub.m] = tan (m[pi]/[2.sup.l]). (26)

Substituting Eq. (23) in the nth branch-cut integral [I.sup.(n)] and closing the integration contour with the arc of upper infinite semi-circle comprised between its ends, as sketched in Fig. 2(b), allows to obtain a closed-contour integral around the [[lambda].sup.(n).sub.m]'s. In fact, such poles lie on the suppressed branch line. Since the residue of Eq. (23) at [lambda]=[[lambda].sup.(n).sub.m] is equal to

[mathematical expression not reproducible], (27)

one obtains

[mathematical expression not reproducible], (28)

and, after letting [A.sup.(n).sub.m] (z, h) = [r.sub.mn] [G.sup.(n).sub.o] [[[lambda].sup.(n).sub.m], z, h], with n=0,1, expression (14) for [A.sup.s] may be rewritten in compact form as

[A.sup.s] ([rho], z, h) = j[[mu].sub.0]/2 [1.summation over (n=- 1)][[M.sub.l].summation over (m=1)] [A.sup.(n).sub.m] (z, h) [H.sup.(1).sub.0] [[[lambda].sup.(n).sub.m][rho]], (29)

where the first [M.sub.l] poles of the pole set corresponding to n= -1 have been considered. It is understood that the numbering of the [[lambda].sup.(-1).sub.m]'s is arranged so that Im [[lambda].sup.(-1).sub.1] < Im [[lambda].sup.(-1).sub.m] < Im [[lambda].sup.(- 1).sub.3] < ... and so on, which means that higher-order poles correspond to highly attenuated waves. Finally, use of Eqs. (4) and (29) in Eqs. (1)-(3) gives rise to the following expressions for the fields

[E.sub.[rho]] = [E.sub.[rho]0] [|.sub.[zeta]=z-h] - [E.sub.[rho]0] [|.sub.[zeta]=z+h] + [E.sup.s.sub.[rho]], (30)

[E.sub.z] = [E.sub.z0] [|.sub.[zeta]=z-h] - [E.sub.z0] [|.sub.[zeta]=z+h] + [E.sup.s.sub.z], (31)

[H.sub.[phi]] = H[[phi].sub.0][|.sub.[zeta]=z-h] - [H.sub.[phi]0][|.sub.[zeta]=z+h] + [H.sup.s.sub.[phi]], (32)

being

[mathematical expression not reproducible], (33)

[mathematical expression not reproducible], (34)

[mathematical expression not reproducible], (35)

with r=[square root of ([[rho].sup.2]+[[zeta].sup.2])], and

[E.sup.s.sub.[rho]] = - 1/2[omega][[epsilon].sub.0] [1.summation over (n=- 1)][[M.sub.l].summation over (m=1)] [[lambda].sup.(n).sub.m] [d.sup.A(n).sub.m]/dz [H.sup.(1).sub.1] [[[lambda].sup.(n).sub.m][rho]], (36)

[E.sup.s.sub.z] = 1/2[omega][[epsilon].sub.0] [1.summation over (n=-1)] [[M.sub.l].summation over (m=1)] [[[[lambda].sup.(n).sub.m]].sup.2] [A.sup.(n).sub.m][H.sup.(1).sub.0] [[[lambda].sup.(n).sub.m][rho]], (37)

[H.sup.s.sub.[phi]] = j/2 [1.summation over (n=-1)][[M.sub.l].summation over (m=1)] [[lambda].sup.(n).sub.m] [A.sup.(n).sub.m][H.sup.(1).sub.1][[[lambda].sup.(n).sub.m][rho]], (38)

with

[mathematical expression not reproducible] (39)

where account has been taken of the identities [21]

d[H.sup.(1).sub.0]([lambda][rho])/d[rho] = - [[lambda].sup.(1).sub.1]([lambda][rho]), (40)

([[partial derivative].sup.2]/[partial derivative][[rho].sup.2] + 1/[rho] [partial derivative]/[partial derivative][rho]) [H.sup.(1).sub.0]([lambda][rho]) = -[[lambda].sup.2][H.sup.(1).sub.0] ([lambda][rho]). (41)

Each of the expressions (36)-(38) consists of three sequences of residue sums, which converge as l is increased in virtue of convergence of the square root algorithm. The residue sums associated with n=0 and n=1 describe the above-surface ground and lateral waves, respectively. Such contributions, together with the direct field and the ideal reflected field, labelled with a subscript 0 in Eqs. (30)-(32), constitute the space wave radiated by the antenna. The rest of the generated field consists of surface waves, and is described by the residue sums associated with n= -1 in Eqs. (36)-(38). Expressions (30)-(32) allow to overcome the limitations implied by the previously published square-root- algorithm-based solution to the same problem [10, 24]. In fact, they allow to relax the assumption that both source and observation points lie at the air-medium interface. At the same time, as will be clarified in Sec. 3, the derived expressions ensure a high level of accuracy and offer advantages in terms of time savings with respect to standard numerical techniques conventionally used for evaluating Sommerfeld- type integrals.

3. RESULTS AND DISCUSSION

Primary objective of this section is to show the convergence of the sequences in Eqs. (36)-(38) as the number of iterations l of the square root algorithm is increased. To this goal, we use Eq. (31) to compute the amplitude-frequency spectrum of the [E.sub.z]-field generated by a unit- moment straight wire antenna lying on the top surface of a two-layer conductive ground. The field point is at height z=50m above the medium and at radial distance [rho]=300 m from the emitter, while the electromagnetic and geometric parameters of the layers of the medium are assumed to be [[sigma].sub.1]=10mS/m, [[epsilon].sub.1]=10[[epsilon].sub.0], [d.sub.1]=400m, [[sigma].sub.2]=1 S/m, and [[epsilon].sub.2]=5[[epsilon].sub.0]. The results of the simulations, depicted in Fig. 3, are compared with the data arising from numerically evaluating Eq. (2) through the Gauss-Kronrod quadrature rule [25]. As is seen, when l is increased the outcomes from Eq. (31) rapidly converge to the exact numerical data, denoted by points, and convergence is faster at higher frequencies. Moreover, perfect agreement between analytical and numerical data is observed when truncating the sequences of residue sums in Eq. (37) at l=12, that is when each residue sum is composed of [M.sub.12] = [2.sup.11] - 1 = 2047 residues. It should be noted that, in general, the three contributions to Eq. (37) do not have the same impact on the [E.sub.z]-field strength. For instance, if the field point is located at height z=50 m from the top surface of the medium, one would expect that the effects of surface and lateral waves, associated with n= -1 and n=1, respectively, are negligible with respect to that of the ground wave (n=0). This point is clarified by Fig. 4, which depicts the magnitude of [E.sub.z] and those of the relevant ground-wave, lateral-wave, and surface- wave terms against the radial distance [rho] from the emitter. The parameters of the medium are taken to be the same as in the previous example, and the operating frequency is assumed to be 10 kHz. Again, the sequences of the residue sums have been truncated at l=12. From the analysis of the plotted curves it emerges that the amplitude of the [E.sub.z]-field substantially coincides with that of the ground-wave field over the considered interval, and that the surface waves suffer from a mild attenuation with increasing [rho]. Instead, the magnitude of the lateral wave rapidly goes down below the plot scale as [rho] grows up. Finally, it should be observed that direct and ideal reflected fields cannot contribute to the total [E.sub.z]- field, since they cancel out each other because of the position of the emitter at the air-medium interface.

One would ask whether the derived solution is advantageous in terms of computation time over previous solutions and standard numerical integration procedures. This aspect is illustrated in Table 1, which shows the time taken by the proposed method, the Gauss-Kronrod quadrature rule, and King's well-established solution [7], to calculate the radial distribution of the vertical electric field that a straight wire antenna lying on a homogeneous medium produces at z=50 m. Here, the operating frequency is 1 MHz, while the electromagnetic parameters of the medium are taken to be [[sigma].sub.1]=0.1 mS/m and [[epsilon].sub.1]=3 [[epsilon].sub.0].

Table 1 is structured as follows. The second column contains time cost data corresponding to simulations running on a single-core 1.8 GHz PC. The third column shows the ratio of the time taken by the proposed expression for [E.sub.z] and the Gaussian rule to that required by King's solution (speed-up offered by King's approach). Finally, the fourth column contains the root-mean- square (RMS) relative errors resulting from using the proposed and King's solution in place of Gaussian quadrature rule. From a glance at the data in Table 1 it is deduced that, accuracy being equal, using Eq. (37) rather than Gaussian quadrature permits to obtain significant time savings. On the other hand, King's approach implies minimum time cost, but the price to be paid for its usage is the loss of accuracy. This is also confirmed by the contour plots sketched in Fig. 5, which depict the magnitudes of [E.sub.[rho]] (Fig. 5(a)) and [E.sub.z] (Fig. 5(b)) versus [rho] and z (the antenna is still positioned at height h=0). For both the contour plots, a significant discrepancy is observed between the solid isolines associated with the proposed solution and the dashed isolines arising from King's approach. The amount of the discrepancy depends on the conductivity of the medium, and higher values of the conductivity make it possible to reduce the discrepancy. For instance, assuming [[sigma].sub.1]=1 mS/m and [[epsilon].sub.1]=10 [[epsilon].sub.0] leads to the contour plots illustrated in Fig. 6, characterized by solid and dashed isolines almost overlapping. Bad accuracy of King's solution for low values of the conductivity is explained by the failure of the assumptions [absolute value of ([k.sub.1])][rho][much greater than]1 and [k.sup.2.sub.0]/[absolute value of ([k.sup.2.sub.1])] [much less than] underlying its derivation.

4. CONCLUSION

This paper presents an explicit solution, in series form, for the fields from a short current-carrying straight wire oriented perpendicular to a stratified medium. The solution arises from the analytical evaluation of the integral representation for the magnetic vector potential, carried out through contour integration upon substituting an equivalent pole set for each branch cut exhibited by the integrand. Each field component is expressed as the sum of five contributions, namely the direct, the ideal reflected, the ground-wave, the lateral-wave, and the surface-wave fields. The latter three terms are given as convergent sequences of residue sums. Performed numerical simulations confirm the correctness of the developed theory. It has been shown how the proposed formulas are accurate and significantly less time consuming than conventional numerical integration procedures.

REFERENCES

[1.] Sommerfeld, A., "On the propagation of waves in wireless telegraphy," Annals of Physics, Vol. 81, 1135-1153, 1926.

[2.] Banos, A., Dipole Radiation in the Presence of a Conducting Half-space, Pergamon Press, Oxford, NY, 1966.

[3.] Durrani, S. H., "Air to undersea communication with electric dipoles," IRE Transactions on Antennas and Propagation, Vol. 10, 524-528, 1962.

[4.] Moore, R. K. and W. E. Blair, "Dipole radiation in a conducting half space," Journal of Research of the National Bureau of Standards-D Radio Propagation, Vol. 65D, No. 6, 547- 563, 1961.

[5.] Wait, J. R., "The electromagnetic fields of a horizontal dipole in the presence of a conducting half-space," Canadian Journal of Physics, Vol. 39, 1017-1028, 1961.

[6.] Bannister, P. R., "Quasi-static fields of dipole antennas at the earth's surface," Radio Science, Vol. 1, 1321-1330, 1966.

[7.] King, R. W. P., M. Owens, and T. T. Wu, Lateral Electromagnetic Waves: Theory and Applications to Communications, Geophysical Exploration, and Remote Sensing, Springer-Verlag, New York, 1992.

[8.] King, R. W. P., S. S. Sandler, and L. C. Shen, "The electromagnetic field of a vertical electric dipole over the earth or sea," IEEE Transactions on Antennas and Propagation, Vol. 42, 382-389, 1994.

[9.] King, R. W. P., "Electromagnetic field of a vertical dipole over an imperfectly conducting halfspace," Radio Science, Vol. 25, 149-160, 1990.

[10.] Parise, M., "An exact series representation for the EM field from a vertical electric dipole on an imperfectly conducting half-space," Journal of Electromagnetic Waves and Applications, Vol. 28, No. 8, 932-942, 2014.

[11.] Wait, J. R., Electromagnetic Waves in Stratified Media, Pergamon Press, New York, 1970.

[12.] Parise, M., "Second-order formulation for the quasi-static field from a vertical electric dipole on a lossy half-space," Progress In Electromagnetics Research, Vol. 136, 509-521, 2013.

[13.] Kong, J. A., Electromagnetic Wave Theory, John Wiley & Sons, New York, 1986.

[14.] Ward, S. H. and G. W. Hohmann, "Electromagnetic theory for geophysical applications," Electromagnetic Methods in Applied Geophysics, Vol. 1, 130-311, Society of Exploration Geophysicists, Tulsa, OK, 1988.

[15.] Balanis, C. A., Antenna Theory: Analysis and Design, 2nd Edition, John Wiley & Sons, New York, 1997.

[16.] Wait, J. R., "Radiation from a vertical electric dipole over a stratified ground," Transactions of the IRE--Antennas and Propagation, Vol. 1, No. 1, 9-11, 1953.

[17.] Wait, J. R., "Influence of a sub-surface insulating layer on electromagnetic ground wave propagation," IEEE Transactions on Antennas and Propagation, Vol. 14, No. 6, 755-767, 1966.

[18.] Wait, J. R., "Asymptotic theory for dipole radiation in the presence of a lossy slab lying on a conducting half-space," IEEE Transactions on Antennas and Propagation, Vol. 15, No. 5, 645-648, 1967.

[19.] Kong, J. A., L. C. Shen, and L. Tsang, "Field of an antenna submerged in a dissipative dielectric medium," IEEE Transactions on Antennas and Propagation, Vol. 25, No. 6, 887-889, 1977.

[20.] Fuller, J. A. and J. R. Wait, "A pulsed dipole in the earth," Transient Electromagnetic Fields, 237-269, Springer, New York, 1976.

[21.] Parise, M., "An exact series representation for the EM field from a circular loop antenna on a lossy half-space," IEEE Antennas and Wireless Prop. Letters, Vol. 13, 23-26, 2014.

[22.] Parise, M., "Exact EM field excited by a short horizontal wire antenna lying on a conducting soil," AEU-International Journal of Electronics and Communications, Vol. 70, 676-680, 2016.

[23.] Parise, M., "Improved Babylonian square root algorithm-based analytical expressions for the surface-to-surface solution to the Sommerfeld half-space problem," IEEE Transactions on Antennas and Propagation, Vol. 63, No. 12, 5832-5837, 2015.

[24.] Parise, M., "On the surface fields of a small circular loop antenna placed on plane stratified earth," International Journal of Antennas and Propagation, Vol. 2015, 1-8, 2015.

[25.] Patterson, T. N. L., "The optimum addition of points to quadrature formulae," Mathematics of Computation, Vol. 22, 847-856, 1968.

Mauro Parise (1), * and Giulio Antonini (2)

Received 9 April 2017, Accepted 2 June 2017, Scheduled 30 June 2017

* Corresponding author: Mauro Parise (m.parise@unicampus.it).

(1) Unit of Electrical Engineering, University Campus Bio-Medico of Rome, Italy. (2) Department of Industrial Engineering, Computer Science and Economics, University of L'Aquila, Italy.

Caption: Figure 1. Sketch of a short straight wire oriented perpendicular to a stratified medium.

Caption: Figure 2. (a) Deformed integration contour and (b) branch cut equivalent pole set.

Caption: Figure 3. Profiles of [absolute value of ([E.sub.z])] against frequency. Outcomes from the proposed solution are compared with the data arising from numerical integration of the field integral.

Caption: Figure 4. Profile of [absolute value of ([E.sub.p])] against [rho]. Separate contributions of ground wave (gw), lateral wave (lw), and surface waves (sw) to the total field.

Caption: Figure 5. Contour plots of (a) [absolute value of ([E.sub.[rho]])] and (b) [absolute value of ([E.sub.z])] versus [rho] and z, computed for [[sigma].sub.1]=0.1 mS/m and [[epsilon].sub.1]=3 [[epsilon].sub.0] by applying the proposed approach (solid isolines). Outcomes from King's approach are denoted by dashed isolines.

Caption: Figure 6. Contour plots of (a) [absolute value of ([E.sub.[rho]])] and (b) [absolute value of ([E.sub.z])] versus [rho] and z, computed for [[sigma].sub.1]=1 mS/m and [[epsilon].sub.1]=3 [[epsilon].sub.0] by applying the proposed approach (solid isolines). Outcomes from King's approach are denoted by dashed isolines.

Table 1. Computation times and accuracy of the proposed and conventional approaches. Approach Time Speed-up RMS rel. consumption error (%) (s) Gaussian rule 126 21 -- King's approach 6 -- 8.3 Proposed solution 19 3.2 0.16

Printer friendly Cite/link Email Feedback | |

Author: | Parise, Mauro; Antonini, Giulio |
---|---|

Publication: | Progress In Electromagnetics Research |

Article Type: | Report |

Date: | Jun 1, 2017 |

Words: | 4298 |

Previous Article: | Layer-by-Layer Design of Bianisotropic Metamaterial and Its Homogenization. |

Next Article: | Electromagnetic Theory with Discrete Exterior Calculus. |

Topics: |