CN113418675A - 一种热线测量风洞流场扰动模态方法 - Google Patents

一种热线测量风洞流场扰动模态方法 Download PDF

Info

Publication number
CN113418675A
CN113418675A CN202110900462.9A CN202110900462A CN113418675A CN 113418675 A CN113418675 A CN 113418675A CN 202110900462 A CN202110900462 A CN 202110900462A CN 113418675 A CN113418675 A CN 113418675A
Authority
CN
China
Prior art keywords
flow field
hot wire
disturbance
value
flow
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
CN202110900462.9A
Other languages
English (en)
Other versions
CN113418675B (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.)
Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center
Original Assignee
Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center
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 Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center filed Critical Equipment Design and Testing Technology Research Institute of China Aerodynamics Research and Development Center
Priority to CN202110900462.9A priority Critical patent/CN113418675B/zh
Publication of CN113418675A publication Critical patent/CN113418675A/zh
Application granted granted Critical
Publication of CN113418675B publication Critical patent/CN113418675B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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，纵坐标值为θ；
Figure BDA0003199589620000011
其中，
Figure BDA0003199589620000021
r＝FCTA/GCTA
Figure BDA0003199589620000022
Figure BDA0003199589620000023
Figure BDA0003199589620000024
式中，e′是热线过热比下的工作电压值均方根偏差；
Figure BDA0003199589620000025
是热线过热比下的工作电压均值；T′是流场总温脉动均方根偏差；
Figure BDA0003199589620000026
是流场总温脉动均值；
Figure BDA0003199589620000027
是流场总温脉动值；ρ是流场密度，ρ′是流场密度均方根偏差，
Figure BDA0003199589620000028
是流场密度均值；u是流场速度，u′是流场速度均方根偏差，
Figure BDA0003199589620000029
是流场速度均值；
Figure BDA00031995896200000210
是流场质量流脉动值；aw是热丝过热比；A和B是热线校准常系数；a*是热丝的电阻温度系数；R*是热丝参考电阻值；η是热丝恢复系数；Re是热丝未通电加热时在流场中的电阻；Re是流场雷诺数；T0是气流总温；k是热线变过热比函数
Figure BDA00031995896200000211
的斜率；
步骤4：根据θ与r的拟合曲线特征，判断流场扰动模态；
当扰动图中的拟合曲线具备以下特征时，判定流场扰动为熵模态，如图1所示：拟合曲线为直线，直线与横坐标相交于-α点附近；
Figure BDA00031995896200000212
式中，γ是比热比，M是马赫数；
当扰动图中的拟合曲线具备以下特征时，判定流场扰动为涡模态，如图2所示：拟合曲线为V型折线，折线与横坐标相交于β点附近；
Figure BDA0003199589620000031
当扰动图中的拟合曲线具备以下特征时，判定流场扰动为声学模态，如图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，纵坐标值为θ；
Figure FDA0003199589610000011
其中，
Figure FDA0003199589610000012
r＝FCTA/GCTA
Figure FDA0003199589610000013
Figure FDA0003199589610000014
Figure FDA0003199589610000015
式中，e′是热线过热比下的工作电压值均方根偏差；
Figure FDA0003199589610000016
是热线过热比下的工作电压均值；T′是流场总温脉动均方根偏差；
Figure FDA0003199589610000017
是流场总温脉动均值；
Figure FDA0003199589610000018
是流场总温脉动值；ρ是流场密度，ρ′是流场密度均方根偏差，
Figure FDA0003199589610000019
是流场密度均值；u是流场速度，u′是流场速度均方根偏差，
Figure FDA00031995896100000110
是流场速度均值；
Figure FDA00031995896100000111
是流场质量流脉动值；aw是热丝过热比；A和B是热线校准常系数；a*是热丝的电阻温度系数；R*是热丝参考电阻值；η是热丝恢复系数；Re是热丝未通电加热时在流场中的电阻；Re是流场雷诺数；T0是气流总温；k是热线变过热比函数
Figure FDA0003199589610000021
的斜率；
步骤4：根据θ与r的拟合曲线特征，判断流场扰动模态；
当扰动图中的拟合曲线具备以下特征时，判定流场扰动为熵模态：拟合曲线为直线，直线与横坐标相交于-α点附近；
Figure FDA0003199589610000022
式中，γ是比热比，M是马赫数；
当扰动图中的拟合曲线具备以下特征时，判定流场扰动为涡模态：拟合曲线为V型折线，折线与横坐标相交于β点附近；
Figure FDA0003199589610000023
当扰动图中的拟合曲线具备以下特征时，判定流场扰动为声学模态：拟合曲线为直线，斜率为正，斜线与纵坐标交点在坐标系原点附近并为正值，或者斜线与纵坐标相交时穿过坐标原点。
当扰动图中的拟合曲线具备以下特征时，判定流场扰动为混合模态：拟合曲线为双曲线。
CN202110900462.9A 2021-08-06 2021-08-06 一种热线测量风洞流场扰动模态方法 Active CN113418675B (zh)

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 true CN113418675A (zh) 2021-09-21
CN113418675B 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)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114818369A (zh) * 2022-05-19 2022-07-29 中国空气动力研究与发展中心设备设计与测试技术研究所 一种连续式跨声速风洞部段设计方法及系统及装置及介质
CN116358823A (zh) * 2023-05-26 2023-06-30 中国空气动力研究与发展中心高速空气动力研究所 高速风洞自由来流质量流量和总温脉动不确定度评估方法

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040144169A1 (en) * 2001-06-12 2004-07-29 Popielas Frank W. Gasket flow sensing apparatus and method
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
CN101614753A (zh) * 2008-06-26 2009-12-30 清华大学 流场传感器及其制作方法
CN102262166A (zh) * 2011-04-08 2011-11-30 内蒙古农业大学 风洞风速廓线仪及其控制方法
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分解法消除流场壁面麦克风阵列噪声信号的方法
CN109632236A (zh) * 2018-11-28 2019-04-16 中国辐射防护研究院 一种环境风洞中同步自动化流场采样与测量系统和方法
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 中国空气动力研究与发展中心高速空气动力研究所 一种亚声速风洞可压缩流体扰动模态测量结果解算方法
US20210223281A1 (en) * 2019-11-20 2021-07-22 Board Of Regents, The University Of Texas System Velocity Measurements Using a Piezoelectric Sensor

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040144169A1 (en) * 2001-06-12 2004-07-29 Popielas Frank W. Gasket flow sensing apparatus and method
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 中国科学院力学研究所 一种超声速流场中有源声波检测方法
CN101614753A (zh) * 2008-06-26 2009-12-30 清华大学 流场传感器及其制作方法
CN102262166A (zh) * 2011-04-08 2011-11-30 内蒙古农业大学 风洞风速廓线仪及其控制方法
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分解法消除流场壁面麦克风阵列噪声信号的方法
CN109632236A (zh) * 2018-11-28 2019-04-16 中国辐射防护研究院 一种环境风洞中同步自动化流场采样与测量系统和方法
CN110455490A (zh) * 2019-08-21 2019-11-15 中国人民解放军国防科技大学 超声速和高超声速风洞流场湍流度的计算方法及装置
US20210223281A1 (en) * 2019-11-20 2021-07-22 Board Of Regents, The University Of Texas System Velocity Measurements Using a Piezoelectric Sensor
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)

* Cited by examiner, † Cited by third party
Title
CHRISTOF GROMKE: "Wind tunnel model of the forest and its reynoles number sensitivity", 《JOURNAL OF WIND ENGINEERING AND INDUSTRIAL AERODYNAMICS》 *
ROBERT D. GARRICK: "Comparison of new ultrasonic and hot wire thermo-anemometer gas flow meters", 《2011 FUTURE OF INSTRUMENTATION INTERNATIONAL WORKSHOP (FIIW) PROCEEDINGS》 *
余涛: "高超声速风洞来流扰动测量及数据后处理技术研究", 《实验流体力学》 *
杜钰锋: "亚声速风洞可压缩流体扰动模态分析", 《航空学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114818369A (zh) * 2022-05-19 2022-07-29 中国空气动力研究与发展中心设备设计与测试技术研究所 一种连续式跨声速风洞部段设计方法及系统及装置及介质
CN114818369B (zh) * 2022-05-19 2023-03-03 中国空气动力研究与发展中心设备设计与测试技术研究所 一种连续式跨声速风洞部段设计方法及系统及装置及介质
CN116358823A (zh) * 2023-05-26 2023-06-30 中国空气动力研究与发展中心高速空气动力研究所 高速风洞自由来流质量流量和总温脉动不确定度评估方法
CN116358823B (zh) * 2023-05-26 2023-08-18 中国空气动力研究与发展中心高速空气动力研究所 高速风洞自由来流质量流量和总温脉动不确定度评估方法

Also Published As

Publication number Publication date
CN113418675B (zh) 2022-12-02

Similar Documents

Publication Publication Date Title
CN113418675B (zh) 一种热线测量风洞流场扰动模态方法
Goldstein et al. Laminar separation, reattachment, and transition of the flow over a downstream-facing step
US8056424B2 (en) Multi-channel flow sensor with extended flow range and faster response
CN104848904A (zh) 进气道流量测量系统
CN102353410B (zh) 一种利用示踪气体测量风道风量的方法及装置
CN105091956A (zh) 发动机进气流量测量装置及发动机系统
CN108195510A (zh) 一种热膜式剪应力传感器的热风洞校准方法
CN110487338A (zh) 一种多孔平衡孔板设计方法及其评价方法
CN218822484U (zh) 一种测量级间三维全参数高频探针
CN106840284A (zh) 基于多管式皮托管矩阵式流量系统及测量方法
CN113607377B (zh) 一种一维热线探头测量亚音速管流上下游噪声源的方法
CN111498141A (zh) 一种基于微型探针实现气流角度实时监测的方法与装置
CN115435931B (zh) 一种能测量级间气流方向的高频熵探针
CN106126858A (zh) 一种动力电池风冷系统中风扇选型的方法
CN206709908U (zh) 一种弯管流量计
CN108304603A (zh) 一种高速飞行器强制转捩装置验证方法
CN115435930A (zh) 一种测量级间三维全参数高频探针
Wei et al. Scale effect during the dynamic stall of dynamic airfoil
CN201335843Y (zh) 双喉径文丘里风速测量管
CN205808480U (zh) 一种止回阀式计量控制装置
CN106841672B (zh) 阵列小机翼检测锅炉二次分风道风速的装置和方法
CN210741871U (zh) 一种冷却器模拟测试平台
CN208334405U (zh) 台风环境下提供参考静压、风速和风向的实测装置
CN115839817B (zh) 一种用于翼型试验的近壁多点式附面层测量装置及方法
Chondrokostas et al. Calibration of pneumatic five-hole probes in the free-jet wind tunnel

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