大气湍流中光传播的数值模拟.doc

PAGE 2 PAGE 1 大气湍流中光传播的数值模拟收稿日期 收稿日期: 2010年3月14日 收到修改稿日期： 基金项目：教育部科技重点项目（105164） 通信作者: 马保科（1972），男，在读博士，副教授，主要从事随机介质中波传播方面的研究。Email: baokema2006@126.com 马保科1,2， 郭立新1 吴振森1 （1.西安电子科技大学，陕西西安 710071 2.西安工程大学，陕西西安 710048 ） 摘 要 光在大气湍流中传播时，受大气分子、气溶胶等粒子的相互作用，将发生光束扩展、漂移和相干性退化等大气湍流效应，这些因素严重影响了光波的远场特性。文章从大气湍流中光传播的理论研究入手，分析了如何构造较为合理的大气湍流相位屏。进而采用McGlamery算法，对Kolmogorov谱下的大气湍流随机相位屏进行了数值模拟，并分析了光波从发射机经湍流大气传播到达接收机时的远场变化特性。研究表明，大气湍流的存在对光的远场传播质量造成很大的影响，研究结果也为大气湍流中与光传播相关的工程应用及自适应光学技术的完善提供了参考。 关键词 大气湍流；McGlamery算法；相位屏模拟； 大气结构常数； 中图分类号 TP391 文献标识码 A 1 引言 大气湍流是一个相当复杂的随机媒质系统，虽然物理学界对湍流的研究已经历了相当漫长的历史，但因涉及的因素千头万绪，其间的相互作用和关系也错综复杂，人们对其物理本质至今未能做到较为清楚的认识。因此，光在大气湍流中传播问题的研究仍存在理论和实验上的挑战[1,2]。通常，当光在湍流大气中传播时，光束截面内包含着许多的大气漩涡，这些漩涡各自对照射到它的那一部分光束形成衍射作用，可导致光束的强度和相位随机变化，进而表现出光束扩展，大气闪烁和相位起伏等大气湍流效应，从而严重降低了接收机的接收效率。目前，突破大气湍流的影响仍是光在随机介质中传播所要解决的关键问题[3]。早在20世纪中期，苏联的Obukhov便采用Rytov平缓微扰法由实验反演湍流特征。在闪烁的饱和现象被发现之后，物理学界又将Markov近似引入求解光场的统计矩，研究大气湍流下的光场特性[1]。然而，在中等起伏条件下，目前仍没有找到很好的解析处理方法。由于数值模拟能够从光的传播过程出发，较为清楚地反映出所涉及问题的物理本质，因而成为研究湍流效应的主要方法[4]。本文采用McGlamery算法[5]，对Kolmogorov谱下的大气随机相位屏进行了数值模拟，进而结合Huygens-Fresnel原理，模拟了在有无大气湍流的情况下，接收机处光场的变化特性。 2 大气湍流中光的传播 在折射率为的随机媒质中，一束波长，波数为（）的单色波的电场由Maxwell波动方程来描述[1,4] 其中，，（，且）为湍流折射率， GOTOBUTTON ZEqnNum490134 REF ZEqnNum490134 \! (1)式左端最后一项反映了偏振特性，当波长远小于湍流的内尺度（）时，此项可以忽略不计，有 因而，电场的任一分量波动方程为 如果介质的非均匀尺度远大于波长，可认为只存在前向小角度散射而无后向散射。对如图1沿方向的传播进行旁轴近似处理，并将光场写为，得 其中，为横向算符。 GOTOBUTTON ZEqnNum993899 REF ZEqnNum993899 \! (4)式还可进一步化简为 其中，算子。如果折射率不随传播距离而变化（或变化较小），则 GOTOBUTTON ZEqnNum258731 REF ZEqnNum258731 \! (5)式中忽略交换算子项，得到如下的广义抛物型方程 通常，算子不同的近似可得不同的抛粅型方程，而算子最简单的近似即为其Taylor展式 当湍流场的折射率满足时， GOTOBUTTON ZEqnNum333191 REF ZEqnNum333191 \! (7)式化简为 将 GOTOBUTTON ZEqnNum565157 REF ZEqnNum565157 \! (8)式代入 GOTOBUTTON ZEqnNum123377 REF ZEqnNum12