# Radar Cross-Section Formulation of a Shell-Shaped Projectile Using Modified PO Analysis

Authors: Mohammad Asif Zaman; Md. Abdul Matin
Journal: Modelling and Simulation in Engineering
 (2012)

---

## Abstract

A physical optics based method is presented for calculation of monostatic Radar Cross-Section (RCS) of a shell-shaped projectile. The projectile is modeled using differential geometry. The paper presents a detailed analysis procedure for RCS formulation using physical optics (PO) method. The shortcomings of the PO method in predicting accurate surface current density near the shadow boundaries are highlighted. A Fourier transform-based filtering method is proposed to remove the discontinuities in the approximated surface current density. The modified current density is used to formulate the scattered field and RCS. Numerical results are presented comparing the proposed method with conventional PO method. The results are also compared with published results of similar objects and found to be in good agreement.


---

## Body

## 1. Introduction

Prediction and measurement of Radar Cross-Section (RCS) have been a significant area of research for scientists and engineers for many years. The widespread uses of radar technology since the Second World War have demanded accurate and efficient prediction of fields scattered by radar targets. The knowledge of echo characteristics of radar targets is of great importance for the design of high-performance radars, as well as low visibility stealth targets [1, 2].In radar technology, an antenna radiates electromagnetic (EM) energy. When an object is illuminated by the radar EM field, it reflects back some EM energy, which is received by an antenna. In monostatic radar system, the transmitting and reception of the EM energy are done by the same antenna or by multiple antennas located very close to each other [1, 3]. In bistatic radar system, separate antennas are used for transmitting and receiving, and the antennas are usually far away from each other [1]. The radar field reflective nature of an object is specified in terms of RCS. The RCS is the area that a target would have to occupy to produce the amount of reflected power that is detected back at the radar [2, 4]. The RCS of an object depends on the viewing angles, the size, geometry, and composition of the object, frequency and polarization of the radar signal, and so forth [3, 4]. Stealth technology concentrates on reducing the RCS of airplanes and missiles to make them invisible to radars. Conversely, radar engineers are developing more sensitive radars that can detect low RCS targets. In both cases, accurate numerical simulation methods are essential for design purposes.Numerical simulation of RCS of an object requires calculation of the scattered field from the object for a given incident field. Several numerical methods exist for scattered EM field calculations. Pure numerical methods such as Method of Moments (MoM), Finite-Difference Time-Domain (FDTD) method, Fast Multipole Method (FMM), and Transmission-Line Matrix (TLM) have been successfully used in predicting RCS of radar targets [4, 5].  Conformal FDTD-based methods have also appeared in literature [6]. Recently, some variants of MoM have been developed for monostatic RCS formulation [7, 8] and other scattering related problems [9]. These methods do not depend on the geometry and can be used for any arbitrary shaped objects. However, these methods are computationally demanding, and the high simulation time for electrically large objects is not always acceptable for design and optimization problems. High-frequency asymptotic methods, such as Geometrical Optics (GO), Physical Optics (PO), Uniform Geometrical Theory of Diffraction (UTD), and Physical Theory of Diffraction (PTD), have also been used for RCS formulation [4]. These methods are based on local interaction of EM fields. Therefore, they are computationally less demanding than the pure numerical methods, and they require far less simulation time [10]. However, these methods are geometry dependent. For complex shaped objects, the scattered field formulation can be tedious, especially when using UTD method [11]. GO and PO methods do not suffer to the same extent as UTD for complex geometry cases. The GO method is the fastest among the high frequency techniques, but it is relatively less accurate [11]. The PO method gives much more accurate results compared to GO. It is a well-accepted method for formulating scattered field from electrically large objects [12, 13].Because of its relatively high accuracy, the PO method has been widely used for RCS formulation [14, 15]. The PO method is also computationally more efficient than pure numeric methods, making it relatively faster which makes it an essential tool for aerospace designers [16]. Several modified versions of PO have been developed to further increase the speed and accuracy of the PO method [17, 18]. In this paper, a new modified PO method is used to formulate the monostatic RCS of a shell-shaped projectile.This paper presents a detailed description and procedure of RCS calculation method of an object using PO method and modified PO method. The procedure includes the geometrical modeling of the shell-shaped radar object, approximation of the induced surface current on the object, and formulation of the PO radiation integral in parametric space. Although the paper concentrates on RCS formulation of a specific shell-shaped radar target, the analysis procedure is general and can be applied to objects of any geometry. In spite of the presence of many research and review articles [4], a complete description of PO method-based procedure for RCS calculation is rare in literature.  In addition to the complete calculation procedure, this paper also presents a new modified PO method. The PO method approximates the induced surface current on the radar target and uses this current to calculate the scattered field [11, 19]. This approximation leads to discontinuous surface current across the target surface near the shadow boundaries [19]. In the proposed method, a better approximate of the induced surface currents is used to remove the unnatural discontinuities and improve the accuracy in calculation of the scattered field. The proposed method incorporates filtering inspired Fourier transform-based methods to remove the current discontinuity.The paper is arranged as follows: in Section2, the geometrical modeling of the shell-shaped projectile is described. Differential geometry-based definition of the surfaces and normals is derived here. Section 3 contains approximation of the induced surface current on object and formulation of the scattered field and RCS. Numerical results are presented in Section 4. Finally, concluding remarks are given in Section 5.


## 2. Geometrical Modeling

The first step in simulating the RCS of an object is to accurately model the surface of the object. This paper concentrates on the RCS of a shell-shaped object. A shell-shaped object is selected as most projectiles represent this basic shape. A schematic diagram of the object and the three-dimensional coordinate system is shown in Figure1. The object can be modeled by three different canonical surfaces: a half sphere (surface S1), a cylinder (surface S2), and a flat disc (surface S3). Most objects can be similarly modeled by a few common canonical shapes. For this reason, scattering from common shapes such as flat plates, cones, and cylinders has received attention in literature since the 1960s [20, 21].Figure 1
Three dimensional geometry of the shell-shaped projectile, and the coordinate system along with incident field vectors.The origin of the coordinate system is taken at the center of the half sphere, andz axis is taken as the axis of the cylinder. The three surfaces can be expressed by the following equations:
(1)S1:x2+y2+z2=b2,-b≤z≤0,S2:x2+y2=b2,0≤z≤d,S3:x2+y2=ρ2,0≤ρ≤b,z=d.
Here, b= radius of the half-sphere = radius of the cylinder, and d = length of the cylinder. For scattering problems and RCS formulation, it is often more convenient to express the surfaces in differential geometry format rather than coordinate geometry format [11]. To express the surfaces in differential geometry format, the following parameters are used:
(2)S1:x=ρscosϕs,y=ρssinϕs,z=-b2-ρs2,0≤ρs≤b,0≤ϕs≤2π,S2:x=bcosϕs,y=bsinϕs,z=ρs,0≤ρs≤d,0≤ϕs≤2π,S3:x=ρscosϕs,y=ρssinϕs,z=d,0≤ρs≤b,0≤ϕs≤2π.
Here, ρs and ϕs are parameters. Note that the ρs and ϕs for S1 are not the same as the ρs and ϕs for S2 or S3. For each surface the parameters are different and unrelated. Same parameter names are used for simplicity only. Using these parameters, the differential geometric expressions of the surfaces are [22]
(3)S1:r=pscosϕsx^+pssinϕsy^-b2-ρs2z^,S2:r=bcosϕsx^+bsinϕsy^+ρsz^,S3:r=pscosϕsx^+pssinϕsy^+dz^.
The limits of the parameters are given in (2). Using these equations, it is possible to construct a three-dimensional model of the shell-shaped projectile using computer coding. The computer generated model is shown in Figure 2. For computation, b=10cm and d=40 cm are used.Figure 2
Computer-generated three-dimensional model of the shell-shaped projectile.Once the surfaces of the object are defined, it is necessary to define normal vectors on each point of the surface. These normal vectors are necessary for GO-, UTD-, or PO-based scattering formulations [19]. The normal vectors,n^, can be calculated from the differential geometric expression of the surfaces using the following equation [11, 22]:
(4)n^=±(∂r/∂ρs)×(∂r/∂ϕs)|(∂r/∂ρs)×(∂r/∂ϕs)|.
The sign of the normal vectors is selected so that they always point away from the surface. Using (4), the normal vectors on the three surfaces of the object are calculated to be
(5)S1:n^=ρsbcosϕsx^+ρsbsinϕsy^-b2-ρs2bz^,S2:n^=cosϕsx^+sinϕsy^,S3:n^=z^.
Using (5), the normal vectors are plotted over the wire frame of the object using computer coding. The results are shown in Figure 3. From visual inspection it is verified that the normals are perpendicular to the surface and points away from it. This verifies the geometrical modeling performed in this section.Figure 3
Normal vectors (blue lines) on the object surfaces.


## 3. RCS Formulation Using Modified PO Method

To formulate the RCS of the object, the incident field must be defined first. As the radar and the target are usually very far away from each other, the incident field can be modeled as a plane wave, implying that the direction of the wave, the direction of the electric field, and the direction of the magnetic field are perpendicular to each other. With respect to a coordinate system defined at the source point of the wave, if the wave travels atz direction, and the electric field is assumed to be polarized along x direction, then the magnetic field will be polarized along y direction. However, in this case, the coordinate system is defined with respect to the object. So, coordinate transformation must be performed to find the expression of the incident field with respect to the defined coordinate system [23, 24]. For a wave that is incident on the angle (θ, ϕ), an x direction polarized electric field converts to a θ polarized wave [24]. So, the incident electric field can be expressed as [11, 24]
(6)Ei=|Ei|θ^=cosθcosϕx^+cosθsinϕy^+sinθz^.
Here, the amplitude of the incident electric field is assumed to be 1. The direction of the incident ray is given by
(7)s^i=sinθcosϕx^+sinθsinϕy^-cosθz^.
The incident ray direction, along with electric and magnetic field polarizations, is shown in Figure 1.For monostatic RCS calculation, the reflected ray is(8)s^r=-s^i.
PO method uses approximate expression of the surface current density induced on the surface of the object due to the incident field to find the scattered field. The surface current density depends on the incident magnetic field. The incident magnetic field is given by [11, 19]
(9)Hi=1η(s^i×Ei).
Here, η=120π = intrinsic impedance of air. Equation (9) is valid only for plane wave incidence. As the object can be considered to be very far away from the source for most radar applications, this assumption is justified. The PO approximates the surface current density, JPO, as [18]
(10)JPO={2n^×Hi,forilluminatedsurface,0,otherwise.
Equation (10) holds for a perfectly conduction object. For calculating RCS of a metallic object, this approximation is justified. Using (7), (9) along with (5) in (10), it is possible to calculate JPO. However, it is necessary to identify which part of the surface is illuminated by the incident field and which part is not for calculation of (10). This can be accomplished by noticing the angle between s^i and n^. When the vectors are perpendicular to each other, the incident field is tangent on the surface [22, 23]. These surface points indicate the shadow boundary, beyond which the surface points will not be illuminated. So, surface points on which the angle between the normal vector and incident ray vector is greater than 90° are in the shadowed region. This statement can be mathematically expressed as
(11)n^·s^i≤0forilluminatedregion,n^·s^i>0forshadowedregion.
Thus using (10) and (11), the PO approximate of the surface current density can be formulated on S1, S2, and S3 surfaces. It should be mentioned that (11) is true only for objects that only have convex surfaces. If one region of the object creates shadow for another region, then (11) cannot be used to identify the illuminated region and shadowed region. For the relatively simple geometry presented in this paper, (11) is sufficient for calculating the shadow boundaries.FromJPO, the magnetic vector potential, A can be calculated using the radiation integral [12, 19]:
(12)A=∫JPOe-jkRRdS=∫JPOe-jkRRdS1+∫JPOe-jkRRdS2+∫JPOe-jkRRdS3.
Here, R is the distance from the object to the receiver. As monostatic radar cross-section is considered, the source and the receiver are located at the same position.  The position of the source point is expressed in polar coordinates as P(r,θ,ϕ)  with r=1000 meter. R can easily be calculated from the coordinates of the source and the object. The surface integrals can be performed over parametric space. The differential surface element dS can be expressed in terms of the parameters as [22, 23]
(13)dS=|∂r∂ρs×∂r∂ϕs|dρsdϕs.
Evaluating (13) for the three surfaces and using it in (12):
(14)A=∫ϕs=02π∫ρs=0bJPO(ρs,ϕs)e-jkRRρsbb2-ρs2dρsdϕs+∫ϕs=02π∫ρs=0dJPO(ρs,ϕs)e-jkRRbdρsdϕs+∫ϕs=02π∫ρs=0bJPO(ρs,ϕs)e-jkRRρsdρsdϕs.
The scattered field, Es, can be formulated from A as [12]
(15)Es=-j60πλA.
Here, λ = wavelength of the incident field. The RCS, σ of the object is calculated from the scattered field using (6) and (15) [1, 4]:
(16)σ=limR→∞4πR2|Es|2|Ei|2.
The method described so far is the conventional PO method. Although sufficiently accurate, one of the weak points of this method is the approximation of the surface current density, JPO. From (10), it can be seen that JPO abruptly changes to zero in the boundary between illuminated and shadow region. Practical surface current densities do not have this discontinuous nature. This discontinuity arises because PO method does not take into account the creeping waves which exist in the boundary between shadowed region and illuminated region. These waves gradually decrease with distance from the shadow boundaries, and the resulting surface current density decreases gradually.So, to increase the accuracy of the PO method, the discontinuityJPO in must be removed. In this paper, a filtering-based approach is proposed to remove the discontinuity in current distribution. The spatial variation in the current distribution over the surface of the object can be compared with temporal variation of an analog signal in time domain. In Fourier expansion of time signals, the time signal is imagined as superimposition of many sinusoidal signals, termed Fourier spectral components. A sharply varying time signal has high Fourier spectral components. A discontinuous signal has large high Fourier spectral components. Using this analogy, the spatial variation of the surface current in scattering problems can be imagined as superimposition of many surface current components with sinusoidal variation. These components can be isolated using Fourier transformation. It can be imagined that the rapidly varying high components of the surface current are responsible for the sharp discontinuous spatial distribution. The high spectral components add in opposite phase with the other terms in the shadow region to create destructive interference. In absence of these high spectral components, the other components will not completely cancel each other, and therefore there will be an oscillating distribution of surface current in the shadow region. These oscillating distributions can be compared to the creeping waves. Thus filtering out high spectral components should create a smooth distribution of surface current densities which may accurately model the actual current distribution with a higher degree of accuracy.To perform filtering operations to makeJPO continuous, Fourier transform is performed in the parametric space. This produces Fourier domain representation of the PO surface current, JFPO. This can be obtained using the following relation [24]:
(17)JFPO(ρs,ζ)=Fϕs[JPO(ρs,ϕs)]=12π∫-∞∞JPO(ρs,ϕs)e-jζϕsdϕs.
Here, Fourier transform is performed with respect to ϕs, and ζ is another parameter corresponding to the frequency term in conventional Fourier transform of time series data. The complete information of JPO is implicit within JFPO. The discontinuity in JPO with respect to ϕs will result in high JFPO values for large ζ. This is analogous to high frequency terms in discontinuous time series data [25]. As these high ζ components contribute to the discontinuity, removing them should result in a smoother surface current density. How many high ζ value components need to be removed to produce an accurate surface current density cannot be analytically calculated due to the complex nature of the mathematics. Observing the discontinuity in the PO current in many canonical problems found in literature [11] and testing the filter for different parameters, it is found through trial and error that removing 70% of high ζ value components results in a relatively accurate surface current distribution. A modified current density, JSPO, is constructed in ζ space by removing 70% of high ζ value components:
(18)JSPO(ρs,ζ)={0,ζ≥0.7ζmaxJFPO(ρs,ζ),otherwise.Now, the modified smoother current density,JMPO, in ρs-ϕs   space can be obtained by using inverse Fourier transform [24]:
(19)JMPO(ρs,ϕs)=Fζ-1[JSPO(ρs,ζ)]=12π∫-∞∞JPO(ρs,ζ)ejζϕsdζ.
Equation (19) is applicable for canonical geometries and surface of revolution. Now, the discontinuity in JPO with respect to ϕs should no longer be present in  JMPO. However, JMPO may still be discontinuous with respect to ρs.  A similar analysis can be performed to remove this discontinuity by taking Fourier transform of JMPO with respect to ρs, eliminating high components and then taking inverse Fourier transform. Then a smoothed surface current density, JM, can be obtained which is continuous across the surface of the object.It is noted that that as some components of the current densities are filtered out, the energy of the conventional PO current may not be equal to the energy of the modified PO current. This may not be consistent with conservation of energy. To rectify this problem, the filtered current densities must be multiplied by a constant scalar so that the resulting current densities have the same energy as the PO current. The value of the constant can easily be calculated by calculating the energy of the conventional PO current and the unscaled modified PO current. This scaling ensures that no changes have occurred in total energy of the surface currents.Due to harmonic oscillating nature of Fourier transform,JM has oscillating characteristics near the shadow boundary which accurately represents diffraction patterns. Using JM instead of JPO in (12)–(16) will result in a more accurate estimate of RCS.


## 4.  Numerical Results

For numerical analysis, the frequency of the incident field is taken to be 10 GHz. The radius of the sphere and cylinder,b=10 cm, and the length of the cylinder, d=40 cm, are assumed. To perform numerical analysis, the surface of the object must be divided into discrete points. Discrete values of the parameters ρs and ϕs are taken to create discrete points on the surface. Large number of points increase accuracy but also take considerable simulation time. Here, 22 values ρs  per wavelength are considered, and 204 values of ϕs per ρs are considered. These values fall in the range of typical selected values for numerical analysis using PO type method [19].  All numerical results are obtained using computer coding.Figures4 and 5 show the surface current density for two different incident field orientations. These current densities are obtained using conventional PO methods using (10). The discontinuity in the current densities can easily be spotted in both figures. All the current densities are plotted after normalization process. The normalization is performed by dividing the current densities by the maximum value of the current density. Due to this division, the normalized current density is unitless. The resulting normalized current densities are expressed in decibels by taking logarithm.Figure 4
Surface current density using PO method for an incident fieldθ=20°, ϕ=60°.Figure 5
Surface current density using PO method for an incident fieldθ=60°, ϕ=45°.Using the proposed modified PO method, the surface current densities are smoothed. Fourier transform is used with respect to both parameters. The resulting smoothed current densities are shown in Figures6 and 7. The difference in current densities can easily be observed. Comparing Figure 4 with Figure 6 and Figure 5 with Figure 7, it can be seen that the current densities decrease gradually when the proposed method is used whereas the conventional PO method creates sharp discontinuities. The oscillating pattern of the currents near the shadow boundary is created when the proposed method is used. These accurately describe realistic diffracted fields [11].Figure 6
Surface current density using modified PO method for an incident fieldθ=20°, ϕ=60°.Figure 7
Surface current density using modified PO method for an incident fieldθ=60°, ϕ=45°.To verify accuracy of the surface current densities approximated in the proposed filtering method, a simplified canonical problem of scattering from an infinite cylinder is considered. Considering an infinitely long cylinder with axis parallel to thez axis is illuminated by a plane wave. The orientation of the magnetic field of the incident wave is considered along the z axis, and the orientation of the electric field is assumed to be along the x axis. The wave propagates along the y axis. The PO current given by (10) is constant throughout the illuminated region of the cylinder. There will be no variation of the surface current along the z axis. The plane wave will illuminate only half the surface of the cylinder. If an angle, χ, is defined along the xy plane as the angle from the y axis, then for a plane wave propagating along the y axis, the illuminated region spans from χ=-90° to 90°. The PO current abruptly falls to zero outside this angle. The proposed filtering method adjusts this discontinuity and makes the current smooth. The actual angular distribution of the current density can be formulated using a Method of Moments (MoM), which is a well-known benchmark method [19].  The current distribution obtained from the PO method, the proposed modified PO method, and MoM are shown in Figure 8. It can be clearly observed that the proposed method produces a current distribution which resembles the results obtained from MoM much closely compared to the PO method. This higher degree of accuracy in estimating the surface current ensures that the proposed method will generate more accurate RCS results compared to the conventional PO method.Figure 8
Surface current densities on an infinite cylinder surface atf=3GHz.The RCS of the shell-shaped object is formulated using both conventional PO and the modified PO method. The RCS atϕ=0° plane as a function of zenith angle, θ, is shown in Figure 9.Figure 9
Monostatic RCS of the shell-shaped projectile at 10 GHz as a function ofθ.RCS is expressed in dB with respect to 1 m2 area. This unit is represented as dB m2 or dB sm [1, 4].  As the object is circularly symmetric around the z axis, the scattered field and RCS are independent of the azimuth angle,ϕ [11, 26]. It is seen that both methods give similar results in most angular regions. The variation in result comes for angular regions where the spherical surface is illuminated. This is expected as the spherical surfaces are affected by the diffracted rays and creeping waves more than the other shaped surfaces [11]. Due to the better approximate of the surface current, the results obtained from the proposed method are more accurate. It is noted that at θ=180° only the half-sphere surface is illuminated. The RCS at this angle is around −15 dB sm. The projected area of the half sphere =πb2=0.0314m2=-15dBsm, which is equal to the RCS. This is expected for a spherical shaped object. This consistency verifies the numerical analysis.The RCS of the projectile at 5.5 GHz as a functionθ is shown in Figure 10. The results from the proposed method and MoM are shown in the same plot for comparison. It can be observed that the results obtained from the proposed method matches with the results obtained from MoM closely for most values of θ.Figure 10
Monostatic RCS of the shell-shaped projectile at 5.5 GHz as a function ofθ.For the angles where RCS is high, the results of both methods are almost identical. The deviation arises for angles where RCS value is low. The higher degree of accuracy obtained from MoM comes at a cost of complex computation and extensive simulation time. The errors produced by the proposed method are relatively small and acceptable for many applications. In most cases the accuracy from PO method is sufficient, and the proposed method is expected to be more accurate than the PO method.The RCS of the shell-shaped projectile for frequencies 3 GHz, 6 GHz, 10 GHz, and 14 GHz as a functionθ is shown in Figure 11. The obtained RCS pattern is similar to the RCS patterns of objects of similar shape described in [4, 16].Monostatic RCS of the projectile at 3 GHz, 6 GHz, 10 GHz, and 14 GHz as a function ofθ.
(a)
(b)
(c)
(d)The computer simulation was performed on an Intel Core i5-2430M 2.4 GHz CPU with 2.94 GB usable RAM. The simulation time of conventional PO and modified PO is compared in Table1.Table 1
Comparison of CPU simulation time.
Frequency (GHz)
Simulation time for PO method (sec)
Simulation time for modified PO method (sec)
3
16.83
26.09
6
32.09
50.79
10
53.19
87.01
14
74.15
144.23From Table1, it can be seen that for lower frequency, the simulation time of conventional PO and the proposed modified PO is comparable. For higher frequency, the modified PO method takes 80% to 90% longer time to simulation compared to the conventional PO method. This excess simulation time can be expressed by the fact that for high-frequency simulation, the number of discrete surface points selected is larger. As there are 22 values ρs selected per wavelength and the wavelength is smaller for high-frequency simulation, the overall number of point increases. The filtering of the surface currents requires this additional simulation time due to large number of points. However, the proposed method is much faster compared to Method of Moments. For each simulation, the Method of Moments requires over 500 seconds of simulation time.


## 5. Conclusion

RCS formulation procedure using PO method is described in detail in this paper. The paper covers geometrical modeling of shell-shaped projectile using differential geometry, formulation of surface current density, evaluation of the scattered field integral in parametric space, and monostatic RCS formulation. A modified PO method is presented which approximates surface current density more accurately. The additional computational steps require Fourier transform and inverse Fourier transform only. These can be easily incorporated in computer code using well-known Fast Fourier Transform (FFT) algorithm. Thus the accuracy is increased without significant increase in computational complexity. The obtained results using the modified PO method are consistent with similar results found in literature.



---

*Source: 328321-2012-09-18.xml*