Newton-Krylov全隐式时间推进方法及其在叶轮机械内流分析中的应用.pdfVIP

2025年10月推进技术Oct.2025

第46卷第10期JOURNALOFPROPULSIONTECHNOLOGYVol.46No.10

Newton-Krylov全隐式时间推进方法及其

在叶轮机械内流分析中的应用*

1111，231，2

李怡，张倩，赵家资，徐慎忍，乔磊，王丁喜

（1.西北工业大学动力与能源学院，陕西西安710129；

2.西北工业大学先进动力研究院，四川成都610299；

3.西北工业大学无人系统技术研究院，陕西西安710072）

摘要：为将Newton-Krylov全隐式时间推进方法应用于基于雷诺平均纳维斯托克斯（RANS）方程

的叶轮机械内流分析中，以此提高叶轮机械内流分析的收敛效率及鲁棒性，本文针对叶轮机械内流分析

相比飞机外流分析在方法上需要的补充，发展了考虑旋转周期边界、转静交界面及离心力源项的精确雅

可比计算方法。后将现有面向飞机外流CFD分析的Newton-Krylov全隐式时间推进方法拓展至多排叶轮机

械定常流场分析，并基于此方法开发了适用于飞机外流和叶轮机械内流分析的全隐式时间推进求解器。最

后，在一系列典型叶轮机械内流分析算例上验证了所开发的流场求解器的计算精度、收敛效率与鲁棒性。

关键词：Newton-Krylov方法；叶轮机械；动静交界面；鲁棒性；自适应库朗数

中图分类号：V231.1文献标识码：A文章编号：1001-4055（2025）10-202409013-18

DOI：10.3724/1001-4055.202409013

［11］

迭代进行三维黏性流场计算。

1引言

隐式迭代时需要在每个非线性迭代步计算残差

基于雷诺平均纳维斯托克斯（RANS）方程的计向量对流场变量的雅可比矩阵并求解相应的线性方

算流体力学（ComputationalFluidDynamics，简称程组，因而其计算成本和稳定性极度依赖此雅可比

CFD）技术广泛应用于飞机（翼型、机翼、整机等）和航矩阵的近似程度和计算方法、大型稀疏线性方程组

空发动机部件（压气机、涡轮、燃烧室等）的气动性能的求解方法以及非线性外迭代策略。对于这类非对

［1-6］

分析与外形设计。虽然非定常RANS（URANS）、称大型稀疏线性方程组，应用最广泛的是由不完全

大涡模拟（LES）、直接数值模拟（DNS）等非定常模拟LU分解（ILU）进行预处理的Krylov子空间方法，例如

方法目前进展迅速，但由于其较高的计算成本，目前［12］

GeneralizedMinimalResidual（GMRES）方法。为

主要用于压气机、涡轮等在非设计工况的流场分节约内存、简化计算，采用隐式时