低动能来流下背负式进气道非定常流动特性分析

为了改善飞翼布局背负式S弯进气道低动能来流状态下的流动性能,采用改进的延迟分离涡模拟（IDDES）方法对原型及改进型背负式进气道流场进行了数值模拟研究,对比分析了进气道流量特性及内部脉动压力特性。结果表明:低动能来流时背负式进气道上部唇口附近存在很大的气流转折角,导致唇口产生分离涡;原型进气道唇口分离涡强度高,高能量分离涡在进气道顶部破裂产生了大范围旋涡结构,进一步加剧了流动分离,从而引发进气道内产生强烈的压力脉动,声压级最大幅值高达145 dB;改进型进气道唇口分离涡得到了有效控制,强度大幅下降,进气道内部压力脉动幅值也显著降低,声压级降幅达8 dB;改进型进气道分离的抑制使进气道有效流通截面积增大,质量流量增加。同时,流场出口品质提升,进气道出口综合畸变指数降低了9.5%。      

DES与DDES在湍流分离中的原理与性能研究

随着工程上流动结构的日益复杂,兼具雷诺平均Navier-Stokes(RANS)方法高效率与大涡模拟（LES）高精度的分离涡模拟（DES）类混合方法成为现阶段工程中最有效的湍流模拟方法之一。围绕DES类混合方法中的DES与延迟分离涡模拟（DDES）方法开展工作,分析二者开关函数构造上的不同,研究延迟因子作用机理,并考察DES与DDES方法的求解能力。研究表明：DES与DDES方法在模拟表现上存在一定差异,DDES方法通过引入延迟因子,保护RANS求解区域,改善模化应力不足,降低了DDES方法对交界面系数CDES敏感程度;DDES方法在计算过程中容易出现过度保护,导致求解瞬时涡结构能力不如DES方法,分析与延迟因子引入比重及开关函数构造形式有关。     

类X-51A飞行器纵向机动数值虚拟飞行仿真

吸气式高超声速飞行器在机动过程中,由于构型复杂,气动特性呈现强烈的非定常特性。传统基于数据库或气动力模型的飞行仿真不能准确描述机动过程中复杂的气动特性和运动规律。针对这一问题,基于现代软件分布式、模块化的发展趋势,建立了一个高效的数值虚拟飞行仿真平台。利用该平台,对一种类X-51A外形的吸气式高超声速飞行器开展了纵向机动闭环数值仿真,并与工程方法的结果进行了对比。研究发现:对于类X-51A外形的吸气式高超声速飞行器,在纵向拉起时,工程方法给出的结果可能不能完全反映非定常效应的影响。此时,应该采用更为精确的虚拟飞行方法来研究飞行器的闭环响应特性。此外,借助该仿真平台还研究了舵回路时间常数对控制系统的影响,为控制律设计提供了一定的参考。      

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

对旋轴流风机非定常流场数值模拟

利用全隐格式的非稳态压力修正求解方法,分析和研究了对旋叶轮级间流场非定常流动情况,上游尾迹与下游叶片位势作用对流场的影响,并与动态实验结果进行了有益的对比.结果表明:非定常流场的计算更能反映叶轮机实际运行性能,如果能很好地组织对旋叶轮的非定常流动,还可以进一步提高其压比和效率.     

浮空气囊定常及非定常气动特性计算

对于0~20 km范围内不同外形的浮空气囊,采用虚拟压缩方法求解不可压缩非定常N-S(Navier-Stokes)方程,数值模拟了5个高度下气囊的定常和非定常流场.定常状态时,随高度及气囊外形的变化,模拟了气囊气动力及背风区复杂流动分离.结果表明,当气囊厚度率逐渐增大时,其背风区存在复杂的空间流动分离.运用双时间步方法及动网格技术,成功实现了气囊平移或俯仰运动状态下非定常流场的数值模拟.考察了不同的平移方向及速度对气囊非定常气动特性的影响,得到了气囊平移过程中气动力及流场的变化情况.对于气囊的俯仰运动,获得了不同俯仰速度及不同最大俯仰角时气囊的气动力及流场.计算结果表明气囊的非定常运动对其流场和气动力影响显著.     

狭长型封闭舱室内非定常流动模拟的湍流模型比较

针对狭长型封闭舱室内非定常流动模拟的湍流模型选择问题,以飞机舱室为典型环境,使用相似准则为依据的热缩比法搭建了实验平台。将实验结果与RNG k-ε、DES、LES三种湍流模型的数值模拟结果进行对比和分析,评估狭长型封闭舱室内非定常流动特征研究中合适的湍流模型。结果显示,RNG k-ε和DES模型可以定性描述流动变化趋势,但是LES模型在流场非定常性和不稳定性捕捉更为准确,其流场结构更接近实验结果。模拟结果与实验结果对比显示,LES模型能更加真实地反映狭长型封闭舱室非定常流动的情况。     

低雷诺数下翼型前缘流动分离机制的研究

采用高精度有限差分格式,对来流雷诺数为1.0×104,攻角为3°的二维翼型流动进行了直接数值模拟,研究了低雷诺数下翼型前缘流动的分离机制,描述了分离涡系的相互作用规律.计算结果表明:前缘椭圆弧靠近叶身位置存在吸力峰,流动在吸力峰内强逆压梯度的作用下发生分离;翼型上表面形成了包含驻留涡、脱落涡和二次涡的涡系结构,其尺度随时间不断变化,具有强烈的非定常性;表面压力分布曲线可以较好的描述翼型边界层流动.      

流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

汽轮机末级三维非定常流动数值模拟

空冷汽轮机组采用空气作为冷却介质,是一种典型的变工况运行机组.深入了解透平叶片在设计状态和高背压条件下的非定常流动机理,能更好设计空冷汽轮机以及提高叶轮机械的性能和稳定性.主要利用三维粘性非定常数值模拟的方法对设计状态和高背压条件下透平叶片内部的流动进行模拟,并分析了设计状态和高背压条件下非定常流动机理.结果表明,在设计状态动静干涉作用是导致非定常现象产生的主要原因;在高背压条件下动静干涉作用很弱,导致非定常现象不明显.     

喷管分离流动及其侧向载荷

利用商业软件CFX对某液体火箭大面积比喷管地面条件下的分离流动进行了三维数值模拟.计算获得了喷管入口总压从8MPa减小到1MPa时的流场参数分布和侧向力载荷情况.结果表明,随着入口总压的降低,喷管内流场会依次经历自由激波分离和受限激波分离两种分离激波模态.受限激波分离模态下喷管壁面压强具有较大波动,再附着点压强甚至高于环境压强.流动分离情况下,喷管将受到一定侧向载荷作用,载荷方向随机分布.入口总压为4MPa时计算得到的侧向载荷最大,实际侧向载荷峰值可能出现在自由激波分离与受限激波分离转换瞬间.     

前飞旋翼三维湍流场的数值模拟

采用计算流体力学CFD(Computational Fluid Dynamics)方法,分析正常前飞和近地前飞时直升机的流场与性能.在直角坐标系中,采用有限体积法和压力耦合方程的半隐式法SIMPLER(Semi-Implicit Method for Pressure-Linked Equations Revised)算法求解三维定常不可压的湍流N-S(Navier-Stokes)方程.在分析中,旋转的螺桨被描绘成沿螺桨桨叶展向分布、与本地流动参数相关以及时间平均的动量源项.计算方法还包括K-ε湍流模型和壁面函数法等措施.流场分析和性能预测同实验测量数据的良好的一致性表明,这种CFD方法可以有效的分析直升机的具体设计问题.      

使用GAO-YONG湍流方程组 对顶盖驱动方腔流的计算

运用SIMPLER方法、QUICK格式及交错网格技术求解GAO-YONG湍流方程组,对二维顶盖驱动方腔流进行数值模拟,得到方腔内的涡量分布和方腔中心水平剖线上的速度变化曲线,并与DNS(Directly Numerical Simulation)结果进行定性对比,有较好的一致性.结果表明,GAO-YONG湍流方程组中与湍流多尺度现象相对应的机械能方程对此湍流方程组能够合理地模拟出顶盖驱动方腔流中复杂的湍流粘性分布以及湍流能量逆转现象起到了关键作用.研究进一步证明了GAO-YONG湍流方程组能够得到复杂流动的真实粘性场.     

电弧推力器约束通道内流动特性数值模拟

约束通道对电弧推力器的性能有着重要的影响,文章采用基于局域热力学模型（LTE）的数值模拟方法对中等功率电弧推力器内等离子体流动进行了数值模拟,考察了电流、入口压力、约束通道尺寸及不同推进剂对约束通道内等离子体流动的影响,分析了约束通道内非均匀流动现象,最后对推力器的性能、效率等进行了讨论。计算结果表明,随着电流的增加电弧高温区变粗变长,随着入口压强的增加电弧高温区半径减小而长度增加,随着约束通道半径的减小电弧高温区变得细长,随着约束通道长度的增加高温区的长度增长而半径无明显变化,氢气的高温区明显小于氮气和氩气;约束通道内只有小部分气体通过高温区被电离,大部分气体沿着壁面附近的低温区流动;约束通道内焦耳热约占总焦耳热的60%~80%,主要受约束通道长度影响。