非均匀湍流路径上光传播数值模拟的相位屏分布.docVIP

第 卷 第 期 年 月 Q( = $$= = （ ） %$$$C#=$R$$=RQ( $= R’’##C$ 物 理 学 报 S@TS +UV,W@S ,WXW@S ， ， ， 0:8Q( X08= ,6P/6EO64 $$= $$= @95I8+9.N8 ,0K8 非均匀湍流路径上光传播数值模拟 的相位屏分布 ! 钱仙妹! 朱文越 饶瑞中 （中国科学院安徽光学精密机械研究所大气光学中心，合肥 ） #$$#% （ 年 月 日收到； 年 月 日收到修改稿） $$ %% ’ $$( # 针对非均匀湍流路径上光传播的数值模拟，研究了相位屏分布的三种不同设置方案，即均匀分布相位屏 （ ）、等 指数间隔的相位屏（ ）和等 参数间隔的相位屏（ ）结果表明， 方案根据折射率 )*+, -./01 2-+, 34567 23+, 8 2-+, 起伏的强弱设置密疏相间的相位屏，实现了对不同折射率起伏区的充分采样，是非均匀路径上相位屏分布的首选 方案 在满足相位屏间距基本要求的前提下，决定相位屏间隔的 指数间隔的选取具有较大的任意性，统计结 8 -./01 果显示，尽管相位屏间 指数变化超过一个数量级，但光束有效半径、能量 比的变化均小于 ，模拟结 -./01 ,/469: ;% 果与理论结果符合很好，表明了 方案的可靠性和稳定性 2-+, 8 关键词：光传播，非均匀湍流路径，数值模拟，相位屏分布 ： ， ， !## =’$ =’$2 =%$? 位屏数目至关重要 然而，据我们所知，目前尚没有 8 %; 引 言 相关文献对上述问题进行讨论，更没有给出明确的 判据来阐述如何设置非等间隔的相位屏 对于地基 8 在研究随机起伏介质中的光传播规律时，数值 激光大气上行传播而言，如何针对湍流强度具有明 模拟因其独特的优势而成为重要手段 常用的数值 显层状结构、随高度剧烈变化的路径设置相位屏是 8 模拟方法是从光传播的波动方程出发 ［ — % ］，在传播 ’ 当前亟待解决的问题之一 8 路径上设置若干个相位屏，光束传播到一个相位屏 本文采用多层相位屏模拟技术，数值分析了均 时即把对应的相位扰动叠加到光场中，然后真空传 匀 间 隔 相 位 屏 LI5M04E 75N/45OL/67 P9FN6 NK466I 播至下一个相位屏处，如此重复，直至终点为止 因 （ ）、等 指数间隔相位屏 8 )*+, -./01 6HL51F:6I/ -./01 此，在传播路径的恰当位置设置相位屏成为光传播 （ ）和等 参数间隔 5I76GC5I/641F: P9FN6 NK466I 2-+, 34567 数值模拟的关键问题之一 对于折射率起伏均匀的 相位屏 8 2HL51F:6I/ 34567 PF4FE6/64C5I/641F: P9FN6 NK466I 光传播路径，一般设置等间隔的相位屏；而对于非均 （ ）等非均匀路径情况下相位屏分布设置的三 23+, 匀路径下的光传播，如地球大气中的上行或下行传 种方案 通过数值模拟与理论结果的对比分析和可 8 播问题，由于传播路径上的折射率起伏强度不相等， 靠性验证，以确定非均匀随机起伏介质光传播数值 若仍设置等间隔的相位屏，则会导致强起伏区不能 模拟中相位屏设置的最佳方案 8 被充分采样，而弱起伏区被过度采样 这不仅使得数 8 ; 相位屏间隔的选择要求 值模拟不能真实地反映光传播的实际过程，又增大 了不必要的计算量，而且由于设置相位屏数目的主 观性，进一步增加了仿真结果的不确定性 因此，对 数值模拟结果的可靠性直接依赖于模拟参数的 8 于非均匀随机起伏介质中的光传播模拟而言，如何 根据折射率起伏的强弱程度设置非等间隔的相位屏 合理选择，其中相位屏间距 受到抽样原则、折射 ! ! 率起伏特性和光传播效应等要求的限制 ［ ，， — ］ % %% 8 以及根据折射率起伏强度和光传播距离自动确定相 相位屏间的真空传播求解是在傅里叶变换的频 中国科学院知识创新工程重要方向性项目（批准号： ）和国家自然科学基金（批准号： ）资助的课题 ! @ABBCD= $$D$% 8 ： ! 2CEF5: GEH5FI%=(J .F9008 K0E8 KI  @/B @ 物 理 学 报 卷 := 域中进行，这要求相邻相位屏间相邻频率上的相位 ，其中 ! + , 为光束的真空波数， , # % !! 差应小于 ，计算得到 ! ）等 / LMJ6N 参数间隔相位屏 （ - ! + + ） ——— * ， （ ） ! ! ! # # $ 其中 为网格间距， 为网格数目， 为激光波 ! # 长 同时，为了能以离散的相位屏代替连续的随机介 % 同理，另选常数 KLHI% 相位屏间湍流大气的 根据