Wind Shear Target Echo Modeling and Simulation

Wind shear is a dangerous atmospheric phenomenon in aviation. Wind shear is defined as a sudden change of speed or direction of the wind. In order to analyze the influence of wind shear on the efficiency of the airplane, this paper proposes a mathematical model of point target rain echo and weather target signal echo based on Doppler effect. The wind field model is developed in this paper, and the antenna model is also studied by using Bessel function. The spectrum distribution of symmetric and asymmetric wind fields is researched by using the mathematical model proposed in this paper. The simulation results are in accordance with radial velocity component, and the simulation results also confirm the correctness of the established model of antenna.


Introduction
As for airborne weather radar, wind shear is a unique atmosphere phenomenon.Microburst is a main form in wind shear.Wind shear often refers to the wind speed or size changing suddenly.Wind shear especially low-level wind shear can cause the maximum damage to airplane.In the aeronautical meteorology, according to the structure of the wind field, wind shear can mainly consist of three kinds of basic situation: vertical wind shear of horizontal wind, horizontal wind shear of horizontal wind, and wind shear of vertical wind.In the actual atmosphere, these three kinds of wind shear can affect flight.According to the airplane relative to the wind vector, the wind shear can be divided into the following: wind shear, cross wind shear, partial wind shear, and head wind shear four forms [1][2][3][4][5].
The scale and strength of low-level wind shear are closely related to wind shear weather system and environmental conditions.Because the phenomenon of wind shear belongs to a small probability event, the existence time of this unique atmosphere phenomenon is only a few minutes [6][7][8][9].It is not repeated.If we rely on the actual test method, not only the cost is very high, but also the risk is quite large, so it is necessary for us to study real wind shear weather change rule of the simulation method [10].
When people realize the serious harm of wind shear, people have invented many equipment which can detect wind field and wind shear, including ground anemometer theodolite, radiosonde wind profile, line radar laser radar, Doppler weather radar, Doppler sound radar, and airborne sensors.The malignant plane crash accident is mainly caused by wind shear [11].Therefore, it is the most effective method for detecting low-level wind shear in time, so that the airplane has enough time and space to avoid it [12].
The testing technology of low-level wind shear is mainly by means of signal processing algorithm, analyzing radar echo signal of wind shear [13].The echo power, wind speed of Doppler, and spectrum width parameters can be extracted by using signal processing algorithm.The wind speed of Doppler is the most important parameter among all the parameters, because it reflects the wind field characteristics of the movement [14][15][16].
The paper is organized as follows.Section 2 contains mathematical model of point target echo, and the signal model of target echo is established in Section 3. Wind field and antenna model are proposed in Section 4, followed by the analysis of simulation results in Section 5. Finally, Section 6 contains the conclusions.

The Proposed Mathematical Model of Point Target Echo
There is a relative motion between airborne weather radar and the meteorological scattering target.So the Doppler frequency shift consists of two parts: the speed of an airplane and meteorological scattering target.For airborne pulsed Doppler radar detecting meteorological target of wind shear, we assume that the speed of airplane is , the average radial velocity among meteorological scattering target is   , the initial distance between the airplane and the meteorological scattering target is , and the distance of airplane away from meteorological scattering unit is ().The relationship between () and time  is as follows: The time rate of the phase  = −4()/ +  is where   is the Doppler frequency shift caused by the airplane and   is the Doppler frequency shift caused by the meteorological scattering target.Equation (2) indicates that the Doppler frequency shift   can be extracted which is caused by the target movement of wind shear rain echo.
The average power of wind shear target rain echo can be expressed as where   is the transmit power,  is the antenna gain,  is the wavelength of airborne weather radar, the scattering rate of meteorological target is , system loss is , and the distance between target and meteorological target is .
We assume that the body reflectivity of meteorological target is ;   represents the target volume unit.
As for airborne weather radar, we assume that the target is the ideal point target; the effective scattering cross section of the target is   .So where the body reflectivity of meteorological target  can be expressed as where  is the wavelength of airborne weather radar. represents the diameter of the target sphere, || 2 is the dielectric constant of scattering particle, and  ( =  6 ) denotes the reflection factor.We assume that the radar transmitted signal is a narrow signal.It can be expressed as where The  is carrier frequency.So, echo signal can be regarded as a signal delay form of transmitted signal.Its amplitude is multiplied by a scale factor .
We assume that the complex reflection constant is .We also defined the scale factor : The echo signal of the mathematical formula can be expressed as where time delay  can be expressed as Putting ( 2) into (10), we can get Then by putting ( 11) into ( 9), we can get where 2/ =  is time delay of target echo.So, for a single wind shear meteorological target scattering body, echo signal of radar target can be expressed as where 2 is the total path length and  is the random phase caused by the scatter body. is the amplitude of echo signal. can be expressed as where  denotes the average power of the wind shear rain echo.

Signal Model of Target Echo
We assume that   denotes the speed component of target scatter body;   represents the speed component of airplane along the radial. = 2/ is constant.We defined Discrete Dynamics in Nature and Society The in-phase component  and quadrature  component of the target echo signal can be expressed, respectively, as where   is the random phase of the scatter target and Δ represents the transmitted phase error.(  ) is the receiver noise. expresses the number of pulses.  is the pulse time interval.
Signal amplitude   can be obtained by radar equation and reflectivity factor.  can be expressed as where where Rae is the constant of radar equation and  loss represents the receiver loss.Vol indicates the volume of the scatter body. shows the multipath fading factor. is the antenna gain.  is wind shear wind field reflectivity.Signal phase is mainly decided by the Doppler frequency shift.Doppler frequency shift consists of two parts,   and   .
The total rain echo signal phase can be expressed as The principle diagram of wind shear target echo simulation is shown in Figure 1.

Wind Field Model of Wind Shear.
In order to simulate the whole space wind field within the scope of the -- three direction of the wind speed, we build a microburst wind shear of wind field model, and the simulation process is as follows: (1) The speed of the vertical flow (a) The calculation of the radial distance where  and  are the airplane's  and  position coordinates and  and  are the center of wind field in  and  position.RC is the radial distance of airplane.The minimum value of RC is set to 2.0.(b) We define a parameter RA; RA is related to the position of airplane, wind shear, and distortion factor where where GX and GY are the distortion factors in  and  position.The minimum value of GR is set to 0.002.The minimum value of RA is set to 2.0.
(2) The radial velocity of horizontal flow The variable of radial velocity distribution is defined by VRR: where  is the height of airplane, VZO is the initial reference velocity, HT is the height limit of horizontal airflow, and GVZ is the gain factor of wind field.
(3) The wind shear model According to the changing rule of the horizontal wind   and vertical wind   , the following wind shear model is established:  In the formula,  and  denote the size of horizontal wind and vertical wind, respectively. 0 = 2/ 0 . 0 shows the total time of the plane through the microburst. 0 expresses the initial phase.To solve the derivative  0 , we also can get

Antenna Model.
The antenna model is made by Bessel function.The gain of main lobe is 34.2 dB.The 3 dB beam width is 2.8 ∘ and the center of first side lobe is 4.4 ∘ and the antenna beam pattern is shown in Figure 2.

Simulation Results Analysis
5.1.The Main Parameters Settings.The main simulation parameter settings are shown in Table 1.

The Condition of No Wind Shear.
The wind speed is set to 0, and the reflectivity is set to 0 dBZ, get rid of the plane, the airplane of the speed along the antenna beam radial velocity takes 0. The speed of Doppler frequency is 0, because there is no the signal of wind shear.The signal of radar received will exist noise and only considers clutter, which is showed in Figure 3.

The Simulation of Symmetric Wind Field.
Considering symmetry wind field, the reflectivity is set to 30 dBZ, because the intensity of rain echo is related to the intensity of the wind field.The center of wind field is large according to the data of wind field model of microburst wind shear, and the intensity of rain echo is large; by contrast wind field edge rain echo is weak, as is shown in Figures 4 and 5. Figures 4 and 5 depict the spectrum distribution of rain echo and range bin change in the first scan line.Figures 6 and 7 depict the spectrum distribution of rain echo and range bin change in the thirtieth scan line.
As can be seen from Figures 4 and 5. the wind speed is almost perpendicular to the radial in microburst wind field edge.The radial velocity component is almost equal to 0. It means that it almost does not produce Doppler frequency shift, and as shown in Figures 6 and 7     relative to the first scan line, the signal strength of rain echo is greater, and the radial velocity component is obvious.

The Simulation of Asymmetric Wind Field.
Considering asymmetry wind field, we assume that the reflectivity is 30 dBZ.The settings of other parameters are showed in Table 1. Figure 8 depicts the spectrum distribution under the head hind, while Figure 9 depicts the spectrum distribution under the partial hind.
Figures 8 and 9 show that the spectrum profile is related to rain echo signal and range bin under the head wind and partial wind in the thirtieth scan line.As can be seen from Figures 8 and 9, the simulation results can reflect the actual rain echo distribution.

Conclusion
In this paper, a novel mathematical model of wind field in wind shear is established.The mathematical model of rain echo point target echo and antenna model are proposed.The Rain Echo Spectrum Distribution under different scan lines is studied in this paper.The simulation realized the symmetric wind field and asymmetric wind field, in a given antenna model and system parameters.It is explained that the data in different coordinate systems can be transformed and an example having analytical solution is given.
The simulation results show that the wind field data of the established model can reflect the basic characteristics of the wind shear in perfect.The distribution of speed spectrum of rain target echo signal can reflect the wind shear condition.The simulation results are in accordance with radial velocity component.

Figure 1 :
Figure 1: The principle diagram of wind shear target echo simulation.

Figure 2 :
Figure 2: The pattern of antenna beam.

Figure 3 :Figure 4 :
Figure 3: Rain Echo Spectrum Distribution of no Wind Shear.

Figure 5 :
Figure 5: Rain Echo Spectrum Distribution of the first scan line (Contour Map).

Figure 6 :Figure 7 :
Figure 6: Rain Echo Spectrum Distribution of the thirtieth scan line (three-dimensional diagram).

Figure 8 :
Figure 8: Rain Echo Spectrum Distribution of asymmetric wind field for head wind (Contour Map).

Figure 9 :
Figure 9: Rain Echo Spectrum Distribution of asymmetric wind field for partial wind (Contour Map).

Table 1 :
The parameters settings of simulation.