^{1}

^{1}

^{2}

^{1}

^{1}

A novel modified four path method (FPM) is presented for calculating coupling field of super-low altitude aircraft target. Based on the hybrid method PO + MEC (Physical Optics and Method of Equivalent Currents), the antenna radiation pattern is introduced to consider the multipath interference from side lobe of seeker. The modified FPM is used to calculate the coupling field from super-low altitude aircraft target above different terrestrial environments. The curves of scattering coefficient are analyzed. The influences of height of target, root mean square (RMS), and incident angle on coupling field characteristics are discussed. The simulation results can be used for reference in detection for super-low altitude target and optimization for radar system.

The scattering characteristics of aircraft target above terrestrial environments are complex and variable. The scattering characteristics often changes with many factors such as direction of detection, position of target and parameters of environment. The composite scattering problems [

This paper attempts to establish a kind of scattering model for super-low altitude aircraft target. Based on the traditional FPM, the Gauss antenna radiation pattern is introduced to modify the FPM model. The hybrid method PO + MEC is applied to calculate the scattering from facets and edges of target. Some examples are given by using the modified FPM, and some curves of coupling field with some factors are analyzed. From the numerical results, the characteristics of coupling field can be concluded which can be used for reference in detection for super-low altitude target and optimization for radar system.

The detection for super-low altitude aircraft target above environment belongs to the composite scattering problem of computational electromagnetics. The FPM proposed by J. T. Johnson is used to analyze the scattering from simple target above plane. However, traditional FPM cannot count the multipath interference which enters the radar/seeker through side lobe. The scattering energy of multipath interference which the radar/seeker receives through main lobe can have a stronger effect than that through side lobe. Therefore, aimed at counting the effect of antenna radiation pattern, the Gauss antenna radiation pattern is introduced. We set one main lobe and two side lobes beside it respectively. And we set the gain of side lobe as g_{3}. The function of Gauss antenna radiation pattern is taken in the form:

G M B ( θ ) = exp ( − 2 ln 2 θ 2 / θ B 2 ) , | θ | ≤ μ (1)

where θ B denotes half-power beam point width; μ denotes azimuth angle when the main lobe’s gain is set as g 3 which is defined as follow:

μ = θ B ln g 3 / ( − 2 ln 2 ) (2)

The antenna radiation pattern of first side lobe can be described by the Gauss function at central position of ± 1.5 μ , which is taken in the form:

G B 1 ( θ ) = g 1 exp ( − 2 ln 2 ( θ + 1.5 μ ) 2 θ B 1 2 ) , μ ≤ | θ | ≤ 2 μ (3)

where g_{1} is the peak gain of the first side lobe; θ B 1 is taken in the form:

θ B 1 = 0.5 μ ( − 2 ln 2 ) / ln ( g 3 / g 1 ) (4)

The antenna radiation pattern of second side lobe can be described by the Gauss function at central position of ± 2.5 μ , which is taken in the form:

G B 2 ( θ ) = g 2 exp ( − 2 ln 2 ( θ ± 2.5 μ ) 2 θ B 2 2 ) , 2 μ < | θ | ≤ 3 μ (5)

where

θ B 2 = 0.5 μ ( − 2 ln 2 ) / ln ( g 3 / g 2 ) (6)

In this paper, he peak gain of the first side lobe is G 1 = 10 log 10 g 1 = − 20 dB ; The peak gain of the second side lobe is G 2 = 10 log 10 g 2 = − 25 dB ; The peak gain of other side lobe is G r = 10 log 10 g 3 = − 30 dB ; The corresponding gain is g 1 = 0.01 , g 2 = 10 − 2.5 = 0.00316 and g 3 = 0.001 respectively.

In

{ h = R cos θ 1 + H D = R sin θ 1 θ 2 = θ 3 = tan − 1 [ D / ( H + h ) ] θ 4 = θ 1 (7)

In Tradition FPM, the surface with small roughness is simplified as infinite smooth plane, and the reflection energy is approximate to concentrating in the direction of specular reflection. However, when surface becomes rougher, the diffuse reflection will become dominant. Hence, we use the modified complex reflection coefficient to describe the effect of environment which is taken in the form:

ρ ( θ i s u r ) = Γ ( θ i s u r ) ρ s ( θ i s u r ) (8)

Here θ i s u r is the incident angle illuminating the environmental surface; Γ denotes the Fresnel reflection coefficient; ρ s denotes the reflection factor of rough surface. We set the polarization mode as Vertical-to-Vertical (VV) in this paper.

Set E ( θ i , θ s ) and E ′ ( θ i , θ s ) as scattering field from target and mirror target respectively. θ i and θ s denote incident angle and scattering angle respectively. According to

1) The first path is “seeker → target → seeker”. The scattering field of target is E 1 ( θ 1 , θ 1 ) .

2) The second path is “seeker → target → environment → seeker”. The equivalent path is “mirror seeker → mirror target → environment → seeker”. The scattering field of mirror target is ρ ( θ 2 ) E ′ 2 ( π − θ 4 , θ 2 ) .

3) The third path is “seeker → environment → target → seeker”. The equivalent path is “mirror seeker → environment → target → seeker”. The scattering field of target is ρ ( θ 2 ) E 3 ( π − θ 3 , θ 1 ) .

4) The fourth path is “seeker → environment → target → environment → seeker”. The equivalent path is “seeker → environment → mirror target → environment → seeker”. The scattering field of mirror target is ρ 2 ( θ 2 ) E ′ 4 ( θ 2 , θ 2 ) .

According to the scattering mechanism in

ξ ( Δ θ ) = G ( Δ θ ) G M B ( θ ) | θ = 0 ∘ = { G M B ( Δ θ ) , | Δ θ | ≤ μ G B 1 ( Δ θ ) , μ ≤ | Δ θ | ≤ 2 μ G B 2 ( Δ θ ) , 2 μ ≤ | Δ θ | ≤ 3 μ (9)

Based on traditional FPM model in

E s = E 1 ( θ 1 , θ 1 ) + ξ ( Δ θ ) ρ ( θ 2 ) E ′ 2 ( π − θ 4 , θ 2 ) + ρ ( θ 2 ) E 3 ( π − θ 3 , θ 1 ) + ξ ( Δ θ ) ρ 2 ( θ 2 ) E ′ 4 ( θ 2 , θ 2 ) (10)

Here E s denotes the electric field vector which the seeker receives in total. Except for Path 1, the sum of electric field vector for other three paths is called as multipath interference, which has great effect on detection for super-low altitude target.

We establish a simulation model for a kind of aircraft and subdivide it into many triangle facets by using Hypermesh software. The PO method [

E facet ≈ i k 0 4 π e − j k 0 r r [ s ^ × ( M + η 0 s ^ × J ) ] ⋅ Δ A ⋅ I (11)

Here J and M denotes electric current and magnetic current of the facet respectively; Δ A denotes the area of triangle facet. Based on the Gardon integral, I is taken in the form:

I = 1 j k 0 | n ^ × w | Δ A ∑ m = 1 3 ( n ^ × w ) ⋅ a m exp ( j k 0 r m ⋅ w ) sin c ( 1 2 k 0 a m ⋅ w ) (12)

In Equation (12), w = s ^ − i ^ , sin c ( x ) = sin x / x ; a m is a vector of the m-th edge structure; r m is position vector for midpoint of the m-th edge structure; i ^ and s ^ denotes unit vector of incident and scattering direction respectively; k 0 and η 0 are wave number and wave impedance in free space respectively; When | n ^ × w | = 0 , Equation (12) can be written as follow:

I = exp ( j k r 0 ⋅ w ) (13)

In Equation (13), r 0 denotes position vector of arbitrary point belonging to the facet i. In this paper, we set that point as mass center of the facet i. Finally, the far field of all facets can be obtained by adding every facet’s far field.

The aircraft model contains several edge structures. The scattering from these edge structures cannot be neglected. We introduce the theory of equivalent electric current [

E e d g e = i k 0 4 π e − j k 0 r r ∫ C [ η 0 s ^ × ( s ^ × J ) + s ^ × M ] e i k 0 s ^ ⋅ r ′ d l (14)

In which J and M present equivalent electric current and equivalent magnetic current of edge C respectively; s ^ denotes unit vector of observation direction; r ′ denotes a vector from original point to certain point located in edge C; dl denotes arc length increment of edge C.

In this section, the characteristics of coupling field from super-low altitude aircraft target above different terrestrial environments are analyzed by modified FPM. The influences of height of target, root mean square (RMS), and incident angle on coupling field characteristics are discussed.

Firstly, we explain the simulation conditions of following computation examples: working frequency is set as 10 GHz; Polarization mode of illuminating wave is VV polarization; Relative dielectric constant of concrete surface equals to 3.0−j0.00036; Relative dielectric constant of dry land surface equals to 4.0−j0.00018; Relative dielectric constant of wet land surface equals to 30.0−j0.36; The root mean height (RMS) of rough surface is 0.1λ in all computation examples except for example 4.3; The target model is set as perfect conductor aircraft; The nose of aircraft is facing the positive direction of x axis; The height of aircraft is 15 min all computation examples except for example 4.4; The spatial distance between seeker and target is set as 5 km; The half-power beam width of seeker is 5˚; The observation angle ranges from 0˚ to 90˚ (corresponding grazing angle ϕ s ranges from 90˚ to 0˚).

Based on the modified FPM, the curves of coupling field from super-low altitude aircraft target above different terrestrial environments are obtained as follow.

In

In this computation example, the azimuth distributions of coupling field with different incident angle are calculated. Set terrestrial surface as concrete surface. Set the incident angles θ_{i} as 45˚, 60˚ and 75˚ (corresponding grazing angles ϕ in are 45˚, 30˚ and 15˚). The curves of azimuth distribution with different incident angle are given as shown in

In this computation example, the influences of RMS on coupling field are discussed. Set terrestrial surface as concrete surface. Set the RMS h_{rms} as 0.1λ, 0.3λ and 0.5λ. The curves of coupling field with different RMS are as shown as follow.

In

In this computation example, the influences of target height on coupling field are discussed. Set terrestrial surface as concrete surface. Set the target height as 15 m, 30 m, 45 m, 60 m, 75 m and 90 m. The curves of coupling field with different target height are given as follow.

In

This paper attempts to establish a kind of scattering model for super-low altitude aircraft target. Based on the traditional FPM, the Gauss antenna radiation pattern is introduced to modify the FPM model. The hybrid method PO + MEC is applied to calculate the scattering from facets and edges of target. The characteristics of coupling field from super-low altitude aircraft target above different terrestrial environments are analyzed by modified FPM. The influences of height of target, RMS and incident angle on coupling field characteristics are discussed. Four important conclusions are obtained: firstly, if the seeker detects target at the Brewster angle, the multipath interference will be weakened in azimuth distribution; secondly, with the increase of RMS, the coupling field of almost angles decreases obviously; thirdly, with the increase of target height, the coupling field of almost angles decreases gradually; fourthly, with the increase of target height, the corresponding angle of minimum value for coupling field decreases. The simulation results can be used for reference in detection for super-low altitude target and optimization for radar system.

The authors declare no conflicts of interest regarding the publication of this paper.

Zou, G.X., Tong, C.M., Yu, B.L., Sun, H.L. and Song, T. (2018) A Modified Four Path Method for Calculating Coupling Field of Super-Low Altitude Aircraft Target. Journal of Computer and Communications, 6, 64-73. https://doi.org/10.4236/jcc.2018.612006