(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111386641.1 (22)申请日 2021.11.22 (71)申请人 北京航空航天大 学 地址 102206 北京市海淀区学院路37号 (72)发明人 潘天宇　李秋实　武文倩　郑孟宗　 (74)专利代理 机构 北京鼎承知识产权代理有限 公司 11551 代理人 顾可嘉　夏华栋 (51)Int.Cl. G06F 30/28(2020.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (54)发明名称 针对进发匹配的压气机稳定边界预测方法、 设备及介质 (57)摘要 本发明提供了一种针对进发匹配的压气机 稳定边界预测方法、 设备及介质， 包括叶片域采 用所述分布式力源体积力模型计算， 非叶片域采 用雷诺平均方法进行计算， 从而获得压气机负荷 的空间分布。 本发明基于对压气机负荷空间分布 及临界负荷的认识对流动稳定性进行分析和预 测， 认为当压气机存在区域其负荷达到当地临界 负荷时会触发失稳， 具有较高的预测精度。 此外 由于采用基于分布式力源体积力模型的计算方 法， 能够实现进发匹配情况下压气机负荷空间分 布的快速获取， 因此能够大幅提高稳定边界预测 工作的效率。 权利要求书2页 说明书7页 附图4页 CN 114117959 A 2022.03.01 CN 114117959 A 1.一种针对进发匹配的压气机稳定边界预测方法， 其特 征在于包括以下步骤： S1:在压气机均匀进气情况 下， 采用定常雷诺平均方法计算获得压气机流场参数； S2:根据所述压气机流场参数及叶片参数， 构建分布式力源体积力模型， 输出分布式力 源体积力； S3:叶片域采用所述分布式力源体积力模型， 非叶片域采用雷诺平均方法进行非定常 计算， 从而获得压气机负荷的空间分布； S4:根据所述压气机负荷的空间分布实施判断， 若压气机整体或任意局部区域的负荷 均未达到当地临界负荷， 则 返回S3且提高出口平均静压再次计算； 若压气机整体或存在局 部区域的负荷达 到当地临界负荷， 则压气机失稳， 输出最后一个稳定点。 2.根据权利要求1所述的一种针对进发匹配的压气机稳定边界预测方法， 其特征在于： 所述压气 机流场参数及叶片参数包括叶片的当地半径、 叶片周向角度坐标、 叶片表面压力、 气流当地速度沿子午面的投影速度、 叶片的当地切线速度、 气流的当地切线速度、 当地静温 度、 当地气流密度、 当地相对坐标系下 的气流速度、 沿流向的熵梯度、 当地流向速度对应的 马赫数。 3.根据权利要求1所述的一种针对进发匹配的压气机稳定边界预测方法， 其特征在于： 所述非叶片域包括进气道、 进气道出口至第一排叶片域前缘、 相 邻叶片域之 间、 最后排叶片 域尾缘至 压气机出口 的全部区域。 4.根据权利要求1所述的一种针对进发匹配的压气机稳定边界预测方法， 其特征在于： 所述当地临界负荷为压气机局部发生失稳时的当地负荷。 5.根据权利要求1所述的一种针对进发匹配的压气机稳定边界预测方法， 其特征在于： 还包括， S5:获得最后一个稳定点和第一个失稳点， 当这两点间的节流步长足够小时， 以最后一 个稳定点作为该进气道情况 下的压气机预测稳定边界； 否则返回S3 。 6.根据权利要求1所述的一种针对进发匹配的压气机稳定边界预测方法， 其特征在于： 所述分布式力源体积力包括叶片对流体的切向体积力， 沿等流向坐标面体积力和沿等叶高 面体积力。 7.根据权利要求1所述的一种针对进发匹配的压气机稳定边界预测方法， 其特征在于： 所述叶片域与非叶片域的交界面为叶片前尾缘在子午面的投影沿周向旋转一周后形成的 曲面。 8.根据权利要求1所述的一种针对进发匹配的压气机稳定边界预测方法， 其特征在于： 进气道进口边界条件给定总压、 总温以及速度方向， 压气机出口边界条件给定平均静压条 件， 并通过减小出 口静压实现对压气机节流过程的模拟， 此外所有实壁边界均为无滑移边 界条件。 9.一种可读存储介质， 其特征在于， 其上具有可执行指令， 当可执行指令被执行时， 使 得计算机执 行如权利要求1 ‑8的任一项所包括的操作。 10.一种计算设备， 其特 征在于， 包括： 一个或多个处 理器； 存储器； 以及 一个或多个程序， 其中所述一个或多个程序存储在所述存储器中并被配置为由所述一权　利　要　求　书 1/2 页 2 CN 114117959 A 2个或多个处 理器执行如权利要求1 ‑8中的任一项所包括的操作。权　利　要　求　书 2/2 页 3 CN 114117959 A 3