Analysis of an Iterative Approach to Determine the Current on the Straight Infinite Wire Above Ground
- Otto von Guericke University, Magdeburg, Germany
Correspondence: Felix Middelstaedt (firstname.lastname@example.org)
An iterative approach which was recently applied to approximate the reflection and scattering coefficients of transmission line ports is analyzed. The iterative solution for the current on an infinite wire above ground is compared to the exact solution. The example is chosen since it is one of the few problems where an exact solution exists. The wire is excited by a lumped voltage source or a plane wave. The convergence of the iterative approach is shown. It can be concluded that the zeroth iteration, which is the classical transmission line solution, coincides with the general transverse electromagnetic mode. Furthermore, it is shown that the first iteration is a very good approximation of the radiation and leaky modes, that occur in the close neighborhood around the lumped source.
Transmission lines play an important part in many branches of electrical engineering. Hence, it is of great interest to understand the propagation of waves along transmission lines and develop analytic methods to determine the currents on wires.
Recently, an iterative approach was developed for thin wires above ground and applied to many examples (see e.g. Rachidi and Tkachenko, 2008; Middelstaedt et al., 2015, 2018). The iterative approach can be seen as an extension to the classical transmission line approximation, which is restricted to uniform transmission lines and to wavelengths much larger than the transverse dimensions of the transmission lines. After only one iteration the iterative approach already yields accurate results for infinite and semi-infinite problems respecting the non-uniformity of the wire.
Due to the success of the iterative approach, questions that arise are:
Does the iterative approach converge if more iterations are used?
What is the interpretation of the solution of each iteration step?
A first step in answering these fundamental questions is done in this paper. The questions are answered using the example of the infinite wire above ground excited by a lumped voltage source and a plane wave. This specific problem is chosen since its exact analytic solution is available for comparison.
Figure 1 shows the thin wire above an infinite, planar, perfectly electrically conducting ground. It is assumed that the radius a of the wire is much smaller than its height h above ground. Hence, the thin wire approximation can be applied. The wire is either excited by a lumped voltage source with amplitude V0 or by a plane wave with wave number and electric field vector E0. The angle of incident is denoted by θ and is defined as shown in Fig. 1.
The current density vector, that arises due to the sources, is defined to be pointing in the positive z-direction. The corresponding current is denoted by I. The current I and the electric scalar potential Φ can both be determined using the so call mixed potential integral equations (MPIE). The MPIE are derived in many publications e.g. in Nitsch and Tkachenko (2010),
The operator 𝒢 is defined as a convolution integral
with the Green's function for the half space
The following relation holds for the angular frequency ω, the permittivity ε0, and the permeability μ0
The tangential electric field on the surface of the wire Etan(z) depends on the kind of excitation. For the two considered cases it is given as
The delta distribution is denoted by δ.
Usually, the goal is to solve equations similar to Eq. (1) for the current I as in Leviatan and Adams (1982), Rachidi and Tkachenko (2008), and Middelstaedt et al. (2015, 2018). It is possible to find the exact solution for the straight infinite wire above ground. For more practical cases there are no exact analytic solutions but only approximations. The iterative approach yields one of these approximations.
In the following the exact solution and the iterative one are derived for the described example. The convergence of the iterative approach is analyzed. Furthermore, the iterative solution is interpreted by comparing it to the exact one.
The problem of the infinite wire above ground is similar to the two wire transmission line in free space, which is dealt with in Schelkunoff (1956), Marin (1975), and Leviatan and Adams (1982). In Marin (1975) the focus lies on the determination of the electromagnetic fields that propagate along the wires in time domain. In Leviatan and Adams (1982) the authors focus only on the transverse electromagnetic (TEM) and leaky modes and ignore the radiation mode. Therefore, the solution for the current along the wire is derived here again with the focus on the current and all modes.
To solve Eq. (1) it is convenient to apply the Fourier transform. The Fourier transform and its inverse, which is marked with a tilde, are defined as
where Km denotes the modified Bessel function of the second kind and order m.
After some mathematical manipulations the current I and its transform are found to be
With Eq. (12) an explicit expression is given for the current that arises due to the excitation Etan(z). Depending on the Fourier transform of the excitation the integral can be solved. The transforms of the excitations given in Eq. (5) are
Hence, the solution resulting from the plane wave excitation is simply
2.1 Lumped Excitation
The procedure to obtain the solution for the lumped excitation is more advanced. The difficulty lies in the solution of the integral in Eq. (12), namely
The integral in Eq. (16) is analyzed using the residue theorem as in Marin (1975), Leviatan and Adams (1982). The variable kz is assumed to be complex. Figure 2 shows the integration contour for Eq. (16) qualitatively in the complex plane for z>0. Since the problem is symmetric it is sufficient to only consider the solution for z>0, which will be done from now on.
The integrand has two branch cuts due to the square root in Eq. (10). The branch cuts in Fig. 2 connect the branch points at via kz=∞ (ref. Riemann sphere). The choice of the paths of the branch cuts originating from the fixed branch points is arbitrary. Here these specific branch cuts are chosen to allow a quick numerical integration along the branch cuts for later considerations.
The residue theorem states that
where Res(kn) denotes the residue of the integrand in Eq. (16) for the poles kn inside the contour (see Fig. 2). The integration along the outer circle vanishes for R→∞ mainly due to the exponential function in the numerator and z>0.
The integral around the branch point at results in the TEM mode ITEM(z). Each residue in Eq. (19) gives rise to leaky modes ILeaky(z). The integration along the branch cut (in both directions) is connected to the radiation mode IRad(z). The total current can be expressed as a sum of these modes
for z>0 with
The poles kn (the roots of Eq. 10) that contribute to the solution are all located between the branch cut and the imaginary axis for this particularly chosen branch cut. There are poles, whose real part is smaller than −k. But they are part of a different Riemann surface. Hence, they are not part of the total current.
The radiation modes arise due to the integration along the branch cuts, specifically along the paths C− and C+.
With the definition of the branch cut (see Fig. 2) follows
The integral can be modified with the substitution , dkz=j dx.
A solution for the latter integral is not known. But the integral is suitable for a fast numerical integration, especially for large z.
The treatment of the radiation mode is omitted in Leviatan and Adams (1982). In Marin (1975) the branch cut is chosen differently and hence, a different representation of the radiation mode is obtained.
The obtained exact solution is later used to interpret the iterative solution.
for the iteration start and
for each iteration step with and
The total current and total scalar potential are
Equations (28) coincide with the classical telegrapher's equations. The equations for each iteration step (Eq. 29) have a similar form. But the source is the electric field emitted by the previously determined current I(n−1).
The equations can be solved using the Fourier transform as for the exact solution. The results are
Hence, the resulting total current is
The geometric series converges if and only if
A parameter study found that Eq. (39) is always true as long as the thin wire approximation is applicable, meaning a≪h and . Further details are illustrated in Appendix B. Hence, the total current from the iterative approach coincides with the exact result in Eq. (11) and therefore, converges to the correct solution.
For general wires it is not so straight forward to obtain the total current. Usually the zeroth and first iteration are determined analytically and the total current is approximated by the sum of those two currents. Therefore, this procedure is applied to this problem and the solutions are compared to the exact modes.
3.1 Lumped Excitation
To find the current of the zeroth and first iteration for the lumped excitation it is more convenient to not apply the Fourier transform. Thus, the complicated integral from the inverse Fourier transform is avoided.
It can be shown that
by applying integration by parts. Decoupling Eqs. (28) and (29) then yields
where Id denotes the identity operator.
To present the solution for each iteration in a concise form it is convenient to define a linear operator 𝒦 by
The latter equations hold for any excitation. For the lumped voltage source follows with Eq. (5)
This section deals with the comparison of the zeroth and first iteration with the exact solution.
4.1 Lumped Excitation
The zeroth iteration , which is the classical transmission line solution, coincides with the exact TEM solution for the lumped excitation (Eq. 21). The first iteration current and the sum of the radiation and leaky modes are illustrated in Figs. 3 and 4. The currents are normalized with the TEM mode. It can be seen that even for large kh the first iteration is a very good approximation of the radiation and leaky modes.
In Appendix A it is shown that the asymptotic expansion of the radiation mode and the first iteration coincide. Due to the exponential damping of the leaky modes, their asymptotic expansion is zero.
This explains the improved accuracy of the approximation with increasing distance to the source.
4.2 Plane Wave Excitation
It is notable that for θ→0 follows
Hence, the iterative approach immediately converges to the exact solution when θ→0. Furthermore, the solution in Eq. (64) coincides with the TEM solution from the lumped excitation with an equivalent voltage . This equivalent behavior is already mentioned in Nitsch and Tkachenko (2002b).
The manuscript presents the exact and iterative solution for the current of an infinite wire above ground excited by a lumped and distributed source. The iterative approach approximates the system behavior independently of the excitation. It is shown that the iterative approach converges for reasonable wire dimensions. When there is a TEM mode present, it coincides with the zeroth iteration solution, which is know from classical transmission line theory. The first iteration, which can also be determined for arbitrary wires, approximates the radiation and leaky modes. They occur when a lumped disturbance, e. g. a lumped voltage source, is present on the wire.
Even though only an infinite wire was analyzed the results help to understand the propagation of currents along practical wires with discrete disturbances (lumped elements, sharp bend, … ). Furthermore, the manuscript presents another example where the iterative approach yields a very good approximation without much effort.
No data sets were used in this article.
The first iteration current is given in Eq. (57). As described in Abramowitz and Stegun (1972) the amplitude of the exponential integral E1 vanishes for large arguments. For small arguments a series exists.
where γ denotes the Euler-Mascheroni constant. Inserting the latter results into Eq. (57) yields
The radiation mode consists of an integral of the form
After some lengthy but straight forward mathematical manipulation it can be shown that
Inserting the latter results into Eq. (27) yields
In the following some results of the parameter study are illustrated. As stated in Sect. 3 the function
needs to be smaller than 1 for the iterative approach to be convergent.
Figure B1 shows the function q(kz) for different k. It can be seen that the convergence condition Eq. (39) holds as long as k≤klim. The upper bound klim is chosen such that , where λ denotes the wavelength. This results in
If k increases, the thin wire approximation is not valid anymore, q(kz) is larger than 1 for some kz and the iterative approach diverges.
A lengthy but straight forward analytic analysis shows that q(kz) is always smaller than 1 for all kz>k.
The authors declare that they have no conflict of interest.
This article is part of the special issue “Kleinheubacher Berichte 2018”. It is a result of the Kleinheubacher Tagung 2018, Miltenberg, Germany, 24–26 September 2018.
This paper was edited by Frank Gronwald and reviewed by two anonymous referees.
Abramowitz, M. and Stegun, I. A. (Eds.): Handbook of Mathematical Functions, With Formulas, Graphs and Mathematical Tables, National Bureau of Standards, 10 edn., U.S. Government Printing Office Washinton D.C., USA, 1972. a, b
Haase, H.: Full-Wave Field Interactions of Nonuniform Transmission Lines, Phd thesis, Otto von Guericke University, Magdeburg, Germany, 2005. a
Leviatan, Y. and Adams, A.: The response of a two-wire transmission line to incident field and voltage excitation, including the effects of higher order modes, IEEE T. Antenn. Propag., 30, 998–1003, https://doi.org/10.1109/TAP.1982.1142893, 1982. a, b, c, d, e, f
Middelstaedt, F., Tkachenko, S., Vick, R., and Rambousky, R.: Analytic approximation of natural frequencies of bent wire structures above ground, in: 2015 IEEE International Symposium on Electromagnetic Compatibility (EMC), 812–817, Dresden, Germany, 16–22 August 2015, https://doi.org/10.1109/ISEMC.2015.7256268, 2015. a, b
Middelstaedt, F., Tkachenko, S. V., Rambousky, R., and Vick, R.: High-Frequency Electromagnetic Field Coupling to a Long, Finite Wire With Vertical Risers Above Ground, IEEE T. Electromagn. C., 58, 1169–1175, https://doi.org/10.1109/TEMC.2016.2544110, 2016. a
Middelstaedt, F., Tkachenko, S. V., and Vick, R.: Transmission Line Reflection Coefficient Including High-Frequency Effects, IEEE T. Antenn. Propag., 66, 4115–4122, https://doi.org/10.1109/TAP.2018.2839914, 2018. a, b, c
Nitsch, J. and Tkachenko, S.: Complex-Valued Transmission-Line Parameters and their Relation to the Radiation Resistance, Interaction Note 573, Air Force Weapons Laboratory, Kirtland Air Force Base, Albuquerque, NM, USA, available at: http://ece-research.unm.edu/summa/notes/In/0573.pdf (last access: 11 April 2019), 2002a. a
Nitsch, J. and Tkachenko, S.: Source Dependent Transmission Line Parameters – Plane Wave vs. TEM Excitation, Interaction Note 577, Air Force Weapons Laboratory, Kirtland Air Force Base, Albuquerque, NM, USA, available at: http://ece-research.unm.edu/summa/notes/In/0577.pdf (last access: 11 April 2019), 2002b. a
Rachidi, F. and Tkachenko, S. V. (Eds.): Electromagnetic Field Interaction with Transmission Lines: From Classical Theory to HF Radiation Effects, WIT Press, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK, 2008. a, b, c, d
Schelkunoff, S. A.: Electromagnetic Waves, D. Van Norstrand Company, Inc., New York, NY, 8th printing, 1956. a
Weinstein, L. A.: The Theory of Diffraction and the Factorization Method, The Golem Press, Boulder, Colorado, 1969. a
Wong, R.: Asymptotic Approximations of Integrals, SIAM, Philadelphia, PA, 2002. a