Next Article in Journal
Enzyme Immobilization on Polymer Membranes: A Quantum and Molecular Mechanics Study
Previous Article in Journal
Fuzzy C-Means Based Clustering and Rule Formation Approach for Classification of Bearing Faults Using Discrete Wavelet Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computation and Experiment on Linearly and Circularly Polarized Electromagnetic Wave Backscattering by Corner Reflectors in an Anechoic Chamber

1
Center for Environmental Remote Sensing, Graduate School of Science and Engineering, Chiba University, Chiba 263-8522, Japan
2
Department of Informatics, Universitas Pembangunan Jaya, Tangerang Selatan 15413, Indonesia
*
Author to whom correspondence should be addressed.
Submission received: 22 August 2019 / Revised: 18 September 2019 / Accepted: 21 September 2019 / Published: 24 September 2019
(This article belongs to the Section Computational Engineering)

Abstract

:
Electromagnetic wave backscattering by corner reflectors in an anechoic chamber is studied using our developed computational tool. The tool applies the Finite-Difference Time-Domain (FDTD) method to simulate the propagation of the wave’s electric and magnetic fields. Experimental measurement in an anechoic chamber is also carried out as a comparison. The two results show agreement, including the finding that the backscatter intensity variation amongst the four circularly polarized modes is significantly smaller than the variation amongst the four linearly polarization modes.

1. Introduction

Understanding electromagnetic (EM) wave backscattering by corner reflectors is of interest to radar researchers. In particular, the backscattering of circularly polarized (CP) EM waves is of great interest. To study these CP EM waves, computations can be useful.
The articles [1,2,3,4,5,6,7,8,9] report previously done computations on EM wave distribution in an anechoic chamber, but not in the context of backscattering by corner reflectors or in the context of CP EM waves. These works contributed to the development and performance improvement of the anechoic chamber. The articles [10,11,12,13,14,15,16] report computations and analyses on EM wave backscattering by corner reflectors. However, most of the reported computations for EM waves were done for non-polarized and linearly polarized waves. Until today, only [10] reports a computation on circularly polarized EM waves. It used the Method of Moments, and the report did not present any technical details about how the circular polarization was realized.
To respond to this situation, our research aims to study circularly polarized EM wave backscattering by corner reflectors in an anechoic chamber through a computation using the Finite-Difference Time-Domain (FDTD) method. There has been no such work reported thus far. As an addition, we did the same computation for linearly polarized (LP) EM waves. An equivalent experiment has also been carried out to obtain a comparison.
The simulation was done using our developed computational tool, called Anechoic Chamber Simulation Tool Version 1 (AST_V1). Both the simulation and the experiment were done for four target types, i.e., triangular trihedral, horizontal plate bar, vertical plate bar, and sphere, and for all linear polarization and all circular polarization modes.
The remainder of this article is arranged as follows: Section 2 addresses the method, Section 3 reports the results, Section 4 discusses the results, and Section 5 summarizes the content of all sections.

2. Method

2.1. The Measurement Plan

The measurement plan is depicted in Figure 1. This was the scenario to be followed by both the computation and the experiment. It can be described as follows:
  • The cuboid anechoic chamber size shall be 6.5 (length) × 4.0 (width) × 2.4 m (height).
  • The distance between the transmitter (TX) and the receiver (RX) antennas shall be 0.8 m.
  • The distance between the target (corner reflector) and antennas in the orthogonal direction shall be 3.0 m.
  • The vector network analyzer (VNA)’s port 1 shall feed signals to the TX antenna.
  • The RX antenna shall feed signals to the VNA’s port 2.
  • The turntable shall be controlled by the turntable controller.
  • Both the VNA and the turntable controller can be driven by the Personal Computer (PC).
  • Phase shifters shall be provided for the TX and RX antennas in the case of circular polarization and shall be installed correctly for each mode, i.e., RX right-handed TX right-handed (RR), RX left-handed TX right-handed (LR), RX left-handed TX right-handed (LR), or RX left-handed TX left-handed (LL). Phase shifters will not be needed for linear polarization modes, i.e., RX horizontal TX horizontal (HH), RX vertical TX horizontal (VH), RX horizontal TX vertical (HV), and RX vertical TX vertical (VV).
  • The transmit power shall be 5 dBm for all measurement cases.
  • The target shall be placed on a turntable in such a way so that elevation angle θ = 0°, and be rotated in such a way so that the aspect angle ϕ changes from −90° to 90° (negative and positive signs respectively represent the left and the right sides of the target) with 5° increments.
  • The EM wave frequency shall be chosen in such a way so that 5λ < target size < 10λ.
Definitions of the elevation angle θ and the aspect angle ϕ are explained in Figure 2. The targets were corner reflectors of four different shapes: Trihedral, horizontal bar, vertical bar, and sphere. To make it realistic for practices, every target was chosen to be large enough compared to the wavelength but not too large to be carried in a radar mission. Here we use the term “medium size” to specify that 5λ < target size < 10λ. The targets and their descriptions are listed in Table 1. The horizontal and the vertical plate bars are simply the same object, only having different standing orientations.

2.2. The Method of Computation

The computation was done systematically to fulfill the measurement plan specified in Section 2.1. A project definition, software requirement review, preliminary design review, and critical design review were carried out before coding and implementation [17]. The critical design review (CDR) of the AST_V1 is presented in Figure 3 and Figure 4. The tool requires default input data, including the permittivity and the permeability of free space, and the relative permittivity (dielectric constant) and the relative permeability of the air (because the transmission medium in the experiments is the air). However, these values can be approximated by 1. The tool needs simulation case-related data, including the anechoic chamber size; the type, size, and material properties of corner reflector; the type and the material properties of the antennas; the distance between the TX and RX antennas; the distance between the antennas and the corner reflector; the polarization mode (either HH, VH, HV, VV, RR, LR, RL, or LL); and the Nyquist factor. The larger the Nyquist factor chosen by the user, the smaller the FDTD cell size configured by the tool will be, and thus the more accurate the simulation results can be expected to be while more computer resources will be consumed.
The field distribution was simulated using the FDTD method by applying the Yee-grid configuration [18,19], employing an electric- and magnetic-fields update technique based on Maxwell equations [18,19,20,21,22], and realizing the energy-absorbing walls of the anechoic chamber with perfectly matched layers (PMLs) [20]. The real size of the anechoic chamber is 6.5 (length) × 4.0 (width) × 2.4 m (height). However, to minimize the usage of computer resources, the computation space domain was reduced to 4.2 (length) × 2.2 (width) × 2.0 m (height). The cuboid space domain was chosen in such a way that it is small enough but still contains all the objects of interest, i.e., the TX antenna, the RX antenna, and the corner reflector. This arrangement is valid because EM wave paths between the TX antenna and the corner reflector, and between the corner reflector and the RX antenna are “line of sight”. This computation space domain was realized by a Yee-grid consisting of 337 × 177 × 161 (9.6 million) cells, where the cell size was 1.25 × 1.25 × 1.25 cm. With a lambda of 10 cm (f = 3 GHz), this cell size corresponds to a Nyquist rate of 8, which is well above the minimum Nyquist rate, i.e., 2.
In the computation, before the FDTD iteration, several procedures were carried out. The TX and RX antennas were synthesized, and targets were also synthesized and rotated according to the requested aspect angle ϕ. The transmitted power was converted to the transmitted voltage, transmitted wave intensity, and transmitted electric field. The transmitted electric field was represented in sinusoidal waves and fed carefully to the TX antenna by considering the requested mode of polarization.
For each mode of linear polarization, i.e., RX Horizontal–TX Horizontal (HH), RX Vertical–TX Horizontal (VH), RX Horizontal–TX Vertical (HV), or RX Vertical–TX Vertical (VV), either the TX horizontal or TX vertical antenna was fed with the sinusoidal waves. For each mode of circular polarization, i.e., RR, LR, RL, or LL, more attention was required since both the TX horizontal and TX vertical antenna needed to be fed simultaneously with the same train of sinusoidal waves, but of different phases.
During the iteration of the FDTD calculation, the value of the electric field arriving at the RX antenna was stored at every step. Upon completion of the calculation iteration, the set of stored electric field values was converted to the received wave intensity and optionally received power. The backscattered wave intensity in dB was calculated by:
BS   Intensity   ( dB )   =   10   log ( I r I t )
where It is the intensity of the transmitted wave and Ir is the intensity of the received wave. Since the intensity equals the voltage amplitude squared, while the voltage is proportional to the electric field, the backscattered wave intensity can be alternatively obtained by:
BS   Intensity   ( dB )   =   20   log ( E r E t )
where Et is the electric field transmitted by the TX antenna and Er is the electric field received by the RX antenna.
The software has now become a graphical user interface (GUI) tool, namely AST_V1. The tool is currently capable of handling simulations on EM wave backscattering in an anechoic chamber with a size range of between 2 × 2 × 2 m to 20 × 20 × 20m in the frequency range of 3 to 8 GHz, with a full circle rotation turntable for the device under test (DUT), which has two horizontal dipole antennas and two vertical dipole antennas.
The AST_V1 user interface’s appearance is shown in Figure 5. The four target types synthesized by the tool are shown in Figure 6. Some visualization frames during a simulation are presented in Figure 7, where target type = trihedral, f = 3 GHz, polarization mode = RR, θ = 0°, and ϕ = 0°. In this figure, the target structure, the antenna structure, and the electric field distribution within the cuboid 3D space are visualized by applying one horizontal cutting plane and one or two vertical cutting plane(s).
It is necessary to mention that every cuboid shown in the figures visualizes a 3D space that is smaller than the anechoic chamber, because it represents only the FDTD space domain. This space domain reduction technique is applied to reduce computation time. The computation was done with a personal computer that uses an Intel CORE i7 processor and a set of 32 GB random access memory (RAM).
For future operations, in addition to the usage of dipole antennas, the tool can be equipped with more antenna choices, e.g., the circularly polarized microstrip antenna [23] that can handle a frequency range of 4.83 to 5.71 GHz. If the tool is modified to handle lower-frequency operations, it can be equipped with lower-frequency antennas, e.g., the circularly polarized microstrip antenna [24] that can handle frequencies of around 2.2 GHz. If it is intended to do simulations for higher frequencies, it can be equipped with, for example, the circularly polarized microstrip [25] that can operate in frequencies of around 9.4 GHz.

2.3. The Method of Experiment

The experiment was carried out thoroughly to fulfill the measurement plan specified in Section 2.1. It was done in the following sequence:
  • The target, the TX antenna, and the RX antenna were placed properly according to Figure 1. Each TX and RX antenna set consisted of a horizontal antenna and a vertical antenna.
  • Two phase shifters were provided to be used in the cases of circular polarization only (RR, LR, LR, LL). One phase shifter was used for the TX antenna, and another phase shifter was used for the RX antenna.
  • Before the measurements, the VNA conducted a calibration involving the measurement of the loss and time delay contributed by extended coax cables and connectors.
  • The VNA conducted the whole measurement cases, where
    • The transmit power was always set to 5 dBm.
    • The measurement mode was always set to Backscattering/Inverse Synthetic Aperture Radar (ISAR).
    • Equal measurements were done for the corner reflectors with shapes of sphere, trihedral, horizontal bar, and vertical bar.
    • For every target, measurements for all linear polarization modes (HH, VH, HV, VV) were done consecutively.
    • For every target and every linear polarization mode, measurements were done consecutively for turntable angle = –90° until 90° with 5° increments, where for each of these turntable angles, measurements were done for f = 3 GHz until f = 9 GHz with 0.01 GHz increments.
    • For every target, measurements for all circular polarization modes (RR, LR, RL, LL) were done consecutively. For these cases, the two phase shifters were employed accordingly.
    • For every target and every circular polarization mode, measurements were done consecutively for turntable angle = –90° until 90° with 5° increments, where for each of these turntable angles, measurements were done for f = 3 GHz until f = 9 GHz with 0.01 GHz increments.
The EM wave feeding from port 1 of the VNA to the TX antennas and from the RX antenna to port 2 of the VNA, and the usage of the phase shifters for every polarization mode are described in Table 2 (theoretical) and Table 3 (practical). Et represents the amplitude of the transmitted electric field (at the TX antennas), whereas Er represents the amplitude of the received electric field (at the RX antennas). Here, the usage of Et, Er, and the sine function is only to describe the arrangements of phase shifting in the cases of circular polarization. The values of Et of different measurement cases were not necessarily the same, and the values of Er of different measurement cases were not necessarily the same either.
Please notice that a right-handed circularly polarized (RHCP) TX antenna can be realized by feeding the TX horizontal antenna with a train of sine waves and at the same time feeding the TX vertical antenna with the 90° forward-shifted version of the same wave train. Then, upon backscattering by a target, when expecting at the RX antenna, an echo with the same circular direction (RHCP) will be sensed properly if its vertical part (the field sensed by the RX vertical antenna) is shifted 90° backward before the summation of the horizontal and vertical parts. This is the reason for the phase-shifting direction of the RX vertical antenna being opposite to the TX vertical antenna in the RR case. In the case of LR, the echo received at the RX antenna is expected to have the opposite direction to the transmitted wave, and thus at the RX vertical antenna, the echo is shifted in the same direction as the phase shifting of the TX vertical antenna. This is the main message of Table 2.
Unfortunately, the theoretical phase shifting described in Table 2 requires two types of phase shifters, i.e., +90° and −90°, for realization. When only the same type phase shifters are available for the experiment, then all CP modes can be realized by the technique described in Table 3.
Some activities that took place during the experiment are presented in Figure 8, whereas the physical appearances of the targets are depicted in Figure 9.

3. Results

Results of both simulation and experiment are presented here. Through the following graphs, the backscattered EM wave intensity is depicted for each target and each polarization mode. Every figure presents (Figure 10, Figure 11, Figure 12 and Figure 13) the measurement results of a target where the results of linear polarization and circular polarization can be compared. From the same figure, the simulation results and experimental results of the same target can be compared too. In all figures, the dotted lines are scatter graphs, whereas the solid lines are their second-order polynomial trendlines.

4. Discussion

The measurement was done for 64 cases, i.e., two types of test (simulation and experiment), four different targets, and eight different polarization modes. Furthermore, every case was done for 37 different angles. As shown in the figures, for every measurement case, a second-order polynomial trendline was drawn. Despite the opportunity to address many possible comparisons amongst the cases, the obvious trends are discussed here.

4.1. Relationship between Backscatter Intensity and Aspect Angle

A relationship between backscatter intensity and aspect angle happens for the trihedral cases, as shown in [13] and as indicated by the studies in [10,14]. The study in [13] showed proof that the strongest backscatter intensity happens when the EM wave comes from the direction that the trihedral faces (aspect angle = 0°). The studies in [10,14] investigated the effects of trihedral plate damage to the backscatter intensity, but implicitly also showed that, for an ideal trihedral, the strongest backscatter intensity occurs when the EM wave comes from the direction that the trihedral faces direction. The result of our simulation on the trihedral showed agreement with this phenomenon, where the backscatter intensity of most polarization modes returned to the maximum when the aspect angle was around 0°. The simulation resulted in the same trends for the horizontal bar and the vertical bar, where the maximum backscatter intensities happened when the aspect angle was 0°, and the minimum backscatter intensities occurred when the aspect angle was 90°. The same trend was demonstrated consistently for both linear polarization and circular polarization. A different trend was shown in the case of the sphere, where the backscatter intensity graph was flat. This can be explained by the fact that in the simulation we used a perfect sphere model, and this geometry was exactly the same for all aspect angles.
The experiment showed the same trend for most linear polarization modes of the trihedral, horizontal bar, and vertical bar, although with less curving degrees compared to the simulation results. For the trihedral, the convex curve (maximum at ϕ = 0° and minimum at ϕ = 90°) occurred for the VH, HV, and VV cases. For the horizontal bar, the convex curve occurred for the VH and VV cases. For the vertical bar, a less obvious convex curve occurred for the VH and VV cases. The sphere resulted in flat backscatter intensity lines for all linear polarization cases, which was the same as the simulation results.
The results of the experiments for the circular polarization modes unexpectedly differed. All of these cases resulted in a flat curve, meaning that the backscatter intensity value was not related to the aspect angle. The reason for this anomaly is yet unknown and should be subject to further investigation.

4.2. Standard Deviations of CPs and LPs

It is in our interest to understand the behaviors of circularly polarized EM waves when hitting a target. Here, their backscatter intensities were observed at an aspect angle of 0°. It was interesting to find that the backscatter intensities of CP modes varied less significantly compared to the backscatter intensities of LP modes. This phenomenon is presented in Table 4.
According to the simulation, the standard deviations of CP modes of the trihedral, the horizontal bar, the vertical bar, and the sphere were 1.13, 0.11, 0.66, and 0.33, respectively, which were significantly smaller than the standard deviations of the LP modes, i.e., 3.00, 1.55, 2.59, and 2.09, respectively. The same trend happened in the experiment, where the standard deviations of CP modes of each target were 1.30, 1.50, 0.91, and 0.87, which were significantly smaller than the standard deviations of the LP modes, i.e., 8.11, 8.88, 8.38, and 2.19. This means that the simulation and the experiment show agreement in returning standard deviations for CP modes that are significantly smaller than the standard deviations of LP modes. This phenomenon, in our understanding, can be explained as follows:
  • The RX sensor of an LP system employs either a horizontally or vertically polarized antenna, so it receives only one component of the backscattered waves (either horizontal or vertical). In this way, the values of backscattered wave intensity for each possible mode, i.e., HH, VH, HV, and VV, are potentially significantly different from each other.
  • The RX sensor of a CP system employs both horizontally and vertically polarized antennas plus a phase shifter, where it receives both horizontal and vertical components of the backscattered waves and sums them up (after shifting the phase of one of the two). Due to this summation mechanism, the resulting values of backscattered wave intensity for each possible mode, i.e., RR, LR, RL, and LL, are less significantly different from each other compared to the LP case. Please be reminded that the backscattered waves off a radar object consist of both horizontal and vertical components (likely with different intensity). The differences of the backscattered wave intensities of RR, LR, RL, and LL modes are merely determined by the phase-shifting manner.

4.3. Backscatter Intensity Differences amongst LP Modes

As discussed in Section 4.2, both the simulation and the experiment showed that the standard deviation of the backscatter intensity of linearly polarized EM waves was relatively large. It means that both the simulation and the experiment could differentiate between the backscatter intensities of each HH, VH, HV, and VV case. The details of the variations shown by LP modes (ϕ = 0°) are presented in Table 5, and the details of the CP modes are also presented as a comparison.
The results of the simulation can be described as the following. For the trihedral, HH returned the largest value (−81.15 dB) whereas VH returned the smallest value (−88.85 dB). For the horizontal bar, HH returned the largest value (−81.08 dB) whereas VH returned the smallest value (−85.38 dB). For the vertical bar, the VV case resulted in the largest value (−79.56 dB) whereas VH resulted in the smallest value (−85.29 dB). For the sphere, HV returned the largest value (−74.95 dB) whereas HH returned the smallest value (−80.54 dB).
The experimental results can be described as the following. For the trihedral, HH returned the largest value (−65.39 dB) whereas HV returned the smallest value (−85.56 dB). For the horizontal bar, HH resulted in the largest value (−65.54 dB) whereas VH resulted in the smallest value (−89.53 dB). For the vertical bar, VV returned the largest value (−65.47 dB) whereas VH returned the smallest value (−88.13 dB). For the sphere, HV returned the largest value (−70.87 dB) whereas VV returned the smallest value (−78.50 dB).
It is interesting to find that the simulation results agreed with the experimental results for the largest backscatter intensity values of every target (shaded in green). HH, HH, VV, and HV were the polarization modes where the trihedral, the horizontal bar, the vertical bar, and the sphere resulted in its largest backscatter intensities, respectively, in both the simulation and the experiment. This conveys important information that (i) the medium-sized sphere tends to reverse EM wave polarization, (ii) the medium-sized trihedral tends to retain the polarization, theoretically because of the double reflection (reverse then reverse again), (iii) the medium-sized horizontal bar tends to strongly backscatter horizontally polarized EM waves while retaining their polarity, and (iv) the medium-sized vertical bar tends to strongly backscatter vertically polarized EM waves while retaining their polarity.

5. Summary

The work reported in this article can be summarized as follows:
The measurement was done for four different medium-sized targets, eight different polarization modes, and 37 different angles, by simulation and experiment, where, in this study, medium-sized was defined as 5λ < target size < 10λ.
The simulation was done using our developed tool, AST_V1, that simulates electromagnetic field propagation using the FDTD method. The experiment was done in our anechoic chamber.
The simulation results mostly agreed with previous studies that the backscatter intensity relates to the target’s aspect angle, where the largest backscatter intensity was obtained when the aspect angle is 0° and the least backscatter intensity was obtained when the aspect angle is 90°.
The simulation agreed with the experiment in the observation of circularly polarized EM waves, where the circular polarization resulted in significantly smaller backscatter intensity standard deviation compared to the linear polarization. This can also mean that the circular polarization provided more stable backscattering.
The simulation to some extent agreed with the experiment in the observation of linearly polarized EM waves, where HH, HH, VV, and HV were the polarization modes where the trihedral, the horizontal bar, the vertical bar, and the sphere resulted in the largest backscatter intensity, respectively. This finding confirms the understanding that:
  • A sphere tends to reverse EM wave polarization,
  • A trihedral tends to retain the EM wave polarization,
  • A horizontal bar tends to strongly backscatter horizontally polarized EM waves while retaining their polarity, and
  • A vertical bar tends to strongly backscatter vertically polarized EM waves while retaining their polarity.
  • However, please note that other polarization modes also exist, with less intensity.
The computation was done with the FDTD cell size of 1.25 × 1.25 × 1.25 cm. The accuracy of the computation results can potentially be improved by reducing the cell size.

Supplementary Materials

Supplementary File 1

Author Contributions

Conceptualization, M.N.; data curation, M.N., C.E.S., P.S., A.H.W., and Y.Y.; formal analysis, M.N.; funding acquisition, J.T.S.S.; investigation, M.N. and J.W.; methodology, M.N., J.T.S.S., P.S., Y.Y., and J.W.; project administration, J.T.S.S.; resources, J.T.S.S., C.E.S., P.S., and A.H.W.; software, M.N.; supervision, J.T.S.S.; validation, P.S.; visualization, M.N.; writing—original draft, M.N.; writing—review & editing, J.T.S.S. and P.S.

Funding

This research has been funded by the Chiba University Strategic Priority Research Promotion Program FY2016-FY2018.

Acknowledgments

The authors would like to convey gratitude to the Mitsubishi Corporation International Student Scholarship as well as the Ministry of Research, Technology, and Higher Education of the Republic of Indonesia (Kemenristekdikti) for their support to the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Herrera, J.F.; Moreno, P. Calculation of reflection losses in a small anechoic chamber. IEEE Lat. Am. Trans. 2015, 13, 1258–1264. [Google Scholar] [CrossRef]
  2. Munteanu, I.; Kakerow, R. Simulation methodology for the assessment of field uniformity in a large anechoic chamber. IEEE Trans. Magn. 2014, 50, 213–216. [Google Scholar] [CrossRef]
  3. Chen, J.; Wen, Y. Simulation of the field uniformity of anechoic chamber. In Proceedings of the 2007 International Symposium on Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications, Hangzhou, China, 16–17 August 2006; pp. 1349–1352. [Google Scholar]
  4. Holloway, C.L.; Mckenna, P.M.; Dalke, R.A.; Perala, R.A.; Devor, C.L. Time-Domain Modeling, Characterization, and Measurements of Anechoic and Semi-Anechoic Electromagnetic Test Chambers. IEEE Trans. Electromagn. Compat. 2002, 44, 102–118. [Google Scholar] [CrossRef]
  5. Bornkessel, C.; Wiesbeck, W. Numerical analysis and optimization of anechoic chambers for EMC testing. IEEE Trans. Electromagn. Compat. 1996, 38, 499–506. [Google Scholar] [CrossRef]
  6. Garn, H.; Miillner, W.; Kremser, R.; Seibersdorf, C.; Probe, F. A Theoretical and Experimental Investigation of Field Homogeneity in Absorber-lined Chambers for Radiated-susceptibility Tests According to IEC 801-3 Second Edition. In Proceedings of the IEEE Symposium on Electromagnetic Compatibility, Chicago, IL, USA, 22–26 August 1994; pp. 138–143. [Google Scholar]
  7. Propper, T.; Daniel, U. Computational Simulation of Electromagnetic Fields in a Large-scale Anechoic Chamber. In Proceedings of the IEEE Symposium on Electromagnetic Compatibility, Chicago, IL, USA, 22–26 August 1994; pp. 308–313. [Google Scholar]
  8. Luebbers, K.K.R.; Steich, D.; Ryan, D. Analysis of Compact Electromagnetic Anechoic Chamber Performance Using Finite Difference Time Domain. In Proceedings of the SOUTHEASTCON ′91, Williamsburg, VA, USA, 7–10 April 1991. [Google Scholar]
  9. Rappaport, C.M.; Guerel, T. Reducing the Computational Domain for FDTD Scattering Simulation using Sawtooth Anechoic Chamber ABC. IEEE Trans. Magn. 1995, 31, 1546–1549. [Google Scholar] [CrossRef]
  10. Hänninen, I.; Pitkonen, M.; Nikoskinen, K.I.; Sarvas, J. Method of moments analysis of the backscattering properties of a corrugated trihedral corner reflector. IEEE Trans. Antennas Propag. 2006, 54, 1167–1173. [Google Scholar] [CrossRef]
  11. Andrade, L.A.; Nohara, E.L.; Peixoto, G.G.; Rezende, M.C.; Martin, I.M. Backscattering analysis of flat plate and dihedral corner reflectors using PO and comparison with RCS measurements in anechoic chamber. In Proceedings of the 2003 SBMO/IEEE MTT-S International Microwave and Optoelectronics Conference—IMOC 2003, Foz do Iguacu, Brazil, 20–23 September 2003; Volume 3, pp. 719–724. [Google Scholar]
  12. Gennarelli, C.; Pelosi, G.; Riccio, G. Physical optics analysis of the field backscattered by a depolarising trihedral corner reflector. IEEE Proc. Microw. Antennas Propag. 1998, 145, 213. [Google Scholar] [CrossRef]
  13. Corona, P.; Ferrara, G.; D’Agostino, F.; Gennarelli, C.; Riccio, G. An improved physical optics model for the evaluation of the field backscattered by triangular trihedral corner reflectors. In Proceedings of the 8th Mediterranean Electrotechnical Conference on Industrial Applications in Power Systems, Computer Science and Telecommunications (MELECON 96), Bari, Italy, 16 May 1996; pp. 534–537. [Google Scholar]
  14. Ferrara, F.; Mattia, G.; Posa, F. Backscattering Study on Non-orthogonal Trihedral Corner Reflectors. IEEE Proc. Microw. Antennas Propag. 1995, 142, 140–143. [Google Scholar] [CrossRef]
  15. Griesser, T.; Balanis, C.A. Backscatter Analysis of Dihedral Corner Reflectors Using Physical Optics and the Physical Theory of Diffraction. IEEE Trans. Antennas Propag. 1987, 35, 1137–1147. [Google Scholar] [CrossRef]
  16. Anderson, W.C. Consequences of Nonorthogonality on the Scattering Properties of Dihedral Reflectors. IEEE Trans. Antennas Propag. 1987, 35, 1154–1159. [Google Scholar] [CrossRef]
  17. Schmidt, R.F. Software Engineering Fundamentals. In Software Engineering Architecture-Driven Software Development; Elsevier Inc.: Morgan Kaufmann, MA, USA, 2013; p. 16. [Google Scholar]
  18. Yee, K.S. Numerical Solution of Initial Boundary Value Problems Involving Maxwell’s Equations in Isotropic Media. IEEE Trans. Antenna Propag. 1966, 14, 302–307. [Google Scholar]
  19. Taflove, A.; Hagness, S.C. Introduction to Maxwell’s Equations and the Yee Algorithm. In Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed.; Artech House: Norwood, MA, USA, 2005; pp. 51–105. [Google Scholar]
  20. Bérenger, J.-P. Perfectly Matched Layer (PML) for Computational Electromagnetics; Morgan & Claypool: San Rafeal, CA, USA, 2007; Volume 2, pp. 1–117. [Google Scholar]
  21. Elsherbeni, A.Z.; Demir, V. The Finite-Difference Time-Domain Method for Electromagnetics with MATLAB Simulations, 2nd ed.; SciTech: Princeton, NJ, USA, 2015. [Google Scholar]
  22. Sullivan, D.M. Electromagnetic Simulation Using the FDTF Method, 1st ed.; IEEE Press: New York, NY, USA, 2000. [Google Scholar]
  23. Santosa, C.E.; Sumantyo, J.T.S.; Chua, M.Y.; Urata, K.; Ito, K.; Gao, S. Subarray Design for C-Band Circularly-Polarized Synthetic Aperture Radar Antenna Onboard Airborne. Prog. Electromagn. Res. 2018, 163, 107–117. [Google Scholar] [CrossRef]
  24. Sitompul, P.P.; Sumantyo, J.T.S.; Kurniawan, F.; Santosa, E.C.; Manik, T.; Hattori, K.; Gao, S.; Liu, J.-Y. A Circularly Polarized Circularly-Slotted-Patch Antenna with Two Asymmetrical Rectangular Truncations for Nanosatellite Antenna. Prog. Electromagn. Res. C 2018, 90, 225–236. [Google Scholar] [CrossRef]
  25. Kurniawan, F.; Sumantyo, J.T.S.; Ito, K.; Kuze, H.; Gao, S. Patch Antenna Using Rectangular Centre Slot and Circular Ground Slot for Circularly Polarized Synthetic Aperture Radar (Cp-Sar) Application. Prog. Electromagn. Res. 2018, 160, 51–61. [Google Scholar] [CrossRef]
Figure 1. The measurement plan. The pyramidal absorbers shall be placed evenly at the inner sides of all six chamber’s walls.
Figure 1. The measurement plan. The pyramidal absorbers shall be placed evenly at the inner sides of all six chamber’s walls.
Computation 07 00055 g001
Figure 2. Definitions of angles used in this work. The elevation angle θ is the angle between the target’s facing direction and the horizontal plane. The aspect angle ϕ is the electromagnetic (EM) wave incident angle in the azimuth point of view (in other words, the angle between the target’s facing direction and the vertical plane). Every measurement case in this study shall be done with θ = 0° and –90° ≤ ϕ ≤ 90°.
Figure 2. Definitions of angles used in this work. The elevation angle θ is the angle between the target’s facing direction and the horizontal plane. The aspect angle ϕ is the electromagnetic (EM) wave incident angle in the azimuth point of view (in other words, the angle between the target’s facing direction and the vertical plane). Every measurement case in this study shall be done with θ = 0° and –90° ≤ ϕ ≤ 90°.
Computation 07 00055 g002
Figure 3. Critical design review of Anechoic Chamber Simulation Tool Version 1 (AST_V1) presented in a flow chart (part 1).
Figure 3. Critical design review of Anechoic Chamber Simulation Tool Version 1 (AST_V1) presented in a flow chart (part 1).
Computation 07 00055 g003
Figure 4. Critical design review of AST_V1 presented in a flow chart (part 2).
Figure 4. Critical design review of AST_V1 presented in a flow chart (part 2).
Computation 07 00055 g004
Figure 5. The appearance of the AST_V1 user interface.
Figure 5. The appearance of the AST_V1 user interface.
Computation 07 00055 g005
Figure 6. The four target types synthesized by the tool. In this figure, the elevation angle is 0° and aspect angle is also 0°.
Figure 6. The four target types synthesized by the tool. In this figure, the elevation angle is 0° and aspect angle is also 0°.
Computation 07 00055 g006
Figure 7. Visualization during a simulation (trihedral corner reflector, 3 GHz, RR, θ = 0°, ϕ = 0°). The 3D volume is visualized by applying a horizontal cutting plane and two vertical cutting planes: (a) Antennas and target visualization before the Finite-Difference Time-Domain (FDTD) calculation started, where each TX and RX antenna set consisted of a horizontal and a vertical dipole antenna. (bd) show the FDTD states respectively at steps 180, 519, and 1154 of 1154.
Figure 7. Visualization during a simulation (trihedral corner reflector, 3 GHz, RR, θ = 0°, ϕ = 0°). The 3D volume is visualized by applying a horizontal cutting plane and two vertical cutting planes: (a) Antennas and target visualization before the Finite-Difference Time-Domain (FDTD) calculation started, where each TX and RX antenna set consisted of a horizontal and a vertical dipole antenna. (bd) show the FDTD states respectively at steps 180, 519, and 1154 of 1154.
Computation 07 00055 g007
Figure 8. Some activities that took place during the experiment: (a) Measuring the target material’s dielectric constant, (b) after the antenna and the phase shifter installation, (c) the phase shifter, and (d) antenna–target horizontal alignment using a laser liner.
Figure 8. Some activities that took place during the experiment: (a) Measuring the target material’s dielectric constant, (b) after the antenna and the phase shifter installation, (c) the phase shifter, and (d) antenna–target horizontal alignment using a laser liner.
Computation 07 00055 g008
Figure 9. Targets available for experiments. The cylinder is not part of this report. During the measurement, the bar can be positioned horizontally or vertically.
Figure 9. Targets available for experiments. The cylinder is not part of this report. During the measurement, the bar can be positioned horizontally or vertically.
Computation 07 00055 g009
Figure 10. Backscatter intensity of a medium-sized trihedral (w/λ = 6.4): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Figure 10. Backscatter intensity of a medium-sized trihedral (w/λ = 6.4): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Computation 07 00055 g010
Figure 11. Backscatter intensity of a medium-sized horizontal bar ((l/λ = 6.3, w/λ = 0.4): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Figure 11. Backscatter intensity of a medium-sized horizontal bar ((l/λ = 6.3, w/λ = 0.4): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Computation 07 00055 g011
Figure 12. Backscatter intensity of a medium-sized vertical bar ((l/λ = 6.3, w/λ = 0.4): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Figure 12. Backscatter intensity of a medium-sized vertical bar ((l/λ = 6.3, w/λ = 0.4): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Computation 07 00055 g012
Figure 13. Backscatter intensity of a medium-sized sphere (d/λ = 6.8): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Figure 13. Backscatter intensity of a medium-sized sphere (d/λ = 6.8): (a) Simulation result of LP, (b) simulation result of CP, (c) experimental result of LP, and (d) experimental result of CP.
Computation 07 00055 g013
Table 1. List of targets and their descriptions.
Table 1. List of targets and their descriptions.
No.TargetDescription
1SphereA conducting thin-plate sphere having a diameter of 26 cm
2TrihedralA conducting thin-plate triangular trihedral having a width of 64 cm
3Horizontal barA conducting thin-plate bar of 4 × 63 cm
4Vertical barThe same as target no. 3 but vertically positioned
Table 2. Theoretical wave feeding, wave sensing and phase shifting for linearly polarized (LP) and circularly polarized (CP) modes.
Table 2. Theoretical wave feeding, wave sensing and phase shifting for linearly polarized (LP) and circularly polarized (CP) modes.
Polarization (RX/TX)Electric Field Sensing at Electric Field Feeding to
Horizontal
RX Antenna
Vertical
RX Antenna
Horizontal
TX Antenna
Vertical
TX Antenna
HHe(t) = Erin(2πft)-e(t)= Etsin(2πft)-
VH-e(t) = Ersin(2πft)e(t)= Etsin(2πft)-
HVe(t) = Ersin(2πft)--e(t)= Etsin(2πft)
VV-e(t) = Ersin(2πft)-e(t)= Etsin(2πft)
RRe(t) = Ersin(2πft)e(t) = Ersin(2πft − 0.5π)e(t)= Etsin(2πft)e(t)= Etsin(2πft + 0.5π)
LRe(t) = Ersin(2πft)e(t) = Ersin(2πft + 0.5π)e(t)= Etsin(2πft)e(t)= Etsin(2πft + 0.5π)
RLe(t) = Ersin(2πft)e(t) = Ersin(2πft – 0.5π)e(t)= Etsin(2πft)e(t)= Etsin(2πft − 0.5π)
LLe(t) = Ersin(2πft)e(t) = Ersin(2πft + 0.5π)e(t)= Etsin(2πft)e(t)= Etsin(2πft − 0.5π)
Table 3. Practical wave feeding, wave sensing, and phase shifting for LP and CP modes. Please notice the different arrangements of phase shifting for the receiver (RX) right-handed transmitter (TX) right-handed (RR), RX left-handed TX right-handed (LR), RX left-handed TX right-handed (LR), or RX left-handed TX left-handed (LL) modes compared to Table 3.
Table 3. Practical wave feeding, wave sensing, and phase shifting for LP and CP modes. Please notice the different arrangements of phase shifting for the receiver (RX) right-handed transmitter (TX) right-handed (RR), RX left-handed TX right-handed (LR), RX left-handed TX right-handed (LR), or RX left-handed TX left-handed (LL) modes compared to Table 3.
Polarization (RX/TX)Electric Field Sensing at Electric Field Feeding to
Horizontal
RX Antenna
Vertical
RX Antenna
Horizontal
TX Antenna
Vertical
TX Antenna
HHe(t) = Ersin(2πft)-e(t) = Etsin(2πft)-
VH-e(t) = Ersin(2πft)e(t) = Etsin(2πft)-
HVe(t) = Ersin(2πft)--e(t) = Etsin(2πft)
VV-e(t) = Ersin(2πft)-e(t) = Etsin(2πft)
RRe(t) = Etsin(2πft + 0.5π)e(t) = Ersin(2πft)e(t) = Etsin(2πft)e(t) = Etsin(2πft + 0.5π)
LRe(t) = Ersin(2πft)e(t) = Etsin(2πft + 0.5π)e(t) = Etsin(2πft)e(t) = Etsin(2πft + 0.5π)
RLe(t) = Etsin(2πft + 0.5π)e(t) = Ersin(2πft)e(t) = Etsin(2πft + 0.5π)e(t) = Etsin(2πft)
LLe(t) = Ersin(2πft)e(t) = Etsin(2πft + 0.5π)e(t) = Etsin(2πft + 0.5π)e(t) = Etsin(2πft)
Table 4. Comparison of the standard deviations of CPs and LPs for every target: (a) According to the simulation, and (b) according to the experiment.
Table 4. Comparison of the standard deviations of CPs and LPs for every target: (a) According to the simulation, and (b) according to the experiment.
(a) Simulation: CP’s vs. LP’s St. Dev. (ϕ = 0°)(b) Experiment: CP’s vs. LP’s St. Dev. (ϕ = 0°)
TargetBackscatter Intensity
Standard Deviation (σ)
TargetBackscatter Intensity
Standard Deviation (σ)
Of All LP ModesOf All CP ModesOf All LP ModesOf All CP Modes
Trihedral3.001.13Trihedral8.111.30
Horizontal Bar1.550.11Horizontal Bar8.881.50
Vertical Bar2.590.66Vertical Bar8.380.91
Sphere2.090.33Sphere2.190.87
Table 5. Backscatter intensity differences amongst LP modes compared to CP modes: (a) Simulation and (b) experiment. Please notice that the backscatter intensities of CP modes varied less significantly compared to LP modes. The ones shaded in green are the largest values of each target for the LP modes.
Table 5. Backscatter intensity differences amongst LP modes compared to CP modes: (a) Simulation and (b) experiment. Please notice that the backscatter intensities of CP modes varied less significantly compared to LP modes. The ones shaded in green are the largest values of each target for the LP modes.
(a) Simulation: Backscatter Intensity Differences amongst All Polarization Modes (ϕ = 0°)
Backscatter of LP Modes [dB]Backscatter of CP Modes [dB]
HHVHHVVVSt. Dev.RRLRRLLLSt. Dev.
Trihedral−81.15−88.85−86.22−82.743.00−83.04−82.95−85.37−85.131.13
Horizontal Bar−81.08−85.38−83.45−83.921.55−85.62−85.70−85.41−85.490.11
Vertical Bar−84.43−85.29−79.89−79.562.59−80.85−81.04−79.55−79.740.66
Sphere−80.54−76.64−74.95−78.502.09−89.56−89.61−88.91−88.940.33
(b) Experiment: (a) Backscatter Intensity Differences amongst All Polarization Modes (ϕ = 0°)
Backscatter of LP Modes [dB]Backscatter of CP Modes [dB]
HHVHHVVVSt. Dev.RRLRRLLLSt. Dev.
Trihedral−65.39−84.75−85.56−76.798.11−66.53−65.71−68.41−64.931.30
Horizontal Bar−65.54−89.53−83.98−80.138.88−65.20−65.16−68.67−65.231.50
Vertical Bar−78.80−88.13−82.77−65.478.38−64.53−66.06−66.37−64.290.91
Sphere−75.44−72.92−70.87−76.482.19−86.55−84.60−84.36−85.590.87

Share and Cite

MDPI and ACS Style

Nasucha, M.; Sri Sumantyo, J.T.; Santosa, C.E.; Sitompul, P.; Wahyudi, A.H.; Yu, Y.; Widodo, J. Computation and Experiment on Linearly and Circularly Polarized Electromagnetic Wave Backscattering by Corner Reflectors in an Anechoic Chamber. Computation 2019, 7, 55. https://0-doi-org.brum.beds.ac.uk/10.3390/computation7040055

AMA Style

Nasucha M, Sri Sumantyo JT, Santosa CE, Sitompul P, Wahyudi AH, Yu Y, Widodo J. Computation and Experiment on Linearly and Circularly Polarized Electromagnetic Wave Backscattering by Corner Reflectors in an Anechoic Chamber. Computation. 2019; 7(4):55. https://0-doi-org.brum.beds.ac.uk/10.3390/computation7040055

Chicago/Turabian Style

Nasucha, Mohammad, Josaphat T. Sri Sumantyo, Cahya E. Santosa, Peberlin Sitompul, Agus H. Wahyudi, Yang Yu, and Joko Widodo. 2019. "Computation and Experiment on Linearly and Circularly Polarized Electromagnetic Wave Backscattering by Corner Reflectors in an Anechoic Chamber" Computation 7, no. 4: 55. https://0-doi-org.brum.beds.ac.uk/10.3390/computation7040055

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