Capacitance matrix values for the structure for the square conductors with a circular shield (see Figure 6).
Abstract
The research described in this chapter analyses twodimensional potential problems for the multibody systems, transverse electromagnetic wave propagation along multiconductor transmission lines and twodimensional plane wave scattering by various arrays. All conductors may be of arbitrary crosssections; the only restriction on the system geometry is a smooth parameterization. These problems are mathematically modelled by Dirichlet boundary value problems for either the Laplace or the Helmholtz equation, with the classical integral representation of the solutions in the form of singlelayer potential. The analyticalnumerical algorithm presented here is based on the method of analytical regularization. The key idea behind this technique is an analytical transformation of the initial illposed integral equations to a wellconditioned Fredholm second kind matrix equation. The resulting system of infinite linear algebraic equations is effectively solved using the truncation method: the solution of the truncated system converges to the solution of the infinite system with the guaranteed accuracy that only depends on the truncation number and thus may be prespecified. The solution obtained is applied to the accurate analysis of 2D electrostatic and electrodynamicfield problems for multiconductor systems with arbitrary profiled conductors. Examples of some conceptual shielded transmission lines incorporating various configurations of conductors and scattering problems for the arrays of thick strips establish the utility of our method and its reliability in various situations
Keywords
 Scattering
 propagation
 analytical regularization
 Laplace equation
 Helmholtz equation
1. Introduction
In electrostatic and electromagnetic studies of highly elongated cylinders, and ensembles or arrays of such cylinders, it is wellknown [1] that the most important effects can be treated by replacing the threedimensional structure by the corresponding crosssectional twodimensional (2D) profile. Such cylinders are described by
The longstanding interest in the investigation of the electrostatic field in periodic structures continues because of numerous applications. One example is the analysis of the propagation of the transverse electromagnetic (TEM) wave in open and shielded multiconductor transmission lines [3, 4]. When the contour of a conductor coincides with the coordinate surface of one of the coordinate systems in which the Laplace’s equation is separable, the Fourier method (method of separation of variables) is used. More generally, a variety of potential problems have been solved by the conformal mapping method. These results are described in many classical handbooks and monographs. The number of such solved problems is highly restricted. Nowadays, the need for simulation of devices used in practice requires development of more universal methods to tackle problems with objects of various finitewidth shapes. One of such numerous examples is the capacitance calculation for thick electrodes [5] where a physically reasonable meaning of ‘edge capacitance’ arises only because an accurate charge distribution over the whole electrode could not be accurately calculated. Though solutions obtained for single objects may adequately describe the real situation, most practical problems deal with a finite number of objects (say, conductors). Even when a conductor is of canonical shape (circular or elliptic cylinder), the solution of an electrostatic multiconductor problem for an assembly of cylinders of different radii is a very bulky and lengthy procedure. Solving this problem as a classical boundary value problem for Laplace’s equation and enforcing the preassigned boundary conditions at the surface of each conductor, it is necessary to make multiple reexpansions of the eigen functions of the Laplace’s equation in each local coordinate system in terms of that chosen to satisfy the boundary conditions.
Electromagnetic and acoustic problems described by the Helmholtz equation can also be considered in two dimensions. A great variety of publications consider the problem involving infinite gratings. They are often used in antenna applications as polarizers and filters. In [6], a vector diffraction formulation for the analysis of perfectly conducting gratings of finite width and thickness is presented. The grating is assumed to have a finite number of infinitely long arbitrarily shaped rods, and is illuminated by an arbitrary plane wave. Electric and magnetic field integral equations are used to numerically solve the corresponding TM and TE electromagnetic problems. Periodic structures in the millimetre wave range are considered in [7]: this paper studies single and double periodic devices using a semianalytical modematching technique. Diffraction of the TMpolarized Gaussian beam by
It should be noted that the long history of solving the Laplace and Helmholtz equations is marked by the development of many numerical methods which are useful in simulation of practical devices. Such methods include the finite difference technique, extrapolation [14], pointmatching method [15], boundary element method [16], spectralspace domain method [17], finite element method [1820], transverse modal analysis [21] and modematching method [22]. A numerical integral equation approach is used in [23] to explore planewave scattering from a nonplanar surface with a sinusoidal height profile for the case of the magnetic field parallel to the surface ridges (TM polarization). In spite of effectiveness of these methods in many cases and flexibility in geometrical representation of the structures, modelling of ridges still have some substantial drawbacks. Most of the methods require large resources in terms of computational time and storage. Often the solutions obtained with such purely numerical methods need to be verified through comparison to other results: accuracy generally cannot be guaranteed for a greater number of iterations or largerscale computations. This problem becomes pronounced in some topologically complex configurations. In electromagnetics, the corresponding class of numerical solutions is applicable in the low to intermediate frequency range. Resonant systems behaviour cannot be reliably analysed with the numerical techniques (see [24]). Analyticalnumerical methods such as those based on the method of analytical regularization (MAR) are designed to overcome these drawbacks in the resonant regime. A comparative analysis of the MAR and other methods was conducted in [25], and in [26] the distinctive features of each of the discussed methods are clearly described. The abovementioned methods are mostly suited for analysis of a single or very few conductors. In the case of a significant number of conductors with individual profiles, the effectiveness of such purely numerical methods is highly problematic because of the rapidly growing scale of computations.
In order to address these difficulties, we present here a semianalytical approach to the analysis of 2D electrostatic and electrodynamic field problems for multiconductor systems. The problems to be solved are treated as the classical Dirichlet boundary value problems for the Laplace and Helmholtz equations. It is wellknown [27, 28] that solutions to the Laplace and Helmholtz equations can be represented as a singlelayer potential at points exterior to the body of a single conductor with contour
where
This equation may be classified as a first kind Fredholm equation with a singular kernel; it is illposed [29]. Nevertheless, this problem has been tackled by many authors who used direct numerical schemes for solving its discrete analogue in a form of a first kind algebraic equation. Theoretically, any numerical method applied to solve this equation is unable to guarantee uniform convergence, or predetermined computational accuracy.
The only way to avoid these shortcomings is to transform the initial equation into a second kind Fredholm equation, discretization of which guarantees uniform convergence and any predetermined accuracy of the numerical solution depending on truncation number. We employ the MAR, in particular, described in [30, 31]. An accurate solution to wave scattering by a single infinitely long cylinder of arbitrary crosssection by the MAR was obtained in [32]. The details of the algorithm for cylinders of closed arbitrary profile are presented in [26, 31]. In this chapter, we generalize the MAR for a multiconductor potential problem where each body is an arbitrary profiled cylinder.
2. Regularization of the electrostatic problem: MAR
2.1. Problem statement
Consider (
The problem is to find electrostatic potential
with boundary conditions of the potentials
To employ the regularization procedure, all contours
2.2. Problem solution
The main challenge of this problem is that all the conductors are arbitraryshaped and the classical separation of variables method is not applicable here. We use a more general approach based on an integral representation. Using the superposition principle, we seek the solution for the total field potential
where
where
Applying boundary conditions to (5), one can arrive at the coupled system of integral equations for the unknowns
Equation (7) represents a system of first kind Fredholm integral equations that is generally illposed.
The contours of the conductors’ crosssections should be smooth. Thus for the analysis of the rectangular and square conductors, corners should be smoothed. The two most common types of parameterization are by angle and by arc length. Here, we use parameterization by angle. After parameterization of the contours
we obtain the system of
The described approach permits us to consider a broader set of possible boundary conditions than simply a constant, though in the application to be described, a constant is deployed on the RHS of (9).
For the kernels
Now we can determine
The function
where
Now we can redefine function
Using the wellknown Fourier expansion, we can formulate an expression for the singular part of the Green’s function:
As the function
Also the unknown function
After substitution of all expansions into (9), one can arrive at the system of
Using orthogonal properties and completeness of the functions
Following the steps suggested in [33], it can be shown that coefficient matrix in (18) is square summable:
Thus the infinite system (18) is of a second Fredholm kind and can now be effectively solved by a truncation method. The solution of the truncated system monotonically and rapidly converges to the exact solution. The above solution automatically incorporates the reciprocal influence of all charged cylinders, allowing accurate calculation of the line charge densities on the boundaries and then the field potentials at any point of the space between the conductors. Fourier expansions in (18) are calculated numerically as all functions are regular.
3. Regularization of the Dirichlet problem for Helmholtz equation: MAR technique for the N body multiple scattering
3.1. Statement of the problem
In this section, we consider the scattering problem for the structure which consists of
The scattered electromagnetic field
where point
Here we consider incident fields in the form of a plane wave. We focus on a transverse magnetic (TM) wave polarization of the incident field (
The field should also satisfy the Sommerfeld radiation condition:
As 
3.2. Regularized solution
Solutions to the Laplace equation can be represented as a singlelayer potential at points exterior to the body. Using the superposition principle, we seek the solution as the sum of singlelayer potentials contributed by each cylinder:
Here
Applying boundary conditions to (23), we obtain the system of
After parameterization of the contours
The following notation is used in (25):
The kernel of the integral equation (25) contains a singularity only in the terms
so that the regular part of Green's function is
The function
where
We expand the singular part of Green's functions in the same way as in Section 2.2, and perform the double Fourier series expansion for the regular function
After substitution of all expansions into (25), one can arrive at the system of
where the following notations are used:
The infinite systems (31) can be effectively solved by a truncation method. The solution of the truncated system steadily and rapidly converges to the exact solution [34]. There are no limitations on the number of cylinders with arbitrary smooth crosssections.
4. Numerical results
The numerical code was validated by comparing obtained results with known analytical solutions for the coaxial line with a centred inner conductor [35] and the coaxial line with a shifted inner conductor [36]. Results obtained coincide for up to 16 decimal places with the published solutions starting with
As an illustration of the effectiveness of the obtained solution, we calculate the capacitance matrix for the assembly of arbitrary profiled cylinders located inside the grounded shield. There are no limitations on the number of cylinders with arbitrary smooth crosssections. The high efficiency of the code is also the result of employment of the discrete Fast Fourier Transform. This makes filling of the matrix very fast routine procedure. For example, the computation time for a problem with the four inner cylinders and truncation number
Efficiency of the developed method is also illustrated by the behaviour of the normalized truncation error versus truncation number (see [30]) calculated in the maximum norm sense as:
where
Figure 3 shows the condition number behaviour in the same case. The results are quite accurate and stable: for a simple structure like this, the condition number has reached a stable value even for small values of the truncation number.
In these examples, the ellipse is parameterized by the angle as a parameter. Fewer number of points on the sides of a slender ellipse results in decreasing accuracy for smaller
Various shapes of conductor will be considered in this chapter. For all system configurations here and below, the inner conductors’ potentials are set to be 1; the shield is grounded. The profile of each interior conductor is described by the superellipse equation (32), where function
In this equation,
This parameterization is infinitely differentiable which gives us a great advantage in accuracy. To demonstrate this property, comparison of two different parameterizations used in the solution to the Helmholtz equation for a single rectangle with roundedoff corners is presented in Figure 4.
Parameterization 1 stands for a superellipse formula; straight lines with a combination of quarter circles are used for the Parameterization 2. The superellipse parameterization uses
4.1. Electrostatic problems
4.1.1. Multiconductor transmission lines
Here the power of the method is illustrated by the analysis of multiconductor transmission lines. Other possible applications of our method for the problems modelled by the Laplace equation include impedance calculations for the transmission lines with adjustable inner conductor, published in [39] and capacitance calculations for the capacitance microscope [40].
The distribution of the electrostatic field for a conceptual configuration of a shielded threeconductor transmission line is shown in Figure 5.
It is worth noting that apparently sharp edges are in fact not sharp but have a very small radius of curvature at some points due to parameterization of the contours.
Next we present the calculations for the capacitance matrix (Table 1) calculated by formula




Two configurations are examined  a symmetrical one and another obtained by the translation of one conductor, as indicated in Figure 6. In each case, accuracy was ensured by examining error estimates as a function of truncation number as explained above.
4.1.2. Transmission lines with the closely spaced conductors
Another example demonstrating the effectiveness of the developed algorithm is a study of the closely spaced conductors case (Figure 7). In Table 2, capacitance values for a circular conductor are shown depending on a distance between the inner conductor of a radius 0.1 and a shield of radius 1. Truncation numbers are chosen to ensure capacitance values and are stable to four decimal places.



1  64  5.1270 
0.5  128  6.9491 
0.1  512  14.9620 
0.05  1024  21.0534 
The analysis shows that reliable results are obtained when the condition
4.2. Scattering of a plane wave by an array of thick strips
Arrays, which are composed of a finite number of strips, are probably the most common periodic structures. They are employed in various electromagnetic radiating and waveguiding devices. For example, a simple but effective leakywave antenna can be designed by placing a microstrip grating above a ground plane, as first proposed by Honey in the 1950s [41] and then studied by different authors with many variations [42]. In addition, periodic structures in the millimetre wave range with high precision requirements must be planar structures, for fabrication reasons [7]. Also, in other applications the strip grating is often used as a circular polarizer [43]. The list of applications can be continued. Nowadays, a lot of attention is paid to more realistic models of the strip gratings: finite gratings, excited by compact directional sources [8]; gratings with thick strips [44, 45]; special elemental positioning [46, 47], etc.
In this section, we consider scattering of an obliquely incident Epolarized plane wave by a finite array of metallic thick strips which is relevant to the problems examined in the papers mentioned above. The case of the array of circular cylinders, including resonant effects, was considered in [48].
4.2.1. Linear array of horizontal thick strips
We consider the scattering of the Epolarized plane wave obliquely incident the linear array composed of the metallic thick strips, as shown in Figure 8. The elemental thick strip is described by its width
The radar crosssection (RCS) is determined from the scattered field as
The dependence of the RCS (
In our calculations, we fix the number of the strips
Because of the symmetrical (for normal wave incidence) location of the elements in the array, the distribution
The scattered pattern in the direction
The distribution of the scattered field in the farfield zone is shown in Figures 1113.
Figure 13 demonstrates the SP for the 3, 5 and 9element array with the same geometrical parameters, as in Figure 11, but for a different incidence angle,
Next we consider the frequency dependence
It is worth noting that extremely high values of the function
The next graph in Figure 15 illustrates the effect of perturbing the periodicity of the array on the RCS: the central strip is moved towards to the neighbouring strip by a distance
Now let us investigate how the thickness of the strips in an array impacts the RCS. In fact, we will consider a more general problem. Usually the term thick strip refers to the strip with the width
Here we consider normal wave incidence
Surprisingly the RCS is fairly insensitive to substantial thickening of the original horizontal thin strip (
4.2.2. Inclined arrays of thick strips
It was shown previously in this chapter that different positioning of the strips affects the reflection properties. Let us consider another geometry, arranging the strips as a 2D truncated corner reflector (see Figure 20).
Here we only consider normal wave incidence, so we normalize the function
The frequency dependence of the RCS in the range for the array of three truncated corner reflectors is shown in Figure 22. The RCS for the single truncated corner reflector is of a regular oscillatory character for each parameter
5. Conclusion
In this chapter, a rigorous approach to the numerical analysis of the multiconductor problems in electrostatics and multiple wave scattering for metallic cylinders is presented. The problems are treated as a classical Dirichlet boundary value problem for the Laplace and Helmholtz equations. All conductors may be of arbitrary crosssections; the only restriction on the system geometry is a smooth parameterization of the boundaries.
The 2D multiconductor problems for the Laplace and Helmholtz equations are rigorously solved by the MAR. The problem is transformed to a numerical analysis of an infinite system of linear algebraic equations of the second kind. This explains its fast convergence and guaranteed computational accuracy, depending only upon truncation number
References
 1.
Jackson JD. Classical Electrodynamics . New York: Wiley; 1999. p. 808.  2.
Shi F, Ramesh P and Mukherjee S. On the Application of 2D Potential Theory to Electrostatic Simulation. Communications in Numerical Methods in Engineering 1995; 11(8): 691701.  3.
Paul CR. Analysis of Multiconductor Transmission Lines . New York: Wiley; 1994.  4.
Brando Faria JA. Multiconductor Transmission Lines. Encyclopedia of RF and Microwave Engineering ; 2005.  5.
Kamchouchi HE and Zaky AA. A Direct Method for the Calculation of the Edge Capacitance of Thick Electrodes. Journal of Physics D 1975; 8: 13651371.  6.
Moreira FJS and Prata A. VectorDiffraction Analysis of Finite Perfectly Conducting Gratings with Arbitrary Profiles. Journal of Microwaves and Optoelectronics 1999; 1(5): 1113.  7.
Filtz M and Henke H. Planar Millimeter Wave Structures. Electrical Engineering (Archiv fur Elektrotechnik) 2004; 86(2): 147155.  8.
MataMendez O, Avendao J and ChavezRivas F. Rigorous Theory of the Diffraction of Gaussian Beams by Finite Gratings: TM Polarization. Journal of Optical Society of America A 2006; 23(8): 18891896.  9.
Borghi R, Gori F and Santarsiero M. PlaneWave Scattering by a Set of Perfectly Conducting Circular Cylinders in the Presence of a Plane Surface. Journal of Optical Society of America A 1996; 13(12): 24412452.  10.
Argatov II and Sabina FJ. Acoustic Scattering by a Cluster of Small SoundSoft Obstacles. Wave Motion 2010; 47(8): 537551.  11.
Yang J, Li LW, Yasumoto K and Liang CH. TwoDimensional Scattering of a Gaussian Beam by a Periodic Array of Circular Cylinders. IEEE Transactions on Geoscience and Remote Sensing 2005; 43(2): 280285.  12.
Yasumoto K, Toyama H and Kushta T. Accurate Analysis of TwoDimensional Electromagnetic Scattering from Multilayered Periodic Arrays of Circular Cylinders Using Lattice Sums Technique. IEEE Transactions on Antennas and Propagation 2004; 52(10): 26032611.  13.
Lucido M, Panariello G and Schettino F. Electromagnetic Scattering by Multiple Perfectly Conducting Arbitrary Polygon Cylinders. IEEE Transactions on Antennas and Propagation 2008; 56(2): 425436.  14.
Judd SV, Whiteley I, Clowes RJ and Rickard DC. An Analytical Method for Calculating Microstrip Transmission Line Parameters. IEEE Transactions on Microwave Theory and Techniques 1970; 18(2): 7887.  15.
Zhu NH, Qiu W, Pun EYB and Chung PS. QuasiStatic Analysis of Shielded Microstrip Transmission Lines with Thick Electrodes. IEEE Transactions on Microwave Theory and Techniques 1997; 45(2): 288290.  16.
Chang T and Tan C. Analysis of a Shielded Microstrip Line with Finite Metallization Thickness by the Boundary Element Method. IEEE Transactions on Microwave Theory and Techniques 1990; 38(8): 11301132.  17.
Gentili GG and Macchiarella G. QuasiStatic Analysis of Shielded Planar Transmission Lines with Finite Metallization Thickness by a Mixed SpectralSpace Domain Method. IEEE Transactions on Microwave Theory and Techniques 1994; 42(2): 249255.  18.
Khebir A, Kouki AB and Mittra RM. Higher Order Asymptotic Boundary Condition for Finite Element Modeling of TwoDimensional Transmission Line Structures. IEEE Transactions on Microwave Theory and Techniques 1990; 38(10): 14331438.  19.
Slade GW and Webb KJ. Computation of Characteristic Impedance for Multiple Microstrip Transmission Lines Using a Vector Finite Element Method. IEEE Transactions on Microwave Theory and Techniques 1992; 40(1): 3440.  20.
Alam MS, Hirayama K, Hayashi Y and Koshiba M. Analysis of Shielded Microstrip Lines with Arbitrary Metallization Cross Section Using a Vector Finite Element Method. IEEE Transactions on Microwave Theory and Techniques 1994; 42(11): 21122117.  21.
Yee HY and Wu K. Printed Circuit TransmissionLine Characteristic Impedance by Transverse Modal Analysis, IEEE Transactions on Microwave Theory and Techniques 1986; 34(11): 11571163.  22.
Hong IP, Yoon N, Park SK and Park HK. Investigation of MetalPenetrating Depth in Shielded Microstrip Line. Microwave and Optical Technology Letters 1988; 19(6): 396398.  23.
Zaki KA and Neureuther AR. Scattering from a Perfectly Conducting Surface with a Sinusoidal Height Profile: TM Polarization. IEEE Transactions on Antennas and Propagation 1971; 19(6): 747751.  24.
Booker SM and Vinogradov SS. A Two Body Benchmark Problem for Transient and Ultrawideband Scattering. Journal of Electromagnetic Waves and Applications 2002; 16(3): 365384.  25.
Nosich AI. The Method of Analytical Regularization in WaveScattering and Eigenvalue Problems: Foundations and Review of Solutions. IEEE Transactions on Antennas and Propagation 1999; 41(3): 3448.  26.
Smith PD. Recent Advances in Regularization Techniques for Scattering and Diffraction. Radio Science 2007; 42(6): doi: 10.1029/2007RS003703.  27.
Kellogg O. Foundations of Potential Theory . BerlinHeidelbergNew York: SpringerVerlag; 1967.  28.
Colton D and Kress R. Integral Equation Methods in Scattering Theory . New York: John Wiley & Sons; 1983.  29.
Kress R. Linear Integral Equations, ser. Applied Mathematical Sciences . Berlin: SpringerVerlag (Volume 82); 1989.  30.
Vinogradov SS, Smith PD and Vinogradova ED. Canonical Problems in Scattering and Potential Theory, Part 1: Canonical Structures in Potential Theory . Boca Raton: Chapman & Hall/CRC; 2001.  31.
Vinogradov SS, Vinogradova ED, Wilson C, Sharp I and Tuchkin YA. Scattering of an EPolarized Plane Wave by TwoDimensional Airfoils. Electromagnetics 2009; 29(3): 268282.  32.
Vinogradova ED and Smith PD. Diffraction from 2D Arbitrary Structures with Cavities or Edges: A Rigorous Approach: Proceedings of the XXVIIth General Assembly of URSI , URSI GASS, 2002, The Netherlands, Maastricht.  33.
Shestopalov VP, Tuchkin YA, Poedinchuk AE and Sirenko YK. New Methods of Solving Direct and Inverse Problems of Diffraction Theory . Kharkiv: Osnova (In Russian); 1997.  34.
Kantorovich LV and Akilov GP. Functional Analysis in Normed Spaces . New York: Pergamon Press; 1974.  35.
Collin RE. Field Theory of Guided Waves . Oxford: Oxford University Press; 1991.  36.
Yakover YM. General Analysis of Shielded Transmission Lines with Cylindrical Geometry Conductors. 27th European Microwave Conference 1997; 2: 955958.  37.
Malits P, Haridim M and Chattah Y. Characteristic Impedance of a Shielded Transmission Line with NonConcentric ArcShaped Conductors. Electromagnetics 2009; 29: 337352.  38.
Jielis J. A Generic Geometric Transformation That Unifies a Wide Range of Natural and Abstract Shapes. American Journal of Botany 2003; 90(3): 333338.  39.
Safonova G, Vinogradova ED and Tuchkin Y. Accurate Impedance Calculations for Transmission Lines with Adjustable Inner Conductor. Microwave and Optical Technology Letters 2014; 56(9): 20662070.  40.
Safonova G and Vinogradova ED. Rigorous Approach to Calculation of Capacitance Images of Metallic Samples for Use in Capacitance Microscopy. Progress in Electromagnetics Research B 2013; 55: 383399.  41.
Honey R. A FlushMounted LeakyWave Antenna with Predictable Patterns. IEEE Transactions on Antennas and Propagation 1959; 7(4): 320329.  42.
Klohn KL, Horn RE, Jacobs H and Freibergs E. Silicon Waveguide Frequency Scanning Linear Array Antenna. IEEE Transactions on Microwave Theory and Techniques 1978; 26(10): 764773.  43.
Brand GF. The Strip Grating as a Circular Polarizer. American Journal of Physics 2003; 71(5): 452.  44.
Kobayashi K and Miura K. Diffraction of a Plane Wave by a Thick Strip Grating. IEEE Transactions on Antennas and Propagation 1989; 37(4): 459470.  45.
Peng ST and Shiao CM. Scattering of Plane Waves by Metallic Gratings. Microwave Symposium Digest. IEEE MTTS International 2327 May 1994; 2: 879882, doi: 10.1109/MWSYM.1994.335217.  46.
Filipovic D and Dlabac T. Low Frequency TM PlaneWave Scattering from a Twolayer Doublestrip Grating with Equal Gaps. Series: Electrical Energy 2006; 19(3): http://www.doiserbia.nb.rs/img/doi/03533670/2006/035336700603379F.pdf (accessed 9 September 2014).  47.
Shimoda M and Itakura T. Scattering of Electromagnetic Plane Waves by an Inclined Parallel Strip Grating. Electronics and Communications in Japan (Part II: Electronics) 1990; 73(8): 5665.  48.
Safonova G and Vinogradova ED. Accurate Backscattering Studies for DoubleLayered Array of Infinitely Long Cylinders. Microwave and Optical Technology Letters 2013; 55: 21272131.  49.
Kalhor HA. Plane Metallic Gratings of Finite Number of Strips. IEEE Transactions on Antennas and Propagation 1989; 37(3): 406407.