留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Validity of beam propagation characteristics through oceanic turbulence simulated by phase screen method

Niu Chaojun Wang Xiaobin Lu Fang Han Xiang’e

牛超君, 王晓斌, 卢芳, 韩香娥. 相位屏法仿真海洋湍流激光传输特性有效性[J]. 红外与激光工程, 2020, 49(7): 20190452. doi: 10.3788/IRLA20190452
引用本文: 牛超君, 王晓斌, 卢芳, 韩香娥. 相位屏法仿真海洋湍流激光传输特性有效性[J]. 红外与激光工程, 2020, 49(7): 20190452. doi: 10.3788/IRLA20190452
Niu Chaojun, Wang Xiaobin, Lu Fang, Han Xiang’e. Validity of beam propagation characteristics through oceanic turbulence simulated by phase screen method[J]. Infrared and Laser Engineering, 2020, 49(7): 20190452. doi: 10.3788/IRLA20190452
Citation: Niu Chaojun, Wang Xiaobin, Lu Fang, Han Xiang’e. Validity of beam propagation characteristics through oceanic turbulence simulated by phase screen method[J]. Infrared and Laser Engineering, 2020, 49(7): 20190452. doi: 10.3788/IRLA20190452

相位屏法仿真海洋湍流激光传输特性有效性

doi: 10.3788/IRLA20190452
详细信息
  • 中图分类号: TN249

Validity of beam propagation characteristics through oceanic turbulence simulated by phase screen method

More Information
    Author Bio:

    牛超君(1991-),女,博士生,主要从事大气传输特性及自适应光学方面的研究。Email:niuchaojun123@163.com

    韩香娥(1962-),女,教授,博士生导师,博士,主要从事光传输、激光雷达和散射方面的研究。Email:xd_xehan@163.com

图(5)
计量
  • 文章访问数:  386
  • HTML全文浏览量:  85
  • PDF下载量:  56
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-05
  • 修回日期:  2020-03-20
  • 网络出版日期:  2020-07-23
  • 刊出日期:  2020-07-23

Validity of beam propagation characteristics through oceanic turbulence simulated by phase screen method

doi: 10.3788/IRLA20190452
    作者简介:

    牛超君(1991-),女,博士生,主要从事大气传输特性及自适应光学方面的研究。Email:niuchaojun123@163.com

    韩香娥(1962-),女,教授,博士生导师,博士,主要从事光传输、激光雷达和散射方面的研究。Email:xd_xehan@163.com

  • 中图分类号: TN249

摘要: 海洋湍流是制约水下激光通信应用的重要因素之一,相对于大气湍流,海洋湍流在较短距离内即可达到强起伏,针对强湍流传输特性及数值仿真的可靠性验证具有重要意义。基于大气湍流传输特性理论公式,利用大气湍流参数等效表达海洋湍流的方法给出了弱起伏到强起伏条件下海洋湍流传输特性的理论计算公式。针对相位屏法在海洋湍流仿真中的应用,给出了选取相位屏间距、网格尺寸和网格数目的基本要求,并数值仿真了海洋湍流不同参数条件下的传输特性,与理论计算结果进行了对比。结果表明:基于相位屏法得到的光强一阶矩传输特性参量与解析结果较为一致,但是在强闪烁条件下数值仿真光强闪烁特性与理论结果偏差较大。

English Abstract

    • In recent years, underwater laser communication system has attracted more and more attentions. Oceanic turbulence is one of the main concerns for research of laser beams propagating through the ocean[1- 2]. Phase screen method is a simple and effective way to simulate the propagation process of complex beams through random medium[3-4]. On one hand, the randomness of each realization could provide good reference for experiments. On the other hand, the statistical property of multiple realizations also testify the theory formulas. Recently, phase screen method has been used to simulate oceanic turbulence[5-6]. Since oceanic turbulence has more influence factors and usually much stronger than atmospheric turbulence, it is necessary to discuss the principle of parameter settings for oceanic turbulence simulated by phase screen method and compare the results calculated by theoretical formula and numerical simulation to confirm the accuracy.

      In this paper, the equation of beam spreading, beam wander and scintillation of Gaussian beam propagating through oceanic turbulence are derived follow the theory of beam propagation through random media[7]. According to the phase screen method used for atmospheric turbulence simulation, the rule of parameter setting for phase screen method used for oceanic turbulence simulation are given. At last, the theoretical computation results and simulation results for statistical property of Gaussian beam propagating through oceanic turbulence are compared and results testify the validity range of phase screen simulated oceanic turbulence.

    • 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

      Figure 1.  Diagram of phase screen simulation method

      $$\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.

    • The following values of parameters are used in the calculations unless the statement of other values: wavelength $\lambda = 532\;{\rm{nm}}$ , waist width ${w_0} = 10\;{\rm{mm}}$ , inner scale $\eta = 1\;{\rm{mm}}$ , $L = 10\;{\rm{m}}$ , $\varepsilon = {10^{ - 6}}{{\rm{m}}^2}{{\rm{s}}^{ - 3}}$ , $\omega = - 2$ , ${\chi _T} = {10^{ - 6}}{K^2}{s^{ - 1}}$ . According to the example show in Sect.3.2, the grids of phase screen is N=1024 and 2048, leading to the parameter setting $N = 1\;024,x = 0.16\;{\rm{m}}, \Delta x = 0.16\;{\rm{mm}}$ and $N = 2\;048$ , $x = 0.24\;{\rm{m}}$ , $\Delta x = 0.12\;{\rm{mm}}$ , respectively. The distance between screens is 2.5 m. Note that the theory results are calculated by the expressions shown in Sec.2 and the numerical simulation results is obtained by the expressions in Sec.3.3 after the average of 1000 realizations.

      The results presented in Fig.2(a) and 2(b) reveal that using phase screen method can achieve good correspondence between theoretical and numerical calculations of such statistical characteristics as long-exposure beam radius WLE and ${\sigma _c}$ on the basis of appropriate parameter setting. Higher resolution numerical grid (2048 × 2 048 pixels) leads to better accuracy and longer distance weaken the correspondence. However, Fig.2(c) shows that the scintillation index $\sigma _I^2$ , which represents the statistical characteristics of second moment of intensity, only have good agreement between theory and numerical results in short distance propagation and the increasement of grids does not make much change to the results.

      Figure 2.  Propagation characteristics versus different propagation distance, (a) long-exposure beam radius WLE, (b) beam centroid displacement ${\sigma _c}$ , (c) on-axis scintillation index $\sigma _I^2$

      Fig. 3-5 shows the similar tendency as Fig.2 that the analytical formula and phase screen numerical simulation for beam spreading and beam wander of Gaussian beam through oceanic turbulence are in good agreements under different turbulence conditions. The simulation results tends to be more approximate when the grids resolution gets higher. Figures also reveal that ${W_{{{LE}}}}$ and ${\sigma _c}$ increase when turbulence gets stronger (higher dissipation of mean-squared temperature ${\chi _T}$ , lower rate of dissipation of kinetic energy per unit mass of seawater $\varepsilon $ and higher ratio of temperature to salinity contribution $\omega $ ). However, stronger turbulence affects consistency of the simulation results while the inconsistency keeps under 10%, which illustrates that phase screen method is suitable for being used under different turbulence conditions.

      Figure 3.  Propagation characteristics versus different $\omega $ ,(a) long-exposure beam radius WLE, (b) beam centroid displacement ${\sigma _c}$ , (c) on-axis scintillation index $\sigma _I^2$

      Figure 4.  Propagation characteristics versus different $\varepsilon $ ,(a) long-exposure beam radius WLE, (b) beam centroid displacement ${\sigma _c}$ , (c) on-axis scintillation index $\sigma _I^2$

      Figure 5.  Propagation characteristics versus different ${\chi _T}$ , (a) long-exposure beam radius WLE, (b) beam centroid displacement $ {\sigma _c} $ , (c) on-axis scintillation index $ \sigma _I^2 $

      Similar to the results in Fig.2, the on-axis scintillation index $\sigma _I^2$ show good match in weak fluctuation regime, the error keeps increase when the turbulence gets stronger $\sigma _I^2 > 1$ . The similar significant mismatch between numerical simulations and theoretical results for the scintillation index have been reported in earlier works about atmospheric turbulence. The reason is explained as the impact of either inner scale effect[14-15] or beam wander induced scintillations[16]. Ref.[12] proves that the major reason for this mismatch is the irregular appearance of large-amplitude (giant) irradiance spikes whose amplitudes greatly exceed the diffraction-limited peaks.

    • In this paper, the analytical formulas of beam spreading, beam wander and scintillation of Gaussian beam through oceanic turbulence are derived. The principle of parameter setting for phase screen simulated oceanic turbulence is mainly discussed and the steps are given in Sec.2.2. Furthermore, the comparison of the simulation results and theory calculation under different parameter settings and different turbulence conditions are shown in figures.

      Results show that beam spreading and beam wander characteristics simulated by phase screen method are in good accordance with the theory formula under different turbulence conditions while the scintillation index has large deviation under strong fluctuation turbulence conditions. In summary, based on appropriate parameter setting, phase screen method can be a simple and effective way to calculate the statistical characteristics depend on the first moment of intensity, such as beam spreading and beam wander. The statistical characteristics depend on the second moment of intensity calculated by numerical simulation achieve good correspondence with theory in weak fluctuation regime and big mismatch in strong fluctuation regime, which should be further researched in the future.

参考文献 (16)

目录

    /

    返回文章
    返回