 Open Access
 Total Downloads : 4
 Authors : Praveen Sahwal
 Paper ID : IJERTCONV1IS02044
 Volume & Issue : NCEAM – 2013 (Volume 1 – Issue 02)
 Published (First Online): 30072018
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
MICROWAVE DETECTION AND IDENTIFICATION OF BURIED METALLIC OBJECTS
MICROWAVE DETECTION AND
ISBN: 9789383758098
IDENTIFICATION OF BURIED METALLIC OBJECTS
1Praveen Sahwal, Assistant Professor(ECE), Geeta Engineering College, Naultha, (Panipat), Haryana, INDIA.
Abstract–In this paper the idea of using acoustically induced Doppler spectra as a means for metallic object detection and identification is introduced. AcoustoEM wave interaction occurs when an electromagnetic wave scatters from an object under acoustic illumination. The incident acoustic wave causes a boundary deformation within the metallic cylinder. The acoustically vibrating cylinder gives rise to a frequency modulated scattered electromagnetic field which is a function of the cylinders natural resonance frequency and both the electromagnetic and acoustic source parameters. Results indicate that the scattered Doppler frequencies correspond to the mechanical vibration frequencies of the cylinder, and the sidelobe Doppler spectrum level is, to the first order, linearly proportional to the degree of deformation. The prior knowledge of doppler responses and signature of each metal help us to detect and identify the buried metallic cylinder.
Keywords: – Electromagnetic Scattering, Acoustics, Metallic cylinders, Doppler spectrum, MATLAB.
I. INTRODUCTION
analytical approach for scattering solutions and was applied by Rayleigh [5] and Maxwell [6] to certain scalar field problems. These methods have found application in electromagnetic scattering from rough surfaces [7, 8], cylinders [9, 10 and 11] and spheres [12, 13].
Fig. 1: Setup for object detection

ANALYTICAL DEVELOPMENT
Here, the bistatic scattered field from a perfectly conducting circular cylinder with arbitrary deformation and illuminated by an electromagnetic plane wave is formulated. Anejwet time dependence for electromagnetic waves and an ejwat time dependence for acoustic waves is understood in all expressions representing waves and is suppressed.
TM Mode
Consider an incident TM plane wave impinging upon a perfectly conducting, slightly deformed circular cylinder as shown in Fig. 2. For a vibrating cylinder whose cross section shape is varying with time, the scattered field will be timevarying which gives rise to the scattered Doppler spectrum. The perimeter of the cylinders surface can be expressed in polar coordinates as
Fig. 2
p, = a + bf(,, t) (1)
Where a is the unperturbed cylinder radius, f(, t) is a
periodic and smooth function of , and b is a perturbation constant assumed to be much smaller than the wavelength.
n n
I (j)nC J (ka)ejn
n=
2jf(, t)
ISBN: 9789383758098
e
jn
n
n
The standard method for the computation of the scattered field from a circular cylinder is the eigenfunction
= nka I (j)
n=
(7)
H(2)(ka)
expansion of the total field. Following the development in
[13], the incident plane wave propagating along the x direction with a polarization parallel to the cylinders axisNow if f(, t) has a Fourier series expansion with respect to
can be expressed as
Ei = ejkpcos = I (j)nJ (kp)ejn
f(, t) = I A (t)ejm
m
m=
(8)
z n
n=
(2)
and is inserted into (7), the resulting solution for the coefficients, Cn, is
Where k is the wavenumber. When there is no perturbation on the surface (i.e.,b = 0 ), the scattered field is given by
Cn(t) =
2(j)n+1
I (j)p
p
Anp
(t)
(9)
Es = I (j)n
Jn(ka)
H(2)(kp)ejn
nkaJn(ka)
p=
H(2)(ka)
z
n=
H(2)(ka) n
n
(3)
Hence for a given displacement function, f(, t), the Cn s calculated from (9) can be used in (4) to yield the bistatic
If a small perturbation is introduced on the boundary, the scattered field may be expanded in terms of a perturbation series in kb. To the first order in, we may write
scattered field at any angle . In this derivation it is
implicitly assumed that the time rate of variations of f(, t) is much slower than the electromagnetic frequency so that relativistic effects can be neglected. Under this
Es = I (j)n
Jn(ka)
(1 + C kb)H(2)(kp)ejn
assumption, a Fourier transform of the timevarying scattered field will provide the scattered Doppler spectrum.
z
n=
H(2)(ka) n n
n
(4)

ACOUSTIC VIBRATION OF A SOLID
CYLINDER
Where the Cn s are unknown coefficients to be determined
using the boundary condition Ei + ES = 0 on the surface
Mathematical Formulation for Displacement
z z
of the perturbed cylinder. This expression is, of course, a
permissible solution of Maxwells equations since it satisfies both the wave equation and the radiation condition. Using Taylor series expansion, the eigenfunctions in (2) and (4) can also be approximated to the first order in kb, i.e.,
Zn(kp) = Zn(k(a + bf(, t)))
Zn(ka) + kbf(, t)Z (ka)
In order to calculate the deformation of a circular cylinder due to a timeharmonic incident acoustic wave, the acoustic scattering from the cylinder must be evaluated. An interesting aspect of acoustic scattering from solid objects is that mechanical waves within the solid are excited which consist of both transverse shear waves and longitudinal compressional waves. The equation of motion for a solid elastic medium can be written [17]
The wave inside the cylinder will be represented by suitable
(5)
n solution of the equation of motion of a solid elastic
Where Zn represents the Bessel or Hankel function of nth order. Applying the boundary condition at the surface of the cylinder
CX)
medium, which may be written as
( ) ( )
a2u
,1 + 2Âµ Âµ Ã— 2w = p1 at2
(10)
I (j)n[J (ka) + kbf(,, t)J, (ka)]ejn1
E(1 a) E
n
n=CX)
n (. u)
(1 + a)(1 2a)
2(1 + a)
Ã— ( Ã— u)
CX)
J (ka)
a2u
= I (j)n n (1
H(2)(ka)
= p1 at2
(11)
n=CX) n
n
+ Cnkb) [H(2)(ka)
n
+ kbf(,, t)H(2)1 (ka)]ejn1
(6)
Neglecting the (kb)2 term and using the Wronskian relationship for Bessel functions results in the following expression:
Where p1 is the density of the scatterer, u is the displacement, and E and a are Youngs modulus and Poissons ratio, respectively. Solutions to this equation for a solid cylinder of infinite length are of the form [16]
Where
= u (12)
And
2w = Ã— u
From Eq. (14) can be derived the equations,
2 p1
a2
= ( ,1 + 2Âµ) at2 (13)
And
By the conditions of symmetry, [rz] = 0
ISBN: 9789383758098
everywhere.Upon
2(2w)
p1
= Âµ
a2(2w)
at2 (14)
substitution from Eqs. (18), (19), (20), (22), (23), and (21), the boundary condition Eqs. (24), (25), and (26), for the nth
Which define the wave velocities
1 1
mode solving these equations simultaneously for cn. The result is
(,1 + 2Âµ) 2 E(1 a) 2
c1 = p1
And
= p1(1 + a)(1 2a)
(15)
cn = P0En(j)n+1sin1n exp(j1n) ,
(27)
p>2
2
Âµ 1 E 1
Mapping to Radial Displacement
p
c2 =
1
=
2p1(1 + a)
(16)
The displacement of the cylinder is now represented by both a radial and angular displacement in (32) and (33).
Solution of Eq. (10) can be found by assuming that
displacement can be derived from a scalar and a vector potential:
u = + Ã— A.
(17)
The displacements thus can be thought of as the sum of two displacements, one associated with compressional waves and other with shear waves. If we assume that potentials satisfy the equations,
2 1 a2
Recall that the perimeter of the cylinder is expressed from
(1) as p, = a + bf(,, t) with the radial displacement given by bf(,, t). The radial and angular displacement can be mapped into only radial displacement by setting
bf(,, t) = ur(, 0, t) (28)
Where 0 is the angular shift corresponding to the angular component of displacement, u . Frombasic geometry (see
Fig. 4) and an application of Taylor series about
= c2 at2
0
u (, , t)
And
1
(18)
1 a2A
0 =
0
a
u(,, t)
u, (,, t)
0
c
2A = 2
2
at2
(19)
Solving for 0
a a
(29)
the displacement derived from it shall be symmetrical about
=0. Now, by Eq. (17) and (11),
u(,, t)
0 = a + u, (,, t)
nb d
u (,, t)
u = I n J (k r) a J (k r) cosn0
r r n 2 n dr n 1 a
n=0
(20)
(30)
u = I nan J (k r) b d J (k r) sinn0
(21)
0
n=0
r n 1
n dr n 2
The factors cn are the unknown coefficients which must be evaluated.
The following boundary conditions are applied at the surface of the cylinder:
pi + ps = [rr] at r
= a, (24)
ui, r + us, r = ur at r
= a, (25)
And
[r0] = [rz] = 0 at r = a. (26)Fig. 4: Mapping Radial Displacement
Since u, (,, t) a. And substituting (30) into (28) gives
, , u(,, t)
bf( , t) = ur a , t .
(31)
Rewriting ur and u with time dependence explicitly shown
CX)
ur(,, t) = I Ur,n cos(wat + 0p,n) cos (n,)
n=0
CX)
(32)
Fig. 3
u(,, t) = I U,n cos(wat + 0,n) sin (n,)
n=0
(33)
Where Ur,nand U,nare the magnitudes of the coefficients
The doppler
ISBN: 9789383758098
response of the scattered field is plotted
for the nth mode given by the expressions in brackets in (20) and (21), respectively, and 0p,nand 0,nare the respective phases of these coefficients. Substituting (32)and (33) into (31).
When only one mode has significant displacement, however, for the nth mode simplifies to
bf = Ur,nCOS(wat + 0p,n) COS n ,
U,nCOS(wat + 0,n) Sin(n,)
a
(34)
this expression should be expanded in a Fourier series about
. The coefficients, Am(t) , of the expansion can then be used in the TM or TE electromagnetic solution to obtain the timevarying scattered field. Taking b = Ur,n, we can expand f(,, t)in a Fourier series about , .
against acoustic frequency. This is highly sensitive to resonances in the cylinder. Note that the n=2 mode is resonant at the lowest frequency, followed by the n=1 andn=3 modes. Hence, by acoustically vibrating the cylinder over a wide range of frequencies, the Doppler response will have significant Doppler components only at the mechanical resonant frequencies of the cylinder.
Fig. 7 show the scattered doppler spectrum for Steel cylinder of radius 5cm for mode n=1 and 2. The young modulus is 200 Ã— 109 Pascal, possoinratio is 0.28, density is 7.7 g/cm3.
Doppler Response For Steel Cylinder For n=1 Mode
60
Magnitude of First Doppler Harmonic(dB)
70
80
90
100
110

SIMULATIONS
Doppler Response
Fig. 6 shows the calculated Doppler spectrum of the backscattered field with TM incidence for n=2 mode. The
120
130
140
150
1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
x1=Propagation constant * Radius
spectrum is similar to a frequency modulated signal with a low index of modulation. The spectrum has frequency components at harmonics of the incident acoustic frequency of which the first is by far the largest. Other acoustic modes produce a similar Doppler spectrum where the frequency components are at harmonics of the vibration frequency and the first component is still the most significant. When the frequency of the incident acoustic wave is not close to a resonant frequency of the cylinder, it behaves as a rigid cylinder. Consequently, there is negligible displacement at the surface of the cylinder, and the magnitude of the Doppler spectrum is also negligible. However, when the frequency of the incident acoustic wave is close to a resonance of the object, the mode corresponding to this resonance is excited and a measurable Doppler spectrum results.
Excitation at Resonance
Here, the magnitude of the first harmonic of the Doppler spectrum in the backscatter field is calculated as a function of the incident acoustic frequency. This illustrate the effect of resonance on the scattered Doppler spectrum. Fig. 7 show measured and computed first harmonic of the scattered Doppler spectrum for the cylinders of various material in air.
100
Magnitude of Doppler Response(db)
50
0
50
100
150
2.5 2 1.5 1 0.5 0 0.5 1 1.5 2 2.5
frequency / f s
Fig. 6: Fourier transform of backscattered Doppler spectrum
Fig. 7 (a): Doppler Response For Steel Cylinder For n=1 Mode
Doppler Response For Steel Cylinder For n=2 Mode
X: 1.298
Y: 49.81
40
Magnitude of First Doppler Harmonic(dB)
50
60
70
80
90
100
110
120
130
1 1.2 1.4 1.6 1.8 2
x1=Propagation constant * Radius
Fig. 7 (b): Doppler Response For Steel Cylinder For n=2 Mode Here ,the profile of three cylinders of Aluminium, Steel and Brass is presented. If we have some knowledge of their pattren of doppler response and magnitude in advance, then we can easily detect the buried cylinders.
Three different mode of excitation of cylinder of different materials is given in the Table 1. In each mode cylinder has different resonant frequency and magnitude of doppler response. Resonant frequency, magnitude of doppler response and difference between 1st and 2nd harmonic of doppler response are the factors that differentiate the cylinder of different material. All these properties of different material are shown in Table 1.
This undesirable result is offset by realizing that in a practical scenario the buried objects of interest are likely to be metallic or plastic shells which will have significantly larger displacements (and lower resonant frequencies) than the virtually rigid, solid objects considered here.
ISBN: 9789383758098
Doppler Response for Aluminium Cylinder 

M od e No . 
X1=Prop agation Constant *Radius 
Reso nant Freq uenc y (KHz ) 
Magnit ude of Dopple r Respon se(dB) 
Diffe rence Betw een 1st and 2nd Harm onic 
n= 1 
1.43 
2.82 
58 
5dB 
n= 2 
1.174 
2.32 
71 
0.5d B 
n= 3 
1.83 
3.61 
36 
3dB 
Table 1: Doppler Respose for Aluminium Cylinder
IV. CONLUSION
In this work, the analytical Doppler spectrum of an acoustically vibrated circular cylinder is examined and which is strongly dependent upon the cylinders mechanical resonances. An analytical solution for the bistatic scattering from a deformed cylinder is derived using a perturbation method. The Doppler spectrum consist of harmonics of the incident acoustic frequency where the first harmonic is the most significant. Also, the first harmonic only becomes measurable when the cylinder is excited near a mechanical resonance. These results indicatethat acoustically vibrating an object at its resonant frequencies and measuring the Doppler scattered response monostatically or bistatically could provide an effective method for detecting and identifying buried objects of different metals. Afer analysing the signature of different metals, it is concluded that the signature of metal buried under the ground have very sharp response, than as it happens in case of steel that can be property of an alloy. From the signature reflected from the different metals buried under the ground, it obsereved that the profile of that metal including amplitude doppler response for resonant frequency and its harmonic are found different for different metal where it is concluded using the profile of the reflected signature, we can identify the metal of the buried objects.
REFERENCES
[1]. Daniel E. Lawrence and Kamal Sarabandi, Acoustic and Electromagnetic Wave Interaction: Analytical Formulation for AcoustoElectromagnetic Scattering Behavior of a Dielectric Cylinder, IEEE Trans.Antennas Propagation, VOL. 49, NO. 10, OCTOBER 2001
[2]. C. Stewart, Summary of mine detection research,U.S. Army engineering res. and devel. labs, Corps. ofEng., Belvoir, VA, Tech. Rep.1636TR, vol. I, pp. 172179, May 1960.
[3]. G. S. Smith, Summary report: Workshop on new directions for electromagnetic detection of nonmetallic mines, Report for U.S. Army BRDEC and ARO, June 1992. [4]. W. R. Scott, Jr. and J. S. Martin, An experimental model of a acoustoelectromagnetic sensor for detecting land mines, in IEEE APS Int. Symp. Dig., Atlanta, GA, July 1998, pp. 978981. [5]. , Experimental investigation of the acousto electromagnetic sensor for locating land mines, Proc. of SPIEInt. Society Optical Engineering, vol. 3710, no. I, pp. 204214, 1999. [6]. Daniel E. Lawrence and Kamal Sarabandi, Electromagnetic Scattering from a Dielectric Cylinder Buried Beneath a Slightly Rough Surface, IEEE Trans. Antennas Propagat. [7]. S. O. Ogunade, Electromagnetic response of an embedded cylinder for line current excitation, Geophys., vol. 46, no. 1, pp. 4552, Jan. 1981. [8]. R. F. Harrington, TimeHarmonic Electromagnetic Fields. New York: McGrawHill, 1961, pp. 232235. [9]. J. J. Faran, Sound scattering by solid cylinders and spheres, J. Acous. Soc. Amer., vol. 23, no. 4, pp. 405 418, July 1951.