高超声速化学非平衡流动Navier-Stokes/DSMC耦合算法

基于相同的化学反应模型,在已有计算流体力学（CFD）和直接仿真蒙特卡罗（DSMC）方法及程序的基础上,采用Modular Particle-Continuum（MPC）耦合技术,建立了包含化学非平衡Navier-Stokes/DSMC耦合算法。算法结构中DSMC计算区域在CFD计算结果上根据当地克努森数自动选取。发展了适用于流场分区信息交换的亚松弛技术,抑制DSMC方法对CFD计算的影响。把DSMC方法和CFD的应用范围拓展到过渡流区,为复杂飞行器近连续过渡流区高超声速化学非平衡流动数值模拟研究提供了一种工程适用的预测分析手段。通过对二维圆柱高超声速化学非平衡绕流的算例与其他结果的比较研究,表明耦合算法不论在流场结构、流场非平衡现象,还是飞行器表面参数、整体气动力/热特性方面,都能够得到与全DSMC计算吻合的结果,证实了所建立的Navier-Stokes/DSMC耦合计算模型与方法的有效性和可靠性。仿真了某航天器解体碎片在过渡区的化学非平衡流动,得到碎片在过渡区的气动力/热特性,为碎片的陨落计算提供依据。     

短壳绝热面积对液氢贮箱绝热性能影响

基于计算流体动力学(CFD)方法研究了典型5 m直径液氢贮箱在短壳未包裹绝热材料、50%面积及100%面积包裹绝热材料3种情况下对贮箱内液氢蒸发特性的影响。数值计算基于流体体积(VOF)模型计算两相流,基于Lee模型计算气液界面传质率,考虑了短壳包裹泡沫表面及未包裹泡沫的暴露表面结霜对漏热的影响,构建的数值模型及界面传质计算具有清晰的气液界面,准确地捕捉到了液氢液面的变化。结果表明:短壳是液氢贮箱漏热的主要因素,对液氢蒸发率影响起重要作用;相对于短壳未绝热,50%绝热使得液氢贮箱气相平均温度从110 K下降到32 K,绝热面积占比增加到100%时,气相平均温度下降到约23 K,绝热改善效果相对降低;比较短壳绝热面积占比从50%增加到100%与从0增加到50%对相对蒸发率影响,前者差异较小,仅降低24%,而后者差异明显,下降了409%。研究结果指导了液氢贮箱绝热结构的优化设计。     

实验中实现"无限展长后掠翼"流场的方法研究

在风洞实验中,为了获得一个近似无限展长后掠翼的三维流场实验条件,必须采用特殊的整流措施,修正洞壁对流场的干扰。对于小展弦比模型的实验,可以通过在两侧洞壁附加衬垫来修正洞壁干扰。本文基于CFD技术,研究了采用无限展长后掠翼空间流线拟合的方法设计洞壁衬垫。数值模拟结果表明,采用流线拟合设计的衬垫能够有效地修正风洞侧壁对流场展向分布的影响,得到近似无限展长后掠翼的流场。     

旋流杯设计参数对燃烧性能的影响

通过三维数值模拟和试验相结合的方法,针对航空发动机主燃烧室开展了火焰筒头部旋流杯设计参数对燃烧室性能影响的研究,其关键设计参数为副旋流器相对进气量和旋流数,关注的燃烧性能包括出口温度场分布、NOx排放量等.结果表明:数值计算获得的规律与试验的结论基本一致.随着副旋流器相对进气量的增加,主燃烧室出口温度场逐渐恶化,而NOx排放量略有降低;副旋流器旋流数偏离基准设计值时,出口温度场恶化,而NOx排放量随旋流数的增加而增加.     

基于非惯性系的悬停状态旋翼CFD/CSD耦合气动分析

旨在提高先进旋翼气动特性的分析精度,在旋翼高精度CFD分析中耦合气动弹性效应,取代传统方法中的刚性桨叶假设,并考虑悬停状态旋翼流场准定常的特性,在非惯性坐标系下建立了一套适合于悬停状态旋翼气动特性计算的CFD/CSD耦合分析方法。旋翼气动载荷通过求解三维Navier-Stokes方程求得,空间离散及通量计算采用Jameson中心格式,时间方向则选用五步Runge-Kutta迭代求解,湍流模型采用B-L模型;基于Hamilton原理建立了描述旋翼弹性运动的非线性微分方程,针对旋翼悬停状态的工作特点,采用Raphson迭代方法求解获得旋翼桨叶的弹性变形量。在CFD/CSD耦合计算中,旋翼桨叶交接面载荷及变形信息通过CFD与CSD模块进行传递,同时为提高桨叶弹性变形后贴体网格生成的效率和质量,采用基于网格点坐标转换的网格变形方法。在CFD和CSD程序分别验证基础上,采用建立的旋翼CFD/CSD耦合分析方法计算了先进的UH-60A直升机旋翼的表面压强及气动载荷。计算结果表明,与刚性旋翼CFD模拟结果比较,本文建立的CFD/CSD耦合分析模型可以更准确地预估旋翼气动载荷和性能。     

地面效应下的机翼失速数值模拟

采用计算流体力学的方法,对NACA0012二维翼型、三维机翼和带端板机翼在地面效应下的失速特性进行了数值模拟。通过对数值模拟结果的比较分析,揭示了地面效应下机翼的失速特性。研究结果为地效飞行器机翼的设计与优化提供了理论依据。     

遮挡板结构对二元喷管红外辐射特性的影响

在宽高比3.33的圆转矩形二元喷管的出口处安装下遮挡板,通过基于计算流体动力学/红外辐射(CFD/IR)数值模拟的方法来研究遮挡板结构对喷管壁面和尾焰的红外辐射强度的抑制效果.研究表明带狭缝气膜冷却的双层遮挡板对排气系统的红外辐射具有更为显著的抑制效果.     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

结冰对翼型流场影响的研究

本文以NACA23012翼型为研究对象,采用CFD软件模拟了结冰对翼型流场的影响。结果表明：结冰对翼型的气动性能有着明显的影响,尤其是在大攻角时,升力系数显著下降,而阻力系数则明显增加;对于光滑翼型和结冰翼型．计算值与试验值在‘定攻角范围内吻合良好：马赫数对两种翼型的升力系数和阻力系数随攻角的变化影响均比较小,雷诺数也一样。     

弹性飞行器大机动飞行气动特性分析

针对两种不同气动布局的前、后掠翼无人机,采用CFD/CSD松耦合方法,对弹性飞行器承受过载为10和20情况下的气动特性与刚性飞行器进行了对比分析。计算结果表明,与刚性飞行器假设情况相比,考虑过载作用下的弹性前掠翼飞机升力减小、阻力减小;而弹性后掠翼飞机的升力增大、阻力增大;两种气动布局的弹性飞行器与刚性情况相比,升阻比都略微增大,而且俯仰力矩系数明显增大。这对于飞行器纵向气动操纵特性会产生较大影响,同时也说明在大机动飞行过程中计入飞行器的弹性变形是非常必要的。