https://doi.org/10.5194/ars-17-177-2019
https://doi.org/10.5194/ars-17-177-2019

19 Sep 2019

19 Sep 2019

# Singularity Expansion Method for thin wires and the Method of Modal Parameters

Singularity Expansion Method for thin wires and the Method of Modal Parameters
Sergey V. Tkachenko1, Juergen B. Nitsch1, Felix Middelstaedt1, Ronald Rambousky2, Martin Schaarschmidt2, and Ralf Vick1 Sergey V. Tkachenko et al.
• 1Faculty of Electrical Engineering and Information Technology, Otto-von-Guericke University, 39106 Magdeburg, Germany
• 2Wehrwissenschaftliche Institut für Schutztechnologien – ABC-Schutz (WIS), 29633 Münster, Germany

Correspondence: Sergey V. Tkachenko (sergey.tkachenko@ovgu.de)

Abstract

Here, we describe a technique to define the Singularity Expansion Method (SEM) poles for short-circuited thin-wire structures developed using the Method of Modal Parameters (MoMP). The MoMP method consists of in the expansion of the system of mixed-potential integral equations (MPIE) into the Fourier series, including the kernels containing Green's function. Corresponding equations for Fourier modes contain infinite matrices of p.u.l. inductance and capacitance, and the solution for current can be obtained using the infinity matrix of p.u.l. impedance. The SEM poles are given by the zeros of the determinant of this matrix. For the case of the symmetrical circular loop, this equation transforms to one well-know from the literature. Numerical investigation of solutions for the poles of the first layer has shown good agreement with previously obtained analytical and numerical results for different wire configurations.

1 Introduction

Thin-wire transmission lines play an important role in EMC. Thin-wire transmission lines facilitate the transmission of the desired signals between electronic devices of different kinds. On the other hand, they are subject to different kinds of electronic interferences. Induced voltages are often the cause of the failure of electronic devices. Different numerical methods (e.g., MoM, FDTD) can be applied to calculate induced currents and voltages, but these are not very helpful for gaining insight into the physics of coupling phenomena (e.g., to make a qualitative analysis of the system response for various values of system parameters, including the case of a statistically defined system or field parameters).

In contrast, the Singularity Expansion Method (SEM) (see pioneer paper Baum, 1971, and reviews of results in books Baum et al., 2012; Tesche et al., 1997) represents the scattering object as a set of oscillators, in which complex frequencies are poles of the response function, which do not depend on the type of excitation of the system. This set of poles yield the main contribution to the coupling response of the object, both in the frequency and time domains. It also defines the radiation of the system and the scattering amplitude. Moreover, the SEM poles can be used for the identification of the systems. Recently, this method has attracted increasing interest in connection to the problem of target identification (Myers et al., 2011a, b; Giri and Tesche, 2012). However, the main method for obtaining the SEM expansion for a general system is the processing of numerical (e.g., Method of Moments) data (Senior and Pond, 1981; Singaraju et al., 1976).

On the other hand, the application of SEM to thin-wire systems (antennas and transmission lines) has several specific features that simplify the investigation of the SEM expansion, especially when the analytic form of the response function is known. Such response functions were investigated earlier for the straight transmission lines, when the length of the line is essentially larger relative to its height (Tkachenko et al., 2013, 2014, 2016; Middelstaedt et al., 2016, 2018). To do this, the so-called asymptotic approach was used (Tkachenko et al., 2001; Rachidi and Tkachenko, 2008). However, the application of the asymptotic approach is restricted by long straight lines. In the present work, we describe the application of another analytical method, the method of modal parameters (Nitsch and Tkachenko, 2005, 2007) for the investigation of the SEM poles. This analytical technique can be applied for thin wires of arbitrary geometrical form and arbitrary excitations.

The rest of the paper is organized in the following. In Sect. 2, the method of modal parameters will be briefly described. In Sect. 3, we present the analytical approach and a numerical algorithm for finding the SEM poles. In Sect. 4, some numerical results and comparisons with other methods will be described. Section 5 concludes the paper.

2 Method of Modal Parameters for finite wires

We consider a wire of arbitrary geometrical form above a perfect-conducting ground (Fig. 1). The wire is represented by a three-dimensional space curve r(l) parameterized with its natural parameter $l\left(\mathrm{0}\le l\le L/\mathrm{2}$, where L∕2 is the length of the wire). It is assumed that both terminals of the line are connected with the ground. The wire is thin, i.e. the radius of the wire is essentially smaller relative to all other characteristic dimensions of the problem (e.g., wavelength λ, length of the line L∕2). The line is excited by an arbitrary electromagnetic field, which can be external (e.g., an exciting plane wave) or a lumped voltage source (a load also can be consider as a current controlled voltage source). Since the thin-wire approximation is used, it is assumed that: the azimuthal component of the current is neglected; the charge density and the current are concentrated on the wire axis; and the boundary condition for the tangential component of the electric field is satisfied on the surface of the wire. The current I(l) and potential ϕ(l) along the wire are described by the integrodifferential Mixed Potential Integral Equations (MPIE, e.g., Nitsch and Tkachenko, 2007). Usually, the system of MPIE is defined on the “physical” wire in the upper half-plane when the perfectly conducting ground is taking into account by mirroring the corresponding Green's functions. Another way is to add to a physical wire its mirror reflection and to consider the loop closed and completed (Fig. 2). During this procedure, the ground plane is removed. To excite only those modes, as in the initial problem, one has to also to consider a symmetrical external and lumped excitation. Such a procedure yields the MPIE system for the closed loop:

$\begin{array}{}\text{(1a)}& \left\{\begin{array}{l}\frac{\mathrm{d}\phantom{\rule{0.125em}{0ex}}\mathit{\varphi }\left(l\right)}{\mathrm{d}\phantom{\rule{0.125em}{0ex}}l}+j\mathit{\omega }\frac{{\mathit{\mu }}_{\mathrm{0}}}{\mathrm{4}\mathit{\pi }}\underset{\mathit{\left\{}L\mathit{\right\}}}{\oint }{\mathbit{e}}_{l}\left(l\right)\cdot {\mathbit{e}}_{l}\left({l}^{\prime }\right)g\left(l,{l}^{\prime }\right)\phantom{\rule{0.125em}{0ex}}I\left({l}^{\prime }\right)\phantom{\rule{0.125em}{0ex}}\mathrm{d}{l}^{\prime }={E}_{l}^{\mathrm{ex}}\left(l\right)\\ \underset{\mathit{\left\{}L\mathit{\right\}}}{\oint }g\left(l,{l}^{\prime }\right)\phantom{\rule{0.125em}{0ex}}\frac{\mathrm{d}\phantom{\rule{0.125em}{0ex}}I\left({l}^{\prime }\right)}{\mathrm{d}\phantom{\rule{0.125em}{0ex}}{l}^{\prime }}\mathrm{d}{l}^{\prime }+j\mathit{\omega }\phantom{\rule{0.125em}{0ex}}\mathrm{4}\mathit{\pi }{\mathit{\epsilon }}_{\mathrm{0}}\phantom{\rule{0.125em}{0ex}}\mathit{\varphi }\left(l\right)=\mathrm{0}\end{array}\right\;\text{(1b)}& {\mathbit{e}}_{l}\left(l\right):=\partial \mathbit{r}\left(l\right)/\partial l,\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\left|{\mathbit{e}}_{l}\left(l\right)\right|=\mathrm{1};\text{(1c)}& g\left(l,{l}^{\prime }\right)=\frac{{e}^{-jk\sqrt{{\left(\mathbit{r}\left(l\right)-\mathbit{r}\left({l}^{\prime }\right)\right)}^{\mathrm{2}}+{r}_{\mathrm{0}}^{\mathrm{2}}}}}{\sqrt{{\left(\mathbit{r}\left(l\right)-\mathbit{r}\left({l}^{\prime }\right)\right)}^{\mathrm{2}}+{r}_{\mathrm{0}}^{\mathrm{2}}}}\end{array}$

Here L is the length of the complete wire loop, ${E}_{l}^{\mathrm{ex}}\left(l\right)$ is a tangential component of a symmetrized exciting field, el(l) is a tangential unit vector, and $g\left(l,{l}^{\prime }\right)$ is a scalar Green's function defined on the curved wire.

Figure 1Inhomogeneous thin conductor over a conductive ground excited by external field and lumped source.

Figure 2Short-circuited, semi-circular wire.

All functions of the natural parameter l in Eqs. (1a)–(1c) are periodical, i.e., $\mathbit{r}\left(l+L\right)=\mathbit{r}\left(l\right)$ and they can be expanded in a complete orthogonal set of functions exp (−jkml) with ${k}_{m}:=\frac{\mathrm{2}\mathit{\pi }m}{L},\phantom{\rule{0.125em}{0ex}}m=\mathrm{\dots }-\mathrm{1},\mathrm{0},\mathrm{1}\mathrm{\dots }$ They can be represented in the form

$\begin{array}{}\text{(2a)}& {E}_{l}^{\mathrm{ex}}\left(l\right)=\sum _{m=-\mathrm{\infty }}^{\mathrm{\infty }}{E}_{l,m}^{\mathrm{ex}}{e}^{-j{k}_{m}l}=\left[{e}^{-j{k}_{m}l}{\right]}^{T}\cdot {\mathbit{E}}_{l}^{\mathrm{ex}};\text{(2b)}& \mathit{\varphi }\left(l\right)=\sum _{m=-\mathrm{\infty }}^{\mathrm{\infty }}{\mathit{\varphi }}_{m}{e}^{-j{k}_{m}l}=\left[{e}^{-j{k}_{m}l}{\right]}^{T}\cdot \mathbf{\Phi };\text{(2c)}& I\left(l\right)=\sum _{m=-\mathrm{\infty }}^{\mathrm{\infty }}{I}_{m}{e}^{-j{k}_{m}l}=\left[{e}^{-j{k}_{m}l}{\right]}^{T}\cdot \mathbit{I}\end{array}$

Here ${\mathbit{E}}_{l}^{\mathrm{ex}}$, Φ and I are column vectors of the Fourier coefficients of the exciting electric field, potential, and current; and $\left[{e}^{-j{k}_{m}l}\right]:=\left[\mathrm{\dots }{e}^{-j{k}_{-\mathrm{1}}l},\mathrm{1},{e}^{j{k}_{\mathrm{1}}l},\mathrm{\dots }{\right]}^{T}$ is an infinite column vector of basic orthogonal functions.

Then, the MPIE system (Eq. 1a) can be re-written in a modal representation:

$\begin{array}{}\text{(3a)}& \left\{\begin{array}{c}-j\mathbit{k}\cdot \mathbf{\Phi }+j\mathit{\omega }\phantom{\rule{0.125em}{0ex}}{\mathbit{L}}^{\prime }\cdot \mathbit{I}={\mathbit{E}}_{l}^{\mathrm{ex}}\\ -j\mathbit{k}\cdot \mathbit{I}+j\mathit{\omega }\phantom{\rule{0.125em}{0ex}}{\mathbit{C}}^{\prime }\cdot \mathbf{\Phi }=\mathrm{0}\end{array}\right\;\text{(3b)}& \mathbit{k}=\mathrm{diag}\left({k}_{m}\right);\text{(3c)}& {\mathbit{L}}^{\prime }=\frac{{\mathit{\mu }}_{\mathrm{0}}}{\mathrm{2}\mathit{\pi }}{\mathbf{G}}_{L};\text{(3d)}& {\mathbit{C}}^{\prime }=\frac{{\mathit{\mu }}_{\mathrm{0}}}{\mathrm{2}\mathit{\pi }}\mathrm{2}\mathit{\pi }{\mathit{\epsilon }}_{\mathrm{0}}{\mathbf{G}}_{C}^{-\mathrm{1}}\end{array}$

where we introduce infinite square matrixes of inductance per-unit-length L and capacitance per-unit-length C, which can be obtained by expansion of the kernels of the system (Eq. 1a) in a double Fourier series:

$\begin{array}{}\text{(4)}& \begin{array}{rl}{\mathbf{G}}_{L}& =\left[{G}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}}^{L}\right];\\ {G}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}}^{L}& =\underset{\mathit{\left\{}L\mathit{\right\}}}{\oint }\mathrm{d}{l}_{\mathrm{1}}\underset{\mathit{\left\{}L\mathit{\right\}}}{\oint }\mathrm{d}{l}_{\mathrm{2}}{\mathbit{e}}_{l}\left({l}_{\mathrm{1}}\right)\cdot {\mathbit{e}}_{l}\left({l}_{\mathrm{2}}\right)g\left({l}_{\mathrm{1}},{l}_{\mathrm{2}}\right)\\ & \cdot \mathrm{exp}\left(-j{k}_{{m}_{\mathrm{2}}}{l}_{\mathrm{2}}+j{k}_{{m}_{\mathrm{1}}}{l}_{\mathrm{1}}\right)\end{array}\text{(5)}& \begin{array}{rl}{\mathbf{G}}_{C}& =\left[{G}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}}^{C}\right];\\ {G}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}}^{C}& =\underset{\mathit{\left\{}L\mathit{\right\}}}{\oint }\mathrm{d}{l}_{\mathrm{1}}\underset{\mathit{\left\{}L\mathit{\right\}}}{\oint }\mathrm{d}{l}_{\mathrm{2}}g\left({l}_{\mathrm{1}},{l}_{\mathrm{2}}\right)\mathrm{exp}\left(-j{k}_{{m}_{\mathrm{2}}}{l}_{\mathrm{2}}+j{k}_{{m}_{\mathrm{1}}}{l}_{\mathrm{1}}\right)\end{array}\end{array}$

Using the modal representation of MPIE (Eq. 3a) one can easily obtain a formal solution for the column-vectors of the current and potential (Eq. 6), where Z is an infinite modal impedance per-unit-length matrix (Eq. 7a), which defines the scattering field on

$\begin{array}{}\text{(6)}& \left\{\begin{array}{l}\mathbit{I}={\mathbf{Z}}^{-\mathrm{1}}\cdot {\mathbit{E}}_{l}^{\mathrm{ex}}\\ \mathbf{\Phi }=\frac{\mathrm{1}}{j\mathit{\omega }}{\mathbit{C}}^{-\mathrm{1}}\cdot j\mathbit{k}\cdot {\mathbf{Z}}^{-\mathrm{1}}\cdot {\mathbit{E}}_{l}^{\mathrm{ex}}\end{array}\right\;\end{array}$

$\begin{array}{}\text{(7a)}& \begin{array}{rl}{\mathbf{Z}}^{\prime }\left(j\mathit{\omega }\right):& =\frac{\mathbit{k}\cdot {{\mathbit{C}}^{\prime }}^{-\mathrm{1}}\cdot \mathbit{k}}{j\mathit{\omega }}+j\mathit{\omega }{\mathbit{L}}^{\prime }\\ & =\frac{{\mathit{\eta }}_{\mathrm{0}}}{\mathrm{4}\mathit{\pi }\phantom{\rule{0.125em}{0ex}}jk}\left[\mathbit{k}\cdot {\mathbf{G}}_{C}\cdot \mathbit{k}-{k}^{\mathrm{2}}{\mathbf{G}}_{L}\right];\end{array}\text{(7b)}& {\mathbit{E}}_{l}^{\mathrm{sc}}=-{\mathbf{Z}}^{\prime }\cdot \mathbit{I};\text{(7c)}& \stackrel{\mathrm{‾}}{W}=\frac{L}{\mathrm{4}}Re\left\{{\mathbit{I}}^{+}\cdot \mathbf{Z}\cdot \mathbit{I}\right\};\end{array}$

the boundary of the wire (Eq. 7b) and the radiation of the system (Eq. 7c) (Nitsch and Tkachenko, 2007).

The general solution (Eq. 6) is valid for any exciting field, including a plane wave excitation (Krauthauser et al., 2005) or lumped voltage sources (Nitsch and Tkachenko, 2005, 2007; Tkachenko et al., 2011). Moreover, lumped loads can be considered using current-controlled sources with unknown current amplitudes, which can be defined by solution of system of linear equations.

In this paper, however, we restrict our scope to the case of a short-circuited wire with an arbitrary excitation.

For the case of the semi-circular wire (when the closed wire is a circular loop, a structure with high symmetry with constant curvature $\mathbit{K}=\mathrm{1}/R=\mathrm{const}$ and zero torsion T=0), the modal inductance, capacitance and impedance per-unit-length are diagonal matrices and one can obtain a known solution for the induced current in the thin-wire approximation (Wu, 1962)

$\begin{array}{}\text{(8a)}& {L}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}}^{\prime }={L}_{{m}_{\mathrm{1}}}^{\prime }{\mathit{\delta }}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}};\text{(8b)}& {L}_{m}^{\prime }={\mathit{\mu }}_{\mathrm{0}}R\frac{{g}_{m+\mathrm{1}}+{g}_{m-\mathrm{1}}}{\mathrm{4}};\text{(8c)}& {C}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}}={C}_{{m}_{\mathrm{1}}}\cdot {\mathit{\delta }}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}};\text{(8d)}& {C}_{m}=\frac{\mathrm{2}\phantom{\rule{0.125em}{0ex}}{\mathit{\epsilon }}_{\mathrm{0}}}{R{g}_{m}}\end{array}$

$\begin{array}{}\text{(9a)}& {Z}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}}^{\prime }\left(j\mathit{\omega }\right)={Z}_{{m}_{\mathrm{1}}}^{\prime }\left(j\mathit{\omega }\right)\cdot {\mathit{\delta }}_{{m}_{\mathrm{1}},{m}_{\mathrm{2}}};\text{(9b)}& {Z}_{m}^{\prime }\left(j\mathit{\omega }\right)=j\mathit{\omega }\phantom{\rule{0.125em}{0ex}}{L}_{m}+\frac{{k}_{m}^{\mathrm{2}}}{j\mathit{\omega }\phantom{\rule{0.125em}{0ex}}{C}_{m}};\text{(9c)}& {k}_{m}:=\frac{m}{R},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}m=\mathrm{\dots }-\mathrm{1},\mathrm{0},\mathrm{1}\mathrm{\dots }\end{array}$

where

$\begin{array}{}\text{(10)}& {g}_{m}=\underset{\mathrm{0}}{\overset{\mathrm{2}\mathit{\pi }R}{\int }}\phantom{\rule{0.125em}{0ex}}\frac{{e}^{j{k}_{m}l-jk\sqrt{\mathrm{4}{R}^{\mathrm{2}}{\mathrm{sin}}^{\mathrm{2}}\left(l/\mathrm{2}R\right)+{a}^{\mathrm{2}}}}}{\sqrt{\mathrm{4}{R}^{\mathrm{2}}{\mathrm{sin}}^{\mathrm{2}}\left(l/\mathrm{2}R\right)+{r}_{\mathrm{0}}^{\mathrm{2}}}}\phantom{\rule{0.125em}{0ex}}\mathrm{d}l\end{array}$
3 SEM poles and Method of Modal Parameters

According to the definition of the SEM poles they are natural frequencies of the system, i.e., “complex frequencies at which an integral equation representation of the scattering problem has solution with no incident electromagnetic field” (Baum). In the initial papers, “they can be found from the zeros of the determinant of the corresponding moment method matrixes” (Baum). Now, however, the MPIE system is represented in the modal form. Therefore, the SEM poles can also be found from the zeros of the determinant of the impedance per-unit length matrix (Eqs. 11a, b).

$\begin{array}{}\text{(11a)}& \mathrm{det}\left({\mathbf{Z}}^{\prime }\left(j\mathit{\omega }\right)\right)=\mathrm{0}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathrm{or}\text{(11b)}& \mathrm{det}\left[{k}^{\mathrm{2}}\mathbf{1}-{\mathbf{G}}_{L}^{-\mathrm{1}}\left(k\right)\cdot \mathbit{k}\cdot {\mathbf{G}}_{C}\left(k\right)\cdot \mathbit{k}\right]=\mathrm{0}\end{array}$

For the symmetrical circular loop (Eqs. 11a, b) splits and one has separate simple equations (Eq. 11b) for each integer m.

$\begin{array}{}\text{(12a)}& \mathrm{det}\left({Z}_{m}\left(j\mathit{\omega }\right)\right)=\mathrm{0}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathrm{or}\text{(12b)}& {k}^{\mathrm{2}}-\frac{\mathrm{2}\phantom{\rule{0.125em}{0ex}}{k}_{m}^{\mathrm{2}}{g}_{m}\left(k\right)}{{g}_{m+\mathrm{1}}\left(k\right)+{g}_{m-\mathrm{1}}\left(k\right)}=\mathrm{0};\end{array}$

The equation for the poles (Eq. 12b) up to notation coincides with the one used by Blackburn (1976) and Blackburn and Wilton (1978), if one assumes thin-wire approximation. In these works, detailed research of the poles for circular loop, including the case of high layers, was carried out.

It is noteworthy that from Eqs. (6), (7a), and (11a), the SEM poles are independent of the kind of excitation of the system (e.g., distributed excitation, lumped excitation). Below we describe both an approximate analytical method and a numerical algorithm to find the poles.

## 3.1 Analytic approach: perturbation theory for SEM poles

It is possible to show that for the case of infinitely thin wire r0→0, the matrices GL and GC become diagonal and approximately equal to each other (${\mathbf{G}}_{L}\approx {\mathbf{G}}_{C}\approx \mathrm{2}\mathrm{ln}\left(\stackrel{\mathrm{̃}}{L}/{r}_{\mathrm{0}}\right)\cdot \mathbf{1}$, where $\stackrel{\mathrm{̃}}{L}$ is some length parameter of the system, (i.e., for the circle $\stackrel{\mathrm{̃}}{L}\approx \mathrm{2}R\right)$, which leads to ${\mathbf{G}}_{L}^{-\mathrm{1}}\left(k\right)\cdot \mathbit{k}\cdot {\mathbf{G}}_{C}\left(k\right)\cdot \mathbit{k}\to {\mathbit{k}}^{\mathrm{2}}$. Thus, we can obtain the approximate equation for the SEM - poles, which has a simple solution:

$\begin{array}{}\text{(13)}& \begin{array}{rl}{k}^{\mathrm{2}}& \cdot \mathbf{1}-{\mathbit{k}}^{\mathrm{2}}=\mathrm{0}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}⇒\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{k}^{\mathrm{2}}={k}_{m}^{\mathrm{2}}\\ & ⇒\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{k}_{m}^{\left(\mathrm{0}\right)}=±\mathrm{2}\mathit{\pi }m/L,\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}m=\mathrm{\dots }-\mathrm{2},-\mathrm{1},\mathrm{0},\mathrm{1},\mathrm{2}\mathrm{\dots }\end{array}\end{array}$

These values can be used as a zero-order approximation of the iteration (perturbation theory) solution. Any subsequent iteration can be evaluated from the previous iteration as:

$\begin{array}{}\text{(14)}& \begin{array}{rl}\mathrm{det}& \left[\mathrm{diag}\left(\left({k}_{m}^{\left(n+\mathrm{1}\right)}{\right)}^{\mathrm{2}}-{\mathbf{G}}_{L}^{-\mathrm{1}}\left({k}_{m}^{\left(n\right)}\right)\cdot \mathbit{k}\\ & \cdot {\mathbf{G}}_{C}\left({k}_{m}^{\left(n\right)}\right)\cdot \mathbit{k}\right]=\mathrm{0},\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}n=\mathrm{0},\mathrm{1},\mathrm{2},\mathrm{\dots }\end{array}\end{array}$

In particular, for the circle wire, the solution of this equation is:

$\begin{array}{}\text{(15)}& \begin{array}{rl}{k}_{m}^{\left(n+\mathrm{1}\right)}& ={k}_{m}^{\left(n\right)}\sqrt{\frac{\mathrm{2}\phantom{\rule{0.125em}{0ex}}{g}_{m}\left({k}_{m}^{\left(n\right)}\right)}{{g}_{m+\mathrm{1}}\left({k}_{m}^{\left(n\right)}\right)+{g}_{m-\mathrm{1}}\left({k}_{m}^{\left(n\right)}\right)}},\\ & n=\mathrm{0},\mathrm{1},\mathrm{2},\mathrm{\dots };\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}{k}_{m}^{\left(\mathrm{0}\right)}={k}_{m}=m/R\end{array}\end{array}$

Note that, as in the case of the circle wire, the SEM poles in the general case can be classified under the index m, which defines the number of waves of current along the thin wire. The considered iteration solution yields the poles of the so-called first layer, which are nearest to the real axis.

## 3.2 Numerical realization of the Method of Modal Parameters – algorithm of solution equation for the SEM poles

In this sub-section we describe the numerical realization of the method of modal parameters. For the first step, it is necessary to calculate the matrix elements of the capacity-like and inductivity-like G-matrixes (Eqs. 4–5). This is realized by double numerical integration. The first integral is solved numerically using the trapezoidal rule with 10 000 points. An important circumstance here is that the first integration is carrying out with function $g\left(l,{l}^{\prime }\right)$, which has a “singularity” with characteristic dimension r0, when $l\approx {l}^{\prime }$. Therefore, a step of the first integration Δ must describe the singularity of the Green's function: in other words, there must be several points of integration on the radius of the wire (e.g., if the complete length of the loop L is about 20 m, and the radius of the wire is r0=1 cm, we have Δ= 20/10 000 = 2 mm, i.e., $a/\mathrm{\Delta }=\mathrm{5}$). However, the second integration contains a smooth integrand obtained as a result of the first integration, where the characteristic length is about the wavelength or length of non-uniformity of the wire. Tests have shown that 500 intervals are enough for the second integration. Of course, in numerical calculations, we use the finite matrix instead of the infinite matrix. The calculations show that, if we consider a matrix with order ${M}_{max}\left(-{M}_{max},-{M}_{max}+\mathrm{1},\mathrm{\dots }\mathrm{0}\mathrm{\dots }{M}_{max}-\mathrm{1},{M}_{max}$, i.e., with $\left(\mathrm{2}{M}_{max}+\mathrm{1}\right)×\left(\mathrm{2}{M}_{max}+\mathrm{1}\right)$ elements), it issufficient to define the SEM poles for index m up to ${M}_{max}-\mathrm{2}$.

The knowledge of the matrices GL and GC gives a possibility to define a matrix Z. Up to the constant factor this matrix is given by:

$\begin{array}{}\text{(16)}& {\mathbf{Z}}^{\prime }\sim {k}^{\mathrm{2}}\cdot {\mathbf{G}}_{L}\left(k\right)-\mathbit{k}\cdot {\mathbf{G}}_{C}\left(k\right)\cdot \mathbit{k}\end{array}$

Next, the determinant of Z(x) is defined as a function of a normalized complex wave number $x:=k\cdot L/\mathrm{2}\mathit{\pi }$. Then, we look for the solution using the Newton algorithm. As the zero-order approximation for the SEM poles of the first layer for the wiring structure short-circuited at both ends, the next values are used (see Eq. 13)

$\begin{array}{}\text{(17)}& {x}_{m}^{\left(\mathrm{0}\right)}=m\end{array}$

Then we begin an iteration procedure. Because the calculation of the double integrals for all ($\mathrm{2}{M}_{max}+\mathrm{1}×\mathrm{2}{M}_{max}+\mathrm{1}$) matrix elements requires a long time (about 30 min on a HP ENVY 17 Notebook personal computer with processor Intel(R) Core(TM) i7-550U CPU@240 GHz using the Fortran code), we will do this for the function:

$\begin{array}{}\text{(18)}& \begin{array}{rl}F\left(\stackrel{\mathrm{̃}}{x},x\right)& :=\mathrm{det}\left({\stackrel{\mathrm{̃}}{k}}^{\mathrm{2}}\cdot {\mathbf{G}}_{L}\left(k\right)-\mathbit{k}\cdot {\mathbf{G}}_{C}\left(k\right)\cdot \mathbit{k}\right)\\ & \mathrm{where}\phantom{\rule{0.25em}{0ex}}x:=k\cdot L/\mathrm{2}\mathit{\pi },\phantom{\rule{0.25em}{0ex}}\stackrel{\mathrm{̃}}{x}:=\stackrel{\mathrm{̃}}{k}\cdot L/\mathrm{2}\mathit{\pi }\end{array}\end{array}$

The iterations are carried out in two circles. In the internal circle (15 iterations), we use the Newton algorithm only for the value $\stackrel{\mathrm{̃}}{x}$ for a fixed value x

$\begin{array}{}\text{(19)}& \begin{array}{rl}{\stackrel{\mathrm{̃}}{x}}_{m}^{\left(j\right)}& :={\stackrel{\mathrm{̃}}{x}}_{m}^{\left(j-\mathrm{1}\right)}-\frac{F\left({\stackrel{\mathrm{̃}}{x}}_{m}^{\left(j-\mathrm{1}\right)},{x}_{m}^{\left(i\right)}\right)}{{\frac{\partial }{\partial \stackrel{\mathrm{̃}}{x}}F\left(\stackrel{\mathrm{̃}}{x},{x}_{m}^{\left(i\right)}\right)|}_{\stackrel{\mathrm{̃}}{x}={\stackrel{\mathrm{̃}}{x}}_{m}^{\left(j-\mathrm{1}\right)}}},\\ & j=\mathrm{1}\mathrm{\dots }{J}_{max}=\mathrm{15}\end{array}\end{array}$

For the next external circle, we repeat this procedure for the ${x}_{m}^{\left(i+\mathrm{1}\right)}={x}_{m}^{\left({I}_{max}\right)}$, repeating this procedure ${I}_{max}=\mathrm{3}$ times (i=0, 1, 2…3). Both initial values are the same ${x}_{m}^{\left(\mathrm{0}\right)}={\stackrel{\mathrm{̃}}{x}}_{m}^{\left(\mathrm{0}\right)}=m$.

4 Numerical results for the SEM poles

## 4.1 Short-circuited semi-circular wire

To check the method, we consider a vertical, short-circuited, semi-circular wire above a perfect conducting ground (Fig. 2). The parameters of the structure are: R=4 m and r0=1 cm. The matrixes GC(k) and GL(k) are calculated by double integration and partially presented in Figs. 4–6. Using the exact calculations this matrix must be diagonal (see Eqs. 8–10). The relative value of the off-diagonal elements relative to diagonal elements represents the error of the method (about 10−3–10−4). But practically, the matrix Z is diagonal and coincides with the exact result. The determinant of this matrix has minimum approximately for xm1 (see Fig. 3 for the inverse value, which has maximums).

Figure 3Function $\mathrm{1}/|F\left(x,x\right)|$ for the circular wire.

Figure 4The real part of the matrix GL for the circular wire. $k=\mathrm{2.3}/R$.

The comparison of the values of the poles of the first layer obtained using the numerical algorithm described above with the solution of the exact Eq. (12b) obtained by Maple software has shown excellent agreement (Table 1). Note that for corresponding values of the parameters of the circular loop the values of the SEM poles of the first layer obtained by the two methods described above coincide with these one obtained earlier by Umashankar and Wilton (1974).

Figure 5The real part of the matrix GC for the circular wire. $k=\mathrm{2.3}/R$.

Figure 6The real part of the matrix Z for the circular wire. $k=\mathrm{2.3}/R$.

Table 1SEM poles of the first layer for the circular wire.

## 4.2 Short-circuited semi-elliptical wire

As the second example, we consider the semi-elliptical short-circuited wire above a perfectly conducting ground (Fig. 7).

Figure 7Short-circuited semi-elliptical wire.

The parameters of the full ellipse (initial wire plus mirrored wire) are: a=4 m, b=1 m, r0=1 cm.

The matrixes GL and GC are presented in Fig. 9 for real k. For non-symmetrical configurations, these matrices are not diagonal. However, for the smooth wire ($\mathit{\lambda }\mathrm{|}\mathrm{d}R\left(l\right)/\mathrm{d}l\mathrm{|}\ll \mathrm{|}R\left(l\right)\mathrm{|}$, where R(l) is the radius of curvature of the wire) the number of essential non-diagonal elements is not too big. This can be observed in Fig. 9, especially for the real parts. For the imaginary parts, the non-diagonal elements are essential.

Figure 8Function $\mathrm{1}/|F\left(x,x\right)|$ for the elliptical wire.

Figure 9Matrices GL (a, b) and GC (c, d) for semi-elliptical short-circuited wire. $kL/\mathrm{2}\mathit{\pi }=\mathrm{2}/\mathrm{3}$.

Next, we investigate the symmetry properties of the matrices, which can reduce the time needed to make the calculations. If one looks carefully at the numerical data for the matrix GC (see Fig. 9c, d), one can see that:

$\begin{array}{}\text{(20a)}& {\left({\mathbf{G}}_{C}\right)}_{-m\mathrm{1},-m\mathrm{2}}={\left({\mathbf{G}}_{C}\right)}_{m\mathrm{1},m\mathrm{2}};\text{(20b)}& {\left({\mathbf{G}}_{C}\right)}_{m\mathrm{2},m\mathrm{1}}={\left({\mathbf{G}}_{C}\right)}_{m\mathrm{1},m\mathrm{2}}\end{array}$

The same is true for the matrixes GL(k) and Z(k). This can be proven in the general case. From the definition of GC(k) we have

$\begin{array}{}\text{(21)}& \begin{array}{rl}& {\left({\mathbf{G}}_{C}\right)}_{m\mathrm{1},m\mathrm{2}}=\underset{\mathrm{0}}{\overset{L}{\int }}\underset{\mathrm{0}}{\overset{L}{\int }}{e}^{j{k}_{m\mathrm{1}}{l}_{\mathrm{1}}-j{k}_{m\mathrm{2}}{l}_{\mathrm{2}}}{g}_{C}\left({l}_{\mathrm{1}},{l}_{\mathrm{2}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \underset{\begin{array}{l}{l}_{\mathrm{1}}->{l}_{\mathrm{2}}\\ {l}_{\mathrm{2}}->{l}_{\mathrm{1}}\end{array}}{=}\underset{\mathrm{0}}{\overset{L}{\int }}\underset{\mathrm{0}}{\overset{L}{\int }}{e}^{j{k}_{m\mathrm{1}}{l}_{\mathrm{2}}-j{k}_{m\mathrm{2}}{l}_{\mathrm{1}}}{g}_{C}\left({l}_{\mathrm{2}},{l}_{\mathrm{1}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\underset{\mathrm{0}}{\overset{L}{\int }}\underset{\mathrm{0}}{\overset{L}{\int }}{e}^{j{k}_{m\mathrm{1}}{l}_{\mathrm{2}}-j{k}_{m\mathrm{2}}{l}_{\mathrm{1}}}{g}_{C}\left({l}_{\mathrm{1}},{l}_{\mathrm{2}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}={\left({\mathbf{G}}_{C}\right)}_{-m\mathrm{2},-m\mathrm{1}}\end{array}\end{array}$

Here we used the symmetry property of the scalar Green function when rearranging the arguments ${g}_{C}\left({l}_{\mathrm{2}},{l}_{\mathrm{1}},k\right)={g}_{C}\left({l}_{\mathrm{1}},{l}_{\mathrm{2}},k\right)$.

To prove the property (Eq. 20b) one can write:

$\begin{array}{}\text{(22)}& \begin{array}{rl}& {\left({\mathbf{G}}_{C}\right)}_{m\mathrm{1},m\mathrm{2}}:=\underset{\mathrm{0}}{\overset{L}{\int }}\underset{\mathrm{0}}{\overset{L}{\int }}{e}^{j{k}_{m\mathrm{1}}{l}_{\mathrm{1}}-j{k}_{m\mathrm{2}}{l}_{\mathrm{2}}}{g}_{C}\left({l}_{\mathrm{1}},{l}_{\mathrm{2}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \underset{\begin{array}{l}{l}_{\mathrm{1}}->-{l}_{\mathrm{1}}\\ {l}_{\mathrm{2}}->-{l}_{\mathrm{2}}\end{array}}{=}\underset{\mathrm{0}}{\overset{-L}{\int }}\underset{\mathrm{0}}{\overset{-L}{\int }}{e}^{-j{k}_{m\mathrm{1}}{l}_{\mathrm{1}}+j{k}_{m\mathrm{2}}{l}_{\mathrm{2}}}{g}_{C}\left(-{l}_{\mathrm{1}},-{l}_{\mathrm{2}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\underset{\mathrm{0}}{\overset{L}{\int }}\underset{\mathrm{0}}{\overset{L}{\int }}{e}^{-j{k}_{m\mathrm{1}}{l}_{\mathrm{1}}+j{k}_{m\mathrm{2}}{l}_{\mathrm{2}}}{g}_{C}\left({l}_{\mathrm{1}},{l}_{\mathrm{2}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \underset{\begin{array}{l}{l}_{\mathrm{1}}->{l}_{\mathrm{2}}\\ {l}_{\mathrm{2}}->{l}_{\mathrm{1}}\end{array}}{=}\underset{\mathrm{0}}{\overset{L}{\int }}\underset{\mathrm{0}}{\overset{L}{\int }}{e}^{-j{k}_{m\mathrm{1}}{l}_{\mathrm{2}}+j{k}_{m\mathrm{2}}{l}_{\mathrm{1}}}{g}_{C}\left({l}_{\mathrm{2}},{l}_{\mathrm{1}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\underset{\mathrm{0}}{\overset{L}{\int }}\underset{\mathrm{0}}{\overset{L}{\int }}{e}^{-j{k}_{m\mathrm{1}}{l}_{\mathrm{2}}+j{k}_{m\mathrm{2}}{l}_{\mathrm{1}}}{g}_{C}\left({l}_{\mathrm{1}},{l}_{\mathrm{2}},k\right)\mathrm{d}{l}_{\mathrm{1}}\mathrm{d}{l}_{\mathrm{2}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}={\left({\mathbf{G}}_{C}\right)}_{m\mathrm{2},m\mathrm{1}}\end{array}\end{array}$

Here we use also the fact of the even number of integrations and the periodic property of all considered functions.

The results of Eqs. (21) and (22) prove the suggestion Eq. (20).

The matrix Z is also non-diagonal. However, near index ${m}_{\mathrm{eff}}:=kL/\mathrm{2}\mathit{\pi }\approx m\phantom{\rule{0.25em}{0ex}}\left(m=±\mathrm{1},±\mathrm{2},\mathrm{\dots }\right)$ corresponding non-diagonal elements are small. The determinant of this matrix has its minimum for xm (see Fig. 8 for the inverse value, which has its maximums). This value serves as a zero-order iteration in the algorithm described above. Using this algorithm, the SEM poles of the first layer were found (Fig. 10). To check the developed method, we compared the results with the SEM poles obtained from the frequency response function of the line obtained using MoM NEC (Numerical Electromagnetic Code) software. The SEM poles can be extracted from numerical data using a Padé approximation, as described by Senior and Pond (1981). The comparison of the poles obtained by the Method of Modal Parameters and proceeding of the data of NEC calculation is shown in Fig. 10.

Figure 10SEM-poles for semielliptical configuration obtained by modal parameter method and by analysis of NEC data.

Figure 11SEM-poles for semi-circle and semi-elliptical structures. ${L}_{\mathrm{ellips}}={L}_{\mathrm{circle}}=\mathrm{17.1568}$ m.

The comparison of the SEM poles for the semi-elliptical wire with the poles for the semi-circular wire with similar length is shown in Fig. 11. One can see that the imaginary part for the semi-circular wire is essentially larger relative to the semi-elliptical wire. Both sets of poles have a practically linear dependence from the real part. We believe that this is because, in the case of prolate semi-elliptical configuration, the horizontal part of the wire is much closer to the conductive surface, which leads to a weaker radiation, which is ultimately responsible for the appearance of the imaginary part of the SEM poles.

## 4.3 Straight horizontal wire with short-circuited risers

As the third example, we consider another wiring structure, which has many practical applications: the straight horizontal wire with short-circuited risers. We also named this the Vance configuration, because it appears in the book of Edward Vance (1978). The parameters of the line are (Fig. 12): Lhor=10 m, h=0.5 m, r0=1 cm.

Figure 12Straight horizontal wire with short-circuited risers.

Figure 13Function $\mathrm{1}/|F\left(x,x\right)|$ for the Vance configuration.

The matrixes GC and GL for this configuration are presented in Fig. 14. As in the previous example with the prolate elliptic wire, these matrixes contain non-diagonal elements, which, however, are larger than in the previous case, because of the non-smoothness of the line. Nevertheless, the properties of the matrixes (Eqs. 20a, b) are still valid.

Figure 14Matrices GL (a, b) and GC (c, d) for a straight horizontal wire with short-circuited risers. $kL/\mathrm{2}\mathit{\pi }=\mathrm{2}/\mathrm{3}$.

The frequency dependence of the inverse determinant of the matrix Z (which is also non-diagonal) is shown in Fig. 13. As in the previous cases, it has its maximums near ${m}_{\mathrm{eff}}:=kL/\mathrm{2}\mathit{\pi }\approx m\phantom{\rule{0.25em}{0ex}}\left(m=±\mathrm{1},±\mathrm{2},\mathrm{\dots }\right)$, which approximately correspond to the real parts of the SEM poles.

Knowledge of the determinant, as a function of complex wave number k, gives us the possibility to determine the SEM poles using the algorithm described above. To check the modal parameter method we compared the results with the one obtained by two alternative methods: the analysis of the response function obtained by NEC (Senior and Pond, 1981), and the asymptotic method (see Tkachenko et al., 2014; Middelstaedt et al., 2016). The results are presented in Fig. 15.

Analyzing the results, we can draw several conclusions. First, the agreement with the “exact” numerical results, as well as with the approximate analytical method is much better in the case of the elliptic wire. This might be caused by the uniform length division by sub-elements in this case, unlike the case of the ellipsoid. In the last case, the uniform division is applied to the parameter t in the usual parametric representation of the ellipse, which leads to a non-uniform division of the lengths of sections of integration and, as a result, the incorrect integration near the singular point at the given number of sub elements (10 000). However, this deficiency can be easily corrected in future calculations.

Figure 15SEM-poles for the Vance configuration obtained by the modal parameter method, by the asymptotic approach and by analysis of NEC data.

Figure 16Differential (red color) and common (green color) modes for the Vance structure.

Another important observation is related to the excitation of two types of modes: the common mode and differential mode. The differential mode ($m=\sim \mathrm{0}$, 2, 4, …, see Fig. 15) appears for any transmission line structure, which has horizontal components. The common mode ($m=\sim \mathrm{1}$, 3, 5, …, see Fig. 15) appears for the transmission line structure, which has vertical elements and excited, for example, by an external field, which has vertical components. Note that both types of modes appear in the two previous examples: semi-circular wire and semi-elliptical wire above the perfectly conducting ground with a non-symmetrical excitation. In the case of the wire without vertical components (horizontal open-circuit line above ground), or in the case of a symmetrical configuration (all three considered configurations) and excitation (i.e., plane wave with normal incidence or with grazing incidence), we can observe only one mode type (differential modes or common modes).

The third note is connected with the fact that the presented numerical calculations, using the method of modal parameters is quite a time consuming (for example, the calculation of the 7–10 poles requires 5–6 h on the computer used here). It is possible to show that the modal parameter method allows additional simplification in the case of smooth wires when the radius of curvature of the wire axis is substantially greater than the wavelength. In this case one can introduce for each mode m smoothly length-dependent modal inductance per-unit-length and capacitance per-unit length. These matrices are diagonal (Tkachenko and Nitsch, 2005). The equation for the SEM poles corresponding to the mth mode is given by the diagonal matrix element of the modal impedance per-unit length with the modal functions. It is possible to show that this result can also be obtained by the variation method (Myers et al., 2011a, b). These problems will be addressed in future work.

5 Conclusion

The SEM poles of the short-circuited thin wire with arbitrary geometry above a perfect conducting ground were investigated by the method of modal parameters. Unlike the Method of Moments using the localized basic functions, the method of modal parameters uses physical modal basic functions, which are the approximate solution of MPIE for thin wires, which essentially simplifies finding the poles. The SEM poles are complex roots of the modal impedance matrix of the system. The numerical method of solution of the corresponding equation was developed and applied for several examples: semi-circular wire, semi-elliptical wire, and horizontal wire with vertical risers. During this procedure, matrices of modal inductance and capacitance were calculated and their symmetry was investigated. The comparison of the SEM poles with the ones obtained by other analytical and numerical methods yields a good agreement. The causes of the disagreement with the numerical results are analyzed. The main cause is using matrices that are not sufficiently large. In our opinion, it possible to improve results by accelerating calculations using the symmetry properties of modal matrixes (Eqs. 20a, b).

Note, that the MoMP can also be applied for the open-circuited wire if one uses a different modal function (sine and cosine). In the future, the developed method can be generalized for the case of a loaded wire by including the lumped impedances as current controlled voltage sources.

Data availability

Data is available upon request.

Author contributions

SVT was responsible for the development of the method of modal parameters, the idea about application of the MoMP for the calculation of SEM poles, the development of the corresponding analytical method and software, discussion of the results, and writing the text. JBN was responsible for the method of modal parameters, the idea about application of the MoMP for the calculation of SEM poles, carrying out the analytical calculations, discussion of the results, and writing the text. FM was responsible for checking the results both analytically (asymptotic approach and the equations for the reflection coefficients) and numerically (method of moments and the software for the asymptotic approach), discussion of the results, and writing the text. RR had the idea about application of the MoMP for the calculation of SEM poles and contributed to the discussion of the results and writing the text (Introduction and Conclusion). MS had the idea about application of the MoMP for the calculation of SEM poles and contributed to the discussion of the results and writing the text (Introduction and Conclusion). RV had the idea about application of the MoMP for the calculation of SEM poles and contributed to the discussion of the results and writing the whole text.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

Review statement

This paper was edited by Frank Gronwald and reviewed by Andres Gallego and one anonymous referee.

References

Baum, C. E.: On the singularity expansion method for the solution of electromagnetic interaction problems, Interaction Notes, Note 88, available at: http://ece-research.unm.edu/summa/notes/In/0088.pdf (last access: 4 August 2019), 1971.

Baum, C. E., Giri, D. V., and Tesche, F. M. (Eds.): The Singularity Expansion Method in Electromagnetics: A Summary Survey and Open Questions, Lulu Enterprises, North Carolina, USA, 2012.

Blackburn, R. F.: Analysis and synthesis of an impedance-loaded loop antenna using the singularity expansion method, Sensor and Simulation Notes, Note 214, available at: http://ece-research.unm.edu/summa/notes/SSN/note214.pdf (last access: 4 August 2019), 1976.

Blackburn, R. F. and Wilton D. R.: Analysis and synthesis of an impedance-loaded loop antenna using the singularity expansion method, IEEE T. Antenn. Propag., 26, 136–140, https://doi.org/10.1109/TAP.1978.1141787, 1978.

Giri, D. V. and Tesche, F. M.: An of the natural frequencies of a straight wire by various methods, IEEE T. Antenn. Propag., 60, 5859–5866, https://doi.org/10.1109/TAP.2012.2211317, 2012.

Krauthauser, H. G., Nitsch, J., Tkachenko, S., Korovkin, N., and Scheibe, H.-J.: Transfer impedance at high frequencies, in Proceedings of the 2005 International Symposium on Electromagnetic Compatibility, Chicago, IL, USA, 8–12 August 2005, 228–233, https://doi.org/10.1109/ISEMC.2005.1513505, 2005.

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.

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.

Myers, J. M., Sandler, S. S., and Wu, T. T.: Electromagnetic resonances of a straight wire, IEEE T. Antenn. Propag., 59, 129–134, https://doi.org/10.1109/TAP.2010.2090479, 2011a.

Myers, J. M., Sandler, S. S., and Wu, T. T.: Electromagnetic resonances of a straight wire on an earth-air interface, IEEE T. Antenn. Propag., 59, 3154–3164, https://doi.org/10.1109/TAP.2011.2161539, 2011b.

Nitsch, J. B. and Tkachenko, S. V.: Global and modal parameters in the generalized transmission line theory and their physical meaning, URSI Radio Sci. Bull., 312, 21–31, 2005.

Nitsch, J. B. and Tkachenko, S. V.: Propagation of current waves along quasi-periodical thin-wire structures: taking radiation losses into account, URSI Radio Sci. Bull., 322, 19–40, 2007.

Rachidi, F. and Tkachenko, S. (Eds.): Electromagnetic Field Interaction with Transmission Lines: From Classical Theory to HF Radiation Effects, WIT Press, Boston, USA, ISBN 978-1-84564-063-7, 2008.

Singaraju, B. K., Giri, D. V., and Baum, C. E.: Further development in the application of contour integration to the evaluation of the zeros of analytic functions and relevant computer programs, Mathematical Notes, Note 42, available at: http://ece-research.unm.edu/summa/notes/Mathematics/0042.pdf (last access: 4 August 2019), 1976.

Senior, T. B. A. and Pond, J. M.: Pole extraction in the frequency domain, Interaction Notes, Note 411, available at: http://ece-research.unm.edu/summa/notes/In/0411.pdf (last access: 4 August 2019), 1981.

Tesche, F. M., Ianoz, M. V., and Karlsson, T.: EMC Analytical Methods and Computational Models, Willey&Son, NY, USA, ISBN 0-471-15573-X, 1997.

Tkachenko, S. and Nitsch, J.: On the Electromagnetic Field Excitation of Smoothly Curved Wires, in: Proceedings of the VI International Symposium on Electromagnetic Compatibility and Electromagnetic Ecology (EMC-2005), St.-Petersburg, Russia, 21–24 June 2005, 115–120, https://doi.org/10.1109/EMCECO.2005.1513078, 2005.

Tkachenko, S., Rachidi, F., and Ianoz, M.: High-frequency electromagnetic field coupling to long terminated lines, IEEE T. Electromagn. C., 43, 117–129, https://doi.org/10.1109/15.925531, 2001.

Tkachenko, S., Nitsch, J., and Rambousky, R.: Electromagnetic field coupling to transmission lines inside rectangular resonators, Interaction Notes, Note 623, available at: http://www.ece.unm.edu/summa/notes/In/IN623.pdf (last access: 4 August 2019), 2011.

Tkachenko, S. V., Nitsch, J., Vick, R., Rachidi, F., and Poljak, D.: Singularity expansion method (SEM) for long terminated transmission lines, in Proceedings of the 2013 International Conference on Electromagnetics in Advanced Applications (ICEAA), Torino, 9–13 September 2013, 1091–1094, https://doi.org/10.1109/ICEAA.2013.6632411, 2013.

Tkachenko, S., Middelstaedt, F., Nitsch, J., Vick, R., Lugrin, G., and Rachidi, F.: High-frequency electromagnetic field coupling to a long finite line with vertical risers, in Proceedings of the 2014 XXXIth URSI General Assembly and Scientific Symposium (URSI GASS), Beijing, https://doi.org/10.1109/URSIGASS.2014.6929526, 2014.

Tkachenko, S. V., Nitsch, J. B. Middelstaedt, F., Magdowski, M., Hellge-Theune, D., Sheibe, H.-J., Rambousky, R., and Vick, R.: Application of singularity expansion method (SEM) to nonuniform transmission lines, in Proceedings of European Electromagnetics Symposium EUROEM 2016, 11–14 July 2016, London, UK, available at: http://ece-research.unm.edu/summa/notes/AMEREM-EUROEM/EUROEM 2016 Book of Abstracts.pdf (last access: 4 August 2019) (37Mb), 2016.

Umashankar, K. R. and Wilton, D. R.: Transient characterization of circular loop using singularity expansion method, Interaction Notes, Note 259, available at: http://ece-research.unm.edu/summa/notes/In/0259.pdf (last access: 4 August 2019), 1974.

Wu, T. T.: Theory of the thin circular loop antenna, J. Math. Phys., 3, 1301–1304, https://doi.org/10.1063/1.1703875, 1962.

Vance, E.: Coupling to Shielded Cables, Willey&Son, NY, USA, 1978.

The determinant can be a large value in dependence of the order or the matrix. For the ${M}_{max}=\mathrm{7}$ its value can reach 1014–1015