Next Article in Journal
Enhanced Pest Recognition Using Multi-Task Deep Learning with the Discriminative Attention Multi-Network
Previous Article in Journal
AISMSNet: Advanced Image Splicing Manipulation Identification Based on Siamese Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Helmholtz–Galerkin Technique in Dipole Field Scattering from Buried Zero-Thickness Perfectly Electrically Conducting Disk

by
Mario Lucido
1,*,
Giovanni Andrea Casula
2,
Gaetano Chirico
1,
Marco Donald Migliore
1,
Daniele Pinchera
1 and
Fulvio Schettino
1
1
Department of Electrical and Information Engineering, University of Cassino and Southern Lazio, 03043 Cassino, Italy
2
Department of Electrical and Electronic Engineering, University of Cagliari, 09124 Cagliari, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(13), 5544; https://doi.org/10.3390/app14135544
Submission received: 14 May 2024 / Revised: 24 June 2024 / Accepted: 24 June 2024 / Published: 26 June 2024
(This article belongs to the Section Electrical, Electronics and Communications Engineering)

Abstract

:
Non-invasive concealed object detection, identification, and discrimination have been of interest to the research community for decades due to the needs to preserve infrastructures and artifacts, guarantee safe conditions for the detection and location of landmines, etc. A modern approach is based on the use of an unmanned aerial vehicle equipped with ground-penetrating radar, which has the advantage of not requiring direct contact with the ground. Moreover, high-resolution underground images are obtained by coherently combining measurements by using a synthetic aperture radar algorithm. Due to the complexity of the real scenario, numerical analyses have always been welcomed to provide almost real-time information to make the best use of the potential of such kinds of techniques. This paper proposes an analysis of the scattering from a zero-thickness perfectly electrically conducting disk buried in a lossy half-space surrounded by air and illuminated by a field generated by a Hertzian dipole located in the air. It is carried out by means of a generalized form of the analytically regularizing Helmholtz–Galerkin technique, introduced and successfully applied by the authors to analyze the plane-wave scattering from a disk and a holed plane in a homogeneous medium. As clearly shown in the numerical results, the proposed method is very effective and drastically outperforms the commercial software CST Microwave Studio 2023.

1. Introduction

Non-invasive concealed object detection, identification, and discrimination have been of interest to the research community for decades due to the need for preserving infrastructures or ancient ruins and artifacts from invasive excavations, with a further advantage in terms of time and costs reduction; or, in war or post-war scenarios, to guarantee safety conditions for the detection and location of landmines in order to avoid devastation and human tragedy caused by explosions due to accidental contact, to name some examples. Among the wide number of approaches developed, those that deserve to be cited are electromagnetic induction, nuclear quadrupole resonance, thermal imaging, and, above all, the techniques based on ground penetrating radar (GPR), which has been widely used in landmine detection (see [1] and the references therein for an overview).
A modern and very effective means to address this problem is based on the use of an unmanned aerial vehicle, i.e., a drone, equipped with GPR [2,3,4,5], which, unlike classic techniques based on a hand-operated metal detector or GPR mounted on a dragged trolley, has the advantage of not requiring direct contact with the ground. Moreover, high-resolution underground images can be obtained by coherently combining the measurements by means of a synthetic aperture radar algorithm [6], and it can be easily conceived that the adoption of modern array configurations for the transmitting and receiving antennas could further improve system performance [7].
Due to the complexity of real scenarios, suitable modeling of the problem is needed to perform numerical simulations that could efficiently provide useful information to best utilize the potential of such kinds of non-invasive techniques. Despite the huge number of techniques based on the finite element method [8], the finite-difference time-domain method [9,10], and the conjugate gradient method [11,12], to cite the most important ones, when dealing with finite scatterers in unbounded media, integral equation formulations are largely preferred because the unknowns and the boundary conditions from which the equations originate are defined on finite supports [13,14,15,16].
Since no closed form solution is available for the obtained integral equation, discretization of the integral operator is needed. However, due to the singular nature of the integral equation at hand, the existence of the solution cannot be generally stated, and even more, if such a solution exists, the convergence of a general discretization scheme cannot be demonstrated. Several techniques, collectively called methods of analytical regularization [17], have been developed to overcome these problems. The rationale is to convert the integral equation at hand to a Fredholm second-kind equation by suitably inverting the most singular part of the integral operator. Of course, based on the problem at hand, the integral operator to be inverted can be conveniently identified with the static part, the high-frequency contribution, a suitable canonical problem, and so on. The well-behaved nature of the equation thus obtained is demonstrated by the Fredholm theory, which allows to establish the existence of a solution from its uniqueness and, moreover, the asymptotic convergence of a discretization scheme preserving the nature of the equation to the exact solution of the problem [18]. In principle, there are no reasons to think of regularization and discretization as two consecutive steps. As a matter of fact, they can proceed “hand to hand” by adopting, for example, the Nystrom discretization technique [19,20,21], which leads to guaranteed-convergence numerical codes founded on the theorems of the quadrature approximations of singular integrals; or the method of analytical preconditioning (MAP) [22], in which Galerkin projection acts as a perfect preconditioner in case a set of orthonormal eigenfunction of the most singular part of the integral operator is selected as expansion basis, thus leading to a Fredholm second-kind matrix equation.
The aim of this paper is to present, for the first time, an analysis of the field scattered from a zero-thickness perfectly electrically conducting (PEC) disk buried in a lossy half-space, surrounded by air, and illuminated by the field generated by a Hertzian dipole located in the air. It is conducted by means of a suitable generalization of the analytically regularizing Helmholtz–Galerkin technique, previously introduced and successfully applied by the authors to the analysis of the plane-wave scattering from a disk [23] and a holed plane [24] in a homogeneous medium. The problem is formulated as an electric field integral equation (EFIE) in the Fourier transform domain for the unknown spectral counterpart of the surface current density on the disk, obtained by imposing vanishing of the tangential to the disk component of the total electric field, i.e., the superposition of the electric field radiated by the dipole and the one scattered from the disk in a planar layered medium. Due to the revolution symmetry of the problem, all the involved functions can be conveniently expanded in a Fourier series and the EFIE recast as an integral equation in the vector Hankel transform (VHT) domain [25]. As a result, the scattered field can be rewritten as the superposition of a free-space contribution, frequently called a “primary” contribution because it contains the most singular part of the integral operator, and a scattered or “secondary” contribution related to the layering of the medium. Since the VHT transforms of the surface curl-free contribution and the surface divergence-free contribution of the general azimuthal harmonic of the surface current density are scalar functions, such contributions are conveniently selected as new unknowns. MAP is used to discretize the integral equation, thus leading to a Fredholm second-kind matrix equation. It is worth noting that the selected basis functions have some additional properties that deserve to be mentioned. First, the guaranteed convergence is fast because the considered basis functions reconstruct the physical behavior of the surface current density at the edge and around the center of the disk. Second, the convolution integrals resulting from the Galerkin projection reduce to algebraic products since the basis functions admit closed-form spectral domain counterparts, i.e., the matrix coefficients are 1-D improper integrals that, in the worst cases, can be efficiently evaluated by means of an analytical technique developed by the authors in a previous paper [26].
This paper is organized as follows: After this Introduction Section, the formulation of the problem and the proposed solution are presented in Section 2. Section 3 is devoted to showing the fast convergence of the proposed method, to the validation of the developed in-house C++ software code by means of comparisons with CST Microwave Studio 2023 (CST-MWS), and to the reconstruction of the scattered field behavior. The conclusions are summarized in Section 4.

2. Formulation of the Problem and Proposed Solution

2.1. Formulation of the Problem

Figure 1 shows the geometry of the problem: a zero-thickness PEC disk of radius a buried at depth d in a homogeneous and isotropic lossy half-space, called medium 2, surrounded by a half-space filled with air, i.e., medium 1. The symbols ε i and μ i denote, respectively, the dielectric permittivity and the magnetic permeability of the i-th medium. A Cartesian coordinate system is introduced such that the xy plane and the z axis coincide, respectively, with the flat interface and the disk axis. An arbitrarily oriented Hertzian dipole located in air, of current density J _ i n c r _ = J _ 0 δ r _ r _ 0 , where r _ = x , y , z and r _ 0 = x 0 , y 0 , z 0 are the vector positions of the observation point and the dipole, respectively, and δ · is the Dirac delta function, radiates an electromagnetic field, E _ i n c r _ , H _ i n c r _ , thus exciting a surface current density on the disk, J _ x , y , that in turn generates a scattered field, E _ s c r _ , H _ s c r _ . The induced surface current density is related to the jump across the disk of the tangential component of the scattered magnetic field, whereas the incident and scattered electric fields are related to each other by the boundary conditions on the disk surface, i.e.,
J _ x , y = z ^ × H _ s c x , y , d + H _ s c x , y , d ,
z ^ × E _ i n c x , y , d + E _ s c x , y , d = 0 ,  
for x 2 + y 2 a 2 .
Hence, a boundary value problem for the Maxwell equation is established. By further imposing the Silver–Muller radiation condition and the power-boundedness condition, the uniqueness of the solution, if it exists, is guaranteed.
Since both the incident and scattered electric fields can be represented by the convolution integral involving the dyadic Green’s function of the problem at hand and the corresponding current density, the Silver–Muller radiation condition is automatically satisfied by the proper selection of Green’s function. Moreover, because Green’s function of a planar layered medium can be expressed in closed form only in the spectral domain [27,28], the incident and scattered electric fields can be conveniently retrieved by means of the inverse Fourier transform of the corresponding spectral domain counterparts.
In general, the electric field generated in medium i = 1,2 by a current density of the kind J ¯ _ j r _ = J _ j x , y δ z z ¯ j located in medium j = 1,2 can be written as the superposition of a primary or free-space contribution, arising only for i = j , and a secondary or scattered contribution, related to the layering of the medium [29], i.e.,
E _ i , j r _ = E _ i , j P r _ + E _ i , j S r _ ,  
E _ i , j P r _ = + G ~ _ _ i , j P u , v J ~ _ j u , v e j u x + v y + R j u , v z z j d u d v i = j 0 i j ,  
E _ i , j S r _ = + G ~ _ _ i , j S u , v J ~ _ j u , v e j u x + v y + R i u , v z + R j u , v z j d u d v ,
where the following definition of a double Fourier transform has been used:
f ~ u , v = 1 4 π 2 + f x , y e j u x + v y d x d y f x , y = + f ~ u , v e j u x + v y d u d v ,
and
G ~ _ _ i , j T u , v = v 2 A i , j T u , v + u 2 B i , j T u , v u v B i , j T u , v A i , j T u , v u C i , j T u , v u v B i , j T u , v A i , j T u , v u 2 A i , j T u , v + v 2 B i , j T u , v v C i , j T u , v u C ¯ i , j T u , v v C ¯ i , j T u , v D i , j T u , v ,  
A i , j T u , v = ω μ j 2 R j u , v u 2 + v 2 1 T = P T i , j μ u , v T = S ,  
B i , j T u , v = R i u , v 2 ω ε i u 2 + v 2 1 T = P T i , j ε u , v T = S ,
C i , j T u , v = R i u , v 2 ω ε i R j u , v s g n z z j T = P 1 j T i , j ε u , v T = S ,  
C ¯ i , j T u , v = 1 2 ω ε i s g n z z j T = P 1 j T i , j ε u , v T = S ,  
D i , j T u , v = u 2 + v 2 2 ω ε i R j u , v 1 + j 2 R i u , v u 2 + v 2 δ z z j T = P 1 j i T i , j ε u , v T = S ,  
T i , j η u , v = 2 η 3 j R j u , v η 1 R 2 u , v + η 2 R 1 u , v δ i , j ,
R i u , v = k i 2 u 2 v 2 ,  
where s g n · is the signum function, δ m , n is the Kronecker delta, k i = ω ε i μ i is the wavenumber of the i-th medium, and ω is the angular frequency.
By using Equation (2) in Equation (1b), the boundary value problem at hand is recast as an integral equation for the Fourier transform of the surface current density on the disk, i.e.,
+ G ~ _ _ 2,2 P u , v + G ~ _ _ 2,2 S u , v e j 2 R 2 u , v d J ~ _ u , v e j u x + v y d u d v = 1 4 π 2 + G ~ _ _ 2,1 S u , v e j u x x 0 + v y y 0 + R 2 u , v d + R 1 u , v z 0 d u d v J _ 0 ,  
for x 2 + y 2 a 2 , where the dyad G ~ _ _ i , j T · is made up of the first two rows of G ~ _ _ i , j T u , v .
The uniqueness of the solution of (5), if it exists, can be stated by the proper selection of the functional space in which the unknown surface current density has to be searched for. For this purpose, it is worth remembering that, in order to satisfy the power boundedness condition, the components of the surface current density parallel and orthogonal to the edge of a PEC planar scatterer behave at most like t 1 / 2 and t 1 / 2 , respectively, as the distance t from the edge approaches zero [30,31].
On the other hand, Equation (5) can be deeply simplified by taking advantage of the azimuthal symmetry of the geometry at hand. As a matter of fact, the surface current density on the disk can be conveniently expanded in a Fourier series,
J _ x , y = n = + J _ n ρ e j n ϕ ,  
where the polar coordinates ρ , ϕ such that x = ρ c ϕ = ρ cos ϕ and y = ρ s ϕ = ρ sin ϕ are introduced. Hence, by means of simple algebraic manipulations, it is possible to write
J ~ _ u , v = J ~ x u , v J ~ y u , v = 1 2 π j n = + j n c ψ j s ψ s ψ j c ψ J ¯ ~ _ n w e j n ψ ,  
where u = w cos ψ and v = w sin ψ , and the following definition of VHT of the order n (VHTn) of the n -th azimuthal harmonic of the surface current density is introduced:
J ¯ ~ _ n w = J ~ C n w j J ~ D n w = 0 + H _ _ n w ρ J ρ n ρ j J ϕ n ρ ρ d ρ = 0 + H _ _ n w ρ J ¯ _ n ρ ρ d ρ ,  
with
H _ _ n w ρ = J n w ρ n J n w ρ w ρ n J n w ρ w ρ J n w ρ ,  
where J n · and J n · denote the n-th order, first-kind Bessel function and its first derivative. By using (7), Equation (5) can be rewritten as follows:
n = + e j n ϕ 0 + H _ _ n w ρ G ~ _ _ w J ¯ ~ _ n w w d w = 1 4 π 2 c ϕ s ϕ j s ϕ j c ϕ 0 + 0 2 π G ~ _ _ 2,1 S w , ψ e j w ρ c ψ ϕ + w ρ 0 c ψ ϕ 0 + R 2 w d + R 1 w z 0 w d w d ψ J _ 0 ,  
where x 0 = ρ 0 c ϕ 0 , y 0 = ρ 0 s ϕ 0 , and
G ~ _ _ w = w 2 B 2,2 P w + B 2,2 S w e j 2 R 2 w d 0 0 A 2,2 P w + A 2,2 S w e j 2 R 2 w d .

2.2. Proposed Solution

By virtue of Helmholtz decomposition, the surface current density on the disk can be thought of as the superposition of a surface curl-free contribution and a surface divergence-free contribution [32]. As a result, the n -th azimuthal harmonic of the surface current density can be rewritten as follows:
J ¯ _ n ρ = d / d ρ n / ρ Φ C n ρ J ¯ C n ρ j n / ρ d / d ρ Φ D n ρ J ¯ D n ρ ,  
where Φ T n · with T = C , D are suitable scalar potential functions.
Despite what happens to the cylindrical components of J ¯ _ n ρ , the VHTn of J ¯ T n ρ are scalar functions [33], i.e.,
J ~ C n w 0 = 0 + H _ _ n w ρ J ¯ C n ρ ρ d ρ ,  
0 j J ~ D n w = 0 + H _ _ n w ρ J ¯ D n ρ ρ d ρ ,  
or, equivalently,
J ~ T n w = 0 + J n w ρ Φ T n ρ w ρ d ρ ,  
thus suggesting J ¯ T n ρ as new unknowns in the spatial domain in order to manage scalar unknowns in the spectral domain.
The discretization of the integral equation in (10) is provided by using the Galerkin method. For this purpose, the following complete and non-redundant expansions of the scalar unknowns in the spectral domain of the orthonormal eigenfunctions of the most singular part of the integral operator that allow the reconstruction of the physical behavior of the surface current density are considered [34]:
J ~ T n w = h = 1 + δ n , 0 + γ T , h n f ~ T n , h w ,  
where
f ~ T n , h w = 0 + J n w ρ f T n , h ρ w ρ d ρ = 2 n + 2 h + p T + 1 J n + 2 h + p T + 1 a w w p T ,  
and where f T n , h ρ denote the expansion functions of Φ T n · , p C = 3 / 2 and p D = 1 / 2 in order to reconstruct the appropriate edge behavior of the surface current density, and the condition
γ D , 1 n = α n γ C , 1 n ,  
where α n = j a   s g n n 4 n 2 1 , has to be verified to guarantee the vanishing of the reconstructed current for ρ > a .
Simple algebraic manipulations, the Galerkin method, the orthogonality property of the azimuthal harmonics, and the extraction of the most singular part of the integral operator allow the recasting of the integral equation in (10) as the following second-kind matrix equation:
x n + A n x n = c n ,  
where, by setting α ^ n = 1 1 + k 2 α n 2 , the unknown vector is
x n = 1 2 ω ε 2 x C n ω μ 2 2 x D n ,  
x C n = x C , h n h = 1 + δ n , 0 + ,  
x D n = x D , h n h = 0 + ,  
x C , 1 n = γ C , 1 n α ^ n       f o r       n 0 ,  
x T , h n = γ T , h n       f o r       h > 0 ,  
the coefficient matrix is
A n = j 2 ω ε 2 A C , C n 2 j ε 2 μ 2 A C , D n 2 j ε 2 μ 2 A D , C n j 2 ω ε 2 A D , D n ,  
A C , C n = A C , C , k , h n k , h = 1 + δ n , 0 + ,  
A C , D n = A C , D , k , h n k = 1 + δ n , 0 ,   h = 0 + ,  
A D , C n = A D , C , k , h n k = 0 ,   h = 1 + δ n , 0 + ,  
A D , D n = A D , D , k , h n k , h = 0 + ,  
A C , C , 1 , 1 n = α ^ n 2 M C , 1 , 1 n α n 2 M D , 1 , 1 n       f o r       n 0 ,  
A C , C , 1 , h n = A C , C , h , 1 n = α ^ n M C , 1 , h n       f o r       n 0 ,   h 0 ,  
A C , C , k , h n = A C , C , h , k n = M C , k , h n       f o r       k , h 0 ,  
A C , D , 1 , h n = A D , C , h , 1 n = α ^ n α n M D , 1 , h n       f o r       n 0 ,   h 0 ,  
A C , D , k , h n = A D , C , h , k n = 0       f o r       k , h 0 ,  
A D , D , k , h n = A D , D , h , k n = M D , k , h n       f o r       k , h 0 ,  
M C , k , h n = 0 + w 2 B 2,2 P w j w 2 ω ε 2 f ~ C n , k w f ~ C n ,   h w w d w + 0 + w 2 B 2,2 S w e j 2 R 2 w d f ~ C n , k w f ~ C n , 𝒽 w w d w ,  
M D , k , h n = 0 + w 2 A 2,2 P w + j ω μ 2 2 w f ~ D n , k w f ~ D n , h w w d w + 0 + w 2 A 2,2 S w e j 2 R 2 w d f ~ D n , k w f ~ D n , 𝒽 w w d w ,  
and the free term is
c n = j 1 2 ω ε 2   c C m j ω μ 2 2 c D m ,  
c C n = c C , k n k = 1 + δ n , 0 + ,  
c D n = c D , k n k = 0 + ,  
c C , 1 n = α ^ n b C , 1 n α n b D , 1 n       f o r       n 0 ,  
c T , k n = b T , k n       f o r       k > 0 ,  
b C , k n = e j n ϕ 0 2 π 0 + 1 2 w 2 B 2,1 S u , v J 0 x j J 0 y e j ϕ 0 J n 1 w ρ 0 J 0 x + j J 0 y e j ϕ 0 J n + 1 w ρ 0 j w C 2,1 S u , v J 0 z J n w ρ 0 e j R 2 w d + R 1 w z 0 f ~ C n , k w w d w ,  
b D , k n = j e j n ϕ 0 2 π 0 + 1 2 w 2 A 2,1 S u , v J 0 x j J 0 y e j ϕ 0 J n 1 w ρ 0 + J 0 x + j J 0 y e j ϕ 0 J n + 1 w ρ 0 e j R 2 w d + R 1 w z 0 f ~ D n , k w w d w ,  
for each n integer number. Following the reasoning proposed by the authors in [33] and the references therein, both the compactness of the matrix operator A n in l 2 and the l 2 -norm boundedness of the free term c n can be demonstrated. As a result, according to Fredholm theory, the existence of the solution originates from its uniqueness, and the approximate solution obtained by truncating the matrix equation in (18) converges to the exact solution as the truncation order tends to infinity.

3. Numerical Results

Of course, the coefficient matrix and the free term elements of the Fredholm second-kind equation in (18) have to be numerically evaluated. To this purpose, it is worth observing that since the integrands of the second integral in (20l) and (20m) and of the integrals in (21f) and (21g) are asymptotically exponentially decaying functions, the corresponding integrals quickly converge. On the other hand, the first integral in (20l) and the first integral in (20m), which are improper integrals of oscillating functions with an asymptotic algebraic decay, can be rewritten as combinations of fast-converging proper integrals by means of an analytical procedure in the complex plane proposed by the authors in a previous paper [26]. Henceforth, the numerical simulations were provided by means of in-house C++ software code running on a laptop equipped with an Intel Core i7-10510U 1.8 GHz processor and 16 GB of RAM, and the integrals were evaluated by means of an adaptive scheme based on a 24-point Gauss–Legendre quadrature rule.
The fast convergence of the proposed method can be quantified by introducing the following normalized truncation error:
e r r M , N = n = N + 1 N 1 x M + 1 n x M n 2 n = N + 1 N 1 x M n 2 ,    
where the symbols · and x M n denote, respectively, the Euclidean norm and the vector of the first M expansion coefficients for each unknown, whereas 2 N 1 is the number of azimuthal harmonics used. It can be shown that the convergence rate is slightly affected by the disk depth and the position and orientation of the dipole. However, as expected, it is affected by the disk radius in terms of both the number of expansion functions and azimuthal harmonics to be used in order to achieve a predetermined accuracy. Such behavior is shown in Figure 2a,b, where the normalized truncation error is plotted as a function of M for N = 7 and as a function of N for M = 8 for a zero-thickness PEC disk of radius a = 5 ,   10 ,   15 ,   20   c m buried in a lossy half-space with ε r 2 = 3.5 j 0.3 and μ r 2 = 1 at a depth d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located in the upper half-space filled by air, ε r 1 = μ r 1 = 1 , at x 0 = 2   m ,   y 0 = 2   m , z 0 = 2   m . It is worth noting that M = 8 and N = 7 achieved a normalized truncation error less than 10 3 for the worst case examined, i.e., for a = 20   c m , with a computation time of 50 s.
Validation of the developed in-house C++ software code was carried out by means of comparisons with CST-MWS. To this purpose, it must be emphasized that the commercial software was not designed to deal with scattering problems involving layered media, thus forcing an approximation of the infinite layered medium to a finite one. Unfortunately, a fairly accurate solution can be obtained only when a very large, with respect to the wavelength, layered region is involved, with an exorbitant number of mesh-cells and an endless computation time to perform the simulations. For this reason, comparisons are provided only in homogeneous medium, in particular, even when medium 2 in Figure 1 is filled with air, i.e., ε r 1 = μ r 1 = ε r 2 = μ r 2 = 1 . In Figure 3, the amplitude of the components of the surface current density on a zero-thickness PEC disk of radius a = 10   c m , located at d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located at x 0 = 2   m ,   y 0 = 2   m , z 0 = 2   m , is plotted along the x and y axes, i.e., for x / a 1 , y = 0 , and for x = 0 , y / a 1 , and compared with the ones reconstructed by means of CST-MWS. The agreement is excellent. However, the proposed method achieved a normalized truncation error less than 10 3 for M = 4 and N = 5 with a computation time of 15 s, whereas the time-domain solver of CST-MWS required about 30 million mesh-cells with a computation time of about 3 h to accurately reconstruct the solution.
As an example of the application of the proposed method, a configuration was considered that resembled a measurement system of the field scattered from a landmine based on a two drones, one equipped with a transmitting antenna and the other equipped with a receiving antenna. As with the cases examined before, a zero-thickness PEC disk, of radius a = 10   c m buried in a lossy half-space with ε r 2 = 3.5 j 0.3 and μ r 2 = 1 at the depth d = 3   c m , was considered. Now, the transmitting Hertzian dipole, directed along the y axis with J _ 0 = 10 3 A m and operating at 800 MHz, was located in the upper half-space filled by air, ε r 1 = μ r 1 = 1 , at z 0 = 2   m and four different positions in the xy plane: (1) x 0 = 2   m ,   y 0 = 2   m , (2) x 0 = 0   m ,   y 0 = 2   m , (3) x 0 = 2   m ,   y 0 = 0   m , and (4) x 0 = 0   m ,   y 0 = 0   m . The scattered field was evaluated in the horizontal area x , y 2   m , 2   m × 2   m , 2   m for z = 1   m . In Figure 4, two separate plots for the amplitudes of the incidence and the scattered electric fields in the area of interest are provided for the four cases examined. In the figures, the disk circumference is denoted by the black dotted circular line. Clearly, the reflection due to the discontinuity interface between the two involved media perturbs the classical spherical wave behavior of the field radiated by the dipole in a homogeneous medium. On the other hand, the disk “footprint” is clearly observable in the scattered field behavior. As a matter of fact, the bright spot (in yellow) is elliptical with a maximum shifted in the opposite direction with respect to the dipole position, except, of course, for the last case examined, in which the dipole is along the disk axis.

4. Conclusions

In this paper, an analysis of the dipole field scattering from a zero-thickness PEC disk buried in a lossy half-space has been carried out by generalizing the guaranteed-convergence Helmholtz–Galerkin technique proposed by the authors in previous papers for dealing with plane-wave scattering from a disk and a holed plane in a homogeneous medium. The advantage of using this method is twofold. First, convergence of the method to the exact solution of the problem is guaranteed due to the Fredholm second-kind nature of the obtained matrix equation. Second, the convergence is very fast because the unknowns are expanded by means of functions reconstructing the physical behavior of the fields. As a matter of fact, the proposed numerical results demonstrate that the numerical resources, in terms of both computation time and memory occupation, to achieve highly accurate results are paltry and drastically less than those required by CST-MWS. Therefore, the proposed method can be particularly suitable for the implementation of effective inverse scattering iterative algorithms. Moreover, it can serve as a reference for the validation of general-purpose commercial software without a mathematically guaranteed convergence.

Author Contributions

Conceptualization, M.L. and M.D.M.; methodology, M.L.; software, M.L.; validation, G.A.C.; formal analysis, M.L.; data curation, G.C. and D.P.; writing—original draft preparation, M.L., M.D.M. and F.S.; writing—review and editing, M.L., G.A.C., M.D.M. and F.S.; supervision, M.L.; funding acquisition, M.D.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Italian Ministry of University and Research (MUR) under the Program for the Development of Research of National Interest (PRIN) 2022 “TERRAIN”, CUP H53D23007300001, funded by the EU in the NextGenerationEU plan.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study were obtained by means of in-house software code implementing the proposed method and by using the commercial software CST Microwave Studio 2023.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Robledo, L.; Carrasco, M.; Mery, D. A survey of land mine detection technology. Int. J. Remote Sens. 2009, 30, 2399–2410. [Google Scholar] [CrossRef]
  2. Amiri, A.; Tong, K.; Chetty, K. Feasibility study of multi-frequency ground penetrating radar for rotary UAV platforms. In Proceedings of the IET International Conference on Radar Systems, Glasgow, UK, 22–25 October 2012. [Google Scholar]
  3. Colorado, J.; Perez, M.; Mondragon, I.; Mendez, D.; Parra, C.; Devia, C.; Martinez-Moritz, J.; Neira, L. An integrated aerial system for landmine detection: SDR-based ground penetrating radar onboard an autonomous drone. Adv. Robot. 2017, 31, 791–808. [Google Scholar] [CrossRef]
  4. Sipos, D.; Planinsic, P.; Gleich, D. On drone ground penetrating radar for landmine detection. In Proceedings of the IEEE 1st International Conference on Landmine, Detection, Clearance Legislations, Beirut, Lebanon, 26–28 April 2017. [Google Scholar]
  5. Ferrara, V.; Pietrelli, A.; Chicarella, S.; Pajewski, L. GPR/GPS/IMU system as buried objects locator. Measurements 2018, 114, 534–541. [Google Scholar] [CrossRef]
  6. García Fernández, M.; Álvarez López, Y.; Arboleya Arboleya, A.; González Valdés, B.; Rodríguez Vaqueiro, Y.; Las-Heras Andrés, F.; Pino García, A. Synthetic aperture radar imaging system for landmine detection using a Ground Penetrating Radar on board a Unmanned Aerial Vehicle. IEEE Access 2018, 6, 45100–45112. [Google Scholar] [CrossRef]
  7. Pinchera, D.; Migliore, M.D. Efficient synthesis of concentric-rings plane wave generators. IEEE Trans. Antennas Propag. 2023, 71, 4967–4975. [Google Scholar] [CrossRef]
  8. O’Neill, K.; Haider, S.A.; Geimer, S.D.; Paulsen, K.D. Effects of the ground surface on polarimetric features of broadband radar scattering from subsurface metallic objects. IEEE Trans. Geosci. Remote Sens. 2001, 39, 1556–1565. [Google Scholar] [CrossRef]
  9. Demarest, K.; Plumb, R.; Huang, Z. FDTD modeling of scatterers in stratified media. IEEE Trans. Antennas Propag. 1995, 43, 1164–1168. [Google Scholar] [CrossRef]
  10. Huang, Z.; Demarest, K.R.; Plumb, R.G. An FDTD/MoM hydrid technique for modeling complex antennas in the presence of heterogeneous grounds. IEEE Trans. Geosci. Remote Sens. 1999, 37, 2692–2698. [Google Scholar] [CrossRef]
  11. Cui, T.J.; Chew, W.C. Fast algorithm for electromagnetic scattering by buried conducting plates of large size. IEEE Trans. Antennas Propag. 1999, 47, 1116–1118. [Google Scholar] [CrossRef]
  12. Zhuang, L.; He, S.Y.; Ye, X.B.; Hu, W.D.; Yu, W.X.; Zhu, G.Q. The BCGS-FFT method combined with an improved discrete complex image method for EM scattering from electrically large objects in multilayered media. IEEE Trans. Geosci. Remote Sens. 2010, 48, 1180–1185. [Google Scholar] [CrossRef]
  13. Losada, V.; Boix, R.R.; Medina, F. Fast and accurate algorithm for the short-pulse electromagnetic scattering from conducting circular plates buried inside a lossy dispersive half-space. IEEE Trans. Geosci. Remote Sens. 2003, 41, 988–997. [Google Scholar] [CrossRef]
  14. Alcuntu, Y.; Yapar, A.; Akduman, I. On the scattering of electromagnetic waves by bodies buried in a half-space with locally rough interface. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1435–1443. [Google Scholar]
  15. Hongo, K.; Jafri, A.D.U.; Naqvi, Q.A. Scattering of electromagnetic spherical wave by a perfectly conducting disk. Prog. Electromagn. Res. 2012, 129, 315–343. [Google Scholar] [CrossRef]
  16. Balaban, M.V.; Sauleau, R.; Garcia-Vigueras, M.; Nosich, A.I. Terahertz range elementary dipole excitation of a thin dielectric disk sandwiched between two graphene covers: Integral equation analysis. In Proceedings of the IEEE International Conference on Microwaves, Antennas, Communications and Electronic Systems (COMCAS), Tel-Aviv, Israel, 4–6 November 2019. [Google Scholar]
  17. Herasymova, D.O.; Dukhopelnykov, S.V.; Lucido, M.; Nosich, A.I. Optical sensing of electron-beam position with twin silver nanotube antenna tuned to hybrid surface plasmon resonance. IEEE J. Sel. Top. Quantum Electron. 2021, 27, 4601008. [Google Scholar] [CrossRef]
  18. Kantorovich, L.V.; Akilov, G.P. Functional Analysis, 2nd ed.; Pergamon/Oxford-Elmsford: New York, NY, USA, 1982. [Google Scholar]
  19. Tsalamengas, J.L. Exponentially converging Nystrom’s methods for systems of singular integral equations with applications to open/closed strip- or slot-loaded 2-D structures. IEEE Antennas Propag. Mag. 2006, 54, 1549–1558. [Google Scholar] [CrossRef]
  20. Bulygin, V.S.; Nosich, A.I.; Gandel, Y.V. Nystrom-type method in three-dimensional electromagnetic diffraction by a finite PEC rotationally symmetric surface. IEEE Trans. Antennas Propag. 2012, 60, 4710–4718. [Google Scholar] [CrossRef]
  21. Bulygin, V.S.; Benson, T.M.; Gandel, Y.V.; Nosich, A.I. Full-wave analysis and optimization of a TARA-like shield-assisted paraboloidal reflector antenna using Nystrom-type method. IEEE Trans. Antennas Propag. 2013, 61, 4981–4989. [Google Scholar] [CrossRef]
  22. Lucido, M.; Schettino, F.; Migliore, M.D.; Pinchera, D.; Di Murro, F.; Panariello, G. Electromagnetic scattering by a zero-thickness PEC annular ring: A new highly efficient MoM solution. J. Electromagn. Waves Appl. 2017, 31, 405–416. [Google Scholar] [CrossRef]
  23. Lucido, M.; Balaban, M.V.; Nosich, A.I. Terahertz-range plasmon and whispering gallery mode resonances in the plane wave scattering from thin microsize dielectric disk with graphene covers. Proc. Math. Phys. Eng. Sci. 2022, 478, 2262. [Google Scholar] [CrossRef]
  24. Lucido, M.; Nosich, A.I. Analytical regularization approach to plane wave diffraction from circular hole in infinite resistive plane. IEEE Trans. Antennas Propag. 2023, 71, 6878–6892. [Google Scholar] [CrossRef]
  25. Chew, W.C.; Kong, J.A. Resonance of nonaxial symmetric modes in circular microstrip disk antenna. J. Math. Phys. 1980, 21, 2590–2598. [Google Scholar] [CrossRef]
  26. Lucido, M.; Di Murro, F.; Panariello, G. Electromagnetic scattering from a zero-thickness PEC disk: A note on the Helmholtz-Galerkin analytically regularizing procedure. Progr. Electromagn. Res. Lett. 2017, 71, 7–13. [Google Scholar] [CrossRef]
  27. Chew, W.C. Waves and Fields in Inhomogeneous Media; IEEE Press: Piscataway, NJ, USA, 1995. [Google Scholar]
  28. Chew, W.C.; Chen, S.Y. Response of a point source embedded in a layered medium. IEEE Antennas Wirel. Propag. Lett. 2003, 2, 254–258. [Google Scholar] [CrossRef]
  29. Michalski, K.A.; Zheng, D. Rigorous analysis of open microstrip lines of arbitrary cross section in bound and leaky regimes. IEEE Trans. Microw. Theory Tech. 1989, 37, 2005–2010. [Google Scholar] [CrossRef]
  30. Meixner, J. The behavior of electromagnetic fields at edges. IEEE Trans. Antennas Propag. 1972, 20, 442–446. [Google Scholar] [CrossRef]
  31. Makarov, G.I.; Osipov, A.V. Structure of Meixner’s series. Radiophys. Quantum Electron. 1986, 29, 544–549. [Google Scholar] [CrossRef]
  32. Van Bladel, J. A discussion of Helmholtz’ theorem on a surface. Arch. Elektron. Ubertrag. 1993, 47, 131–136. [Google Scholar]
  33. Lucido, M.; Balaban, M.V.; Dukhopelnykov, S.V.; Nosich, A.I. A fast-converging scheme for the electromagnetic scattering from a thin dielectric disk. Electronics 2020, 9, 1451. [Google Scholar] [CrossRef]
  34. Wilkins, J.E. Neumann series of Bessel functions. Trans. Am. Math. Soc. 1948, 64, 359–385. [Google Scholar] [CrossRef]
Figure 1. Geometry of the problem.
Figure 1. Geometry of the problem.
Applsci 14 05544 g001
Figure 2. Normalized truncation error for a zero-thickness PEC disk of radius a = 5 ,   10 ,   15 ,   20   c m buried in a lossy half-space with ε r 2 = 3.5 j 0.3 and μ r 2 = 1 at a depth d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located in air, ε r 1 = μ r 1 = 1 , at x 0 = 2   m ,   y 0 = 2   m , z 0 = 2   m , as a function of (a) M for N = 7 and as a function of (b) N for M = 8 .
Figure 2. Normalized truncation error for a zero-thickness PEC disk of radius a = 5 ,   10 ,   15 ,   20   c m buried in a lossy half-space with ε r 2 = 3.5 j 0.3 and μ r 2 = 1 at a depth d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located in air, ε r 1 = μ r 1 = 1 , at x 0 = 2   m ,   y 0 = 2   m , z 0 = 2   m , as a function of (a) M for N = 7 and as a function of (b) N for M = 8 .
Applsci 14 05544 g002
Figure 3. Amplitude of the components of the surface current density along the x and y axes on a zero-thickness PEC disk of radius a = 10   c m , located at d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located at x 0 = 2   m ,   y 0 = 2   m , z 0 = 2   m , for ε r 1 = μ r 1 = ε r 2 = μ r 2 = 1 . (a) J ρ along the x axis, (b) J ρ along the y axis, (c) J ϕ along the x axis, and (d) J ϕ along the y axis. Blue solid line: this method; red dashed line: CST-MWS.
Figure 3. Amplitude of the components of the surface current density along the x and y axes on a zero-thickness PEC disk of radius a = 10   c m , located at d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located at x 0 = 2   m ,   y 0 = 2   m , z 0 = 2   m , for ε r 1 = μ r 1 = ε r 2 = μ r 2 = 1 . (a) J ρ along the x axis, (b) J ρ along the y axis, (c) J ϕ along the x axis, and (d) J ϕ along the y axis. Blue solid line: this method; red dashed line: CST-MWS.
Applsci 14 05544 g003
Figure 4. Amplitudes of the incident and scattered electric fields in the horizontal area x , y 2   m , 2   m × 2   m , 2   m for z = 1   m in the case of a zero-thickness PEC disk of radius a = 10   c m buried in a lossy half-space with ε r 2 = 3.5 j 0.3 and μ r 2 = 1 at a depth d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located in the upper half-space filled by air, ε r 1 = μ r 1 = 1 , at z 0 = 2   m , and 4 different positions in the xy plane: (a,b) x 0 = 2   m ,   y 0 = 2   m ; (c,d) x 0 = 0   m ,   y 0 = 2   m ; (e,f) x 0 = 2   m ,   y 0 = 0   m ; and (g,h) x 0 = 0   m ,   y 0 = 0   m . The black dotted circular line denotes the disk circumference.
Figure 4. Amplitudes of the incident and scattered electric fields in the horizontal area x , y 2   m , 2   m × 2   m , 2   m for z = 1   m in the case of a zero-thickness PEC disk of radius a = 10   c m buried in a lossy half-space with ε r 2 = 3.5 j 0.3 and μ r 2 = 1 at a depth d = 3   c m , and illuminated by a Hertzian dipole directed along the y axis with J _ 0 = 10 3 A m , operating at 800 MHz and located in the upper half-space filled by air, ε r 1 = μ r 1 = 1 , at z 0 = 2   m , and 4 different positions in the xy plane: (a,b) x 0 = 2   m ,   y 0 = 2   m ; (c,d) x 0 = 0   m ,   y 0 = 2   m ; (e,f) x 0 = 2   m ,   y 0 = 0   m ; and (g,h) x 0 = 0   m ,   y 0 = 0   m . The black dotted circular line denotes the disk circumference.
Applsci 14 05544 g004aApplsci 14 05544 g004b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lucido, M.; Casula, G.A.; Chirico, G.; Migliore, M.D.; Pinchera, D.; Schettino, F. Helmholtz–Galerkin Technique in Dipole Field Scattering from Buried Zero-Thickness Perfectly Electrically Conducting Disk. Appl. Sci. 2024, 14, 5544. https://doi.org/10.3390/app14135544

AMA Style

Lucido M, Casula GA, Chirico G, Migliore MD, Pinchera D, Schettino F. Helmholtz–Galerkin Technique in Dipole Field Scattering from Buried Zero-Thickness Perfectly Electrically Conducting Disk. Applied Sciences. 2024; 14(13):5544. https://doi.org/10.3390/app14135544

Chicago/Turabian Style

Lucido, Mario, Giovanni Andrea Casula, Gaetano Chirico, Marco Donald Migliore, Daniele Pinchera, and Fulvio Schettino. 2024. "Helmholtz–Galerkin Technique in Dipole Field Scattering from Buried Zero-Thickness Perfectly Electrically Conducting Disk" Applied Sciences 14, no. 13: 5544. https://doi.org/10.3390/app14135544

APA Style

Lucido, M., Casula, G. A., Chirico, G., Migliore, M. D., Pinchera, D., & Schettino, F. (2024). Helmholtz–Galerkin Technique in Dipole Field Scattering from Buried Zero-Thickness Perfectly Electrically Conducting Disk. Applied Sciences, 14(13), 5544. https://doi.org/10.3390/app14135544

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop