This paper proposes a method to estimate directly the incoherent scattered intensity and radar cross section (RCS) from the effective permittivity of a random media. The proposed method is derived from the original concept of incoherent scattering. The incoherent scattered field is expressed as a simple formula. Therefore, to reduce computation time, the proposed method can estimate the incoherent scattered intensity and RCS of a random media. To verify the potential of the proposed method for the desired applications, we conducted a Monte-Carlo analysis using the method of moments; we characterized the accuracy of the proposed method using the normalized mean square error (NMSE). In addition, several medium parameters, such as the density of scatterers and analysis volume, were studied to understand their effect on the scattering characteristics of a random media. The results of the Monte-Carlo analysis show good agreement with those of the proposed method, and the NMSE values of the proposed method and Monte-Carlo analysis are relatively small at less than 0.05.
In designing radar and remote sensing systems to detect objects buried in the ground or submerged within a chaff cloud, it is important to accurately grasp the characteristics of wave propagation in a random media. In the case of a chaff cloud, the atmosphere and chaff fibers can be modelled as the host material and inclusions of a random medium, respectively.
In particular, to enhance the detection capability of a radar system, it is essential that the clutter signature from the random media surrounding a target be removed. This clutter signal by the random media should be predicted and reflected in the radar system during the design phase
. Even in cases of obtaining synthetic aperture radar images, the image quality is highly dependent on the accurate estimation of the radar clutter signal by the random media
The methods used to estimate the scattered field from a random medium can be classified as Monte-Carlo analysis and effective permittivity calculation. The former method uses multiple estimations of the wave scattering for a large number of inclusions through the use of low-frequency methods such as the method of moments (MoM) and the T-matrix
. The computed results are usually averaged over the realizations; therefore, this method provides an accurate solution through the use of numerical methods, despite suffering from an enormous amount of computation time and computer resources. The latter method, on the other hand, is the most commonly used one, and includes a mixing formula such as the Maxwell-Garnet, Polder-van Sante, or quasi-crystalline approximation with coherent potential
. The mixing formula is based on the concept of low frequency; thus, the analysis frequency is limited. In particular, this type of formula can be applied to calculate the effective permittivity of a random medium with inclusions, which have a formal shape, a uniformly random position, and orientations in all directions. The generalized equivalent conductor (GEC) method has recently been announced; this method can be used to determine the effective permittivity of a randomly oriented media with a probability distribution function
. The obtained effective permittivity of a random medium defines the characteristics of the medium. In tracing the history of an effective permittivity calculation, in most cases, the effective permittivity is only used to obtain the extinction coefficient (or absorption attenuation) and to estimate the attenuation of the field intensity
The scattered field generated by a random medium with random properties consists of a coherent field and an incoherent field. An incoherent field represents fluctuations in the scattered field, and is therefore usually considered trivial or negligible. In general radar systems, the signatures reflected by a target are expressed as the radar cross section (RCS). The average RCS of a random medium is the sum of the coherent and incoherent RCSs due to the coherent and incoherent fields, respectively
. From the point of the average RCS, consequently, an incoherent field is also a significant factor that cannot be ignored.
To evaluate an incoherent scattered field, in a conventional random media field, the effective propagation constant and amplitude of the coherent plane wave are gained by solving an integral equation; these values are used to construct a coherent transmitted field, which is used to obtain a coherent scattered field. The incoherent scattered intensity (squared amplitude) is then obtained through a distorted Born approximation
. Because the series of procedures for the determination of the incoherent scattered field uses the spherical vector function of the T-matrix approach, obtaining unknown coefficients is cumbersome, and the incoherent scattered intensity can only be obtained by using a coherent transmitted field.
The method for estimating the coherent field from the effective permittivity can be used to provide a clear prediction, which has been verified by other analysis methods
. Few studies on the determination of an incoherent field using the effective permittivity can be found in the literature. In
, calculations of the incoherent scattered intensity from the effective permittivity were conducted. However, the calculations were not derived from electromagnetic theory, and the derivation principles were unclear; moreover, some errors in the calculation results can be observed.
, the calculation of the effective permittivity from both an incoherent and a coherent scattered field is proposed. Based on this method, the present paper proposes a direct method for estimating the incoherent scattered intensity and RCS from the effective permittivity, which is obtained using the GEC method. In comparison with the results of
, the results of the proposed method are derived from the original concept of incoherent scattering. Additionally, since the proposed method uses the effective permittivity estimated with the GEC method, it can estimate the incoherent RCS of scatterers with arbitrary orientation distribution. The potential of the proposed method is verified through comparison with the results from a Monte-Carlo analysis using the MoM. To investigate the characteristics of coherent and incoherent scattering, we also estimate the RCS of a random medium for various inclusion densities.
The remainder of this paper is organized as follows. In Section II, the coherent, incoherent, and average RCSs of a random medium are briefly explained, and the proposed approach for directly calculating an incoherent RCS is presented. The results of the RCS with respect to the depths of a random medium are shown in Section III, using the orientation distribution and the inclusion density. Finally, Section IV summarizes this work and provides some concluding remarks.
- 1. Review of Scattered Field and Average RCS in Random Media
Consider multiple identical scatterers distributed randomly and oriented arbitrarily within volume
of a free space, as shown in
is regarded as a random medium owing to the random characteristics of the multiple scatterers. The magnitude and phase of the scattered field of the random medium fluctuate randomly in space and time; hence, the scattered field is a function of position
, and the number of realizations
. For a total of
realizations, the scattered field by the
th realization is represented by the sum of the coherent and incoherent scattered fields, which can be given as follows:
(1) E s (m) = E coh + E incoh (m) ,
where the coherent scattered field,
, is the average scattered field for each realization. That is, the coherent scattered field is irrelevant to any particular realization. On the other hand, the incoherent scattered field,
E incoh (m)
, indicates fluctuations, and its sum is zero. Therefore, the coherent and incoherent scattered fields can be expressed as
(2) E coh = 1 M ∑ m=1 M E s (m) ,
(3) ∑ m=1 M E incoh (m) =0.
Equalization model of random medium: (a) realizations of random medium with multiple identical scatterers and (b) homogeneous media with volume current and effective permittivity.
If the scatterers are totally random, then the scattered field will have nearly incoherent components, and the coherent scattered field will be diminished. The more the scatterers tend toward having an arbitrary orientation or position, the more the coherent scattered field increases. Based on the same principle, in the line-of-sight problem, the coherent scattered field is almost the only field in the proximity of the transmitter. However, at a farther distance from the transmitter, there is a strengthening of the incoherent scattered field
The coherent scattered field depends on the size and shape of the entire analysis volume containing multiple scatterers, and can be calculated by assuming the entire analysis volume as a homogeneous material with an effective permittivity. On the other hand, an incoherent scattered field is caused by fluctuations in the dielectric constant between the inclusions/spaces and the background of the random medium. For a different arrangement of the scatterers for each realization, a fluctuating scattered field arises from spatial variations in the dielectric constant between the inclusions/spaces and the background
When a unit plane wave is incident on the random medium, the instantaneous RCS of the
th realization becomes
(4) σ (m) =4π E s (m) E s (m) ∗ =4π E coh E coh ∗ +4π[ E coh E coh (m)∗ + E incoh (m) E coh ∗ + | E incoh (m) | 2 ],
where the first term on the right-hand side arises from the coherent scattered field, and the second term arises from both the coherent and incoherent scattered fields. For a total of
realizations, the average RCS is therefore given by
(5) 〈 σ 〉 =4π〈 | E s (m) | 2 〉=4π | E coh | 2 +4π〈 | E incoh (m) | 2 〉 = σ coh + σ incoh .
Contrary to the instantaneous RCS, the average RCS is expressed as the sum of the coherent and incoherent RCSs. To obtain the average RCS of a random medium, the coherent scattered field and incoherent scattered field are both needed.
- 2. Formulation of Incoherent Method
As mentioned above, the fluctuations in a scattered field are caused by spatial variations of the dielectric constant between the inclusions/spaces
) and the background,
. Such variations operate as the volume current density,
. Thus, the wave equation can be summarized as follows
(6) ∇×∇× E → 1 − κ eff 2 E → 1 = ω 2 μ( ε− ε eff ) E → 1 , =jωμ J → f ,
E → 1
is the electric field in the boundary enclosing the random medium, and
denote the permeability of the equivalent volume and extinction coefficient of the incoherent field, respectively. Additionally,
are the permittivity of the inclusions/space and the effective permittivity of the equivalent homogeneous volume, respectively. The volume current density
is the current density induced by the incident field on the equivalent volume. That is,
can be considered an equivalent current source that generates the incoherent scattered field. If the effective permittivity of the random medium is defined using a non-dimensional parameter, such as
, then we have
(7) ε eff = ε 0 ( 1−js ).
Then, the volume current density
(8) J → f =jω( ε eff − ε 0 ) E → 1 =ω ε 0 s E → 1 .
The magnetic vector potential,
, based on the volume current density
(9) A → = μ 4π ∭ V J → f ( x ′ , y ′ , z ′ ) e −jkR R d v ′ ,
is the distance between the source and the observation point. We assume that the observation point is in the far zone from the target, which is located at the origin of the coordinate system. Using this approximation in the integrals for vector potential
, the incoherent scattered field originating from volume current density
can be written as
(10) E incoh s = −jωμ 4πr e −j k 0 r × ∭ V [ J → f ( x ′ , y ′ , z ′ )−( J → f ( x ′ , y ′ , z ′ )⋅ r ^ ) r ^ ] e −jk r ¯ ′ ⋅ r ^ d v ′ +O( 1 r 2 ).
Consider that an equivalent homogeneous volume of width
, and thickness
is illuminated by a uniform plane wave, as shown in
. The media is bounded on both sides by air. The propagation direction of the incident plane wave and the direction of the observation point are assumed to be parallel to the
-axis and are defined by
k ^ i = x ^
k ^ s =− x ^
, respectively. Thus, the incident field is given by
(11) E i = E 0 exp( −j k 1 x ).
The ratio between the wave number of the effective medium and that of free space,
, can be made as close to 1 − j
/2 as desired by making the non-dimensional parameter
sufficiently close to zero
. Substituting (8) and (11) into (10), we obtain
(12) E incoh s ≈ −jωμ e −j k 0 r 4πr ∭ V ( ωεs E 0 e −j k 1 x ′ )⋅ e −j k 0 x ′ d v ′ = −j k 0 2 e −j k 0 x s 4πx ∫ 0 h ∫ 0 w ∫ 0 d E 0 ⋅ e −j( k 1 + k 0 ) x ′ d v ′ = −j k 0 2 e −j k 0 x s 4πx ∫ 0 h ∫ 0 w ∫ 0 d E 0 ⋅ e − k 0 ( 2j+ s 2 ) x ′ d v ′ .
We suppose that a number of scatterers,
, divide the equivalent homogeneous volume into
-directions, as shown in
, such that
equals the total number of total subvolumes (
). The integral in (12) can be expressed in the following Riemann summation form:
(13) E incoh s =− j k 0 2 e −j k 0 x s 4πx ΔxΔyΔz ∑ n h =1 N h ∑ n w =1 N w ∑ n d =1 N d e − k 0 ( 2j+ s 2 )mΔx ,
, and Δ
, respectively. In addition, Δ
is the reciprocal of the scatterers’ density. Then, (13) becomes
(14) E incoh s =− j k 0 2 e jkx s 4πxρ ( N h ⋅ N w ) ∑ n d =1 N d e − k 0 ( 2j+ s 2 )mΔx .
Substituting (14) into (5) leads us to the following incoherent RCS:
(15) σ incoh = 2 π 2 λ sA ρ ( 1− e − k 0 sd ).
Therefore, (15) provides the incoherent RCS of a random medium; this incoherent RCS is defined by the shape of an analysis volume, effective permittivity, and inclusion density at the desired frequency. Equation (15) does not have terms related to the inclusions; the properties of inclusions is reflected in the effective permittivity of the random medium.
III. Simulation Results
To validate the proposed method, simulation results are compared with the incoherent RCS obtained from the Monte-Carlo analysis using the MoM with 50 realizations, employing MATLAB R2013b with a 2.83 GHz Quad CPU and an 8 GB RAM PC. The coherent RCS, meanwhile, is obtained from the GEC method
and the Monte-Carlo analysis. The multiple identical scatterers are assumed to be thin, perfectly conducting wire of a half-wavelength in length. These scatterers are also assumed to be uniformly distributed in volume
), as shown in
. The scatterers have two orientation distributions — a uniformly random orientation for all directions and a horizontal orientation (parallel to the horizontal plane). The cross-sectional area of the random medium is fixed as
= 10λ × 10λ, and its depth,
, varies from 0.1λ to 4λ. Although the total volume
varies, the density of the scatterers is kept at 1.0 [no./λ
]; thus, the total number of scatterers varies.
Geometry of equalized homogeneous media of w × d × h (m3) with N sub-blocks.
As mentioned before, each realization of the Monte-Carlo analysis uses the MoM. For the
th realization, the result when using the MoM is the scattered field of (1). Thus, by substituting the
th scattered field (1) into (4), the instantaneous RCS is obtained for the
th realization. After all realizations, the average RCS is calculated by averaging all the instantaneous RCSs. The coherent field is also calculated by averaging the scattered field over a total of
realizations, as described in (2); then, the coherent RCS is obtained using (5). Therefore, the incoherent RCS can be obtained by subtracting the coherent RCS from the average RCS.
show the backscattering incoherent, coherent, and average RCSs of the random medium normalized to λ
, in terms of the resonance frequency of the scatterers as a function of the depth of the slab. In terms of the values of peak, null, and level, the proposed method shows good agreement with the Monte-Carlo analysis performed using the MoM.
The coherent RCSs shown in
have nulls at points where the depth of media are multiple half-wavelengths, but maintain constant peak value. When a plane wave is incident normally on the random medium, the phase difference between the fields reflected from the front and rear of the equivalent media generates nulls and peaks. These patterns mean that the coherent RCS can be calculated from the equivalent homogeneous medium of the random medium, as mentioned earlier. On the other hand, as the number of scatterers increases, the incoherent RCS also increases. The incoherent RCS even becomes larger than the coherent RCS when there is an increase of the depth of the equivalent homogeneous medium. Therefore, as illustrated in
, the average RCS, which is the sum of the coherent and incoherent RCSs, has an increasing value with oscillating form as a function of depth.
The scatterers in
orient uniformly in all directions, whereas those in
have a horizontal orientation. In other words, all scatterers are parallel to the
-plane. For a uniformly oriented distribution, the
θ ^ θ ^
-polarization RCS is equal to the
θ ^ θ ^
-polarization RCS. Otherwise, for a perfectly horizontal orientation, the
θ ^ θ ^
-polarization average, coherent, and incoherent RCS are about 3 dB larger than those of the scatterers with a uniform orientation, and the
θ ^ θ ^
-polarization RCS has a very small level of less than −600 dBλ
. These results can be also confirmed through the relative effective permittivity of
. Similar to the monostatic RCS,
is identical to
for a uniform orientation. In the case of horizontal orientation,
of 1.0-j0.0 means that the scatterers have no electrical effect on a
Monostatic RCS of scatterers with uniform orientation: (a) incoherent RCS and coherent RCS, and (b) average RCS.
Monostatic RCS of scatterers with horizontal orientation: (a) coherent RCS and incoherent RCS, and (b) average RCS.
Relative effective permittivity with respect to wave polarization and orientation of scatterers.
| ||Wave polarization |
|εθθ ||εϕϕ |
|Uniform orientation ||0.9892 - j0.0195 ||0.9892 - j0.0194 |
|Horizontal orientation ||1.0 - j0.0 ||0.9842 - j0.0285 |
To illustrate the differences between the calculated RCSs and the reference RCSs, the normalized mean square error (NMSE) is often used; it can be defined as
(16) NMSE= ∑ n=1 N | RCS n p − RCS n r | 2 ∑ n=1 N | RCS n r | 2 ,
RCS n p
RCS n r
are the linear scaled RCS estimated using the proposed method and the Monte-Carlo analysis, respectively, at the
shows the NMSE for the proposed method using the results of the Monte-Carlo analysis as the reference values for uniformly oriented scatterers. It can be seen that the proposed method has a very small NMSE and provides a very accurate value. In particular, the proposed method provides a more accurate solution than does the GEC method in terms of the NMSE.
NMSE between proposed method and Monte-Carlo analysis for uniform orientation distribution.
|Method ||Monte-Carlo results as reference |
|Type ||NMSE |
|Proposed ||Incoherent RCS ||0.0331 |
|GEC ||Coherent RCS ||0.1172 |
|Proposed & GEC ||Average RCS ||0.0277 |
shows the backscattering coherent, incoherent, and average RCSs for several inclusion densities. The coherent RCS is calculated using an equivalent homogeneous medium, and is significantly affected by the effective permittivity, which is mainly dependent on the inclusion density. The peak values as a function of the inclusion density vary from −30 dBλ
to 30 dBλ
. When the depth of the medium is increased, the peak level decreases slightly and the null level increases slightly. The cause of these behaviors is the attenuation of the reflected field after multiple bounces within the equivalent homogeneous medium. In the same vein, the attenuation in a random medium with a high inclusion density also increases. Therefore, the ripples of the coherent RCS decrease with increases in the depth and inclusion density. These tendencies are clearly confirmed at a higher density.
The incoherent RCS curves of
are asymptotic, not oscillatory. For a higher inclusion density or greater depth, the incoherent RCS has a tendency to rapidly converge to a certain level, as indicated in (15).
Monostatic RCS of scatterers with uniform orientation with respect to inclusion density: (a) coherent, (b) incoherent, and (c) average RCS.
When the number of scatterers increases, spaces with difference in the effective permittivity and the scatterers’ permittivity widen. This situation implies an increase in the incoherent scattering. If the number of scatterers is even greater, the wider space causes a difference between the effective permittivity and the scatterers’ permittivity. Owing to the high effective permittivity, the difference is lowered. This means a convergence at a particular level. This convergence can be traced to the phenomenon in which a low current is induced to scatterers with small gaps between them. In other words, the problem is whether the predominant scattering is generated by the inter-scatterer coupling or the shielding effect.
From the pattern of the average RCS in
, we can ascertain the dominant term with respect to the scatterer density. When the scatterer density is less than 1.0 [no./λ
], the incoherent RCS is dominant. On the contrary, when the scatterer density is larger than 1.0 [no./λ
], the coherent RCS is dominant. The reference density that determines the dominant factor depends on the shape, orientation, and material properties of the scatterers.
This paper proposed a method for calculating the incoherent scattered intensity and RCS of random media. The proposed method was expressed as a simple formula, which can easily and quickly calculate the incoherent RCS and scattered intensity without multiple realization or complex calculation procedures. Therefore, it is expected that this method will be very useful in scattering analysis and RCS prediction of random media. To verify the proposed method, we compared our results to results obtained using a Monte-Carlo analysis with the MoM. The proposed method is in good agreement with the Monte-Carlo analysis. In addition, because both the coherent RCS and the incoherent RCS are directly estimated based on the effective permittivity of a random medium, the proposed method is efficient in comparison with conventional methods, such as the T-matrix and the radiative transfer method. The proposed method can also be applied to inclusions with not only a fiber shape but also various other shapes. Additionally, the proposed method can be also expanded to multi-environments in which the densities of scatterers are varying spatially, if the effective permittivity and analysis geometry are known.
This work was supported by the MSIP (the Ministry of Science, ICT and Future Planning), Rep. of Korea (No. 10041876).
Dong-Wook Seo received his BS degree in electrical engineering from Kyungpook National University, Daegu, Rep. of Korea, in 2003 and his MS and PhD degrees in electrical engineering from the Korea Advanced Institute of Science and Technology, Daejeon, Rep. of Korea, in 2005 and 2011, respectively. He was a senior researcher at the Defense Agency for Technology and Quality, Daegu, Rep. of Korea, from 2011 to 2012. Since 2012, he has been a senior researcher at ETRI, Daegu, Rep. of Korea. His current research interests include numerical techniques in the area of electromagnetics, radar cross-section analysis, wireless power transfer, biomedical implantable devices, and automotive radar systems.
Corresponding Author email@example.com
Jae-Ho Lee received his BS degree in electronic and electrical engineering from Kyungpook National University, Daegu, Rep. of Korea, in 2002 and his MS degree in electrical and electronic engineering from the Korea Advanced Institute of Science and Technology, Daejeon, Rep. of Korea, in 2004. He was awarded his PhD degree in electrical and electronic engineering from the Tokyo Institute of Technology Japan, in 2010. From 2004 to 2005, he worked for the Mobile Communication PM team, Institute of Information and Technology Assessment, Daejeon, Rep. of Korea. He worked for the Research Center at Samsung Thales, Yongin, Rep. of Korea, from 2010 to 2012. Since 2013, he has been a senior researcher at ETRI, Daegu, Rep. of Korea. His research interests include waveguide arrays, electromagnetic numerical analysis, biomedical implantable devices, automotive radar systems, and antennas.
Hyung Soo Lee received his BS degree in electrical engineering from Kyungpook National University, Daegu, Rep. of Korea, in 1980 and his PhD degree in IT engineering from Sungkyunkwan University, Suwon, Rep. of Korea, in 1996. Since 1983, he has been a principal researcher at ETRI, Daegu, Rep. of Korea. His research interests are spectrum engineering, WPAN system design, and biomedical IT convergence devices.
“Performance Analysis of the Clutter Map CFAR Detector with Noncoherent Integration,”
DOI : 10.4218/etrij.93.0193.0001
“Finite-Difference Time-Domain Simulation of Scattering from Objects in Continuous Random Media,”
IEEE Trans. Geosci. Remote Sens.
DOI : 10.1109/36.981359
“The Application of Recursive Aggregate T-Matrix Algorithm in the Monte Carlo Simulations of the Extinction Rate of Random Distribution of Particles,”
DOI : 10.1029/94RS01567
“Scattering Properties of Dense Media from Monte Carlo Simulations with Application to Active Remote Sensing of Snow,”
DOI : 10.1029/96RS00939
“The Average Backscattering Cross Section of Clouds of Randomized Resonant Dipoles,”
IEEE Trans. Antennas Propag.
DOI : 10.1109/TAP.1974.1140813
“Monte Carlo Simulations of the EM Bistatic Scattering from a Novel Foil Cloud,”
Int. Symp. Antennas Propag. EM Theory
Aug. 15–18, 2000
“Maxwell-Garnett Mixing Rule in the Presence of Multiple Scattering: Derivation and Accuracy,”
Phys. Rev. B
DOI : 10.1103/PhysRevB.72.014205
“Theory of Microwave Remote Sensing,”
New York, USA
“Generalized Equivalent Conductor Method for a Chaff Cloud with an Arbitrary Orientation Distribution,”
Progress Electromagn. Res.
DOI : 10.2528/PIER10051202
“The Effect of Fiber Orientation Distribution on the Effective Permittivity of Fiber Composite Materials,”
J. Electromagn. Waves Appl.
DOI : 10.1163/156939310793675835
“Electromagnetic Wave Propagation through Chaff Clouds,”
IEEE Trans. Antennas Propag.
DOI : 10.1109/TAP.2007.900178
“Multiple Scattering of Waves by a Half-Space of Distributed Discrete Scatterers with Modified T-Matrix Approach,”
J. Phys. D: Appl. Phys.
DOI : 10.1088/0022-3727/27/3/002
“Incoherent Scattering from Dense Clouds of Wire Dipoles,”
IEEE Antennas Propag. Soc. Int. Symp.
San Diego, CA, USA
July 5–11, 2008
“Numerical Scattering Analysis for Two-Dimensional Dense Random Media: Characterization of Effective Permittivity,”
IEEE Trans. Antennas Propag.
DOI : 10.1109/8.575638
“Wave Propagation and Scattering in Random Media,”
New York, USA