The analysis of the TM electromagnetic scattering from perfectly electrically conducting polygonal cross-section cylinders is successfully carried out by means of an electric field integral equation formulation in the spectral domain and the method of analytical preconditioning which leads to a matrix equation at which Fredholm’s theory can be applied. Hence, the convergence of the discretization scheme is guaranteed. Unfortunately, the matrix coefficients are improper integrals involving oscillating and, in the worst cases, slowly decaying functions. Moreover, the classical analytical asymptotic acceleration technique leads to faster decaying integrands without overcoming the most important problem of their oscillating nature. Thus, the computation time rapidly increases as higher is the accuracy required for the solution. The aim of this paper is to show a new analytical technique for the efficient evaluation of such kind of integrals even when high accuracy is required for the solution.
1. Introduction
Spectral domain formulations are particularly suitable for the analysis of a wide class of electromagnetic problems ranging from the propagation in planar guides and waveguides or the radiation by planar antennas to the scattering from cylindrical structures or planar surfaces involving homogeneous or stratified media, just to give some examples. In general, the obtained integral equation in the spectral domain does not admit a closed form solution; hence, numerical schemes have to be adopted. The fast convergence of such methods is a key point. When dealing with polygonal cross-section cylindrical structures or canonical shape planar surfaces, just for examples, a well-posed matrix operator equation can be obtained by means of the method of analytical preconditioning [1]. It consists of the discretization of the integral equation by means of Galerkin’s method with a suitable set of expansion functions leading to a matrix equation at which Fredholm’s or Steinberg’s theorems can be applied [2, 3]. In the literature, it has been widely shown that this goal can be fully reached by selecting expansion functions reconstructing the physical behaviour of the fields on the involved objects with a closed-form spectral domain counterpart [4–15]. With such a choice, few expansion functions are needed to achieve highly accurate results and the convolution integrals are reduced to algebraic products. However, the obtained matrix coefficients are improper integrals of oscillating and, in the worst cases, slowly decaying functions to be numerically evaluated. The classical analytical asymptotic acceleration technique (CAAAT), consisting of the extraction from the kernels of such kind of integrals of their asymptotic behaviour while the slowly converging integrals of the extracted parts are expressed in closed form, allows us to obtain faster decaying integrands without overcoming the most important problem of their oscillating nature. Consequently, the convergence of the accelerated integrals becomes slower and slower as higher is the accuracy required for the solution.
In order to overcome this problem, a novel technique has been proposed for the analysis of the propagation in single and multiple coupled microstrip lines in planarly layered media [16], the propagation in multilayered single and coplanar coupled striplines [17, 18], the scattering form a rectangular plate in homogeneous medium or buried in a lossy half-space [19, 20], the complex resonances of a rectangular patch in a multi-layered medium [21], the scattering from a tilted strip buried in a lossy half-space at oblique incidence [22], and recently the scattering from a circular plate in homogeneous medium [23]. By means of algebraic manipulations and a suitable integration procedure in the complex plane, the matrix coefficients, which are single/double integrals involving products of Bessel functions of the first kind, are expressed as linear combinations of proper integrals and/or improper integrals of nonoscillating functions which can be quickly evaluated. It is interesting to observe that, despite the same line of reasoning, the procedures developed in the quoted papers are, in general, different from problem to problem.
In the papers [13, 15], the analysis of the electromagnetic scattering from perfectly electrically conducting (PEC) polygonal cross-section cylinders when a TM polarized plane wave orthogonally impinges onto the scatterer surface has been successfully approached. The problem has been formulated in terms of electric field integral equation (EFIE) in the spectral domain and discretized by means of Galerkin’s method. Jacobi polynomials multiplied by their orthogonality weight have been used as expansion functions. By means of a suitable choice of the polynomials’ parameters, the physical behaviour of the surface current density on each side of the polygonal cross-section and even on the adjacent wedges has been reconstructed. Moreover, it has been shown that the spectral domain counterpart of the selected expansion functions can be expressed in closed form in terms of confluent hypergeometric functions of the first kind. Due to the reciprocity theorem, it has been demonstrated that it is always possible to reduce the convolution integrals to algebraic products; i.e., the matrix coefficients are single improper integrals involving products of confluent hypergeometric functions of the first kind and complicated arguments numerically evaluated by means of CAAAT.
The aim of this paper is the introduction of a new analytical technique for the efficient evaluation of such kind of integrals. Algebraic manipulations and a suitable integration procedure in the complex plane allow us to reduce each integral to a linear combination of proper integrals and improper integrals of nonoscillating functions involving confluent hypergeometric functions of the first kind and the second kind. As will be shown in the following, the proposed technique outperforms the CAAAT especially when a high accuracy is required for the solution.
2. Formulation and Solution of the Problem2.1. Background
In Figure 1, a polygonal cross-section PEC cylinder is sketched. A coordinate system x,y,z is introduced such that the z axis coincides with the cylinder axis. The L sides of the polygonal cross-section are numbered clockwise and a local coordinate system xi,yi,z is introduced on the i-th side with the origin at the centre of the side itself in the position x=x-i,y=y-i and the yi axis oriented in the outward direction. The angle ϕi∈ ]−π, π] denotes the orientation of the xi axis with respect to the x axis and 2ai denotes the length of the i-th side. A TM polarized plane wave impinges onto the cylinder surface with an angle ϕ with respect to the x axis and orthogonally with respect to the cylinder axis, namely,(1)E_incx,y=E0z^e-jkxx+kyy,where kx=-kcosϕ and ky=-ksinϕ, k=2π/λ=ωεμ is the wavenumber, λ is the wavelength, ε is the dielectric permittivity and μ the magnetic permeability of the external medium, ω is the angular frequency. With such a choice, only TM solutions can be obtained; i.e., the induced current is longitudinal and the electromagnetic field is invariant along the z axis.
Geometry of the problem.
By means of Meixner’s theory [24], the following edge behaviour can be established for the longitudinal current density on the i-th surface:(2)Jizxi=1-xiaiti+1+xiaiti-J-izxiwith i=1,2,…L, where (3)ti±=ψi±-π2π-ψi±ψi±≤3π21ψi±≥3π2,ψi± being the angles of the wedges at the abscissas xi=±ai and the well-behaved function J-izxi belongs to the weighted Hilbert space Li2-ai,ai with inner product f,g=∫-aiai1-xi/aiti+1+xi/aiti-fxig∗xidxi. The edge behaviour in (2) is strictly related to the asymptotic behaviour of the Fourier transform of the current itself with respect to xi (J~izu). Indeed, by means of Watson’s lemma [25] it can be stated that(4)J~izu~u→+∞J~iz∞u=ηi-e-juaiuti-+1+ηi+ejuaiuti++1,where ηi± are suitable parameters depending on the geometry of the problem and the incident field.
Hence, as shown in [13], the following spectral domain representation for the longitudinal electric field can be obtained by invoking the superposition principle:(5)Ezx,y=∑i=1L∫-∞+∞J~izugiu,x,ydu,where (6)giu,x,y=ωμ2je-juxix,y-yix,yu2-k2u2-k2.
By imposing the total electric field to be vanishing on the cylinder surface, an EFIE is obtained(7)Ezx,yyj=0=-E0e-jkxx+kyyyj=0with xj≤aj and j=1,2,…L.
The obtained integral equation is discretized by means of Galerkin’s method. In order to achieve fast convergence, the function J-izxi is expanded in series of Jacobi polynomials Pnti+,ti-xi/ai which form an orthogonal basis in the weighted Hilbert space Li2-ai,ai. Hence, the following expansion can be established for the longitudinal component of the surface current on the i-th side:(8)Jizxi=∑n=0+∞Jniφnti+,ti-xiai,where (9a)φnα,βxa=1-xaα1+xaβPnα,βx/aξnα,βΠxa,(9b)ξnα,β=∫-aa1-xaα1+xaβPnα,βxa2dx=a2α+β+1Γn+α+1Γn+β+1n!2n+α+β+1Γn+α+β+1are suitable normalization quantities, Π· is the unitary rectangular window, and Γ· denotes the Gamma function [26]. Moreover, the Fourier transform of the selected expansion functions can be expressed in closed-form in terms of confluent hypergeometric functions of the first kind F11·;·;· [26], i.e.,(10a)φ~nα,βau=12π∫-∞∞φnα,βxaejuxdx=ξ-nα,β2juane-juaF11n+β+1;2n+α+β+2;2jua,(10b)ξ-nα,β=a2α+βBn+α+1,n+β+1πn!ξnα,β.
As shown in [15], by means of reciprocity, it is always possible to individuate a representation of the matrix coefficients such that the convolution integrals can be always interpreted as the Fourier transform in the complex plane of the expansion functions. Hence, the general element of the scattering matrix can be written as(11)In,mi,j=∫-∞+∞f^n,mi,juduwhere(12a)f^n,mi,ju=giu,x-j,y-jφ~nti+,ti-aiuφ~mtj+,tj--ajGiju∗,(12b)Giju=ucosφi-φj+ju2-k2sgnyix-j,y-jsinφi-φj.
Remembering the asymptotic behaviour of the confluent hypergeometric function of the first kind [26], (13)F11a;b;z~z→+∞e±jπaz-aΓbΓb-a+ezza-bΓbΓa,where z=u+jv is a complex number, the upper sign has to be chosen when -π/2≤arg(z)<3π/2 and the lower sign in the case -3π/2<arg(z)≤π/2, and it is simple to obtain the following oscillating asymptotic behaviour for the integrand in (11):(14)f^n,mi,ju~u→+∞∑t=01e-βti,ju∑s=01cn,m,s,ti,jejαs,ti,juutis+tjt+3,where cn,m,s,ti,j are suitable parameters, tl0=tl+, with l∈i,j,(15a)αs,ti,j=--1sai+-1tajcosφi-φj-xix-i,y-i,(15b)βti,j=-1tajsgnyix-j,y-jsinφi-φj+yix-j,y-j≥0,where sgn· is the signum function.
It is interesting to observe that(16)minβti,j=-ajsinφi-φj+yix-j,y-jis the distance between the ordinate of the vertex of the j-th side closer to the i-th side in the coordinate system of the i-th side. It is not difficult to understand that when minβti,j~0, e.g., the i-th and the j-th sides are adjacent/almost adjacent or coplanar/almost coplanar sides, the integral in (11) is slowly convergent. In order to speed up the convergence of such kind of integrals, CAAAT can be applied; i.e., it is possible to extract the asymptotic behaviour of their integrands and express the integrals of the extracted contribution in closed form. Unfortunately, such a technique allows us to obtain faster decaying integrands without overcoming the most important problem of their oscillating nature.
2.2. Proposed Solution
A new analytical technique is now introduced in order to speed up the numerical evaluation of the slowly converging scattering matrix coefficients.
In the complex plane z=u+jv, the following relation can be established [26]:(17)F11a;b;z=Γbe±jπaΓb-aUa;b;z+Γbeze∓jπb-aΓaUb-a;b;-z,where U·;·;· is the confluent hypergeometric function of the second kind [26] and the upper sign has to be chosen when -π/2≤arg(z)<3π/2 while the lower sign in the case -3π/2<arg(z)≤π/2.
Therefore, by posing(18)ψ~nα,βz=ξ-nα,βΓ2n+α+β+1Γn+α+1e±jπn+β+12jzne-jzUn+β+1;2n+α+β+2;2jz,where the upper sign has to be chosen when 0≤arg(z)<2π while the lower sign in the case -π<arg(z)≤π, it is possible to write(19)φ~nα,βz=ψ~nα,βz+-1nψ~nβ,α-z.
It is interesting to observe that ψ~nα,βz is a discontinuous function for Rz=0. Moreover, since [26](20)Ua,b,z=πsinπbF11a,b,zΓ1+a-bΓb-z1-bF111+a-b,2-b,zΓaΓ2-b,it is simple to state that ψ~nα,βz~z→0dnα,β/z2n+α+β+1, where dnα,β is a suitable parameter. Then, ψ~nα,βz is singular for z=0 (except when α=β=-1/2 and n=0 simultaneously) and it is a many-valued function with a branch-cut for Iz>0 if α+β is a noninteger number.
Therefore, the general matrix coefficient can be rewritten as(21)In,mi,j=∫-r+rf^n,mi,judu+∑w,s,t=01∫r+∞fn,m,w,s,ti,juduwhere (22)fn,m,w,s,ti,ju=-1ns+mtgi-1wu,x-j,y-jψ~ntis,tis-1-1w+saiuψ~mtjt,tjt-1--1tajGij-1wu∗,it is simple to show that the integration path of the improper integrals in (21) does not intersect the branch-cuts of the integrands, and it has been chosen r>k in order to avoid the possible singularities of fn,m,w,s,ti,jz which can only be located on the real axis for -k≤u≤k (see Figure 2, just for an example).
Integration path.
On the other hand, since [26](23)Ua;b;z~z→+∞z-a,for -3π/2<arg(z)<3π/2, the many-valued function fn,m,w,s,ti,jz can be rewritten on its principal sheet as(24)fn,m,w,s,ti,jz=f-n,m,w,s,ti,jzejαs,ti,j-1wz-βti,jz2-k2,where f-n,m,w,s,ti,jz~z→∞c-n,m,w,s,ti,j/zγn,m,w,s,ti,j with γn,m,w,s,ti,j≥2 and c-n,m,w,s,ti,j is a suitable parameter.
Hence, the path in the complex plane along with the exponential function in (24) does not oscillate and, simultaneously, fn,m,w,s,ti,jz is an integrable function can be obtained by posing (25)αu-βIRz=0αv+βRRz≥0,where the positions α=αs,ti,j-1w and β=βti,j have been introduced in order to simplify the notation, while R· and I· denote the real part and the imaginary part of a complex number, respectively.
For α=β=0, formulas (25) are always verified. For α≠0 and β=0 the integration path defined by (25) is u=0, sgnαv≥0, while for α=0 and β>0 the integration path is given by u≥k, v=0.
Now, let us consider the more general case for α≠0 and β>0. By means of algebraic manipulations and supposing to consider u≥0, (25) can be reduced to(26)u2β2-v2α2-k2α2+β2=0sgnαv≥0,which leads to the following expressions of v and u as a function of u and v, respectively:(27a)v=αu2β2-k2α2+β2for u≥βkα2+β2,(27b)u=βv2α2+k2α2+β2for sgnαv≥0.It is worth noting that, from a numerical point of view, (27a) should be preferred for α/β~0 while (27b) is more suitable for α/β~∞. Hence, for α,β>0, just for an example, a possible alternative integration path is the one in Figure 3.
Alternative integration path.
Now, let us consider the closed contour C sketched in Figure 4(a) nonintersecting any of the branch-cuts of fn,m,w,s,ti,jz.
Integration contours in the complex plane. (a) Integration contour nonintersecting the branch-cuts of the integrand. (b) Integration contour intersecting one of the branch-cuts of the integrand.
By means of Cauchy’s integral theorem and Jordan’s lemma, it is possible to write(28)limR→∞∮Cfn,m,w,s,ti,jzdz=0,and then, the general improper integral at the right-hand side of (21) can be rewritten as (29a)∫r+∞fn,m,w,s,ti,judu=j∫0v-fn,m,w,s,ti,jr+jvdv+∫ru-fn,m,w,s,ti,ju+jv-du+I-n,m,w,s,ti,j(29b)I-n,m,w,s,ti,j=∫u-+∞fn,m,w,s,ti,ju+jαu2β2-k2α2+β21+jα/β2uu2/β2-k2/α2+β2dufor αβ≤1∫v-+∞fn,m,w,s,ti,jβv2α2+k2α2+β2+jvβ/α2vv2/α2+k2/α2+β2+jdvfor αβ≥1where sgnαv->0 and u-=βv-2/α2+k2/α2+β2, which is a linear combination of proper integrals and an improper integral of an asymptotically nonoscillating function.
Unfortunately, the integration contour can intersect one of the branch-cuts of fn,m,w,s,ti,jz along the horizontal/vertical path of the alternative integration path (see Figure 4(b)). In such a case, one of the confluent hypergeometric functions of the second kind of the integrand of the integral in (29b) takes values on its upper/lower secondary sheet. However [26],(30)Ua,b,zej2πM=1-e-j2πMbΓ1-bΓ1+a-bF11a,b,z+e-j2πMbUa,b,zfor -π<argz≤π, where M=-1,0,1 identifies the lower secondary sheet, the principal sheet, and the upper secondary sheet of U·;·;·, respectively. Hence, it is not difficult to conclude that the integrand of the integral in (29b) is an asymptotically oscillating function for M=±1. Such a problem is completely overcome by introducing the following generalization of the representation (19)(31)φ~nα,βz=χ~n,M+,M-α,βz+-1nχ~n,M-,M+β,α-z,being(32)χ~n,M+,M-α,βz=ej2πM-n+β+1ψ~nα,βz+1-ej2πM+n+α+1-1nψ~nβ,α-z,where the parameters M±∈-1,0,+1 take into account the possible intersection of the integration contour in Figure 4 with one of the branch-cuts of fn,m,w,s,ti,jz.
3. Results and Discussion
The aim of this section is to show the efficiency of the presented technique even by means of comparisons with the CAAAT.
The following normalized truncation error is introduced:(33)errN=JN+1-JNJN,where · is the usual Euclidean norm and JM is the vector of all the expansion coefficients of the longitudinal currents on all the sides evaluated with M expansion functions on each side. The simulations are performed on a laptop equipped with an Intel Core 2 Duo CPU T9600 2.8-GHz 3-GB RAM, running Windows XP and the integrals evaluated by means of a Gauss-Legendre quadrature routine. In this way, for the proposed examples, about 70 integrals per second can be numerically evaluated.
First of all, the analysis of the scattering from an irregular triangular cross-section cylinder of sides AB¯=a1+1/3, BC¯=2a/3, CA¯=a2 for a=λ/2,λ,2λ when a TM polarized plane wave impinges with ϕ=π/6 and E0=1V/m is performed. Such a configuration involves only adjacent sides for which, as underlined above, the integrands of the integrals of the coefficients’ matrix are slowly decaying functions. In Figure 5(a), the normalized truncation error by varying the number of expansion functions used on each side of the cylinder cross-section (N) is plotted revealing a very fast convergence in all the examined cases. As a matter of fact, a normalized truncation error of 10-1,10-2,10-3 is achieved for N ranging from 4, 6, 13 to 11, 14, and 17, respectively. It is interesting to note that less than 35 seconds are needed to fill the coefficients’ matrix for N=20. For the sake of completeness, in Figure 5(b) the reconstructed surface current densities for all the examined cases are plotted. In Figure 5(c) the CPU time ratio, i.e., the ratio between the computation time needed to reconstruct the solution as obtained by using the CAAAT with respect to the presented technique, is reported for all the examined cases as a function of N. As clearly shown, the presented technique always outperforms the CAAAT. Moreover, the CPU time ratio quickly increases as higher is N, i.e., as higher is the accuracy required for the solution.
TM scattering from an irregular triangular cross-section cylinder with AB¯=a1+1/3, BC¯=2a/3, CA¯=a2, and for a=λ/2,λ,2λ, ϕ=π/6, and E0=1V/m. (a) Normalized truncation error. (b) Surface current density. (c) CPU time ratio.
In a second example, the scattering from two identical irregular triangular cross-section cylinders of sides AB¯=DE¯=λ1+1/3, BC¯=EF¯=2λ/3, and CA¯=FD¯=λ2 with AB¯ and DE¯ coplanar sides and BD¯=λ/4,λ/2,λ when a TM polarized plane wave impinges with ϕ=π/6 and E0=1V/m is examined. This is an even more troublesome configuration because it involves even coplanar sides for which the integrands of the corresponding integrals of the coefficients’ matrix are slowly decaying functions which oscillate faster and faster as greater is the distance between the coplanar sides (see formula (16)). In Figures 6(a) and 6(b), the normalized truncation error and the surface current density are plotted, respectively. The convergence is very fast and substantially independent of the distance between the cylinders. An accuracy for the solution of 10-1,10-2,10-3 is achieved for N=7,9,14, respectively. It is interesting to note that less than 100 seconds are needed to fill the coefficients’ matrix for N=20. In Figure 6(c), the CPU time ratio is reported as a function of N. In all the examined cases, the presented technique outperforms the CAAAT especially as higher is the distance between the cylinders and as higher is the number of expansion functions used.
TM scattering from two identical irregular triangular cross-section cylinders with AB¯=DE¯=λ1+1/3, BC¯=EF¯=2λ/3, CA¯=FD¯=λ2, AB¯ and DE¯ coplanar sides, and for BD¯=λ/4,λ/2,λ, ϕ=π/6, and E0=1V/m. (a) Normalized truncation error. (b) Surface current density. (c) CPU time ratio.
In the last example, the scattering from an irregular pentagonal cross-section cylinder of sides AB¯=DE¯=EA¯=2λ, BC¯=λ3, and CD¯=λ when a TM polarized plane wave impinges with ϕ=π/6 and E0=1V/m is analysed. As shown in Figure 7(a), the convergence is very fast even in such a case, as N=7,9,11 is enough to achieve an accuracy for the solution of 10-1,10-2,10-3, respectively. It is interesting to note that less than 70 seconds are needed to fill the coefficients’ matrix for N=20. Even if the correctness of the integral formulation and the discretization scheme has been widely discussed by the authors in the papers [13, 15], for the sake of completeness, a comparison with the commercial software CST Microwave Studio (CST-MWS) is shown in Figure 7(b) in terms of surface current density revealing a very good agreement. To conclude, the CPU time ratio is reported in Figure 7(c). As can be seen, the presented technique always outperforms the CAAAT. Moreover, the CPU time ratio increases by increasing N.
TM scattering from an irregular pentagonal cross-section cylinder with AB¯=DE¯=EA¯=2λ, BC¯=λ3, CD¯=λ, and for ϕ=π/6 and E0=1V/m. (a) Normalized truncation error. (b) Surface current density. (c) CPU time ratio.
4. Conclusions
In recent papers, the TM scattering from polygonal cross-section PEC cylinders has been addressed by means of an EFIE formulation in the spectral domain and the method of analytical preconditioning leading to the numerical evaluation of improper integrals of asymptotically oscillating and slowly decaying functions. In this paper, a new analytical technique has been introduced in order to represent such kind of integrals as linear combinations of quickly converging integrals. The proposed numerical results show that such a technique is very effective and outperforms the CAAAT.
Data Availability
The data used to support the findings of this study are obtained by implementing the method detailed in [13, 15] and by means of CST-MWS.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
NosichA. I.Method of analytical regularization in computational photonics2016518142114302-s2.0-8498781120710.1002/2016RS006044KolmogorovA. N.FominS. V.1999Dover, EnglandMR0435771SteinbergS.Meromorphic families of compact operators1968/19693137237910.1007/BF00251419MR02332402-s2.0-34250536952HongoK.Diffraction by a Flanged Parallel‐Plate Waveguide19727109559632-s2.0-001541929510.1029/RS007i010p00955HongoK.IshiiG.Diffraction of an Electromagnetic Plane Wave by a Thick Slit19782634944992-s2.0-001797145610.1109/TAP.1978.1141870EswaranK.On the solutions of a class of dual integral equations occurring in diffraction problems19904291877399427MR105796710.1098/rspa.1990.0066VelievE. I.VeremeyV. V.HashimotoM.IdemenM.TretyakovO. A.Numerical-Analytical Approach for the Solution to the Wave Scattering by Polygonal Cylinders and Flat Strip Structures1993TokyoScience HouseDavisA.ScharsteinR.Electromagnetic Plane Wave Excitation of an Open-Ended, Finite-Length Conducting Cylinder20127230131910.1163/156939393X00354HongoK.SerizawaH.Diffraction of electromagnetic plane wave by a rectangular plate and a rectangular hole in the conducting plate19994761029104110.1109/8.7771282-s2.0-0032643192BliznyukN. Y.NosichA. I.KhizhnyakA. N.Accurate computation of a circular-disk printed antenna axisymmetrically excited by an electric dipole20002532112162-s2.0-003460753810.1002/(SICI)1098-2760(20000505)25:3<211::AID-MOP15>3.0.CO;2-DTsalamengasJ. L.Rapidly converging direct singular integral-equation techniques in the analysis of open microstrip lines on layered substrates20014935555592-s2.0-003527465610.1109/22.910563LosadaV.BoixR. R.MedinaF.Fast and accurate algorithm for the short-pulse electromagnetic scattering from conducting circular plates buried inside a lossy dispersive half-space20034159889972-s2.0-3534891169010.1109/TGRS.2003.810678LucidoM.PanarielloG.SchettinoF.Analysis of the electromagnetic scattering by perfectly conducting convex polygonal cylinders2006544122312312-s2.0-3364570628510.1109/TAP.2006.872662HongoK.NaqviQ. A.Diffraction of electromagnetic wave by disk and circular hole in a perfectly conducting plane20076811315010.2528/PIER060731022-s2.0-33947118119LucidoM.PanarielloG.SchettinoF.Electromagnetic scattering by multiple perfectly conducting arbitrary polygonal cylinders20085624254362-s2.0-4054913454410.1109/TAP.2007.915419LucidoM.A new high-efficient spectral-domain analysis of single and multiple coupled microstrip lines in planarly layered media2012607202520342-s2.0-8486348586410.1109/TMTT.2012.2195025LucidoM.An analytical technique to fast evaluate mutual coupling integrals in spectral domain analysis of multilayered coplanar coupled striplines2012544103510392-s2.0-8485725931310.1002/mop.26674LucidoM.An efficient evaluation of the self-contribution integrals in the spectral-domain analysis of multilayered striplines2013123603632-s2.0-8487584841810.1109/LAWP.2013.2252139ColucciniG.LucidoM.A new high efficient analysis of the scattering by a perfectly conducting rectangular plate20136152615262210.1109/TAP.2012.2237533MR3055894Zbl1370.78100LucidoM.Electromagnetic scattering by a perfectly conducting rectangular plate buried in a lossy half-space20145210636863782-s2.0-8490185400610.1109/TGRS.2013.2296353LucidoM.Complex resonances of a rectangular patch in a multilayered medium: A new accurate and efficient analytical technique201414512313210.2528/PIER140202042-s2.0-84896481081LucidoM.Scattering by a tilted strip buried in a lossy half-space at oblique incidence20143751622-s2.0-8490317090110.2528/PIERM14041507LucidoM.Di MurroF.PanarielloG.Electromagnetic scattering from a zero-thickness PEC disk: A note on the Helmholtz-Galerkin analytically regularizing procedure2017717132-s2.0-85031736272MeixnerJ.The Behavior of Electromagnetic Fields at Edges1972442446MR0068445JonesD. S.1964New York, NY, USAPergamon PressMR0161555AbramowitzM.StegunI. A.1984The NetherlandsVerlag Harri DeutschMR0208797