CN106768818A - 一种激波风洞中混合气体来流运行参数获得方法 - Google Patents

一种激波风洞中混合气体来流运行参数获得方法 Download PDF

Info

Publication number
CN106768818A
CN106768818A CN201611245685.1A CN201611245685A CN106768818A CN 106768818 A CN106768818 A CN 106768818A CN 201611245685 A CN201611245685 A CN 201611245685A CN 106768818 A CN106768818 A CN 106768818A
Authority
CN
China
Prior art keywords
specific heat
mixed gas
pressure
heat ratio
gas
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201611245685.1A
Other languages
English (en)
Other versions
CN106768818B (zh
Inventor
王丹
卢洪波
谌君谋
刘明
林键
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Aerospace Aerodynamics CAAA
Original Assignee
China Academy of Aerospace Aerodynamics CAAA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Academy of Aerospace Aerodynamics CAAA filed Critical China Academy of Aerospace Aerodynamics CAAA
Priority to CN201611245685.1A priority Critical patent/CN106768818B/zh
Publication of CN106768818A publication Critical patent/CN106768818A/zh
Application granted granted Critical
Publication of CN106768818B publication Critical patent/CN106768818B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种激波风洞中混合气体来流运行参数获得方法,在激波风洞试验中,来流气体比热比是重要的物理量,对于空气,比热比随温度变化很小,可以获得来流马赫数,而对于混合气体,其难点在于比热比随温度变化可能较大,难以确定来流马赫数。在数值模拟中,一般取气体的比热比为定值,即等效比热比,这个取值对于流场参数计算非常关键,但该值的大小在文献和资料中差别很大,难以确定,本发明用试验和理论分析相结合的方法,获得了试验气体的等效比热比,为特殊气体来流试验研究,尤其是深空探测试验方法提供依据。

Description

一种激波风洞中混合气体来流运行参数获得方法
技术领域
本发明涉及一种激波风洞中混合气体来流运行参数获得方法,用于深空探测领域。
背景技术
进入21世纪,世界各航天国家纷纷推出新的深空探测发展战略和规划,其中包括多个向地外天体发射着陆器的项目。而由于大多数地外天体大气环境与气动加热原理同地球存在较大的差异性(组分、密度不同),因此对再入探测器的表面气动热防护和气动力变轨控制等一系列问题的研究提出了全新的挑战。
我国还未对特殊气体来流气动特性开展过***研究,尚缺乏对该类问题的明确认知,相关的气动力和气动热的试验技术手段也比较匮乏。而开展该方向研究的前提是确定特殊气体及混合气体来流参数,该部分内容国外文献也没有详细说明。
发明内容
本发明所要解决的技术问题:针对现有技术的不足,提出一种激波风洞中混合气体来流运行参数获得方法,获得试验气体的等效比热比,进而获得混合气体来流参数,为深空探测试验提供依据。
本发明的技术解决方案:
一种激波风洞中混合气体来流运行参数获得方法,包括如下步骤:
步骤1:在低压被驱动段前后夹好隔膜,使低压被驱动段形成封闭的空间,将n种混合气体分别n次充入低压被驱动段中,每次充入混合气体的质量分数不同,采集每次的来流总压Pe及皮托压力Ps;
步骤2:根据如下方程,确定n次充入的混合气体对应的马赫数Ma和比热比r1
其中,A为喷管出口面积,Acr为喷管喉道面积;
步骤3:将步骤2中得到的n次充入的混合气体对应的比热比r1代入如下方程,解得n种混合气体中每种气体的比热比ri
其中,cpi为定压比热容,cvi为定容比热容,mi为混合气体中每种气体的质量分数;
步骤4:验证比热比ri的正确性:选用n+1种混合比例的n种气体,返回步骤1,解得n+1种比例的n种混合气体中每种气体的比热比ri′,若ri′与ri的偏差在2%以内,则认为比热比ri准确,否则,调整的比值,返回步骤1重新计算、验证,直至将偏差控制在2%以内为止。
激波风洞包括高压驱动段、低压被驱动段、第一隔膜、第二隔膜、喷管、试验段、测压钯和压力测量传感器,高压驱动段用于充入高压空气,第一隔膜设置在高压驱动段与低压被驱动段之间,第二隔膜设置在低压被驱动段与喷管之间,测压钯设置在试验段内,测压钯上设置多个压力测量传感器,用于测量均匀区的皮托压力Ps,喷管入口处的壁面上安装压力测量传感器,测量来流总压Pe。
混合气体质量分数的确定方法为:混合气体各个组分间不发生化学反应,将低压被驱动段抽吸到真空状态,然后按照公式充入不同分压的各种气体,获得所需质量分数的混合气体,
其中P1为被驱动低压段的压力,Pi为每种组分气体的分压,Xi为每种组分气体的摩尔分数。
所述步骤3中的值按照公式根据空气来流流场校测所得马赫数Ma和空气比热比r反推得到。
本发明的有益效果:
(1)在激波风洞试验中,来流气体比热比是重要的物理量,对于空气,比热比随温度变化很小,可以获得来流马赫数,而对于混合气体,其难点在于比热比随温度变化可能较大,难以确定来流马赫数。在数值模拟中,一般取气体的比热比为定值,即等效比热比,这个取值对于流场参数计算非常关键,但该值的大小在文献和资料中差别很大,难以确定,本发明用试验和理论分析相结合的方法,获得了试验气体的等效比热比,为特殊气体来流试验研究,尤其是深空探测试验方法提供依据;
(2)本发明详细说明了在激波风洞中,以混合气体为来流时,参数获得方法,可以成为地外天体(如火星)再入气动理论、计算方法和再入飞行器气动设计的重要地面试验验证手段。
附图说明
图1是本发明流程图;
图2为本发明激波风洞结构图。
具体实施方式
下面将结合附图和具体实施对本发明做进一步详细的说明。
一种激波风洞中混合气体来流运行参数获得方法,如图1所示,包括如下步骤:
步骤1:在低压被驱动段前后夹好隔膜,使低压被驱动段形成封闭的空间,将n种混合气体分别n次充入低压被驱动段中,每次充入混合气体的质量分数不同,采集每次的来流总压Pe及皮托压力Ps;
步骤2:根据如下方程,确定n次充入的混合气体对应的马赫数Ma和比热比r1
其中,A为喷管出口面积,Acr为喷管喉道面积;
步骤3:将步骤2中得到的n次充入的混合气体对应的比热比r1代入如下方程,解得n种混合气体中每种气体的比热比ri
其中,cpi为定压比热容,cvi为定容比热容,mi为混合气体中每种气体的质量分数;
步骤4:验证比热比ri的正确性:选用n+1种混合比例的n种气体,返回步骤(1),解得n+1种比例的n种混合气体中每种气体的比热比ri′,若ri′与ri的偏差在2%以内,则认为比热比ri准确,否则,调整的比值,返回步骤1重新计算、验证,直至将偏差控制在2%以内为止。
步骤1中混合气体质量分数的确定方法为:混合气体各个组分间不发生化学反应,将低压被驱动段抽吸到真空状态,然后按照公式充入不同分压的各种气体,获得所需质量分数的混合气体,
其中P1为被驱动低压段的压力,Pi为每种组分气体的分压,Xi为每种组分气体的摩尔分数。
所述步骤3中的值按照公式根据空气来流流场校测所得马赫数Ma和空气比热比r反推得到。
如图2所示,激波风洞包括高压驱动段4、低压被驱动段1、第一隔膜、第二隔膜、喷管、试验段、测压钯和压力测量传感器,高压驱动段4用于充入高压空气,第一隔膜设置在高压驱动段4与低压被驱动段1之间,第二隔膜设置在低压被驱动段1与喷管之间,测压钯设置在试验段内,测压钯上设置多个压力测量传感器,用于测量均匀区的皮托压力Ps,喷管入口处的壁面上安装压力测量传感器,测量来流总压Pe。
对于激波风筒,有激波运行方式和炮风洞运行方式两种:
其中,P51为5区压力/1区压力,Ms为激波马赫数,T51为5区温度/1区温度,T5和P5为激波运行方式时的来流总温和总压,Te和Pe为炮风洞运行方式时的总温和总压,
利用Ms,求得P5、T5(T1一般为288.16K)及Pe、Te,即激波运行方式或炮风洞运行方式时的来流总压和总温,进而通过常规计算,可以获得马赫数及雷诺数等参数。
在激波风洞试验中,来流气体比热比是重要的物理量,对于空气,比热比随温度变化很小,可以获得来流马赫数,而对于混合气体,其难点在于比热比随温度变化可能较大,难以确定来流马赫数。在数值模拟中,一般取气体的比热比为定值,即等效比热比,这个取值对于流场参数计算非常关键,但该值的大小在文献和资料中差别很大,难以确定,本发明用试验和理论分析相结合的方法,获得了试验气体的等效比热比,为特殊气体来流试验研究,尤其是深空探测试验方法提供依据。
本发明说明书中未公开内容为本领域技术人员公知常识。

Claims (4)

1.一种激波风洞中混合气体来流运行参数获得方法,其特征在于,包括如下步骤:
步骤1:在低压被驱动段前后夹好隔膜,使低压被驱动段形成封闭的空间,将n种混合气体分别n次充入低压被驱动段中,每次充入混合气体的质量分数不同,采集每次的来流总压Pe及皮托压力Ps;
步骤2:根据如下方程,确定n次充入的混合气体对应的马赫数Ma和比热比r1
P s P e = [ ( r 1 + 1 ) Ma 2 ( r 1 - 1 ) Ma 2 + 2 ] r 1 r 1 - 1 [ r 1 + 1 2 r 1 Ma 2 - ( r 1 - 1 ) ] 1 r 1 - 1 A A c r = 1 M a [ ( 1 + r 1 - 1 2 Ma 2 ) ( 2 r 1 + 1 ) ] r 1 + 1 2 ( r 1 - 1 )
其中,A为喷管出口面积,Acr为喷管喉道面积;
步骤3:将步骤2中得到的n次充入的混合气体对应的比热比r1代入如下方程,解得n种混合气体中每种气体的比热比ri
r 1 = Σ n i = 1 m i c p i Σ n i = 1 m i c v i = Σ n i = 1 r i m i c v i Σ n i = 1 m i c v i
其中,cpi为定压比热容,cvi为定容比热容,mi为混合气体中每种气体的质量分数;
步骤4:验证比热比ri的正确性:选用n+1种混合比例的n种气体,返回步骤1,解得n+1种比例的n种混合气体中每种气体的比热比ri′,若ri′与ri的偏差在2%以内,则认为比热比ri准确,否则,调整的比值,返回步骤1重新计算、验证,直至将偏差控制在2%以内为止。
2.如权利要求1所述的一种激波风洞中混合气体来流运行参数获得方法,其特征在于,激波风洞包括高压驱动段(4)、低压被驱动段(1)、第一隔膜、第二隔膜、喷管、试验段、测压钯和压力测量传感器,高压驱动段(4)用于充入高压空气,第一隔膜设置在高压驱动段(4)与低压被驱动段(1)之间,第二隔膜设置在低压被驱动段(1)与喷管之间,测压钯设置在试验段内,测压钯上设置多个压力测量传感器,用于测量均匀区的皮托压力Ps,喷管入口处的壁面上安装压力测量传感器,测量来流总压Pe。
3.如权利要求1所述的一种激波风洞中混合气体来流运行参数获得方法,其特征在于,混合气体质量分数的确定方法为:混合气体各个组分间不发生化学反应,将低压被驱动段抽吸到真空状态,然后按照公式充入不同分压的各种气体,获得所需质量分数的混合气体,
其中P1为被驱动低压段的压力,Pi为每种组分气体的分压,Xi为每种组分气体的摩尔分数。
4.如权利要求1所述的一种激波风洞中混合气体来流运行参数获得方法,其特征在于,所述步骤3中的值按照公式根据空气来流流场校测所得马赫数Ma和空气比热比r反推得到。
CN201611245685.1A 2016-12-29 2016-12-29 一种激波风洞中混合气体来流运行参数获得方法 Active CN106768818B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611245685.1A CN106768818B (zh) 2016-12-29 2016-12-29 一种激波风洞中混合气体来流运行参数获得方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611245685.1A CN106768818B (zh) 2016-12-29 2016-12-29 一种激波风洞中混合气体来流运行参数获得方法

Publications (2)

Publication Number Publication Date
CN106768818A true CN106768818A (zh) 2017-05-31
CN106768818B CN106768818B (zh) 2019-04-30

Family

ID=58928824

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611245685.1A Active CN106768818B (zh) 2016-12-29 2016-12-29 一种激波风洞中混合气体来流运行参数获得方法

Country Status (1)

Country Link
CN (1) CN106768818B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107806977A (zh) * 2017-11-29 2018-03-16 中国航空工业集团公司沈阳空气动力研究所 一种组合式宽马赫数高焓脉冲风洞管体结构
CN110057537A (zh) * 2019-04-12 2019-07-26 北京空天技术研究所 飞行器气动性能影响预测方法
CN111346528A (zh) * 2020-04-29 2020-06-30 中国空气动力研究与发展中心高速空气动力研究所 一种变比热比冷喷流混合气体的制备方法
CN115828470A (zh) * 2023-01-10 2023-03-21 中国空气动力研究与发展中心计算空气动力研究所 一种喷管型线设计方法、装置、设备及介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104374543A (zh) * 2014-11-06 2015-02-25 沈阳黎明航空发动机(集团)有限责任公司 一种研究气流混合特性的装置
CN205670075U (zh) * 2016-04-28 2016-11-02 贵州航空发动机研究所 一种飞机进气道出口流场测量机构

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104374543A (zh) * 2014-11-06 2015-02-25 沈阳黎明航空发动机(集团)有限责任公司 一种研究气流混合特性的装置
CN205670075U (zh) * 2016-04-28 2016-11-02 贵州航空发动机研究所 一种飞机进气道出口流场测量机构

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
李锋 等: "高超声速气动试验模拟现状与面临的新挑战", 《气体物理》 *
杨肖峰 等: "MSL火星探测器高超声速流场预测及气动性分析", 《宇航学报》 *
杨肖峰: "火星探测器气动力热和传热特性研究", 《中国优秀硕士学位论文全文数据库》 *
汪球 等: "高焓激波风洞有效试验时间的测量研究", 《第十六届全国瘠薄与激波管学术会议论文》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107806977A (zh) * 2017-11-29 2018-03-16 中国航空工业集团公司沈阳空气动力研究所 一种组合式宽马赫数高焓脉冲风洞管体结构
CN107806977B (zh) * 2017-11-29 2024-04-09 中国航空工业集团公司沈阳空气动力研究所 一种组合式宽马赫数高焓脉冲风洞管体结构
CN110057537A (zh) * 2019-04-12 2019-07-26 北京空天技术研究所 飞行器气动性能影响预测方法
CN111346528A (zh) * 2020-04-29 2020-06-30 中国空气动力研究与发展中心高速空气动力研究所 一种变比热比冷喷流混合气体的制备方法
CN115828470A (zh) * 2023-01-10 2023-03-21 中国空气动力研究与发展中心计算空气动力研究所 一种喷管型线设计方法、装置、设备及介质
CN115828470B (zh) * 2023-01-10 2023-04-28 中国空气动力研究与发展中心计算空气动力研究所 一种喷管型线设计方法、装置、设备及介质

Also Published As

Publication number Publication date
CN106768818B (zh) 2019-04-30

Similar Documents

Publication Publication Date Title
CN106768818A (zh) 一种激波风洞中混合气体来流运行参数获得方法
Zhou et al. Experimental measurement of wind-sand flux and sand transport for naturally mixed sands
Wallace Space-time correlations in turbulent flow: A review
CN103575631B (zh) 岩石渗透性测试***及测试方法
Schetz et al. Molecular weight and shock-wave effects on transverse injection in supersonic flow
CN103400035B (zh) 一种高可信度快速预测飞行器滚转动导数的方法
Wang et al. A numerical study of the evolution of the solar wind from Ulysses to Voyager 2
CN105550467A (zh) 高速风洞斜孔壁壁面流动模型的建立方法
Theiss et al. Wake flow instability studies behind discrete roughness elements on a generic re-entry capsule
CN114777855A (zh) 基于加速度检测的涡街湿气分相流量测量方法
CN110186522A (zh) 结合涡街幅值特性的湿气过读补偿与流量测量方法
Anyoji et al. Development of a low-density wind tunnel for simulating martian atmospheric flight
Li et al. Multi-component LBM-LES model of the air and methane flow in tunnels and its validation
CN106844820A (zh) 一种基于深空探测试验的二氧化碳来流参数确定新方法
CN102621276B (zh) 可控校正大气水氢氧稳定同位素比率、梯度及其垂直分馏过程测量的装置
Cartieri et al. Using CFD to calculate support interference effects
Nakiboğlu et al. Stack gas dispersion measurements with large scale-PIV, aspiration probes and light scattering techniques and comparison with CFD
CN106338399A (zh) 一种跨超音速总静压探针测量真值的计算方法
Fedioun et al. Boundary layer transition on the LEA hypersonic vehicle forebody
Laderman Adverse pressure gradient effects on supersonic boundary-layer turbulence
CN114993392A (zh) 基于加速度幅值建模的涡街湿气流量测量方法
Surry et al. Model studies of wind effects-a perspective on the problems of experimental technique and instrumentation
Hill Jr et al. Initial Development of Turbulent, Compressible Free Shear Layers
Xiong et al. Experimental and numerical simulation investigations on particle sampling for high-pressure natural gas
Ge et al. The study of in situ wind and gravity wave determination by the first passive falling-sphere experiment in China's northwest region

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant