|
1,IntroductionLasers play an important role in the fields of optical communications, lidar, directional damage, aerospace and other fields, and have great scientific and practical value in environmental media measurement and target detection [1]. However, due to the scattering and absorption of the droplets, the laser energy will continue to decay, obviously decay and result in the effectiveness weakness of the laser system. Therefore, it is of great value to study the laser transmission in the fog. Wei Hailiang used the fog droplet spectral index distribution model to study the attenuation characteristics of laser in fog[2]. Li Wei compared the applicability of the fog drop error model in sea fog of different intensities [3], and concluded that the greater the fog density, the closer the model value is to the measured value. Leng Kun studied the influence of propagation distance, visibility and initial RMS radius of the laser on the spatial distribution of photons reaching the receiving plane[4]. Rao Ruizhong mainly introduced the optical properties of the atmosphere, atmospheric refraction, molecular absorption and scattering, light scattering of aerosol particles, etc. In the second chapter of the book, he gives the composition, geometric structure and spatial distribution of droplet particles [5]. Dong Feibiao based on the Mie scattering theory, applied the Monte Carlo multiple scattering model to numerically simulate the transmission characteristics of laser pulses in the fog [6]. Yuan Hui calculated the extinction coefficient of 1.064 μm laser radiation under different visibility conditions using experimental data and empirical formulas [7]. Han Yong calculated the atmospheric extinction coefficient at a wavelength of 0.55 μm through visibility [7] [8], and then obtained the atmospheric extinction coefficient at a wavelength of 1.064 μm. Research on the impact of the above documents on laser transmission includes theoretical research, experimental research and numerical simulation research. Theoretical research is based on the theory of laser propagation attenuation in the atmosphere. In many cases, analytical expressions cannot be obtained; experimental research requires real laser systems and field test conditions, and the work is complicated and costly [9]. Therefore, as a numerical simulation method, Monte Carlo method can simulate all different situations, avoid a lot of experiments, and reflect the real situation of laser transmission to a certain extent. Based on the Monte Carlo method, this paper quantitatively analyzes the spatial distribution of photons on the receiving plane and the temporal distribution of photons on the receiving plane, and compares the calculation results with the MODTRAN online calculation results. 2,Theoretical Analysis2.1,Analysis of laser attenuation characteristics in fogGenerally, the scattering generated when the size of the droplet particles is similar to the laser wavelength, which AOPC 2020: Advanced Laser Technology and Application, edited by Zhiyi Wei, Jing Ma, Wei Shi, Xuechun Lin, can be solved by Mie theory. Due to the arbitrariness and randomness of the orientation of the droplet particles in space and the small particle size, it is statistically regarded as an isotropic spherical particle for analysis [2]. According to the Mie scattering theory, the attenuation cross section σt and scattering cross section of a single droplet particle σs are defined as follows [5]: Where an and bn are Mie coefficients, k = 2π/λ,n is the eigenvalue. At present, the commonly used droplet Gamma particle size distribution model is [6]: Where N is the number density of droplet particles and r is the radius of droplet particles. r is in meters. Where V is the visibility of fog (km) and W is the water content (g/m3). To calculate the attenuation coefficient of fog droplets for laser transmission, the relationship between visibility and attenuation coefficient can be approximated [8]. Where μt0.55 is the attenuation coefficient of the most sensitive wavelength of the human eye. The attenuation coefficient and scattering coefficient of fog droplets to laser transmission can also be obtained by using the attenuation cross section and scattering cross section: Where n(r) is the probability density function of the droplet particle size distribution [5], Ni is the total number of particles of component i per unit volume, and the two characteristic parameters rmodN,i and σi are the average radius and standard deviation, respectively. The asymmetry factor g of the particles indicates the degree of asymmetry of the front and back scattering. The calculation formula is as follows [11]: Where, Qs is the scattering efficiency factor, and ω = μs/μt is the single scattering efficiency. Table 1 shows the calculation results of the extinction parameters of 1.064 μm laser with different droplet particles using the above formula, where the visibility is V = 100m and the complex refractive index of the droplet on the infrared laser is n=1.327-2.89*10-6i [6], and the scale parameter of different droplet particles can be found in[10]. Table 1the comparison of the calculation results of the extinction parameters in different fogs
Through theoretical analysis of droplet particles, it is found that the attenuation of fog on a laser with a wavelength of 1.064 μm is closely related to the type of fog. Different types of droplets have different particle radius and different attenuation coefficients for light. At the same time, when the visibility is V = 1000m, the attenuation coefficient of advection fog is 0.002, the single scattering rate is 0.9999, and the asymmetry factor is 0.7214. 2.2,Monte Carlo calculation model2.2.1,Initial beam analysisAccording to the theory of laser resonator cavity, the self-reproducing mode generated by the stabilized cavity follows the gaussian distribution, and the output of the optical cavity is the gaussian beam. In the plane perpendicular to the transmission direction of the laser beam, its light intensity distribution is [12]. Where, γS is the RMS radius of the laser beam, which is used to describe the radial distribution of the beam, x and y are the abscissa and ordinate on the plane perpendicular to the transmission direction of the laser beam, respectively. In the source sampling of the Monte Carlo method, a Gaussian distribution sampling model is introduced to make the statistical characteristics of the emitted photons satisfy the characteristics of the Gaussian beam. 2.2.2Analysis of photon transport in fogThe basic idea of the Monte Carlo method to simulate photon transmission is: photons enter the random medium with a certain direction and weight and initial time, then the next transmission step is given after being scattered by the droplet particles, and the scattering phase function is sampled to obtain new direction of motion after scattering. This process is repeated until the photon is received by the target screen or the weight is less than the threshold (10-6 in this paper)[6]. The propagation process is shown in Figure 1 [4]. The motion state of the photon is represented by the spatial position (x,y,z), the moving direction (θ, φ), and the photon weight W. θ and φ are the scattering angle and azimuth angle, respectively. The geometric path length of each step of the photon is as follows. Where εr1 is a random number evenly distributed between [0,1]. After the step l is determined, the next collision position (x’, y’, z’) of the photon is determined by the current photon position (x,y,z) and the direction cosine (μx,μy,μz). Then determine whether the photon is scattered. Taking a random number εr2 evenly distributed between [0,1]. If εr2 < ω, the photon is scattered; otherwise, the photon is absorbed and the tracking of the photon is terminated [13]. Assuming that the initial weight of the photon is W = 1, after scattering, the weight of the photon becomes W = W·ω. After the photon is scattered, the new propagation direction is determined by the scattering phase function, and the Henyey-Greenstein function can be used to approximate its scattering phase probability function; sometimes in order to obtain a higher accuracy of the scattering phase probability function, the weighted sampling method can also be used find the phase function [14]. In this paper, the Henyey-Greenstein function is selected in article [15], which has the following form. Where g is the asymmetry factor. Assuming that the scattering azimuth φ is evenly distributed in [0, 2π], the sampled value can be expressed as: Where, εr3 is a random number evenly distributed between [0,1]. According to the scattering direction of the photon, the new direction cosine in the global coordinate system after the photon collision can be obtained by coordinate transformation as 3.Simulation results and analysisAs shown in Figure 1, after the laser been transmitted for a distance L, a receiving plane is set perpendicular to the propagation direction to study the spatial distribution of photons. The total number of simulated photons is 2000000, The divergence angle β is 0.001rad. The attenuation coefficient μt, the single scattering rate ω), and the asymmetry factor g are given in [6]. The influence of the propagation distance on the spatial distribution of photons is firstly considered in the receiving plane. The Monte Carlo simulation results are shown in Figure 2, where the visibility V = 100m and the initial beam RMS radius γs =0.5m. Then, the influence of visibility on the spatial distribution of photons in the receiving plane is considered. The Monte Carlo simulation results are shown in Figure 3, where the propagation distance L = 10m and the initial beam root-mean-square radius γs =0.5m. Finally, we consider the effect of the initial RMS radius on the spatial distribution of the photons in the receiving plane is considered. The Monte Carlo simulation results are shown in Figure 4, where the propagation distance L = 20m and visibility V = 100m. It can be seen from Figure 2 that with the increase of the propagation distance, the spatial distribution of photons at the receiving plane has not significantly deviated from the Gaussian distribution, but the number of photons reaching the receiving plane has decreased significantly. This is because the farther the propagation distance, the greater the probability of photons being scattered, resulting in a reduction in the number of photons received. It can be seen from Figure 3 that the number of photons at the receiving plane decreases first and then increases with the increase of visibility, but the spatial distribution of photons has not deviated significantly from the Gaussian distribution. This is because the density of droplet particles decreases as the visibility increases, the probability of photons being scattered decreases, and the number of photons reaching the receiving plane increases. It can be seen from Figure 4 that the larger the initial RMS radius of the laser, the more dispersed the spatial distribution of photons, but it has not deviated significantly from the Gaussian distribution. Figures 2 and 4 are consistent with the conclusions of reference [4]. When the visibility is less than 20m in Figure 3, the law is slightly different from reference [4], showing the randomness of photon propagation in dense fog. To further analyze the quantitative relationship between the spatial distribution of photons and various parameters, Figure 5 shows the correspondence between the number of photons reaching the receiving plane and the propagation distance, visibility, and initial RMS radius. As shown in Figure 5, as the propagation distance increases, the number of photons arriving at the receiving plane decreases approximately linearly; as the visibility increases, the number of photons reaching the receiving plane decreases slightly at first, and then increases rapidly; As the root radius increases, the number of photons reaching the receiving plane increases slowly, but the change is not obvious and is roughly linear. Figure 6 shows the normalized photon time distribution received at different distances. It can be seen from Figure 6 that for a certain propagation distance, the number of normalized photons received at the receiving plane increases rapidly with time and then decreases slowly; while with the increase of the propagation distance, the peak time of receiving photons at the receiving plane gradually increases. In order to further understand the relationship between the photon peak time at the receiving plane and the propagation distance, we extracted the relevant data, and the results are shown in Figure 7. Figure 7 shows the relationship between the peak time of the photon received by the receiving plane and the propagation distance. It can be seen that this relationship basically shows a linear change rule. In addition, the parameter visibility V = 1000m, the propagation distance L = 1000m, and the attenuation coefficient μt1.064 = 0.002 are substituted into the atmospheric transmittance formula T = exp(–μ · L) to calculate the atmospheric transmittance of 0.1353. Take the urban surface albedo 0.15 [16], substitute MODTRAN online calculation results to get T = 0.1253, the relevant conclusions of the two are basically the same. 4.ConclusionThe numerical simulation by Monte Carlo method shows that as the propagation distance increasing, the number of photons reaching the receiving plane decreases approximately linearly, but the spatial distribution does not deviate significantly from the Gaussian distribution. It begins to decrease slowly at first and then increases, but the spatial distribution has not significantly deviated from the Gaussian distribution; with the increase of the initial RMS radius, the number of photons reaching the receiving plane does not change initially and then increases slowly. In addition, for a certain propagation distance, the number of normalized photons received at the receiving plane increases rapidly with time and then decreases slowly; while with the increasement of the propagation distance, the peak time of receiving photons at the receiving plane gradually increases. REFERENCESChen, X. X., Wu, Z. S.,
“Study on the Scattering Polarization Characteristics of Laser in Low Visibility Atmosphere,”
Xidian University. Papers, 1
–98
(2019). Google Scholar
Wei, H. L., Shao, L. M., Li, M. L.,
“Research on propagation attenuation characteristics of laser fog based on Mie theory,”
Ship Electronic Engineering. Papers, 34
(1), 148
–152
(2014). Google Scholar
Li, W., Shao, L. M., Yuan, Q. Z.,
“Experimental Study on the Attenuation of Yellow Sea Advection Fog to 8-12μm Infrared Radiation,”
Infrared And Laser Engineering. Papers, 48
(7), 0704005
(2019). https://doi.org/10.3788/IRLA Google Scholar
Leng, K., Zhang, X.,
“Monte Carlo simulation of laser propagation in marine aerosol,”
Physics and Engineering. Papers, 28
(6), 74
–79
(2018). Google Scholar
Rao, R. Z., [Modern atmospheric optics], 152
–154 Science Press, Beijing
(2012). Google Scholar
Dong, F. B., Lin, W. B.,
“Monte Carlo Simulation of 1.064μm Pulse Laser Transmission in Fog,”
Infrared. Papers, 35
(9), 37
–42
(2014). Google Scholar
Yuan, H.,
“Attenuation characteristics of 1.06μm laser radiation in water mist,”
Infrared And Laser Engineering. Papers, 47
(10), 181
–187
(2018). Google Scholar
Han, Y., Rao, R. Z., Wang, Y. J.,
“Obtaining multi-wavelength aerosol optical properties using atmospheric visibility,”
Infrared And Laser Engineering. Papers, 36
(2), 267
–269
(2007). Google Scholar
Liu, J. R.,
“[High Energy Laser System Test and Evaluation],”
36
–38 National Defense Industry Press, Beijing
(2014). Google Scholar
Liu, R. Y.,
“Observational Study on Macro and Micro Physical Structure of Fog in Winter in Nanjing,”
Nanjing University of Information Science and Technology. Papers, 47
–49
(2011). Google Scholar
Wang, H. X., Zhu, Y. Z.,
“Study on the transmission characteristics of laser in different types of aerosol,”
Acta Phys. Papers, 62
(2), 1
–10
(2013). Google Scholar
Zhang, X.,
“Influence of fluctuating water on laser underwater target detection based on Monte Carlo method,”
CHINESE JOURNAL OF LASERS. Papers, 39
(7), 0702003
(2012). https://doi.org/10.3788/CJL Google Scholar
Ramella-Roman, J. C., Prahl, S. A., Jacques, S. L.,
“Three Monte Carlo programs of polarized light transport into scattering media,”
Opt. Express., 13
(4420),
(2005). Google Scholar
Hu, S., Gao, T. Z., Li, H.,
“Near infrared laser transmission characteristics based on phase function piecewise weighted sampling,”
CHINESE JOURNAL OF LASERS. Papers, 42
(8), 0802008
(2015). https://doi.org/10.3788/CJL Google Scholar
Toublanc, D.,
“Henyey–Greenstein and Mie phase functions in Monte Carlo radiative transfer computations,”
Applied optics. Papers, 35
(18), 3270
–3274
(1996). https://doi.org/10.1364/AO.35.003270 Google Scholar
Liu, Q. Q.,
“Study on Spectral Radiation Albedo of Different Land Use Types in China,”
REMOTE SENSING TECHNOLOGY AND APPLICATION. Papers, 34
(1), 50
–56
(2019). Google Scholar
|