突加不平衡激励碰摩转子反向涡动分析

建立了考虑叶片-机匣碰摩、挤压油膜阻尼的模型,推导了转子发生失稳的判别条件。从复非线性模态角度分析了突加不平衡激励下转子的动力特性,揭示转子反向涡动响应的形成过程及存在条件。通过参数分析获得转子发生反向涡动的敏感参数及其影响规律,并根据突加不平衡激励下转子反向涡动的响应特征分析某航空发动机叶片飞失故障。计算结果表明：转子能够发生反向涡动需要满足两个条件,其一,转子本身存在反进动模态失稳区;其二,冲击载荷使叶盘幅值达到反进动模态阻尼失稳点并进入反进动模态失稳区。实际航空发动机转子中具有因突加不平衡而发生反向涡动的风险,会造成支承结构破坏,严重威胁航空发动机的安全。增加转子阻尼、降低叶尖-机匣摩擦因数、降低静子叶片刚度、采用挤压油膜阻尼结构均有利于降低该风险。     

W型封严环表面粗糙度对间隙泄漏流动的影响

金属封严环密封结构的泄漏是影响发动机金属封严环间隙泄漏通道内气体流动的因素之一。本文采用有限元分析软件建立不同表面粗糙度下的泄漏通道模型;采用CFD 技术对不同结构的间隙内流体流动进行数值模拟,研究表面粗糙度对泄漏流动的影响;计算封严环的泄漏率,通过密封试验验证计算结果的正确性。结果表明：气体在泄漏通道内受到表面粗糙度的影响,靠近粗糙壁面处产生明显的回流现象,降低了通道内流体的泄漏流速;在接触应力不变的条件下,密封间隙的泄漏量随表面粗糙度的增大而增大,提高封严环的加工精度能够降低表面粗糙度从而增强其密封性能。     

基于多核并行计算的超燃冲压发动机一维模型实时性研究

在发动机控制系统设计中,为了缩短设计周期、降低研发成本,需要建立面向控制的、较为精确的、实时性高的超燃冲压发动机性能计算模型,以保证模型精度、提高计算速度为研究目标,基于多核高性能计算仿真平台,开展了面向控制的超燃冲压发动机一维模型实时性优化工作。运用简化计算流程、改进C语言程序、开拓缓存区等方法有效提高了一维模型计算速度。创新性地尝试了计算流体力学并行化方法,对隔离段和燃烧室一维模型进行结构分解。计算网格平衡分配至多个中央处理器,并借助核间数据通讯实现多核并行计算。与串行模型计算结果对比,七核并行计算模型性能参数偏差不超过0.1%,全工况仿真时间小于30ms,计算耗时较优化前缩短了75%以上。实时性优化后的多核并行模型计算精度高、速度快、收敛性好,可以作为超燃冲压发动机控制系统设计和半实物仿真验证平台。     

弱可压缩流体与可压缩流体模型的管路水击研究

为提升液体推进系统模型对可压缩流体的仿真能力,针对供应管路的瞬变流动,采用一维有限体积法,使用Modelica语言基于MWorks平台分别开发了适用于补燃循环发动机仿真模型库的弱可压缩流体与可压缩流体的管路模型,不同于控制方程中显含声速的弱可压缩流体模型,可压缩流体模型通过将守恒方程与真实流体的SRK型状态方程耦合求解来计算压力波的传播。在管路供应系统中,分别以液氢、液氧及液甲烷为工质进行关阀水击的仿真研究,结果表明:使用SRK状态方程避免了可压缩流体模型对内能拟合公式的依赖,减轻了建模难度;可压缩流体模型对液氧及液甲烷的仿真结果较为准确,最大误差不超过2%,对液氢的最大计算误差约9%;相对于可压缩流体模型,弱可压缩流体模型计算得到的水击压力在振幅及频率上均偏大;可压缩流体模型计算得到的流体密度波动不超过3%,因此工程实践中处理水击问题时将液体假设为弱可压具有一定的合理性。     

斜流工况下泵喷推进器推进性能数值分析

斜流入水的现象在推进器的运行过程中广泛存在,会对推进器的水动力性能造成不利影响。为了分析斜流工况对泵喷推进器推进效率及流场造成的影响,采用分块网格技术,对泵喷推进器内外流场进行了结构化网格划分。基于RANS 方程,SST k-ω 湍流模型,对不同转速（n=3000,3600,4200r/min）、不同斜流角（β=0°,5°,10°,15°,20°,25°,30°,40°） 的泵喷推进器进行了数值模拟。数值模拟方法通过Ka4-70/19A导管桨进行了验证,计算值与实验值吻合情况良好。研究结果表明,小角度下,斜流角的改变对泵喷推进器的效率影响很小;斜流角大于20°时,斜流角增大将使推进效率明显减少,这是由于泵喷推进器内流场分布不均,出现了水流聚集区和分散区,降低了泵喷推进器的平衡性以及对流体的加速效果;斜流角越大,叶尖间隙对主流影响越大,造成的能量损失越大。     

基于四面体的三维流线构造

流线是计算流体力学(CFD,Computational Fluid Dynamics)矢量可视化研究的重要技术.基于四面体单元,首先将点云数据场六面体单元分割成四面体单元,存储四面体之间的邻接关系;然后据此逐步将点定位到数据场中四面体单元内;再根据流线上速度矢量方向的变化量和四面体单元网格的内切球半径,自适应确定积分步长,以解决积分步长太小影响流线追踪速度或步长太大影响流线精度问题;最后利用龙格库塔数值积分法求得流线上点,并根据点定位和调整积分步长构造出三维流线.由于该算法避免了传统流线追踪算法中Jacobian矩阵逆变换,因而减小了矩阵转换过程带来的误差,提高了流线追踪的精度和效率.      

降落伞充满稳定阶段流场数值模拟

降落伞系统具有非常复杂的动力学特性,为了较精确的计算降落伞气动力,深入研究其流场特性,采用计算流体动力学方法,对某试验用降落伞充满稳定阶段流场进行数值模拟,计算结果与试验数据吻合较好,表明该数值模型可靠,具有较高精度。     

基于CFD/CSD耦合方法的新型桨尖旋翼气动弹性载荷计算

为提高旋翼非定常气动弹性载荷的分析精度,在刚性旋翼计算流体力学(CFD)方法中引入计算结构动力学(CSD)方法,建立了一套适合于新型桨尖旋翼气动弹性载荷分析的CFD/CSD耦合方法。旋翼流场分析采用NavierStokes/Euler方程作为控制方程,围绕旋翼生成运动嵌套网格。在流场求解中,采用双时间法推进,通量计算采用Jameson中心格式,并采用B-L(Baldwin-Lomax)湍流模型。基于Hamilton变分原理和中等变形梁理论开展桨叶弹性运动变形分析,并发展了一套具有任意转角梁单元的新方法以提高新型桨尖旋翼的动力学分析精度。采用基于代数变换方法的网格变形策略,建立了一套CFD/CSD松耦合方法,桨叶运动变形和旋翼气动力信息通过流固交接面传递。首先分别对CSD和CFD模块进行了验证,然后计算了UH-60A旋翼在高速前飞状态下的气动弹性载荷,并与试验值进行了对比,最后重点对旋翼桨尖形状进行了参数分析。计算结果表明,相比于升力线理论和刚性旋翼CFD方法,CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,并能更准确地反映新型桨尖旋翼的气动弹性耦合效应;同时采用后掠桨尖在桨叶前行侧30°~90°方位角范围可以显著降低激波强度,有利于改善旋翼的气动特性。     

基于高阶谐波平衡的跨声速颤振高效预测方法

跨声速流场激波及其诱导的附面层分离等非线性因素导致跨声速颤振边界很难被准确预测,尤其是目前工程常用的偶极子格网法,在跨声速时该方法的预测精度大幅下降。在雷诺平均Navier-Stokes方程流场求解器的框架内,利用结构模态建立广义结构运动方程,利用径向基函数建立模态振型的插值方法,结合径向基函数和无限插值两种网格变形方法的优点实现高效高鲁棒性网格变形方法,从而实现颤振时间推进分析流程,利用国际颤振标模AGARD445.6机翼验证程序在跨声速颤振边界预测中的可靠性。然而,时域方法在气动/结构反复迭代,需要耗费大量的计算资源和时间。为了提高颤振预测效率,基于高阶谐波平衡（HOHB）方法快速获得广义力影响系数矩阵,利用该矩阵建立频域模态位移和气动力之间的关系,实现高效颤振频域分析方法。通过二维翼型和三维机翼算例进行验证,结果表明：在不对计算精度产生明显影响的前提下,HOHB方法能够提高颤振预测效率约6倍。     

旋翼非定常动态失速的三维效应分析

采用建立的高精度计算流体动力学(CFD)方法,针对旋翼非定常动态失速的三维（3D）效应特性进行研究。以Helishape 7AD旋翼为基准,开展三维效应对旋翼非定常气动特性的影响分析。研究了来流马赫数对旋翼翼型动态失速特性影响。在此基础上,针对三维情形旋翼动态失速非定常涡流动特性及诱导分离特征开展了数值分析,通过与二维情形对比表明:受旋翼旋转、轴向诱导速度等三维效应影响,旋翼桨叶剖面动态失速涡的产生、对流和脱落明显滞后于二维翼型情形,并且涡强度也更弱。越靠近桨叶内段,桨叶剖面非定常动态失速特性与二维旋翼翼型情形的差距越为明显。