CN114486155B - 一种高焓激波风洞参数诊断方法和系统 - Google Patents
一种高焓激波风洞参数诊断方法和系统 Download PDFInfo
- Publication number
- CN114486155B CN114486155B CN202111619055.7A CN202111619055A CN114486155B CN 114486155 B CN114486155 B CN 114486155B CN 202111619055 A CN202111619055 A CN 202111619055A CN 114486155 B CN114486155 B CN 114486155B
- Authority
- CN
- China
- Prior art keywords
- shock
- tube
- flow field
- enthalpy
- high enthalpy
- 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.)
- Active
Links
- 230000035939 shock Effects 0.000 title claims abstract description 160
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000003745 diagnosis Methods 0.000 title claims description 16
- 239000007921 spray Substances 0.000 claims abstract description 56
- 238000012360 testing method Methods 0.000 claims abstract description 43
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 34
- 238000000295 emission spectrum Methods 0.000 claims abstract description 33
- 238000005259 measurement Methods 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims abstract description 9
- 239000000523 sample Substances 0.000 claims description 58
- 238000012545 processing Methods 0.000 claims description 30
- 230000003068 static effect Effects 0.000 claims description 28
- 230000003287 optical effect Effects 0.000 claims description 10
- 230000003595 spectral effect Effects 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 8
- 239000013307 optical fiber Substances 0.000 claims description 7
- 238000004847 absorption spectroscopy Methods 0.000 claims description 5
- 238000004993 emission spectroscopy Methods 0.000 claims description 5
- 238000004020 luminiscence type Methods 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 description 3
- 238000010494 dissociation reaction Methods 0.000 description 3
- 230000005593 dissociations Effects 0.000 description 3
- 230000005284 excitation Effects 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
- G01M9/04—Details
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
- G01M9/065—Measuring arrangements specially adapted for aerodynamic testing dealing with flow
- G01M9/067—Measuring arrangements specially adapted for aerodynamic testing dealing with flow visualisation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/08—Aerodynamic models
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种高焓激波风洞流场参数诊断方法和系统,涉及到高焓激波风洞试验领域。本发明的方法利用高焓激波风洞,采用接触测量技术和非接触光谱测量技术,测量激波管末端驻室参数和喷管自由流参数,诊断高焓激波风洞流场。本发明利用接触测量技术、激光纹影技术、非接触吸收光谱测量技术、非接触发射光谱测量技术和多组分多温度数值模拟技术诊断高焓激波风洞激波管末端驻室参数及喷管自由流参数,不仅仅可以获得流场温度、压力、流场组分类别和组分浓度,还可以获得流场非平衡态信息、风洞有效运行时间。
Description
技术领域
本发明涉及一种高焓激波风洞参数诊断方法和系统,属于高超声速空气动力学试验领域。
背景技术
随着人类航天科技不断进步,近地空间可能存在各类超声速及高超声速再入式飞行物体,当再入体高速重返大气层时,再入体与大气之间存在强烈摩擦,其前端形成很强的脱体激波,在激波强烈压缩作用下,大量的飞行器动能转化为热能,形成高焓流动,引起气体分子振动能激发、离解、复合和电离等复杂物理化学变化。
为研究此类再入体再入现象,通常利用高焓激波风洞进行再入物理地面模拟。再入物理地面模拟中高焓流动特点是具有较强的自发发射特性,空气发生离解甚至电离,信息含量极其丰富,是研究流动的平衡与非平衡性质较珍贵的试验资料。
高焓流动中气体分子振动能激发、离解、复合和电离等复杂物理化学变化过程中发生的能级跃迁,导致能量发生变化,此变化会造成分子具有有特定光谱线,因此利用接触测量技术结合非接触光谱技术来诊断高焓流动气体信息,进而获得高焓激波风洞自由流的平动温度、振动温度、气体组分及流动时间等参数,联合多组分多温度数值模拟技术,为高焓激波风洞驻室及自由流参数诊断提供一种经过试验验证的可行思路。本发明即基于上述背景开展。
发明内容
本发明的技术解决问题是:针对高焓激波风洞产生高温、高压试验气体及有效流场时间短的特点,利用接触测量技术、激光纹影技术、非接触吸收光谱技术、非接触发射光谱技术和多组分多问数值模拟技术来诊断高焓流动气体信息,提出了一种高焓激波风洞流场参数诊断方法和系统。提高了高焓激波风洞试验流场精度,能够获得更加准确的自由流参数。
本发明的技术方案是:一种高焓激波风洞流场参数诊断系统,包括:压电传感器、总压传感器、吸收光谱系统、发射光谱系统、皮托压探针、静压探针、驻点热流探针、第一数据处理系统和第二数据处理系统;
在激波管上安装若干压电传感器,在激波管末端安装总压传感器;在激波管末端总压传感器附近安装非接触吸收光谱系统和发射光谱系统的准直器及探测器,在试验段安装非接触吸收光谱系统和发射光谱系统的准直器及探测器;在试验段中的排架上安装皮托压探针、静压探针、驻点热流探针,用于测量试验段流场的皮托压、静压和驻点热流;非接触吸收光谱系统、发射光谱系统和皮托压探针、静压探针、驻点热流探针所测量的数据通过光纤或者数据线分别发送至第一数据处理系统和第二数据处理系统中进行数据处理。
发射光谱系统的测量波段在0.1μm~6μm,根据激波管末端驻室和高焓喷管自由流参数的光谱信息决定。
非接触吸收光谱系统采用NO和O激光器,测量激波管末端驻室和喷管自由流组分NO和O的浓度;利用近红外探测装置,采用一条光谱谱线测量激波管末端驻室和高焓喷管自由流的温度;在试验段的十字排架上安装近红外探测装置,采用两条谱线进行测速,测速的两条光路夹角为30°~60°。
激光纹影系统穿透高焓喷管出口的自发光的流场,并将自发光过滤,拍摄清晰的球头头部脱体激波距离;单色激光光源穿透强光和高温区域,通过在相机前安装配套的滤光片,能够观察到清晰的流场结构。
第一数据采集系统、第二数据采集系统采用PXI测试平台,采用多套PXIe机箱级联方式实现多通道级联,采集频率大于100kHz。
一种高焓激波风洞流场参数诊断方法,包括以下步骤:
步骤(一)、在激波管上安装若干压电传感器,在激波管末端安装总压传感器;在激波管末端总压传感器附近安装非接触吸收光谱系统和发射光谱系统的准直器及探测器,在试验段安装非接触吸收光谱系统和发射光谱系统的准直器及探测器;在试验段中的排架上安装皮托压探针、静压探针、驻点热流探针,用于测量试验段流场的皮托压、静压和驻点热流;将非接触吸收光谱系统、发射光谱系统和皮托压探针、静压探针、驻点热流探针所测量的数据通过光纤或者数据线分别发送至第一数据处理系统和第二数据处理系统中进行数据处理;
步骤(二)、利用激波管上壁面的压电传感器,测量入射激波经过压电传感器的时间间隔Δt,根据压电传感器彼此间的距离ΔL,计算入射激波速度V=ΔL/Δt;
步骤(三)、利用激波管末端壁面安装总压传感器,测量入射激波反射后激波管末端驻室的总压P0;
步骤(四)、利用非接触吸收光谱系统测量激波管末端驻室气体和高焓喷管(3)自由流气体的平动温度、NO/O组分浓度;
步骤(五)、利用发射光谱系统测量激波管末端驻室气体和喷管自由流气体的振动温度、组分光谱信息;
步骤(六)、在试验段十字排架上安装球头,利用激光纹影系统,测量球头周围流场的结构,获得球头头部激波的距离;
步骤(七)、利用试验段十字排架上安装皮托压探针、静压探针和驻点热流探针,获得流场的皮托压、静压和驻点热流;
步骤(八)、利用步骤(二)和步骤(三)测量获得入射激波速度V和激波管末端驻室的总压P0,采用高温空气的热力学数值模型和准一维激波管理论,计算高温状态条件(2000K~10000K)下的激波管末端驻室参数,将激波管末端驻室参数与步骤(四)中测得的激波管末端驻室的温度和NO/O组分浓度进行对比,若两者数据的偏差大于设定的阈值a,则需要修改热力学数值模型,再次迭代,直到两者数据的偏差≤设定的阈值a;
步骤(九)、将步骤(八)中迭代完成的激波管末端驻室参数作为输入数值,计算高焓喷管流场的初始条件,高焓喷管流场计算利用多组分多温度的热化学非平衡数值模拟方法进行,计算得到的高焓喷管出口参数;将高焓喷管出口参数与步骤(三)~步骤(七)中的测量数据进行对比,若两者数据的偏差大于设定的阈值b,则修改数值模型,再次迭代,直到两者数据的偏差≤设定的阈值b。
所述的一种高焓激波风洞流场参数诊断方法,利用发射光谱系统测量高焓喷管驻室和喷管出口自由流气体组分,根据是否有驱动气体组分,判断风洞有效运行时间t1;
利用非接触吸收光谱系统测量高焓喷管驻室和喷管出口自由流的NO/O浓度大小,判断风洞有效运行时间t2;
利用激光纹影系统拍摄的球头脱体激波距离,根据脱体激波距离,判断风洞有效运行时间t3;取时间t1、时间t2和时间t3中的最小值作为风洞有效运行时间。
激光纹影系统中采用的单色激光光源的波长范围为520~720nm。
本发明与现有技术相比的有益效果是:
本发明利用接触测量技术、激光纹影技术、非接触吸收光谱测量技术、非接触发射光谱测量技术和多组分多温度数值模拟技术等共同诊断高焓激波风洞激波管末端驻室参数及喷管自由流参数,不仅仅可以获得流场温度、压力、流场组分类别和组分浓度,还可以获得流场非平衡态信息、风洞有效运行时间。利用本发明的高焓激波风洞流场参数诊断方法,可以搭建高焓流场参数测量平台,具有非常高的灵敏度,连续的时间分辨率及快速时间响应,是一种较好的优选方法。
附图说明
图1为本发明涉及的高焓激波风洞流场参数诊断方法流程。
图2为本发明涉及的高焓激波风洞流场参数诊断方法示意图。
其中,1.激波管,2.膜片,3.高焓喷管,4.试验段,5.安装传感器的十字排架,6.光学窗口,7.吸收光谱系统,8.数据处理系统,9.发射光谱系统,10.数据处理系统。
具体实施方式
下面通过对本发明进行详细说明,本发明的特点和优点将随着这些说明而变得更为清楚、明确。
在这里专用的词“示例性”意为“用作例子、实施例或说明性”。这里作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。尽管在附图中示出了实施例的各种方面,但是除非特别指出,不必按比例绘制附图。
如图2,一种高焓激波风洞流场参数诊断系统,包括:压电传感器、总压传感器、吸收光谱系统7、发射光谱系统9、皮托压探针、静压探针、驻点热流探针、第一数据处理系统8和第二数据处理系统10;
在激波管1上安装若干压电传感器,在激波管1末端安装总压传感器;在激波管1末端总压传感器附近安装非接触吸收光谱系统7和发射光谱系统9的准直器及探测器,在试验段4安装非接触吸收光谱系统7(吸收光谱系统7的光纤探头设置在试验段4的光学窗口6处)和发射光谱系统9的准直器及探测器;在试验段4中的排架5上安装皮托压探针、静压探针、驻点热流探针,用于测量试验段流场的皮托压、静压和驻点热流;非接触吸收光谱系统7、发射光谱系统9和皮托压探针、静压探针、驻点热流探针所测量的数据通过光纤或者数据线分别发送至第一数据处理系统8和第二数据处理系统10中进行数据处理。
发射光谱系统9的测量波段在0.1μm~6μm,根据激波管1末端驻室和高焓喷管3自由流参数的光谱信息决定。
非接触吸收光谱系统7采用NO和O激光器,测量激波管1末端驻室和喷管自由流组分NO和O的浓度;利用近红外探测装置,采用一条光谱谱线测量激波管1末端驻室和高焓喷管3自由流的温度;在试验段的十字排架5上安装近红外探测装置,采用两条谱线进行测速,测速的两条光路夹角为30°~60°。
激光纹影系统穿透高焓喷管3出口的自发光的流场,并将自发光过滤,拍摄清晰的球头头部脱体激波距离;单色激光光源穿透强光和高温区域,通过在相机前安装配套的滤光片,能够观察到清晰的流场结构。
第一数据采集系统8、第二数据采集系统10采用PXI测试平台,采用多套PXIe机箱级联方式实现多通道级联,采集频率大于100kHz。
如图1,本发明给出一种高焓激波风洞流场参数诊断方法,为高焓流场参数诊断提供一种经过试验验证的可行思路,关键点利用利用接触测量技术、激光纹影技术、非接触吸收光谱测量技术、非接触发射光谱测量技术和多组分多温度数值模拟技术等共同诊断高焓激波风洞激波管1末端驻室参数及喷管自由流参数,当验证结果满足流场要求,即试验测量结果与数值模拟结果相匹配时,完成高焓激波风洞流场参数诊断。本发明不仅仅可以获得流场温度、压力、流场组分类别和组分浓度,还可以获得流场非平衡态信息、风洞有效运行时间。
步骤一、在激波管1上安装一系列压电传感器,在激波管1末端安装总压传感器。在激波管1末端总压传感器附近安装非接触吸收光谱系统7和发射光谱系统9的准直器及探测器,同时在试验段4也安装非接触吸收光谱系统7和发射光谱系统9的准直器及探测器。在试验段4中的十字排架5上安装皮托压探针、静压探针、驻点热流探针,将非接触吸收光谱系统7、发射光谱系统9和皮托压探针、静压探针、驻点热流探针所测量的数据通过光纤或者数据线分别发送至第一数据处理系统8和第二数据处理系统10中进行数据处理;激波管1和高焓喷管3之间设置膜片2。
步骤二、利用激波管1上壁面的压电传感器,测量入射激波经过压电传感器的时间间隔Δt,根据压电传感器彼此间的距离ΔL,计算入射激波速度V=ΔL/Δt。
步骤三、利用激波管1末端壁面安装总压传感器,测量入射激波反射后激波管1末端驻室的总压P0,测量获得入射激波速度V和激波管1末端驻室的总压P0,采用高温空气的热力学数据,计算高温状态条件下的驻室参数。
步骤四、利用非接触吸收光谱系统7测量激波管1末端驻室气体和高焓喷管3自由流气体的平动温度、NO/O组分浓度等。
步骤五、在激波管1末端靠近总压传感器位置以及试验段都安装发射光谱系统9,测量激波管1末端驻室气体和喷管3自由流气体的振动温度、组分光谱信息。
步骤六、在试验段4流场均匀区安装球头,利用激光纹影系统,测量球头周围流场的结构,获得球头头部激波的距离。
步骤七、在试验段4流场均匀区安装皮托压探针、静压探针和驻点热流探针,获得流场的皮托压、静压和驻点热流。
步骤八、利用步骤二和步骤三测量获得入射激波速度V和激波管1末端驻室的总压P0,采用高温空气的热力学数值模型和准一维激波管理论,计算高温状态条件(2000K~10000K)下的激波管1末端驻室参数,将激波管1末端驻室参数与测量激波管1末端驻室的温度和NO/O组分浓度进行对比,若两者数据的偏差大于4%,则需要修改热力学数值模型,再次迭代,直到两者数据的偏差≤4%;
步骤九、将步骤八中迭代完成的激波管1末端驻室参数作为输入数值,计算高焓喷管3流场的初始条件,高焓喷管3流场计算利用多组分多温度的热化学非平衡数值模拟方法进行,计算得到的高焓喷管3出口参数;将高焓喷管3出口参数与步骤三~步骤七中的测量数据进行对比,若两者数据的偏差大于6%,则修改数值模型,再次迭代,直到两者数据的偏差≤6%。
利用发射光谱系统9测量高焓喷管3驻室和喷管出口自由流气体组分,根据是否有驱动气体组分,判断风洞有效运行时间t1;
利用非接触吸收光谱系统7测量高焓喷管3驻室和喷管出口自由流的NO/O浓度大小,判断风洞有效运行时间t2;
利用激光纹影系统拍摄的球头脱体激波距离,根据脱体激波距离,判断风洞有效运行时间t3;取时间t1、时间t2和时间t3中的最小值作为风洞有效运行时间。
本发明利用接触测量技术、激光纹影技术、非接触吸收光谱测量技术、非接触发射光谱测量技术和多组分多温度数值模拟技术等共同诊断高焓激波风洞激波管1末端驻室参数及喷管自由流参数,不仅仅可以获得流场温度、压力、流场组分类别和组分浓度,还可以获得流场非平衡态信息、风洞有效运行时间。利用本发明的高焓激波风洞流场参数诊断方法,可以搭建高焓流场参数测量平台,具有非常高的灵敏度,连续的时间分辨率及快速时间响应,是一种较好的优选方法。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (13)
1.一种高焓激波风洞流场参数诊断方法,其特征在于,包括以下步骤:
步骤(一)、在激波管(1)上安装若干压电传感器,在激波管(1)末端安装总压传感器;在激波管(1)末端总压传感器附近安装非接触吸收光谱系统(7)和发射光谱系统(9)的准直器及探测器,在试验段(4)安装非接触吸收光谱系统(7)和发射光谱系统(9)的准直器及探测器;在试验段(4)中的排架(5)上安装皮托压探针、静压探针、驻点热流探针,用于测量试验段流场的皮托压、静压和驻点热流;将非接触吸收光谱系统(7)、发射光谱系统(9)和皮托压探针、静压探针、驻点热流探针所测量的数据通过光纤或者数据线分别发送至第一数据处理系统(8)和第二数据处理系统(10)中进行数据处理;
步骤(二)、利用激波管(1)上壁面的压电传感器,测量入射激波经过压电传感器的时间间隔Δt,根据压电传感器彼此间的距离ΔL,计算入射激波速度V=ΔL/Δt;
步骤(三)、利用激波管(1)末端壁面安装总压传感器,测量入射激波反射后激波管(1)末端驻室的总压P0;
步骤(四)、利用非接触吸收光谱系统(7)测量激波管(1)末端驻室气体和高焓喷管(3)自由流气体的平动温度、NO/O组分浓度;
步骤(五)、利用发射光谱系统(9)测量激波管(1)末端驻室气体和喷管(3)自由流气体的振动温度、组分光谱信息;
步骤(六)、在试验段十字排架(5)上安装球头,利用激光纹影系统,测量球头周围流场的结构,获得球头头部激波的距离;
步骤(七)、利用试验段十字排架(5)上安装皮托压探针、静压探针和驻点热流探针,获得流场的皮托压、静压和驻点热流;
步骤(八)、利用步骤(二)和步骤(三)测量获得入射激波速度V和激波管(1)末端驻室的总压P0,采用高温空气的热力学数值模型和准一维激波管理论,计算高温状态条件下的激波管(1)末端驻室参数,将激波管(1)末端驻室参数与步骤(四)中测得的激波管(1)末端驻室的温度和NO/O组分浓度进行对比,若两者数据的偏差大于设定的阈值a,则需要修改热力学数值模型,再次迭代,直到两者数据的偏差≤设定的阈值a;
步骤(九)、将步骤(八)中迭代完成的激波管(1)末端驻室参数作为输入数值,计算高焓喷管(3)流场的初始条件,高焓喷管(3)流场计算利用多组分多温度的热化学非平衡数值模拟方法进行,计算得到的高焓喷管(3)出口参数;将高焓喷管(3)出口参数与步骤(三)~步骤(七)中的测量数据进行对比,若两者数据的偏差大于设定的阈值b,则修改数值模型,再次迭代,直到两者数据的偏差≤设定的阈值b。
2.根据权利要求1所述的一种高焓激波风洞流场参数诊断方法,其特征在于:发射光谱系统(9)的测量波段在0.1μm~6μm,根据激波管(1)末端驻室和高焓喷管(3)自由流参数的光谱信息决定。
3.根据权利要求2所述的一种高焓激波风洞流场参数诊断方法,其特征在于:非接触吸收光谱系统(7)采用NO和O激光器,测量激波管(1)末端驻室和喷管自由流组分NO和O的浓度;利用近红外探测装置,采用一条光谱谱线测量激波管(1)末端驻室和高焓喷管(3)自由流的温度;在试验段的十字排架(5)上安装近红外探测装置,采用两条谱线进行测速,测速的两条光路夹角为30°~60°。
4.根据权利要求3所述的一种高焓激波风洞流场参数诊断方法,其特征在于:激光纹影系统穿透高焓喷管(3)出口的自发光的流场,并将自发光过滤,拍摄清晰的球头头部脱体激波距离;单色激光光源穿透强光和高温区域,通过在相机前安装配套的滤光片,能够观察到清晰的流场结构。
5.根据权利要求4所述的一种高焓激波风洞流场参数诊断方法,其特征在于:激光纹影系统中采用的单色激光光源的波长范围为520~720nm。
6.根据权利要求5所述的一种高焓激波风洞流场参数诊断方法,其特征在于:第一数据处理系统(8)、第二数据处理系统(10)采用PXI测试平台,采用多套PXIe机箱级联方式实现多通道级联,采集频率大于100kHz。
7.根据权利要求1所述的一种高焓激波风洞流场参数诊断方法,其特征在于:利用发射光谱系统(9)测量高焓喷管(3)驻室和喷管出口自由流气体组分,根据是否有驱动气体组分,判断风洞有效运行时间t1;
利用非接触吸收光谱系统(7)测量高焓喷管(3)驻室和喷管出口自由流的NO/O浓度大小,判断风洞有效运行时间t2;
利用激光纹影系统拍摄的球头脱体激波距离,根据脱体激波距离,判断风洞有效运行时间t3;取时间t1、时间t2和时间t3中的最小值作为风洞有效运行时间。
8.一种高焓激波风洞流场参数诊断系统,其特征在于,包括:压电传感器、总压传感器、吸收光谱系统(7)、发射光谱系统(9)、皮托压探针、静压探针、驻点热流探针、第一数据处理系统(8)和第二数据处理系统(10);
在激波管(1)上安装若干压电传感器,在激波管(1)末端安装总压传感器;在激波管(1)末端总压传感器附近安装非接触吸收光谱系统(7)和发射光谱系统(9)的准直器及探测器,在试验段(4)安装非接触吸收光谱系统(7)和发射光谱系统(9)的准直器及探测器;在试验段(4)中的排架(5)上安装皮托压探针、静压探针、驻点热流探针,用于测量试验段流场的皮托压、静压和驻点热流;非接触吸收光谱系统(7)、发射光谱系统(9)和皮托压探针、静压探针、驻点热流探针所测量的数据通过光纤或者数据线分别发送至第一数据处理系统(8)和第二数据处理系统(10)中进行数据处理;
根据所述高焓激波风洞流场参数诊断系统进行高焓激波风洞流场参数诊断方法如下:
步骤一、利用激波管(1)上壁面的压电传感器,测量入射激波经过压电传感器的时间间隔Δt,根据压电传感器彼此间的距离ΔL,计算入射激波速度V=ΔL/Δt;
步骤二、利用激波管(1)末端壁面安装总压传感器,测量入射激波反射后激波管(1)末端驻室的总压P0;
步骤三、利用非接触吸收光谱系统(7)测量激波管(1)末端驻室气体和高焓喷管(3)自由流气体的平动温度、NO/O组分浓度;
步骤四、利用发射光谱系统(9)测量激波管(1)末端驻室气体和喷管(3)自由流气体的振动温度、组分光谱信息;
步骤五、在试验段十字排架(5)上安装球头,利用激光纹影系统,测量球头周围流场的结构,获得球头头部激波的距离;
步骤六、利用试验段十字排架(5)上安装皮托压探针、静压探针和驻点热流探针,获得流场的皮托压、静压和驻点热流;
步骤七、利用步骤一和步骤二测量获得入射激波速度V和激波管(1)末端驻室的总压P0,采用高温空气的热力学数值模型和准一维激波管理论,计算高温状态条件下的激波管(1)末端驻室参数,将激波管(1)末端驻室参数与步骤三中测得的激波管(1)末端驻室的温度和NO/O组分浓度进行对比,若两者数据的偏差大于设定的阈值a,则需要修改热力学数值模型,再次迭代,直到两者数据的偏差≤设定的阈值a;
步骤八、将步骤七中迭代完成的激波管(1)末端驻室参数作为输入数值,计算高焓喷管(3)流场的初始条件,高焓喷管(3)流场计算利用多组分多温度的热化学非平衡数值模拟方法进行,计算得到的高焓喷管(3)出口参数;将高焓喷管(3)出口参数与步骤二~步骤六中的测量数据进行对比,若两者数据的偏差大于设定的阈值b,则修改数值模型,再次迭代,直到两者数据的偏差≤设定的阈值b。
9.根据权利要求8所述的一种高焓激波风洞流场参数诊断系统,其特征在于:发射光谱系统(9)的测量波段在0.1μm~6μm,根据激波管(1)末端驻室和高焓喷管(3)自由流参数的光谱信息决定。
10.根据权利要求9所述的一种高焓激波风洞流场参数诊断系统,其特征在于:非接触吸收光谱系统(7)采用NO和O激光器,测量激波管(1)末端驻室和喷管自由流组分NO和O的浓度;利用近红外探测装置,采用一条光谱谱线测量激波管(1)末端驻室和高焓喷管(3)自由流的温度;在试验段的十字排架(5)上安装近红外探测装置,采用两条谱线进行测速,测速的两条光路夹角为30°~60°。
11.根据权利要求10所述的一种高焓激波风洞流场参数诊断系统,其特征在于:激光纹影系统穿透高焓喷管(3)出口的自发光的流场,并将自发光过滤,拍摄清晰的球头头部脱体激波距离;单色激光光源穿透强光和高温区域,通过在相机前安装配套的滤光片,能够观察到清晰的流场结构。
12.根据权利要求11所述的一种高焓激波风洞流场参数诊断系统,其特征在于:激光纹影系统中采用的单色激光光源的波长范围为520~720nm。
13.根据权利要求12所述的一种高焓激波风洞流场参数诊断系统,其特征在于:第一数据处理系统(8)、第二数据处理系统(10)采用PXI测试平台,采用多套PXIe机箱级联方式实现多通道级联,采集频率大于100kHz。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111619055.7A CN114486155B (zh) | 2021-12-27 | 2021-12-27 | 一种高焓激波风洞参数诊断方法和系统 |
US18/023,679 US20240310240A1 (en) | 2021-12-27 | 2021-12-30 | Method and system for diagnosing high-enthalpy shock tunnel parameters |
PCT/CN2021/142954 WO2023123180A1 (zh) | 2021-12-27 | 2021-12-30 | 一种高焓激波风洞参数诊断方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111619055.7A CN114486155B (zh) | 2021-12-27 | 2021-12-27 | 一种高焓激波风洞参数诊断方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114486155A CN114486155A (zh) | 2022-05-13 |
CN114486155B true CN114486155B (zh) | 2022-10-14 |
Family
ID=81496526
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111619055.7A Active CN114486155B (zh) | 2021-12-27 | 2021-12-27 | 一种高焓激波风洞参数诊断方法和系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20240310240A1 (zh) |
CN (1) | CN114486155B (zh) |
WO (1) | WO2023123180A1 (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115343016B (zh) * | 2022-07-18 | 2024-08-06 | 中国航天空气动力技术研究院 | 一种高焓激波风洞流场动压的测量方法及系统 |
CN115618168B (zh) * | 2022-11-30 | 2023-03-14 | 中国空气动力研究与发展中心空天技术研究所 | 高马赫数高焓风洞中皮托压数据处理方法 |
CN116296227B (zh) * | 2023-05-26 | 2023-07-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法 |
CN118376379B (zh) * | 2024-06-21 | 2024-08-16 | 中国空气动力研究与发展中心超高速空气动力研究所 | 高焓瞬时流场光辐射多角度多波段的同时测量装置及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107976295A (zh) * | 2017-12-27 | 2018-05-01 | 中国航天空气动力技术研究院 | 一种2m量级自由活塞驱动的高焓激波风洞 |
CN207703439U (zh) * | 2017-12-27 | 2018-08-07 | 中国航天空气动力技术研究院 | 一种2m量级自由活塞驱动的高焓激波风洞 |
CN109655227A (zh) * | 2018-12-07 | 2019-04-19 | 中国航天空气动力技术研究院 | 一种低焓电弧加热器气流焓值诊断系统及诊断方法 |
CN112649172A (zh) * | 2020-12-21 | 2021-04-13 | 中国航天空气动力技术研究院 | 静压探针及高焓激波风洞静压测量方法 |
CN113280996A (zh) * | 2021-04-25 | 2021-08-20 | 中国航天空气动力技术研究院 | 一种高焓流场自由流的速度测量方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06258177A (ja) * | 1993-03-08 | 1994-09-16 | Mitsubishi Heavy Ind Ltd | 高エンタルピ風洞 |
JP3409083B2 (ja) * | 2001-05-11 | 2003-05-19 | 独立行政法人航空宇宙技術研究所 | 自由噴流式極超音速風洞試験装置 |
US8833153B2 (en) * | 2012-09-20 | 2014-09-16 | The Boeing Company | Correction of pressure signals measured during supersonic wind tunnel testing |
CN102937655B (zh) * | 2012-10-31 | 2015-06-03 | 中国科学院力学研究所 | 一种测量激波速度的系统和方法 |
CN107806977B (zh) * | 2017-11-29 | 2024-04-09 | 中国航空工业集团公司沈阳空气动力研究所 | 一种组合式宽马赫数高焓脉冲风洞管体结构 |
CN111024357B (zh) * | 2019-12-11 | 2021-07-13 | 中国航天空气动力技术研究院 | 一种大尺寸自由活塞高焓激波风洞模拟飞行环境的方法 |
CN112067241B (zh) * | 2020-08-19 | 2022-08-12 | 中国航天空气动力技术研究院 | 一种高焓激波风洞总温测量方法 |
CN112067243A (zh) * | 2020-08-25 | 2020-12-11 | 中国航天空气动力技术研究院 | 一种用于高焓激波风洞的流场温度测量方法 |
-
2021
- 2021-12-27 CN CN202111619055.7A patent/CN114486155B/zh active Active
- 2021-12-30 US US18/023,679 patent/US20240310240A1/en active Pending
- 2021-12-30 WO PCT/CN2021/142954 patent/WO2023123180A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107976295A (zh) * | 2017-12-27 | 2018-05-01 | 中国航天空气动力技术研究院 | 一种2m量级自由活塞驱动的高焓激波风洞 |
CN207703439U (zh) * | 2017-12-27 | 2018-08-07 | 中国航天空气动力技术研究院 | 一种2m量级自由活塞驱动的高焓激波风洞 |
CN109655227A (zh) * | 2018-12-07 | 2019-04-19 | 中国航天空气动力技术研究院 | 一种低焓电弧加热器气流焓值诊断系统及诊断方法 |
CN112649172A (zh) * | 2020-12-21 | 2021-04-13 | 中国航天空气动力技术研究院 | 静压探针及高焓激波风洞静压测量方法 |
CN113280996A (zh) * | 2021-04-25 | 2021-08-20 | 中国航天空气动力技术研究院 | 一种高焓流场自由流的速度测量方法 |
Non-Patent Citations (2)
Title |
---|
Diagnosis of iso-octane combustion in a shock tube by emission spectroscopy;Changhua Zhang等;《Chemical Physics letters》;20130131;第13-17页 * |
高焓激波风洞试验技术综述;谌君谋等;《空气动力学学报》;20180815;第36卷(第04期);第543-552页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2023123180A1 (zh) | 2023-07-06 |
CN114486155A (zh) | 2022-05-13 |
US20240310240A1 (en) | 2024-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114486155B (zh) | 一种高焓激波风洞参数诊断方法和系统 | |
CN107884392B (zh) | 电弧加热器气流焓值光谱测量系统 | |
Marineau et al. | Investigation of Mach 10 boundary layer stability of sharp cones at angle-of-attack, part 1: experiments | |
Watkins et al. | Measuring surface pressures on rotor blades using pressure-sensitive paint | |
JP2010117360A (ja) | 補正した空気流量の直接的非貫入型測定のためのシステム及び方法 | |
Kennedy et al. | Visualizations of hypersonic boundary-layer transition on a variable bluntness cone | |
Burns et al. | Application of FLEET velocimetry in the NASA Langley 0.3-meter transonic cryogenic tunnel | |
Bounitch et al. | Improved measurements of" Tunnel Noise" pressure fluctuations in the AEDC hypervelocity wind tunnel No. 9 | |
Dhanagopal et al. | High-Speed Pressure Sensitive Paint Measurements of the Initial Concept 3. X Vehicle at Mach 7 | |
Holden | Experimental studies in the LENS shock tunnel and expansion tunnel to examine real-gas effects in hypervelocity flows | |
MacDonald et al. | Temperature and CO number density measurements in shocked CO and CO2 via tunable diode laser absorption spectroscopy | |
Allen | Nonequilibrium shock front rotational, vibrational and electronic temperature measurements | |
Eitelberg | First results of calibration and use of the HEG | |
Hoberg et al. | Freestream Characterization and Condensation Detection in the AFOSR--Notre Dame Large Mach-6 Quiet Tunnel | |
CN106844820B (zh) | 一种基于深空探测试验的二氧化碳来流参数确定方法 | |
CN114184348B (zh) | 一种高焓流场光电特性辨识装置及方法 | |
Xu et al. | A tunable diode laser absorption based velocity sensor for local field in hypersonic flows | |
James et al. | Detecting shock arrival in expansion tubes and shock tunnels using high-frequency photodiodes | |
Schneider et al. | Progress in the operation of the Boeing/AFOSR Mach-6 quiet tunnel | |
Dhanagopal et al. | Global Surface Temperature Distribution on the Initial Concept 3. X Vehicle at Mach 7 | |
Wadhams et al. | Recent experimental studies of high speed boundary layer transition in LENS facilities to further the development of predictive tools for boundary layer transition in flight | |
Reda et al. | Aerothermodynamic testing of ablative reentry vehicle nosetip materials in hypersonic ballistic-range environments | |
Tanno et al. | Measurement of hypersonic high-enthalpy boundary layer transition on a 7º cone model | |
Wadhams et al. | Continuing experimental studies of high speed boundary layer transition in lens facilities to further the development of predictive tools for boundary layer transition in flight | |
Sharma | Post-shock thermochemistry in hypervelocity CO2 and air flow |
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 |