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

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

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

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

垂直/短矩起降飞机机翼内埋式风扇布局气动特性分析

以垂直/短距起降飞机过渡飞行状态为背景,针对机翼内埋式风扇布局的自由来流/风扇喷流混合型流动,基于结构/非结构混合网格使用CFD方法进行了非定常数值模拟和分析.首先使用滑移网格技术对NASA涵道螺旋桨进行算例验证,其时均计算结果与实验值的误差为5.3%,证明了计算方法的可靠性和准确性,然后数值模拟了机翼内埋式风扇布局在不同迎角下的气动性能.结果表明:风扇喷流在机翼上产生了特有的“抽吸”和“堵塞”效应,引起了机翼总升阻力的显著增加,升力最大增量达到干净机翼升力的2.6倍,阻力最大增量为干净机翼阻力的3.2倍,混合流场在机翼后缘引起了升力损失并卷起对涡.      

塞锥气膜冷却对二元塞式喷管红外特征的影响

为了研究二元塞式喷管塞锥壁面多斜孔气膜冷却对喷管红外特征的影响,设计了带冷却结构的喷管试验模型,测量了吹风比从0增加到1时的塞锥壁面温度、喷管出口截面处喷流温度以及喷管红外辐射强度.结果表明：随着吹风比的增加,气膜冷却效率逐渐提高,最大达到0.68;喷管出口截面在中心区域喷流的温度逐渐降低,最大降幅为14%;喷管红外辐射强度逐渐降低,在0°方位角上最大降低52.8%,在90°方位角上最大降低13%.     

针栓式喷注器锥形液膜破碎特性试验

采用高速摄影获得了针栓式喷注器在不同喷注压降和结构参数下的表面波破碎图像,测量了锥形液膜的破碎长度和破碎时间,研究了变工况时液膜破碎长度和破碎时间的变化规律.试验结果表明:在喷注压降不变的条件下,针栓式喷注器能够实现流量的线性调节.针栓式喷注器设计时,在合理的推进剂动量比范围内,狭缝宽度应尽量取小.液膜在低工况时破碎得更快.液膜破碎长度和破碎时间均随喷注压降的增加而减小.      

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

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

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

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

分开排气式喷管喷流噪声预测及试验研究

为揭示涡扇发动机分开排气式喷管喷流噪声的频谱分布和指向特性,对喷管缩比模型远声场喷流噪声进行了预测及试验研究,并开展了全尺寸喷管基准型及锯齿型喷流噪声的预测工作。结果表明:采用Tam & Auriault方法可以预测出喷流噪声的声压级值和空间指向分布,喷流噪声具有明显的指向性。与基准型喷管相比,锯齿型喷管具有良好的降噪效果;内外涵锯齿型喷管在中低频段的降噪效果优于内涵锯齿型喷管,在高频段的降噪效果差于内涵锯齿型喷管;在本文研究参数范围内,指向角小于88°时内外涵锯齿型喷管总声压级(OASPL)值高于内涵锯齿型喷管。     

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

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

带反馈控制的基本带翼外形的耦合气动/运动数值模拟研究

介绍了一种先进的、时间精度的多学科数值模拟平台的开发和应用,该平台可用于预测飞行器机动过程中的非定常气动特性和飞行力学特性.通过耦合求解非定常Navier-Stokes方程、刚体动力学方程及控制律,对基本带翼外形进行了虚拟飞行数值模拟研究.求解器采用基于结构网格的有限体积法,时间方法采用双时间步方法.使用反馈控制器获得舵偏控制指令,重叠网格方法用来模拟物体间的相对运动.模拟结果显示控制器的增益必须正确选择才能获得正确的响应.仿真结果表明本文发展的虚拟飞行数值模拟方法能够预测包含动边界的复杂非定常流动,具有较强的工程适用性.