CN113901591B - 热线单元干功率计算方法及液态水含量计算方法 - Google Patents
热线单元干功率计算方法及液态水含量计算方法 Download PDFInfo
- Publication number
- CN113901591B CN113901591B CN202111492134.6A CN202111492134A CN113901591B CN 113901591 B CN113901591 B CN 113901591B CN 202111492134 A CN202111492134 A CN 202111492134A CN 113901591 B CN113901591 B CN 113901591B
- Authority
- CN
- China
- Prior art keywords
- hot wire
- wire unit
- dry air
- liquid water
- cloud
- 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
- 239000007788 liquid Substances 0.000 title claims abstract description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 44
- 238000004364 calculation method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 30
- 239000003595 mist Substances 0.000 claims abstract description 25
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 230000008020 evaporation Effects 0.000 claims description 11
- 238000001704 evaporation Methods 0.000 claims description 11
- 230000003068 static effect Effects 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 abstract description 15
- 238000005507 spraying Methods 0.000 description 6
- 239000002245 particle Substances 0.000 description 4
- 239000008276 ice cloud Substances 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 230000003116 impacting effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64D—EQUIPMENT FOR FITTING IN OR TO AIRCRAFT; FLIGHT SUITS; PARACHUTES; ARRANGEMENT OR MOUNTING OF POWER PLANTS OR PROPULSION TRANSMISSIONS IN AIRCRAFT
- B64D15/00—De-icing or preventing icing on exterior surfaces of aircraft
- B64D15/20—Means for detecting icing or initiating de-icing
-
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/56—Investigating or analyzing materials by the use of thermal means by investigating moisture content
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Automation & Control Theory (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Fluid Mechanics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明适用于结冰风洞试验技术领域,提供了一种热线单元干功率计算方法及基于此的液态水含量计算方法,本发明先在干空气条件下标定,得到热线单元干空气努塞尔数与热线雷诺数、干空气普朗特数的拟合关系,利用该关系可以计算结冰云雾条件下热线单元的干功率,再基于干功率进一步计算液态水含量。本发明的这种方法可以实现结冰风洞中结冰云雾液态水含量的精细化计算,有效地克服了现有技术测量精度较低的技术问题,为大型结冰风洞结冰云雾精细化测量提供了重要技术支撑。
Description
技术领域
本发明涉及结冰风洞试验技术领域,尤其是涉及一种热线单元干功率计算方法及基于此的液态水含量计算方法。
背景技术
当飞机在云层中飞行时,云层中的过冷水滴(即温度低于冰点的液态水滴)会不断撞击飞机迎风面,导致其表面发生结冰现象。飞机结冰广泛存在于飞行实践中,并严重威胁飞行安全。
结冰风洞是开展飞机结冰研究及飞机部件防除冰系统验证的重要地面试验设备,其在飞机结冰适航审定中扮演着重要角色。结冰云雾环境模拟能力是结冰风洞性能的核心内容。结冰云雾微物理特征的准确测量与评估是结冰风洞准确模拟结冰云雾环境的关键,其中云雾水含量(定义为结冰云雾中,单位体积内云雾颗粒质量)是结冰云雾重要的微物理特征之一。
基于恒温热线原理的云雾液态水含量测量方法是目前大型结冰风洞中重要的液态水含量测量方法。该原理将热线单元保持在恒定温度(一般在100℃以上),进而使撞击到热线单元表面的液滴完全蒸发,通过测量此时热线单元所需的功率,计算得到云雾液态水含量。热线单元功率由两部分组成,分别为补偿干空气对流换热损耗的功率(干功率)和补偿云雾蒸发损耗的功率(湿功率),其中湿功率直接用于计算云雾液态水含量。由此可见,准确获得热线单元干功率是准确计算云雾液态水含量的关键。
目前,通常采用两种方法获得热线单元干功率:(1)去皮法。测量喷雾前(无云雾状态)热线单元功率,同时假设喷雾后,干功率不发生变化,将此功率作为整个试验过程中的热线干功率。(2)补偿法。该方法通常在测量热线单元后部增加一根补偿热线单元,通过巧妙设计补偿探头位置,使得云雾颗粒不会撞击到补偿热线表面,使得该热线仅能测量干空气功率,进而标定获得测量热线与补偿热线干功率间的函数关系,最终在实验过程中,通过补偿热线单元测量功率计算测量热线单元的干功率。
但是,以上两种方法均存在明显缺陷。对于皮法,实验过程气流参数难免会发生变化,因此该方法核心假设难以完全成立,显著影响测量精度。对于补偿法,尽管通过增加补偿热线,可以有效捕捉到实验过程中流场参数变化造成的干功率的变化,但是实验过程中,测量热线单元的温度以及其与补偿热线间的标定关系通常会发生一定变化,这些变化会降低该方法得到的干功率可靠性。
因此,为系统解决大型结冰风洞中复杂结冰云雾液态水含量的精细化测量问题,突破热线水含量测量仪干功率及液态水含量精细化计算瓶颈,亟需发展一种全面考虑热线单元性能参数以及流场条件的干空气功率计算方法,从而提高云雾液态水含量计算精准度。
发明内容
为解决现有技术的不足,本发明针对大型结冰风洞中结冰云雾,建立了基于牛顿冷却方程的干功率标定和计算方法,提出了一种基于热线原理的液态水含量精细化计算方法。
本发明是这样实现的,一种热线单元干功率计算方法,其特征在于,包括以下步骤:
S10. 实验采集热线单元干空气条件标定数据
采集干空气条件下不同工况的热线单元功率P、热线温度Tw、气流速度Va、气流静温Ts;
S20. 计算热线单元干空气努塞尔数Nu:
其中,k为干空气热传导系数,lw为热线长度,π为圆周率;
S30. 获得热线单元干空气努塞尔数的拟合系数
S31. 计算热线雷诺数Rew和干空气普朗特数Pr:
S32. 设定热线单元干空气努塞尔数Nu的拟合关系式为:
其中,A和B为拟合系数;
S33. 将步骤S20计算的热线单元干空气努塞尔数Nu和步骤S31计算获得的Rew和Pr带入步骤S32设定的拟合关系式进行拟合,计算得到参数A=A1,B=B1,从而得到热线单元干空气努塞尔数Nu的表达式:
S30. 根据下式计算云雾场中热线单元干功率:
其中,Tw1为云雾场中的热线温度,Ts1为云雾场中的气流静温。
进一步地,步骤S31中,热线雷诺数Rew和干空气普朗特数Pr的计算公式为:
其中,ρa为气流密度,dw为热线直径,μ为干空气力学粘性系数,Cpa为干空气定压比热。
本发明还提供一种液态水含量的计算方法,其特征在于,基于前述计算的干功率来计算,采用下式计算结冰云雾液态水含量LWC:
其中,Ptotal是在云雾场中测得的热线单元功率,Ec为热线单元的液态水滴撞击系数,Rc为热线单元的液态水滴蒸发比,Le为液滴蒸发潜热,Cl为液态水比热,Te为液滴蒸发温度,SV为热线单元的采样体积。
进一步地,SV的计算式为:
其中,Va1为云雾场中的气流速度,dw为热线直径。
采用本发明的干功率计算方法以及基于此的液态水含量计算方法,相对于现有技术,至少具有以下有益效果:
本发明创造性地发现了热线单元干空气努塞尔数Nu与热线雷诺数Rew和干空气普朗特数Pr有很好的拟合关系,基于此可以利用该拟合关系式计算云雾条件下的干功率,进一步计算结冰云雾场的液态水含量。
该方法可以实现结冰风洞中结冰云雾液态水含量的精细化计算,有效地克服了现有技术测量精度较低的技术问题,为大型结冰风洞结冰云雾精细化测量提供了重要技术支撑。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面所描述的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例的一种液态水含量的计算方法流程图;
图2是本发明实施例1的干空气条件中各气流速度条件下热线单元功率测量结果;
图3是发明实施例1的干空气条件中各气流速度条件下热线单元干空气努塞尔数Nu计算结果;
图4是本发明实施例1的干空气条件中各气流速度条件下热线雷诺数Rew计算结果;
图5是本发明实施例1的干空气条件中各气流速度条件下干空气普朗特数Pr计算结果;
图6是本发明实施例1的干空气条件下Nu拟合曲线;
图7是本发明实施例1的热线单元干功率Pdry计算结果;
图8是本发明实施例2的液态水含量计算结果随时间的变化曲线。
具体实施方式
以下的说明提供了许多不同的实施例、或是例子,用来实施本发明的不同特征。以下特定例子所描述的元件和排列方式,仅用来精简的表达本发明,其仅作为例子,而并非用以限制本发明。
实施例1
本实施例1提供一种热线单元干功率计算方法,包括以下步骤:
S10. 实验采集热线单元干空气条件标定数据
采集干空气条件下不同工况的热线单元功率P、热线温度Tw、气流速度Va、气流静温Ts;本实施例中选用典型试验工况条件:气流速度Va=40m/s、50m/s、60m/s、70m/s、80m/s、90m/s、100m/s,气流静温Ts=-3℃,气流静压Ps=1atm。其中,所采集的各速度条件下热线单元功率如图2所示。
值得说明的是,干空气条件是指干空气流场中仅包含空气,不含有云雾颗粒。
S20. 计算热线单元干空气努塞尔数Nu:
其中,k为干空气热传导系数,lw为热线长度,π为圆周率;图3给出了各速度条件下热线单元空气努塞尔数计算结果。
S30. 获得热线单元干空气努塞尔数的拟合系数
S31. 根据以下公式计算热线雷诺数Rew和干空气普朗特数Pr:
其中,ρa为气流密度,dw为热线直径,μ为干空气力学粘性系数,Cpa为干空气定压比热,计算结果如图4和图5所示。
S32. 设定热线单元干空气努塞尔数Nu的拟合关系式为:
其中,A和B为拟合系数;
S33. 将步骤S20计算的热线单元干空气努塞尔数Nu和步骤S31计算获得的Rew和Pr带入步骤S32设定的拟合关系式进行拟合,本实施例中,采用最小二乘法进行拟合,拟合得到参数A=0.535,B=0.491,拟合曲线如图6所示,从而得到热线单元干空气努塞尔数Nu的表达式:
S30. 根据下式计算云雾场中热线单元的干功率:
其中,Tw1为云雾场中的热线温度,Ts1为云雾场中的气流静温。
本实施例中,选取典型云雾试验工况:气流速度Va=80m/s,气流静温Ts1=-3℃,气流静压Ps=1atm,喷雾时间t=120s。试验过程中,记录整个喷雾过程中的热线温度Tw1,使用标定的热线单元干空气努塞尔数Nu的表达式,计算整个喷雾过程中热线单元干功率Pdry,结果如图7所示。
实施例2
本实施例提供一种液态水含量的计算方法,在云雾试验过程中,采集整个喷雾过程中热线单元功率Ptotal,采用实施例1所计算的干功率计算液态水含量,如图1所示:
采用下式计算结冰云雾液态水含量LWC:
其中,Ptotal是在云雾场中测得的热线单元功率,Ec为热线单元的液态水滴撞击系数,Rc为热线单元的液态水滴蒸发比,Le为液滴蒸发潜热,Cl为液态水比热,Te为液滴蒸发温度,SV为热线单元的采样体积;
Va1为云雾场中的气流速度,dw为热线直径,液态水含量计算结果随时间的变化曲线如图8所示。
本实施例中,在结冰云雾中,热线单元为保持恒定的热线温度,需耗费一定的能量,即仪器测量得到的热线单元功率Ptotal,该功率由两部分组成,分别为补偿干空气对流换热热损耗的功率和补偿云雾蒸发损耗的功率,分别定义为干功率Pdry和湿功率Pwet,表示为:
而结冰云雾的液态水含量是由湿功率决定的:
式中,Td为液滴初始温度,一般情况下,液滴温度和颗粒温度与气流静温Ts一致,可以得到液态水含量基于干功率的计算公式为:
本发明基于创造性地发现了热线单元干空气努塞尔数Nu与热线雷诺数Rew和干空气普朗特数Pr有很好的拟合关系,可以得到干空气努塞尔数Nu的关系式,用来计算结冰云雾条件下的干功率,并进一步基于干功率来计算结冰云雾场中的液态水含量。该方法可以实现结冰风洞中结冰云雾液态水含量的精细化计算,为大型结冰风洞结冰云雾精细化测量提供了重要技术支撑。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种液态水含量的计算方法,其特征在于,包括以下步骤:
S10. 实验采集热线单元干空气条件标定数据
采集干空气条件下不同工况的热线单元功率P、热线温度Tw、气流速度Va、气流静温Ts;
S20. 计算热线单元干空气努塞尔数Nu:
其中,k为干空气热传导系数,lw为热线长度,π为圆周率;
S30. 获得热线单元干空气努塞尔数的拟合系数
S31. 计算热线雷诺数Rew和干空气普朗特数Pr:
S32. 设定热线单元干空气努塞尔数Nu的拟合关系式为:
其中,A和B为拟合系数;
S33. 将步骤S20计算的热线单元干空气努塞尔数Nu和步骤S31计算获得的Rew和Pr带入步骤S32设定的拟合关系式进行拟合,计算得到参数A=A1,B=B1,从而得到热线单元干空气努塞尔数Nu的表达式:
S40. 根据下式计算云雾场中热线单元干功率:
其中,Tw1为云雾场中的热线温度,Ts1为云雾场中的气流静温;
采用下式计算结冰云雾液态水含量LWC:
其中,Ptotal是在云雾场中测得的热线单元功率,Ec为热线单元的液态水滴撞击系数,Rc为热线单元的液态水滴蒸发比,Le为液滴蒸发潜热,Cl为液态水比热,Te为液滴蒸发温度,SV为热线单元的采样体积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111492134.6A CN113901591B (zh) | 2021-12-08 | 2021-12-08 | 热线单元干功率计算方法及液态水含量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111492134.6A CN113901591B (zh) | 2021-12-08 | 2021-12-08 | 热线单元干功率计算方法及液态水含量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113901591A CN113901591A (zh) | 2022-01-07 |
CN113901591B true CN113901591B (zh) | 2022-03-01 |
Family
ID=79025751
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111492134.6A Active CN113901591B (zh) | 2021-12-08 | 2021-12-08 | 热线单元干功率计算方法及液态水含量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113901591B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114858849B (zh) * | 2022-07-11 | 2022-10-04 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种动态冰的导热系数获取方法 |
CN116895039B (zh) * | 2023-09-11 | 2023-11-17 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种结冰云雾伪颗粒图像识别及特征参数测量方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112069462A (zh) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种大型结冰风洞高压空气资源消耗计算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2263776B (en) * | 1992-01-28 | 1995-05-17 | Endress & Hauser Ltd | Fluid mass flowmeter |
CN109815514B (zh) * | 2017-11-22 | 2021-06-08 | 广州汽车集团股份有限公司 | 散热器换热性能仿真的方法和系统 |
EP3537124B1 (de) * | 2018-03-08 | 2021-01-27 | ABB Schweiz AG | Verfahren und system zur nicht-intrusiven ermittlung einer temperatur eines durch einen leitungsabschnitt strömenden fluids |
CN108509738A (zh) * | 2018-04-09 | 2018-09-07 | 湖南大学 | 大功率海岛特种电源模块散热装置及其温度控制方法 |
CN110309591A (zh) * | 2019-07-01 | 2019-10-08 | 山东科技大学 | 一种平翅片换热器空气侧层流状态下换热与阻力计算方法 |
-
2021
- 2021-12-08 CN CN202111492134.6A patent/CN113901591B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112069462A (zh) * | 2020-08-28 | 2020-12-11 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种大型结冰风洞高压空气资源消耗计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113901591A (zh) | 2022-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113901591B (zh) | 热线单元干功率计算方法及液态水含量计算方法 | |
CN113899663B (zh) | 一种大型结冰风洞中复杂结冰云雾mvd测量方法 | |
King et al. | A hot-wire liquid water device having fully calculable response characteristics | |
Rogers et al. | A continuous-flow diffusion chamber for airborne measurements of ice nuclei | |
Kind et al. | Experimental and computational simulation of in-flight icing phenomena | |
Mazin et al. | Thermodynamics of icing cylinder for measurements of liquid water content in supercooled clouds | |
US11161629B2 (en) | System for numerical simulation and test verification of icing characteristics of an aerostat | |
CN111307406B (zh) | 一种结冰风洞液态水含量测量方法 | |
CN113899784B (zh) | 一种大型结冰风洞中复杂结冰云雾水含量测量方法 | |
CN114676560B (zh) | 基于光纤结冰探测的云雾参数计算及冰形预测方法、介质 | |
CN112678189B (zh) | 一种改进的结冰传感器安装位置确定方法 | |
Steen et al. | An assessment of the icing blade and the SEA multi-element sensor for liquid water content calibration of the NASA GRC Icing Research Tunnel | |
CN112985347A (zh) | 一种飞机结冰表面粗糙度计算方法 | |
Haman et al. | A new ultrafast thermometer for airborne measurements in clouds | |
Orchard et al. | Measurement of liquid water content for supercooled large drop conditions in the NRC’s altitude icing wind tunnel | |
Guo et al. | Liquid water content measurement with SEA multi-element sensor in CARDC icing wind tunnel: Calibration and performance | |
Fraser et al. | Thermodynamic limitations of ice accretion instruments | |
Feller | Investigation of Equilibrium Temperatures and Average Laminar Heat-Transfer Coefficients for the Front Half of Swept Circular Cylinders at a Mach Number of 6.9 | |
Lilie et al. | Description and results for a simple ice crystal detection system for airborne applications | |
Neel et al. | The calculated and measured performance characteristics of a heated-wire liquid-water-content meter for measuring icing severity | |
Juttula et al. | Study of the aerodynamic sampling effects of a holographic cloud droplet instrument | |
Clark et al. | A Guide Creating SAE AS5562 Ice Crystal, Mixed Phase and Rain Conditions in a Wind Tunnel Environment | |
Zhivoglotov | Estimation of liquid water content effects on the air temperature measurements in the clouds based on the wind tunnel experiments | |
Bhalwankar et al. | The evaporation of the charged and uncharged water drops suspended in a wind tunnel | |
CN115649455B (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 |