重气体介质中跨声速非定常气动特性

刘永平; 欧阳炎; 孟少飞; 寇西平; 夏洪亚

doi:10.13224/j.cnki.jasp.20240470

重气体介质中跨声速非定常气动特性

doi: 10.13224/j.cnki.jasp.20240470

刘永平^1,,
欧阳炎¹,
孟少飞¹,
寇西平^{1, 2},
夏洪亚^1,*, ,

1.
中国空气动力研究与发展中心高速空气动力研究所，四川绵阳 621000
2.
西北工业大学航空学院，西安 710072

基金项目: 四川省自然科学基金（2023NSFSC0400）

详细信息

作者简介:
刘永平（1991-），男，博士生，主要从事飞行器气动弹性力学研究。E-mail：liuyp917@163.com

通讯作者:
夏洪亚（1981-），女，工程师，硕士，主要从事飞行器气动弹性试验研究。E-mail：421988039@qq.com

中图分类号: V211.24
计量
- 文章访问数: 554
- HTML浏览量: 270
- PDF量: 47
- 被引次数: 0
出版历程
- 收稿日期: 2024-07-12
- 网络出版日期: 2025-02-18

Transonic unsteady aerodynamic characteristics in heavy gas medium

LIU Yongping^1
,,
OUYANG Yan¹,
MENG Shaofei¹,
KOU Xiping^{1, 2},
XIA Hongya^{1,*
, ,}

1.
High Speed Aerodynamics Institute，China Aerodynamics Research and Development Center，Mianyang Sichuan 621000，China
2.
School of Aeronautics，Northwestern Polytechnical University，Xi’an 710072，China

摘要

摘要:
采用雷诺平均的Navier-Stokes方程（RANS）和Spalart-Allmaras一方程湍流模型（SA）模拟了NACA64A010翼型在空气和重气体介质R-134a中的非定常跨声速气动特性。在相同马赫数、雷诺数及减缩频率条件下，俯仰简谐运动的计算结果表明：重气体介质中翼型表面非定常压力系数的模值和相位的分布在激波位置与空气中差异明显，非定常升力系数与空气介质中差异不大，俯仰力矩系数在幅值和相位上都与空气中存在一定差异，将跨声速相似律应用于非定常气动力修正，俯仰力矩系数的幅值和相位实现了到空气的相似性转化，但随着减缩频率的增大，跨声速相似律的修正效果变差。通过俯仰力矩做功分析表明，如果对重气体中的非定常气动力不进行相似性修正，将使空气和重气体中颤振特性出现差异，影响重气体风洞颤振特性的评估。该研究为后续重气体介质中飞行器颤振特性研究及修正方法的发展提供了基础支持。
- 重气体介质 /
- 跨声速 /
- 非定常气动力 /
- 俯仰简谐运动 /
- 颤振
Abstract:
The unsteady transonic aerodynamic characteristics of a NACA64A010 airfoil in air and heavy gas medium R-134a were simulated by using the Reynolds-averaged Navier-Stokes equations (RANS) and the Spalart-Allmaras (SA) one-equation turbulence model. Under the same Mach number, Reynolds number and reduced frequency, the calculation results of simple harmonic motion in pitch showed that the distribution of modulus and phase of unsteady pressure coefficient in heavy gas medium was obviously different from that in air, the unsteady lift coefficient was not different from that in air medium, and the amplitude and phase of pitch moment coefficient were different from those in air medium. The transonic similarity law was applied to the unsteady aerodynamic force correction, and the amplitude and phase of pitch moment coefficient were transformed into air similarity, but with the increase of reduced frequency, the correction effect of transonic similarity law became worse. Through analysis of pitching moment work, it was shown that if the unsteady aerodynamic force in heavy gas was not corrected by similarity, the flutter characteristics in air and heavy gas could be different, thus affecting the evaluation of flutter characteristics in heavy gas wind tunnel. This study can provide a basic support for the follow-up research on flutter characteristics of aircraft in heavy gas media and the development of correction methods.
- heavy gas medium /
- transonic /
- unsteady aerodynamic force /
- pitching simple harmonic motion /
- flutter

HTML

图 1 模型和边界条件示意图

Figure 1. Diagram of model and boundary conditions

下载: 全尺寸图片幻灯片

图 2 翼型周边网格分布

Figure 2. Grid distribution around airfoil

下载: 全尺寸图片幻灯片

图 3 不同网格定常压力分布对比

Figure 3. Comparison of steady pressure distribution with different grids

下载: 全尺寸图片幻灯片

图 4 不同网格上表面非定常压力分布对比

Figure 4. Comparison of unsteady pressure distribution on the upper surface with different grids

下载: 全尺寸图片幻灯片

图 5 不同网格非定常气动力对比

Figure 5. Comparison of unsteady aerodynamic forces with different grids

下载: 全尺寸图片幻灯片

图 6 文献[25]C_m-α结果

Figure 6. C_m-α results from Ref.[25]

下载: 全尺寸图片幻灯片

图 7 不同介质中定常压力分布对比

Figure 7. Comparison of steady pressure distribution in different media

下载: 全尺寸图片幻灯片

图 8 不同介质中上表面非定常压力分布对比

Figure 8. Comparison of unsteady pressure distribution on the upper surface in different media

下载: 全尺寸图片幻灯片

图 9 不同介质中非定常气动力对比

Figure 9. Comparison of unsteady aerodynamic forces in different media

下载: 全尺寸图片幻灯片

图 10 减缩频率k=0.100的非定常气动力

Figure 10. Unsteady aerodynamic force with reduced frequency k=0.100

下载: 全尺寸图片幻灯片

图 11 减缩频率k=0.202的非定常气动力

Figure 11. Unsteady aerodynamic force with reduced frequency k=0.202

下载: 全尺寸图片幻灯片

图 12 减缩频率k=0.300的非定常气动力

Figure 12. Unsteady aerodynamic force with reduced frequency k=0.300

下载: 全尺寸图片幻灯片

表 1 计算状态

Table 1. Calculation status

介质	Ma	α₀/（°）	α_m/（°）	Re/10⁶	X_a/c
空气	0.7960	−0.21	1.01	25	0.248
R-134a修正前	0.7960
R-134a修正后	0.8095

下载: 导出CSV

表 2 俯仰力矩系数特性和做功

Table 2. Pitching moment coefficient characteristics and work

k	介质	f	φ/rad	W
0.100	空气	0.0067	−2.55	−0.01186
	R-134a修正前	0.0055	−2.22	−0.01390
	R-134a修正后	0.0069	−2.59	−0.01147
0.202	空气	0.0108	−2.29	−0.02578
	R-134a修正前	0.0101	−2.03	−0.02873
	R-134a修正后	0.0109	−2.35	−0.02461
0.300	空气	0.0145	−2.33	−0.03337
	R-134a修正前	0.0150	−2.07	−0.04179
	R-134a修正后	0.0140	−2.41	−0.02968

下载: 导出CSV

参考文献(31)

[1]	LIVNE E. Future of airplane aeroelasticity[J]. Journal of Aircraft, 2003, 40(6): 1066-1092. doi: 10.2514/2.7218
[2]	DOWELL E H. Unsteady transonic aerodynamics and aeroelasticity[J]. Applied Mechanics Reviews, 1988, 41(8): 299-319. doi: 10.1115/1.3151909
[3]	COLE S, GARCIA J. Past, present, and future capabilities of the Transonic Dynamics Tunnel from an aeroelasticity perspective[R]. AIAA 2000-1767, 2000.
[4]	罗务揆, 谭申刚, 谢怀强, 等. 确定颤振模型设计参数的方法研究[J]. 航空学报, 2013, 34(10): 2383-2390. LUO Wukui, TAN Shengang, XIE Huaiqiang, et al. Research on methods used to determine flutter model design factors[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(10): 2383-2390. (in Chinese LUO Wukui, TAN Shengang, XIE Huaiqiang, et al. Research on methods used to determine flutter model design factors[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(10): 2383-2390. (in Chinese)
[5]	钱卫, 张桂江, 刘钟坤. 飞机全动平尾颤振特性风洞试验[J]. 航空学报, 2015, 36(4): 1093-1102. QIAN Wei, ZHANG Guijiang, LIU Zhongkun. Flutter characteristics for aircraft all-movable horizontal tail through wind tunnel test[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(4): 1093-1102. (in Chinese QIAN Wei, ZHANG Guijiang, LIU Zhongkun. Flutter characteristics for aircraft all-movable horizontal tail through wind tunnel test[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(4): 1093-1102. (in Chinese)
[6]	IVANCO T G. Unique testing capabilities of the NASA langley transonic dynamics tunnel, an exercise in aeroelastic scaling[R]. AIAA 2013-2625, 2013.
[7]	YEAGER W T, KVATERNIK R G. A historical overview of aeroelasticity branch and transonic dynamics tunnel contributions to rotorcraft technology and development[R]. NASA Technical Memorandum 2001-211054, 2001.
[8]	YEAGER W Jr, KVATERNIK R. Contributions of the Langley Transonic Dynamics Tunnel to rotorcraft technology and development[R]. AIAA 2000-1771, 2000.
[9]	COLE S, KELLER D, PIATAK D. Contributions of the NASA langley transonic dynamics tunnel to launch vehicle and spacecraft development[R]. AIAA 2000-1772, 2000.
[10]	PERR B, NOLL T, SCOTT R. Contributions of the transonic dynamics tunnel to the testing of active control of aeroelastic response[R]. AIAA 2000-1769, 2000.
[11]	HUBER P W. Use of Freon-12 as a fluid for aerodynamic testing[R]. NACA-TN-1024, 1946.
[12]	VON DOENHOFF A E, BRASLOW A L, SCHWARTZBERG M A. Studies of the use of Freon-12 as a wind tunnel testing medium[R]. NACA-TN-3000, 1953.
[13]	WELLER W. Comparison of aerodynamic data measured in air and Freon-12 wind-tunnel test mediums[R]. NASA-TM-78671, 1978.
[14]	ANDERSON W. A numerical study on the use of sulfur hexafluoride as a test gas for wind tunnels[R]. AIAA 1990-1421, 1990.
[15]	ANDERS J, ANDERSON W, MURTHY A. The use of heavy gas for increased Reynolds numbers in transonic wind tunnels[R]. AIAA 1998-2882, 1998.
[16]	刘永平, 寇西平, 查俊, 等. 重气体介质的等熵流动特性[J]. 航空动力学报, 2025, 40(9): 20230421. LIU Yongping, KOU Xiping, ZHA Jun, et al. The isentropic flow characteristics of heavy gas medium[J]. Journal of Aero-space Power, 2025, 40(9): 20230421. (in Chinese LIU Yongping, KOU Xiping, ZHA Jun, et al. The isentropic flow characteristics of heavy gas medium[J]. Journal of Aero-space Power, 2025, 40（9）: 20230421. (in Chinese)
[17]	张汇卓. 重气体环境机翼气动特性数值模拟研究[D]. 北京: 军事科学院, 2021. ZHANG Huizhuo. Numerical simulation study on aerodynamic characterristics of wing in heavy gas environment[D]. Beijing: Academy of Military Sciences, 2021. (in Chinese ZHANG Huizhuo. Numerical simulation study on aerodynamic characterristics of wing in heavy gas environment[D]. Beijing: Academy of Military Sciences, 2021. (in Chinese)
[18]	查俊, 曾开春, 寇西平, 等. 重气体介质中超临界翼型跨声速流动特性[J]. 航空动力学报, 2021, 36(9): 1894-1905. ZHA Jun, ZENG Kaichun, KOU Xiping, et al. Transonic flow characteristics of supercritical airfoil in heavy gas medium[J]. Journal of Aerospace Power, 2021, 36(9): 1894-1905. (in Chinese ZHA Jun, ZENG Kaichun, KOU Xiping, et al. Transonic flow characteristics of supercritical airfoil in heavy gas medium[J]. Journal of Aerospace Power, 2021, 36(9): 1894-1905. (in Chinese)
[19]	刘永平, 查俊, 胡哲, 等. 重气体介质中气动特性的相关性修正[J]. 航空动力学报, 2025, 40(9): 20230748. LIU Yongping, KOU Xiping, ZHA Jun, et al. The isentropic flow characteristics of heavy gas medium[J]. Journal of Aerospace Power, 2025, 40(9): 20230748. (in Chinese LIU Yongping, KOU Xiping, ZHA Jun, et al. The isentropic flow characteristics of heavy gas medium[J]. Journal of Aerospace Power, 2025, 40（9）: 20230748. (in Chinese)
[20]	BENDIKSEN O. Role of shock dynamics in transonic flutter[R]. AIAA 1992-2121, 1992.
[21]	DAVIS S S, MALCOLM G N. Transonic shock-wave/boundary-layer interactions on an oscillating airfoil[J]. AIAA Journal, 1980, 18(11): 1306-1312. doi: 10.2514/3.50886
[22]	BENDIKSEN O O. Review of unsteady transonic aerodynamics: theory and applications[J]. Progress in Aerospace Sciences, 2011, 47(2): 135-167. doi: 10.1016/j.paerosci.2010.07.001
[23]	AGARD. Compendium of unsteady aerodynamic measurements[R]. AGARD Report No. 702, 1982.
[24]	DAVIS S S. NACA64A010 (NACA Ames Model) oscillatory pitching[R]. AGARD Report No. 702, 1982.
[25]	LIU Feng, JI Shanhong. Unsteady flow calculations with a multigrid Navier-Stokes method[J]. AIAA Journal, 1996, 34(10): 2047-2053. doi: 10.2514/3.13351
[26]	WANG Baoyuan, ZHA Gecheng. Numerical simulation of transonic limit cycle oscillations using high-order low-diffusion schemes[J]. Journal of Fluids and Structures, 2010, 26(4): 579-601. doi: 10.1016/j.jfluidstructs.2010.02.003
[27]	YUAN Weixing, SANDHU R, DIAS de MATOS Jr O, et al. Methodology development for coupled aeroelastic analysis of wing flutter[R]. AIAA 2016-1550, 2016.
[28]	刘南. 机翼跨声速非线性颤振及高效分析方法研究[D]. 西安: 西北工业大学, 2016. LIU Nan. Study on transonic nonlinear flutter and efficient analysis method of wing[D]. Xi’an: Northwestern Polytechnical University, 2016. (in Chinese LIU Nan. Study on transonic nonlinear flutter and efficient analysis method of wing[D]. Xi’an: Northwestern Polytechnical University, 2016. (in Chinese)
[29]	REED. Correlation with flight of some aeroelastic model studies in the NASA Langley transonic dynamics tunnel[R]. NASA-SP-415, 1976.
[30]	WANG Shengyi, INGHAM D B, MA Lin, et al. Numerical investigations on dynamic stall of low Reynolds number flow around oscillating airfoils[J]. Computers & Fluids, 2010, 39(9): 1529-1541.
[31]	LU K, XIE Y H, ZHANG D, et al. Numerical investigations into the asymmetric effects on the aerodynamic response of a pitching airfoil[J]. Journal of Fluids and Structures, 2013, 39: 76-86. doi: 10.1016/j.jfluidstructs.2013.02.001