地外天体起飞过程真空羽流导流力热效应研究

为研究地外天体起飞真空羽流对探测器分离产生的力热扰动,使用计算流体力学-直接模拟蒙特卡洛(CFD-DSMC)耦合计算模型对锥面导流的真空羽流场进行了计算。采用组分输运模型计算三维连续流场,并获取当地的克努森数作为判断连续流和离散流的依据。使用基于分子动力学的直接模拟蒙特卡洛方法(DSMC)计算离散流场,采用可变软球(VSS)碰撞模型和Larsen and Borgnakke传能模型计算离散流分子间的能量传递,将计算结果与试验进行了对比,验证了计算方法的可靠性。研究结果表明,A器受到的侧向干扰力矩为62N·m,底部受到的最大压力为100Pa,最大热流密度为100k W/m2;B器受到的侧向干扰力矩为558 N·m,表面最大压力为8k Pa,最大热流密度为600k W/m2,喷口与导流装置顶面距离为400mm时,激波已移出喷管内部。     

应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

飞机机头附近液态水含量分布规律的数值研究

在飞行器飞行过程中,安装在机头上的各类传感器结冰会对飞机的飞行安全产生不利影响。而液态水含量是影响飞机结冰的重要参数,为了更好地进行传感器结构及布局设计,对机头附近液态水含量分布规律进行数值研究。采用计算流体力学方法对飞机飞行时的空气-过冷水滴的稀疏两相流场进行数值模拟,通过求解N-S方程获得空气流场,采用欧拉法求解水滴运动场得到水滴运动轨迹。分析在不同来流马赫数、来流迎角及液滴直径等条件下,机头关键截面上广义水滴遮蔽高度的变化及机头关键位置上液态水含量的变化规律,研究不同状态参数对液态水含量分布的影响。结果表明:水滴遮蔽高度及浓度增加区范围随马赫数增大而增加,迎角的变化对机头上下表面液态水含量分布产生相反的影响。     

以柴油为燃料的燃烧室燃烧效率实验及分析

地面动力常常需要改进利用航空燃气轮机技术,为了研究燃气轮机燃烧室在低工况下使用柴油燃料的燃烧效率问题,采用燃气分析试验方法,对燃气轮机燃烧室在不同进口总温、空气流量和油气比等工况下的燃烧效率进行了研究。结果表明:油气比较低时,燃烧区温度低,燃烧不完全导致燃烧效率急剧下降,随着燃烧室油气比的增加,燃烧效率逐渐接近100%;进口空气流量以及进口总温的增加会提高空气雾化喷嘴的雾化能力以及燃烧温度,燃烧更加充分完全,燃烧效率由96%左右提高至99%以上;总结归纳得到了适用于柴油燃料的燃烧效率预估经验关系式。     

固体火箭发动机流体喉部推力矢量特性

针对采用水作为二次流工质的流体喉部进行了冷流实验及数值模拟研究.研究了该种固体火箭发动机流体喉部的一般规律,包括不同二次流射流方式,不同二次流流量下流体喉部的扼流性能,推力偏角及推力效率,数值模拟及实验结果吻合较好.结果表明:扼流性能与二次流的注射位置、注射角度及流量比有关,且随二次流/主流流量比的增大而增大.喉部二次流喷射能有效的调节有效喉部面积进而调节推力大小,当流量比为0.4时,最大有效喉部面积比为0.8;扩张段二次流喷射能有效调节推力方向,当流量比为0.4时,最大推力偏角为20°;喉部二次流与扩张段二次流入射位置存在相位差可有效降低喉部与扩张段二次流干扰.      

北京首都机场畅通滑行时间研究

随着机场规模不断扩大和流量持续增长,准确把握拥挤态势成为运行管理中的难题。美国联邦航空局(FAA)和欧洲空管局(EUROCONTROL)采用畅通滑行时间(unimpeded taxi time)作为量化分析机场运行效能的标杆。国内笼统地规定了少数机场的平均滑行时间作为机场运行指标,但无法反映不同机场、不同停机位、不同起飞跑道造成的运行差异。因此可借鉴机场畅通滑行时间这个指标来考核机场效能,但是欧美机场的该指标数值并不完全适合国内机场,因此在欧美计算方法的基础上提出了改进的畅通滑行时间计算方法,并构建多元线性回归模型,运用北京首都机场实际运行数据进行检验,证明了该方法的可行性。     

冲击距与气膜孔方位角对旋流气膜冷却性能影响

为了获得冲击距Hi与气膜孔方位角α对旋流气膜冷却性能的影响规律,以六边形供气腔圆形气膜孔平板气膜冷却结构为研究对象,对五种冲击距参数(0.74D,1.14D,1.54D,1.94D,2.34D)(D为气膜孔直径)和五种气膜孔方位角参数(0°,10°,15°,20°,25°)进行了三维数值计算研究,得到了绝热壁面气膜冷却效率、展向平均气膜冷却效率、流场空间无量纲浓度分布等随冲击距与气膜孔方位角的变化规律,分析了肾形涡对旋流气膜冷却性能的影响机理。结果表明:冲击距对绝热壁面气膜冷却效率展向分布规律影响不大,而方位角增加能够明显提高绝热壁面气膜冷却效率及展向气膜覆盖面积,方位角0°模型展向气膜冷却效率最大值为0.42,方位角25°模型的最大值为0.48,相比前者增加14.3%;绝热壁面同一流向位置的展向平均气膜冷却效率随冲击距的增加而增大,随方位角的增加而增大,Hi=2.34D时的展向平均气膜冷却效率最佳,α=20°时的展向平均气膜冷却效率曲线最佳。方位角的增加能够明显破坏流场中存在的肾形涡结构。     

软式空中加油受油机头波数值仿真分析

针对软式空中加油对接过程中受油机头波加剧软管锥套飘摆、降低对接成功率的问题,建立了受油机三维模型,利用Delaunay方法自下向上生成非结构网格.通过计算ONERA M6三维翼型表面压力分布和压力系数,验证了三维Navier-Stokes方程和SST k-ω湍流模型的有效性.在此基础上对不同飞行条件下受油机头波进行静态数值计算,分析了头波随马赫数和迎角的变化规律,并通过仿真解释了头波影响下软管锥套的飘摆原因,预测了对接时锥套的运动趋势.     

气膜孔局部堵塞对叶片压力面冲击-扰流柱-气膜结构综合冷却效率的影响

运用数值模拟方法研究了气膜孔局部堵塞对叶片压力面上射流冲击-扰流柱-气膜结构综合冷却效率的影响,重点分析了堵塞位置和堵塞比的影响。研究结果表明：无论是气膜孔内无堵塞还是存在局部堵塞情形,随着吹风比增大,综合冷却效率均呈现逐渐增大的趋势;在低的吹风比下,气膜孔出口-尾缘局部堵塞的综合冷却效率略低于无堵塞气膜孔,而在气膜孔进口-前缘和气膜孔出口-前缘的局部堵塞则导致综合冷却效率有较为轻微的上升;在高吹风比下,位于气膜孔出口-前缘和气膜孔进口-前缘的局部堵塞能够抑制气膜射流与主流相互作用所形成的卵形涡,从而冷却效率下降较少,而气膜孔出口-尾缘的局部堵塞则导致综合冷却效率降低较多;堵塞比对壁面沿程综合冷却效率的影响呈现非单调的变化趋势,这是由于冲击-扰流柱-气膜整体式冷却结构的冷却效果取决于内部强化传热和外部气膜防护的多重作用机制。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。