高速飞行器热流固耦合光传输分析论文

Word第第页高速飞行器热流固耦合光传输分析论文高速飞行器在大气层内高速飞行时，其光学头罩与来流互相作用产生大量的热，形成严峻的气动光学效应。一方面，光学头罩外空气受到突然压缩和摩擦，其温度、压力、密度和化学成分发生改变，形成高速冗杂流场;另一方面，光学头罩的温度和应变量会随着飞行时间发生改变，进而引起晶体材料折射率的转变。来自目标的光线到达探测器将受到流场与光学头罩的双重影响。出射波面发生畸变，对成像探测系统造成光传输干扰，使目标图像消失模糊、抖动等现象，严峻影响探测结果。因此，除了需事先对红外成像系统的性能进行测试，还需对高速飞行器气动光学效应进行全面的讨论，为系统的矫正供应数据基础。

目前，对气动环境下的光传输影响分析主要是针对高速流场或光学头罩其中一个方面进行的，对其耦合光传输的影响分析较少，本文针对这一状况，通过建立气动环境下超群声速飞行器在飞行过程中流场与头罩的耦合参数改变模型，计算CFD网格数据，采纳光线追迹计算热流固耦合光传输影响。

1飞行过程的热流固耦合分析

1.1绕流流场的数值模拟

湍流的形成是由于流体中的惯性力对流体的影响占到主导地位。流体的流淌受到物理守恒定律的支配。由于流场结构冗杂，无法用统一的方程来描述，为便利讨论，建立CFD网格，在流场中任取正六面体微元作为讨论对象，该微元和四周流体之间存在互相作用力以及质量和能量的交换。基于所取微元的任意性，对该微元建立可代表流体的运动方程，称为Navier-Stokes方程组:

1.2窗口的热流固耦合分析

高速飞行器在飞行时与四周的空气猛烈摩擦，动能转化为热能，导致流场的温度急剧上升，并加热头罩。在热流固耦合数值模拟中，绕流流场作用于头罩，使后者产生温度和应力应变的动态改变。因此，光学头罩在气动环境下的加热为持续的非稳态过程，在计算过程中应将流场与头罩进行直接耦合计算。

1.2.1外表压力场的计算依据流体力学模型可以假设，当气体质点与头罩外表碰撞后，气体沿头罩外表的法向动量全部损失掉，形成施加于物体上的力，通过计算动量的改变，可以计算出作用在头罩外表的压力:

1.2.2温度场分析

为保证流场与头罩之间的气动对流换热，头罩内部传热主要以传导和辐射方式进行。在直角坐标系中，温度分布的基本方程为

2耦合光传输过程

2.1折射率场的计算

光线通过非匀称的折射率场会发生偏折，并产生附加相位，影响光线传输。精确计算折射率场对光线追迹计算光传输效应具有重要的意义。依据之前对绕流流场数值模拟可以计算出流场的压力、密度等信息，通过物质原子理论和Gladston-Dale定律计算流场的折射率场分布:

气动热环境下的光学窗口由于受到非匀称的温度场和应力场影响，材料的密度和极化率发生转变，进而引起折射率的改变。主要表现为热光效应和弹光效应。式(10)中的其次项说明光学材料的折射率是随着温度的改变而产生的改变量，此过程称为热光效应。第三项说明由应力应变的作用引起窗口材料的光学性质转变，主要是指折射率发生改变，称为弹光效应。转变后的折射率表示为:

在热流固耦合条件下，不同飞行时间下z=0处的对称截面上流场与头罩的折射率场的分布如下图。

2.2光线追迹法

由于流场与头罩的折射率场状况冗杂并且分布无规律，本文通过数值方法求解光线方程实现光线追迹:

为保证计算精度和运算速度满意工程需要，本文采纳四阶Runge-Kutta法计算此光线方程。通过引入外推参量，将方程转化为一阶方程组，并以三维表示如下:

2.3交界面数据处理

基于应力场与温度场的存在，头罩在飞行过程中发生形变，对头罩外表进行光线追迹时采纳最小二乘法对发生形变后的外表网格进行曲面拟合。由于流体的可压缩性大，没有固定的样子，不能用计算固体材料应变理论计算流体的样子分布，导致近头罩处的流场与头罩膨胀部分的坐标数据消失交叠，这在计算耦合光传输时是不合理的，因此，需要对该部分的CFD网格数据做进一步计算处理。

如今从理论上解决该问题的方法分为两种:一是动网格技术，二是依据仿真结果对数据做插值拟合处理。动网格技术可用来模拟流体的样子随时间和边界运动而发生转变的状况，更加真实地模拟出物体运动过程的改变，目前被用于鱼雷等低速武器上。对于高速和超群声速飞行器的仿真，由于缺乏有力的运动模型支持，迭代仿真参数难以收敛。本文采纳其次种方法，通过数值手段，对交叠的网格数据进行插值，使计算后的流场能够紧密包裹住光学头罩。

2.4交界面光传输处理

当光线在两种折射率不同的.透亮介质中传播时，依据入射状况不同，会发生折射、反射、全反射等现象。除此之外，在气动热环境下，光学头罩会发生不行忽视的形变，因此，精确计算光线在流场与头罩交界面处的传播状况是分析耦合光传输效应的重要内容。

2.5像质评价指标

通过合理的光学质量评价指标可以直观的描述非匀称折射率场对光束质量的影响，为图像畸变校正供应根据，本文通过选择光学传递函数、能量集中度以及斯特列尔比等三项指标来进行像质评价，详细计算方法如下:

3仿真结果分析

设来流的静压为101325Pa，静温为300K。飞行器采纳半球形硫化锌头罩，头罩厚度为5mm，飞行高度为20km，飞行速度分别为3Ma和5Ma，攻角都为0°。其余的光学特性参数如所示。利用光线追迹计算得到不同时刻的耦合光学传递函数和能量集中度状况如所示。由可以看出，飞行时间较短时，系统的传递函数较好，此时间学系统与抱负状况比较接近，失真程度较小。当飞行速度到达5Ma，飞行时间为120s时，4Hz的目标传递函数已经衰减到接近10%，系统几乎失效。

而从能量集中度可以看出，随着飞行速度的增大以及飞行时间变长，到达90%能量时的弥散半径越来越大，系统失真程度增大。

不同飞行条件下的Strehl比同样反映图像失真状况，如所示，随着飞行条件越苛刻，其比值越小，能量分布越分散，探测器接受到的图像成像质量越低。

4结论

本文分析了高速飞行器在飞行过程中流场与头罩的热流固耦合状况，计算相应条件下的折射率场分布状况。通过该耦合折射率场光程差状况得出相应的光学传递函数、能量集中度以及斯特列尔比，依据计算结果可知:

随着飞行速度增大，耦合折射率场对光