高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

空心光束在非Kolmogorov湍流传输路径上的区域分割

陈晓文 魏小琴 汤明玥 邓涵凌

引用本文:
Citation:

空心光束在非Kolmogorov湍流传输路径上的区域分割

    作者简介: 陈晓文(1984-), 男, 副教授, 主要从事激光传输与变换方面的研究。E-mail: 2098647523@qq.com.
  • 基金项目:

    南充市科学技术局市校合作科研专项资金资助项目 18SXHZ0163

    南充市科学技术局市校合作科研专项资金资助项目 19SXHZ0239

    四川省教育厅资助项目 17ZB0163

  • 中图分类号: O436;TN012

Region division of hollow beams in non-Kolmogorov turbulent path

  • CLC number: O436;TN012

  • 摘要: 为了研究空心光束在非Kolmogorov湍流传输路径上的区域范围与各参量之间的关系及不同区域内光束的扩展情况, 采用广义惠更斯-菲涅耳原理推导了空心光束传输于非Kolmogorov湍流中的二阶矩宽度、瑞利区间及湍流距离的解析式, 并利用湍流距离把传输路径分割为3个区域进行数值分析。结果表明, 区域Ⅰ、区域Ⅱ的长度及区域Ⅲ的起始点都随湍流广义指数α的增大而先减小再增大(当α=3.11时出现一个极小值), 且随遮拦比η和光束阶数M(及N)的增大而增大; M(及N)取值较小时(M(及N) < 3), 湍流在瑞利区间内对光束扩展造成的影响不能忽略, M(及N)和η越大, 越容易忽略湍流在瑞利区间内对光束扩展所构成的影响; 光束在传输路径上依次进入区域Ⅰ、区域Ⅱ及区域Ⅲ, 其光束扩展逐渐变得更加剧烈, 且随着M(及N)和η的增加, 区域Ⅱ长度和区域Ⅲ的起始点相较于区域Ⅰ的长度增加更为显著。该研究结果为空心光束传输于湍流中的相关应用提供了参考。
  • Figure 1.  Contour of intensity distribution of hollow beams

    Figure 2.  zt, 1, zR versus α

    Figure 3.  zt, 1/zR versus η with different value of versus M(and N)

    Figure 4.  zt, 1 versus η with different value of M(and N)

    Figure 5.  Δzt versus α with different value of η

    Figure 6.  Δzt versus α with different value of M(and N)

    Figure 7.  zt, 2 versus α

    Figure 8.  zt, 2 versus η with different value

    Figure 9.  w(z) versus z with different value of M(and N)

    Figure 10.  w(z) versus z with different value of η

  • [1]

    WANG S C H, PLONUS M A. Optical beam propagation for a partially coherent source in the turbulent atmosphere[J]. Journal of the Optical Society of America, 1979, A69(9): 1297-1304.
    [2]

    GBUR G, WOLF E. Spreading of partially coherent beams in random media[J]. Journal of the Optical Society of America, 2002, A19(8): 1592-1598. 
    [3]

    CHEN X W, JI X L. Influence of turbulence on spatial correlation properties of partially coherent annular beams[J]. Chinese Journal of Lasers, 2009, 36(9): 2319-2325(in Chinese). doi: 10.3788/CJL20093609.2319
    [4]

    TANG M Y, CHEN X W, JI X L. Influence of turbulence on propagation and far-field beam quality of M×N cosh-Gaussian bams[J]. Acta Photonica Sinica, 2009, 38(3): 713-718 (in Chinese).
    [5]

    CHEN X W, JI X L. Consistency of the directionality of partially coherent beams in turbulence expressed in terms of the angular spread and the far-field average intensity[J]. Chinese Physics, 2010, B19(2): 024203. 
    [6]

    TOSELLI I, ANDREWS L C, PHILLIPS R L, et al. Free space optical system performance for laser beam propagation through non Kolmogorov turbulence[J]. Optical Engineering, 2008, 47(2): 023003. 
    [7]

    DENG J P, JI X L, LU L. Propagation of polychromatic partially coherent decentred laser beams propagating in non-Kolmogorov turbulence[J]. Acta Physica Sinica, 2013, 62(14): 144211(in Chin-ese). doi: 10.7498/aps.62.144211
    [8]

    XU H F, CUI Z F, QU J. Propagation of elegant Laguerre-Gaussian beam in non-Kolmogorov turbulence[J]. Optics Express, 2011, 19(22): 21163-21173. doi: 10.1364/OE.19.021163
    [9]

    ZHOU Zh L, YUAN Y Sh, SHU J, et al. Beam wander of a partially coherent-like beam in non-Kolmogorov turbulence[J]. Laser Techno-logy, 2019, 43(4): 143-148 (in Chinese). 
    [10]

    LIU D J, ZHONG H Y, WANG G Q, et al. Propagation of a radial phase-locked partially coherent elegant Laguerre-Gaussian beam array in non-Kolmogorov medium[J]. Applied Physics, 2019, B 125(3): 52.1-52.14. doi: 10.1007/s00340-019-7161-8
    [11]

    LUO Ch K, LU F, YIN Ch X, et al. Numerical study on transmi-ssion performance of Laguerre-Gaussian beam in non-Kolmogorov turbulence[J]. Laser & Optoelectronics Progress, 2018, 55(12): 120101 (in Chinese).
    [12]

    CHENG M J, GUO L X, LI J T. Influence of moderate-to-strong a-nisotropic non-Kolmogorov turbulence on intensity fluctuations of a Gaussian-Schell model beam in marine atmosphere[J]. Chinese Physics, 2018, B27(5): 054203. 
    [13]

    WU G H, GUO H, YU S, et al. Spreading and direction of Gaussian-Schell model beam through a non-Kolmogorov turbulence[J]. Optics Letters, 2010, 35(5): 715-717. doi: 10.1364/OL.35.000715
    [14]

    JI X L, ZHONG E T, LV B D. Spreading of partially coherent flattened Gaussian beams propagating through turbulent media[J]. Journal of Modern Optics, 2006, 53(12): 1753-1763. doi: 10.1080/09500340600590588
    [15]

    TANG M Y, LI B Zh. Turbulence distance of cosh-Gaussian beams in non-Kolmogorov turbulence[J]. Laser Technology, 2015, 39(4): 581-584 (in Chinese). 
    [16]

    CHEN X W, TANG M Y, HE C, et al. Spreading of apertured partially coherent beams in turbulent media[J]. Optik, 2014, 125(20): 6032-6036. doi: 10.1016/j.ijleo.2014.06.098
    [17]

    YANG J, WANG H, ZHANG X. Influence of atmospheric turbulence on the Rayleigh range of partially coherent laser[J]. Laser Technology, 2016, 40(3): 456-460 (in Chinese).
    [18]

    YANG T, JI X L. Regions of spreading of laser beams propagating through non-Kolmogorov atmospheric turbulence[J]. Acta Optica Sinica, 2015, 35(5): 0501001(in Chinese). doi: 10.3788/AOS201535.0501001
    [19]

    LIU D J, LUO X X, YIN H M. Effect of optical system and turbulent atmosphere on the average intensity of partially coherent flat-topped vortex hollow beam[J]. Optik, 2017, 130(C): 227-236. 
    [20]

    HONG X R, ZHENG Y L, TANG R A, et al. Propagation characteristics of a hollow Gaussian laser beam in a tapered plasma channel[J]. Physics of Plasmas, 2020, 27(4): 043109. doi: 10.1063/1.5145348
    [21]

    ALAVINEJAD M, TAHERABADI G, HADILOU N, et al. Changes in the coherence properties of partially coherent dark hollow beam propagating through atmospheric turbulence[J]. Optics Communications, 2013, 288(1): 1-6. 
    [22]

    XU Y G, YANG T, DAN Y Q, et al. Average intensity and spreading of partially coherent dark hollow beam through the atmospheric turbulence along a slant path[J]. Optik, 2016, 127(19): 7794-7802. doi: 10.1016/j.ijleo.2016.05.081
  • [1] 汤明玥李宾中 . 双曲余弦高斯光束在非Kolmogorov湍流中的湍流距离. 激光技术, 2015, 39(4): 581-584. doi: 10.7510/jgjs.issn.1001-3806.2015.04.034
    [2] 赵琦郝红宇樊红英李建欣蒋泽伟肖星 . 部分相干cosh-Gaussian光束通过大气湍流后的聚焦特性. 激光技术, 2016, 40(5): 750-755. doi: 10.7510/jgjs.issn.1001-3806.2016.05.028
    [3] 杨盛凯王晓艳赵亮徐勇根 . 反常涡旋光束在各向异性大气湍流中的漂移. 激光技术, 2021, 45(2): 252-258. doi: 10.7510/jgjs.issn.1001-3806.2021.02.022
    [4] 杨军王慧章曦 . 大气湍流对部分相干激光瑞利区间影响的研究. 激光技术, 2016, 40(3): 456-460. doi: 10.7510/jgjs.issn.1001-3806.2016.03.033
    [5] 姜楠李晓英牛春晖刘鑫 . 大气湍流对激光空间传输特性影响的实验研究. 激光技术, 2022, 46(5): 708-712. doi: 10.7510/jgjs.issn.1001-3806.2022.05.022
    [6] 束杰屈军 . EEMGSM光束在各向异性湍流中的光束质量因子研究. 激光技术, 2019, 43(6): 834-840. doi: 10.7510/jgjs.issn.1001-3806.2019.06.019
    [7] 周正兰袁扬胜束杰徐翔屈军 . 部分相干月牙形光束在非Kolmogorov谱中的漂移. 激光技术, 2019, 43(4): 579-584. doi: 10.7510/jgjs.issn.1001-3806.2019.04.025
    [8] 高金全周正兰徐华锋吴彬屈军 . 部分相干Airy涡旋光束在非Kolmogorov谱中的模态强度. 激光技术, 2021, 45(4): 522-529. doi: 10.7510/jgjs.issn.1001-3806.2021.04.018
    [9] 刘兵王巨胜杨泽后李晓锋樊冬任鹏李斌罗雄冯力天 . 基于非视线红外激光大气散射通信技术研究. 激光技术, 2014, 38(6): 854-858. doi: 10.7510/jgjs.issn.1001-3806.2014.06.028
    [10] 赵琦樊红英李轶国蒋泽伟胡绍云赖庚辛黄燕琳耿旭 . 大能量中空光束大气传输的仿真与实验比对研究. 激光技术, 2014, 38(4): 542-545. doi: 10.7510/jgjs.issn.1001-3806.2014.04.021
    [11] 刘玉丽 . 探测边界层大气温度的转动喇曼激光雷达. 激光技术, 2018, 42(4): 541-544. doi: 10.7510/jgjs.issn.1001-3806.2018.04.020
    [12] 刘琳琳杨杰黄见苑克娥尹凯欣胡顺星 . 激光雷达观测淮南大气SO2和NO2浓度廓线实例分析. 激光技术, 2019, 43(3): 353-358. doi: 10.7510/jgjs.issn.1001-3806.2019.03.012
    [13] 石玉立杨丰恺曹念文 . MPL反演南京北郊气溶胶光学厚度准确度的研究. 激光技术, 2015, 39(3): 372-376. doi: 10.7510/jgjs.issn.1001-3806.2015.03.020
    [14] 钱天陆健唐懋张冲张宏超 . 飞秒激光致双液滴光学击穿和等离子体分布研究. 激光技术, 2023, 47(2): 193-199. doi: 10.7510/jgjs.issn.1001-3806.2023.02.006
    [15] 王玉诏 . 一种用于激光雷达信号的自适应分段平滑算法. 激光技术, 2014, 38(5): 599-602. doi: 10.7510/jgjs.issn.1001-3806.2014.05.005
    [16] 杨建清罗积军徐军孙琦云王婷 . 基于雨滴谱分布数据库的近红外辐射衰减研究. 激光技术, 2018, 42(2): 161-165. doi: 10.7510/jgjs.issn.1001-3806.2018.02.004
    [17] 陈莉英姜洪波沈利沣赵志刚 . 脉冲腔衰荡技术探测气溶胶消光系数的研究. 激光技术, 2016, 40(2): 241-244. doi: 10.7510/jgjs.issn.1001-3806.2016.02.019
    [18] 刘秋武陈亚峰王杰王晓宾曹开法黄见胡顺星 . 差分吸收NO2激光雷达光源的设计与实现. 激光技术, 2018, 42(4): 433-439. doi: 10.7510/jgjs.issn.1001-3806.2018.04.001
    [19] 刘杰辉李鑫 . 基于ZEMAX二氧化硫荧光采集光路的设计仿真. 激光技术, 2020, 44(2): 221-225. doi: 10.7510/jgjs.issn.1001-3806.2020.02.015
    [20] 闫顺生李海花 . 利用喇曼-米方法确定激光雷达几何因子. 激光技术, 2013, 37(4): 511-514. doi: 10.7510/jgjs.issn.1001-3806.2013.04.021
  • 加载中
图(10)
计量
  • 文章访问数:  5202
  • HTML全文浏览量:  3701
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-07-08
  • 录用日期:  2020-08-04
  • 刊出日期:  2021-05-25

空心光束在非Kolmogorov湍流传输路径上的区域分割

    作者简介: 陈晓文(1984-), 男, 副教授, 主要从事激光传输与变换方面的研究。E-mail: 2098647523@qq.com
  • 1. 川北医学院 医学影像学院, 南充 637000
  • 2. 川北医学院 基础医学院, 南充 637000
基金项目:  南充市科学技术局市校合作科研专项资金资助项目 18SXHZ0163南充市科学技术局市校合作科研专项资金资助项目 19SXHZ0239四川省教育厅资助项目 17ZB0163

摘要: 为了研究空心光束在非Kolmogorov湍流传输路径上的区域范围与各参量之间的关系及不同区域内光束的扩展情况, 采用广义惠更斯-菲涅耳原理推导了空心光束传输于非Kolmogorov湍流中的二阶矩宽度、瑞利区间及湍流距离的解析式, 并利用湍流距离把传输路径分割为3个区域进行数值分析。结果表明, 区域Ⅰ、区域Ⅱ的长度及区域Ⅲ的起始点都随湍流广义指数α的增大而先减小再增大(当α=3.11时出现一个极小值), 且随遮拦比η和光束阶数M(及N)的增大而增大; M(及N)取值较小时(M(及N) < 3), 湍流在瑞利区间内对光束扩展造成的影响不能忽略, M(及N)和η越大, 越容易忽略湍流在瑞利区间内对光束扩展所构成的影响; 光束在传输路径上依次进入区域Ⅰ、区域Ⅱ及区域Ⅲ, 其光束扩展逐渐变得更加剧烈, 且随着M(及N)和η的增加, 区域Ⅱ长度和区域Ⅲ的起始点相较于区域Ⅰ的长度增加更为显著。该研究结果为空心光束传输于湍流中的相关应用提供了参考。

English Abstract

    • 激光束在大气中的传输研究对光通讯、卫星遥感等实际应用有一定指导意义[1],值得科研工作者对其展开相关研究。然而,复杂的大气环境需要找到一个合适的湍流模型来进行描述,前期科研工作者们均采用常规的Kolmogorov功率谱来研究湍流[2-5],但实际湍流中含有各向异性统计特征的区域,由湍流广义指数和广义幅度因子描述的非Kolmogorov湍流功率谱[6]所得结果与试验数据更相符,具有一定的普适性。因此,国内外同行们对激光束在非Kolmogorov湍流模型中的传输特性研究产生了浓厚的兴趣,并于近期展开了一些相关研究[7-12]。研究发现,光束在湍流中传输的光强分布、相干度、扩展等均对湍流广义指数表现出非单调的变化。光束扩展是激光束传输特性中值得重点关注的一个方面,目前已有大量的文献[13-18]中对激光束在湍流中的扩展展开了相关研究,并得到了一些有意义的结论。空间衍射及湍流介质是造成光束扩展的两个主要因素,那么值得进一步讨论的是,可以考虑把光束的传输路径分割为3个区域:主要因空间衍射导致光束扩展的区域、空间衍射和湍流共同影响光束扩展的区域及主要由湍流导致光束扩展的区域。因此,本研究中利用湍流距离将光束的传输路径分割为3个区域,并对3个区域的范围及不同区域内光束的扩展进行详尽的讨论。

      另一方面,空心光束的中心光强分布为零,具有特殊的性质和实际应用[19-20]。空心光束在生物医学、激光加工及微观粒子的冷却与囚禁等方面应用广泛,关于空心光束的产生、特性仍是目前研究的热点。为此,科研工作者们就湍流对空心光束传输特性的影响开展了大量研究[21-22]。然而,迄今为止,关于空心光束在湍流中的区域分割研究尚未涉及。

      综上所述,本文中依据湍流距离,将空心光束在湍流中的传输路径分割为3个区域,详尽讨论了3个区域的范围与光束参量和湍流参量之间的关系,同时直观地展示了光束在3个区域内的扩展情况。此外,详细讨论了光束参量的选取与湍流是否在瑞利区间内就对光束扩展构成影响之间的关系。

    • 空心光束的理论模型可由两束不同束宽(w0w0η分别为两束光束的束宽,且0 < η < 1)的平顶光束之差构成。那么,空心光束的复振幅可写为:

      $ \begin{array}{c} U_{0}\left(\boldsymbol{r}^{\prime}, 0\right)=\sum\limits_{m=1}^{M} \alpha_{m} \exp \left[-\left(m p_{m} \frac{x^{\prime 2}}{w_{0}{}^{2}}\right)\right] \sum\limits_{n=1}^{N} \alpha_{n} \times \\ \exp \left[-\left(n p_{n} \frac{y^{\prime 2}}{w_{0}^{2}}\right)\right]-\sum\limits_{m=1}^{M} \alpha_{m} \exp \left[-\left(m p_{m} \frac{x^{\prime 2}}{\eta^{2} w_{0}{}^{2}}\right)\right] \times \\ \sum\limits_{n=1}^{N} \alpha_{n} \exp \left[-\left(n p_{n} \frac{y^{\prime 2}}{\eta^{2} w_{0}{}^{2}}\right)\right] \end{array} $

      (1)

      式中,η为光束遮拦比;M, N为光束阶数; αm, αn, pmpn的表达式详见参考文献[3]

      则空心光束在z=0处的光强表达式为[3]

      $ \begin{array}{c} I_{0}\left(\boldsymbol{r}_{1}{ }^{\prime}, \boldsymbol{r}_{2}, 0\right)=U_{0}{ }^{*}\left(\boldsymbol{r}_{1}{ }^{\prime}, 0\right) U_{0}\left(\boldsymbol{r}_{2}{ }^{\prime}, 0\right)= \\ I_{01}\left(\boldsymbol{r}_{1}{ }^{\prime}, \boldsymbol{r}_{2}{ }^{\prime}, 0\right)+I_{02}\left(\boldsymbol{r}_{1}{ }^{\prime}, \boldsymbol{r}_{2}{ }^{\prime}, 0\right)- \\ I_{03}\left(\boldsymbol{r}_{1}{ }^{\prime}, \boldsymbol{r}_{2}{ }^{\prime}, 0\right)-I_{04}\left(\boldsymbol{r}_{1}{ }^{\prime}, \boldsymbol{r}_{2}{ }^{\prime}, 0\right) \end{array} $

      (2)

      式中,U0是环状光束的复振幅,U0*U0的共轭;r1′和r2′分别是U0函数、U0*函数中z=0处垂直于传输方向平面上的空间位置。其中,

      $ \begin{array}{c} I_{01}\left(\boldsymbol{r}_{1}{ }^{\prime}, \boldsymbol{r}_{2}{ }^{\prime}, 0\right)=\sum\limits_{m=1}^{M} \sum\limits_{m^{\prime}=1}^{M} \alpha_{m} \alpha_{m^{\prime}} \times \\ \exp \left[-\left(m p_{m} \frac{x_{1}{ }^{\prime 2}}{w_{0}{ }^{2}}+m^{\prime} p_{m^{\prime}} \frac{x_{2}{ }^{\prime 2}}{w_{0}{ }^{2}}\right)\right] \times \\ \sum\limits_{n=1}^{N} \sum\limits_{n^{\prime}=1}^{N} \alpha_{n} \alpha_{n^{\prime}} \exp \left[-\left(n p_{n} \frac{y_{1}{ }^{\prime 2}}{w_{0}{ }^{2}}+n^{\prime} p_{n^{\prime}} \frac{y_{2}{ }^{\prime 2}}{w_{0}{ }^{2}}\right)\right] \end{array} $

      (3)

      $ \begin{array}{c} I_{03}\left(\boldsymbol{r}_{1}{ }^{\prime}, \boldsymbol{r}_{2}{ }^{\prime}, 0\right)=\sum\limits_{m=1}^{M} \sum\limits_{m^{\prime}=1}^{M} \alpha_{m} \alpha_{m^{\prime}} \times \\ \exp \left[-\left(m p_{m} \frac{x_{1}{ }^{\prime 2}}{w_{0}{ }^{2}}+m^{\prime} p_{m^{\prime}} \frac{x_{2}{ }^{\prime 2}}{\eta^{2} w_{0}{ }^{2}}\right)\right] \times \\ \sum\limits_{n=1}^{N} \sum\limits_{n^{\prime}=1}^{N} \alpha_{n} \alpha_{n^{\prime}} \exp \left[-\left(n p_{n} \frac{y_{1}{ }^{\prime 2}}{w_{0}{ }^{2}}+n^{\prime} p_{n^{\prime}} \frac{y_{2}{ }^{\prime 2}}{\eta^{2} w_{0}{ }^{2}}\right)\right] \end{array} $

      (4)

      式中,x1′和y1′是r1′中的直角坐标;x2′和y2′是r2′中的直角坐标;m′,n′是分别对应x2′和y2′的取值阶数。将(3)式中的w0换为ηw0,得到I02(r1′, r2′, 0);将(4)式中的w0ηw0互换,得到I04(r1′, r2′, 0)。

      为展示空心光束在z=0处光强分布等高线图,可令(3)式中r1′=r2′,并选取光束参量w0=0.02m,η=0.6,M=N=4, 见图 1

      Figure 1.  Contour of intensity distribution of hollow beams

    • 在光束的传输路径上,光束的二阶矩宽度定义为[3]

      $ w^{2}(z)=\frac{2 \iint \boldsymbol{r}^{2} I(\boldsymbol{r}, z) \mathrm{d}^{2} \boldsymbol{r}}{\iint I(\boldsymbol{r}, z) \mathrm{d}^{2} \boldsymbol{r}} $

      (5)

      式中,I(r, z)为空心光束传输路径上z处的光强分布,r为传输路径上的空间位置。利用广义惠更斯-菲涅耳公式,I(r, z)可表示为:

      $ \begin{array}{*{20}{c}} {I(\mathit{\boldsymbol{r}}, z) = {{\left( {\frac{k}{{2{\rm{ \mathit{ π} }}z}}} \right)}^2}\iint {{{\rm{d}}^2}} {\mathit{\boldsymbol{r}}_1}^\prime \iint {{{\rm{d}}^2}} {\mathit{\boldsymbol{r}}_2}^\prime {I_0}\left( {{\mathit{\boldsymbol{r}}_1}^\prime , {\mathit{\boldsymbol{r}}_2}^\prime , 0} \right) \times }\\ {\exp \left\{ {\left( {\frac{{{\rm{i}}k}}{{2z}}} \right)\left[ {\left( {{\mathit{\boldsymbol{r}}_1}^{\prime 2} - {\mathit{\boldsymbol{r}}_2}^{\prime 2}} \right) - 2\left( {\mathit{\boldsymbol{r}} \cdot {\mathit{\boldsymbol{r}}_1}^\prime - \mathit{\boldsymbol{r}} \cdot {\mathit{\boldsymbol{r}}_2}^\prime } \right)} \right]} \right\} \times }\\ {{{\left\langle {\exp \left[ {{\psi ^*}\left( {\mathit{\boldsymbol{r}}, {\mathit{\boldsymbol{r}}_1}^\prime , z} \right) + \psi \left( {\mathit{\boldsymbol{r}}, {\mathit{\boldsymbol{r}}_2}^\prime , z} \right)} \right]} \right\rangle }_{\rm{m}}}} \end{array} $

      (6)

      式中,波数k=2π/λ; ψ(r, r2′, z)是湍流介质的复相位结构函数,ψ*(rr1′, z)是其共轭函数。湍流系综平均〈〉m表示为[4-5]

      $ \begin{array}{c} \left\langle\exp \left[\psi^{*}\left(\boldsymbol{r}, \boldsymbol{r}_{1}{}^{\prime}, z\right)+\psi\left(\boldsymbol{r}, \boldsymbol{r}_{2}{}^{\prime}, z\right)\right]\right\rangle_{\mathrm{m}}= \\ \exp \left\{-4 {\rm{ \mathit{ π}}}^{2} k^{2} z \int_{0}^{1} \int_{0}^{\infty} k \varPhi_{n}(\kappa, \alpha) \times\right. \\ \left.\left[1-\mathrm{J}_{0}\left(\kappa \xi\left|\boldsymbol{r}_{2}{}^{\prime}-\boldsymbol{r}_{1}{}^{\prime}\right|\right)\right] \mathrm{d} \kappa \mathrm{d} \xi\right\} \end{array} $

      (7)

      式中,J0()为零阶贝塞尔函数,κ表示空间频率,ξ为传输路径参量,κm为低空间截止频率, κ0为高空间截止频率。

      根据非Kolmogorov湍流统计[6], Φn(κ,α)=$F(\alpha )\tilde C_n^2\exp \left[ { - \left( {\frac{{{\kappa ^2}}}{{{\kappa _m}^2}}} \right)} \right]/{\left( {{\kappa ^2} + {\kappa _0}^2} \right)^{\alpha /2}}$,其中, $\tilde C_n^2$为折射率结构参量;高空间截止频率${\kappa _0} = \frac{{2{\rm{ \mathsf{ π} }}}}{{{L_0}}}$,低空间截止频率${\kappa _m} = \frac{{t(\alpha )}}{{{l_0}}}$(L0l0分别为湍流的外尺度和内尺度);F(α)=Γ(α-1)×$\frac{{\cos \left( {\alpha {\rm{ \mathsf{ π} /2}}} \right)}}{{4{{\rm{ \mathsf{ π} }}^2}}}$,$t(\alpha ) = {\left[ {\Gamma \left( {5 - \frac{\alpha }{2}} \right) \cdot F(\alpha ) \cdot \frac{2}{{3{\rm{ \mathsf{ π} }}}}} \right]^{1/(\alpha - 5)}}$,其中Γ(·)为伽马函数,α表示湍流广义指数。

      综合(2)式、(6)式及(7)式,并代入(5)式,采用积分变换法,并经繁琐的积分运算,可求得空心光束传输于非Kolmogorov湍流中的二阶矩宽度解析式:

      $ w^{2}(z)=w_{1}{ }^{2}+w_{2}{ }^{2} z^{2}+T z^{3} $

      (8)

      其中,

      $ \begin{aligned} w_{1}{ }^{2}=& w_{0}{ }^{2}\left\{\sum _ { m = 1 } ^ { M } \sum _ { m ^ { \prime } = 1 } ^ { M } \sum _ { n = 1 } ^ { N } \sum _ { n ^ { \prime } = 1 } ^ { N } \left[\frac{P}{\sqrt{Q_{1} Q_{2}}}\left(\frac{1}{Q_{1}}+\frac{1}{Q_{2}}\right)\left(1+\eta^{4}\right)-\right.\right.\\ &\left.\left.\frac{P \eta^{4}}{\sqrt{Q_{3} Q_{4}}}\left(\frac{1}{Q_{3}}+\frac{1}{Q_{4}}\right)-\frac{P \eta^{4}}{\sqrt{Q_{5} Q_{6}}}\left(\frac{1}{Q_{5}}+\frac{1}{Q_{6}}\right)\right]\right\} C^{-1} \end{aligned} $

      (9)

      $ \begin{array}{c} w_{2}{}^{2}=\left\{\sum \limits_ { m = 1 } ^ { M } \sum \limits_ { m ^ { \prime } = 1 } ^ { M } \sum \limits_ { n = 1 } ^ { N } \sum \limits_ { n ^ { \prime } = 1 } ^ { N } \frac { 4 } { w _ { 0 } ^ { 2 } k ^ { 2 } } \left[\frac{2 P}{\sqrt{Q_{1} Q_{2}}}\left(\frac{R_{1}}{Q_{1}}+\frac{R_{2}}{Q_{2}}\right)-\right.\right. \\ \left.\left.\frac{P \eta^{2}}{\sqrt{Q_{3} Q_{4}}}\left(\frac{R_{1}}{Q_{3}}+\frac{R_{2}}{Q_{4}}\right)-\frac{P \eta^{2}}{\sqrt{Q_{5} Q_{6}}}\left(\frac{R_{1}}{Q_{5}}+\frac{R_{2}}{Q_{6}}\right)\right]\right\} C^{-1} \end{array} $

      (10)

      $ \begin{array}{c} C=\sum\limits_{m=1}^{M} \sum\limits_{m^{\prime}=1}^{M} \sum\limits_{n=1}^{N} \sum\limits_{n^{\prime}=1}^{N}\left\{\frac{P}{\sqrt{Q_{1} Q_{2}}}\left(1+\eta^{2}\right)-\right. \\ \left.\left(\frac{P}{\sqrt{Q_{3} Q_{4}}}+\frac{P}{\sqrt{Q_{5} Q_{6}}}\right) \eta^{2}\right\} \end{array} $

      (11)

      $ T=\frac{2 {\rm{ \mathit{ π}}}^{2} z^{3}}{3} \int_{0}^{\infty} \kappa^{3} \varPhi_{n}(\kappa, \alpha) \mathrm{d} \kappa $

      (12)

      式中,P=αmαm′αnαn′Q1=mpm+m′pm′Q3=ηmpm+m′pm′Q5=mpm+ηm′pm′R1=(mpm)(m′pm′)。若将mm′分别改为nn′,即得Q2, Q4, Q6R2。(8)式中,湍流项T与湍流广义指数α、内尺度l0及外尺度L0相关,即α, l0L0的取值将影响空心光束的二阶矩宽度。

      瑞利区间的定义:光束的横截面积扩展达z=0处两倍时的传输距离[17],即:

      $ w^{2}\left(z_{\mathrm{R}}\right)=w_{1}{ }^{2}+w_{2}{ }^{2} z^{2}+T z^{3}=2 w_{1}{ }^{2} $

      (13)

      求解可得:

      $ z_{\mathrm{R}}=\frac{1}{3 T}\left(D+\frac{w_{2}{ }^{4}}{D}-w_{2}{ }^{2}\right) $

      (14)

      式中,$D = {\left[ {\frac{3}{2}T{{\left( {81w_1^4{T^2} - 12w_1^2w_2^6} \right)}^{1/2}} - w_2^6 + \frac{{27}}{2}w_1^2{T^2}} \right]^{1/3}}$。

      另一方面,湍流距离表示为光束的横截面积因受湍流影响而扩展10%及90%的传输距离[18],分别用zt, 1zt, 2表示:

      $ \frac{w^{2}\left(z_{\mathrm{t}}\right)-w_{\mathrm{f}}{ }^{2}\left(z_{\mathrm{t}}\right)}{w^{2}\left(z_{\mathrm{t}}\right)}=\mu $

      (15)

      式中,wf为自由空间中的束宽,wf2=w12+w22z2。求解三次方程,可得到空心光束传输于非Kolmogorov湍流中的湍流距离zt解析表达式:

      $ z_{\mathrm{t}}=\frac{\mu^{2} w_{2}{ }^{4}+\mu w_{2}{ }^{2} J+J^{2}}{3(1-\mu) J T} $

      (16)

      式中,

      $ \begin{array}{c} J=\left\{3(1- \mu) T \frac{\left[81 \mu^{2}(1-\mu)^{2} w_{1}{ }^{4} T^{2}+12 w_{1}{ }^{2} w_{2}{ }^{6}\right]^{1 / 2}}{2}+\right.\\ \left.\frac{27 \mu(1-\mu)^{2} w_{1}{ }^{2} T^{2}}{2}+\mu^{3} w_{2}{ }^{6}\right\}^{1 / 3} \end{array} $

      (17)

      显然,当μ=10%时,(16)式表示zt, 1;当μ=90%时,(16)式表示zt, 2

      应当指出的是,光束的传输路径可被zt,1zt, 2划分为3个区域:(1)区域Ⅰ:0 < zzt, 1,湍流介质对光束扩展所构成的影响甚小,可忽略,扩展主要由空间衍射造成;(2)区域Ⅱ:zt, 1 < z < zt, 2,空间衍射和湍流均对光束的扩展造成一定影响;(3)区域Ⅲ:zzt, 2,湍流为光束扩展的主要影响因素,空间衍射对其的影响可以忽略。

    • 基于空心光束二阶矩宽度表达式(8)式、瑞利区间表达式(14)式及湍流距离解析式(16)式,利用MATHEMATICA进行数值模拟,以分析空心光束在传输路径中的区域分割及在不同区域中的光束扩展。计算参量如下(不含特别指定的参量取值):$\tilde C_n^2$=10-14m3-α, w0=0.02m, L0=10m, l0=0.01m, λ=1.06μm, η=0.6, M=N=4。图 2为瑞利区间zR及湍流距离zt,1关于湍流广义指数α的变化情况。可以发现,zt, 1zR都随α的增大而先减小再增大,且当α=3.11时出现一个极小值。究其原因是此处湍流为极大值,zt,1zR受极大湍流的影响而成为极小值。从图 2还可以看出,当α>3.72时,zt, 1>zR,说明瑞利区间在区域Ⅰ之内,即湍流在瑞利区间内对光束扩展的影响可忽略不计;而α < 3.72时,zt, 1 < zR,瑞利区间已到达区域Ⅱ,说明湍流在瑞利区间内对光束扩展造成了一定的影响。为进一步定量比较zt, 1zR图 3中给出了zt, 1/zR在不同光束阶数M(及N)下随遮拦比η的变化情况。3条曲线均随M(及N)、η的增大而变大。M(及N)取值较小时(M(及N) < 3),不论η取何值,zt, 1均小于zR,即湍流在瑞利区间内对光束扩展构成的影响不能忽略;随M(及N)增大,曲线与虚线出现交点,且交点对应的M(及N)越大,η则越小。说明了M(及N)和η越大,越容易忽略湍流在瑞利区间内对光束扩展造成的影响。

      Figure 2.  zt, 1, zR versus α

      Figure 3.  zt, 1/zR versus η with different value of versus M(and N)

      图 4为区域Ⅰ长度zt, 1在不同M(及N)下随η的变化情况。可以发现,zt, 1M(及N)和η的增大而增大,其物理解释为:M(及N)和η越大,光束的原有扩展越大,湍流对其扩展的影响越小,则区域Ⅰ的长度越大。

      Figure 4.  zt, 1 versus η with different value of M(and N)

      区域Ⅱ的长度Δzt=zt, 2-zt, 1在不同ηM(及N)下随α的变化见图 5图 6。同样,Δzt也随α的增大而先减小再增大(当α=3.11时出现一个极小值),且ΔztηM(及N)的增大而增大。此外,α取较小值时,各条曲线差异较小,随α的增大,Δzt差异更为明显。

      Figure 5.  Δzt versus α with different value of η

      Figure 6.  Δzt versus α with different value of M(and N)

      图 7图 8中分别为区域Ⅲ的起点zt, 2αη的变化情况。由图可看出,zt, 2同样随α的增大而先减小再增大(当α=3.11时出现一个极小值),且随M(及N)和η的增大而增大。为了更加直观展示空心光束在湍流中的扩展,图 9图 10中给出了光束二阶矩宽度w(z)在不同M(及N)和η下随传输路径z的变化情况,同时给出了相应的zt, 1zt, 2的具体数值。由图可知,光束在传输路径上依次进入区域Ⅰ、区域Ⅱ及区域Ⅲ,其二阶矩宽度w(z)的增大变得更加剧烈,且zt, 1zt, 2M(及N)和η的增大而增大,这与上述研究结果一致;同时,发现随着参量M(及N)和η的增加,zt, 1缓慢增加,而zt, 2的增大却十分显著。这说明光束参量M(及N)和η对区域Ⅱ长度和区域Ⅲ的起始点影响更大。

      Figure 7.  zt, 2 versus α

      Figure 8.  zt, 2 versus η with different value

      Figure 9.  w(z) versus z with different value of M(and N)

      Figure 10.  w(z) versus z with different value of η

    • 建立空心光束理论模型,并推导其传输于湍流中的二阶矩宽度、瑞利区间及湍流距离解析式。同时,利用湍流距离把传输路径分割为3个区域,并对扩展区域的分割范围及不同区域内光束的扩展进行了数值理论分析。

      (1) 区域Ⅰ、区域Ⅱ的长度及区域Ⅲ的起始点都随α的增大而先减小再增大(当α=3.11时出现一个极小值),且随M(及N)和η的增大而增大。

      (2) M(及N)取值较小时(M(及N) < 3),不论η取何值,湍流在瑞利区间内对光束扩展的影响都不能忽略;M(及N)和η越大,越容易忽略湍流在瑞利区间内对光束扩展造成的影响。

      (3) α取较大值时,不同M(及N)和η下区域Ⅱ长度的差异更为显著。

      (4) 传输路径上,光束依次进入区域Ⅰ、区域Ⅱ及区域Ⅲ, 其光束扩展依次变得更加剧烈;随着参量M(及N)和η的增加,区域Ⅰ的长度缓慢增加,而区域Ⅱ长度和区域Ⅲ的起始点增加更为显著。应当特别指出:各类激光束在湍流中传输,都可以采用湍流距离将光束传输路径分割为3个区域,3个扩展区域的范围均与光束参量和湍流参量相关。本研究中所指出的区域Ⅰ、区域Ⅱ的长度及区域Ⅲ的起始点都随α的增大而先减小再增大,这一结论也能适用于其它各类光束。究其原因:本研究小组[4]曾在讨论相干光束在湍流中的扩展及方向性时指出,任意相干光束在湍流中二阶矩宽度的解析式均可表达为本文中的(8)式,即二阶矩宽度均由3项组成,前两项与光束参量相关,第3项湍流项是独立存在,因而各类光束的扩展受湍流影响的变化趋势是一致的,结合(7)式和(16)式可知,湍流广义指数对空心光束湍流距离(区域Ⅰ、区域Ⅱ的长度及区域Ⅲ的起始点)的影响也能适用于其它光束。本研究结果对于空心光束传输于湍流中的相关应用有指导意义。

参考文献 (22)

目录

    /

    返回文章
    返回