-
随着以美国太空探索技术公司(SpaceX)、维珍银河公司(Virgin Galaxy)和蓝色起源公司(Blue Origin)为代表的太空公司的兴起,太空旅游、商业发射等产业蓬勃发展,世界对低成本、可重复使用运载器的发射需求进一步增长[1]。
“猎鹰9号”(Falcon-9)火箭是SpaceX公司研制的可重复使用运载火箭,采用栅格舵作为一子级垂直返回过程中的稳定控制面。在一子级回收过程中,栅格舵打开,且进行摆动,用于一子级回收过程中调整姿态,实现落区控制[2-6]。蓝色起源公司的“新格伦”(New Glenn)重复使用火箭以第一位完成绕地球轨道飞行的美国宇航员约翰•格伦(John Glenn)的名字命名,其运载能力将超过SpaceX的“重型猎鹰”(Falcon Heavy)和联合发射联盟(United Launch Alliance,ULA)的“德尔塔Ⅳ 型”重型火箭(Delta Ⅳ Heavy)”[7-10]。不同于SpaceX采用栅格舵进行回收控制的外形,“新格伦”火箭第一级采用两片滑翔翼和4片空气动力侧翼,以满足再入返回阶段气动减速和控制的需求[11]。
本文主要对基于栅格舵和滑翔翼的两种垂直起降火箭一子级再入返回段的气动特性进行计算,分析其不同攻角再入的气动特性规律,研究其性能的优劣和任务适用模式。
-
“猎鹰9号”火箭是SpaceX公司研制的中型、低成本、两级构型的全液体运载火箭。自2010年“猎鹰9号”火箭v1.0型首飞成功开始,该系列火箭经历大量技术上的升级换代,期间推出“猎鹰9 v1.1型”“猎鹰9 v1.2-FT型”(包含Block 1-4亚型),最终定型在“猎鹰9 v1.2-Block 5”。栅格舵是一种由许多栅格组成的具有蜂窝式结构的新型弹翼,用作火箭垂直返回的控制面,如图1所示。
“新格伦”火箭第一级具有空气动力滑翔翼,如图2所示。当一级分离后调转方向并再次点火,空气动力滑翔翼可减缓下降速度,并在着陆前展开6个着陆支架,最终将垂直降落在海上的一艘无人平台上。“新格伦”火箭尾部的两片滑翔翼和级间段的4片空气动力侧翼,可控制箭体飞向着陆点并实施着陆,同时取消了再入减速点火,节省了箭体返回所需的大量推进剂。
采用有限体积(Finite Volume)方法,通过数值求解三维坐标系下的Navier-Stokes方程[12]进行数值模拟。采用时间相关法求解流场,通过有限体积法离散方程,将流场划分成一系列有限大小的体积单元,每一个体积单元上都有一个节点作为代表,将各方程进行分离求解。离散方程的统一形式可表示为
$$ \frac{\rm d}{{{\rm d}t}}{(\rho \phi V)_{{P_0}}} + {a_0}{\phi _{{P_0}}} = \sum\limits_{j = 1}^N {{a_j}} {\phi _{{P_j}}} + {b_0} $$ (1) 对于稳态过程,
$\dfrac{\rm d}{{\rm d} t}(\rho \phi V)_{P_0}=0$ 。式(1)中:ρ为大气密度;ϕ为广义变量,代表流场内的待求参数(如速度、温度等);$ P_{0} $ 为控制容积中心节点;$ P_{j} $ 为其周围节点;$ V $ 为控制容积体积;$j$ 为控制容积中界面的角标;$ N $ 为控制容积的界面数;$ a_{0} $ 和$ a_{j} $ 分别为离散方程中对应于$ P_{0} $ 和$ P_{j} $ 的算子系数;$ b_{0} $ 为离散方程源项。采用Spalart-Allmaras(SA)一方程湍流模型[13]来计算湍流效应,该模型适用于无限制流和附着流动,目前已经广泛应用于工程实际中,尤其是在火箭、飞机机翼、翼身组合体等绕流的数值模拟中[14]。SA模型的控制方程如下
$$ \frac{\rm d}{{ {\rm d}t}}(\rho \tilde v) = {G_v} + {D_v} - {Y_v} $$ (2) 其中:
$ \tilde v $ 为模型因变量;$ {G_v} $ 、$ {D_v} $ 和$ {Y_v} $ 分别为湍流黏性生成项、黏性扩散项和黏性耗散项,具体如下$$ {G_v} = {C_{b1}}\rho \tilde S\tilde v $$ (3) $$ {D_v} = \frac{1}{{{\sigma _{\tilde v}}}}\Bigg[\frac{\partial }{{\partial x{}_j}}\Bigg((\mu + \rho \tilde v)\frac{{\partial \tilde v}}{{\partial {x_j}}}\Bigg) + {C_{b2}}{\Bigg(\frac{{\partial \tilde v}}{{\partial {x_j}}}\Bigg)^2}\Bigg] $$ (4) $$ {Y_v} = {C_{w1}}\rho {f_w}{\Bigg(\frac{{\tilde v}}{d}\Bigg)^2} $$ (5) 湍流黏性系数可以由下式得出
$$ {v_t} = \tilde v{f_{v1}} $$ (6) 式(3)~(6)中,各系数如下
$$ {f_{v1}} = \frac{{{\chi ^3}}}{{{\chi ^3} + C_{v1}^3}} \chi = \frac{{\tilde v}}{v} $$ $$ \tilde S = S + \frac{{\tilde v}}{{{\kappa ^2}{d^2}}}{f_{v2}} S = \sqrt {2{\varOmega _{ij}}{\varOmega _{ij}}} $$ $$ {f_{v2}} = 1 - \frac{\chi }{{1 + \chi {f_{v1}}}} {\varOmega _{ij}} = \frac{1}{2}\Bigg(\frac{{\partial {u_i}}}{{\partial {x_j}}} - \frac{{\partial {u_j}}}{{\partial {x_i}}}\Bigg) $$ $$ {f_w} = g{\left[\frac{{1 + C_{w3}^6}}{{{g^6} + C_{w3}^6}}\right]^{\frac{1}{6}}} g = r + {C_{w2}}({r^6} - r) $$ 其中:
$ v $ 、$ \mu $ 分别为分子运动黏性系数和分析动力黏性系数;$ {f_{v1}} $ 为阻尼函数;$ d $ 为该控制单元的物面距;$\varOmega$ 为涡量的模;$ u $ 为速度。出现的常数如下:${C}_{b1}=0.133\;5$ ,$ {C_{b2}} = 0.622 $ ,${\sigma _{\tilde v}} = {2/3}$ ,$ {C_{v1}} = 7.1 $ ,$ {C_{w1}} = \dfrac{{{C_{b1}}}}{{{\kappa ^2}}} + \dfrac{{1 + {C_{b2}}}}{{{\sigma _{\tilde v}}}} $ ,$ {C_{w2}} = 0.3 $ ,$ {C_{w3}} = 2.0 $ ,$\kappa =0.418\;7$ 。计算时采用的空间格式为高分辨率的Roe通量差分裂(Flux Difference Splitting,FDS)格式,时间格式为高可靠性、高效率的LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。
-
栅格舵构型一子级外形示意如图3所示。一子级全长17 m,直径3 m,质心距火箭原点的相对距离为0.6,力矩参考点为(0,0,0)。建模时对发动机喷管里的复杂结构进行简化处理,栅格舵与一子级连接部位的结构也做了简化。对栅格舵构型一子级周围流体区域进行网格划分,栅格翼网格采用分区对接网格技术,以O型网格为主。区与区之间满足光滑过渡,在边界层和栅格舵前缘部分适当加密。其中亚音速流场域取为15倍的参考长度,超音速流场域选用3倍的参考长度,计算网格量1 000 万。箭体远场网格示意如图4所示。
滑翔翼构型一子级外形示意如图5所示,一子级全长20 m,直径3.5 m,质心距火箭原点的相对距离为0.6,力矩参考点为(0,0,0)。建模时对发动机喷管里的复杂结构进行简化处理,对滑翔翼构型一子级周围流体区域进行网格划分,采用结构化全模网格。在靠近箭体区域进行了必要的加密,边界网格y+<100。由于主要计算状态为超音速,计算区域可适当缩小,流场域选用3倍的参考长度,全箭采用网格数约为273万,箭体远场网格示意如图6所示。
为了验证网格无关性,以带栅格舵一子级为例,对一子级生成粗、中、密3套网格,网格量及第一层网格高度见表1。以马赫数(Ma)4、攻角30°,侧滑角0°的计算条件对3套网格进行了数值模拟,获得一子级轴向力系数CA并进行了误差分析。从表1中统计的相对误差可以看出,3套网格采用的网格量及第一层网格高度满足气动力计算的网格无关性要求。选取第一层网格高度为0.01 mm用于气动力的计算,此时网格量适中,可兼顾计算精度和效率需求。
表 13种计算网格情况
Table 1.Three different kinds of mesh
参数 粗 中 密 总网格量/万 675 1 000 1 500 第一层网格高度/mm 0.015 0.01 0.007 轴向力系数CA 1.862 9 1.860 1 1.879 8 CA相对误差/% 0.15 — 1.06 -
计算区域外边界采用远场边界条件,箭体壁面边界采用无滑移、绝热壁面边界条件。远场条件根据需要计算的工况详细给定来流压力、马赫数和方向、温度以及空气的物性参数(包括:比热、层流黏性系数、导热系数等);采用理想气体模型,不同飞行状态对应不同的高度和大气密度。
层流黏性系数和导热系数取为定值;边界条件设定如图7所示;理想气体的比热容可以表示成温度的函数
$$ {c_{\rm p}} = {a_0} + {a_1}T + {a_2}{T^2} + {a_3}{T^3} $$ (7) 式(7)中系数由实验测定,对于理想气体,通常取
${a}_{0}=1\;161.482$ ,${a}_{1}=-2.368\;8$ ,${a}_{2}=0.014\;9$ ,${a}_{3}=-5.034\;909 \times {10^{-5}}$ 。 -
对栅格舵构型一子级和滑翔翼构型一子级再入返回段10°和30°攻角状态下的气动特性进行数值仿真计算,得到不同马赫数下的流场结构,如图8~10所示。
图 8攻角10°、Ma= 0.8工况栅格舵构型一子级流场分布图
Figure 8.The velocity flow distribution of the first sub-stage with gridrudder forα= 10°,Ma= 0.8
图 9攻角10°、Ma= 2.0工况栅格舵构型一子级流场分布图
Figure 9.The velocity flow distribution of the first sub-stage with gridrudder forα= 10°,Ma= 2.0
对栅格舵构型一子级,亚音速状态下,来流在一子级头部开始减速,压力升高。经过喷管后,靠近壁面处,因壁面的黏性影响和攻角的影响,上壁面流动因绕过一子级而形成了较大范围的分离区,下壁面低速区范围较小;远离壁面处气流绕过一子级肩部加速,然后向下游流经栅格舵,3号栅格舵处于靠近壁面的低速流场区域,1号栅格舵则处于较为接近自由来流的流场区域中。
超音速状态下,来流在一子级头部形成激波,气流减速增压,由于攻角的影响,上激波远离壁面,下激波靠近壁面。经过一子级肩部后膨胀加速,向下游流经栅格舵,3号栅格舵处于靠近壁面的低速流场区域,1号栅格舵则处于较为接近自由来流的流场区域中。
亚声速时随着来流马赫数增加,一子级肩部的分离区先增大然后减小,其形态呈左右非对称分布,至超音速状态,分离区逐渐减小,其形态逐渐向对称面靠近形成对称的形态。对比马赫云图可以发现,随着来流马赫数的增大,一子级头部逐渐形成激波,马赫数越大,激波越贴体。
栅格舵构型一子级的关键控制部件是栅格舵,栅格舵流场较为复杂,图11为不同马赫数下单独栅格舵流场图。Ma= 0.8时,栅格舵在亚音速条件下,来流过栅格舵前缘绕流膨胀再压缩,在壁面的边界层影响下,栅格内里的有效横截面面积减少,栅格内里形成超声速流动区域,并在尾部区域存在结尾激波,流动重新减速为亚声速流动。
图 10攻角10°、Ma= 7.0工况栅格舵构型一子级流场分布图
Figure 10.The velocity flow distribution of the first sub-stage with gridrudder forα= 10°,Ma= 7
图 11不同马赫数下栅格舵流场图
Figure 11.The velocity flow distribution of the gridrudder for different Mach number
Ma= 1.0与Ma= 0.8流动类似,随马赫数的增加,栅格舵前缘亚声速来流减速,通道内气流出现壅塞,部分流量从栅格外侧溢流。栅格后部的超声速流动区域增大,轴向力系数明显升高,在Ma= 1左右达到峰值,如图11(b)所示。
Ma= 1. 5时,来流马赫数大于1,需通过激波改变入口流速和压力来调节栅格内流量,从而在栅格入口前方形成一道很强的弓形激波,由于波后仍为亚音速流动,栅格内里流动与Ma= 1.0时流动近似。
Ma= 3.0时,随着马赫数提高,弓形脱体激波逐渐向栅格舵靠近,当大于第二临界马赫数时,超声速气流能通过收缩通道,此时弓形激波变成斜激波且靠近舵前缘,相邻栅格产生的两道斜激波大约在上下两片隔片中间相交,相交区域后方形成一个高压区。
Ma= 5.0时,激波相交的区域继续向栅格出口方向移动,波后高压区对隔片下表面压力影响减小,导致栅格法向力降低,如图11(e)所示。
对滑翔翼构型一子级再入返回段Ma= 2.0 ~ 9.0,10°攻角和30°攻角状态下的气动特性进行数值仿真计算,得到了各工况的气动计算结果。图12展示了攻角10°、Ma= 2.0工况“新格伦”火箭一子级流场分布图,在Ma= 2.0的超音速状态下,来流在一子级尾部的喷管处形成激波,气流减速增压,经过一子级肩部后膨胀加速,向下游流经滑翔翼和侧翼时又产生激波,在一子级级间段处形成了一个分离区。图13展示了攻角10°、Ma= 9.0工况“新格伦”火箭一子级流场分布图,在Ma= 9.0的超音速状态下,激波角大幅减小,在一子级级间段处也形成了一个分离区。图14展示了攻角30°、Ma= 9.0工况“新格伦”火箭一子级流场分布图,在30°大攻角状态下,来流经过一子级尾部的喷管后,在背风面产生大的分离区,在迎风面产生贴体激波,且激波发展至侧翼后,产生相互作用。
图 12攻角10°、Ma= 2.0工况滑翔翼构型一子级流场分布图
Figure 12.The velocity flow distribution of the first sub-stage with glider forα= 10°,Ma= 2.0
-
对于火箭来说,压心和质心之间的相对位置对飞行中的稳定性起决定性作用,如果压心位于质心之后,则火箭具有静稳定性,反之具有静不稳定性。在设计中,一般引入静稳定裕度(
$ \eta $ )的概念,如式(8)所示$$ \eta = ({x_{\rm cp}} - {x_{\rm cg}}) \times 100 \text% $$ (8) 其中:
${x_{\rm cp}}$ 为压心距火箭原点的相对距离(压心系数);${x_{\rm cg}}$ 为质心距火箭原点的相对距离(质心系数)。如果压心和质心的相对位置在飞行中变化很大,将对控制系统设计带来困难,因压心和质心的相对位置控制在一定的范围内。用$ \Delta \eta $ 表征火箭飞行期间质心的移动程度,一般$ \Delta \eta $ 不超过8%,当质心在压心之后,且质心与压心相对位置变化较小时,则为静稳定性较优。图15展示了栅格舵构型一子级10°和30°攻角下气动力系数随马赫数的变化规律。轴向力系数CA随马赫数增大先增大后减小,在Ma=2左右达到峰值。法向力系数CN的变化规律与轴向力系数CA相似,也是随马赫数增大呈现先增大后减小的规律。10°攻角时,俯仰力Cmz大体上是随着马赫数基本不变。
图 15栅格舵构型一子级气动力系数变化曲线
Figure 15.Aerodynamic force and moment coefficient at different Machs of the first sub-stage with gridrudder
压心系数xcp维持在距离原点0.6 ~ 0.75的位置,在质心之后,且随着马赫数变化较小,维持在质心位置附近,因此静稳定性较好。
图16给出了滑翔翼构型一子级10°和30°攻角时的轴向力系数、法向力系数、俯仰力矩系数和压心系数随马赫数的变化规律,其中俯仰力矩系数和压心系数相对于坐标原点。由图可知,在
马赫数大于2之后,轴向力系数CA随 马赫数没有明显变化,10°攻角的轴向力系数高于30°攻角。法向力系数CN随 马赫数增大而逐渐减小,在马赫数大于5后呈线性减小趋势,30°攻角的法向力系数要远高于10°攻角。10°攻角时俯仰力矩系数Cmz随马赫数变化不大,压心系数xcp维持在距离原点0.35 ~ 0.45的位置,在质心位置(xcg=0.6)之前,是静不稳定的;30°攻角时压心靠后,维持在距离原点0.80左右的位置,在质心之后,是静稳定的。 图 16滑翔翼构型一子级气动力系数变化曲线
Figure 16.Aerodynamic force and moment coefficient at different Machs of the first sub-stage with glider
通过对比基于栅格舵构型和滑翔翼构型一子级的压心系数范围可知,在10°和30°攻角下,基于栅格舵构型的一子级可以维持较大范围的静稳定状态,且在马赫数和攻角变化时,压心在一定范围内变化,降低了对舵面控制的要求;而基于滑翔翼构型的一子级在10°攻角下是静不稳定的状态,对舵面控制效率提出较高的要求。
图17给出了栅格舵构型一子级和滑翔翼构型一子级升阻比随马赫数的变化规律。栅格舵构型一子级的10°攻角的升阻比在0.05~0.3之间,随马赫数的增大先增大再减小;30°攻角的升阻比在0.5~0.8之间,随马赫数的增大先增大再减小。滑翔翼构型一子级10°攻角时升阻比在Ma= 2时为0.75,之后随马赫数的增大而逐渐减小,在马赫数大于5后呈线性减小趋势。30°攻角时升阻比维持在1.0左右。在滑翔阶段,一子级采用10°以内的小攻角飞行,而10°攻角条件下,滑翔翼构型一子级的升阻比是栅格舵构型的2~5倍,升力特性更优越,有利于滑翔段的增升减速。
通过两种构型一子级的气动特性对比可知,栅格舵构型一子级的静稳定性优于滑翔翼构型一子级,适合于需要精确控制落点的任务,可以在指定位置定点着陆。滑翔翼构型一子级的升阻比较大,更适合需要远距离滑翔的任务,机动性能更强,可以在更大范围内选择着陆点。
-
本文对基于栅格舵和滑翔翼的两种构型火箭一子级再入返回段的气动特性进行了计算和分析,从气动特性规律、静稳定性和升阻特性3个方面进行了对比,结论如下。
1)在气动特性分布规律方面,基于栅格舵构型和滑翔翼构型的火箭一子级,在小攻角下的轴向力系数CA、法向力系数CN、俯仰力矩系数Cmz随Ma变化规律相似;在大攻角下俯仰力矩系数Cmz规律不同。
2)在静稳定性方面,小攻角下,基于栅格舵构型一子级的压心系数位于质心之后,而滑翔翼位于质心之前,基于栅格舵构型一子级的静稳定性优于滑翔翼构型一子级,更适合于需要精确控制落点的任务。
3)在升阻特性方面,基于滑翔翼构型一子级的升阻比高于栅格舵构型,更适合需要远距离滑翔的任务。
An Analysis of Aerodynamic Characteristics of Reusable Rocket’s First Sub-Stage with Grid Rudder and Glider
-
摘要:针对基于栅格舵和滑翔翼的两种垂直起降火箭一子级再入返回段的气动特性进行了仿真分析,得到了两种构型一子级的气动特性,分析了其流动结构特征,并研究了气动特性随马赫数( Ma)和攻角的变化规律。结果显示:基于栅格舵构型的一子级轴向力系数和法向力系数随马赫数先增大后减小;基于滑翔翼构型的一子级在小攻角下轴向力系数和法向力系数在马赫数大于2后随马赫数增大而逐渐减小。通过对比可知,基于栅格舵构型一子级的静稳定性优于滑翔翼构型一子级,更适合于需要精确控制落点的重复使用任务;基于滑翔翼构型一子级的升阻比高于栅格舵构型,更适合于需要远距离滑翔的重复使用任务。Abstract:The aerodynamic characteristics of the reentry stage of the first sub-stage of two vertical take-off and landing rockets based on grid rudder and glider are simulated and analyzed. The aerodynamic characteristics of the first sub-stage of two configurations are obtained. The flow structure characteristics of the sub-stage are analyzed, and the variation of aerodynamic characteristics with Ma number and angle of attack is studied. The results show that the axial and the normal force coefficient of the sub-stage with grid rudder first increase and then decrease with the Ma number. The axial and normal force coefficient of the sub-stage based on glider configuration decrease gradually with the increase of Ma number when Ma number is above 2 at small angle of attack. By comparison, the static stability of the sub-stage with grid rudder configuration is better than the sub-stage with glide configuration, ,which means that the sub-stage with grid rudder was more suitable for a mission requiring accurate control of the landing point. Meanwhile, the lift-drag ratio of the sub-stage with glider configuration is higher than that of the grid rudder configuration, indicating that the sub-stage with glider was more applicable for long distance gliding.
-
Key words:
- grid rudder/
- glider/
- sub-stage/
- reusable rocket
Highlights● The aerodynamic characteristics of grid rudder and glider for reusable rocket first sub-stage are compared. ● The sub-stage with grid rudder is more suitable for a mission requiring accurate control of the landing point. ● The sub-stage with glider is more applicable for long distance gliding. -
表 13种计算网格情况
Table 1Three different kinds of mesh
参数 粗 中 密 总网格量/万 675 1 000 1 500 第一层网格高度/mm 0.015 0.01 0.007 轴向力系数CA 1.862 9 1.860 1 1.879 8 CA相对误差/% 0.15 — 1.06 -
[1] MARSHALL T. Anover view of the NASA aeronautics test program strategic plan[C]//27th AIAA Aerodynamic Measurement Technology andGround Testing Conference. Reston: AIAA, 2010. [2] 郑雄,杨勇,姚世东,等. 法尔肯9 可重复使用火箭发展综述[J]. 导弹与航天运载技术,2016(2):30-46.ZHENG X,YANG Y,YAO S D,et al. Survey and review on development of Falcon 9 reusable rocket[J]. Missiles and Space Vehicles,2016(2):30-46. [3] 牟宇,孙冀伟,秦旭东. 猎鹰9火箭Block5构型首次飞行任务解析[J]. 宇航总体技术,2018,2(5):1-7.MOU Y,SUN J W,QIN X D. The analysis of Falcon 9 Block5’s maiden flight[J]. Astronautical Systems Engineering Technology,2018,2(5):1-7. [4] CLARK S. Falcon 9 launch timeline with Bangabandhu1[EB/OL]. (2018-05-11)[2022-07-06]. http: //spaceflightnow.com/2018/05/11/falcon-9-launch-timeline-with-Bangabandhu1. [5] 胡炜. SpaceX, 一时传奇or一世传奇?[M]. 北京: 中国宇航出版社, 2015. [6] 龙雪丹,魏威. SpaceX商业载人运输系统首次无人试飞解析[J]. 飞航导弹,2019(3):5-11,44.doi:10.16338/j.issn.1009-1319.20190817 [7] 廖小刚,王岩松. 2018年国外载人航天发展综述[J]. 载人航天,2019,25(1):122-127.doi:10.3969/j.issn.1674-5825.2019.01.018LIAO X G,WANG Y S. Review of human spaceflight development abroad in 2018[J]. Manned Spaceflight,2019,25(1):122-127.doi:10.3969/j.issn.1674-5825.2019.01.018 [8] 何慧东,肖武平,苑艺,等. 2020年全球载人航天发展综述[J]. 国际太空,2021(2):24-37. [9] 廖小刚,王岩松. 2020年国外载人航天发展综述[J]. 载人航天,2021,27(1):127-134.doi:10.16329/j.cnki.zrht.2021.01.020LIAO X G,WANG Y S. Review of human spaceflight development abroad in 2020[J]. Manned Spaceflight,2021,27(1):127-134.doi:10.16329/j.cnki.zrht.2021.01.020 [10] 王辰等. 可重复使用运载火箭发展研究[J]. 飞航导弹,2018(9):18-26.doi:10.16338/j.issn.1009-1319.20180017 [11] SPISZT S, TAYLOR J C, GIBSON D, ETAL. Processing in frared imagery of the SpaceX Falcon first stage reentry during CRS-4 mission[C]//AIAA SPACE and Astronautics Forum and Exposition. Reston: AIAA, 2017. [12] BLAZEK J. Computational fluid dynamics: principles and applications[M]. Oxford: Elsevier Science Ltd, 2001. [13] SPALART P R, ALLMARAS S R. A one-equation turbulence model for aerodynamic flows: AIAA Paper 92-0439[R]. [S. l. ]: AIAA, 1992. [14] 杨凌,李大春,钟兢军. 不同湍流模型对平面叶栅气动性能预测的研究[J]. 中国科技论文,2017,12(23):2690-2694.doi:10.3969/j.issn.2095-2783.2017.23.009LI L,LI D C,ZHONG J J,et al. Research on predicting accuracy of turbulence models for transonic plane cascade aerodynamic[J]. China Sciencepaper,2017,12(23):2690-2694.doi:10.3969/j.issn.2095-2783.2017.23.009 -