Design and actualization of the variable density plane cascade wind tunnel
-
摘要: 为满足先进涡扇发动机对变雷诺数平面叶栅试验的需求，设计了亚/跨/超声速来流高效变换、雷诺数和马赫数独立调节、压气机和涡轮平面叶栅试验为一体、换热与冷却试验能力兼具的变密度平面叶栅风洞，提出了风洞的总体设计方案。文章详细介绍了风洞引射器、半柔壁喷管及试验舱等部件设计问题，分析了流场调试及典型叶栅试验结果。调试结果表明：采用的部件设计技术实现了变密度平面叶栅风洞的主要功能，试验雷诺数可低至3.1×105 m–1，具备开展低雷诺数平面叶栅试验的能力。风洞流场调试结果满足《低速风洞和高速风洞流场品质要求》（GJB 1179A—2012），为研究亚/跨/超声速压气机和涡轮叶栅低雷诺数流动问题提供了重要试验平台。Abstract: In order to meet the need of the variable Reynolds number plane cascade tests for the advanced turbofan engine, the variable density plane cascade wind tunnel was designed to change the subsonic, transonic and supersonic flow efficiently, adjust the Mach number （Ma） and Reynolds number （Re） independently, combine the compressor and turbine cascade tests and possess the ability of heat transfer or cooling experiments. The general design project of the wind tunnel was put forward, and the ejector, the flexible nozzle as well as the test chamber design problems were introduced in details. The results of the flow field debugging and the typical cascade tests were analysed. The research reveals that the design technologies of the components satisfy the main functions of the variable density plane cascade wind tunnel. The Reynolds number can be as low as 3.1×105 m–1 so that the low Reynolds number experiments can be made easyly with the facility. The flow field debugging results satisfy the National military standard GJB 1179A—2012 of the specification for flow quality of high and low speed wind tunnels. It provides a key test platform for the study of the transonic, supersonic Mach number and low Reynolds number flow problems of the compressor and turbine cascades.
-
Keywords:
- plane cascade /
- wind tunnel /
- turbomachinery /
- Reynolds number /
- Mach number
-
0 引言
目前，要求地面试验设备完全模拟高超声速飞行环境是非常困难的，在一个试验设备上进行所有环境的模拟试验更是不可能的。因此，需要研制多种试验设备，以满足高超声速试验需求[1-5]。
常规高超声速风洞主要模拟参数是马赫数和雷诺数[6]，其他一些参数则主要通过高焓试验设备来模拟[7-8]。
在高超声速风洞运行中，气源不仅要维持风洞所必须达到的压力比，还需要满足雷诺数模拟的要求。一般来说，风洞的总压要足够高。如果风洞连续工作，高的压力比和高总压将使风洞消耗很大的动力，因而高超声速风洞多为暂冲式。
暂冲式常规高超声速风洞的运行方式是：气罐中的压缩空气经过加热系统达到所需温度，然后通过型面喷管，在试验段形成所需的高超声速流场，最后经由超声速扩散段升压后进入引射排气系统，排入大气，或经过冷却器进入真空系统，由真空泵系统排入大气。常规高超声速风洞一般加热温度在800K以下，主要是为了防止气流冷凝，不能真实模拟实际总温(马赫数6时达到1800K)，因此，开展发动机及飞行器一体化带动力试验还需要满足总温要求的高焓设备。
燃烧加热风洞是目前高超声速飞行器地面模拟试验尤其是带动力试验的主要设备。在过去的十几年间，中国空气动力研究与发展中心研制了不同尺度的脉冲燃烧风洞，探索了一体化飞行器设计、计算与试验技术[7-9]。在此基础上，发展了基于脉冲燃烧风洞的大尺度飞行器带动力一体化试验技术[10-11]，提出了一种基于一体化试验直接测力结果的飞行器和发动机性能评估方法。
由于采用燃烧加热方式，无论是采用氢燃料还是碳氢燃料，都会产生水蒸气(H2O)、二氧化碳(CO2)等燃烧产物，即形成所谓“污染”。污染组分将造成风洞试验气体物理化学属性与真实空气存在一定差异，导致地面试验难以完全模拟真实飞行状态下的所有来流参数，且给地面试验结果向真实飞行状态的外推带来不确定性。为了尽可能降低这种不确定性，地面试验一般需要慎重考虑试验来流与模拟对象环境之间的参数匹配问题：即通过有选择地调整污染来流的某些状态参数，使之逼近对应的真实飞行环境，同时放弃一些无法兼顾的非关键参数，以尽可能达到飞行器气动与推进性能的可靠模拟。
目前主要有两种模拟方式：对于气动力试验，一般采用静温、静压、马赫数模拟；对于发动机性能试验，也可以采用总焓、动压、马赫数模拟。
中国科技大学杨基明、罗喜胜等[12-13]通过典型升力体飞行器的试验和计算研究，证明污染组分使得斜激波波后参数相对纯净空气有一定变化，从而造成模型表面压力及模型的气动力有一定变化，但变化量较小。相对于纯净空气来流，污染空气来流时发动机推力、单位推力、比冲均有所下降。脉冲燃烧风洞采用氢氧燃烧加热，污染空气来流时推力性能下降相对较小，与纯净空气来流时较为接近。
西北工业大学宋文艳等[14]研究了H2O/CO2污染煤油燃料对超声速燃烧室的影响，认为污染对点火和超声速燃烧具有一定的抑制作用，污染组分的存在会导致燃烧室模态转换点发生变化。
谭宇等[15]在酒精燃烧和氢气燃烧两种加热方式的风洞设备上开展了匹配方案对超燃冲压发动机性能影响的试验研究，比较了两种目前较常用的气流参数匹配方案，结果表明：对于采用氢气燃烧加热方式的风洞设备，总焓动压马赫数(h0QM)匹配比静温静压马赫数(TPM)匹配能够获得更高的壁面静压和推力收益；对于采用酒精燃烧加热方式的风洞设备，两种匹配方案表现相当。
壁温比通常定义为模型壁面温度与恢复温度的比值。常规高超声速风洞来流温度较低，气流壁面恢复温度与壁面温度相差不大，接近于绝热状态；而燃烧加热脉冲风洞来流温度较高，风洞试验时间很短(300~600ms)，壁面温度基本保持常温，接近于等温状态，气流壁面恢复温度与试验模型表面温差较大，气流对壁面的加热效应明显。壁面温度条件会影响到高超声速边界层内的流动参数，进而影响到流场的波系结构和气动性能。壁温降低引起摩擦系数的增加和边界层厚度减小，使得激波边界层作用区域变小；同时，冷壁使得边界层内亚声速部分声速较低，马赫数更高。
为了评估真实飞行条件下的飞行器气动性能，需重点考虑壁温对燃烧加热脉冲风洞试验结果的影响。
本文采用不通气标模，在常规高超声速风洞以及两个不同尺度(Φ2.4m和Φ600mm)的脉冲燃烧风洞中开展对比测力试验，结合数值计算，研究脉冲燃烧风洞水凝结、雷诺数及壁温比对模型气动性能的影响规律。
1 试验模型
试验模型采用不通气标模。对应Φ2.4m脉冲燃烧风洞，试验模型为大不通气标模，采用背部支撑，如图 1所示；对应常规高超声速风洞和Φ600mm脉冲燃烧风洞，试验模型为1/5缩比的小不通气标模，采用尾部支撑，如图 2所示。
2 脉冲燃烧风洞水凝结影响研究
为了获得水凝结对脉冲燃烧风洞试验数据的影响，采用小不通气标模，在Φ600mm脉冲燃烧风洞开展不同总温条件下的对比试验研究，分析测力结果。
2.1 对比试验参数
对Φ600mm脉冲燃烧风洞喷管进行了试验配套改造，在来流总温、组分不同的条件下，保持风洞出口马赫数一致。以现有的Ma6喷管为基础(总温1500K)，设计制造总温为1200和1800K的喉道段，喷管出口马赫数与现有喷管一致。表 1为不通气高超标模试验参数(p0、T0分别为总压、总温，p、T分别为静压、静温)。
表 1 不同总温高超标模试验参数Table 1 Flow parameters of different T0 testsMa p0/MPa T0/K p/Pa T/K Re 6.0 4.5 1200 2368 162 7.08×106 6.0 5.1 1500 2398 210 4.98×106 6.0 5.8 1800 2374 263 3.61×106 2.2 试验结果
图 3、4给出了不通气高超标模的轴向力系数CA、法向力系数CN随迎角α变化的试验结果曲线。
![]()
图 3 高超标模不同总温试验CA~α曲线Fig. 3 CA~α graph of different T0 tests for the typical hypersonic model![]()
图 4 高超标模不同总温试验CN~α曲线Fig. 4 CN~α graph of different T0 tests for the typical hypersonic model试验表明：在总温1200~1800K试验状态下，最大轴向力系数差别为5.4%，说明水凝结对不通气高超标模的轴向力影响很小，约为5%左右。在-2°迎角时，总温1200K(有水凝结)试验状态的法向力系数较总温1500K减小了约40%；在0°~6°迎角下，法向力系数减小12%~5%，说明随着迎角的增加，模型下表面的凝结水蒸发，导致水凝结对法向力的影响减小。在-2°迎角时，总温1800K试验状态的法向力系数较总温1500K增加了约30%，说明总温1500K、-2°迎角状态在模型的下表面也发生了水凝结，导致模型下表面流场压力升高，法向力系数减小；而在其他迎角状态，总温1500K时并没有发生水凝结现象，因此法向力系数差别较小。
3 脉冲燃烧风洞雷诺数影响研究
3.1 对比试验参数
为了获得雷诺数对试验结果的影响，用不同尺度的不通气标模在Φ2.4m脉冲燃烧风洞和Φ600mm脉冲燃烧风洞进行对比试验，试验参数见表 2(CHIF：脉冲燃烧风洞)。
表 2 试验参数Table 2 Test parametersFacility Ma T0/K p0/MPa Φ2.4m CHIF 6 1575 4.76 Φ600mm CHIF 6 1500 4.61 3.2 试验和计算结果
图 5和6给出了大、小不通气标模在两个脉冲燃烧风洞的轴向力系数CA、法向力系数CN试验结果，同时也给出了两个模型按照对应风洞来流条件计算获得的气动力系数结果。
![]()
图 5 不通气标模轴向力系数计算和试验结果Fig. 5 Numerical and experimental axis force coefficients of the test model![]()
图 6 不通气标模法向力系数计算和试验结果Fig. 6 Numerical and experimental normal force coefficients of the test model计算结果表明：模型尺度对不通气标模的轴向力有一定影响，小不通气标模的轴向力略大于大不通气标模，相差约2%~5%。由于两个脉冲燃烧风洞模拟条件基本相同，试验时间也接近，因此这个差异可以解释为模型尺度(雷诺数)造成的：小不通气标模的雷诺数小于大不通气标模，因此轴向力略大。试验并未能够反映出雷诺数的影响，这可能是由于两个模型支撑方式不同以及试验测量误差掩盖了差异。
两个模型的法向力结果不一致，主要是由于支撑方式不同造成的。前期研究表明，腹部支撑对轴向力影响不大，但对法向力影响显著。
计算与试验结果符合较好，使计算和试验得到了相互验证。
4 脉冲燃烧风洞壁温比影响研究
采用小不通气标模，通过Φ600mm脉冲燃烧风洞及常规高超声速风洞的试验对比，结合数值计算，研究壁温比对脉冲燃烧风洞试验结果的影响。
4.1 对比试验参数
Φ600mm脉冲燃烧风洞与Φ1m常规高超声速风洞的试验流场存在总温、雷诺数、壁温比等方面的差异。假设壁温300K，风洞流场参数对比如表 3所示(HWT：常规高超声速风洞；Tratio为壁温比)。
表 3 Φ1m常规高超声速风洞与Φ600mm脉冲燃烧风洞流场参数Table 3 Test parameters of Φ1m hypersonic wind tunnel and Φ600mm combustion heated impulse facilityFacility Ma p0/MPa T0/K γ Re Tratio Φ1m HWT 6.0 3.00 500 1.40 2.62×107 5.17 Φ600mm CHIF 6.0 4.65 1560 1.36 4.10×106 1.49 4.2 试验和计算结果
脉冲燃烧风洞试验与计算结果表明：采用等温壁300K条件计算获得的阻力系数与试验数据比较接近，误差在6%以内；绝热壁计算阻力系数较等温壁300K计算小5%(迎角6°)~15%(迎角-2°)；不同壁温条件对升力系数的计算结果影响不大，与试验结果的误差分别在6%和8%以内；采用等温壁300K条件计算获得的结果与试验基本一致。图 7和8给出了脉冲燃烧风洞试验和计算的阻力系数CD及升力系数CL的对比曲线。
![]()
图 7 CD~α试验与计算结果对比(脉冲风洞，总温1500K)Fig. 7 CD~α comparison of numerical and experimental results from Φ600mm combustion heated impulse facility![]()
图 8 CL~α试验与计算结果对比(脉冲风洞，总温1500K)Fig. 8 CL~α comparison of numerical and experimental results from Φ600mm combustion heated impulse facility常规高超声速风洞试验和计算结果表明：不同壁温条件对计算结果影响不大，计算结果与试验符合较好，阻力误差在4%以内，升力误差在3%以内。图 9和10给出了相应的阻力系数和升力系数比较曲线。
![]()
图 9 CD~α试验与计算结果对比(常规高超声速风洞)Fig. 9 CD~α comparison of numerical and experimental results from Φ1m hypersonic wind tunnel![]()
图 10 CL~α试验与计算结果对比(常规高超声速风洞)Fig. 10 CL~α comparison of numerical and experimental results from Φ1m hypersonic wind tunnel通过试验验证，表明数值计算具有较高的可信度。为了进一步分析规律，采用数值计算分析了脉冲燃烧风洞与常规高超声速风洞阻力系数的数据相关性，如图 11所示。
![]()
图 11 CD~α计算与试验结果对比(脉冲燃烧风洞和常规高超声速风洞)Fig. 11 CD~α comparison of numerical and experimental results from combustion heated impulse facility and hypersonic wind tunnel(1) 脉冲燃烧风洞等温壁计算得到的阻力系数比常规高超声速风洞大7%(迎角6°)~17%(迎角-2°)，迎角越小，差异越大。这个差异由多方面因素导致，主要有雷诺数、比热比、壁温比等。
(2) 针对同一个模型，分别采用脉冲燃烧风洞和常规高超声速风洞的模拟参数计算(壁面都取绝热壁条件)，此时得到的结果差异应该排除了壁温比的影响，而反映了两个风洞模拟条件(如雷诺数、比热比)不同引起的变化量，计算获得的这一部分的影响量约为3%(迎角6°)~6%(迎角-2°)。
(3) 除去两个风洞模拟条件(如雷诺数、比热比)不同引起的阻力变化量，可得到壁温比的影响量约为4%(迎角6°)~12%(迎角-2°)，迎角越小，差异越大。
5 结论
通过脉冲燃烧风洞与常规高超声速风洞不通气高超标模对比试验和计算，得到如下结论：
(1) 脉冲燃烧风洞获得的气动力系数变化规律与常规高超声速风洞一致；
(2) 雷诺数、壁温比对阻力系数均有影响，其中壁温比影响显著，脉冲燃烧风洞获得的阻力系数明显大于常规高超声速风洞；
(3) 水凝结对不通气高超标模的轴向力影响较小，对法向力影响较大，且不同迎角状态影响量不同。
-
![]()
图 8 试验段不同宽度位置流场马赫数分布（Ma=0.8）
Fig. 8 Flow field mach number distributions in different width positions
![]()
图 9 不同宽度位置流场气流角分布（Ma=0.8）
Fig. 9 Flow field flow angle distributions in different width positions
表 1 风洞主要技术参数
Table 1 Main technical parameters of the wind tunnel
序号 参 数 数 值 1 试验段截面尺寸 190 mm×445 mm 2 试验马赫数 0.3～1.8 3 试验雷诺数 3.1×105～4.5×107 m−1 4 次流温度 170～373 K 5 稳定段总压 5～300 kPa 6 稳定段总温 常温 7 叶片数 ≥7 8 典型弦长 75～120 mm 9 气流角调节范围 0°～180° 
下载: 导出CSV
表 2 三级引射器主要技术参数
Table 2 Main technical parameters of the three stages injector
级号 py/MPa m/（kg﹒s–1） May pb/kPa 1 0.25 3.96 3.6 93.6 2 0.36 16.38 3.0 80.0 3 0.96 79.60 2.8 20.0
下载: 导出CSV
-
[1] 肖洪,吴虎,廉筱纯. 雷诺数对涡扇发动机性能及稳定性影响[J]. 航空动力学报,2005,20(3):394-398. DOI: 10.3969/j.issn.1000-8055.2005.03.010 XIAO H,WU H,LIAN X C. Reynolds number effects on performance and aerodynamic stabilities of the turbofan engines[J]. Journal of Aerospace Power,2005,20(3):394-398. doi: 10.3969/j.issn.1000-8055.2005.03.010
[2] 王进,骆广琦,陶增元. 雷诺数对压气机特性及发动机稳定性影响的计算和分析[J]. 航空动力学报,2003,18(1):20-23. DOI: 10.13224/j.cnki.jasp.2003.01.004 WANG J,LUO G Q,TAO Z Y. Effects of Reynolds number on compressor performance and engine stability[J]. Journal of Aerospace Power,2003,18(1):20-23. doi: 10.13224/j.cnki.jasp.2003.01.004
[3] MAYLE R E. The 1991 IGTI scholar lecture: the role of laminar-turbulent transition in gas turbine engines[J]. Journal of Turbomachinery,1991,113(4):509-536. doi: 10.1115/1.2929110
[4] RHODEN H G. Effects of Reynolds number on the flow of air through a cascade of compressor blades[R]. RM 2919, 1956.
[5] ROBERTS W B. The effect of Reynolds number and laminar separation on axial cascade performance[J]. Journal of Engineering for Power,1975,97(2):261-273. doi: 10.1115/1.3445978
[6] IRVING A J and ROBERT O B. Aerodynamic design of axial-flow compressors[R]. NASA SP-36, 1965.
[7] GOSTELOW J P. Cascade aerodynamics[M]. Sydney: Pergamon Press, 1984: 68-69.
[8] BRUNNER S, FOTTNER L, SCHIFFER H P. Comparison of two highly loaded low pressure turbine cascades under the influence of wake-induced transition[C]//Proc of ASME Turbo Expo 2000: Power for Land, Sea, and Air. 2014. doi: 10.1115/2000-GT-0268
[9] CITAVY J,NORBURY J F. The effect of Reynolds number and turbulence intensity on the performance of a compressor cascade with prescribed velocity distribution[J]. Journal of Mechanical Engineering Science,1977,19(3):93-100. doi: 10.1243/jmes_jour_1977_019_022_02
[10] LAKE J, KING P, RIVIR R. Reduction of separation losses on a turbine blade with low Reynolds numbers[C]//Proc of the 37th Aerospace Sciences Meeting and Exhibit. 1999. doi: 10.2514/6.1999-242
[11] WUNDERWALD D, FOTTNER L. Experimental investigation of boundary layer transition and turbulence structures on a highly loaded compressor cascade[R]. 95-GT-129, 1995.
[12] KOCH C C,SMITH L H Jr. Loss sources and magnitudes in axial-flow compressors[J]. Journal of Engineering for Power,1976,98(3):411-424. doi: 10.1115/1.3446202
[13] CRAIG H R M,COX H J A. Performance estimation of axial flow turbines[J]. Proceedings of the Institution of Mechanical Engineers,1970,185(1):407-424. doi: 10.1243/pime_proc_1970_185_048_02
[14] 刘永泉,刘太秋,季路成. 航空发动机风扇/压气机技术发展的若干问题与思考[J]. 航空学报,2015,36(8):2563-2576. DOI: 10.7527/S1000-6893.2015.0078 LIU Y Q,LIU T Q,JI L C. Some problems and thoughts in the development of aero-engine fan/compressor[J]. Acta Aeronautica et Astronautica Sinica,2015,36(8):2563-2576. doi: 10.7527/S1000-6893.2015.0078
[15] SCHREIBER H A,STEINERT W,KÜSTERS B. Effects of Reynolds number and free-stream turbulence on boundary layer transition in a compressor cascade[J]. Journal of Turbomachinery,2002,124(1):1-9. doi: 10.1115/1.1413471
[16] SONODA T,SCHREIBER H A. Aerodynamic characteristics of supercritical outlet guide vanes at low Reynolds number conditions[J]. Journal of Turbomachinery,2007,129(4):694-704. doi: 10.1115/1.2720868
[17] GOMES R A,STOTZ S,BLAIM F,et al. Hot-film measurements on a low pressure turbine linear cascade with bypass transition[J]. Journal of Turbomachinery,2015,137(9):091007. doi: 10.1115/1.4029967
[18] 王文涛,王子楠,张宏武,等. 压气机静叶栅层流分离泡转捩与角区分离数值模拟与实验[J]. 航空动力学报,2017,32(9):2273-2282. WANG W T,WANG Z N,ZHANG H W,et al. Numerical simulation and experiment of laminar separation bubble transition and corner separation of compressor stator cascade[J]. Journal of Aerospace Power,2017,32(9):2273-2282.
[19] 凌代军,王晖,马昌友. 低雷诺数亚声速扩压平面叶栅试验[J]. 航空动力学报,2013,28(1):171-179. DOI: 10.13224/j.cnki.jasp.2013.01.024 LING D J,WANG H,MA C Y. Subsonic compressor plane cascade experiment at low Reynolds number[J]. Journal of Aerospace Power,2013,28(1):171-179. doi: 10.13224/j.cnki.jasp.2013.01.024
[20] 林岳峥,祝利,王海. 全球鹰无人侦察机的技术特点与应用趋势[J]. 飞航导弹,2011(9):21-24,32. [21] 王治敏,徐让书,赵长宇. 某叶栅风洞栅前流场的分析与改进[J]. 沈阳航空航天大学学报,2016,33(5):12-17. DOI: 10.3969/j.issn.2095-1248.2016.05.003 WANG Z M,XU R S,ZHAO C Y. Analysis and improvement of flow field in front of a cascade wind tunnel[J]. Journal of Shenyang Aerospace Ace University,2016,33(5):12-17. doi: 10.3969/j.issn.2095-1248.2016.05.003
[22] 廖达雄,陈吉明,彭强,等. 连续式跨声速风洞设计关键技术[J]. 实验流体力学,2011,25(4):74-78. DOI: 10.3969/j.issn.1672-9897.2011.04.014 LIAO D X,CHEN J M,PENG Q,et al. Key design techniques of the low noise continuous transonic wind tunnel[J]. Journal of Experiments in Fluid Mechanics,2011,25(4):74-78. doi: 10.3969/j.issn.1672-9897.2011.04.014
[23] 国防科学技术工业委员会. 低速风洞和高速风洞流场品质要求: GJB 1179A—2012 [S]. 北京: 总装备部标准出版发行部, 2012. [24] DUNKER R,RECHTER H,STARKEN H,et al. Redesign and performance analysis of a transonic axial compressor stator and equivalent plane cascades with subsonic controlled diffusion blades[J]. Journal of Engineering for Gas Turbines and Power,1984,106(2):279-287. doi: 10.1115/1.3239560
[25] KIOCK R,LEHTHAUS F,BAINES N C,et al. The transonic flow through a plane turbine cascade as measured in four European wind tunnels[J]. Journal of Engineering for Gas Turbines and Power,1986,108(2):277-284. doi: 10.1115/1.3239900
计量
- 文章访问数: 572
- HTML全文浏览量: 204
- PDF下载量: 72
