高超声速飞行器气体辐射噪声计算方法

高超声速飞行器周围的激波层内高温气体会发生剧烈的物理化学变化,伴随强烈的光辐射过程,直接影响红外导引头的光学成像效果。采用流体力学Navier-Stokes方程和热化学非平衡模型模拟高温非平衡流动,考虑电子跃迁和振转跃迁以窄带法求解气体辐射特性参数,基于有限体积法离散辐射传输方程,在分别验证流场解算、辐射参数求解和辐射传输计算的基础上,进行了飞行器高速飞行的流动和辐射模拟。数值求解得到了飞行器流场特征和粒子数空间分布。计算的选定波长范围内的气体辐射发射系数空间分布显示其与激波形状和波后气体温度分布相似。通过传输得到的飞行器光学窗口视线路径上的气体辐射噪声成轴对称分布,发现辐射噪声和飞行速度、气体成分等密切相关,马赫数增加时气体辐射噪声显著增强。     

高超声速飞行器气体辐射噪声计算方法

高超声速飞行器周围的激波层内高温气体会发生剧烈的物理化学变化,伴随强烈的光辐射过程,直接影响红外导引头的光学成像效果。采用流体力学Navier-Stokes方程和热化学非平衡模型模拟高温非平衡流动,考虑电子跃迁和振转跃迁以窄带法求解气体辐射特性参数,基于有限体积法离散辐射传输方程,在分别验证流场解算、辐射参数求解和辐射传输计算的基础上,进行了飞行器高速飞行的流动和辐射模拟。数值求解得到了飞行器流场特征和粒子数空间分布。计算的选定波长范围内的气体辐射发射系数空间分布显示其与激波形状和波后气体温度分布相似。通过传输得到的飞行器光学窗口视线路径上的气体辐射噪声成轴对称分布,发现辐射噪声和飞行速度、气体成分等密切相关,马赫数增加时气体辐射噪声显著增强。     

高超声速飞行目标尾焰红外辐射测量

为了研究超高声速飞行器发动机尾焰喷射高温高速气流的辐射特性,对尾焰成分CO₂及H₂O分子在4.3和2.7 μm大气窗口红外辐射波段进行了测量.利用高温燃气激波风洞模拟产生超高声速飞行器尾焰喷流,喷流速度M=5.5.实验中选用单元型InSb红外探测器,并利用黑体进行原位定标.测量距离为0.7 m,采用单透镜成像加光阑的方法收集光信号.实验中分别沿喷流方向喷流垂直方向进行了多点测量,通过定标结果反演得到尾焰在4.3和2.7 μm分别沿喷流方向和喷流垂直方向的光谱辐亮度和波段辐亮度分布,测量结果表明4.3 μm辐射强度及稳定性均高于2.7 μm.      

一种靶机红外辐射特性增强技术研究

为了在外场试验中检验红外制导武器系统的工作性能,需要使用靶机红外增强装置实现对真实目标红外辐射特性的模拟。介绍了一种利用光学方法实现靶机红外辐射增强和目标辐射特性模拟的技术。使用高亮度红外光源作辐射源,通过对光学系统的控制,在实现红外辐射增强的同时,可模拟不同距离和方向的红外辐射强度特性。通过建模对所模拟目标的红外辐射强度及辐射方向特性进行了分析,并以此为设计依据设计了光学系统。实现了对最大红外辐射强度为1600W/sr的目标辐射方向特性的模拟。     

高温二氧化碳气体红外辐射实验研究

主要介绍了依托高频等离子体风洞建立的高温气体辐射测量平台,并在此平台上开展了高温二氧化碳气体红外辐射实验测量。介绍了高频等离子体风洞的运行原理、流场特性及工作介质;介绍了实验测量的条件、装置、标定、数据处理方法和结果分析;通过自建的高温气体发射光谱测量平台实验测量了二氧化碳气体在1 500～3 000 K范围内4个温度点的红外发射光谱;介绍了Abel变换在测量二氧化碳气体红外辐射空间分布中的应用,通过Abel变换获得了高温下二氧化碳气体红外辐射的空间分布结果; 分析了高温二氧化碳气体在4.3 μm附近的红外辐射的强度及其中心波长随温度变化的分布,得到了发射峰中心波长随温度的升高向长波方向展开的结果,并与文献结果进行了对比分析。     

三维高超声速无粘定常绕流的数值模拟

本文采用一种简单有效的通量分裂结合一种二阶TVD格式的数值通量的方法,提出一种隐式的迎风有限体积格式,并利用这种格式,从气体动力学非定常Euler方程组出发,数值模拟了三维不对称物体的高超声速无粘定常绕流。数值结果表明此格式具有分辨率较高和收敛速度较快的优点。     

飞艇红外探测系统探测高超声速目标性能研究

以X-51A为例,研究了飞艇红外探测系统对临近空间高超声速目标的探测性能。首先,根据飞行器的飞行状态和飞行高度建立了临近空间高超声速目标不同波段的红外辐射特性模型,以及随高度变化的目标背景红外辐射强度模型;其次,综合考虑飞行器与飞艇高度、地球曲率及红外辐射在大气中传播的波段选择性等因素,建立了红外辐射在临近空间大气中传播的透过率模型;在此基础上,建立了飞艇红外探测系统对高超声速目标的探测距离模型。通过仿真得到了临近空间高超声速目标在不同飞行状态下3个波段的红外辐射强度随目标飞行高度变化的曲线,以及飞艇红外探测系统对飞行器在不同飞行状态下3个红外辐射波段的探测能力。研究结果表明:飞艇红外探测系统对高超声速目标的有效探测距离可以达到百公里量级;当飞行器飞行状态一定时,随着飞行器飞行高度的增加,系统对目标的探测距离先增大后减小;与长波波段相比,中短波波段的探测距离更大,并给出了临近空间飞艇应尽量布置在海拔高度大于18 km的高空中的部署建议。     

迎头方向隐身导弹红外辐射特性分析北大核心CSCD

隐身导弹应用红外隐身技术将降低导弹被红外探测系统的发现概率。通过以AGM-158C导弹为参考,对迎头观测方向目标的红外辐射特性,包括目标蒙皮、目标蒙皮反射周围背景和目标尾焰的红外辐射强度进行估算,并与典型亚音速导弹红外辐射特性的数据进行对比分析。分析结果表明:红外隐身技术可以明显降低导弹的中、长波红外辐射,尤其是中波波段的红外辐射;红外探测系统设计时尾焰探测通常选用中波红外探测器,弹体探测通常选用长波红外探测器;隐身导弹采用红外隐身技术将提高导弹的突防能力。     

新颖的高温流场红外辐射物性参数计算方法

为解决高速飞行器红外辐射特性研究中高温气体辐射参数计算问题,在转动平滑谱带（smeared rotational band, SRB）模型的假设下,首次提出基于贝叶斯估计的跳转马尔科夫链蒙特卡罗（reversible jump Markov chain monte carlo, RJMCMC）方法来计算高温流场中气体组分红外光谱物性参数。仿真计算和实验测量表明,仅需很少的先验知识,且不需要相关跃迁系数,该方法可以很好的估计主要振动跃迁带的位置,所提取的SRB模型参数的变化规律与SRB理论模型结果一致,利用拟合参数计算得到不同温度下谱带积分发射和吸收系数与逐线计算的结果非常近似。     

高超声速飞行器前缘流固耦合计算方法研究

对高超声速流场和结构温度场进行了耦合计算分析, 同时基于准静态假设对结构应力进行了分析. 流场部分采用基于非定常Navier-Stokes (N-S)方程的有限体积法, 湍流模型采用SST k-ω 模型, 固体部分采用基于非稳态热传导方程的有限元法, 同时基于准静态假设对固体结构的应力应变进行了分析. 在流固交界面处, 高速流体从固体结构得到温度边界条件, 固体结构从高速流体得到热流边界条件, 从而实现了流场和固体温度场的紧耦合计算.通过与超声速无限长圆管绕流试验结果进行对比, 验证了该方法的可靠性. 同时对二维圆管结构在气动加热过程中的温度、应力等的变化进行了比较详细的分析. 研究结果表明: 随着气动加热时间的推进, 由于圆管结构的高温区在不断扩大, 导致了结构的热变形在不断地增大; 圆管最小变形区出现在θ为60°处; 同时研究发现在计算时间内圆管热变形对外部流场的影响可以忽略不计.     

高超声速流场成像质量特性及校正方法

光学成像制导飞行器在大气层中以高超声速飞行时，其成像窗口附近强压缩流场的气动光学效应会导致成像过程出现抖动、偏移和模糊，影响制导精度.为研究该问题，搭建了基于高超声速（Ma = 6.0）炮风洞的气动光学地面模拟平台.利用高速摄像机获取了多种喷流压比状态下光学头罩成像图片，研究了成像特性.基于背景纹影技术（background oriented schlieren，BOS）直接获取气动光学畸变的点扩散函数信息，结合Wiener滤波方法对地面模拟平台获取的成像畸变结果进行了校正，并结合灰度分布、峰值信噪比（peak signal-to-noise ratio，PSNR）和结构相似度（structural similarity，SSIM）对校正结果进行了定性和定量评价.成像结果表明，头罩无冷却喷流时成像质量最好，在压力匹配附近头罩成像质量相对于欠压喷流和过压喷流成像质量较好.图像校正结果表明，风洞运行过程中采集的时间序列图像在校正之后所对应的灰度分布情况、PSNR和SSIM都得到提高.      

A Numerical Tackling on Sakiadis Flow with Thermal Radiation

Momentum and energy laminar boundary layers of an incompressible fluid with thermal radiation about a moving plate in a quiescent ambient fluid are investigated numerically. Also, it has been underlined that the analysis of the roles of both velocity and temperature gradient at infinity is of key relevance for our results.     

隐身飞机尾焰的红外辐射特性

隐身飞机尾焰的红外辐射是隐身飞机探测的主要辐射源.本文提出了一种新的隐身飞机尾焰红外辐射特性计算模型.该模型以普通飞机尾焰红外辐射特性计算模型为基础,进而考虑红外隐身措施的影响,间接实现隐身飞机尾焰的红外辐射特性的计算.计算模型分别考虑了隐身飞机的二元喷口、引射技术、红外遮蔽云以及遮挡对尾焰辐射的影响.计算结果得出,添...     

A New Numerical Method for Solving the Acoustic Radiation Problem

A numerical method of solving the problem of acoustic wave radiation in the presence of a rigid scatterer is described. It combines the finite element method and the boundary algebraic equation one. In the proposed method, the exterior domain around the scatterer is discretized, so that there appear an infinite domain with regular discretization and a relatively small layer with irregular mesh. For the infinite regular mesh, the boundary algebraic equation method is used with spurious resonance suppression according to Burton and Miller. In the thin layer with irregular mesh, the finite element method is used. The proposed method is characterized by simple implementation, fair accuracy, and absence of spurious resonances.     

隐身飞机尾焰的红外辐射特性

隐身飞机尾焰的红外辐射是隐身飞机探测的主要辐射源.本文提出了一种新的隐身飞机尾焰红外辐射特性计算模型.该模型以普通飞机尾焰红外辐射特性计算模型为基础,进而考虑红外隐身措施的影响,间接实现隐身飞机尾焰的红外辐射特性的计算.计算模型分别考虑了隐身飞机的二元喷口、引射技术、红外遮蔽云以及遮挡对尾焰辐射的影响.计算结果得出,添加隐身措施后尾焰辐射强度仅为添加前辐射强度的5.8%.针对隐身飞机尾焰红外辐射特性很难获取的问题,将计算结果与喷灯燃烧航空煤油的光谱峰位数据进行了比较,实验结果显示隐身前后辐射能量量级变化与国外文献相同,表明该模型可以用于隐身飞机尾焰红外辐射特性计算.     

土卫六探测器高超声速化学非平衡特性数值计算

为了研究土卫六探测器在地球大气和土卫六大气中飞行时气动热环境的差异,采用数值模拟的方法分别对地球大气成分和土卫六大气成分下探测器前方流场内气动热环境进行了研究,得到如下结论:与地球大气条件相比,土卫六大气条件下探测器表面的热流密度较高,其峰值可达到地球大气条件的1倍以上.探测器表面的热流密度先随着高度的减小而增加,而后...     

A Numerical Method on Eulerian Grids for Two-Phase Compressible Flow

We develop a numerical method to simulate a two-phase compressible flow with sharp phase interface on Eulerian grids. The scheme makes use of a level set to depict the phase interface numerically. The overall scheme is basically a finite volume scheme. By approximately solving a two-phase Riemann problem on the phase interface, the normal phase interface velocity and the pressure are obtained, which is used to update the phase interface and calculate the numerical flux between the flows of two different phases. We adopt an aggregation algorithm to build cell patches around the phase interface to remove the numerical instability due to the breakdown of the CFL constraint by the cell fragments given by the phase interface depicted using the level set function. The proposed scheme can handle problems with tangential sliping on the phase interface, topological change of the phase interface and extreme contrast in material parameters in a natural way. Though the perfect conservation of the mass, momentum and energy in global is not achieved, it can be quantitatively identified in what extent the global conservation is spoiled. Some numerical examples are presented to validate the numerical method developed.     

任意四边形网格上解辐射输运方程的一种数值方法

解二维辐射输运问题时,对于矩形网格,用Sn菱形(diamond)格式处理十分简单有效,但该格式不能直接用于非矩形网格,而辐射流体力学计算往往要求在非矩形网格上进行.本文讨论任意四边形网格的辐射输运数值方法,给出一种全正的子网格平衡计算格式,并进行数值计算,与矩形网格Sn菱形格式结果比较,得出本方法可行性结论.     

Features of Hypersonic Flow Past a Yawed Plate

Doklady Physics - Flow past a yawed plate in the strong interaction mode is considered in the case when the pressure at its trailing edge varies along the transverse coordinate instead of being...