大气光学中湍流微扰计算讨论.docVIP

大气光学中湍流微扰计算讨论 　　摘要：光在湍流大气中传输，理论分析湍流大气对光束性质的影响时，需要考虑湍流大气的微扰结构函数，本文首先采用微扰湍流谱的二阶统计平均分析光束在接受面的性质，再采用微扰湍流谱的波结构函数分析接受面处的光束性质，分析表明，采用两种不同的微扰湍流谱结构，可以得到一致的结果，即满足kξv=1条件时，湍流微扰谱密度统计平均的二阶近似和波结构函数表示是等价的。 　　关键词：大气光学；微扰计算；统计平均；结构函数 　　中图分类号：G642 文献标志码：A 文章编号：1674-9324（2014）24-0114-02 　　在大学光学教学中，一般都会讲到光传输问题。光传输指光在真空中（自由空间）、光学系统或随机介质中传输。光在随机介质中传输时，最常见的是光在湍流大气中传输，我们每天见到的阳光就是光穿越大气照射到地球上，正在发展中的激光无线光通讯也是需要光在大气中传输，光束在传输时受不均匀大气引起的折射、衍射，微小粒子散射等影响，使其在传播过程中强度、相位在时间和空间上呈现随机起伏，产生波前扭曲变形、光束抖动、光强闪烁、光束弯曲和漂移、光束扩展以及接收端光斑发生畸变等现象，严重影响信息在大气中传递和接收。如何处理大气扰动是必须解决的问题，在具体的教学或分析文献时，不同的书籍或文献对湍流大气的处理方式不同，同学们容易混淆，这里分析常见的一种情况。光在大气中传输受到大气中湍流的影响，一般采用惠更斯-菲涅尔原理加上湍流微扰项处理，即： 　　式中s=（x，y），r=（x，y）分别表示入射面和接收面处的位置矢量，U（s，0），U（r，L），分别表示入射面和接受面处光波的电场，i是虚数单位，d2s=dxdy，L是传输距离，k是波数，f是频率，ψ（s，r，L）是湍流微扰项。由方程（1）可以看出，如果知道入射光束的电场U（s，0），频率f和湍流微扰项ψ（s，r，L），就可以根据方程得到接收面处光束的电场，当接收面处的电场知道，可以分析光束的强度分布，光束的相干态、偏振态的变化以及光束的发散情况等；还可以分析光束在传输过程中光束的闪烁，漂移和信噪比等参数的变化情况。可是在实际的分析过程中，湍流微扰项给出的是大气折射率起伏的结构函数，不能直接计算，因此在求解方程（1）时通常采用谱密度〈U（r1，L）U*（r2，L）〉m的统计平均来表示，其中*表示取复共轭，〈〉m表示统计平均。 　　一种情况下，微扰项的统计平均用二阶近似表示： 　　式中φn（k）表示大气折射率指数起伏的空间能谱，具体的谱表示可以参考文献[2]，k是空间波数，单位是m-1。方程（2）给出了湍流微扰项的计算方式，而对于方程（1），需要对接收面处的电场做统计平均，即：〈W（r1，r2，L）〉=〈U（r1，L）U*（r2，L）〉，当r1=r2时，表示接收面处的光强，即〈I（r，L）〉=〈W（r，r，L）〉。从方程（2）可以看出，在计算接收面处光束的电场时，接收面处的点r和入射面处的点s有耦合项（r1-r2）（s1-s2），这在计算过程中增加了一定的难度。 　　另一种情况下，微扰项的统计平均用波结构函数表示： 　　式中J0表示零阶贝塞尔函数，从方程（3）中可以看出，接收面处的点r和入射面处的点s没有耦合项，在理论计算过程中降低了一定的难度，并且方程（3）在表示湍流扰动方面与实际情况符合也比较好。方程（2）和方程（3）都是基于科尔莫戈罗夫（Kolmogorov）描述湍流的三分之一标度律，描述的是简单的，理想化的均匀的，各项同性的湍流。 　　在计算模拟过程中，方程（2）和方程（3）都是正确的描述了湍流的扰动，那么在什么情况下使用方程（2），什么情况下使用方程（3），没有绝对的要求，具体使用哪一个方程比较好，根据实际情况判断。 　　在满足|kξv|=1条件时，式J0（kξv）：1-■（kξv）2成立，此时可以对方程（3）进行化简，即把贝塞尔函数展开成指数函数形式，具体计算如下： 　　-4π2k2L■dξ■[1-J0（k|s1-s2）ξ+（1-ξ）（r1-r2）|] 　　φn（k）kdk 　　=-π2k2L■φn（k）kdk■dξ×（|k2（s1-s2）2ξ2+k2（1-ξ）2（r1-r2）2+2ξ（1-ξ）k2（r1-r2）（s1-s2）|）2 　　=-π2k2L■φn（k）k3dk■dξ×（|（s1-s2）2ξ2+（1-ξ）2（r1-r2）2+2ξ（1-ξ）（r1-r2）（s1-s2）|）2 　　=-■■k3φn（k）dk[（r1-r2）2+（r1-r2）（s1-s2）+（s1-s2）2] （4） 　　方程（4）表明，在满足近似条件下，方程（3）可以展开成方程（2）的形式。但是方程（2）不可以变换成方程（3），因为方程（2）是在近似展开的条件下由方程（3）获得，因此，在实际计算过程中