International Journal of Antennas and Propagation

Volume 2013 (2013), Article ID 584260, 9 pages

http://dx.doi.org/10.1155/2013/584260

## Time-Domain Physical Optics Method for the Analysis of Wide-Band EM Scattering from Two-Dimensional Conducting Rough Surface

School of Science, Xidian University, Xi’an 710071, China

Received 2 March 2013; Accepted 18 August 2013

Academic Editor: Daniel S. Weile

Copyright © 2013 Jia Chungang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Time-domain physical optics (TDPO) method is extended to investigate electromagnetic (EM) scattering from two-dimensional (2D) perfectly electrically conducting (PEC) rough surface in both time domain and frequency domain. The scheme requires relatively small amounts of computer memory and CPU time, and has advantage over the Kirchhoff Approximation (KA) method in obtaining transient response of rough surface by a program run. The 2D Gaussian randomly rough surface is generated by Monte Carlo method and then is partitioned into small triangle facets through the meshing preprocess. The accuracy of TDPO is validated by comparing the numerical results with those obtained by the KA method in both backward and specular directions. The transient response and its frequency distribution of radar cross section (RCS) from rough surface is shown, respectively. The scattering results from rough surface with different size in the specular direction are given. The influence of the root mean square height (*σ*) and correlation length () on electromagnetic scattering from PEC rough surface is discussed in detail. Finally, the comparisons of backscattering results at different incident angles are presented and analyzed.

#### 1. Introduction

Nowadays, the problems of EM scattering from randomly rough surface have been widely applied in fields of remote sensing, target identification, radar detection, and so on. A lot of investigations are carried out to deal with this random problem in terms of roughness and slope of the rough surface. Take the approximate analytical methods, for example, the KA method [1], which is considerably valid approach to solve the problem that rough surface is relatively smooth, the small-perturbation method (SPM) in which standard derivation of rough surface is small enough compared with the wavelength [2], the small-slope approximation (SSA) that is applied to small and large surface roughness [3], and the two-scale method (TSM) [4], which considers that the rough surface can be separated into large- and small-scale surfaces, allowing use of physical optics (PO) and SPM methods. These techniques are limited by the roughness, the incident angle, and so forth. Some numerical methods are also employed, which include the Monte Carlo method [5], the finite difference time-domain (FDTD) [6, 7] method, method of moments (MOM) [8, 9], the finite element method (FEM) [10], and the forward back method (FBM) [11], the fast multipole method (FMM) [12].

However, these approaches and accelerated algorithms can hardly handle the model that is considerably electrically large due to complex computation. High frequency technologies are still extremely valid to solve the problem that the electrical size of rough surface is large, but high frequency methods in time domain rarely appeared in the previous literatures about investigating on scattering from the randomly rough surface, and transient scattering has important significance in synthetic aperture radar (SAR) imaging or inverse SAR (ISAR) imaging. Up to now, to our knowledge, few works have been reported for the wide-band scattering of rough surface by TDPO. In this paper, emphasis is put on discussing this situation in the time domain. The TDPO method is the time version of physical optics, which reserves the advantage of fast computing speed and less memory demand, and can obtain wide-band property by a program running. This paper is devoted to the wide-band scattering characteristics of the 2D rough surface by utilizing TDPO method, which can deal with the electrically large rough surface problem in the time domain. The earlier TDPO method is not carried out in the time domain directly. Each frequency component response is obtained by the frequency method, and the time-domain response of the origin signal is then achieved by inverse Fourier transform (IFT). TDPO is firstly presented by En-Yuan Sun, and the method is used to analyze the scattering field of paraboloidal reflector and hyperboloidal reflector [13]. In [14], the scheme is employed to deal with the large-size structure of the combinative objects.

In this paper, we employ TDPO to investigate scattering from the 2D PEC randomly rough surface. The method is performed to obtain the transient response reflected from rough surface, and then the property of wide-band frequency is obtained by fast Fourier transform (FFT). The formulation of TDPO is presented in [13], which derives the far field by transforming the equivalent electric current density and far-field expression in the frequency domain to those in the time domain. The model of 2D rough surface is presented by the Monte Carlo method. The comparison between the wide-band RCS obtained by TDPO and calculated by the sweep-frequency KA method verifies the validation of TDPO. Finally, the presented scheme is utilized to analyze scattering properties of 2D randomly rough surface.

The paper remaining is organized as follows. In Section 2, the formulation of TDPO is presented, as well as the theoretical formulae of 2D randomly rough surface which is of Gaussian type are given. Several examples of 2D rough surface generated with different parameters of root mean square height (rms height) and correlation length are demonstrated. Section 3 provides the validity of TDPO by comparing the numerical results with those obtained by sweep-frequency KA. The influence of size of rough surface, correlation length, rms height, and the incident angle on the scattering from rough surface is discussed. Finally, some concluding remarks are addressed in Section 4.

#### 2. Theory and Formulation

##### 2.1. Randomly Rough Surface Generation

Firstly, 2D randomly rough surface is modeled in order to determine its scattering characteristic. Using the spectrum density, rough surface is simulated by Monte Carlo method [15] by which the power spectrum is filtered in frequency domain, and then the height of rough surface can be obtained by inverse fast Fourier transform (IFFT). The height distribution function of rough surface is written as where and are the lengths of the rough surface along the -axis and -axis, respectively. The numbers of discrete points are and . and are the spatial frequency at the corresponding point. is a random number following normal distributions function. is the power spectrum density function, which is Gaussian type in this paper and is given as follows: where is rms height, and and are the correlation lengths along the - and -direction, by which the profile of the rough surface is determined. Examples of 2D randomly rough surfaces are simulated and depicted with different correlation length and rms height in Figures 1(a)–1(d).

##### 2.2. Time-Domain Physical Optics Method for 2D PEC Rough Surface

Based on the generation of the rough surface profiles, the surface is divided into small triangle facets though mesh processing. For the conducting rough surface, the surface-current-density distribution by the physical optics method in the frequency domain approximates as [16] where is the incident magnetic field, is the unit vector normal to the surface. The far scattering field in the frequency domain is approximately derived as In (5), is the characteristic impedance in free space. is the distant observing point, and is the scatterer integration point, where is the surface current density. is the surface area of the lighten region. In order to derive the scattering field in the time domain, the IFT is needed for (5), and the electric field in time domain can be obtained as follows where is time delay from the integration point to the distant observing point. The relation of surface-current density between the time domain and frequency domain satisfies the Fourier transform: Substituting (4) to (7), could also be written as By taking IFT of the incident magnetic field , one can obtain its formulation in time domain: Based on (6), (8), and (9), the scattering field in the time domain is derived as where is time delay in the incident direction. For the derivation in detail, one can refer to [12]. From the equation above, integration over the whole scattering area only refers to the incident field with time delay () in the lit region and is not related to interaction between the other points on the surface. For each small triangle facet in lit region, (10) is implemented to calculate scattered electric fields in time domain at the point . The wide-band RCS is obtained through carrying out FFT on transient response. And the normalized RCS in the far zone is defined as [17] where is the illuminated area of 2D rough surface.

#### 3. Numerical Results and Discussions

##### 3.1. Validation of TDPO

In this section, the TDPO method is utilized to investigate the problem of transient scattering from 2D Gaussian randomly rough surface. It is noted that rough surface with finite length is concerned, similar to that in [18]; no tapering or window is introduced. Here, Gaussian pulse is chosen as the incident source. The proposed method is employed to calculate wide-band scattering from rough surface in the specular and backward direction with different incident angles, respectively, and its accuracy against KA is verified by comparing the and polarized results. In this paper, the proposed TDPO method is exhibited by averaging 15 Monte Carlo realizations.

Geometry of the scattering problem is illustrated and defined in Figure 2, where an incident wave impinges on the surface in the direction of , which makes angle relevant to the -axis and relative to the *-*axis. The scattered angle is , and the scattered azimuthal is . The polarization angle is defined as .

Because TDPO is valid only within the high frequency approximation, the incident source is a modulated Gaussian pulse, which is written as follows: where is the modulation circular frequency, is the modulation frequency, is the time delay, and is the pulse width.

In order to ensure the validity of the algorithm presented in this paper, we firstly calculate the RCS of 2D PEC rough surface using the TDPO method and KA, respectively. The size of rough surface generated is , the rms height is , and the correlation length is . The frequency band width and the pulse width of the modulated Gaussian pulse are = 1~4 GHz and , respectively. In Figure 3(a), the incident angle is , the azimuth angle is , and the polarization angle is . In Figure 3(b), the incident parameters are set to , , and .

Figure 3(a) compares the analysis results obtained by TDPO and those by KA in the backward directions with polarization. Figure 3(b) presents the results of polarization in the specular direction using above two mentioned methods. It is obvious that the scattering results by TDPO are coincident with those of KA method, which verifies the validation of the proposed TDPO. In addition, the curve of RCS from the rough surface in the specular direction is much smoother than that in the backward direction.

##### 3.2. Discussion on the Results of EM Scattering from Rough Surface

In this section, the proposed TDPO method is employed to analyze rough surface with different sizes and different scale of rough surface in terms of correlation length and rms height both in the time and frequency domain. Moreover, the backscattering results for different incident angles are discussed.

In Figure 4, the scattering results from 2D PEC rough surface in the specular direction are shown for polarization discussed with different sizes of surface (, ), where the rms height and the correlation length are and , respectively. The incident angle is , the incident azimuth angle is , and the incident frequency is = 1~4 GHz. Figures 4(a)-4(b) investigate the transient response rough surface, where it is obvious that the magnitude of the electric field gets larger when the size of the rough surface increases. We attribute this phenomenon to the fact that the scattering in the specular direction becomes stronger when the size of rough surface is larger. Figure 4(c) plots the wide-band RCS from the rough surface, where one can find that the RCS increases with increasing the surface size for the whole frequency bands. In addition, it is also seen that the variation tendency of the curve with parameter is similar to that with . The reason for this phenomenon is that the parameters of the two kinds of rough surface are the same in terms of the correlation length and the rms height, which determine the profile of surface.

Figures 5(a)–5(c) plot the polarized electric fields in the time domain from the rough surface in the specular direction, where the parameters of rough surface are , , and for the case of , and , respectively. The incident angle is , and the incident azimuth angle is . It is found that the first pulse signal has little difference, but it is obviously observed that the second pulse increases with increasing .

Figure 5(d) shows the analysis result of wide-band characteristics through transforming transient response into frequency. It is obvious that RCS in the specular direction increases with increasing correlation length , especially at the high frequency band. The reason for this is that by keeping the rms height constant and by increasing the correlation length, the electromagnetic roughness is constant, but the rms slope decreases, which leads to an increase of the scattered energy in the coherent scattering direction.

The results of transient response and wide-band RCS from 2D rough surface in specular direction for different root square mean heights (, , , and ) are presented in Figure 6, where the size of rough surface is , the incident angle is , and the azimuth angle is . The polarization is discussed. Figures 6(a)–6(c) show the transient response from rough surface, where the magnitude of the second signal decreases with increasing rms height . In Figure 6(d), wide-band RCS decreases obviously with larger over the whole frequency range; the primary reason for this is because the roughness increases with the increase of , which results in the decrease of the scattering in the specular direction with increasing of the surface roughness.

To further explore the important scattering characteristics of rough surface, the backscattering results ( polarization) for different incident angles are compared in Figure 7, where the size of Gaussian rough surface is and the rms height and the correlation length are and , respectively. The frequency band width of the Gaussian pulse is = 1~4 GHz.

Figures 7(a)-7(b) illustrate the backscattered electric fields in time domain, where the magnitude of electric field is visually smaller with increasing incident angles. The wide-band RCS of 2D PEC rough surface is depicted in Figure 7(c). One can observe that the RCS becomes smaller with the increase of the incident angles over almost all the frequency band of 1–4 GHz, which is caused by the fact that the smaller angle leads to stronger scattering in the backward direction.

#### 4. Conclusion

In this paper, a time-domain high frequency TDPO method is presented and adopted to investigate the scattering characteristics of 2D PEC Gaussian rough surface. Firstly, the wide-band RCS from 2D rough surface in both backward and specular directions calculated by the TDPO method are compared with the results obtained by the KA method. Good agreement is achieved for the two cases, which verifies the validation of the presented TDPO method. Then, the scattering characteristics of rough surface with different size are analyzed. The transient response in time domain and the wide-band RCS obtained by FFT in frequency domain are presented to examine the effect of the correlation lengths and rms height on scattering properties in detail. Finally, we will put our focus on investigating the dielectric and lossy surfaces by employing MECA [19] in the time domain (TDMECA). Furthermore, the backscattering from rough surface is present and analyzed for the different incident angles. Similar to KA method [20], the TDPO method presented in this paper is invalid for the case of larger incident angle or the larger roughness due to the neglect of the multiple scattering. It should be pointed out that the future investigation on this topic will include the scattering from 2D conducting rough surface at low grazing angle illumination with larger roughness and scattering from the dielectric rough surfaces.

#### Acknowledgments

This work was supported by the National Science Foundation for Distinguished Young Scholars of China (Grant no. 61225002), the Specialized Research Fund for the Doctoral Program of Higher Education (Grant no. 20100203110016), and the Fundamental Research Funds for the Central Universities (Grant no. K50510070001).

#### References

- D. Holliday, “Resolution of a controversy surrounding the Kirchhoff approach and the small perturbation method in rough surface scattering theory,”
*IEEE Transactions on Antennas and Propagation*, vol. 35, no. 1, pp. 120–122, 1987. View at Google Scholar · View at Scopus - S. O. Rice, “Reflection of electromagnetic waves from slightly rough surfaces,” in
*Theory of Electromagnetic Waves*, M. Kline, Ed., pp. 351–378, Wiley, New York, NY, USA, 1951. View at Google Scholar - G. Berginc, “Small-slope approximation method: a further study of vector wave scattering from two-dimensional surfaces and comparison with experimental data,”
*Progress in Electromagnetics Research*, vol. 37, pp. 251–287, 2002. View at Google Scholar - S. L. Durden and J. F. Vesecky, “Numerical study of the separation wavenumber in the two-scale scattering approximation,”
*IEEE Transactions on Geoscience and Remote Sensing*, vol. 28, no. 2, pp. 271–272, 1990. View at Publisher · View at Google Scholar · View at Scopus - R. M. Axline and A. K. Fung, “Numerical computations of scattering from a perfectly conducting random surface,”
*IEEE Transactions on Antennas and Propagation*, vol. 26, no. 3, pp. 482–488, 1978. View at Google Scholar · View at Scopus - C. H. Chan, S. H. Lou, L. Tsang, and J. A. Kong, “Electromagnetic scattering of waves by random rough surface: a finite-difference time-domain approach,”
*Microwave and Optical Technology Letters*, vol. 4, no. 9, pp. 355–359, 1991. View at Google Scholar · View at Scopus - J. Li, L.-X. Guo, and H. Zeng, “FDTD investigation on bistatic scattering from two-dimensional rough surface with UPML absorbing condition,”
*Waves in Random and Complex Media*, vol. 19, no. 3, pp. 418–429, 2009. View at Publisher · View at Google Scholar · View at Scopus - R. R. Lentz, “A numerical study of electromagnetic scattering from ocean-like surfaces,”
*Radio Science*, vol. 9, no. 12, pp. 1139–1146, 1974. View at Google Scholar · View at Scopus - R. T. Marchand, “On the use of finite surfaces in the numerical prediction of rough surface scattering,”
*IEEE Transactions on Antennas and Propagation*, vol. 47, no. 4, pp. 600–604, 1999. View at Publisher · View at Google Scholar · View at Scopus - S. H. Lou, L. Tsang, and C. H. Chan, “Application of the finite element method to Monte Carlo simulations of scattering of waves by random rough surfaces: penetrable case,”
*Waves in Random Media*, vol. 1, no. 4, article 006, pp. 287–307, 1991. View at Publisher · View at Google Scholar · View at Scopus - A. Iodice, “Forward-backward method for scattering from dielectric rough surfaces,”
*IEEE Transactions on Antennas and Propagation*, vol. 50, no. 7, pp. 901–911, 2002. View at Publisher · View at Google Scholar · View at Scopus - V. Jandhyala, E. Michielssen, S. Balasubramaniam, and W. C. Chew, “A combined steepest descent-fast multipole algorithm for the fast analysis of three-dimensional scattering by rough surfaces,”
*IEEE Transactions on Geoscience and Remote Sensing*, vol. 36, no. 3, pp. 738–748, 1998. View at Publisher · View at Google Scholar · View at Scopus - E.-Y. Sun and W. V. T. Rusch, “Time-domain physical-optics,”
*IEEE Transactions on Antennas and Propagation*, vol. 42, no. 1, pp. 9–15, 1994. View at Publisher · View at Google Scholar · View at Scopus - L.-X. Yang, D.-B. Ge, and B. Wei, “FDTD/TDPO hybrid approach for analysis of the EM scattering of combinative objects,”
*Progress in Electromagnetics Research*, vol. 76, pp. 275–284, 2007. View at Publisher · View at Google Scholar · View at Scopus - L. Tsang and J. A. Kong,
*Scattering of Electromagnetic Waves- Numerical Simulations*, Wiley, New York, NY, USA, 2000. - C. Scott,
*Modern Methods of Reflector Antenna Analysis and Design*, Artech H, Boston, Mass, USA, 1990. - J. Li, L.-X. Guo, and H. Zeng, “FDTD method investigation on the polarimetric scattering from 2-D rough surface,”
*Progress in Electromagnetics Research*, vol. 101, pp. 173–188, 2010. View at Google Scholar · View at Scopus - J. Li, B. Wei, Q. He, L. Guo, and D. Ge, “Time-domain iterative physical optics method for analysis of EM scattering from the target half buried in rough surface: PEC case,”
*Progress in Electromagnetics Research*, vol. 121, pp. 391–408, 2011. View at Google Scholar · View at Scopus - J. G. Meana, J. Á. Martínez-Lorenzo, F. Las-Heras, and C. Rappaport, “Wave scattering by dielectric and lossy materials using the Modified Equivalent Current Approximation (MECA),”
*IEEE Transactions on Antennas and Propagation*, vol. 58, no. 11, pp. 3757–3761, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Collaro, G. Franceschetti, M. Migliaccio, and D. Riccio, “Gaussian rough surfaces and Kirchhoff approximation,”
*IEEE Transactions on Antennas and Propagation*, vol. 47, no. 2, pp. 392–398, 1999. View at Publisher · View at Google Scholar · View at Scopus