# Line Source Scattering by Buried Perfectly Conducting Circular Cylinders

Authors: F. Frezza; L. Pajewski; C. Ponti; G. Schettini
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

A two-dimensional scattering problem of a line source by a set of perfectly conducting circular cylinders buried in a semi-infinite medium is solved, in both TE and TM polarization. A cylindrical-wave approach is used and applied to both the field emitted by the source and the field scattered by the buried objects. Reflection and transmission of such fields through the planar interface are evaluated making use of the plane-wave spectrum of a cylindrical wave. Numerical results are presented, with checks confirming the validity of the method.


---

## Body

## 1. Introduction

Analytical and numerical techniques, to solve direct scattering by buried objects, offer useful tools in postprocessing of measured data from environmental sensing, detection of land mines, or of buried utilities, inspection of roads’ and buildings’ anomalies. Several methods are developed in the literature to solve such a scattering problem with cylinders buried below a planar interface [1–4] or embedded in a dielectric slab [5–9]. In the mentioned works, the excitation is given by a plane wave, which, propagating in a vacuum half-space, impinges on the interface with the medium hosting the objects. Less attention has been devoted to alternative kinds of excitations, with a more direct link to sources used in practical applications. As far as ground-penetrating radar investigation is considered, in most cases the excitation field is emitted by a dipole antenna [10]. Therefore, modelling such a field with a point or a line source is a more realistic assumption. In the work of D’Yakonov [11], an analytical solution to the scattering of a time-varying electromagnetic field from a point source, by a cylindrical inhomogeneity below the surface of the earth, is developed. It is assumed that the field produced by the source is parallel to the axis of the cylinder and it is constant along it. D’Yakonov’s exact solution has been extended by Ogunade [12] to get numerical results for the case of an embedded circular cylinder excited by a line current. The problem of scattering of a subterranean cylindrical inhomogeneity from a line source is solved by Howard [13] with a modified mode-matching method in conjunction with an integral equation.In [4], the cylindrical wave approach has been used to solve scattering of a plane wave by perfectly conducting cylindrical objects buried in a dielectric half-space. The method, which is based on expansions into cylindrical waves of the fields scattered by the buried cylinders, allows to take into account in an accurate way the scatterers’ geometry. Moreover, making use of the plane-wave spectrum of a cylindrical wave [14], reflection and transmission through the planar interface between air and ground of the scattered fields, can also be evaluated.In this paper, in reply to a need of realistic description of the scattering scenario, the interaction of a line source with buried perfectly conducting cylinders is dealt with. As the geometry of the problem is two-dimensional, the field radiated by the source can be described by a cylindrical wave. Therefore, the CWA presented in [4] is here suitably extended to evaluate reflection in the source half-space, the transmission in the ground, and the interaction with the buried objects. In particular, reflected and transmitted cylindrical waves of zeroth order are defined, which have been numerically evaluated in order to implement the theoretical analysis into a Fortran code.The paper is organized as follows: in Section2, the theoretical analysis is presented in detail, while in Section 3 some hints on the numerical implementation are given, and numerical results both in near- and far-field regions are reported.


## 2. Theoretical Analysis

The geometry of the problem is depicted in Figure1: N perfectly conducting circular cylinders of infinite extension along the y-axis are buried in a dielectric half-space of permittivity ε1=ε0n12, placed below a half-space of permittivity ε0. In the upper medium, a source given by a line of constant electric or magnetic current, of infinite extension and parallel to the planar interface and the cylinders’ axes, is placed. A main reference frame MRF (O,ξ,ζ) with normalized coordinates ξ=k0x and ζ=k0z, where k0=2π/λ0 is the propagation constant in medium 0, is introduced. The source has axis centred in (χL,ηL) in MRF; with (OL,ξL,ζL) standing for the reference frame RFL centered on the line source, the following change of coordinates applies
(1)ξ=ξL+χL,ζ=ζL+ηL,
where χL<0.Figure 1
Geometry of the problem.Moreover, a reference frameRFq centered on the axis of the qth cylinder (q=1,…,N) is employed, both rectangular (O,ξq,ζq) and polar (O,ρq,θq), with ξq=k0xq=ξ-χq,ζq=k0zq=ζ-ηq, and  ρq=k0rq, with (χq,ηq) being the coordinates in MRF of the center of the qth cylinder, with radius αq=k0aq.The problem is solved by means of a scalar functionV(ξ, ζ), which stands for the total field component parallel to the cylinders’ axes, that is, the electric field Ey or the magnetic field Hy, in TM and TE polarization, respectively. Moreover, the field V(ξ, ζ) is decomposed in each medium into a certain number of terms, that is, the excitation field Vi(ξ, ζ) from the line source, and the terms following its interaction with the interface and the cylinders. In particular, the interaction of the field Vi(ξ, ζ) with the planar interface, in absence of the cylinders, gives a reflected field Vr(ξ, ζ) in medium 0, and a transmitted field Vt(ξ, ζ) in medium 1, as shown in Figure 2(a). Taking into account the geometry of the source, the incident field Vi is here described by a cylindrical wave centered on (χL,ηL), and the reflected and transmitted fields are determined by means of the plane-wave spectrum of a cylindrical wave. From the interaction of the transmitted field with the buried cylinders, a field Vs(ξ, ζ)  is scattered in medium 1, and the scattered-reflected and scattered-transmitted fields Vsr(ξ, ζ) and Vst(ξ, ζ) are excited, due to reflection and transmission of Vs(ξ, ζ) through the interface, respectively (Figure 2(b)). Scattered field contributions are expressed as a sum of cylindrical waves, and the plane-wave spectrum is used to cope with reflection and transmission of such waves.Decomposition of the total field: (a) incident, reflected, and transmitted cylindrical waves; (b) scattered fields.
(a)
(b)The line source radiates ay-directed field expressed through a Hankel function of first kind and zeroth order, with argument proportional to the distance from the source to the observation point. If the source is placed as in Figure 1, the field generated in absence of the planar interfaces and the cylinders can be written as [15]
(2)Vi(ξ,ζ)=-V0H0(1)(n0(ξ-χL)2+(ζ-ηL)2),
where n0 is the refraction index of the upper medium. In order to obtain the reflected and transmitted fields Vr(ξ,ζ) and Vt(ξ,ζ), respectively, the cylindrical function CW0=H0(1) can be expressed as a spectrum of plane waves [14]:
(3)CW0(ξL,ζL)=12π∫-∞+∞F0(ξL,n∥)ein∥ζLdn∥,
with n∥ being the parallel component with respect to the z-axis of the generic plane wave of the spectrum.In (3), the spectrum F0 evaluated in ξL>0 is
(4)F0(ξL,n∥)=2eiξL1-n∥21-n∥2.The reflected fieldVr(ξ,ζ)  is proportional to a reflected cylindrical function of zeroth order RW0(ξ,ζ):
(5)Vr(ξ,ζ)=-V0RW0(ξ,ζ),
which is derived from (3), evaluating the reflected plane wave on each plane wave of the spectrum:
(6)RW0(ξ,ζ)=12π∫-∞+∞Γ01(n∥)F0(-ξ-χL,n∥)ein∥(ζ-ηL)dn∥.In a similar way, the transmitted fieldVt(ξ,ζ) is expressed through a transmitted cylindrical function of zeroth order TW0(ξ,ζ):
(7)Vt(ξ,ζ)=-V0TW0(ξ,ζ),
obtained by evaluating the transmission on each plane wave of the spectrum (3):
(8)TW0(ξ,ζ;-χL)=12π∫-∞+∞T01(n∥)F0(-χL,n∥)×ein11-(n∥/n1)2ξein∥(ζ-ηL)dn∥.An expression of the transmitted field in (7) as a function of polar coordinates centred on the pth cylinder’s axis has to be employed, in order to impose boundary conditions on the surface of each cylinder in an easier way. For this purpose, (7) can be written in the following form:
(9)TW0(ξp,ζp;-χL)=12π∫-∞+∞T01(n∥)F0(-χL,n∥)×ein11-(n∥/n1)2χPein∥(ηP-ηL)×ein1(n⊥ξP+n∥ζP)dn∥,
where the exponential ein1(n⊥ξP+n∥ζP)is a transmitted  plane-wave in RFP, which can be expanded into a series of Bessel functions:
(10)ein1(n⊥ξP+n∥ζP)=∑l=-∞+∞ile-ilφtJl(n1ρP)eilθP.In (10), φt=arctan[n∥/1-(n∥/n1)2] is the angle of transmission of the generic plane wave of the spectrum, and the final expression of the transmitted field is
(11)Vt(ξp,ζp)=-V0∑l=-∞+∞ilJl(n1ρP)eilθP×12π∫-∞+∞T01(n∥)F0(-χL,n∥)ein11-(n∥/n1)2χP×ein∥(ηP-ηL)e-ilarctan[n∥/1-(n∥/n1)2]dn∥.The field in (11) can also be written in the following more compact form:
(12)Vt(ξp,ζp)=-V0∑l=-∞+∞ilJl(n1ρP)eilθP×TW0,l(χP,ηP-ηL,-χL).As to the scattered fieldsVs(ξ,ζ), Vsr(ξ,ζ), and Vst(ξ,ζ), the same expressions derived in [4] are used. The scattered field Vs(ξ,ζ) is an expansion into cylindrical functions with unknown coefficients cqm:
(13)Vs(ξ,ζ)=V0∑q=1N∑m=-∞+∞cqmCWm(n1ξq,n1ζq),
with CWm(n1ξq,n1ζq)=Hm(1)(n1ρq)eiℓθq, where Hm(1)(n1ρq) is the first-kind Hankel function of order m. Making use of the addition theorem of Hankel functions, the field in (13) can be written in RFq, as follows:
(14)Vs(ξq,ζq)=V0∑l=-∞∞Jl(n1ρp)eilθp∑q=1N∑m=-∞+∞imcqm×[Hl(1)(n1ρp)Jl(n1ρp)CWm-l(n1ξqp,n1ζqp)(1-δqp)××+Hl(1)(n1ρp)Jl(n1ρp)δqpδlm].To express the scattered fieldsVsr(ξ,ζ) and Vst(ξ,ζ), the cylindrical function CWm is expanded into a Fourier spectrum Fm:
(15)CWm(ξ,ζ)=12π∫-∞+∞Fm(ξ,n∥)ein∥ζdn∥,
and the explicit expression of  the spectrum is given by
(16)Fm(ξ,n∥)=2ei|ξ|1-n∥21-n∥2{e-imarccosn∥,ξ≥0,eimarccosn∥,ξ≤0.The scattered-reflected fieldVsr(ξ,ζ) is the sum of the fields scattered by each cylinder and reflected by the interface. In particular, reflected cylindrical waves are defined as basis functions of the scattered-reflected field, which are obtained by evaluating the reflection on each plane wave of the spectrum (16):
(17)RWm(n1ξ,n1ζ)=12π∫-∞+∞Γ10(n∥)Fm(n1ξ,n∥)ein1n∥ζdn∥.Thus, the scattered-reflected field is(18)Vsr⁡(ξ,ζ)=V0∑q=1N∑m=-∞+∞cqmRWm[n1(-2χq-ξq),n1ζq].An expression of the scattered field in medium 1, which is a function of polar coordinates centred on thepth cylinder, is used when boundary conditions are imposed:
(19)Vsr⁡(ξ,ζ)=V0∑l=-∞+∞Jl(n1ρp)eilθp∑q=1N∑m=-∞+∞cqmRWm+l×[-n1(χq+χp),n1(ηq-ηp)].Equation (19) is obtained from (18), making use of (16) and (17), and of the expansion of a plane wave into Bessel functions.Transmitted cylindrical waves are defined as basis functions of the scattered-transmitted field:(20)TWm(ξ,ζ,χ)=12π∫-∞+∞T10(n∥)Fm(-n1χ,n∥)×e-i1-(n1n∥)2(ξ+χ)ein1n∥ζdn∥,
which is defined as follows:
(21)Vst(ξ,ζ)=V0∑q=1N∑m=-∞+∞cqmTWm(ξq,ζq,χq).Under the hypothesis of perfectly conducting cylinders, boundary conditions of zero electric field on the cylinders’ surface are imposed, which are(22)Vt+Vs+Vsr⁡|ρp=αp=0,
in TM polarization, and
(23)∂∂ρp(Vt+Vs+Vsr⁡)|ρp=αp=0,
in TE polarization.From boundary conditions (22) and (23), a linear system in the unknown expansion coefficients cqm is derived:
(24)∑q=1N∑m=-∞+∞Aqplm(TM,TE)cqm=Bpl(TM,TE)p=1,…,N,l=0,±1,…,±∞,
with
(25)Aqplm(TM,TE)=i-lGl(TM,TE)(n1ρp)×{δqpδlmGl(TM,TE)(n1ρp)CWm-l(n1ξqp,n1ζqp)(1-δqp)××+RWm+l[-n1(χq+χp),n1(ηq-ηp)]××+δqpδlmGl(TM,TE)(n1ρp)},Bpl(TM,TE)=-Gl(TM,TE)(n1ρp)TW0,l(χp,ηp-ηL,-χL),
where Gℓ(TM)(·)=Jℓ(·)/Hℓ(1)(·) and Gℓ(TE)(·)=Jℓ′(·)/Hℓ′(1)(·).


## 3. Numerical Results

The numerical implementation of the analytical method presented in Section2 has been performed according to criteria of accuracy. As to solve system (24), a truncation on the cylindrical orders given by m and ℓ to a finite number of terms is needed, and the rule Mt=⌊3n1α⌋ [16], with α being the radius of the largest cylinder, has been employed, with a good compromise between accuracy and computational heaviness. A numerical evaluation of Cylindrical Functions relevant to the spectral integrals (6), (8), (17), and (20) is also needed. Suitable algorithms have been developed in [17] for the evaluation of Reflected Cylindrical Functions RWm defined in (18), and they have been generalized in [18] to the evaluation of Transmitted Cylindrical Functions TWm reported in (20). The developed algorithms take into account the infinite extension of the integration domain, and they detect where the spectrum of each wave function is oscillating and evanescent, applying a different integration procedure. On oscillating integrands, a further decomposition of the integration domain is applied making use of adaptive techniques, and Gaussian quadrature rules are employed. As to the evanescent part of the spectrum, it is decomposed into a suitable number of subintervals, where fixed low-order Gauss-Legendre rules are applied.The Reflected Cylindrical Function of zeroth orderRW0 in (6) can be easily evaluated applying the algorithm developed in [17], while a modified integration algorithm has been developed for the evaluation of the Transmitted Cylindrical Functions of zeroth order TW0 in (8). With the transmitted field expressed in a cylindrical frame centred on the axis of the pth cylinder, as in (11), the TW0 is expressed as follows:
(26)TW0(χp,ηp-ηL;-χL)=12π∫-∞+∞T01(n∥)F0(-χL,n∥)×ein11-(n∥/n1)2χPein∥(ηP-ηL)×e-ilarctan[n∥/1-(n∥/n1)2]dn∥.Comparing (26) to (20), the additional term e-iℓarctan[n||/1-(n∥/n1)2] appears in (26). Therefore, a different decomposition and evaluation have been applied, following the general rules of the algorithm in [18].Some checks of convergence have been performed on the implemented numerical code, as the one reported in Figure3, where the scattered far field from a plane-wave excitation is compared to the one from a line source, which is placed at an increasing distance from the interface. A perfectly conducting cylinder of normalized radius α=1, centered in χ=2.57 and η=0, has been simulated. The refractive index of the lower medium is n1=2, and the polarization state is TM. The line source has been centred on ηL=0, and placed at three different distances χL = −0.5, −2, −5 from the interface. A good convergence to the scattered field by a plane-wave excitation can be observed, as the cylindrical wave front is moved away from the interface. In particular, in the cases χL=-2 and χL=-5 the curves are superimposed to the plane wave results. Accuracy of results is shown in Figure 4, for the case of scattered field due to a line source placed in χL=-2. Truncation on cylindrical functions is applied according to the rule Mt=⌊3n1α⌋, and with the given parameters it yields Mt=6. Scattered-transmitted field is then evaluated increasing the order Mt to Mt+1=7 and Mt+2=8. In Figure 4, the absolute value of the difference between scattered-transmitted fields evaluated for two consecutive value of Mt is displayed, confirming the accuracy of the rule Mt=⌊3n1α⌋.Figure 3
Scattered far-field diagram of a perfectly conducting cylinder (α=1, χ=2.57, η=0, n1=2, TM polarization) by a line source centred on ηL=0, and placed at different positions χL; results are compared to the scattered field by a plane-wave excitation.Figure 4
Check of accuracy on results in Figure3, in the case of line source in χL=-2, increasing the truncation order on cylindrical functions Mt.The check of convergence to plane-wave excitation is proposed also in the near field, in Figure5. The simulated parameters are the same as in Figure 3, and the scattered field is now observed along a line parallel to the interface, placed at the normalized distance ξ=-1. A good convergence of the curves with χL=-2 and χL=-5 to the plane wave one can be observed.Figure 5
Scattered near-field diagram of a perfectly conducting cylinder (α=1, χ=2.57, η=0, n1=2, TM polarization) by a line source centred on ηL=0 and placed at different positions χL. The field is evaluated at ξ=-1, and results are compared to the scattered field by a plane-wave excitation.A case of interaction between two identical cylinders of radiusα1=α2=1, buried at a symmetric distance of η1,2=±1.5 from the ξ-axis, at the depth χ1=χ2=2.57, in a medium with refraction index n1=2, is also reported. The scattered near-field diagram, evaluated at the distance ξ=-1, is given in Figure 6, in the TM polarization state, with a source centred in ηL=0 and χL=-0.5.Figure 6
Scattered near-field diagram of two perfectly conducting cylinders (α1=α2=1, χ1=χ2=2.57, η1=1.5,η2=-1.5, n1=2, TM polarization) by a line source centred on ηL=0 and χL=-0.5. Results are evaluated at the distance ξ=-1.


## 4. Conclusions

A theoretical method to solve the two-dimensional scattering of a line source by a set of perfectly conducting cylinders buried in a dielectric half-space has been presented. The CWA has been used to express both the source radiation and its interaction with the interface as well as the scattered fields by the cylinders. The method can deal with TM or TE polarization states and gives solutions both in the near and far field. Results of its numerical implementation have been presented, which have been obtained through an accurate evaluation of spectral integrals.As future work, the presented analysis, carried out in the spectral domain, will be extended to a time-domain solution following an approach analogous to the one developed in [19] for pulsed plane-wave excitation. A possible extension of the method may also consider a lossy dielectric half-space, extending the approach developed in [20, 21] for monochromatic plane-wave excitation.  The scattered field would be represented as an expansion of Hankel’s functions of complex argument. As to Reflected and Transmitted Cylindrical Functions, used both to express reflection and transmission of the incident wave and as basis functions of the scattered-reflected and scattered-transmitted fields, the general definitions presented in the lossless case are still applied, apart from replacing the real refraction index of the lower medium with a complex one. The algorithm developed in [17] for the evaluation of the spectral integrals in the lossless case may still be applied, although the complex reflection index would lead to a different behaviour of the integrand both on the evanescent and homogeneous parts of the spectrum.



---

*Source: 261818-2012-08-07.xml*