-
The spatial power spectrum of the refractive-index fluctuations for the oceanic turbulence can be described as[8]
$$\begin{array}{l} {\varPhi _n}(\kappa ) = 0.388 \times {10^{ - 8}}{\varepsilon ^{ - 1/3}}{\kappa ^{ - 11/3}}[1 + 2.35{(\kappa \eta )^{2/3}}] \times \\ \dfrac{{{\chi _T}}}{{{\omega ^2}}}[{\omega ^2}\exp( - {A_T}{\delta _{TS}}) + \exp ( - {A_S}{\delta _{TS}}) - 2\omega \exp ( - {A_{TS}}{\delta _{TS}})] \\ \end{array} $$ (1) Here,
$\varepsilon $ is the rate of dissipation of kinetic energy per unit mass of fluid ranging from 10–1 − 10–10 m2s–3;${\chi _T}$ is the rate of dissipation of mean-squared temperature with the range${10^{ - 4}}-{10^{ - 10}}{{\rm{K}}^{\rm{2}}}{{\rm{s}}^{{\rm{ - 1}}}}$ ;$\omega $ is the ratio of temperature and salinity contributions to the refractive index spectrum varying in –5-0 corresponding to dominating temperature-induced and salinity-induced turbulence, respectively;$\eta $ is the Kolmogorov micro scale.${A_{TS}} = 9.41 \times {10^{ - 3}}$ ,${A_S} = 1.9 \times {10^{ - 4}} $ ,${A_T} = 1.863 \times {10^{ - 2}}$ ,${\delta _{TS}} = 8.284{\left( {\kappa \eta } \right)^{4/3}} + 12.987{\left( {\kappa \eta } \right)^2}$ .Ref.[5] introduces the basic expressions of Gaussian beam propagate through random media. The long-term spot radius
${W_{{{LE}}}}$ through turbulence is defined as$${W_{{{LE}}}} = W\sqrt {1 + T} $$ (2) $$\begin{split} T &= 4{\pi ^2}{k^2}L\int_0^1 {d\xi } \int_0^\infty {d\kappa \kappa {\varPhi _n}\left( \kappa \right)}\times \\ &\left[ {1 - \exp \left( { - \dfrac{{\varLambda L{\kappa ^2}{\xi ^2}}}{k}} \right)} \right] \\ \end{split} $$ (3) Here,
$W = {{{w_0}} / {\sqrt {{\Theta ^2} + {\varLambda ^2}} }}$ , w0 is the initial spot size. Based on the general model, beam wander through turbulence can be expressed as$$\begin{split} \left\langle {r_c^{\rm{2}}} \right\rangle & = 4{\pi ^2}{k^2}L\int_0^1 {d\xi } \int_0^\infty {d\kappa \kappa {\varPhi _n}\left( \kappa \right)} \\ &{H_{LS}}\left[ {1 - \exp \left( { - \dfrac{{\varLambda L{\kappa ^2}{\xi ^2}}}{k}} \right)} \right] \end{split} $$ (4) Large-scale filter function HLS is
$${H_{LS}} = \exp \left\{ { - {\kappa ^2}W_0^2\left[ {{{\left( {{\Theta _0} + {\Theta _0}\xi } \right)}^2} + \varLambda _0^2{{\left( {1 - \xi } \right)}^2}} \right]} \right\}$$ (5) The on-axis scintillation of Gaussian beam through turbulence is
$$\begin{split} \sigma _I^2\left( {0,L} \right)& = 8{\pi ^2}{k^2}L\int_0^1 {d\xi } \int_0^\infty {d\kappa \kappa {\varPhi _n}\left( \kappa \right)} \exp \left( { - \dfrac{{\varLambda L{\kappa ^2}{\xi ^2}}}{k}} \right) \times \\ & \left\{ {1 - \cos \left[ {\frac{{L{\kappa ^2}}}{k}\xi \left( {1 - \left( {1 - \Theta } \right)\xi } \right)} \right]} \right\} \end{split} $$ (6) Following the approach derived in Ref.[7], Λ0 and Θ0 are defined to characterize the unperturbed beams in terms of the initial phase front radii of curvature F0 and spot size w0:
$${\varLambda _0} = \dfrac{{2L}}{{kw_0^2}},{\Theta _0} = 1 - \dfrac{L}{{{F_0}}}$$ (7) Collimated Gaussian beam with
${F_0} = \infty $ is considered in this paper. Also, k represents the wavenumber, and L is the propagation path length. Λ and Θ can also be used to simplify the mathematical expressions$$\varLambda = \dfrac{{{\varLambda _0}}}{{\Theta _0^2 + \varLambda _0^2}},\Theta = \dfrac{{{\Theta _0}}}{{\Theta _0^2 + \varLambda _0^2}}$$ (8) Since oceanic turbulence could achieve strong fluctuation in a short distance, the theory formula should apply to from weak to strong fluctuation conditions. The method of effective beam parameters maybe used to extend these following formulas into the strong regime by replacing Λ and Θ by their effective beam parameters
$${\varLambda _e} = \dfrac{\varLambda }{{1 + 4q\varLambda /3}},{\Theta _e} = \dfrac{{\Theta - 2q\varLambda /3}}{{1 + 4q\varLambda /3}}$$ (9) Substituting Eq.(1) into Eq.(2)−(6), the long-exposure beam radius
${W_{{{LE}}}}$ , beam centroid displacement${\sigma _c}$ and on-axis scintillation index$\sigma _I^2$ for Gaussian beam propagating through oceanic turbulence can be derived. -
The principle of phase screen method [7]is considering the beam propagate through vacuum and phase screens along the transmission path alternately as shown inFig.1. The propagation distance L is divided into many parts, turbulence effects of each part is considered to be a phase screen, the light field between two adjacent phase screen can be expressed as
$$\begin{split} u({{r}},{z_{j + 1}})& = {F^{ - 1}}\{ F\left[ {u({{r}},{z_j})\exp \left[ {{\rm{i}}\varphi (x,y)} \right]} \right]\times \\ & \exp ( - {\rm{i}}\dfrac{{\kappa _x^2 + \kappa _y^2}}{{2k}}\Delta {z_{j + 1}})\} \\ \end{split} $$ (10) Here,
$\Delta {z_{j + 1}} = {z_{j + 1}} - {z_j}$ is the distance between two phase screen located at${z_{j + 1}}$ and${z_j}$ ;${\kappa _x}$ and${\kappa _y}$ denotes the spatial frequency, the grids are divided to be$N \times N$ and the resolution is$\Delta x$ ,$F$ and${F^{ - 1}}$ represents Fourier transform and Fourier inverse transform, respectively;$\varphi \left( {x,y} \right)$ reveals the phase fluctuation caused by turbulence which can be expressed as [7]$$\begin{split} \varphi \left( {x,\;y} \right) &= C\sum\limits_{{K_x}} {\sum\limits_{{K_Y}} {a\left( {{\kappa _x}{\rm{,}}\;{\kappa _y}} \right)\sqrt {{\varPhi _\theta }\left( {{\kappa _x}{\rm{,}}\;{\kappa _y}} \right)} } }\times \\ & \exp \left[ {i\left( {{\kappa _x}x + {\kappa _y}y} \right)} \right] \end{split} $$ (11) Here,
$x = m\Delta x$ ,$y = m\Delta y$ represents the spatial domain,$\Delta x$ ,$\Delta y$ is the sample interval, m, n are integers;${\kappa _x} = {m^{'}}\Delta {\kappa _x}$ ,${\kappa _y} = {n^{'}}\Delta {\kappa _y}$ represents the frequency domain,$\Delta {\kappa _x}$ ,$\Delta {\kappa _y}$ is the sample interval,$m'$ ,$n'$ are integers;$a\left( {{\kappa _x},{\kappa _y}} \right)$ is the Fourier transform of the Gaussian random matrix with the mean value 0 and variance 1;${\varPhi _\theta }\left( {{\kappa _x}{\rm{,}}\;{\kappa _y}} \right)$ represents the phase power spectrum, which is related to the power spectrum as$${\varPhi _\theta }({\kappa _x},{\kappa _y}) = 2\pi {k^2}\int\limits_z^{\Delta z} {{\varPhi _n}({\kappa _x},{\kappa _y},\xi )} {\rm{d}}\xi $$ (12) -
The parameter of phase screen mainly include the grid size
$\Delta x$ , screen size x and the number of phase screen along the propagation path NPS. Since FFT method is usually used to calculate the beam propagation for program, the sample rule request is derived by Ref.[9]$$\Delta x = \dfrac{{\lambda L}}{x}$$ (13) Here,
$\lambda $ is wavelength;L is the propagation distance. Considering that$x = \sqrt {\lambda LN} $ ,$\Delta x = \sqrt {\lambda L/N} $ , the relationships between grid size$\Delta x$ , screen size x and grids number N is decided. Usually the grids number N is decided firstly, and with higher N, the resolution of sample is more approximate while the computation complexity is higher. In this paper, the influence of different N is shown and compared.Meanwhile, the description of beam should also be considered. For Gaussian beam, the beam width
${w_0}$ should be much bigger than$\Delta x$ . Also, consider about the characteristics of turbulence, the sampling should be less than half of the coherence length$\Delta x \leqslant {r_0}/2$ . Here the oceanic turbulence coherence length${r_0} \approx 2.1{\rho _0}$ [10], spatial coherence scale${\rho _0}$ for oceanic turbulence is derived by Ref.[11].$$\begin{split} &{\rho _0} = {[3.603 \times {10^{ - 7}}{k^2}\Delta z{\varepsilon ^{ - 1/3}}\dfrac{{{\chi _T}}}{{2{\omega ^2}}}(16.958{\omega ^2} - 44.175\omega + 118.923)]^{ - 1/2}}\;\;\;\;({\rho _0} < < \eta ) \\ &{\rho _0} = {[3.603 \times {10^{ - 7}}{k^2}\Delta z{\varepsilon ^{ - 1/3}}\dfrac{{{\chi _T}}}{{2{\omega ^2}}}(1.116{\omega ^2} - 2.235\omega + 1.119)]^{ - 3/5}}\;\;\;\;({\rho _0} > > \eta ) \\ \end{split} $$ (14) The number of phase screen along the path is decided by the turbulence strength. According to Ref.[12], the fluctuation between two phase screens should be small enough and meet the requirements of Rytov variance
$\sigma _R^2\left( {\Delta z} \right) < 0.1$ , then the requirement of phase screen number could be expressed as$${N_{PS}} > {\left[ {10\sigma _R^2\left( L \right)} \right]^{6/11}}$$ (15) Rytov variance
$\sigma _R^2$ can be approximately expressed as[13]:$$\begin{split} \sigma _R^2 &= 3.063 \times {10^{ - 7}}{k^{7/6}}{L^{11/6}}{\varepsilon ^{ - 1/3}}{\chi _T}× \\ &\left( {0.358{\omega ^2} - 0.725\omega + 0.367} \right)/{\omega ^2} \end{split} $$ (16) The example of parameter settings is as follows:
(1) Consider about the strongest turbulence condition as an example, L=50 m,
$\varepsilon = {10^{ - 6}}{{\rm{m}}^2}{{\rm{s}}^{ - 3}}$ ,${\chi _T} = {10^{ - 6}}{K^2}{s^{ - 1}}$ ,$\omega = - 2$ . Substituting into Eq.(14),${r_0} \approx 0.7\;{\rm{mm}}$ is obtained. Consider that$\Delta x \leqslant {r_0}/2$ ,$\Delta x < < {w_0}$ , the grid size should fulfill$\Delta x \leqslant 0.35\;{\rm{mm}}$ , grids number should fulfill$N \geqslant 238$ .(2) Use Eq.(15) to calculate the phase screen number, and the distance between adjacent phase screens are should meet
$\Delta z \leqslant 2.5$ . During the simulation in this paper, the distance between screens is 2.5 m, which is much lower than the atmospheric turbulence simulation conditions. -
Assuming that the long exposure intensity distribution is Gaussian, the statistical properties calculation could be expressed as [7]. Long exposure beam footprint radius:
$$W_{{{LE}}}^2 = {{2\int\limits_{ - \infty }^\infty {{r^2}\left\langle {I\left( {{{r}},L} \right)} \right\rangle {\rm{d}}{{r}}} } / {\int\limits_{ - \infty }^\infty {\left\langle {I\left( {{{r}},L} \right)} \right\rangle {\rm{d}}{{r}}} }}$$ (17) Standard deviation of the fluctuation of the instantaneous center of the beam:
$${\sigma _c} = \sqrt {\left\langle {{{r}}_c^{\rm{2}}} \right\rangle } $$ (18) $${{{r}}_c} = {{2\int\limits_{ - \infty }^\infty {{{r}}I\left( {{{r}},L} \right){\rm{d}}{{r}}} } / {\int\limits_{ - \infty }^\infty {\left\langle {I\left( {{{r}},L} \right)} \right\rangle {\rm{d}}{{r}}} }}$$ (19) On-axis scintillation index:
$$\sigma _I^2 = {{\left\langle {{{\left[ {I(0,L) - \left\langle {I(0,L)} \right\rangle } \right]}^2}} \right\rangle }/ {{{\left\langle {I(0,L)} \right\rangle }^2}}}$$ (20) Here, < > represents the average and
${{{r}}_c}$ is the position of beam center.
Validity of beam propagation characteristics through oceanic turbulence simulated by phase screen method
-
摘要: 海洋湍流是制约水下激光通信应用的重要因素之一,相对于大气湍流,海洋湍流在较短距离内即可达到强起伏,针对强湍流传输特性及数值仿真的可靠性验证具有重要意义。基于大气湍流传输特性理论公式,利用大气湍流参数等效表达海洋湍流的方法给出了弱起伏到强起伏条件下海洋湍流传输特性的理论计算公式。针对相位屏法在海洋湍流仿真中的应用,给出了选取相位屏间距、网格尺寸和网格数目的基本要求,并数值仿真了海洋湍流不同参数条件下的传输特性,与理论计算结果进行了对比。结果表明:基于相位屏法得到的光强一阶矩传输特性参量与解析结果较为一致,但是在强闪烁条件下数值仿真光强闪烁特性与理论结果偏差较大。Abstract: Oceanic turbulence is an important factor to restrict the application of underwater optical communication. Phase screen method is a simple and effective way to simulate the propagation process of complex beams through turbulence. The constraints of parameter setting for phase screen simulated oceanic turbulence based on the sampling principle and turbulence effects were firstly discussed here. Furthermore, the theoretical expressions of propagation characteristics of Gaussian beam through oceanic turbulence from weak to strong fluctuation regime were derived. Our goal in this research was to testify the validity of phase screen method in oceanic turbulence by comparison of major statistical characteristics of Gaussian beam propagating in oceanic turbulence simulated by phase screen method and the theoretical expressions derived. Results show good match between simulation results and theory formulas for long exposure beam radius and centroid displacement under different turbulence conditions, as well as the scintillation index under weak fluctuation regime. However, results show significant mismatch between numerically estimated and theoretically predicted values for the on-axis scintillation index in strong fluctuation regime.
-
Key words:
- oceanic turbulence /
- phase screen method /
- beam wander /
- scintillation index
-
-
[1] Li Y, Yu L, Zhang Y. Influence of anisotropic turbulence on the orbital angular momentum modes of Hermite-Gaussian vortex beam in the ocean [J]. Optics Express, 2017, 25(11): 12203. doi: 10.1364/OE.25.012203 [2] Li D, Deng D, Zhao F, et al. Propagation properties of spatiotemporal chirped Airy Gaussian vortex wave packets in a quadratic index medium [J]. Optics Express, 2017, 25(12): 13527. doi: 10.1364/OE.25.013527 [3] Xu Chenlu,Hao Shiqi,Zhang Dai, et al. Design of atmospheric turbulence phase screen set under the influence of combined oblique propagation and beam propagation [J]. Infrared and Laser Engineering, 2019, 48(4): 0404003. (in Chinese) doi: 10.3788/IRLA201948.0404003 [4] Li Dun, Ning Yu, Wu Wuming, et al. Numerical simulation and verification of dynamic atmospheric turbulence with rotating phase screens [J]. Infrared and Laser Engineering, 2017, 46(12): 1211003. (in Chinese) doi: 10.3788/IRLA201746.1211003 [5] Yang T, Zhao S. Study on the Model of Ocean Turbulence Random Phase Screen [J]. Acta Optica Sinica, 2017, 37(12): 1201001. [6] Niu Chaojun, Lu Fang, Han Xiang'e. Propagation properties of Gaussian array beams transmitted in oceanic turbulence simulated by phase screen method [J]. Acta Optica Sinica, 2018, 38(6): 0601004. doi: 10.3788/AOS201838.0601004 [7] Andrews L, Phillips R. Laser Beam Propagation Through Random Media[M]. Bellingham: SPIE Press, 2005. [8] Nikishov V V, Nikishov V I. Spectrum of turbulent fluctuations of the sea-water refraction index [J]. Fluid Mechanics Research, 2000, 27(1): 82−98. doi: 10.1615/InterJFluidMechRes.v27.i1.70 [9] Voelz D, Roggemann M. Digital simulation of scalar optical diffraction: revisiting chirp function sampling criteria and consequences [J]. Applied Optics, 2009, 48(32): 6132−6142. doi: 10.1364/AO.48.006132 [10] Fried D. Optical resolution through a randomly inhomogeneous medium for very long and very short exposures [J]. JOSA A, 1966, 56(10): 1372−1379. doi: 10.1364/JOSA.56.001372 [11] Lu L, Ji X, Baykal Y. Wave structure function and spatial coherence radius of plane and spherical waves propagating through oceanic turbulence [J]. Optics Express, 2014, 22(22): 27112−27122. doi: 10.1364/OE.22.027112 [12] Lachinova S, Vorontsov M. Giant irradiance spikes in laser beam propagation in volume turbulence: analysis and impact [J]. Journal of Optics, 2016, 18(2): 025608. doi: 10.1088/2040-8978/18/2/025608 [13] Lu Lu. Influence of oceanic turbulence on propagation of laser beams[D]. Hefei:University of Science and Technology of China, 2016: 52-56.(in Chinese) [14] Flatte S M, Wang G Y, Martin J. Irradiance variance of optical waves through atmospheric turbulence by numerical simulation and comparison with experiment [J]. JOSA A, 1993, 10(10): 2363−2370. [15] Consortini A, Cochetti F, Churnside J H, et al. Inner-scale effect on irradiance variance measured for weak-to-strong atmospheric scintillation [J]. JOSA A, 1993, 10(11): 2354−2362. doi: 10.1364/JOSAA.10.002354 [16] Andrews L, Phillips R, Ronald L. Strehl ratio and scintillation theory for uplink Gaussian-beam waves: beam wander effects [J]. Optical Engineering, 2006, 45(7): 076001. doi: 10.1117/1.2219470