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

Authors: Jia Chungang; Guo Lixin; Yang Pengju
Journal: International Journal of Antennas and Propagation
 (2013)

---

## 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 (l) 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.


---

## Body

## 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 Section2, 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
(1)f(xm,yn)=1LxLy∑mk=-M/2+1M/2∑nk=-N/2+1N/2F(kmk,knk)×exp[i(kmkxm+knkyn)],
where
(2)F(kmk,knk)=2π[LxLyS(kmk,knk)]1/2·{[N(0,1)+iN(0,1)]2,mk≠0,M2,andnk≠0,N2N(0,1),mk≠0,M2,ornk≠0,N2,Lx and Ly are the lengths of the rough surface along the x-axis and y-axis, respectively. The numbers of discrete points are M and N. kmk and knk are the spatial frequency at the corresponding point. N(0,1) is a random number following normal distributions function. S(kmk,knk) is the power spectrum density function, which is Gaussian type in this paper and is given as follows:
(3)S(kmk,knk)=σ2lxly4πexp(-kmk2lx2+knk2ly24),
where σ is rms height, and lx and ly are the correlation lengths along the x- and y-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).Simulation of 2D Gaussian randomly rough surface (the unit of bar is meter).
(a)
σ
=
0
.
03
m, lx=ly=0.1m
(b)
σ
=
0
.
03
m, lx=ly=0.2m
(c)
σ
=
0
.
02
m, lx=ly=0.15m
(d)
σ
=
0
.
04
m, lx=ly=0.15m


### 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]
(4)J→s={2n^×H→inc(r→′,ω),litregion0,otherwise,
where H→inc(r→′,ω) is the incident magnetic field, n^ is the unit vector normal to the surface. The far scattering field in the frequency domain is approximately derived as
(5)E→s(r→,ω)=jkz0exp(-jkr)4πr×∬sr^×[r^×J→s(r→′,ω)]exp(jk→·r→′)ds′.
In (5), z0 is the characteristic impedance in free space. r→ is the distant observing point, and r→′ is the scatterer integration point, where J→s(r→′,ω) is the surface current density. s 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
(6)E→s(r→,t)=η4πrc∬sr^×[r^×∂∂tJ→s(r→′,t-τ1)]ds′,
where τ1=r→·(r→-r→′)/c 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:
(7)J→s(r→′,t)=12π∫J→s(r→′,ω)exp(jωt)dω.
Substituting (4) to (7), J→s could also be written as
(8)J→s(r→′,t)=12π∫2n^×H→inc(r→′,ω)exp(jωt)dω.
By taking IFT of the incident magnetic field H→inc(r→′,ω), one can obtain its formulation in time domain:
(9)H→inc(r→′,t)=12π∫H→inc(r→′,ω)exp(jωt)dω.
Based on (6), (8), and (9), the scattering field in the time domain is derived as
(10)E→s(r→,t)=z02πrc∬sr^×[r^×(n^×∂∂tH→inc(r→′,t-τ1-τ2))]ds′,
where τ2=k^i·r→′/c 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 (τ1+τ2) 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 r→. 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]
(11)RCS=limR→∞4πR2S|E→s|2|E→i|2,
where S is the illuminated area of 2D rough surface.


## 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
(1)f(xm,yn)=1LxLy∑mk=-M/2+1M/2∑nk=-N/2+1N/2F(kmk,knk)×exp[i(kmkxm+knkyn)],
where
(2)F(kmk,knk)=2π[LxLyS(kmk,knk)]1/2·{[N(0,1)+iN(0,1)]2,mk≠0,M2,andnk≠0,N2N(0,1),mk≠0,M2,ornk≠0,N2,Lx and Ly are the lengths of the rough surface along the x-axis and y-axis, respectively. The numbers of discrete points are M and N. kmk and knk are the spatial frequency at the corresponding point. N(0,1) is a random number following normal distributions function. S(kmk,knk) is the power spectrum density function, which is Gaussian type in this paper and is given as follows:
(3)S(kmk,knk)=σ2lxly4πexp(-kmk2lx2+knk2ly24),
where σ is rms height, and lx and ly are the correlation lengths along the x- and y-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).Simulation of 2D Gaussian randomly rough surface (the unit of bar is meter).
(a)
σ
=
0
.
03
m, lx=ly=0.1m
(b)
σ
=
0
.
03
m, lx=ly=0.2m
(c)
σ
=
0
.
02
m, lx=ly=0.15m
(d)
σ
=
0
.
04
m, lx=ly=0.15m


## 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]
(4)J→s={2n^×H→inc(r→′,ω),litregion0,otherwise,
where H→inc(r→′,ω) is the incident magnetic field, n^ is the unit vector normal to the surface. The far scattering field in the frequency domain is approximately derived as
(5)E→s(r→,ω)=jkz0exp(-jkr)4πr×∬sr^×[r^×J→s(r→′,ω)]exp(jk→·r→′)ds′.
In (5), z0 is the characteristic impedance in free space. r→ is the distant observing point, and r→′ is the scatterer integration point, where J→s(r→′,ω) is the surface current density. s 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
(6)E→s(r→,t)=η4πrc∬sr^×[r^×∂∂tJ→s(r→′,t-τ1)]ds′,
where τ1=r→·(r→-r→′)/c 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:
(7)J→s(r→′,t)=12π∫J→s(r→′,ω)exp(jωt)dω.
Substituting (4) to (7), J→s could also be written as
(8)J→s(r→′,t)=12π∫2n^×H→inc(r→′,ω)exp(jωt)dω.
By taking IFT of the incident magnetic field H→inc(r→′,ω), one can obtain its formulation in time domain:
(9)H→inc(r→′,t)=12π∫H→inc(r→′,ω)exp(jωt)dω.
Based on (6), (8), and (9), the scattering field in the time domain is derived as
(10)E→s(r→,t)=z02πrc∬sr^×[r^×(n^×∂∂tH→inc(r→′,t-τ1-τ2))]ds′,
where τ2=k^i·r→′/c 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 (τ1+τ2) 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 r→. 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]
(11)RCS=limR→∞4πR2S|E→s|2|E→i|2,
where S 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 Figure2, where an incident wave impinges on the surface in the direction of k^i, which makes angle θi relevant to the z-axis and ϕi relative to the x-axis. The scattered angle is θs, and the scattered azimuthal is ϕs. The polarization angle is defined as η.Figure 2
Definition of scattering geometry with relevant parameters in the coordinatexyz.Because TDPO is valid only within the high frequency approximation, the incident source is a modulated Gaussian pulse, which is written as follows:(12)ei(t)=-cos(ωt)exp[-4π(t-t0)τ2],
where ω=2πf0 is the modulation circular frequency, f0 is the modulation frequency, t0 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 isS=2.4m×2.4m, the rms height is σ=0.01m, and the correlation length is lx=ly=0.1m. The frequency band width and the pulse width of the modulated Gaussian pulse are f = 1~4 GHz and τ=1.875ns, respectively. In Figure 3(a), the incident angle is θi=30∘, the azimuth angle is ϕi=0∘, and the polarization angle is η=0∘. In Figure 3(b), the incident parameters are set to θi=45∘, ϕi=0∘, and η=90∘.TDPO and KA methods for the wide-band RCS from 2D PEC rough surface.
(a)
θθ polarization in backward direction (θi=30∘)
(b)
φ
φ polarization in specular direction (θi=45∘)Figure3(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 Figure4, the scattering results from 2D PEC rough surface in the specular direction are shown for φφ polarization discussed with different sizes of surface (S=1.2m×1.2m, S=2.4m×2.4m), where the rms height and the correlation length are σ=0.01m and lx=ly=0.1m, respectively. The incident angle is θi=30∘, the incident azimuth angle is ϕi=0∘, and the incident frequency is f = 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 S for the whole frequency bands. In addition, it is also seen that the variation tendency of the curve with parameter S=1.2m×1.2m is similar to that with S=2.4m×2.4m. 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.Scattering from 2D rough surface in specular direction for different sizes of rough surface.
(a)
Transient response (S=1.2m×1.2m)
(b)
Transient response (S=2.4m×2.4m)
(c)
Wide-band RCSFigures5(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 S=2.4m×2.4m, σ=0.01m, and for the case of lx=ly=0.1m, lx=ly=0.2m and lx=ly=0.3m, respectively. The incident angle is θi=30∘, and the incident azimuth angle is ϕ=0∘. It is found that the first pulse signal has little difference, but it is obviously observed that the second pulse increases with increasing l.Scattering from the 2D PEC rough surface in the specular direction for different correlation lengths.
(a)
Transient response (lx=ly=0.1m)
(b)
Transient response (lx=ly=0.2m)
(c)
Transient response (lx=ly=0.3m)
(d)
Wide-band RCSFigure5(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 l, 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 (lx=ly=0.2m, σ=0.006m, 0.008m, and 0.01m) are presented in Figure 6, where the size of rough surface is S=2.4m×2.4m, the incident angle is θ=30∘, and the azimuth angle is ϕi=0∘. 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.Scattering from the 2D PEC rough surface in the specular direction for different root square mean heights.
(a)
Transient response (σ=0.006m)
(b)
Transient response (σ=0.008m)
(c)
Transient response (σ=0.01m)
(d)
Wide-band RCSTo 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 S=2.4m×2.4m and the rms height and the correlation length are σ=0.01m and lx=ly=0.01m, respectively. The frequency band width of the Gaussian pulse is f = 1~4 GHz.Backscattering from the 2D PEC rough surface for different incident angles.
(a)
Transient response (θi=30∘)
(b)
Transient response (θi=45∘)
(c)
Wide-band RCS in the backward directionFigures7(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.


## 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 Figure2, where an incident wave impinges on the surface in the direction of k^i, which makes angle θi relevant to the z-axis and ϕi relative to the x-axis. The scattered angle is θs, and the scattered azimuthal is ϕs. The polarization angle is defined as η.Figure 2
Definition of scattering geometry with relevant parameters in the coordinatexyz.Because TDPO is valid only within the high frequency approximation, the incident source is a modulated Gaussian pulse, which is written as follows:(12)ei(t)=-cos(ωt)exp[-4π(t-t0)τ2],
where ω=2πf0 is the modulation circular frequency, f0 is the modulation frequency, t0 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 isS=2.4m×2.4m, the rms height is σ=0.01m, and the correlation length is lx=ly=0.1m. The frequency band width and the pulse width of the modulated Gaussian pulse are f = 1~4 GHz and τ=1.875ns, respectively. In Figure 3(a), the incident angle is θi=30∘, the azimuth angle is ϕi=0∘, and the polarization angle is η=0∘. In Figure 3(b), the incident parameters are set to θi=45∘, ϕi=0∘, and η=90∘.TDPO and KA methods for the wide-band RCS from 2D PEC rough surface.
(a)
θθ polarization in backward direction (θi=30∘)
(b)
φ
φ polarization in specular direction (θi=45∘)Figure3(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 Figure4, the scattering results from 2D PEC rough surface in the specular direction are shown for φφ polarization discussed with different sizes of surface (S=1.2m×1.2m, S=2.4m×2.4m), where the rms height and the correlation length are σ=0.01m and lx=ly=0.1m, respectively. The incident angle is θi=30∘, the incident azimuth angle is ϕi=0∘, and the incident frequency is f = 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 S for the whole frequency bands. In addition, it is also seen that the variation tendency of the curve with parameter S=1.2m×1.2m is similar to that with S=2.4m×2.4m. 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.Scattering from 2D rough surface in specular direction for different sizes of rough surface.
(a)
Transient response (S=1.2m×1.2m)
(b)
Transient response (S=2.4m×2.4m)
(c)
Wide-band RCSFigures5(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 S=2.4m×2.4m, σ=0.01m, and for the case of lx=ly=0.1m, lx=ly=0.2m and lx=ly=0.3m, respectively. The incident angle is θi=30∘, and the incident azimuth angle is ϕ=0∘. It is found that the first pulse signal has little difference, but it is obviously observed that the second pulse increases with increasing l.Scattering from the 2D PEC rough surface in the specular direction for different correlation lengths.
(a)
Transient response (lx=ly=0.1m)
(b)
Transient response (lx=ly=0.2m)
(c)
Transient response (lx=ly=0.3m)
(d)
Wide-band RCSFigure5(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 l, 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 (lx=ly=0.2m, σ=0.006m, 0.008m, and 0.01m) are presented in Figure 6, where the size of rough surface is S=2.4m×2.4m, the incident angle is θ=30∘, and the azimuth angle is ϕi=0∘. 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.Scattering from the 2D PEC rough surface in the specular direction for different root square mean heights.
(a)
Transient response (σ=0.006m)
(b)
Transient response (σ=0.008m)
(c)
Transient response (σ=0.01m)
(d)
Wide-band RCSTo 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 S=2.4m×2.4m and the rms height and the correlation length are σ=0.01m and lx=ly=0.01m, respectively. The frequency band width of the Gaussian pulse is f = 1~4 GHz.Backscattering from the 2D PEC rough surface for different incident angles.
(a)
Transient response (θi=30∘)
(b)
Transient response (θi=45∘)
(c)
Wide-band RCS in the backward directionFigures7(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.



---

*Source: 584260-2013-09-19.xml*