RCS Measurement and Analysis of Rectangular and Circular Cross-section Cavities

Similar documents
Lecture 7: Introduction to HFSS-IE

Progress In Electromagnetics Research M, Vol. 20, 29 42, 2011

HFSS Hybrid Finite Element and Integral Equation Solver for Large Scale Electromagnetic Design and Simulation

Progress In Electromagnetics Research, PIER 43, , 2003

Accurate Radar Cross Section Modelling of Jet Inlets & Engines

Aspects of RF Simulation and Analysis Software Methods. David Carpenter. Remcom. B = t. D t. Remcom (Europe)

A Graphical User Interface (GUI) for Two-Dimensional Electromagnetic Scattering Problems

Comparison of iteration convergences of SIE and VSIE for solving electromagnetic scattering problems for coated objects

A MATLAB PHYSICAL OPTICS RCS PREDICTION CODE

RCS CALCUATIONS USING THE PHYSICAL OPTICS CODES

Simulation Advances for RF, Microwave and Antenna Applications

Numerical Simulation of RCS for A Spherical Convergent Flap Nozzle with A Non-rectangular Divergent Duct

For functionality and CAD/EDA import filter, see technical specifications of the CST STUDIO SUITE

Comparison of TLM and FDTD Methods in RCS Estimation

ABSORBER FOAM CHARACTERIZATION FOR PREDICTING OVERALL ANECHOIC CHAMBER PERFORMANCE

EM Software & Systems GmbH

A Magnetic Field Integral Equation Based Iterative Solver for Scattered Field Prediction

2

Acceleration of the Iterative Physical Optics Using Graphic Processing Unit. Abstract

IMPLEMENTATION OF ANALYTICAL (MATLAB) AND NUMERICAL (HFSS) SOLUTIONS ADVANCED ELECTROMAGNETIC THEORY SOHAIB SAADAT AFRIDI HAMMAD BUTT ZUNNURAIN AHMAD

Simulation Advances. Antenna Applications

Design and Analysis of an Euler Transformation Algorithm Applied to Full-Polarimetric ISAR Imagery

ELECTROMAGNETIC diffraction by perfectly conducting

Advanced Surface Based MoM Techniques for Packaging and Interconnect Analysis

Defense Technical Information Center Compilation Part Notice

Modeling the Effects of Wind Turbines on Radar Returns

Simulating Reflector Antenna Performance with GRASP9

CHAPTER 6 MICROSTRIP RECTANGULAR PATCH ARRAY WITH FINITE GROUND PLANE EFFECTS

Full-Wave and Multi-GTD Analysis of the Ice Cloud Imager for MetOp-SG

HFSS Ansys ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary

Reflection, Refraction and Polarization of Light

Development and Applications of an Interferometric Ground-Based SAR System

Physics 11. Unit 8 Geometric Optics Part 1

Terahertz backscattering behavior of various absorbing materials

Bistatic SAR coherence improvement through spatially variant polarimetry

Modeling the Radar Signature of Thin Metallic Objects with the AFDTD Software

Insights into EMC Chamber Design:

Stable Laser Resonator Modeling: Mesh Parameter Determination and Empty Cavity Modeling

newfasant PO Multiple effects analysis

25 The vibration spiral

A systematic approach to the analysis of polarisation-sensitive reflectors

D&S Technical Note 09-2 D&S A Proposed Correction to Reflectance Measurements of Profiled Surfaces. Introduction

Experiment 8 Wave Optics

Exterior Orientation Parameters

AUTOMATIC TARGET RECOGNITION IN HIGH RESOLUTION SAR IMAGE BASED ON BACKSCATTERING MODEL

The sources must be coherent. This means they emit waves with a constant phase with respect to each other.

High Frequency Wave Propagation and Discrete Geodesics

INTRODUCTION TO The Uniform Geometrical Theory of Diffraction

An Efficient Analysis on the Fitting Error Caused by the Deformation of Metal Pylon in the RCS Measurement

Lecture 2: Introduction

MODELLING OF THE SCATTERING BY A SMOOTH DIELECTRIC CYLINDER: STUDY OF THE COMPLEX SCATTERING MATRIX

Reflection and Refraction

REDESIGN AND OPTIMIZATION OF THE PAVING AL- GORITHM APPLIED TO ELECTROMAGNETIC TOOLS (INVITED PAPER)

ERMES: NUMERICAL TOOL FOR SAR COMPUTATIONS

A fast direct solver for high frequency scattering from a cavity in two dimensions

SCATTERING AND IMAGE SIMULATION FOR RECON- STRUCTION OF 3D PEC OBJECTS CONCEALED IN A CLOSED DIELECTRIC BOX

Advanced modelling of gratings in VirtualLab software. Site Zhang, development engineer Lignt Trans

IMAGING WITH SYNTHETIC APERTURE RADAR

M-SPURT Compressing the target characterization for a fast monostatic RCS simulation

On Wire-Grid Representation of Solid Metallic Surfaces

E x Direction of Propagation. y B y

Diffraction. Single-slit diffraction. Diffraction by a circular aperture. Chapter 38. In the forward direction, the intensity is maximal.

Modeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics

Design of Electromagnetic Test Sites

Simulation of Transition Radiation from a flat target using CST particle studio.

Enhanced Characteristic Basis Function Method for Solving the Monostatic Radar Cross Section of Conducting Targets

Study of the Effects of Target Geometry on Synthetic Aperture Radar Images using Simulation Studies

HFSS 14 Update for SI and RF Applications Markus Kopp Product Manager, Electronics ANSYS, Inc.

Development and assessment of a complete ATR algorithm based on ISAR Euler imagery

Validation of aspects of BeamTool

Memorandum. Clint Slatton Prof. Brian Evans Term project idea for Multidimensional Signal Processing (EE381k)

Mie scattering off plasmonic nanoparticle

EO-1 Stray Light Analysis Report No. 3

Polarization Properties of ACS

Condenser Optics for Dark Field X-Ray Microscopy

Introduction to the FEKO Suite

DESIGNER S NOTEBOOK Proximity Detection and Link Budget By Tom Dunn July 2011

Three-Dimensional Laser Scanner. Field Evaluation Specifications

ARTIFICIAL SCATTERERS FOR S.A.R. INTERFEROMETRY

THE concept of using a lossy material to absorb an

Guidelines for proper use of Plate elements

Spherical Geometry Selection Used for Error Evaluation

Designing Horn Antenna utilizing FEM Symmetry Boundary Conditions

NEW BENCHMARK RADAR TARGETS FOR SCATTERING ANALYSIS AND ELECTROMAGNETIC SOFTWARE VALIDATION

GPR DATA ANALYSIS OF AIRFIELD RUNWAY USING MIGRATION VELOCITY SIMULATION ALGORITHM

Unique Two-Way Field Probe Concept Utilizing a Geodesic Sphere and Quad-Rotor

Using the Discrete Ordinates Radiation Model

Aradar target is characterized by its radar cross section

Iterative methods for use with the Fast Multipole Method

Deterministic Approach to Achieve Broadband Polarization-Independent Diffusive Scatterings Based on Metasurfaces

Electromagnetics. R14 Update. Greg Pitner ANSYS, Inc. February 24, 2012

Mu lt i s p e c t r a l

Can Xia 1, 2,WanqingYou 2, and Yufa Sun 2, *

Plane Wave Imaging Using Phased Array Arno Volker 1

High-Frequency Algorithmic Advances in EM Tools for Signal Integrity Part 1. electromagnetic. (EM) simulation. tool of the practic-

HFSS PO Hybrid Region

HFSS - Antennas, Arrays and FSS's. David Perry Applications Engineer Ansoft Corporation

SESSION 5: INVESTIGATING LIGHT. Key Concepts. X-planation. Physical Sciences Grade In this session we:

3 - SYNTHETIC APERTURE RADAR (SAR) SUMMARY David Sandwell, SIO 239, January, 2008

Modeling the Transmission Loss of Passthroughs in Sound Package using Foam Finite Elements

Transcription:

RCS Measurement and Analysis of Rectangular and Circular Cross-section Cavities Abhinav Bharat, M L Meena, S. Sunil Kumar, Neha Sangwa, Shyam Rankawat Defence Laboratory, DRDO Jodhpur, India-342011 Abstract The RCS of re-entrant structures, like cavities, is difficult to estimate either using High-Frequency Computational Electromagnetics (CEM) solvers with multibounce capability, or with exact solvers. Results from High- Frequency solvers are not consistent amongst themselves, or with measured data. Exact solvers, such as MoM, require heavy compute resources, and cavity structures make the convergence slower. In this work, the RCS of a simple circular crosssection duct and a simple rectangular cross-section duct has been investigated using both computations and measurements. The ducts have been fabricated and embedded into a low-rcs testbody, so that only the ductcontributions are measured. The results of the measurements, carried out at DRDO's NFDRF RCS Range, are presented. The RCS of the ducts were computed using an SBR Solver, and using an MoM solver. The measured data is compared with the results of the computations. It is shown that the High-Frequency Solver is inadequate to correctly predict the RCS behavior of ducts, even for very simple duct shapes. Keywords- RCS, SBR, MOM I. INTRODUCTION In real-world electromagnetic problems, like scattering, radiation, etc. for electrically large objects with cavity are not analytically calculable. Cavity structures, such as air-intakes and exhausts of an aircraft, give rise to large backscatter echo over a wide angular sector, and therefore, require critical attention in the design of low-observable platforms. The measurement of radar echo from real target is costly and complex equipments are involved in measurement setup. For such setup, precise positioning of antennas, fairly clean (low background clutter and multipath) measurement facility including LO pylon, and accurate angular scan become mandatory. In addition, the microwave instrumentation should be configured to illuminate the object with sufficient energy so that the returned signal from test target has sufficient isolation w.r.t., the background clutter. The measurement of complex cavity shapes further requires a suitable setup such that only the interior portion of cavity is irradiated with Electromagnetic interaction. So, it becomes necessary to hide the cavity outer surface from EM interaction. This can be performed either by appropriate placement of absorbers around the cavity periphery or by blending the cavity structure on a low RCS testbody, e.g., Almond testbody. Thus, an accurate characterization of the complex cavity structures radar scattering is a rather challenging task. CEM solvers have been in use for modeling the interaction of electromagnetic fields with physical objects and related scattering. Several numerical techniques are present which can derive the results by solving the closed form of Maxwell s equations. An effort into analysis of cavity shapes has been taken up by applying various numerical techniques for such analysis. The techniques though very well suited for single bounce analysis have been put to test for the multi-bounce phenomenon involved in our scenario. This paper presents the monostatic RCS study in azimuth plane for Rectangular and Circular crosssection cavities. Both the computational and measurement techniques have been used for RCS analysis, exclusively at an incident frequency of 10GHz. The computation results have been generated for the study performed using a High Frequency solver (PO- PTD), a full-wave solver (MLFMM) and through Nearfield RCS measurements. II. SIMULATION AND MEASUREMENT TECHNIQUES UTILIZED The solution to the Electromagnetic scattering problems requires an appropriate solution of the Maxwell's Equations. For the computational techniques involving surface meshes, the surface of the object is discretized into small elements (triangles or quadrilaterals). This mesh is then utilized for the computation of induced Electric and Magnetic fields (or

induced charges and currents) w.r.t. the incident Electromagnetic waves along with application of suitable boundary conditions at the interfaces. Here, we have used two solvers, (i) Asymptotic solver based on PO- PTD with multi-bounce capability, (ii) Full-wave solver based on MLFMM technique. Asymptotic Solver (PO-PTD with multi-bounce): The Asymptotic computation technique is an analysis in the frequency domain based on use of different high frequency approximation methods like GO, GO-GTD, PO, PO-PTD, hybrid SBR. The incident and reflected rays can be considered either independent or bundled together within the ray tubes. This solver is typically used for scattering computations on electrically very large objects which are difficult to handle by other EM solution methods. The mesh generation is very robust and quite insensitive to the quality of the CAD model. The approximation involved is the use of ray optics approach to estimate the induced surface currents and hence, to calculate the scattered fields. Full-wave Solver (MLFMM/ MoM): The full-wave solver based on integral equation solution to Maxwell's Equations, is of special interest for electrically large objects. Here, the discretization of the calculation domain is reduced to the object boundaries (surface mesh), which leads to the solution of a linear equation system with less unknowns than the volumetric techniques. The integral equation solver based on MoM technique creates a dense matrix, as each element is considered coupled to all other discretized elements of the object. The other technique, Fast Multipole Method (FMM), further uses the sub-division of computation domain to boxes so as to combine the corresponding coupling sparse matrices. Another technique, Multi- Level Fast Multipole Method (MLFMM), further uses a recursive sub-division scheme to speed up this matrix vector multiplication process. So, the MLFMM solver is more suitable for analysis of very large electrical size problems as the standard MoM approach becomes cumbersome. It speeds up the matrix vector multiplications, which are dominating operations in the iterative solver and also enhances the memory efficiency of computation platform. The solver time for MLFMM is estimated as N iter. N.log(N) while for MoM it is approximately N 3. Similarly, the memory requirement for MLFMM is approximately N.log(N) as compared to N 2 for MoM. Measurement Method: EM waves are transmitted towards the target and the received signal can be used to determine the shape and composition of target. SAR and ISAR techniques are used for collection of this backscattered echo. In our case, at NFDRF range, DRALA complex, DL Jodhpur, ISAR technique has been employed for RCS data acquisition through measurements. In ISAR technique, the target to be analyzed is rotated over the angular sector of interest while the radar is kept stationary. The reflectivity data is collected at each angular sample over the whole frequency band. The target of interest can be considered as comprising of various reflecting centers called scatterers and the total backscattered echo at the receiver T field is the vector (phase and magnitude) sum of the reflectivity data received from these scatterers: j2π ( ) ( 2d( x, y) ) T field = κ ρ x, y exp dxdy (3) λ area where, κ = a proportionality constant, ρ(x,y) = 2D spatial distribution of reflectivity, d(x,y) = distance of reflecting centre to antenna aperture, λ = wavelength of incident wave, The total field thus depends on the wavelength of incident wave, reflectivity of point scatterers and the distance between observation point and the point scatterer. This field data at the receiver is then calibrated to estimate the expected RCS levels. III. TESTBED AND CAVITY SHAPES An Almond shape (low RCS body) [1] of length 125cm was designed and developed as a testbed for performing RCS computation and measurement studies on rectangular and circular cross-section cavities as shown in Figure 1. This testbed was identified as having very low signature in the frontal sector and selected for performing measurements of backscatter echo from the interior of the cavities, while removing scattering from its exterior surface. The only source of backscatter, that was observed from the testbed (within X-Band), is the tip-diffraction effect from its nose. Figure 1 Fabricated Almond-shaped testbed The open ended rectangular and circular crosssection cavities (dimensions described in Figure 2 &

Figure 3, respectively) were then blended on to the low RCS testbed (Figure 4 & Figure 5). The cavities are placed inside the almond testbed with their axis placed horizontally w.r.t. the ground. The cavity openings were also blended with the surface of almond testbed. The fabrication of the measurement models was done using FRP, the surface of which was subsequently metallized using conductive paint coating. Figure 4 Rectangular cavity blended in low-rcs testbed Figure 2 The dimension of the rectangular cross section is 15 cm x 25 cm. The cavity lower and upper faces length is about 74 cm and 43 cm. These lengths are chosen so as to fit in a 125 cm almond shaped testbed Figure 3 The diameter is 150 cm. The maximum length is about 83 cm and minimum length is about 49 cm RCS analysis of the models was carried out at X-Band (10 GHz) by means of both computations and measurements. At this frequency, the electrical size of the rectangular cavity aperture (distance between parallel plates) for azimuth plane scan is approximately 10λ. So, the ray bouncing approach is expected to yield adequate similarity to the full-wave solver results for this high frequency analysis. But, the diameter of circular crosssection cavity aperture is approximately 5λ. So, ray optical behavior assumption may not be suitable for analysis of this cavity shape, as suggested by the available literature. Figure 5 Cylindrical cavity blended in low-rcs testbed IV. TESTBED VERIFICATION A verification study of the blending of rectangular and circular cross-section cavities within the low-rcs Almond testbed is as shown in Figure 6. The plot shows measured RCS levels vs. aspect angle in azimuth plane at 10 GHz incident frequency, horizontal polarization, for almond testbed and almond with rectangular & circular cross-section cavities. This plot clearly suggests that the RCS levels of almond shaped testbed is quite low (below -20 dbsm) in the angular window of Nose-on ± 40 deg. So the cavity shapes can be blended on such platform for RCS analysis within this angular window. Figure 6 RCS plot of Almond testbed and blending with Rectangular and Cylindrical Cavity- Measured results, 10 GHz,

V. RESULTS & OBSERVATIONS A comprehensive study was performed to analyze the RCS response as a function of aspect angle in azimuth plane for three models: (i) simple almond shaped testbed of length 125cm, (ii) rectangular cross-section cavity blended on the almond shaped testbed, (iii) circular cross-section cavity blended on the almond shaped testbed. The study shown here, is mainly focused on analysis of the scattering response at 10 GHz and horizontal polarization. Data was collected from computations carried out using commercially available Asymptotic solver and MLFMM solver and also through measurements at NFDRF range. The results of analysis of low RCS testbed, rectangular and circular cross-section cavity blended on testbed using different solvers and measurements are shown in Figure 7, Figure 8 and Figure 9, respectively. Figure 8 RCS comparison plot of rectangular cross-section cavity blended on Almond shaped testbed for Asymptotic, Integral Equation(MLFMM) solver and measurement results, 10 GHz, Figure 7 RCS comparison plot of Almond shaped testbed for Asymptotic, Integral Equation(MLFMM) solver and measurement results, 10 GHz, Figure 8 shows similar RCS comparison results for rectangular cross-section cavity blended on almond shaped testbed. For Nose-on incidence, if only the flat plate present at the cavity backing is assumed dominant, then the observed RCS level should be about +12.93 dbsm. However, such Nose-on incidence RCS level is observed only for computation results of Asymptotic solver. The results for full-wave solver computations and measurements are slightly lower in magnitude and have good correlation between them. This reflects the approximate nature of Asymptotic solver, whose accuracy decreases with the decrease in electrical size of the cavity aperture. However, considering the good correlation observed with other results, the solver can still be used as a time efficient alternative. As observed in Figure 7, a low-rcs angular window is observed for the almond shape testbed within Nose-on ± 40 deg span for all the results. The measured results are observed to be slightly higher as compared to the computation results. The computation results from the full-wave solver are more closer to the measurement results. This difference in the RCS levels can be attributed due to: (i) dominant tip diffraction effects in the testbed within the low-rcs angular window, (ii) approximate nature of the Asymptotic solver, (iii) manufacturing deformities.

The RCS study was then performed on these models by means of both measurements and computational tools: both approximate and exact solvers. Figure 9 RCS comparison plot of circular cross-section cavity blended on Almond shaped testbed for Asymptotic, Integral Equation(MLFMM) solver and measurement results, 10 GHz, Similarly, Figure 9 shows the RCS comparison plots for circular cross-section cavity blended on almond shaped testbed. For Nose-on incidence, if only the flat plate present at the cavity backing is assumed dominant, then the observed RCS level in this case should be about +6.40 dbsm. However, only the Asymptotic solver is seen to show such Nose-on incidence RCS level. The results for full-wave solver computations and measurements are slightly lower in magnitude and have better correlation between them as compared to the Asymptotic solver. This RCS plot demonstrates the effect of reduction in RCS levels due to partial blocking of Electromagnetic waves at the cavity aperture. The asymptotic solver is seen not able to properly capture the RCS response over the angular sector. The full-wave solver is although observed more close to the measurement results. However, slight deviations are expected from measurement results possibly due to inaccurate positioning of the almond as well as cavity back-plate in both pitch and tilt aspects. VI. CONCLUSIONS The cavity structures are observed as prime back-scatterers on aircraft platforms while present as airintakes and exhaust cavities. They show a significant increase in the observed RCS levels over a wide angle sector both in azimuth and pitch planes when blended on such platforms. In order to assess the RCS contribution of typical cavity shapes due to interior irradiation with Electromagnetic waves while ignoring the contribution from its outer surface, the cavity shapes were blended on a specially designed models of low-rcs testbed, which were also subsequently fabricated. The observations are as mentioned below: i. The low-rcs almond testbed can be used to study the cavity RCS within the Nose-on ± 40 deg angular sector in azimuth plane ii. The full-wave solver simulations should obviously be the choice for cavity RCS estimation however, owing to multi-bounce phenomenon involved, they show lesser rate of convergence and hence are not time efficient. iii. For the cavity aperture diameter exceeding 10λ size, the asymptotic solver is expected to give good correlation with the actual results. So they can be preferred for analysis of electrically very large cavity apertures as they are time efficient. iv. The positioning of the almond shaped testbed is very crucial for both pitch and tilt aspects in order to correctly predict the RCS of blended cavity shape. ACKNOWLEDGEMENTS The authors would like to acknowledge ADA, Bangalore team for their support to this activity. The support received from Director, DLJ, Dr SR Vadera, and the kind guidance received from Group Head, Ms Anjali Bhatia, is acknowledged. We further extend our sincere gratitude to the support staff at NFDRF Range, DL Jodhpur. REFERENCES 1. 'US Patent No. 4809003 for Almond Test Body', Dominek et al, dated 28 Feb 1989 2. High Frequency RCS of Open Cavities with Rectangular and Circular Cross Sections, Ling, Lee, Chou, IEEE Trans. AP-37, May 1989 3. 'Shooting and Bouncing Rays: Calculating the RCS of an Arbitrarily Shaped Cavity', 'Hao Ling, Ri Chee Chou, Shung Wu Lee, IEEE Trans. AP-37, Feb 1989' 4. RCS Study of Air Intake Configuration of AURA Interim Technical Progress Report, Defence Lab Jodhpur, Feb 2012