可压缩混合层光学传输效应理论分析与实验研究

利用量级分析和风洞实验研究了可压缩混合层流动第二发展阶段气动光学效应的规律性.理论分析主要针对二维大尺度结构存在时视线误差(boresight error,BSE)与混合层流场及其特征参数之间的关系进行了讨论.研究结果表明:在混合层发展的第二阶段,时均BSE与对流马赫数呈现出复杂的非线性关系;同时还发现流场中的湍动能和混合层界面处的雷诺应力分布也是影响时均BSE的重要因素;采用细光束穿越混合层流场的风洞试验结果主要证实了时均BSE与对流马赫数之间的非线性关系.     

可压缩混合层气动光学效应研究

建立具有普适性的流场光学传输效应计算方法,针对脉动流场的光学传输计算方法可以预测空间非均匀分布的脉动流场引起的光学传输效应.以高阶紧致差分格式求解三维可压缩滤波N-S方程,对可压缩混合层流动进行空间模式的大涡数值模拟,得到流动失稳转捩直至完全湍流的空间发展全过程,采用时间模式的直接数值模拟进行对比验证.分析流动不同阶段对光学传输效应的影响,初步讨论流场特性与光学传输效应之间的联系,计算结果表明,在转捩区域流场引起的成像畸变最为严重.     

不同对流马赫数下光束畸变的实验研究

利用双曝光CCD相机成像技术,研究了小直径光束穿越可压缩混合层流场后所引起的气动光学效应.实验研究结果表明:在不同对流马赫数(Mc=0.17,0.45)下,光束投影(点扩散函数,PSF)发生了不同程度的变形和偏移.此外在不同对流马赫数下细光束投影都存在扩束、缩束和面积几乎不变的现象,但光束投影出现扩、缩束的概率不同.     

混合层流场中涡结构对流速度的特性

基于大涡模拟和光线追踪方法, 对光线穿越流场后的光程分布与混合层流场中涡结构之间的关系进行了分析, 提出了一种基于涡核位置提取的涡结构瞬时对流速度定量计算方法, 并使用直接几何测量数据进行了验证. 通过对不同尺寸的涡结构、涡-涡配对及融合过程中的涡结构和强压缩性流场中涡结构瞬时对流速度的定量数值计算, 揭示了混合层流场中涡结构对流速度的特性: 对单个涡结构而言, 其瞬时对流速度具有脉动特性, 且脉动幅度随涡结构尺寸和流场压缩性而变化; 在涡-涡配对及融合过程中, 涡对中各个涡结构的瞬时对流速度都表现出类似正弦波动的特点. 针对混合层流场中涡结构对流速度的特性, 给出了其背后的物理原因.     

光学窗口凹腔流场的光学传输效应研究

当带有光学成像探测制导系统的高速飞行器在大气层内飞行时,光学窗口与来流之间会形成复杂的凹腔绕流流场,产生气动光学效应。建立了求解超声速流场的高精度LES/RANS混合算法模型,研究了超声速流动条件下的Settles三维凹腔流动;在计算得到精确流场数据的基础上,研究了凹腔剪切层区域的光学传输效应。结果表明,凹腔流动的剪切层结构将引起光线抖动,导致严重的波面畸变,明显降低光强斯特雷尔(Strehl)比,严重影响光学传输性能,进而大大降低光学制导精度。     

超高声速飞行器光学窗口气动光学效应分析

超高声速飞行器在大气中飞行时,由于气动热和气动力的作用,光学窗口会产生严重的气动光学效应,使目标图像发生像偏移、抖动、模糊和能量衰减。利用有限元分析方法对光学窗口的热光效应、弹光效应和窗口热变形进行研究,计算了由热光学效应、光学窗口热变形引起的点扩散函数峰值大小及峰值位置(像偏移)随时间的变化趋势。根据对温度场的分布分析,计算了温度梯度对透过率的影响以及透过率与窗口热辐射随时间变化的趋势,可为光学窗口的设计、材料的选择及后期的图像处理提供依据。     

超音速湍流脉动流场气动光学效应分析

对高超音速飞行器在大气中飞行时所产生的湍流脉动气动光学效应进行了理论分析与计算.根据CFD计算流场时所使用的湍流模型及其控制方程,推导出流场的折射率脉动方差控制方程.用统计方法,求出该脉动流场的系综平均光学传递函数及相位均方差,从不同角度描述了湍流脉动的气动光学效应.计算结果表明,气动光学传输函数的幅度响应函数具有低通特征,使所获得的图像发生像模糊,而相位响应函数则导致红外成像相位非线性偏移.此外,在相同飞行高度下,马赫数越高,图像模糊越严重.电弧风洞实验结果验证了本文理论分析的正确性.     

超音速湍流脉动流场气动光学效应分析

对高超音速飞行器在大气中飞行时所产生的湍流脉动气动光学效应进行了理论分析与计算.根据CFD计算流场时所使用的湍流模型及其控制方程,推导出流场的折射率脉动方差控制方程.用统计方法,求出该脉动流场的系综平均光学传递函数及相位均方差,从不同角度描述了湍流脉动的气动光学效应.计算结果表明,气动光学传输函数的幅度响应函数具有低通特征,使所获得的图像发生像模糊,而相位响应函数则导致红外成像相位非线性偏移.此外,在相同飞行高度下,马赫数越高,图像模糊越严重.电弧风洞实验结果验证了本文理论分析的正确性.     

气动光学平均流场效应时间特性研究

根据气动光学效应定义了绕流流场厚度,分析了气动光学平均流场效应的时间特性与绕流流场厚度的关系,给出特征时间的粗估公式。使用大型流体计算软件SAED模拟计算了机载光学平台绕流流场随时间的变化,计算了绕流流场厚度,研究了气动光学平均场效应随时间的变化规律。计算结果验证了粗估公式的可靠性,表明与湍流效应的特征时间相似,气动光学平均流场效应的特征时间,正比于绕流流场厚度与飞机速度的比值。在给定的计算模型中,特征时间约为5 ms量级,特征频率低于300 Hz;平均流场气动光学效应引起的光程差均方根随时间围绕平均值变化,方差与平均值之比约为8%。     

高速光学头罩气动光学效应研究进展

伴随着飞行器飞行速度的不断增加,高速光学头罩技术具有极大的应用潜力和前景.只是,高速光学头罩在大气层中飞行,受到气动光学效应的影响,成像目标会出现畸变、抖动、模糊和能量衰减,影响成像精度.通过对数十年来高速光学头罩流体力学和气动光学效应的总结和梳理,初步构建了高速光学头罩气动光学效应相似准则模型.归纳了气动光学效应抑制的主要技术路线以及目前已经开展的一些工作.最后总结和讨论了高速光学头罩气动光学效应关键技术难点及未来的发展趋势,可为今后的气动光学效应研究提供一些参考和帮助.      

Experimental evidence for non-linear growth in compressible mixing layer

An experimental study of compressible mixing layers(CMLs)was conducted using planar laser Mie scattering(PLMS)visualizations from condensed ethanol droplets in the flow.Large ensembles of digital images were collected for two flow conditions at convective Mach numbers Mc=0.11 and 0.47.The coherent vortices,braids and eruptions in the mixing zone were observed,interpreted as evidence of multi-scale,three-dimensional structures at a high Reynolds number.The mixing layers with a large visualized range present two stages along the streamwise direction,corresponding to the initial mixing and the well-developed stage.A new method,the gray level ensemble average method(GLEAM),by virtue of the similarity of the mixing layer,was applied to measure the growth rate of the CML thickness.New evidence for a nonlinear growth of CML is reported,providing an interpretation of previous observations of the scattering of the growth rate.     

An experimental study of aero-optical aberration and dithering of supersonic mixing layer via BOS

The optical performance of supersonic mixing layer is heavily deteriorated by the aero-optical aberration and dithering of coherent structures, but current measuring methods limit the spatiotemporal resolution in relevant studies. A high resolution whole-field aero-optical aberration and dithering measuring method based on the Background Orient Schlieren (BOS) technique was studied. The systematic structure, sensitivity and resolution of BOS are analyzed in this paper. The aero-optical aberration and dither...     

The fractal measurement of experimental images of supersonic turbulent mixing layer

Flow visualization of supersonic mixing layer has been studied based on the high spatiotemporal resolution Nano-based Planar Laser Scattering(NPLS) method in SML-1 wind tunnel. The corresponding images distinctly reproduced the flow structure of laminar,transitional and turbulent region,with which the fractal measurement can be implemented. Two methods of measuring fractal dimension were introduced and compared. The fractal dimension of the transitional region and the fully developing turbulence region of supersonic mixing layer were measured based on the box-counting method. In the transitional region,the fractal dimension will increase with turbulent intensity. In the fully developing turbulent region,the fractal dimension will not vary apparently for different flow structures,which em-bodies the self-similarity of supersonic turbulence.     

DNS of compressible turbulent boundary layer over a blunt wedge

Wall turbulence is more complicated than free turbulence, and the direct numerical simulation (DNS) of wall turbulence is more difficult. In recent years, most of DNS cases for wall turbulence are simplified by using temporal mode, where streamwise pe- riodic boundary condition is imposed. In temporal mode, spatial transition will be con- sidered as an analogue of time-evolving transition. For the channel turbulence, an equivalent body force can substitute for the mean gradient of pressure, …     

Interaction between strain and vorticity in compressible turbulent boundary layer

The interaction of strain and vorticity in compressible turbulent boundary layers at Mach number 2.0 and 4.9 is studied by direct numerical simulation(DNS)of the compressible Navier-Stokes equations.Some fundamental characteristics have been studied for both the enstrophy producing and destroying regions.It is found that large enstrophy production is associated with high dissipation and high enstrophy,while large enstrophy destruction with moderate ones.The enstrophy production and destruction are also correlated with the dissipation production and destruction.Moreover,the enstrophy producing region has a distinct tendency to be‘sheet-like’structures and the enstrophy destroying region tends to be‘tube-like’in the inner layer.Correspondingly,the tendency to be‘sheet-like’or‘tube-like’structures is no longer obvious in the outer layer.Further,the alignment between the vorticity vector and the strain-rate eigenvector is analyzed in the flow topologies.It is noticed that the enstrophy production rate depends mainly on the alignment between the vorticity vector and the intermediate eigenvector in the inner layer,and the enstrophy production(destruction)mainly on the alignment between the vorticity vector and the extensive(compressive)eigenvector in the outer layer.     

Evaluation of a flamelet/progress variable approach for pulverized coal combustion in a turbulent mixing layer

A steady flamelet/progress variable (FPV) approach for pulverized coal flames is employed to simulate coal particle burning in a turbulent shear and mixing layer. The configuration consists of a carrier-gas stream of air laden with coal particles that mixes with an oxidizer stream of hot products from lean combustion. Carrier-phase DNS (CP-DNS) are performed, where the turbulent flow field is fully resolved, whereas the coal is represented by Lagrangian point particles. CP-DNS with direct chemistry integration is performed first and provides state-of-the-art validation data for FPV modeling. In a second step the control variables for FPV are extracted from the CP-DNS and used to test if the tabulated manifold can correctly describe the reacting flow (a priorianalysis). Finally a fully coupled a posteriori FPV simulation is performed, where only the FPV control variables are transported, and the chemical state is retrieved from the table and fed back to the flow solver. The a priori results show that the FPV approach is suitable for modeling the complex reacting multiphase flow considered here. The a posteriori data is similarly in good agreement with the reference CP-DNS, although stronger deviations than a priori can be observed. These discrepancies mainly appear in the upper flame (of the present DNS), where premixing and highly unsteady extinction and re-ignition effects play a role, which are difficult to capture by steady non-premixed FPV modeling. However, the present FPV model accurately captures the lower, more stable flame that burns in non-premixed mode.