小尺度流场内涡旋光束的气动光学特性

张永棠

(1.广东东软学院 计算机科学与技术系，广东 佛山，528225； 2.广东省大数据分析与处理重点实验室，广东 广州，510006； 3.南昌工程学院 江西省协同感知与先进计算技术研究所，江西 南昌，330003)

涡旋光束因其特殊的螺旋形相位和自身携带的轨道角动量(OAM)由Vig JR等人[1]提出后就引起了学术界的广泛讨论[2-3]。尤其是OAM可以改善光通信容量[4]，吸引了众多学者对涡旋光束在大气中的传播特性进行研究[5-6]。然而在一些应用环境中，涡旋光束需要经过空气密度高频变化的气动区域，其传输特性在现有文献中提及较少，因此对涡旋光束气动光学效应的研究是涡旋光束在光通讯领域内应用的进一步拓展。

目前在气动光学的研究中，例如Garciar等人[7]采用沿着光束行进路径积分计算OPD的方式，对经过具有平型窗口凸台周围流场的光束的气动光学效应进行了探究；Katos等人[8]通过搭建经过超音速流场的实验光路并以OPD为计算结果，研究了超音速流场中边界层处的气动光学效应；杨海马等人[9]对光束经过超音速湍流边界层后光畸变的统计特性进行了探讨；Song X等人[10]研究了剪切层与孔径对气动光学效应的影响，以上对于气动光学效应的研究，除了被研究光源均为高斯光束以外，光场计算的方法也均为基于费马原理的光线追迹法。

由于其他含有OAM的光束在通信时理论上可以看作是具有不同径向指数和相同OAM的LG光束的线性叠加[11]，并且LG光束是涡旋光束中最为基本的一种涡旋光，因此文中将研究具有3种不同拓扑荷数的LG光束，分别以攻角、海拔、速度作为流场唯一变化条件时，经过小尺度流场后的气动光学效应。构建了流场的气动模型和光场的传输模型，采用的是基于亥姆霍兹方程的数值仿真算法，光场的传输计算方法则有别于传统的光线追迹理论。在该小尺度流场模型上，对LG光束和高斯光束在不同流场环境下经过远场衍射后的SR和成像偏移进行研究，得出LG光束在小尺度流场内所表现出的气动光学效应规律性的结论以及与高斯光束有别的特性。

1 理论分析

1.1 气动模型

为了便于研究，文中构建了如图1所示的启动模型，LG光束所经过的流场为一个由微小凸起(以下简称凸起)造成的流场。凸起为一圆柱体加上半球体构成，置于平面中心，即凸起底面圆心位置为(0,0)。截取凸起正上方2×2×1.73 mm3的流场作为光束要经过的小尺度流场。在LG光束传播1.73mm后，可认为光束进入到折射率均匀的大气之中，一般而言，2mm左右的传播距离已经能够使得光束表现出气动光学应。选用小尺度流场则是因为小尺度时流场的雷诺数较小，湍流激烈程度也较小，在流场仿真中更易形成稳定可收敛的流场。

采用计算流体动力学(CFD)仿真时，设置飞行速度为0.6Ma至4.6Ma，飞行高度为6km至10km，攻角为0°至20°，流场仿真方法选用大涡模拟。此处的攻角定义为来流与凸起横轴之间的夹角，光束沿着凸起的纵轴垂直出射，观察方向沿着光束出射方向的反方向。

图1 流场尺寸示意图Fig.1 Flow field parameters

1.2 传输模型

在非均匀介质中，麦克斯韦方程可以简化为如下的亥姆霍兹方程：

2Ur+k2n2rUr=0

(1)

将光场复振幅Ur展开为抛物型表达式Ur=Ax,y,zeikz代入(1)式中，在满足旁轴近似∂zzA≪k∂zA时，(1)式可以演变为振幅A对传播方向z的隐式微分方程[12]：

(2)

(3)

光束出气动区域后，使用离散化的夫琅禾费衍射积分公式对远场光振幅进行计算[13]：

(4)

式中z0为衍射距离，2a、2b分别为y、x两方向的衍射面长度。Δx1=Δy1=0.01，Δz=0.001。

光束穿过小尺度流场后，其远场成像质量采用SR和成像位移两种参数分别从相位稳定性和相对于光信号接收器件而言的振幅稳定性两方面进行评价。SR通常采用的计算公式为[14-15]

SR=exp-k×Δrms2

(5)

其中Δrms为OPD的均方差。成像位移的计算公式为

(6)

其中，高斯光束在初始面处的振幅分布如下：

(7)

(8)

以上对光场的传输计算和对评价参数的计算，采用的软件均为Mathematica。

2 振幅和相位稳定性

为了更好的分析光束的振幅和相位稳定性，并且有较直观的仿真视觉图像，参照了文献[7,10,13]中的实验参数，本研究采取攻角取值为[0°,5°,10°,15°,20°]，海拔取值为[6km,7km,8km,9km,10km]，速度取值为[0.6Ma,1.6Ma,2.6Ma,3.6Ma,4.6Ma]的条件进行仿真。采用径向指数p=0，拓扑荷数l=1,2,3的LG光束(为方便描述，分别简记为L1、L2、L3)以及高斯光束作为仿真对象，分别探究攻角、海拔、马赫数的变化对LG光束光强和OPD的影响。

图2展示了上述实验参数下的初始光强分布和相位分布图像。其中图2(a)-(d)为光强分布、图2(e)-(h)相位分布。

从图2可以看出，不同拓扑荷数的LG光束以及高斯光束经过气动区域后，不同大小的光斑会携带不同位置的OPD数值，这些OPD数值会体现在实部和虚部的变化当中，变化了的实部虚部会在远场衍射中进一步影响光强的分布，使其充分地显示出气动光学效应。因此在计算光束的SR和成像偏移，可以研究光束的相位稳定性和振幅稳定性。

2.1 光束的相位稳定性

由(5)式可知，SR与光束的OPD均方差呈负相关关系，因此，本研究采用SR来描述光束的相位稳定性。图3是高斯光束和L1、L2、L3光束的SR值随攻角变化的曲线图。

图3 SR随攻角变化曲线图Fig.3 SR at different angles of attack

由图3可知，随着攻角的增大，各光束的SR均逐渐下降。对于LG光束而言，在同一攻角下，随着拓扑荷数的增加，SR也在增加，可见当攻角作为流场唯一变化条件时，拓扑荷数越大，光束的相位稳定性越好，并且LG光束的SR总是大于高斯光束的SR，故LG光束的相位稳定性高于高斯光束。此外，随着攻角的增大，不同光束间的SR差值在缩小，可见在大攻角的情况下，拓扑荷数对相位稳定性的影响不大，而在小攻角处，拓扑荷数大的LG光束具有明显的相位稳定优势。

图4展示了在0°攻角、2.6Ma的飞行速度下，高斯光束和L1、L2、L3光束SR值随海拔增加的变化曲线。

图4 SR随海拔变化曲线图Fig.4 SR at different altitudes

由图4可知，随着海拔的增大，各光束的SR值均在逐渐增大，各光束间SR的差值也在扩大，且LG光束的SR值依旧高于高斯光束。对于LG光束而言，在同一海拔下，随着拓扑荷数的增加，SR也在增加。所以在海拔作为流场唯一变化条件时，拓扑荷数越大，LG光束的相位稳定性越好，且相互之间的差距随着海拔的增大而扩大。

图5展示了在0°攻角，10km飞行高度下，高斯光束和L1、L2、L3光束的SR值随马赫数的变化曲线。

图5 SR随马赫数变化曲线图Fig.5 SR at different Mach Numbers

由图5可知，在马赫数由亚音速增加到超音速之时[15]，各光束的SR值都有一个大的跌落，但是对于L2和L3两束LG光束而言，当马赫数继续增加时，其SR值几乎维持不变，而L1光束与之相比，则依旧呈现出下降的趋势。可见当马赫数作为流场唯一变化条件时，LG光束的拓扑荷数越大，光束相位稳定性越好，且相互之间差距较大。当马赫数大于2Ma时，马赫数对拓扑荷数大于1的LG光束相位稳定性的影响不大。

图5中显示，高斯光束于3.6Ma处的SR值较其在2.6Ma时的SR值略大，约在2×10-22，而于4.6Ma又有所下降。该现象的可能原因是3.6Ma下的仿真流场，在光束路径所经过的流场部分，较大的空气密度浮动与2.6Ma相比处在所计算的截面较为外侧的区域，而相同束腰下高斯光束的光斑最小，所以于3.6Ma处高斯光束的SR又略有上升。

就相位稳定性而言，LG光束的拓扑荷数越大，相位稳定性越好，且均优于高斯光束。其中拓扑荷数大于1的LG光束在经过高马赫数飞行产生的流场后其相位稳定的优势非常明显。当且仅当攻角作为流场唯一变化条件时，随着攻角的增加，LG光束间的相位稳定性差距在减小。

2.2 光束的振幅稳定性

图6 成像偏移随着攻角的变化Fig.6 Imaging displacements at different angles of attack

图7 成像偏移随着海拔的变化Fig.7 Imaging displacements at different altitudes

当海拔、马赫数、攻角分别作为流场唯一变量时，采用成像偏移(以下简称偏移)来查看气动光学效应对光束质心位置产生的影响。高斯光束以及L1、L2、L3三种LG光束在10km、2.6Ma飞行速度下偏移量随攻角增大而变化的曲线图如图6所示。

由图6可以看出，随着攻角的增加，偏移也在逐渐增大。对于高斯光束和L1、L2两束LG光束而言，在攻角小于15°的时候，偏移变化幅度很小，而当攻角大于15°后，这三束光的偏移量均大大增加，可见对于这三束光，在一定的攻角范围内，均能较好的保持成像位置。并且发现L1光束的偏移量一直保持在高斯光束之下，可见其比高斯光束有着更好的振幅稳定性。而对于L3光束而言，随着攻角的增加，其偏移大幅增加。因此，对LG光束而言，拓扑荷数越大，其成像偏移越大，当拓扑荷数大于2时，振幅稳定性极差。在攻角为0°至15°的范围内时，对高斯光束和拓扑荷数小于3的LG光束而言，攻角变化对其偏移的影响并不明显。

高斯光束以及L1、L2、L3三种LG光束在0°攻角、2.6Ma飞行速度下偏移量随海拔增加而变化的曲线图如图7所示。

由图7可以看到，随着海拔的增加各光束的偏移量均呈现下降的趋势。L1、L2光束偏移量均明显小于L3，且三者之差越来越小。可见海拔越高，光束的振幅稳定性越好，虽然各光束之间的振幅稳定性差距随着海拔的增加而减小，但相互之间的差距仍然很大。

图8展示了高斯光束以及L1、L2、L3三种LG光束在0°攻角、10km飞行高度下随马赫数的增大，偏移量的变化曲线。

图8 成像偏移与马赫数的关系Fig.8 Imaging displacement at different Mach numbers

由图8可知，随着马赫数的增加，成像偏移量有增加有减少。与之前的规律相似，高斯光束和L1、L2两束LG光束的偏移量均明显小于L3光束，其中L1光束又小于高斯光束，并且在0.6Ma、1.6Ma处，L2的偏移量也小于高斯光束。对于高斯光束和L1、L2两束LG光束而言，在马赫数由亚音速升至超音速的过程中，偏移量有个极大的增加。而在速度由1.6升至4.6Ma的过程中，偏移量经历了下降、上升、下降的过程。对于L3光束而言，当速度从0.6Ma升至3.6Ma时，其偏移量一直在增加，只是增加的幅度不一样，但在4.6Ma处，偏移量有所下降。该现象的可能原因是当马赫数较大时，激波进一步靠近头部，雷诺数也在增加，因此边界层变薄[2]，故图像偏移量下降。所以在其它飞行条件不变时，各光束的振幅稳定性在马赫数大于音速后都有极大的降低，而当马赫数继续增加时，振幅稳定性则产生波动，没有稳定的走向。

就振幅稳定性而言，拓扑荷数越大，振幅稳定性越差，其中L3光束的振幅稳定性远远小于L1、L2光束，且仅有L1光束，其振幅稳定性总是优于高斯光束。在攻角为0°到15°时，攻角的变化对拓扑荷数小于3的LG光束的振幅稳定性没有明显影响。

3 结论

采用大涡模拟对小尺度流场进行了仿真，得到了流场区域内的空气密度数据，再由旁轴近似的数值算法，仿真计算了不同拓扑荷数的LG光束在经过小尺度流场后，又经过2m远场衍射后的OPD和光强分布，并以SR和成像偏移对远场光束质量进行了评价。在对不同拓扑荷数的LG光束远场成像质量的比较中发现，LG光束的拓扑荷数越大，其穿过小尺度气动区域后的相位稳定性就越好，且总是优于高斯光束；在高马赫数时，拓扑荷数大的LG光束相位稳定的优势相当明显；当且仅当攻角作为小尺度流场唯一变化条件增大时，不同拓扑荷数的LG光束其相位稳定性的差距是在减小的。而拓扑荷数与振幅稳定性则呈现出负相关的关系，尤其是拓扑荷数大于2的LG光束，其振幅稳定性远远小于L1、L2光束，此时有且只有L1光束的振幅稳定性优于高斯光束，即L1光束无论在振幅稳定性还是相位稳定性方面，其在小尺度流场内的表现都优于高斯光束；当攻角在0°到15°之间变化时，拓扑荷数小于3的LG光束振幅稳定性几乎不受影响。通过对不同拓扑荷数的LG光束与高斯光束经过小尺度流场后的远场成像质量的比较，对涡旋光束在气动环境下的气动光学效应探究提供了方法上的借鉴。