Volume 4, Issue 1, January 2016, Pages: 1-7

Nonlinear Time-Varying Turbulence Signal Processing and Simulation Under Curvilinear Coordinate System

Xiaoyang Liu1, 2, Wanping Liu1, Chao Liu1, Xiaoping Zeng2

1School of Computer Science and Engineering, Chongqing University of Technology, Chongqing, China

2Postdoctoral Research Station of Information and Communication, Engineering, Chongqing University, Chongqing, China

Email address:

(Xiaoyang Liu)

To cite this article:

Xiaoyang Liu, Wanping Liu, Chao Liu, Xiaoping Zeng. Nonlinear Time-Varying Turbulence Signal Processing and Simulation Under Curvilinear Coordinate System. Communications. Vol. 4, No. 1, 2016, pp. 1-7. doi: 10.11648/

Received: April 16, 2016; Accepted: April 25, 2016; Published: May 13, 2016

Abstract: In order to analyse the nature of the turbulence, the curvilinear coordinate system is established according to the relations of cartesian coordinate system. Aiming at this turbulence mechanism, the turbulence boundary equations were established, then all grid points of the coordinates can be got by the established equations. The dangerous degree measure (risk factor) was put forward. The mathematical model of risk factor was established as to quantitatively describe the danger level of turbulence. A turbulence signal processing algorithm is proposed combined with the Doppler effect and the computing grid. The simulation results show that the curvilinear coordinate system and three-dimensional turbulence field can reflect the characteristics of turbulence, and risk factor can reflect the danger level of turbulence, the proposed turbulence signal processing algorithm can detect and forecast the turbulence effectively.

Keywords: Turbulence, Curvilinear Coordinate System,Grid Generation, Spectrum Width, Wind Velocity

1. Introduction

Turbulence has been recognized as a cause of aircraft crashes for some time.

Turbulence refers to the wind speed fluctuations in a short period time. Evoked a lot of In order to understand the nature of the proposed turbulence model, a description of microbursts is necessary. The necessary description and the characteristics of the JAWS (Jamming and Warning System) data sets are included in this paper.

These descriptions contain a discussion of some turbulence models. In anticipation of the need for adding turbulence to the JAWS data for the purpose of flight simulation, relevant aspects of turbulence theory are presented [1, 2, 3, 4]. The rule of turbulence motion is irregularity and instability. As for the weather radar, turbulence is the meteorological goal of large particles which speed suddenly change greatly.

The wind loading of a tracking radar antenna, for example, results from a meanvelocity component that varies with time plus superimposed random gusts. In the turbulence region, because of the great variation of river air velocity along the direction of change is quite sharp. Cases like this not only create dangers for the airplane but in danger of life [5, 6, 7].

The classic turbulence models of Von Karman and Dryden has been attracting the people’s attention. But these models do not consider the influence of asymmetric wind field. Kolmogorov and other people have realized the no characteristic scales of turbulence [8, 9, 10]. They think that the number of Reynold is large, turbulence is made up of all kinds of different scales of vortex. Maelstrom is produced by boundary condition and instability. Maelstrom is broken into smaller spiral by inertia and energy delivery. The difficulty of turbulence research is that it contains large and small scales of turbulence vortex.

Due to the scale self-similarity of the wavelet transform. Bezerra and Ashrafian used the method of wavelet transform to study the turbulence, its main idea is to use of the domain is transform the space domain into  domain by using Fourier transform. Looking for  value when the signal is the strongest, and to determine the main ingredients or cycle of the signal. So, they can decide the signal is strong or not, or in which scale the signal is the strongest [11, 12]. This approach does not play the advantages of wavelet transform. Because the wavelet transform is related to the time factor, but Fourier transform doesn't have the time factor.

Turbulence is a highly complicated nonlinear dissipative system. The MLSD (Maximum-likelihood Sequence Detection) is realized in lognormal statistical model of the atmospheric turbulence by Susumu. G and Philippe R [13, 14]. But MLSD cannot track the time change of the multiplicative noise statistics. The influences of turbulence detection is researched with PPP (pulse pair processing) algorithm by Y. Li and V. Raj Mohan [15, 16]. References [17, 18, 19, 20] researched the autocorrelation characteristics of the Von Karman model, but they have not considered the effect of turbulence scale factors. References [21, 22, 23] mainly studied the numerical simulation of atmospheric turbulence, Gaussian noise sequence is dealt by adopting a double random exchange of minimization method. The method will vastly improve the degree of bleaching, but there is no discussion the wind speed estimation of turbulence flow field. Those studies have been important for turbulence detecting, and in particular their the frequency spectrum analysis of turbulence [24, 25].

Depending upon the root mean square of turbulence Doppler velocity , turbulence level classification can be defined in Table 1.

Table 1. Turbulence level classification.

turbulence intensity negligible slight moderate strong

0-1.5 1.5-3 3-4.5 >4.5

The sketch of turbulence observation under airborne weather radar is shown in Figure 1.

Figure 1. The sketch of turbulence observation.

In Figure 1. The plane fly at a constant velocity  along X axis. The antenna angle of view is . from the scatterer S return of the received signal The Doppler frequency shift of the received signal from the scatterer S is described by


Where  is the wavelength of airborne weather radar.  is the azimuth of Scatterer S.

2. Curvilinear Coordinate System

For convenience, cartesian coordinate system  is written as , and cylindrical coordinate system  is written as . The mutual conversion relationship of the two coordinate systems can be expressed as



It can be seen that the relations between the two set of functions is single-valued continuous function. The two Jacobian determinant of coordinate transformation are not zero, in addition to the coordinate. In a certain space area, Jacobian determinant can be described by


In the region D, there is one-to-one correspondence between a dotand an array . Take  as the coordinate of the dot. So, a curvilinear coordinate system is established [26, 27, 28], which is shown in Figure 2.

Figure 2. Curvilinear coordinate system in region D.

3. The Generation of Turbulence Calculation Grid

Now make  represents the quality of the fluid micelle, shows the average density of the fluid micelle, expresses the volume of the fluid micelle.  is absolute velocity.

Based on the mass conservation law, the energy conservation law and the definition of divergence. When the velocity field is continuously and differentiable, we have


The motion equations of fluid micelle is described by




In the above formulas,  is the mass force of the unit mass fluid and inertial force.  expresses the stress tensor.  denotes the viscous force.  is a proportionality constant associated with the viscosity of fluid.  is a unit vector.

Under the  curve coordinate system, a two-dimensional curve coordinate system can be derived according to the fluid mechanics in combination with Lame coefficient.


A set of equations can be worked out of meeting with some formulas and other boundary conditions.


The coordinates of all grid points  can be attained associated with . Then connect these grid points by the smooth curve, the grid lines of curve coordinate can be got by connecting these grid points with the smooth curve.

4. Turbulence Signal Processing and Detecting

4.1. The Risk Factor Model

To conduct quantitative analysis on the risk of turbulence, a risk factor for danger judgment is proposed. According to the theorem of quality continuity, a risk factor can be defined when the vertical wind speed is proportional to the horizontal divergence in the given the height.


In the above formula,  is the ground speed of aircraft. .  expresses the horizontal wind velocity.  is the rate of change of horizontal wind velocity. denotes the height of aircraft.  is the distance from the airplane to the center of the micro air flow. The positive value of risk factor expresses the descended performance of airplane. The negative value of risk factor denotes the increased performance of airplane.

4.2. Three-Dimensional Turbulence Field

In order to meet the demands of numerical simulation. The turbulence field in the space is divided into  grid points. The Three-dimensional turbulence field is shown in Figure 3.

Figure 3. Three-dimensional turbulence field.

According to the nature of FFT transform, the two relations between the space sampling rate of time domain and the sampling rate of frequency domain can be described by



In the above formulas,  is the sampling rate of time domain in the direction .  is the number of grid points in the direction .  is the grid spacings of frequency domain in the direction .  is the grid spacings of time domain in the direction .

Figure 4. The generation of turbulence field.

In Figure 4,  is the three road of Gaussian white noise sequence, .  denotes the three dimensional formed filter function. IDFT expresses the inverse transform of three dimensional discrete Fourier.  shows the isotropic turbulence.

4.3. Turbulence Signal Processing

The three dimensional discrete Fourier transform and its inverse transform can be expressed as



In the above formulas, ; .  is the numerical sequence of time domain of three dimensional Fourier transformation.  is the numerical sequence of frequency domain of three dimensional Fourier transformation.

The time domain noise variance of the three dimensional zero mean is . So the frequency domain noise variance of the three dimensional zero mean can be described by


The following formulas can be got by the FFT of three-dimensional symmetrical characteristic.





In turbulence wind field, the echo signal sampling sequence of a certain range bin can be described by


and represent the in-phase and quadrature components of pulse sampling signal, respectively.

The power spectrum of echo signal can be got


In the above formula. is the length of sampling sequence. denotes the sampling interval of echo signal.  expresses the amplitude.

The average Doppler speed  and spectrum width  of echo signal can be described by



The echo signal of turbulence processing is shown in Figure 5.

Figure 5. The echo signal of turbulence processing flow chart.

5. Simulation Results Analysis

The echo signal of three-dimensional power spectrum Doppler velocity of turbulence echo signal is shown in Figure 6. It can be seen that the proposed algorithm can estimate the wind speed distribution of the turbulence wind field effectively.

Figure 6. Three-dimensional power spectrum analysis.

In order to verify the validity of the algorithm, we compare the real wind velocity and the proposed algorithm. As can be seen from the Figure 7, the wind speed can be accurately estimated by the proposed algorithm.

Figure 7. The Doppler velocity in different range bin.

When the turbulence signal is processed by airborne weather radar, firstly estimate the spectrum width of echo signal in wind filed with different range bins, in order to determine whether the target echo is turbulence or not. Figure 8 and Figure 9 are the risk factor analysis of turbulence in different range bins.

Figure 8. The risk factor analysis of turbulence (have ground clutter).

Figure 9. The risk factor analysis of turbulence (no ground clutter).

As can be seen from Figure 8 and Figure 9, the risk factor is different in different range bins, the risk factor reaches the maximum between the 20th rang bin and the 25th range bin. The performance of no ground clutter is better than the situation of have ground clutter through the Figure 8 and Figure 9.

Figure 10. Three dimensional risk factor analysis.

The diagram in Figure 10 shows the three dimensional risk factor in different scanning lines. The Figure11 expresses the risk factor in different range bin. According to risk factor definition, the threshold value of risk factor can be set to 0.2. It is dangerous in some areas as can be seen from Figure 11.

Figure 11. The risk factor in different range bin.

6. Conclusions

In this paper, a curvilinear coordinate system is established in this paper. The mathematical model of turbulence calculation grid is set up combined with the continuous equation and motion equation of turbulence. The grid computing model is good for subdivision turbulence field. The grid could enable detailed simulations of turbulence. In order to be able to quantitative analysis the harm of turbulence, the risk factor of turbulence is set up which reflects hazard of turbulence. A turbulence signal algorithm is proposed combined with Doppler effect. This algorithm can reflect the characteristics of turbulence and can carry on the detection and prediction of turbulence.

This turbulence signal processing model and mathematical model of risk factor can be used as the reference model of microburst wind shear. Considering the nonlinear of turbulence, the space-time adaptive processing and deconvolution method can be used for turbulence research and other characteristics.


The paper are supported by Chongqing Postdoctoral Researchers Research Project Special Funded (Grant No. Xm2015029), Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant No. KJ1500926, KJ1500904, KJ1500920), Chongqing Science and Technology Committee Fund of Foundation and Cutting-edge Research Plan (Grant No. cstc2014jcyjA40007), Natural Science Foundation of China (Grant No. 61502064, 61503052, 91438104, 11547148).


  1. Oreifej. O, Xin-Li, Shah. M, Simultaneous Video Stabilization and Moving Object Detection in Turbulence, IEEE Transactions on Pattern Analysis and Machine Intelligence, 2013, 35 (2): 450-462.
  2. Treib. M, Burger. K, Reichl. F, Meneveau. C, Szalay. A, Westermann. R,Turbulence Visualization at the Terascale on Desktop PCs, IEEE Transactions on Visualization and Computer Graphics, 2012, 18 (2): 2169-2177.
  3. Gil. E, Laguna. P, Martinez. J. P, Barquero-Perez. O, Garcia-Alberola. A, Sornmo. L,Heart Rate Turbulence Analysis Based on Photoplethysmography, IEEE Transactions on Biomedical Engineering, 2013, 60 (11): 3149-3155.
  4. Gibson. K. B, Nguyen. T. Q, An Analysis and Method for Contrast Enhancement Turbulence Mitigation, IEEE Transactions on Image Processing, 2014, 23 (7): 3179-3190.
  5. Xuan Tang, Zhaocheng Wang, Zhengyuan Xu, Ghassemlooy. Z,Multihop Free-Space Optical Communications Over Turbulence Channels with Pointing Errors using Heterodyne Detection, Journal of Lightwave Technology, 2014, 32 (15): 2597-2604.
  6. Rajbhandari. S, Ghassemlooy. Z, Haigh. P. A, Kanesan. T, Xuan Tang, Experimental Error Performance of Modulation Schemes Under a Controlled Laboratory Turbulence FSO Channel,Journal of Lightwave Technology, 2015, 33 (1): 244-250.
  7. Pham. H. T. T, Dang. N. T, Pham. A. T, Effects of atmospheric turbulence and misalignment fading on performance of serial-relaying M-ary pulse-position modulation free-space optical systems with partially coherent Gaussian beam, IET Communications, 2014, 8 (10): 1762-1768.
  8. Puryear. A. L, Shapiro. J. H, Parenti. R. R, Reciprocity-enhanced optical communication through atmospheric turbulence — Part II: Communication architectures and performance, IEEE/OSA Journal of Optical Communications and Networking, 2013, 5 (8): 888-900.
  9. Kaur. P, Jain. V. K, Kar. S, Performance Analysis of FSO Array Receivers in Presence of Atmospheric Turbulence, IEEE Photonics Technology Letters, 2014, 26 (2): 1165-1168.
  10. Xuegui Song, Fan Yang, Julian Cheng, Alouini. M. S, BER of Subcarrier MPSK and MDPSK Systems in Atmospheric Turbulence, Journal of Lightwave, 2015, 33 (1): 161-170.
  11. Vu. B. T, Dang. N. T, Thang. T. C, Pham. A. T, Bit error rate analysis of rectangular QAM/FSO systems using an APD receiver over atmospheric turbulence channels, IEEE/OSA Journal of Optical Communications and Networking, 2013, 5 (5): 437-446.
  12. Katherine Mc Caffrey, Baylor Fox-Kemper, Peter E, Hamlington. Jim Thomson, Characterization of turbulence anisotropy, coherence, and intermittency at a prospective tidal energy site: Observational data analysis,Renewable Energy, 2015, 76 (4): 441-453.
  13. Susumu. G, J. C. Vassilicos,Energy dissipation and flux laws for unsteady turbulence, Physics Letters A, 2015, 379 (16): 1144-1148.
  14. Philippe R, Spalart. Philosophies and fallacies in turbulence modeling,Progress in Aerospace Sciences, 2015, 74: 1-15.
  15. Y. Li, A. M. Castro, T. Sinokrot, W. Prescott, P. M. Carrica, Coupled multi-body dynamics and CFD for wind turbine simulation including explicit wind turbulence, Renewable Energy, 2015, 76, (4): 338-361.
  16. V. Raj Mohan, D. C. Haworth, Turbulence–chemistry interactions in a heavy-duty compression–ignition engine, Proceedings of the Combustion Institute, 2015, 35 (3): 3053-3060.
  17. J. Shinjo, J. Xia, A. Umemura, Droplet/ligament modulation of local small-scale turbulence and scalar mixing in a dense fuel spray, Proceedings of the Combustion Institute, 2015, 35 (2): 1595-1602.
  18. K. Prabu, D. Sriram. Kumar, MIMO free-space optical communication employing coherent BPOLSK modulation in atmospheric optical turbulence channel with pointing errors, Optics Communications, 2015, 343 (15): 188-194.
  19. Erhan Pulat, Hıfzı Arda Ersan, Numerical simulation of turbulent airflow in a ventilated room: Inlet turbulence parameters and solution multiplicity, Energy and Buildings, 2015, 93 (15): 227-235.
  20. Andreas Engelen, Susanne Schmidt, Michael Buchsteiner, The Simultaneous Influence of National Culture and Market Turbulence on Entrepreneurial Orientation: A Nine-country Study, Journal of International Management, 2015, 21 (7): 18-30.
  21. Xiaoyang Liu, Yong Li, Chengyu Feng, Simulation and analysis of turbulence signals in airborne pulse Doppler radar, Systems engineering and electronics, 2015, 21 (1): 18-30.
  22. Xiaoyang Liu, Yong Li, Turbulence signal processing in the airborne weather radar, International Journal of Advancements in Computing Technology, 2013, 34 (5): 816-824.
  23. Yong Li, Xiaoyang Liu, Chengyu Feng, Three dimensional turbulent flow formation and simulation analysis in airborne radar, Systems engineering and electronics, 2013, 35 (6): 1193-1198.
  24. Paolo Orlandi, Sergio Pirozzoli, Matteo Bernardin, George F. Carnevale, A minimal flow unit for the study of turbulence with passive scalars,Journal of Turbulence, 2014, 15 (2): 731-751.
  25. L. Djenidi, S. F. Tardu, R. A. Antonia, L. Danaila, Breakdown of Kolmogorov’s first similarity hypothesisin grid turbulence, Journal of Turbulence, 2014, 15 (3): 596-610.
  26. Xiaoyang Liu, Chao Liu. Wanping Liu. Wind Shear Target Echo Modeling and Simulation. Discrete Dynamics in Nature and Society, 2015 (4): 1-6.
  27. Ugur Cakir, Ertugrul Kargi, Hakan Sarman, Cengiz Isik. Impact of Diabetic Foot on Selected Psychological or Social Characteristics.Journal of Diabetes Research, 2014, 8: 2-9.
  28. Detlef Lohse, Siegfried Grossmann. Intermittency in turbulence. Physica A: Statistical Mechanics and its Applications, 1993, 194 (2): 519-531.

Article Tools
Follow on us
Science Publishing Group
NEW YORK, NY 10018
Tel: (001)347-688-8931