CN116296227A - 一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法 - Google Patents
一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法 Download PDFInfo
- Publication number
- CN116296227A CN116296227A CN202310604823.4A CN202310604823A CN116296227A CN 116296227 A CN116296227 A CN 116296227A CN 202310604823 A CN202310604823 A CN 202310604823A CN 116296227 A CN116296227 A CN 116296227A
- Authority
- CN
- China
- Prior art keywords
- heat flow
- test
- radiation
- wind tunnel
- model
- 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
Links
- 238000010998 test method Methods 0.000 title claims abstract description 28
- 230000005855 radiation Effects 0.000 claims abstract description 152
- 238000012360 testing method Methods 0.000 claims abstract description 120
- 230000001052 transient effect Effects 0.000 claims abstract description 73
- 238000010438 heat treatment Methods 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000004519 manufacturing process Methods 0.000 claims abstract description 5
- 239000005304 optical glass Substances 0.000 claims description 30
- 238000000034 method Methods 0.000 claims description 27
- 230000035939 shock Effects 0.000 claims description 25
- 239000000463 material Substances 0.000 claims description 18
- 238000010521 absorption reaction Methods 0.000 claims description 11
- 239000010408 film Substances 0.000 claims description 9
- 238000009434 installation Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 7
- 239000000758 substrate Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000001595 flow curve Methods 0.000 claims description 5
- 238000002834 transmittance Methods 0.000 claims description 5
- 230000005457 Black-body radiation Effects 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 4
- 229910000838 Al alloy Inorganic materials 0.000 claims description 3
- 239000000853 adhesive Substances 0.000 claims description 3
- 230000001070 adhesive effect Effects 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 3
- 239000007921 spray Substances 0.000 claims description 3
- 239000010935 stainless steel Substances 0.000 claims description 3
- 229910001220 stainless steel Inorganic materials 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 239000010409 thin film Substances 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 9
- 238000011160 research Methods 0.000 abstract description 4
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 6
- 229910052697 platinum Inorganic materials 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- DYRBFMPPJATHRF-UHFFFAOYSA-N chromium silicon Chemical compound [Si].[Cr] DYRBFMPPJATHRF-UHFFFAOYSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000010963 304 stainless steel Substances 0.000 description 1
- 229910000589 SAE 304 stainless steel Inorganic materials 0.000 description 1
- 230000037237 body shape Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012546 transfer 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
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K17/00—Measuring quantity of heat
- G01K17/06—Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device
- G01K17/08—Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device based upon measurement of temperature difference or of a temperature
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明属于超高速风洞试验技术领域，公开了一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法。试验方法利用高焓膨胀风洞产生的极短时超高速流场，模拟超高速飞行器高温绕流，通过隔离气流对试验模型的对流加热，测量高温绕流对试验模型表面形成的辐射热流。试验方法首先根据飞行工况确定高焓膨胀风洞的模拟试验条件，再根据模拟试验条件估算辐射主导波段及辐射热流范围后设计制作瞬态辐射热流传感器，然后对瞬态辐射热流传感器的测量信号进行数据处理，得到试验模型表面测点处的辐射热流。本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法为高焓气动热环境研究提供了技术支持。
Description
技术领域
本发明属于超高速风洞试验技术领域，具体涉及一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法。
背景技术
超高速飞行器飞行过程中，头部的激波层气体发生离解甚至电离，受激发的气体粒子会以辐射的形式对超高速飞行器表面加热。在飞行速度9km/s时，辐射加热能够达到总气动加热的30%左右；在更高速度下，辐射加热甚至在总气动加热中占主导地位。因此，辐射热流的预测是超高速飞行过程的关键问题之一。辐射加热的计算预测具有较大的难度和不确定度，需要开展相应的地面模拟试验研究。
高焓膨胀风洞的高温气流速度高达12km/s，具有超高速气流模拟能力，模型头部的激波层温度最高超过10000K，具备了飞行器辐射加热效应的模拟能力，但是，目前还缺乏相应的辐射热流测量技术。
当前，亟需发展一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法。
发明内容
本发明所要解决的技术问题是提供一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法。
本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，包括以下步骤：
S1.确定高焓膨胀风洞试验条件；
S11.根据高焓膨胀风洞试验需要模拟的飞行高度和速度，确定高焓膨胀风洞的运行状态，并根据高焓膨胀风洞流场调试校测获得的流场参数，确定试验模型尺寸；
S12.采用CFD计算试验模型头部的激波层参数，估算试验模型表面的辐射热流范围，设置瞬态辐射热流传感器的位置，即测点位置；
根据激波层气流温度，假设测点的光辐射亮度
符合黑体辐射理论，按下式预估绕流的主要辐射波段：
式中，
为黑体温度，/>
为辐射波长，/>
为第一辐射常数，/>
为第二辐射常数；
计算辐射与波长的对应曲线，确定试验工况下，激波层气流对试验模型表面的辐射主导波段；
试验模型表面的辐射热流q r按下式估算：
其中，
，为常数；r n为模型头部等效球头半径，单位为m；/>
为气流密度，单位为kg/m3；V为气流速度，单位为m/s；a、b及f(V)按下式计算：
辐射热流q r的单位为W/cm2；a max表示a的最大值；
S13.根据高温绕流流场温度所对应的主要辐射波段及模型表面辐射热流，确定瞬态辐射热流传感器表面的柱形光学玻璃及内部的对流热流计的参数及选型；
S2.安装试验装置；
S21.制作瞬态辐射热流传感器；瞬态辐射热流传感器为从外至内同中心轴线的双层柱体结构，包括位于外层的安装套、位于内层的瞬态辐射热流传感器测量体，瞬态辐射热流传感器测量体包括从上至下顺序固定的柱形光学玻璃和对流热流计；对流热流计的上表面设置辐射吸收层；对流热流计的下表面设置信号线，信号线穿出瞬态辐射热流传感器测量体，与数据采集设备连接；
S22.在高焓膨胀风洞的试验段舱体内，通过模型支杆将试验模型安装在攻角机构上，试验模型位于风洞喷管的出口、试验段的流场均匀区内，在试验模型的头部驻点安装瞬态辐射热流传感器，瞬态辐射热流传感器与试验模型的表面平齐或者平滑过渡；瞬态辐射热流传感器将高温流场的辐射光信号转变成电压信号，电压信号经数据信号线传输至数据采集设备，通过控制及数据处理计算机的数据处理软件处理电压信号，得到试验模型表面的辐射热流q r；
S23.在高焓膨胀风洞的洞体外，安装数据采集设备的触发器，触发器以高焓膨胀风洞的激波管管壁的压电传感器发送的电压信号为触发信号；
S24.启动数据采集设备和触发器，使得试验装置处于等待状态；
S3.进行高焓膨胀风洞试验；
启动高焓膨胀风洞，高焓膨胀风洞上游激波管中形成的入射激波作用于压电传感器，压电传感器输出电压信号触发数据采集设备，数据采集设备开始采集，高焓膨胀风洞试验结束后，停止采集；
S4.进行试验数据处理；
S41.计算对流热流计表面热流曲线：
瞬态辐射热流传感器输出的信号为对流热流计感受到加热后输出的电压变化信号，将电压信号通过对流热流计的电阻温度系数转换为温度后，采用基于一维半无限假设的热流计算方法，获得热流随时间变化曲线：
式中，
表示瞬态辐射热流传感器的对流热流计直接测量的对流热流计表面热流；/>
为对流热流计基底材料热物性综合参数；其中，/>
为基底材料的密度，单位为kg/m3；c为基底材料的比热，单位为J/（kg·K），k为基底材料的热传导系数，单位为w/（m·K）；T为温度；t为时间；
将热流随时间的变化曲线，结合高焓膨胀风洞流场的有效试验时间区间，以试验气流到达前的数据平均值作为基线值，有效试验时间的热流平均值作为测量平均值，测量平均值与基线值的之差为有效热流平均值；
S42.将对流热流计表面热流换算为试验模型表面辐射热流；
根据瞬态辐射热流传感器的结构和辐射光的传输过程，将对流热流计表面热流换算为试验模型表面辐射热流：
式中，
为试验模型表面辐射热流，单位为W/m2；/>
为对流热流计表面热流，单位为W/m2；/>
为对流热流计的敏感元件面积，单位为m2；/>
为瞬态辐射热流传感器的柱形光学玻璃的窗口面积，单位为m2；/>
为瞬态辐射热流传感器的视角因子；/>
为对流热流计前端面的辐射吸收层的吸收率；/>
为瞬态辐射热流传感器的柱形光学玻璃的透射率；
其中，视角因子按下式计算：
式中，
为对流热流计的敏感元件能够看到的柱形光学玻璃的窗口半径，单位为m；/>
为对流热流计的敏感元件到柱形光学玻璃的窗口外表面的高度，单位为m。
进一步地，所述的瞬态辐射热流传感器安装到试验模型壳体上预留的安装通孔内，采用胶合剂从试验模型内壁面一侧粘结固定。
进一步地，所述的柱形光学玻璃采用MgF2玻璃。
进一步地，所述的对流热流计为薄膜热流计或同轴热电偶。
进一步地，所述的辐射吸收层的感应波长范围200nm～2500nm。
进一步地，所述的安装套的材质与试验模型的壳体的材质相同；或者安装套的材质采用铝合金或者不锈钢中的一种；安装套内壁涂敷有反射膜。
由于高焓膨胀风洞的有效试验气流持续时间极短，通常在几百微秒以内，因此，瞬态辐射热流传感器必须能够快速响应；同时，瞬态辐射热流传感器还必须能够隔离高温气流的对流加热，并考虑辐射场的光谱特性及辐射吸收率等影响因素。
本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法使用的瞬态辐射热流传感器根据模拟试验条件预估试验模型的绕流温度、主导辐射波段和辐射热流值，确定瞬态辐射热流传感器的柱形光学玻璃及对流热流计的参数和选型，标定柱形光学玻璃的透过率，辐射吸收层的吸收率，以及对流热流计的灵敏度系数。在对流热流计的上方采用柱形光学玻璃作为辐射透射窗口，让辐射光通过，柱形光学玻璃还能够隔绝气流以隔离对流加热，采用快速响应的对流热流计作为感热器件，在对流热流计的上表面即感应面增加辐射吸收层以吸收光辐射，并设计相应的封装结构，最终获得响应快、结构紧凑、安装方便的瞬态辐射热流传感器。
本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法根据瞬态辐射热流传感器的结构和辐射光传输特性，建立了辐射热流数据处理方法。对流热流计的输出信号是在试验模型测点当地，进入到瞬态辐射热流传感器内部的对流热流计上表面，即感应面的辐射加热引起的热流计表面温升-电压信号；首先将此温升-电压信号经一维半无限假设传热理论处理后得到对流热流计表面热流，然后考虑试验模型外形和瞬态辐射热流传感器构型，根据辐射传输过程，再通过对流热流计热流计算与辐射传输换算相结合，最后将对流热流计表面热流换算为试验模型表面辐射热流。
简而言之，本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法利用高焓膨胀风洞产生的极短时超高速流场，模拟超高速飞行器高温绕流，通过隔离气流对试验模型的对流加热，测量高温绕流对试验模型表面形成的辐射热流。试验方法首先根据飞行工况确定高焓膨胀风洞的模拟试验条件，再根据模拟试验条件估算辐射主导波段及辐射热流范围后设计制作瞬态辐射热流传感器，然后对瞬态辐射热流传感器的测量信号进行数据处理，得到试验模型表面测点处的辐射热流。本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法为高焓气动热环境研究提供了技术支持。
附图说明
图1为本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法使用的瞬态辐射热流传感器示意图；
图2为本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法使用的试验装置示意图；
图3为本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法获得的热流计表面热流曲线。
图中，1.试验段舱体；2.风洞喷管；3.瞬态辐射热流传感器；4.试验模型；5.模型支杆；6.攻角机构；7.数据信号线；8.数据采集设备；9.控制及数据处理计算机；10.柱形光学玻璃；11.辐射吸收层；12.安装套；13.对流热流计；14.信号线。
具体实施方式
下面结合附图和实施例详细说明本发明。
本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，包括以下步骤：
S1.确定高焓膨胀风洞试验条件；
S11.根据高焓膨胀风洞试验需要模拟的飞行高度和速度，确定高焓膨胀风洞的运行状态，并根据高焓膨胀风洞流场调试校测获得的流场参数，确定试验模型4尺寸；
S12.采用CFD计算试验模型4头部的激波层参数，估算试验模型4表面的辐射热流范围，设置瞬态辐射热流传感器3的位置，即测点位置；
根据激波层气流温度，假设测点的光辐射亮度
符合黑体辐射理论，按下式预估绕流的主要辐射波段：
式中，
为黑体温度，/>
为辐射波长，/>
为第一辐射常数，
为第二辐射常数；
计算辐射与波长的对应曲线，确定试验工况下，激波层气流对试验模型4表面的辐射主导波段；
试验模型4表面的辐射热流q r按下式估算：
其中，
，为常数；r n为模型头部等效球头半径，单位为m；/>
为气流密度，单位为kg/m3；V为气流速度，单位为m/s；a、b及f(V)按下式计算：
辐射热流q r的单位为W/cm2；a max表示a的最大值；
S13.根据高温绕流流场温度所对应的主要辐射波段及模型表面辐射热流，确定瞬态辐射热流传感器3表面的柱形光学玻璃10及内部的对流热流计13的参数及选型；
S2.安装试验装置；
S21.如图1所示，制作瞬态辐射热流传感器3；瞬态辐射热流传感器3为从外至内同中心轴线的双层柱体结构，包括位于外层的安装套12、位于内层的瞬态辐射热流传感器测量体，瞬态辐射热流传感器测量体包括从上至下顺序固定的柱形光学玻璃10和对流热流计13；对流热流计13的上表面设置辐射吸收层11；对流热流计13的下表面设置信号线14，信号线14穿出瞬态辐射热流传感器测量体，与数据采集设备8连接；
S22.如图2所示，在高焓膨胀风洞的试验段舱体1内，通过模型支杆5将试验模型4安装在攻角机构6上，试验模型4位于风洞喷管2的出口、试验段的流场均匀区内，在试验模型4的头部驻点安装瞬态辐射热流传感器3，瞬态辐射热流传感器3与试验模型4的表面平齐或者平滑过渡；瞬态辐射热流传感器3将高温流场的辐射光信号转变成电压信号，电压信号经数据信号线7传输至数据采集设备8，通过控制及数据处理计算机9的数据处理软件处理电压信号，得到试验模型4表面的辐射热流q r；
S23.在高焓膨胀风洞的洞体外，安装数据采集设备8的触发器，触发器以高焓膨胀风洞的激波管管壁的压电传感器发送的电压信号为触发信号；
S24.启动数据采集设备8和触发器，使得试验装置处于等待状态；
S3.进行高焓膨胀风洞试验；
启动高焓膨胀风洞，高焓膨胀风洞上游激波管中形成的入射激波作用于压电传感器，压电传感器输出电压信号触发数据采集设备8，数据采集设备8开始采集，高焓膨胀风洞试验结束后，停止采集；
S4.进行试验数据处理；
S41.计算对流热流计13表面热流曲线：
瞬态辐射热流传感器3输出的信号为对流热流计13感受到加热后输出的电压变化信号，将电压信号通过对流热流计13的电阻温度系数转换为温度后，采用基于一维半无限假设的热流计算方法，获得热流随时间变化曲线：
式中，
表示瞬态辐射热流传感器3的对流热流计13直接测量的对流热流计13表面热流；/>
为对流热流计13基底材料热物性综合参数；其中，/>
为基底材料的密度，单位为kg/m3；c为基底材料的比热，单位为J/（kg·K），k为基底材料的热传导系数，单位为w/（m·K）；T为温度；t为时间；
将热流随时间的变化曲线，结合高焓膨胀风洞流场的有效试验时间区间，以试验气流到达前的数据平均值作为基线值，有效试验时间的热流平均值作为测量平均值，测量平均值与基线值的之差为有效热流平均值；
S42.将对流热流计13表面热流换算为试验模型4表面辐射热流；
根据瞬态辐射热流传感器3的结构和辐射光的传输过程，将对流热流计13表面热流换算为试验模型4表面辐射热流：
式中，
为试验模型4表面辐射热流，单位为W/m2；/>
为对流热流计13表面热流，单位为W/m2；/>
为对流热流计13的敏感元件面积，单位为m2；/>
为瞬态辐射热流传感器3的柱形光学玻璃10的窗口面积，单位为m2；/>
为瞬态辐射热流传感器3的视角因子；/>
为对流热流计13前端面的辐射吸收层11的吸收率；/>
为瞬态辐射热流传感器3的柱形光学玻璃10的透射率；
其中，视角因子按下式计算：
式中，
为对流热流计13的敏感元件能够看到的柱形光学玻璃10的窗口半径，单位为m；/>
为对流热流计13的敏感元件到柱形光学玻璃10的窗口外表面的高度，单位为m。
进一步地，所述的瞬态辐射热流传感器3安装到试验模型4壳体上预留的安装通孔内，采用胶合剂从试验模型4内壁面一侧粘结固定。
进一步地，所述的柱形光学玻璃10采用MgF2玻璃。
进一步地，所述的对流热流计13为薄膜热流计或同轴热电偶。
进一步地，所述的辐射吸收层11的感应波长范围200nm～2500nm。
进一步地，所述的安装套12的材质与试验模型4的壳体的材质相同；或者安装套12的材质采用铝合金或者不锈钢中的一种；安装套12内壁涂敷有反射膜。
实施例1：
本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，在试验前根据高焓膨胀风洞的高焓试验工况预测高温气体辐射主导波段范围及预估辐射热流值，作为瞬态辐射热流传感器3设计参考。假设，针对高焓膨胀风洞气流速度10km/s的典型高焓状态，通过跨驻激波的理论计算方法或CFD仿真方法，可知钝头体形状的试验模型4头部弓形激波的激波层温度约10000K，基于黑体辐射假设，预估该温度下激波层主导辐射在近紫外-可见光-近红外波段，因此选择在该波段范围内具有稳定高透过率的MgF2作为柱形光学玻璃10的材料，并选择该波段范围内吸收率相对较高的铬硅膜作为对流热流计13上表面即感应面的辐射吸收层11；然后根据预估的辐射热流值范围以及流场有效试验时间确定对流热流计13的选型，在辐射热流较小且有效试验时间较短时，选择灵敏度较高、响应速度快的铂薄膜热流计作为对流热流计13。
本实施例的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，具体参数如下：
高焓膨胀风洞的来流速度约11.5km/s，试验模型4为类似阿波罗返回舱的钝头体，钝头体的头部具有弓形激波，能够参考前述的典型高焓状态，开展本实施例的研究工作。
瞬态辐射热流传感器3安装在钝头体迎风大底表面中心点及附近区域，与瞬态辐射热流传感器3匹配的钝头体的安装通孔直径为4.1±0.05mm，瞬态辐射热流传感器3的安装套12与钝头体的表面平齐。
瞬态辐射热流传感器3的安装套12外径4mm，长12mm，材质为304不锈钢，安装套12与内部的瞬态辐射热流传感器测量体的各部件之间使用环氧树脂密封和固定；柱形光学玻璃10为MgF2玻璃，直径3mm，厚度1mm，辐射透过率约90%；对流热流计13为直径2mm的柱状铂薄膜热流计，上表面的敏感元为S型铂薄膜，厚度0.1μm、面积
约0.805mm2；对流热流计13的上表面还镀制一层厚度0.2μm的铬硅膜作为辐射吸收层11，辐射吸收层11的辐射吸收系数取250nm～1000nm波段的初步标定结果，约为75%；柱形光学玻璃10的下表面与辐射吸收层11的上表面之间的间隙约0.15mm，柱形光学玻璃10的上表面与辐射吸收层11的上表面之间的间隙约1.05mm。
瞬态辐射热流传感器3的视角因子f v =0.0765，瞬态辐射热流传感器3的对流热流计13直接测量的热流
与试验模型4表面辐射热流/>
的换算关系为/>
。
本实施例获得了如图3所示的热流计表面热流曲线，图中的竖线Ⅰ、竖线Ⅱ之间的83μs为有效试验时间区间，该区间内的平均热流为36W/cm2，按本发明的用于高焓膨胀风洞模型表面辐射热流测量的试验方法的数据处理方法，得到试验模型4表面辐射热流
为80W/cm2。
尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域。对于熟悉本领域的人员而言，在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。
Claims (6)
1.一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法，其特征在于，包括以下步骤：
S1.确定高焓膨胀风洞试验条件；
S11.根据高焓膨胀风洞试验需要模拟的飞行高度和速度，确定高焓膨胀风洞的运行状态，并根据高焓膨胀风洞流场调试校测获得的流场参数，确定试验模型（4）尺寸；
S12.采用CFD计算试验模型（4）头部的激波层参数，估算试验模型（4）表面的辐射热流范围，设置瞬态辐射热流传感器（3）的位置，即测点位置；
根据激波层气流温度，假设测点的光辐射亮度
符合黑体辐射理论，按下式预估绕流的主要辐射波段：
式中，
为黑体温度，/>
为辐射波长，/>
为第一辐射常数，
为第二辐射常数；
计算辐射与波长的对应曲线，确定试验工况下，激波层气流对试验模型（4）表面的辐射主导波段；
试验模型（4）表面的辐射热流q r按下式估算：
其中，
，为常数；r n为模型头部等效球头半径，单位为m；/>
为气流密度，单位为kg/m3；V为气流速度，单位为m/s；a、b及f(V)按下式计算：
辐射热流q r的单位为W/cm2；a max表示a的最大值；
S13.根据高温绕流流场温度所对应的主要辐射波段及模型表面辐射热流，确定瞬态辐射热流传感器（3）表面的柱形光学玻璃（10）及内部的对流热流计（13）的参数及选型；
S2.安装试验装置；
S21.制作瞬态辐射热流传感器（3）；瞬态辐射热流传感器（3）为从外至内同中心轴线的双层柱体结构，包括位于外层的安装套（12）、位于内层的瞬态辐射热流传感器测量体，瞬态辐射热流传感器测量体包括从上至下顺序固定的柱形光学玻璃（10）和对流热流计（13）；对流热流计（13）的上表面设置辐射吸收层（11）；对流热流计（13）的下表面设置信号线（14），信号线（14）穿出瞬态辐射热流传感器测量体，与数据采集设备（8）连接；
S22.在高焓膨胀风洞的试验段舱体（1）内，通过模型支杆（5）将试验模型（4）安装在攻角机构（6）上，试验模型（4）位于风洞喷管（2）的出口、试验段的流场均匀区内，在试验模型（4）的头部驻点安装瞬态辐射热流传感器（3），瞬态辐射热流传感器（3）与试验模型（4）的表面平齐或者平滑过渡；瞬态辐射热流传感器（3）将高温流场的辐射光信号转变成电压信号，电压信号经数据信号线（7）传输至数据采集设备（8），通过控制及数据处理计算机（9）的数据处理软件处理电压信号，得到试验模型（4）表面的辐射热流q r；
S23.在高焓膨胀风洞的洞体外，安装数据采集设备（8）的触发器，触发器以高焓膨胀风洞的激波管管壁的压电传感器发送的电压信号为触发信号；
S24.启动数据采集设备（8）和触发器，使得试验装置处于等待状态；
S3.进行高焓膨胀风洞试验；
启动高焓膨胀风洞，高焓膨胀风洞上游激波管中形成的入射激波作用于压电传感器，压电传感器输出电压信号触发数据采集设备（8），数据采集设备（8）开始采集，高焓膨胀风洞试验结束后，停止采集；
S4.进行试验数据处理；
S41.计算对流热流计（13）表面热流曲线：
瞬态辐射热流传感器（3）输出的信号为对流热流计（13）感受到加热后输出的电压变化信号，将电压信号通过对流热流计（13）的电阻温度系数转换为温度后，采用基于一维半无限假设的热流计算方法，获得热流随时间变化曲线：
式中，
表示瞬态辐射热流传感器（3）的对流热流计（13）直接测量的对流热流计（13）表面热流；/>
为对流热流计（13）基底材料热物性综合参数；其中，/>
为基底材料的密度，单位为kg/m3；c为基底材料的比热，单位为J/（kg·K），k为基底材料的热传导系数，单位为w/（m·K）；T为温度；t为时间；
将热流随时间的变化曲线，结合高焓膨胀风洞流场的有效试验时间区间，以试验气流到达前的数据平均值作为基线值，有效试验时间的热流平均值作为测量平均值，测量平均值与基线值的之差为有效热流平均值；
S42.将对流热流计（13）表面热流换算为试验模型（4）表面辐射热流；
根据瞬态辐射热流传感器（3）的结构和辐射光的传输过程，将对流热流计（13）表面热流换算为试验模型（4）表面辐射热流：
式中，
为试验模型（4）表面辐射热流，单位为W/m2；/>
为对流热流计（13）表面热流，单位为W/m2；/>
为对流热流计（13）的敏感元件面积，单位为m2；/>
为瞬态辐射热流传感器（3）的柱形光学玻璃（10）的窗口面积，单位为m2；/>
为瞬态辐射热流传感器（3）的视角因子；/>
为对流热流计（13）前端面的辐射吸收层（11）的吸收率；/>
为瞬态辐射热流传感器（3）的柱形光学玻璃（10）的透射率；
其中，视角因子按下式计算：
式中，
为对流热流计（13）的敏感元件能够看到的柱形光学玻璃（10）的窗口半径，单位为m；/>
为对流热流计（13）的敏感元件到柱形光学玻璃（10）的窗口外表面的高度，单位为m。
2.根据权利要求1所述的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，其特征在于，所述的瞬态辐射热流传感器（3）安装到试验模型（4）壳体上预留的安装通孔内，采用胶合剂从试验模型（4）内壁面一侧粘结固定。
3.根据权利要求1所述的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，其特征在于，所述的柱形光学玻璃（10）采用MgF2玻璃。
4.根据权利要求1所述的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，其特征在于，所述的对流热流计（13）为薄膜热流计或同轴热电偶。
5.根据权利要求1所述的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，其特征在于，所述的辐射吸收层（11）的感应波长范围200nm～2500nm。
6.根据权利要求1所述的用于高焓膨胀风洞模型表面辐射热流测量的试验方法，其特征在于，所述的安装套（12）的材质与试验模型（4）的壳体的材质相同；或者安装套（12）的材质采用铝合金或者不锈钢中的一种；安装套（12）内壁涂敷有反射膜。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310604823.4A CN116296227B (zh) | 2023-05-26 | 2023-05-26 | 一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310604823.4A CN116296227B (zh) | 2023-05-26 | 2023-05-26 | 一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116296227A true CN116296227A (zh) | 2023-06-23 |
| CN116296227B CN116296227B (zh) | 2023-07-25 |
Family
ID=86822587
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310604823.4A Active CN116296227B (zh) | 2023-05-26 | 2023-05-26 | 一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116296227B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116698199A (zh) * | 2023-08-07 | 2023-09-05 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于原子层热电堆的辐射热流传感器及其安装试验方法 |
| CN117782515A (zh) * | 2024-02-28 | 2024-03-29 | 中国空气动力研究与发展中心计算空气动力研究所 | 激波风洞来流参数影响的气动热数据不确定度评估方法 |
| CN118376379A (zh) * | 2024-06-21 | 2024-07-23 | 中国空气动力研究与发展中心超高速空气动力研究所 | 高焓瞬时流场光辐射多角度多波段的同时测量装置及方法 |
| CN118533420A (zh) * | 2024-07-26 | 2024-08-23 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于光纤的嵌入式模型绕流光辐射测量装置与测量方法 |
| CN118583430A (zh) * | 2024-07-31 | 2024-09-03 | 中国航空工业集团公司沈阳空气动力研究所 | 一种常规高超声速风洞的瞬态热流测量装置及试验方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107843405A (zh) * | 2016-09-21 | 2018-03-27 | 北京空天技术研究所 | 试验件和发动机燃气对飞行器底部辐射热流的获取方法 |
| CN109029907A (zh) * | 2018-07-18 | 2018-12-18 | 大连理工大学 | 一种气动热环境试验模拟条件的参数相似方法 |
| AU2019205004B1 (en) * | 2019-07-11 | 2020-10-01 | Iyinomen, Daniel Odion DR | A Novel Plasma Preheating Test Device for Replicating Planetary Reentry Surface Temperatures in Hypersonic Impulse Facilities |
| CN112067241A (zh) * | 2020-08-19 | 2020-12-11 | 中国航天空气动力技术研究院 | 一种高焓激波风洞总温测量方法 |
| CN113945354A (zh) * | 2021-12-14 | 2022-01-18 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于辨识膨胀风洞加速段流动分区特性的试验方法 |
| CN114486155A (zh) * | 2021-12-27 | 2022-05-13 | 中国航天空气动力技术研究院 | 一种高焓激波风洞参数诊断方法和系统 |
| CN115452180A (zh) * | 2022-09-23 | 2022-12-09 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种高焓气流恢复温度测量方法及测量装置 |
-
2023
- 2023-05-26 CN CN202310604823.4A patent/CN116296227B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107843405A (zh) * | 2016-09-21 | 2018-03-27 | 北京空天技术研究所 | 试验件和发动机燃气对飞行器底部辐射热流的获取方法 |
| CN109029907A (zh) * | 2018-07-18 | 2018-12-18 | 大连理工大学 | 一种气动热环境试验模拟条件的参数相似方法 |
| US20200217749A1 (en) * | 2018-07-18 | 2020-07-09 | Dalian University Of Technology | Parmeter similarity method for test simulation conditions of aerodynamic heating environment |
| AU2019205004B1 (en) * | 2019-07-11 | 2020-10-01 | Iyinomen, Daniel Odion DR | A Novel Plasma Preheating Test Device for Replicating Planetary Reentry Surface Temperatures in Hypersonic Impulse Facilities |
| CN112067241A (zh) * | 2020-08-19 | 2020-12-11 | 中国航天空气动力技术研究院 | 一种高焓激波风洞总温测量方法 |
| CN113945354A (zh) * | 2021-12-14 | 2022-01-18 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于辨识膨胀风洞加速段流动分区特性的试验方法 |
| CN114486155A (zh) * | 2021-12-27 | 2022-05-13 | 中国航天空气动力技术研究院 | 一种高焓激波风洞参数诊断方法和系统 |
| CN115452180A (zh) * | 2022-09-23 | 2022-12-09 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种高焓气流恢复温度测量方法及测量装置 |
Non-Patent Citations (2)
| Title |
|---|
| 唐贵明, 袁生学, 曾远金: "激波风洞高焓流动及其驻点对流和辐射热流测量", 流体力学实验与测量, no. 01, pages 51 - 56 * |
| 朱广生;聂春生;曹占伟;袁野;: "气动热环境试验及测量技术研究进展", 实验流体力学, no. 02, pages 3 - 12 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116698199A (zh) * | 2023-08-07 | 2023-09-05 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于原子层热电堆的辐射热流传感器及其安装试验方法 |
| CN116698199B (zh) * | 2023-08-07 | 2023-10-03 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于原子层热电堆的辐射热流传感器及其安装试验方法 |
| CN117782515A (zh) * | 2024-02-28 | 2024-03-29 | 中国空气动力研究与发展中心计算空气动力研究所 | 激波风洞来流参数影响的气动热数据不确定度评估方法 |
| CN117782515B (zh) * | 2024-02-28 | 2024-05-07 | 中国空气动力研究与发展中心计算空气动力研究所 | 激波风洞来流参数影响的气动热数据不确定度评估方法 |
| CN118376379A (zh) * | 2024-06-21 | 2024-07-23 | 中国空气动力研究与发展中心超高速空气动力研究所 | 高焓瞬时流场光辐射多角度多波段的同时测量装置及方法 |
| CN118376379B (zh) * | 2024-06-21 | 2024-08-16 | 中国空气动力研究与发展中心超高速空气动力研究所 | 高焓瞬时流场光辐射多角度多波段的同时测量装置及方法 |
| CN118533420A (zh) * | 2024-07-26 | 2024-08-23 | 中国空气动力研究与发展中心超高速空气动力研究所 | 基于光纤的嵌入式模型绕流光辐射测量装置与测量方法 |
| CN118583430A (zh) * | 2024-07-31 | 2024-09-03 | 中国航空工业集团公司沈阳空气动力研究所 | 一种常规高超声速风洞的瞬态热流测量装置及试验方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116296227B (zh) | 2023-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116296227B (zh) | 一种用于高焓膨胀风洞模型表面辐射热流测量的试验方法 | |
| CN106908352B (zh) | 基于分布式温敏光纤的翼型表面边界层转捩位置测量方法 | |
| CN109632867A (zh) | 一种用于考核材料高超声速抗烧蚀性能的试验系统及方法 | |
| WO2024148899A1 (zh) | 激波风洞等离子鞘套与尾迹光辐射特性测量系统及方法 | |
| CN110174527A (zh) | 一种热式光纤光栅风速风向传感器及检测方法 | |
| CN109387349A (zh) | 一种不规则天线窗表面热流密度和压力测量装置 | |
| CN117871026A (zh) | 一种高焓流场辐射热流试验系统和方法 | |
| CN113390602B (zh) | 一种基于内嵌式温敏漆的气动热全场测量系统及方法 | |
| CN116380402B (zh) | 用于高焓膨胀风洞的瞬态辐射热流传感器的安装制作方法 | |
| CN211978276U (zh) | 一种用于高焓激波风洞燃烧场的热环境测量装置 | |
| CN109870406B (zh) | 一种材料表面涂层附着力测试方法及系统 | |
| CN114184348B (zh) | 一种高焓流场光电特性辨识装置及方法 | |
| Yang et al. | Multiphysical field analysis of a temperature sensor for meteorological observation | |
| CN206292249U (zh) | 一种热式风速传感单元及传感器 | |
| JPH09236494A (ja) | 放射誤差除去のための方法及び温度センサー装置 | |
| CN209524882U (zh) | 一种光纤式结冰探测仪 | |
| CN209102139U (zh) | 一种脱落插座热流密度和压力测量装置 | |
| CN116698199B (zh) | 基于原子层热电堆的辐射热流传感器及其安装试验方法 | |
| CN209102329U (zh) | 一种不规则天线窗表面热流密度和压力测量装置 | |
| Mani et al. | Experimental investigation of blunt cone model at hypersonic Mach number 7.25 | |
| Lee et al. | Compensation of solar radiation and ventilation effects on the temperature measurement of radiosondes using dual thermistors | |
| Wanqing et al. | Research and development of surface heat flux sensor for high speed aircraft | |
| CN110186572A (zh) | 一种具有探头预热功能的耳温枪 | |
| Schülein | Experimental investigation of laminar flow control on a supersonic swept wing by suction | |
| CN113758599B (zh) | 一种用于动态总温测量的光纤法-珀总温探针及其制作方法 |
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 |