Wind field characteristics on a bridge site under complex mountain terrain
-
摘要: 通过地貌模型风洞试验研究复杂群山情况下某桥址的风场特性,分析平均风速、风攻角、湍流强度、湍流积分尺度等随风向角和测点位置的变化特性,获得了复杂群山环境下典型位置测点脉动风速功率谱的变化情况。研究表明:复杂群山环境下,桥址的平均风速均小于梯度风高度的风速;爬坡效应使得顺山谷方向来流(顺风向,即风从谷口吹入)产生显著的正攻角,最大值为+35.3°;顺风向时,桥址各测点的纵向湍流强度和横向湍流强度达到最小,约为10%,其他风向下的湍流强度较大;顺风向时,桥址的湍流积分尺度较大,且随测点高度上升而增大,桥址各测点功率谱较来流功率谱发生明显变化,功率谱高频段能量显著增大、单峰特征降低。Abstract: In order to obtain wind field characteristics on a bridge site on a complex mountain terrain, small-scale topographic wind-tunnel model tests are employed. The variation characteristics of the mean wind speed, wind attack angle, turbulent intensity and turbulence integral scale with respect to the wind azimuth and measuring point position are analyzed. The changes of the wind speed spectrum of typical measuring points on the complex mountain terrain are also studied. The results show that the mean wind speeds of all the measuring points on the bridge site are less than that of the gradient height. When wind blows along the valley, significant positive attack angles are produced due to the wind climbing effect, with the maximum value reaching +35.3°. The along-wind and crosswind turbulent intensities in this wind direction reach their minimum values as low as nearly 10%, and become larger in other wind directions. The turbulence integral scale in this direction is much larger than that in other directions, and increases with the increasing height of the measuring point. The wind speed power spectrum in the direction has significant differences compared with the incoming wind speed spectrum, of which the energy in the high frequency section increases significantly and the feature of the single peak diminishes.
-
Keywords:
- terrain model /
- wind tunnel test /
- mean wind speed /
- turbulent intensity /
- turbulence integral scale
-
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随迎角α变化的试验结果曲线。
试验表明:在总温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试验结果,同时也给出了两个模型按照对应风洞来流条件计算获得的气动力系数结果。
计算结果表明:模型尺度对不通气标模的轴向力有一定影响,小不通气标模的轴向力略大于大不通气标模,相差约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的对比曲线。
常规高超声速风洞试验和计算结果表明:不同壁温条件对计算结果影响不大,计算结果与试验符合较好,阻力误差在4%以内,升力误差在3%以内。图 9和10给出了相应的阻力系数和升力系数比较曲线。
通过试验验证,表明数值计算具有较高的可信度。为了进一步分析规律,采用数值计算分析了脉冲燃烧风洞与常规高超声速风洞阻力系数的数据相关性,如图 11所示。
(1) 脉冲燃烧风洞等温壁计算得到的阻力系数比常规高超声速风洞大7%(迎角6°)~17%(迎角-2°),迎角越小,差异越大。这个差异由多方面因素导致,主要有雷诺数、比热比、壁温比等。
(2) 针对同一个模型,分别采用脉冲燃烧风洞和常规高超声速风洞的模拟参数计算(壁面都取绝热壁条件),此时得到的结果差异应该排除了壁温比的影响,而反映了两个风洞模拟条件(如雷诺数、比热比)不同引起的变化量,计算获得的这一部分的影响量约为3%(迎角6°)~6%(迎角-2°)。
(3) 除去两个风洞模拟条件(如雷诺数、比热比)不同引起的阻力变化量,可得到壁温比的影响量约为4%(迎角6°)~12%(迎角-2°),迎角越小,差异越大。
5 结论
通过脉冲燃烧风洞与常规高超声速风洞不通气高超标模对比试验和计算,得到如下结论:
(1) 脉冲燃烧风洞获得的气动力系数变化规律与常规高超声速风洞一致;
(2) 雷诺数、壁温比对阻力系数均有影响,其中壁温比影响显著,脉冲燃烧风洞获得的阻力系数明显大于常规高超声速风洞;
(3) 水凝结对不通气高超标模的轴向力影响较小,对法向力影响较大,且不同迎角状态影响量不同。
-
[1] 中华人民共和国交通运输部. 公路桥梁抗风设计规范: JTG/T 3360-01—2018[S]. 北京: 人民交通出版社, 2018. Ministry of Transport of the People’s Republic of China. Wind-resistant design specification for highway bridges: JTG/T 3360-01—2018[S]. Beijing: China Communications Press, 2018.
[2] CHOCK G Y K,COCHRAN L. Modeling of topographic wind speed effects in Hawaii[J]. Journal of Wind Engineering and Industrial Aerodynamics,2005,93(8):623-638. doi: 10.1016/j.jweia.2005.06.002
[3] MERONEY R N. Wind-tunnel simulation of the flow over hills and complex terrain[J]. Journal of Wind Engineering and Industrial Aerodynamics,1980,5(3-4):297-321. doi: 10.1016/0167-6105(80)90039-2
[4] TEUNISSEN H W,SHOKR M E,BOWEN A J,et al. The Askervein Hill Project: Wind-tunnel simulations at three length scales[J]. Boundary-Layer Meteorology,1987,40(1-2):1-29. doi: 10.1007/BF00140067
[5] LUBITZ W D,WHITE B R. Wind-tunnel and field investigation of the effect of local wind direction on speed-up over hills[J]. Journal of Wind Engineering and Industrial Aerodynamics,2007,95(8):639-661. doi: 10.1016/j.jweia.2006.09.001
[6] CARPENTER P,LOCKE N. Investigation of wind speeds over multiple two-dimensional hills[J]. Journal of Wind Engineering and Industrial Aerodynamics,1999,83(1-3):109-120. doi: 10.1016/S0167-6105(99)00065-3
[7] 沈国辉,翁文涛,王轶文,等. 某复杂山体的三维风场特征研究[J]. 振动与冲击,2020,39(4):75-80. DOI: 10.13465/j.cnki.jvs.2020.04.009 SHEN G H,WENG W T,WANG Y W,et al. A study on three-dimensional wind field characteristics of a complex hill[J]. Journal of Vibration and Shock,2020,39(4):75-80. doi: 10.13465/j.cnki.jvs.2020.04.009
[8] 陈政清,李春光,张志田,等. 山区峡谷地带大跨度桥梁风场特性试验[J]. 实验流体力学,2008,22(3):54-59, 67. CHEN Z Q,LI C G,ZHANG Z T,et al. Model test study of wind field characteristics of long-span bridge site in mountainous valley terrain[J]. Journal of Experiments in Fluid Mechanics,2008,22(3):54-59, 67.
[9] 庞加斌,宋锦忠,林志兴. 山区峡谷桥梁抗风设计风速的确定方法[J]. 中国公路学报,2008,21(5):39-44. DOI: 10.19721/j.cnki.1001-7372.2008.05.008 PANG J B,SONG J Z,LIN Z X. Determination method for wind-resistant design wind speed of mountainous-valley bridge[J]. China Journal of Highway and Transport,2008,21(5):39-44. doi: 10.19721/j.cnki.1001-7372.2008.05.008
[10] 刘黎阳,张志田,汪志雄,等. 不同规模地形模型对某山区桥梁设计风特性确定的影响[J]. 实验流体力学,2018,32(6):49-54. DOI: 10.11729/syltlx20170140 LIU L Y,ZHANG Z T,WANG Z X,et al. Scope effects of terrain models on wind properties design of a bridge located at mountainous area[J]. Journal of Experiments in Fluid Mechanics,2018,32(6):49-54. doi: 10.11729/syltlx20170140
[11] 张宏杰,赵金飞,蔡达章,等. 垭口地貌要素对风速分布规律影响的风洞试验研究[J]. 实验流体力学,2014,28(4):25-30. DOI: 10.11729/syltlx20130044 ZHANG H J,ZHAO J F,CAI D Z,et al. Wind tunnel test on the influence of col features on wind speed distribution[J]. Journal of Experiments in Fluid Mechanics,2014,28(4):25-30. doi: 10.11729/syltlx20130044
[12] 庞加斌,宋锦忠,林志兴. 四渡河峡谷大桥桥位风的湍流特性实测分析[J]. 中国公路学报,2010,23(3):42-47. DOI: 10.19721/j.cnki.1001-7372.2010.03.007 PANG J B,SONG J Z,LIN Z X. Field measurement analysis of wind turbulence characteristics of Sidu river valley bridge site[J]. China Journal of Highway and Transport,2010,23(3):42-47. doi: 10.19721/j.cnki.1001-7372.2010.03.007
[13] HUI M C H,LARSEN A,XIANG H F. Wind turbulence characteristics study at the Stonecutters Bridge site: Part I—Mean wind and turbulence intensities[J]. Journal of Wind Engineering and Industrial Aerodynamics,2009,97(1):22-36. doi: 10.1016/j.jweia.2008.11.002
[14] HUI M C H,LARSEN A,XIANG H F. Wind turbulence characteristics study at the Stonecutters Bridge site: Part Ⅱ: Wind power spectra, integral length scales and coherences[J]. Journal of Wind Engineering and Industrial Aerodynamics,2009,97(1):48-59. doi: 10.1016/j.jweia.2008.11.003
[15] 于舰涵,李明水,廖海黎. 山区地形对桥位风场影响的数值模拟[J]. 西南交通大学学报,2016,51(4):654-662. DOI: 10.3969/j.issn.0258-2724.2016.04.008 YU J H,LI M S,LIAO H L. Numerical simulation of effect of mountainous topography on wind field at bridge site[J]. Journal of Southwest Jiaotong University,2016,51(4):654-662. doi: 10.3969/j.issn.0258-2724.2016.04.008
[16] HUANG G Q,CHENG X,PENG L L,et al. Aerodynamic shape of transition curve for truncated mountainous terrain model in wind field simulation[J]. Journal of Wind Engineering and Industrial Aerodynamics,2018,178:80-90. doi: 10.1016/j.jweia.2018.05.008
[17] 胡朋,李永乐,廖海黎. 山区峡谷桥址区地形模型边界过渡段形式研究[J]. 空气动力学学报,2013,31(2):231-238. DOI: 10.7638/kqdlxxb-2011.0184 HU P,LI Y L,LIAO H L. Shape of boundary transition section for mountains-gorge bridge site terrain model[J]. Acta Aerodynamica Sinica,2013,31(2):231-238. doi: 10.7638/kqdlxxb-2011.0184
[18] 余世策,陈勇,李庆祥,等. 建筑风环境风洞试验中风速探头的研制与应用[J]. 实验流体力学,2013,27(4):83-87. DOI: 10.3969/j.issn.1672-9897.2013.04.015 YU S C,CHEN Y,LI Q X,et al. Development and application of wind speed probe for wind tunnel test of wind environment around buildings[J]. Journal of Experiments in Fluid Mechanics,2013,27(4):83-87. doi: 10.3969/j.issn.1672-9897.2013.04.015
[19] 中华人民共和国住房和城乡建设部. 建筑工程风洞试验方法标准: JGJ/T 338—2014[S]. 北京: 中国建筑工业出版社, 2015. Ministry of Housing and Urban-Rural Development of the People's Republic of China. Standard for wind tunnel test of buildings and structures: JGJ/T 338—2014[S]. Beijing: China Architecture & Building Press, 2015.
[20] SIMIU E, SCANLAN R H. Wind effects on structure[M]. 3rd ed. New York: Wiley InterScience, 1996.
[21] FLAY R G J,STEVENSON D C. Integral length scales in strong winds below 20 m[J]. Journal of Wind Engineering and Industrial Aerodynamics,1988,28(1-3):21-30. doi: 10.1016/0167-6105(88)90098-0