基于CFD/RBD方法的旋转弹轨迹仿真

旋转弹在飞行过程中,由于表面气动载荷分布不对称,可能会出现耦合共振、灾难偏航等现象,使导弹飞行失常偏离目标。针对此类问题,采用非定常N-S方程和刚体动力学方程耦合求解的方法,仿真旋转弹飞行轨迹。非定常气动力计算使用非结构混合网格N-S方程求解程序HUNS3D,刚体六自由度运动方程采用改进的四阶Adams预估校正法求解。将这种耦合求解方法应用于美国ARL旋转弹试验模型的轨迹仿真,并将计算得到的旋转弹姿态角、空间位置和气动力与相关文献进行对比分析,结果表明:本文方法能较精确地仿真旋转弹的飞行轨迹。     

固液混合火箭发动机喷焰红外辐射特性分析

为了研究固液混合火箭发动机喷焰在不同飞行状态的红外辐射特性,建立了固液两相状态下喷焰的红外辐射计算模型,获得了喷焰在2~12μm波段的红外辐射.将喷焰的温度、压力、摩尔组分等作为输入条件,米氏散射(Mie)模型计算颗粒相的散射,采用有限体积法(FVM)离散辐射传递方程,利用中分辨率大气辐射传输(MODTRAN)模型计算不同观测距离的大气透过率,编程计算得到不同飞行高度,不同观测方向下喷焰的红外辐射.结果表明喷焰红外辐射亮度与温度分布类似,在2.7~3.0μm和4.2~4.5μm波段范围内的辐射较强,高温粒子在近红外比远红外对总辐射的贡献大,太空与地面观测不同飞行高度的发动机喷焰辐射率因大气衰减的不同而有所差异,可为目标跟踪与发动机设计提供参考.      

空间目标红外辐射特性分析

研究空间目标的红外辐射特性,对于探测识别空间目标、监视空间碎片、保障卫星运行安全等问题有重要意义.本文以空间卫星为目标,分析了其红外辐射特性形成的原理,依据空间目标的热平衡方程,应用节点网络法求解目标表面的温度分布,并结合普朗克定律推导空间目标红外辐射度的计算方法,在分析总结相关资料的基础上得到了求解空间目标红外辐射特...     

飞行器表面红外辐射的模型以及红外特征分析(英文)

隐身最重要的技术之一是评估飞行器被红外探测器接收到的红外辐射强度。本文针对飞行器表面红外辐射特征计算的问题提出一个方法。首先通过分析飞行器表面气动加热、发动机热部件等主要热源,建立了CFD计算模型并获得了飞行表面的温度场。并结合反向蒙特卡洛方法建立求解飞行器表面复杂几何外形的红外辐射模型。通过数值模拟,比较了飞行器主要部件在水平和侧向上红外辐射强度所占的比重;并针对飞行器主要部件冷却对红外辐射强度的影响效果进行了研究。计算结果表明:在水平方向降低飞行器头部、垂尾、机翼部件表面温度10K能降低红外辐射强度大于8%;而在侧方向降低飞行器头部、垂尾、机翼部件表面温度10K能降低红外辐射强度小于8%。计算结果可为隐身设计提供参考。     

矢量喷管内燃气辐射与壁面温度的耦合计算

选取内、外调节片和隔热屏建立几何模型,基于封闭腔净辐射模型和壁面热平衡模型建立了燃气辐射与喷管壁面温度的耦合算法.波段为1~5μm的气体辐射采用窄波带模型计算,其他波段不考虑气体辐射,建立辐射净热流密度-有效辐射亮度-壁面温度的关联式求解燃气与壁面的辐射换热,采用牛顿-拉斐尔森迭代法求解壁面热平衡方程计算其温度.对某轴对称矢量喷管(偏转20°),计算了喷管壁面的红外光谱辐射和辐射净热流密度,以及各部分结构的温度.作为验证,还计算了文献中某液体火箭发动机轴对称矢量喷管壁面的辐射净热流密度,与文献的结果进行对比一致性较好.研究表明:轴对称矢量喷管偏转段沿周向的辐射热流密度和温度差异很大,沿偏转方向部位壁面的温度和辐射热流密度都较低,偏转方向壁面的温度比相反方向大约低10%,辐射热流密度大约低50%.      

某型飞机短舱辐射换热计算

针对某飞机地面慢车工况,进行了短舱内复合换热数值计算研究,分别运用球谐函数法、离散坐标法和蒙特卡洛法计算辐射换热.结果显示:不同计算方法所得结果总体上趋于一致,热辐射与自然对流耦合作用明显.与单纯自然对流的计算结果相比,舱内温度显著升高,各区域温度值上升7%~30%.考虑空气介质对不同波段热辐射的选择性吸收后,短舱内各计算点温度进一步上升.综合分析表明:地面条件下飞机短舱换热分析过程中热辐射的影响不可忽略.      

用温度修正的k-ε模型研究排气系统红外辐射特性

为研究温度修正的k-ε湍流模型对排气系统红外辐射特性计算精度的改善程度,计算了某轴对称收敛喷管在3~5μm波段的空腔-喷流组合红外辐射特性,并与试验结果对比.喷管流场及温度场采用有限体积法(FVM)求解Navier-Stokes(N-S)方程得到,湍流模型分别采用标准k-ε模型和温度修正的k-ε模型.红外辐射特性的计算采用有限体积法求解吸收-发射性介质条件下的三维辐射传输方程,并考虑大气衰减作用.计算结果表明:温度修正的k-ε模型可有效预测排气系统热喷流的核心区长度,从而获得与试验结果一致的红外辐射强度.与标准k-ε模型相比,采用温度修正的k-ε模型计算得到的红外辐射强度的最大相对误差由64%降至了4.3%,说明温度修正的k-ε模型适用于精确模拟排气系统的红外辐射特性.      

两种涡扇发动机排气系统红外辐射特性的比较

针对弹用涡扇发动机，计算和比较了3μm-5μm波段轴对称分开排气喷管和轴对称混合排气喷管的红外辐射特性，并与模型实验的结果进行了比较。结果表明，两种类型的排气喷管在红外辐射的空间分布和光谱分布特性方面具有相似的特征，但是红外辐射强度的数值有较大差别；内、外涵混合排气，可以显著降低喷管空腔、尾喷流和整个排气系统的红外辐射。     

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

十一组元高温平衡、非平衡空气光谱吸收系数计算

飞行器以很高的马赫数再入大气层时 ,头部激波层、尾迹的气体辐射产生的紫外、可见、红外特征信号 ,是地面监测和反导制导系统探测、识别的主要依据 ;同时 ,头部气体层热辐射也是再入体壁面加热热流的重要源项。激波层内高温气体吸收系数是目标紫外、可见光、红外辐射特性计算、热防护计算和气动流场辐射场耦合计算的基本参数。本文分析了再入过程中高超声速稀薄气体流的非平衡现象 ,采用三温度模型 (电子能温度 Te、振动能温度 TV、转动和平动能温度 TR)来表征激波层内非平衡态的气体组份各个能级占有数分布 ,并由原子分子辐射理论直接计算 N2 ,N+ 2 ,N,N+ ,O2 ,O+ 2 ,O,O+ ,NO,NO+和 e-等 1 1种主要空气组元各种能级跃迁对辐射的贡献 ,最终得到平稀和非平衡高温空气吸收系数计算模型。计算结果可见光区与试验结果符合较好 ,紫外、红外区结果稍差     

高速飞行器瞬态气动热试验模拟系统

为得到高速飞行器表面各部分的热应力，应变，结构膨胀量等高温力学性能参数，须建立高速飞行器瞬态气动热试验模拟系统。针对高速飞行器气动模拟试验瞬态热控过程所具有变化复杂，高度非线性，瞬变，强耦合的特点，将模拟控制方法应用于瞬态气动热模拟控制系统。该系统能够按照高速飞行器飞行过程中表面热流和温度的瞬态连续变化对气动模拟加热过程实施快速，准确的动态控制。     

基于再入轨迹和气动热环境的返回舱烧蚀研究简

 针对再入全过程合理预测热防护罩表面材料烧蚀深度和温度的动态变化问题,提出融合再入轨迹、气动热以及Newton-Raphson和三对角矩阵算法（TDMA）构建动态烧蚀的方法。该方法建立直入式和跳跃式三自由度再入轨迹,应用修正的牛顿流体理论估算气动参数,以及修正的Fay-Riddell和Sutton-Grave理论计算驻点区域的热流密度,利用一维非线性热传导方程模拟了热防护材料的烧蚀过程。仿真结果表明：此方法实现了再入全过程热防护材料烧蚀深度和温度连续动态变化的预测,同样适用于更为复杂结构飞行器的动态烧蚀预测,与热平衡积分法（HBI）相比其结果可靠合理,为进一步优化热防护系统（TPS）提供了一定的参考依据。     

基于再入轨迹和气动热环境的返回舱烧蚀研究

针对再入全过程合理预测热防护罩表面材料烧蚀深度和温度的动态变化问题,提出融合再入轨迹、气动热以及Newton-Raphson和三对角矩阵算法（TDMA）构建动态烧蚀的方法。该方法建立直入式和跳跃式三自由度再入轨迹,应用修正的牛顿流体理论估算气动参数,以及修正的Fay-Riddell和Sutton-Grave理论计算驻点区域的热流密度,利用一维非线性热传导方程模拟了热防护材料的烧蚀过程。仿真结果表明：此方法实现了再入全过程热防护材料烧蚀深度和温度连续动态变化的预测,同样适用于更为复杂结构飞行器的动态烧蚀预测,与热平衡积分法（HBI）相比其结果可靠合理,为进一步优化热防护系统（TPS）提供了一定的参考依据。     

气动加热下光学窗口传热特性的数值模拟

针对气动加热条件下光学窗口辐射导热耦合传热开展研究,以了解光学窗口的热辐射发射规律及热防护性能。采用有限体积法、蒙特卡罗法、谱带模型建立了具有镜反射界面的光学窗口辐射传输及辐射导热耦合传热计算模型,对具有两层SiO2玻璃的光学窗口传热特性进行了模拟。结果表明处于高温区玻璃向外发射热辐射,低温区玻璃吸收热辐射;除了加热面附近外的玻璃区域,辐射热流密度远大于导热热流密度,在2.9～4.2μm谱带内的辐射热流密度大于0.4～2.9μm谱带内的;在0.4～2.9μm谱带内,可见光波长内的出射辐射热流密度最小,中红外波长范围内出射辐射热流密度最大。在光学窗口总厚度一定的情况下,降低玻璃层厚度和玻璃夹层气体热导率能够有效降低光学窗口非加热面稳态温度。     

典型尖楔结构燃气流热试验工况确定方法

在传统冷壁热流模拟方法的基础上,进一步提出以热壁温度及热流密度的时序变化曲线为控制目标的燃气流热试验工况确定方法,即利用壁温控制目标与实测值的偏差对热壁热流控制目标做一定修正,以尽可能消除和弥补前期试验误差,同时利用300K冷壁边界热流密度数据库插值迭代方法,快速确定一定气动热模拟所需燃气流温度,解决了沿飞行轨迹瞬态热试验技术难题之一。利用CFD数值模拟方法,建立了典型尖楔结构高/中温双路燃气流组合热试验300K冷壁边界热流密度数据库,并针对典型尖楔结构沿某飞行轨迹9个典型状态气动热模拟需求,确定相应双路燃气流热模拟参数。相关数值计算结果显示,驻点区域热流密度平均模拟偏差为4.5%,平板区热流密度平均模拟偏差为4.6%,两者最大模拟偏差均不大于8%,满足工程试验精度要求。同时,瞬态热分析结果显示第45s时,距驻点1mm处最大温度梯度达到21K/mm,距驻点10.1mm处最大温度梯度达到18K/mm,满足气动热大温度梯度效应需求。      

高超声速典型弹道下的壁板热气动弹性动力学分析

对高超声速流中带有热防护系统（TPS）的二维壁板进行了热气动弹性的双向耦合建模与分析,采用三阶活塞理论计算气动力,通过参考焓法获得气动热流,在有限差分法的基础上进行结构热传导计算,拟合了结构材料特性随温度退化的曲线,最后将气动热模块、气动弹性模块进行双向耦合以考虑气动热与结构形变之间的相互反馈,并在2种典型弹道状态下进行热气动弹性响应分析。结果表明,在X-34A的设计弹道下,双向耦合分析会引起更加剧烈的热应力与热弯矩的变化与较长的瞬态混沌过程。在FALCON弹道下,双向耦合得到的结果加热更为剧烈,而温度下降的过程也更快。对比2种弹道发现,长时间的高超声速飞行更容易引发颤振,而机动性较强的弹道面临的主要问题则是屈曲,需要考虑材料的应力及强度特性。同时说明了双向耦合策略对于现代飞行器在弹道状态下工作的热气弹响应分析的必要性。     

“Ⅳ型”激波干扰中流-热-固耦合问题一体化计算分析

针对高超声速进气道前缘"Ⅳ型"激波干扰产生的气动加热与结构传热多物理场耦合计算问题，发展了一种基于有限体积法的流-热-固一体化计算方法。该方法采用一体化控制方程组统一离散求解外部高速流场与内部结构温度场，规避了传统分区耦合算法在时间域内交替迭代的繁琐数据交换策略。另外，提出一种新的双温阻模型计算流-固交界面的物性参数以保证计算准确性，采用LU-SGS隐式时间迭代和自适应时间步长以提高计算效率。采用经典高超声速二维圆管流-热-固耦合算例对该一体化方法进行验证，计算结果与试验值和参考文献数据吻合较好，证明了该方法的可靠性和正确性。利用一体化方法对高超声速前缘"Ⅳ型"激波干扰流-热-固耦合问题进行定常/非定常计算与分析，给出了温度与热流的时变特性，计算结果表明，激波干扰作用产生的超声速"喷流"不断冲击壁面，使得壁面最大压力系数增大约9倍，壁面最大热流增大约4.7倍，给高速飞行器的热防护设计与选材带来严峻挑战。同时，也表明了一体化计算方法可以较好地用于长航时飞行条件下与复杂飞行环境下的高超声速热防护系统的热环境特性分析与综合性能评估。     

固定翼双旋弹动力学分岔特性分析

针对一种滚转稳定的固定翼双旋弹,对其非线性动力学进行了分岔特性分析,并在此基础上研究了各系统参数对其动力学分岔特性的影响。根据固定翼双旋弹非对称的特点,通过数值计算方法研究其飞行过程中平衡点随同向鸭翼安装角的变化规律,通过系统的分岔图得知系统具有三组稳定平衡点,其中只有一组平衡点为理想可行的稳定平衡点,因此需限定同向鸭翼安装角的范围以使固定翼双旋弹保持稳定飞行。在此基础上针对固定翼双旋弹弹道修正组件周期旋转和转角固定两种工作模式,通过各系统参数下的系统分岔图总结了固定翼双旋弹结构及气动力参数对其动力学系统分岔特性的影响。仿真结果表明,固定翼双旋弹的各气动力参数及飞行速度均对系统的分岔特性具有较大影响,应合理选定这些系统参数以使其具有良好的气动特性。