Waveform Distortion of Gaussian Beam in Atmospheric Turbulence Simulated by Phase Screen Method

Yang, Zhi-Qiang; Yang, Li-hong; Gong, Lei; Wang, Liguo; Wang, Xinyi

doi:https://doi.org/10.1155/2022/5939293

Journal of Mathematics

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Fuzzy Approximation and its Applications in the Stabilization of Dynamic Systems 2021

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 5939293 | https://doi.org/10.1155/2022/5939293

Waveform Distortion of Gaussian Beam in Atmospheric Turbulence Simulated by Phase Screen Method

Zhi-Qiang Yang,¹Li-hong Yang,¹Lei Gong,¹Liguo Wang,¹and Xinyi Wang²

Academic Editor: Heng Liu

Received30 Nov 2021

Accepted20 Jan 2022

Published21 Feb 2022

Abstract

The atmospheric turbulence phase screen is generated based on the power spectrum inversion method, and the multiple transmission processes are statistically averaged. The waveform distortion of the Gaussian beam in the atmospheric turbulence is analyzed; the simulation results show that the property of Gaussian beam has been destroyed after its passing through the atmospheric turbulence. When the beam waist radius is close to a certain radius, the degree of change becomes larger. When it approaches the critical value, the wave surface no longer changes sharply, and as the turbulence intensity increases, the phase fluctuation becomes more and more severe, the coherence of the beam is destroyed, and the spot may be split into several pieces. Finally, the relationship of intensity fluctuation, amplitude fluctuation, and bit error rate with distance is analyzed.

1. Introduction

The particularity of light propagation in turbulent media lies in the random fluctuation of the refractive index of the medium; the key of numerical simulation for laser atmospheric propagation is to construct a phase screen that can correctly reflect the statistical characteristics of atmospheric turbulence [1–4]. The earliest phase screen simulation turbulence method is proposed by Fleck, who used the “multiphase screen method” in 1976 to simulate the influence of atmospheric turbulence in free space on the transmitted beam. This method is widely used to simulate the atmospheric turbulence effect of the beam [5]. The current methods of generating phase screens are mainly to filter Gaussian white noise by using the refractive index power spectrum of the turbulence, by which the disturbed wavefront phase can be obtained through Fourier transform; this method is called “power spectrum inversion method.” Since Mcglamery proposed the power spectrum inversion method in 1967, this method has been widely used to simulate turbulence [6]. In 1976, Noll proposed to use orthogonal and complete two-dimensional basis functions (such as Zenike polynomials and K-L polynomials) to directly obtain the distorted wavefront of turbulence [7]. In 1983, Walner proposed a method of using the covariance of the phase structure function to generate the Kolmogorov turbulence phase screen [8], but this method is computationally expensive. In 1990, Bowman proposed an algorithm for reconstructing the turbulent wavefront using wavelet transform [9]. Among many methods, the power spectrum inversion method is applicable to a wide range of atmospheric turbulence spectrum models [10–12].

Based on the Gaussian beam horizontal propagation theory, this paper uses MATLAB to carry out detailed simulation on the power spectrum inversion method and analyzes the simulation results. The propagation characteristics of the Gaussian beam under the atmospheric turbulence spectrum are discussed.

2. Theoretical Model of Laser Transmission in Turbulent Atmosphere

Assuming that the monochromatic electromagnetic wave is transmitted along the Z direction, the wave field equation iswhere ω is the frequency of monochromatic electromagnetic waves. When it propagates in the atmosphere (nonmagnetic medium μ = 1) without free charge (ρ₀ = 0) and free current (J₀ = 0), its wave equation can be expressed aswhere is a function of which represents the local refractive index at this point. It is very difficult to accurately solve this basic equation of light transmission in atmosphere.

By removing the vector operator, the equation can be turned into a scalar equation:where is a slowly varying function of z; it can only change at the distance of random medium scale l, so as long as l >> λ; then, , and approximately replaces ( represents the fluctuation of atmospheric refractive index at position ). By simplifying formula (3), the following formula can be obtained:

Equation (4) is called the parabolic equation or quasi-optical approximation relationship, which is suitable for the propagation of narrow-angle expanded beams. k is equal to 2π/λ which represents wave number in free space, and . The infinitely extending random medium can be divided into many sections with thickness of ; the influence of turbulence in each section on the phase of light waves can be concentrated on a thicknessless phase screen, but this phase screen has no effect on the amplitude of the light wave. The change of light wave amplitude is the cumulative result of free-space diffraction at the distance between many two-phase screens.

In order to specify each random phase screen, the correlation function and power spectrum of the phase fluctuation on each phase screen can be calculated. Assuming that there is no diffraction on each phase screen, from the perspective of geometric optics, the correlation function of the phase fluctuation can be written:where is the vector in the coordinate plane perpendicular to the z direction (the xy plane), is the correlation function of the phase fluctuation, and is the correlation function of refractive index fluctuation. Assuming ( is the coherence length of the turbulence), formula (5) can be simplified:

So, the relationship between the power spectrum function and the power spectrum of the turbulent refractive index of each phase screen can be obtained:

Because the transmission of the light wave field between two adjacent phase screens can be approximated as the transmission in free space, which represents the fluctuation of the refractive index , it can be obtained from equation (4):

It can be obtained by Fourier transform of two dimensional on both sides at the same time:where is the wave vector perpendicular to the transmission direction; by integral in the Z direction on both sides of equation (9), solution can be obtained:

Equation (10) shows the free-space diffraction between the two-phase screens.

As shown in Figure 1, after replacing the infinitely extending turbulent medium with a large number of phase screens, the transmission of light waves through the turbulent atmosphere is equivalent to the transmission of light waves between phase screens after passing through the phase screens.(a)The contribution of the phase screen to the phase of the light wave field is as follows: where represents the right side of the ith phase screen and represents the left side of the ith phase screen.(b)Perform Fourier transform on to get .(c)Through the free-space transmission with a distance of ,(d)By performing Fourier transform on , can be obtained:

Then, the transmission of light waves in the entire infinitely extending turbulent medium can be realized by repeatedly using equation (13) until the transmission of the last phase screen is completed.

3. Numerical Simulation Method

The parabolic equation is used to express the light field when the beam propagates in the turbulent atmosphere. The propagation path of the beam can be divided into two parts, the free space and the thin phase screen dispersed in it. When the light beam is transmitted to the phase screen, the wave front of the light is added with the corresponding disturbance part, and it passes through a section of free space to reach the next phase screen. When the beam reaches the receiving end, the entire simulation process ends. From equation (13), the light field at any adjacent phase screen can be expressed by the following formula [13, 14]:

Equation (14) represents the turbulent phase screen, and k_x and k_y represent the wave number. The light wave transmission distance is divided into N_Z sections, and the jth phase screen is in z_j, is the phase screen spacing, and each phase screen is divided into N × N square grids whose grid spacing is .

Power spectrum of the atmospheric turbulence is used to filter the complex Gauss random number matrix; by inverse Fourier transform, the atmospheric disturbance phase will be obtained, which is the core of the power spectrum inversion method.

MATLAB is used to generate a complex Gaussian random number matrix which has a mean value of 0 and a variance of 1, and then, the atmospheric refractive index power spectrum is introduced; the modified Von Karman spectrum is as follows:

The Kolmogorov spectrum is as follows:which is used to produce a continuous phase screen by filtering matrix a with the atmospheric refraction power spectrum [14]:where R represents the spatial spectrum of random function with normal distribution with mean value of 0 and variance of 1 and F represents atmospheric turbulence power spectrum. and are the sampling interval, and the constant C comes from the scale factor, which is used to adjust the phase screen variance. Fresnel diffraction is adopted in the vacuum transmission part, and the light field after the beam transmission distance z can be expressed as [15]where f is the spatial frequency; by Fourier transform of , the transmission in the vacuum part can be realized.

4. Numerical Simulation Results

The simulation parameters are set as follows: wavelength is 1.55 μm, number of sampling points is 256 × 256, phase screen width is 0.5 m, phase screen spacing is 300 m, and 1000 phase screens are constructed.

For the modified Von Karman spectrum, atmospheric turbulence phase screen for the different atmospheric refractive index structure constant and distortion of Gaussian beam after passing through the phase screen are simulated, and the results are as follows.

Figure 2 is the atmospheric turbulence phase screen generated by the power spectrum inversion method. Figures 2(a) to 2(f) are the atmospheric turbulence phase screens when the atmospheric refractive index structure constant takes different values. From figures, it is can be seen that the turbulence intensity increases gradually with the increase of atmospheric refractive index structure constant .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3 shows the light intensity distribution diagram of the initial beam with the beam waist radius . Figure 3 shows Gaussian beam after passing through turbulence simulated by the phase screen theory, where Figures 4(a)–4(f) represent the phase fluctuation and light intensity distribution of the distorted Gaussian beam after passing through phase screens of different intensities.

(a)

(b)

Comparing Figures 4(a), 4(c), and 4(e) with Figure 3(a), after the Gaussian beam passes through atmospheric turbulence, the Gaussian property of the beam is destroyed. Compared with the initial Gaussian beam, the wave surface becomes rugged, the faster the change within the beam waist radius, but when approaching a certain distance from the center, the wave surface no longer changes sharply but up and down changes. As the atmospheric refractive index structure constant increases, that is, the turbulence intensity increases, in fact, the three images correspond to weak turbulence, medium turbulence, and strong turbulence. When the beam waist radius is close to a certain radius, its degree of change keeps getting bigger. When it approaches the critical value, the wave surface no longer changes sharply, but as the distance is close to the center point, the greater the turbulence intensity, the greater the ups and downs of the wave surface. By comparing Figures 4(b), 4(d), and 4(f) with Figure 3(b), it can be obtained after the Gaussian beam passing through the atmospheric turbulence phase screen which is generated by the spectrum inversion method; the spot becomes blurred and scattered, and the focusing ability becomes worse, which indicates the phase of the Gaussian beam fluctuates. With the increase of , the intensity of atmospheric turbulence continues to increase, the phase fluctuation becomes more and more severe, the coherence of the beam is destroyed, and the light spot may be split into several pieces. Compared with the initial Gaussian beam, the intensity of the beam attenuates significantly, and the larger the , the greater the attenuation of the light intensity.

(a)

(b)

(c)

(d)

(e)

(f)

5. The Effect of Turbulent Atmosphere on Free-Space Optical Communications

In Free-Space Optical Communications (FSO) [16, 17], the laser from the transmitting end can be treated as a plane wave after being collimated by an optical lens [18, 19]; for the plane wave,

From this, the relationship between bit error rate and light intensity fluctuation can be obtained [20, 21]:

Under different fluctuation conditions, the variation of logarithmic amplitude fluctuation and logarithmic intensity fluctuation with transmission distance are shown in Figures 5 and 6, and the relationship between bit error rate and logarithmic light intensity fluctuation variance is shown in Figure 7. In addition, the influence of transmission distance on bit error rate is also given, as shown in Figure Figure 8.

Figures 5 and 6 show the variation of logarithmic amplitude fluctuation and logarithmic intensity fluctuation with transmission distance under different intensity turbulence conditions. It can be clearly seen from the figures that when the transmission distance is constant, the stronger the turbulence intensity, the greater the fluctuation of logarithmic amplitude and logarithmic intensity.

In Free-Space Optical communication (FSO) system with a rate below 2.5 Gbps, the system’s bit error rate is generally required to be below 10⁻⁹. It can be seen from Figure 7 that, under weak fluctuation conditions, for the atmospheric laser communication system’s bit error rate less than the requirement, the light intensity fluctuation should be less than 0.67. It can be seen from Figures 7 and 8 that when the transmission distance is constant, with the increase of turbulence intensity, the bit error rate rises quickly. For the bit error rate 10⁻⁹, when is equal to , the effective distance of communication is close to 4.6 km. With the increase of , the effective distance of transmission is getting smaller and smaller; when it reaches , the effective distance of transmission reduced to about 0.37 km.

6. Conclusion

By comparing the phase screens under atmospheric turbulence of different intensities, it is found that atmospheric turbulence of different intensities affects the propagation of Gaussian beams in terms of light intensity and phase fluctuations. The greater the intensity of the atmospheric turbulence is, the greater the light intensity of the Gaussian beam is affected, the more obvious the fluctuation of the light intensity, and the more obvious the tendency of the light spot to split into several pieces. In terms of phase fluctuations, the greater the intensity of the turbulence, the more intense the phase fluctuations of the Gaussian beam, the larger the difference of the light spot in the figure, and the more obvious the chromatic aberration, which puts forward higher and more specific requirements in the accuracy and reliability of laser transmission under atmospheric turbulence of different intensities. With the increase of turbulence intensity, the bit error rate rises quickly. With the increase of C_n², the effective distance of transmission is getting smaller and smaller; when it reaches , the effective distance of transmission reduced to about 0.37 km.

Data Availability

All datasets generated for this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Key R&D Program-Key Special Project of Inter-governmental S&T Innovation Cooperation (2018YFE0199200), Key R & D Program of Shaanxi Province (2019GY-081), and Key projects of Shaanxi Provincial Department of Education (20JY030).

References

R. Frehlich, “Simulation of laser propagation in a turbulent atmosphere,” Applied Optics, vol. 39, no. 3, pp. 393–397, 2000.
View at: Publisher Site | Google Scholar
H. Liu, Y. Pan, J. Cao, H. Wang, and Y. Zhou, “Adaptive neural network backstepping control of fractional-order nonlinear systems with actuator faults,” IEEE Transactions on Neural Networks and Learning Systems, vol. 99, pp. 1–12, 2020.
View at: Publisher Site | Google Scholar
S. A. Ponomarenko, J. J. Greffet, and E. Wolf, “The diffusion of partially coherent beams in turbulent media,” Optics Communications, vol. 208, no. 1, 2002.
View at: Publisher Site | Google Scholar
Y. Yang, W. Li, and H. Zhang, “A seven-parameter BRDF model with double-peak characteristic suitable for sandy soil,” Mathematical Problems in Engineering, vol. 10, no. 3, Article ID 9398608, 8 pages, 2018.
View at: Publisher Site | Google Scholar
J. A. Fleck, J. R. Morris, and M. D. Feit, “Time-dependent propagation of high energy laser beams through the atmosphere,” Applied Physics, vol. 10, no. 2, pp. 129–160, 1976.
View at: Publisher Site | Google Scholar
B. L. McGlamery, “Restoration of turbulence-degraded images ,” Journal of the Optical Society of America, vol. 57, no. 3, pp. 293–297, 1967.
View at: Publisher Site | Google Scholar
R. J. Noll, “Zernike polynomials and atmospheric turbulence ,” Journal of the Optical Society of America, vol. 66, no. 3, pp. 207–211, 1976.
View at: Publisher Site | Google Scholar
E. P. Walner, “Optimal wave-front correction using slope measurements,” Journal of the Optical Society of America, vol. 73, no. 12, pp. 1771–1776, 1983.
View at: Google Scholar
K. W. Bowman and W. T. Rhodes, “Application of wavelets to wavefront reconstruction in adaptive optical systems,” Adaptive Optics and Applications. International Society for Optics and Photonics, vol. 3126, pp. 288–299, 1997.
View at: Publisher Site | Google Scholar
L. C. Andrews and R. L. Phillips, Laserbeam Propagation Though Random Media, SPIE Optical Engineering Press, Bellingham, WA, USA, 2001.
R. L. Fante, “Electromagnetic beam propagation in turbulent media,” Proceedings of the IEEE, vol. 63, no. 12, pp. 1669–1692, 1975.
View at: Publisher Site | Google Scholar
Y. Jiang, J. Ma, L. Tan, S. Yu, and W. Du, “Measurement of optical intensity fluctuation over an 11.8 km turbulent path,” Optics Express, vol. 16, no. 10, 2008.
View at: Publisher Site | Google Scholar
V. P. Aksenov, V. V. Dudorov, V. V. Kolosov, and V. Y. Venediktov, “Probability distribution of intensity fluctuations of arbitrary-type laser beams in the turbulent atmosphere,” Optics Express, vol. 27, no. 17, 2019.
View at: Publisher Site | Google Scholar
R. Z. Rao, Light Propagation in the Turbulent Atmosphere, Anhui Science and Technology Press, Hefei, China, 2005.
H. M. Zhang, The Research of Simulation of Laser Transmit in the Atmosphere Turbulence, Institute of Optics and Electronics of Chinese Academy of Sciences, Chengdu, China, 2005.
D. Chen, Ke Xi-zheng, and Lu Zhang, “Laser intermodulation distortion and characteristic under the turbulence channel,” Acta Photonica Sinica, vol. 45, no. 2, Article ID 0206007, pp. 1–5, 2016.
View at: Publisher Site | Google Scholar
Y. Baykal, “Signal-to-noise ratio reduction due to oceanic turbulence in oceanic wireless optical communication links,” Optics Communications, vol. 427, pp. 44–47, 2018.
View at: Publisher Site | Google Scholar
B. Dima, “Free-space optical channel simulator for weak-turbulence conditions,” Applied Optics, vol. 31, 2015.
View at: Google Scholar
Y. Xiang and M. Yao, “Free-space communications over exponentiated Weibull turbulence channels with nonzero boresight pointing errors,” Optics Express, vol. 23, no. 3, 2015.
View at: Publisher Site | Google Scholar
J. C. Ricklin, S. M. Hammel, D. Frank, and S. Lachinova, “Atmospheric channel effects on free-space laser communication,” Journal of Optical and Fiber Communications Reports, vol. 3, no. 2, 2006.
View at: Publisher Site | Google Scholar
A. K. Majumdar, “Free-space laser communication performance in the atmospheric channel,” Journal of Optical and Fiber Communications Reports, vol. 2, no. 4, 2005.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Zhi-Qiang Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies