CN113418675B - 一种热线测量风洞流场扰动模态方法 - Google Patents
一种热线测量风洞流场扰动模态方法 Download PDFInfo
- Publication number
- CN113418675B CN113418675B CN202110900462.9A CN202110900462A CN113418675B CN 113418675 B CN113418675 B CN 113418675B CN 202110900462 A CN202110900462 A CN 202110900462A CN 113418675 B CN113418675 B CN 113418675B
- Authority
- CN
- China
- Prior art keywords
- flow field
- hot wire
- disturbance
- value
- mode
- 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
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005259 measurement Methods 0.000 title abstract description 16
- 230000010349 pulsation Effects 0.000 claims abstract description 30
- 239000000523 sample Substances 0.000 claims abstract description 6
- 230000035945 sensitivity Effects 0.000 claims abstract description 4
- 238000010586 diagram Methods 0.000 claims description 13
- 238000013021 overheating Methods 0.000 claims description 6
- 238000011084 recovery Methods 0.000 claims description 2
- 238000000691 measurement method Methods 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 230000003068 static effect Effects 0.000 description 5
- 238000000926 separation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000000917 particle-image velocimetry Methods 0.000 description 1
- 230000000007 visual effect Effects 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/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
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种热线测量风洞流场扰动模态方法,属于流场测试技术领域。本发明根据热丝工作电压的无量纲脉动值θ与热丝灵敏度系数无量纲值r之间的拟合关系曲线特征,判断流场扰动的四种模态:涡、熵、声学和混合模态。本发明的一种热线测量风洞流场扰动模态方法,探头安装简单易行,测量方法明确,可以准确测量流场的扰动模态,为流场扰动模态评估和流场品质优化提供依据。
Description
技术领域
本发明属于流场测试技术领域,具体涉及一种热线测量风洞流场扰动模态方法。
背景技术
风洞流场在不同的流速、温度和洞体结构条件下,存在以下四种扰动模态:当扰动以温度脉动为主导时,称为熵模态;当扰动以速度脉动为主导时,称为涡模态;当扰动以静压脉动为主导时,称为声学模态;当无明显突出脉动量时,称为混合模态。测量风洞流场扰动模态,可以掌握风洞流场扰动的主导因素和流动机理,进而对风洞流场品质进行评估和针对性优化。
现有技术主要采用脉动压力传感器或者麦克风测量风洞静压脉动,采用总温探头测量风洞温度脉动,采用粒子图像测速法等技术测量亚声速流场涡量,但是,无统一方法同时测量以上四种风洞流场扰动模态,无法测量流场的主导脉动量和模态变化特征。
基于热线的风洞流场扰动模态测量方法,采用连续变热线风速仪(简称热线)过热比方法,测量获得流场扰动图,根据流场扰动图特征,可分析获得亚跨超声速流场的四种扰动模态、主导脉动量和模态变化特征,测量方法直观可靠。
发明内容
本发明要解决的技术问题是提供一种热线测量风洞流场扰动模态方法,探头安装时,使热丝垂直于来流方向,采用连续变热线过热比测量方法,获得热丝在连续变过热比下的输出电压值信号,根据这一组数据作出电压值与热丝灵敏度系数之间拟合曲线,根据拟合曲线特征确定流场的扰动模态。
本发明解决上述技术问题采用的技术方案:一种热线测量风洞流场扰动模态方法,该方法包括:
步骤1:将一维热线探头的热丝垂直于来流方向安装于流场核心流中;
步骤2:在流场中连续改变热丝过热比获得一组热线仪输出电压数据,热丝过热比a=(Rw-R0)/Rw,其中,Rw为热线在流场中加热时的工作电阻,R0为热线在流场中的非工作电阻;
步骤3:根据公式1获得热丝工作电压的无量纲脉动值θ与热丝灵敏度系数无量纲值r之间的拟合关系曲线,横坐标值为r,纵坐标值为θ;
其中,
r=FCTA/GCTA
式中,e′是热线过热比下的工作电压值均方根偏差;是热线过热比下的工作电压均值;T′是流场总温脉动均方根偏差;是流场总温脉动均值;是流场总温脉动值;ρ是流场密度,ρ′是流场密度均方根偏差,是流场密度均值;u是流场速度,u′是流场速度均方根偏差,是流场速度均值;是流场质量流脉动值;aw是热丝过热比;A和B是热线校准常系数;a*是热丝的电阻温度系数;R*是热丝参考电阻值;η是热丝恢复系数;Re是热丝未通电加热时在流场中的电阻;Re是流场雷诺数;T0是气流总温;k是热线变过热比函数的斜率;
步骤4:根据θ与r的拟合曲线特征,判断流场扰动模态;
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为熵模态,如图1所示:拟合曲线为直线,直线与横坐标相交于-α点附近;
式中,γ是比热比,M是马赫数;
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为涡模态,如图2所示:拟合曲线为V型折线,折线与横坐标相交于β点附近;
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为声学模态,如图3所示:拟合曲线为直线,斜率为正,斜线与纵坐标交点在坐标系原点附近并为正值,或者斜线与纵坐标相交时穿过坐标原点。
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为混合模态,如图4所示:拟合曲线为双曲线。
本发明的有益效果:
本发明的一种热线测量风洞流场扰动模态方法,探头安装简单易行,测量方法明确,可以有效解决流场扰动模态测量问题,为流场扰动模态评估和流场品质优化提供依据。
本发明为开展热线测试技术在流场扰动模态测量中的应用提供了设计参考。
附图说明
图1流场熵模态示意图;
图2流场涡模态示意图;
图3流场声学模态示意图;
图4流场混合模态示意图;
图5流场熵模态实测图;
图6流场涡模态实测图;
图7流场声学模态实测图;
图8流场混合模态实测图。
具体实施方式
图5是流场熵模态实测图,图6是流场涡模态实测图,图7是流场声学模态实测图,图8是流场混合模态实测图。
案例1:在某低温风洞,流场速度为马赫数0.5,总温为110K,测量过程中按照0.7、0.8、0.9、1.0、1.2、1.4、1.6、1.8过热比采集8组过热比数据,根据公式1拟合出横坐标值为r,纵坐标值为θ的脉动特征曲线图。测量获得图5所示的流场熵模态图,该图为拟合曲线为直线,直线与横坐标相交于-α点附近,流场扰动为温度脉动主导。为提高流场质量,重点对流场温度均匀性进行优化。
案例2:在某跨声速风洞试验段,流场速度为马赫数0.8,测量过程中按照0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8过热比采集8组过热比数据,根据公式1拟合出横坐标值为r,纵坐标值为θ的脉动特征曲线图。测量获得图6所示的流场涡模态图,该图为典型V折线,可判断流场扰动为速度脉动主导。为提高流场质量,重点抑制涡扰动,对试验段入口整流装置进行优化,抑制流动分离涡形成。
案例3:在某低速低湍流开口风洞,流场速度为80米/秒,测量过程中按照0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8过热比采集8组过热比数据,根据公式1拟合出横坐标值为r,纵坐标值为θ的脉动特征曲线图。测量获得图7所示的流场声学模态图,该图为接近原点的直线,可判断流场速度脉动很低,扰动为静压脉动主导。为提高流场质量,重点抑制低频压力脉动,采用导流片对风洞喷口整流方式进行优化。
案例4:在某跨声速风洞,流场速度为马赫数0.7,测量过程中按照0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8过热比采集8组过热比数据,根据公式1拟合出横坐标值为r,纵坐标值为θ的脉动特征曲线图。测量获得图8所示的流场混合模态图,该图的拟合曲线为双曲线,流场扰动无突出脉动量,流场有两种以上量级相当的扰动量。为提高流场质量,需同时考虑抑制流场的流动分离、静压脉动及两者的耦合。
Claims (1)
1.一种热线测量风洞流场扰动模态方法,该方法包括:
步骤1:将一维热线探头的热丝垂直于来流方向安装于流场核心流中;
步骤2:在流场中连续改变热丝过热比获得一组热线仪输出电压数据,热丝过热比a=(Rw-R0)/Rw,其中,Rw为热线在流场中加热时的工作电阻,R0为热线在流场中的非工作电阻;
步骤3:根据公式1获得热丝工作电压的无量纲脉动值θ与热丝灵敏度系数无量纲值r之间的拟合关系曲线,横坐标值为r,纵坐标值为θ;
其中,
r=FCTA/GCTA
式中,e′是热线过热比下的工作电压值均方根偏差;是热线过热比下的工作电压均值;T′是流场总温脉动均方根偏差;是流场总温脉动均值;是流场总温脉动值;ρ是流场密度,ρ′是流场密度均方根偏差,是流场密度均值;u是流场速度,u′是流场速度均方根偏差,是流场速度均值;是流场质量流脉动值;aw是热丝过热比;A和B是热线校准常系数;a*是热丝的电阻温度系数;R*是热丝参考电阻值;η是热丝恢复系数;Re是热丝未通电加热时在流场中的电阻;Re是流场雷诺数;T0是气流总温;k是热线变过热比函数的斜率;
步骤4:根据θ与r的拟合曲线特征,判断流场扰动模态;
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为熵模态:拟合曲线为直线,直线与横坐标相交于-α点附近;
式中,γ是比热比,M是马赫数;
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为涡模态:拟合曲线为V型折线,折线与横坐标相交于β点附近;
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为声学模态:拟合曲线为直线,斜率为正,斜线与纵坐标交点在坐标系原点附近并为正值,或者斜线与纵坐标相交时穿过坐标原点;
当扰动图中的拟合曲线具备以下特征时,判定流场扰动为混合模态:拟合曲线为双曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110900462.9A CN113418675B (zh) | 2021-08-06 | 2021-08-06 | 一种热线测量风洞流场扰动模态方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110900462.9A CN113418675B (zh) | 2021-08-06 | 2021-08-06 | 一种热线测量风洞流场扰动模态方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113418675A CN113418675A (zh) | 2021-09-21 |
CN113418675B true CN113418675B (zh) | 2022-12-02 |
Family
ID=77719784
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110900462.9A Active CN113418675B (zh) | 2021-08-06 | 2021-08-06 | 一种热线测量风洞流场扰动模态方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113418675B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114818369B (zh) * | 2022-05-19 | 2023-03-03 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种连续式跨声速风洞部段设计方法及系统及装置及介质 |
CN116358823B (zh) * | 2023-05-26 | 2023-08-18 | 中国空气动力研究与发展中心高速空气动力研究所 | 高速风洞自由来流质量流量和总温脉动不确定度评估方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1670491A (zh) * | 2005-05-16 | 2005-09-21 | 中国科学院力学研究所 | 一种超声速流场中有源声波检测方法 |
WO2006076252A2 (en) * | 2005-01-10 | 2006-07-20 | Semih Olcmen | Spark-plug ldv, lif, and lii probe for engine flow and combustion analysis |
WO2013052932A1 (en) * | 2011-10-06 | 2013-04-11 | Ofs Fitel, Llc | Broadband fiber sensor array |
WO2014076419A1 (fr) * | 2012-11-19 | 2014-05-22 | Eveon | Dispositif de mesure de débit d'un fluide dans un canal, dispositif de transfert d'un fluide et procédé de mesure de débit d'un fluide dans un canal |
CN104048808A (zh) * | 2013-03-14 | 2014-09-17 | 中国科学院工程热物理研究所 | 一种动态熵探针 |
CN106289506A (zh) * | 2016-09-06 | 2017-01-04 | 大连理工大学 | 一种使用pod分解法消除流场壁面麦克风阵列噪声信号的方法 |
CN110455490A (zh) * | 2019-08-21 | 2019-11-15 | 中国人民解放军国防科技大学 | 超声速和高超声速风洞流场湍流度的计算方法及装置 |
CN111624363A (zh) * | 2020-05-27 | 2020-09-04 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种考虑流体压缩性影响的热线风速仪校准方法 |
CN111707439A (zh) * | 2020-07-10 | 2020-09-25 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种可压缩流体湍流度测量试验数据的双曲线拟合方法 |
CN111766039A (zh) * | 2020-07-10 | 2020-10-13 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种亚声速风洞可压缩流体扰动模态测量结果解算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6948361B2 (en) * | 2001-06-12 | 2005-09-27 | Dana Corporation | Gasket flow sensing apparatus and method |
CN101614753A (zh) * | 2008-06-26 | 2009-12-30 | 清华大学 | 流场传感器及其制作方法 |
CN102262166B (zh) * | 2011-04-08 | 2013-01-09 | 内蒙古农业大学 | 风洞风速廓线仪及其控制方法 |
CN109632236A (zh) * | 2018-11-28 | 2019-04-16 | 中国辐射防护研究院 | 一种环境风洞中同步自动化流场采样与测量系统和方法 |
US20210223281A1 (en) * | 2019-11-20 | 2021-07-22 | Board Of Regents, The University Of Texas System | Velocity Measurements Using a Piezoelectric Sensor |
-
2021
- 2021-08-06 CN CN202110900462.9A patent/CN113418675B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006076252A2 (en) * | 2005-01-10 | 2006-07-20 | Semih Olcmen | Spark-plug ldv, lif, and lii probe for engine flow and combustion analysis |
CN1670491A (zh) * | 2005-05-16 | 2005-09-21 | 中国科学院力学研究所 | 一种超声速流场中有源声波检测方法 |
WO2013052932A1 (en) * | 2011-10-06 | 2013-04-11 | Ofs Fitel, Llc | Broadband fiber sensor array |
WO2014076419A1 (fr) * | 2012-11-19 | 2014-05-22 | Eveon | Dispositif de mesure de débit d'un fluide dans un canal, dispositif de transfert d'un fluide et procédé de mesure de débit d'un fluide dans un canal |
CN104048808A (zh) * | 2013-03-14 | 2014-09-17 | 中国科学院工程热物理研究所 | 一种动态熵探针 |
CN106289506A (zh) * | 2016-09-06 | 2017-01-04 | 大连理工大学 | 一种使用pod分解法消除流场壁面麦克风阵列噪声信号的方法 |
CN110455490A (zh) * | 2019-08-21 | 2019-11-15 | 中国人民解放军国防科技大学 | 超声速和高超声速风洞流场湍流度的计算方法及装置 |
CN111624363A (zh) * | 2020-05-27 | 2020-09-04 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种考虑流体压缩性影响的热线风速仪校准方法 |
CN111707439A (zh) * | 2020-07-10 | 2020-09-25 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种可压缩流体湍流度测量试验数据的双曲线拟合方法 |
CN111766039A (zh) * | 2020-07-10 | 2020-10-13 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种亚声速风洞可压缩流体扰动模态测量结果解算方法 |
Non-Patent Citations (4)
Title |
---|
Comparison of new ultrasonic and hot wire thermo-anemometer gas flow meters;Robert D. Garrick;《2011 Future of Instrumentation International Workshop (FIIW) Proceedings》;20111231;全文 * |
Wind tunnel model of the forest and its reynoles number sensitivity;Christof Gromke;《Journal of wind engineering and industrial aerodynamics》;20180430;全文 * |
亚声速风洞可压缩流体扰动模态分析;杜钰锋;《航空学报》;20210630;全文 * |
高超声速风洞来流扰动测量及数据后处理技术研究;余涛;《实验流体力学》;20191030;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113418675A (zh) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113418675B (zh) | 一种热线测量风洞流场扰动模态方法 | |
CN105157948B (zh) | 一种适用于超声速/高超声速流道的流量测试系统及测试方法 | |
Goldstein et al. | Laminar separation, reattachment, and transition of the flow over a downstream-facing step | |
CN102353410B (zh) | 一种利用示踪气体测量风道风量的方法及装置 | |
CN105987773B (zh) | 阻滞式总温传感器 | |
US20140208755A1 (en) | Gas Turbine Air Mass Flow Measuring System and Methods for Measuring Air Mass Flow in a Gas Turbine Inlet Duct | |
JPS5882121A (ja) | 流体速度測定装置 | |
CN104819806B (zh) | 一种高精度传感器校验装置 | |
CN105465860A (zh) | 一种吸油烟机实际流量测量装置及其测量方法 | |
CN108195510A (zh) | 一种热膜式剪应力传感器的热风洞校准方法 | |
CN115435930A (zh) | 一种测量级间三维全参数高频探针 | |
CN218822484U (zh) | 一种测量级间三维全参数高频探针 | |
CN115950493A (zh) | 一种适用于亚声速流道的流量测试系统及测试方法 | |
CN106840284A (zh) | 基于多管式皮托管矩阵式流量系统及测量方法 | |
CN113607377B (zh) | 一种一维热线探头测量亚音速管流上下游噪声源的方法 | |
CN109238738B (zh) | 基于可测试性格栅的特种车辆冷却风量测试装置及方法 | |
CN108304603A (zh) | 一种高速飞行器强制转捩装置验证方法 | |
CN106126858A (zh) | 一种动力电池风冷系统中风扇选型的方法 | |
CN216116185U (zh) | 一种电磁差压式质量流量计 | |
CN206709908U (zh) | 一种弯管流量计 | |
CN201335843Y (zh) | 双喉径文丘里风速测量管 | |
CN210741871U (zh) | 一种冷却器模拟测试平台 | |
CN111677683B (zh) | 一种基于流量补偿法的微型风扇气动性能测试方法及装置 | |
CN209385257U (zh) | 一种发动机egr孔板流量计 | |
CN106841672B (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 |