超声速平板边界层旁路转捩直接数值模拟

采用具有七阶精度的有限差分格式和基于Euler方程组的特征边界条件,通过直接数值模拟研究了由固定相位周期性吹吸扰动诱导的、Ma∞=2.25、Re∞=2.5×107/m(635 000/inch)空间发展平板边界层的旁路转捩过程,并建立了相应的直接数值模拟数据库。数值结果表明:扰动相位的随机性并不是转捩发生的必要条件,固定相位周期性吹吸扰动即可诱发边界层的转捩;当采用固定相位时,随着吹吸频率的增加,转捩位置向下游移动;当扰动频率足够大时,周期性吹吸无法诱导出边界层的转捩。随后,在所得数据库中提取了速度及热力学变量的平均值,与已有的数据库和理论分析吻合良好;速度脉动均方根与公认的不可压缩平板边界层实验基本吻合;热力学变量脉动均方根和雷诺应力与已有数据库的统计结果及实验数据具有相同的分布特征。对比发现,不同转捩的诱导因素对平均转捩位置和热力学变量均方根、雷诺应力等高阶统计量的影响较大。     

热声载荷下薄壁开孔结构振动响应与寿命预估

金属薄壁开孔结构在强热声载荷下,表现出大挠度强非线性响应特性,疲劳寿命缩短。基于时域的基础上,利用有限元法(FEM)结合降阶模态法(ROM)获取4边固支开孔薄板在不同热、声载荷组合下的位移和应力的动态响应,并进行响应时间历程和功率谱(PSD)的统计分析。采用雨流计数法和Morrow平均应力模型,结合Miner线性累积损伤理论对结构进行疲劳寿命的预估和分析。结果表明:屈曲后的位移响应由热声载荷的相对强弱决定。此外,屈曲前结构随热、声载荷的增加,寿命缩短。屈曲后,持续跳变使得结构的寿命缩至最短;进入间歇跳变区域,间歇时间短的跳变要比间歇时间长的跳变造成的结构损伤大,即快频跳变引起的损伤更大。跳变结束后,随着温度的升高,寿命先延长后缩短。     

一种超声速状态下机载推进系统模型的建模方法

为了满足机载航空推进系统模型对精度、实时性及数据存储量的苛刻要求,提出了一种机载模型建模方法.对传统的机载稳态模型建立过程中使用的相似准则进行了充分讨论,得出了一种具备更高建模精度的相似准则.在新的相似准则基础上,建立了包含泰勒展开式中非线性2阶项的高精度稳态变量模型,有效解决了分段线性稳态变量模型在大包线、变状态条件下精度难以保证的问题;为了解决超声速状态下安装推力与发动机净推力差别较大的问题,建立了考虑外流特性的简化进气道模型,给出了一种计算安装推力的方法,即由简化发动机模型计算安装推力中符合相似准则的部分,由进气道模型计算溢流阻力、放气阻力及发动机推力中不符合相似准则的部分.仿真结果表明:基于新的相似准则所建立的机载模型具有大包线、变状态适应性,输出参数误差均在0.5%以内.      

叶尖泄漏流和引射效应对跨声速串列转子失速机制的影响

为探究跨声速串列转子的失速机制,利用数值模拟的手段对不同间隙情况的某跨声速串列转子的叶尖流场进行分析.研究发现:前排二次泄漏区会对失速产生重要影响;叶尖泄漏流和引射效应是前排二次泄漏区流场结构的决定因素,而间隙大小直接影响着两者的强弱.当前排大间隙时,叶尖泄漏流强,是前排二次泄漏区流场的主导因素,这时失速首先发生在前排;当前排间隙减小时,引射效应增强,并逐渐成为主导因素,这时前排二次泄漏区的堵塞情况得到改善,失速区域转移到后排.改变前排间隙对该跨声速转子的性能有更显著的影响,当前排为最先失速排时,改变后排间隙对串列转子性能没有显著影响.      

基于多岛遗传-模拟退火算法的气动外形优化

针对超声速翼型滑翔机的气动外形优化问题,提出了一种多岛遗传算法与模拟退火算法相结合的混合优化算法。首先,通过多岛遗传算法产生一个随机的初始种群,找到全局最优点附近区域的一个次优解;然后,将此次优解作为模拟退火算法的初值启动退火进程,缩小设计空间的范围,找到全局最优解。优化结果表明,所提混合算法可以有效地解决飞行器外形优化问题;在数量巨大的设计空间中,滑翔机的最优气动外形能够以较低的计算资源代价快速得到;双弧形翼型的气动特性较六边形翼型更有优势,叉形尾翼的静稳定性和控制效率高于十字形;优化后的外形极大地增加了滑翔机的滑翔距离。     

各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

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

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

基于拟动力学的航空发动机主轴球轴承热弹流润滑分析

以航空发动机主轴球轴承高速高温重载典型工况为算例,基于轴承拟动力学分析,得到滚动体与套圈之间的微接触区运动和受力状态,分别用Hamrock-Dowson(H-D)拟合公式和翟文杰热修正公式计算了最小膜厚(MFT),同时结合热弹流润滑分析得到最小膜厚分布,对比了试验测试、H-D拟合、翟文杰热修正和热弹流润滑分析4种算法的最小膜厚,并根据膜厚比判定了轴承的弹流润滑状态.结果表明:H-D拟合和翟文杰热修正最小膜厚与试验测试结果相差较远,热弹流润滑分析结果与试验数据吻合性较好,误差10%以内;探讨了不同工况下的轴承弹流润滑性能.结果表明:存在径向载荷时,润滑膜整体压力增大,膜厚减小;随着转速增大,压力减小,膜厚增大,但是增大的趋势减缓.      

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

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

周向槽轴向位置影响机匣处理扩稳能力的机理

为了揭示周向槽机匣处理轴向位置影响压气机稳定性的流动机理,采用试验与数值方法研究了三种不同轴向位置对机匣处理扩稳能力的影响。试验与数值结果均表明机匣处理轴向位置不变时扩稳能力最强,轴向位置前移次之,轴向位置后移最弱,对应的试验综合裕度改进量分别为10.3%,7.82%及5.65%。通过详细分析压气机叶顶流场表明,周向槽机匣处理轴向位置后移使叶顶前缘区域没有受到到周向槽的作用,叶尖部分间隙泄漏流在叶顶前缘形成低速带并造成进口堵塞。轴向位置前移使叶顶后部区域没有受到周向槽的作用,在叶顶通道出口部分叶高范围内形成低速堵塞区,流通能力弱。轴向位置不变的机匣处理对应的叶顶通道流通能力最强,因此扩稳效果最好。