高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

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

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

90°多弯管阻力系数计算方法

为准确获得90°组合多弯管阻力系数,采用数值模拟方法对S弯管管内流动进行模拟,并通过分析归纳S弯管局部阻力的相邻影响机制及规律,提出了90°组合弯管阻力系数计算方法.研究结果表明:弯管阻力系数受2弯连接段长度影响的规律与文献试验结果一致,数值计算方法准确有效;当连接段长度大于4倍管径时,S弯管阻力系数呈线性增长;提出的90°组合弯管阻力系数计算方法考虑了多弯管间相邻影响,且采用提出的方法计算了3弯管和多弯管的阻力系数,并将结果与数值模拟结果进行了对比,二者之间的差距均在3％之内,适合工程应用.     

侧载及加热方位对槽道内临界热流密度影响

为了研究侧载以及加热方位对矩形窄缝槽道内临界热流密度特性的影响,在旋转平台上进行了两相流的超重力实验.以蒸馏水为工质,通过改变质量流速、入口温度、侧载大小以及加热方位,获得了发生临界换热现象时的质量流速、实验段压降和壁温的变化趋势,考察了侧载和加热方位对临界热流密度的影响规律,并对侧载下两相流不稳定性进行了简析.实验表明:发生临界换热时,壁温迅速上升,有效加热热流开始减少,关闭加热电源后,实验段压降下降、质量流速回升较明显;侧载以及不同的加热方位对临界热流密度有明显影响;临界热流密度发生前后,蒸汽回流导致进口温度上升,增加了流动的不稳定性.      

涡流发生器对高负荷扩压叶栅性能影响的机理分析

为探明涡流发生器流动控制技术对高负荷扩压叶栅性能影响及作用机理,根据高负荷扩压叶栅的流动特点,提出了在叶栅入口端壁处加涡流发生器的流动控制方案,通过计算研究了采用涡流发生器前后叶栅气动性能、附面层及主要旋涡结构的变化。研究结果表明:采用涡流发生器后,叶栅正攻角下的气动性能显著提升,总压损失减小,静压升增大,稳定工作最大正攻角从3°增加至5°,其中在3°攻角下总压损失系数下降0.028,静压系数提升0.033;涡流发生器生成的尾涡阻挡端壁附面层由压力面向吸力面的横向迁移,使吸力面/端壁区域聚集的低能流体减少,改善了角区流动;采用涡流发生器后,通道涡、集中脱落涡和壁角涡减弱,角区分离得到抑制。     

启动过程对大膨胀比喷管外压失稳的影响研究

针对上面级发动机大膨胀比喷管高模试车时发生的外压失稳现象,进行了高模试车时的稳定工作和启动过程数值仿真和实验研究。结果表明:对于室压4.5MPa,环境压力87k Pa,喷管面积比70的上面级发动机,稳定工作时,喷管承受的是内压载荷,不会发生外压失稳;不预抽真空启动时,喷管内的流场建立过程所需的时间很短约0.4s,而试验舱的压力只能依靠发动机的引射而降低,从0.4s开始喷管承受较大的外压载荷,直至10s左右试验舱的压力低于喷管内壁压力,在较长时间内喷管一直承受外压载荷,很容易发生外压失稳;预抽真空启动时喷管在0.14s到0.42s时间段承受外压载荷,且只在0.14s到0.25s内承受的外压载荷较大,喷管承受外压载荷的时间很短,不会发生外压失稳;启动前抽真空是避免喷管发生外压失稳的一种有效手段。     

充气机翼的褶皱与失效行为研究

在经典工程梁理论的基础上,结合张力薄膜的应力状态分析,提出充气机翼褶皱失稳的判据。计入表面薄膜褶皱引起的刚度退化效应,将机翼等效处理为一个变截面刚度的梁,建立了充气悬臂机翼的等效梁模型,并采用微分求积法进行充气机翼弯曲变形分析。计算结果与充气机翼的静力弯曲试验结果相吻合,验证了充气机翼弯曲变形分析方法的有效性。应用片条理论引入气动力模型,并与所建立的等效梁模型相耦合,建立充气机翼的静气动弹性耦合模型,并用迭代算法进行求解。考虑起皱和失稳两种判据,并计算获取了试验机翼的起皱动压和皱褶失稳动压形式,计算结果与风洞试验结果一致。根据所建立的充气机翼静气动弹性分析方法,可以预测充气机翼表面褶皱区的扩展和弯曲变形,进而绘制充气机翼的静气弹许用包线,为充气机翼的设计提供必要的安全边界参考依据。     

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

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

基于IDDES方法的翼型结冰失速分离流动数值模拟

应用基于k-ωSST湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,就失速点附近翼型前缘典型双角状积冰导致的复杂分离流动进行了数值模拟研究.通过与风洞试验结果进行对比,表明对于此类分离流动问题,IDDES方法能够在壁面附近取得良好的速度预测结果,有效解析分离区域内的中小尺度湍流结构,较为准确地描述大尺度时均分离泡的再附位置和形态特征,适用于翼型结冰后复杂流动的精细分析.同时计算结果显示当此带冰翼型位于失速点附近时,角状冰后方脱落剪切层内部的旋涡不稳定析出和输运过程促进了外部流动与回流区域流动间的掺混,将导致流动发生非定常再附现象.     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。