CN112197721B - 一种水膜厚度测量方法 - Google Patents

一种水膜厚度测量方法 Download PDF

Info

Publication number
CN112197721B
CN112197721B CN202011413320.1A CN202011413320A CN112197721B CN 112197721 B CN112197721 B CN 112197721B CN 202011413320 A CN202011413320 A CN 202011413320A CN 112197721 B CN112197721 B CN 112197721B
Authority
CN
China
Prior art keywords
water
thickness
test
water film
ultrasonic probe
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
Application number
CN202011413320.1A
Other languages
English (en)
Other versions
CN112197721A (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.)
Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Low Speed Aerodynamics 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 Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202011413320.1A priority Critical patent/CN112197721B/zh
Publication of CN112197721A publication Critical patent/CN112197721A/zh
Application granted granted Critical
Publication of CN112197721B publication Critical patent/CN112197721B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
    • G01B17/02—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring thickness

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明适用于风洞测量技术领域，提供了一种水膜厚度测量方法，其包括如下步骤：将测量装置布置于结冰风洞试验段内；将结冰风洞试验段内的气流速度和气体温度调节到预设值，通过冷却盘管将试验平板的温度调节到预设值；活塞不断上移使储水容器内的水上升至排水管的出口；测量超声探头对应位置处的水膜厚度和结冰厚度；测量出试验平板的倾斜角度；修正测量的水膜厚度和结冰厚度。本发明在测量水膜厚度时，并没有将结冰层的厚度计入水膜厚度，且通过试验平板的倾斜角对水膜厚度进行了修正，从而使得水膜厚度的测量更加精确。

Description

一种水膜厚度测量方法
技术领域
本发明属于风洞测量技术领域，尤其涉及一种水膜厚度测量方法。
背景技术
当飞机穿越含有过冷水滴的云层时，其迎风部件表面会发生结冰现象，进而严重威胁飞行安全。试验研究发现，当环境温度比较高时，撞击在结冰表面的过冷水滴不会完全冻结，而会在剪切力、重力及表面张力的作用下向后流动，在机翼结冰表面与空气之间形成一层流动的水膜。这层水膜具有流动速度小、厚度薄且厚度随冻结和蒸发作用而不断变化等特点。水膜流动冻结是飞机结冰及防除冰过程中广泛存在的物理现象。
水膜在冷板表面流动时，在空间上是沿液固接触面向上推移、沿水膜流动方向向后推移，在时间上是伴冻结、再辉、冷却的非平衡相变过程。然而现有对于水膜流动的研究过于简化，如将水膜视为剪切力驱动的Couette流，忽略了水膜流动与冻结的耦合作用，使得基于水膜模型的计算结果与实验结果仍存在一定偏差。
发明内容
本发明的目的在于提供一种水膜厚度测量方法，旨在解决现有技术中的水膜测量精度不高的技术问题。
本发明提供了一种水膜厚度测量方法，其包括如下步骤：
步骤S10:将测量装置布置于结冰风洞试验段内，所述测量装置包括试验组件、供水组件、积水组件，所述供水组件包括储水容器、活塞、排水管，所述活塞可沿储水容器滑动，所述排水管与储水容器的出口连接；所述试验组件包括试验平板、超声探头组件、冷却盘管、保温材料，所述超声探头组件、冷却盘管位于试验平板和所述保温材料之间；所述积水组件包括排水箱和位于排水箱内的吸水材料；
步骤S20:将结冰风洞试验段内的气流速度和气体温度调节到预设值，通过冷却盘管将试验平板的温度调节到预设值；其中，所述试验平板的长度为l，宽度为w，所述超声探头组件由m×n个超声探头组成，其中，m为超声探头的总行数，n为超声探头的总列数，超声探头表示为i x,y ，其中，x为超声探头的行号，y为超声探头的列号，1≤x≤m，1≤y≤n；
步骤S30:活塞不断上移使储水容器内的水上升至排水管的出口；
步骤S40：测量m×n个超声探头对应位置处的水膜厚度和结冰厚度，其中，超声探头i x,y 对应位置的水膜厚度表示为H water,x,y 和超声探头i x,y 对应位置的结冰厚度表示为H ice,x,y ；
步骤S50：测量出试验平板的倾斜角度α；
步骤S60：根据试验平板的倾斜角度α，修正步骤S40中测量的水膜厚度和结冰厚 度，其中，水膜厚度H water,x,y 修正后表示为
Figure 765867DEST_PATH_IMAGE001
，结冰厚度H ice,x,y 修正后表示为
Figure 135669DEST_PATH_IMAGE002
，其中，
Figure 337849DEST_PATH_IMAGE003
。
进一步地，所述步骤S50包括如下步骤：
步骤S51：计算试验平板上结冰的总体积V 1 ；
步骤S52：计算所述超声探头组件探测到的结冰的总体积V 2 ；
步骤S53：计算倾斜角度α，其中
Figure 134903DEST_PATH_IMAGE004
，其中，H AB 通过以下方程联 立计算得到：
Figure 743739DEST_PATH_IMAGE005
，
Figure 979548DEST_PATH_IMAGE006
表示 位于试验平板的起始端下方的超声探头测得的结冰厚度。
进一步地，所述步骤S51包括如下步骤：
步骤S511：计算流向试验平板的水的总质量
Figure 990230DEST_PATH_IMAGE007
，计算流出试验平板的水 的总质量
Figure 771235DEST_PATH_IMAGE008
，计算从试验平板表面蒸发的水的质量
Figure 867367DEST_PATH_IMAGE009
；
步骤S512：计算试验平板上结冰的总质量
Figure 641288DEST_PATH_IMAGE010
，其中
Figure 506476DEST_PATH_IMAGE011
，并根据
Figure 645333DEST_PATH_IMAGE012
计算V 1 ，其中，
Figure 540346DEST_PATH_IMAGE013
，
Figure 852378DEST_PATH_IMAGE014
表示结冰密度。
进一步地，流向试验平板的水的总质量
Figure 650701DEST_PATH_IMAGE015
通过以下公式计算得到：
Figure 22777DEST_PATH_IMAGE016
，其中，p为活塞上移的行程，A为活塞的截面积，
Figure 405085DEST_PATH_IMAGE017
为水密度。
进一步地，从试验平板表面蒸发的水的质量
Figure 458492DEST_PATH_IMAGE018
通过以下公式计算得到：
Figure 298272DEST_PATH_IMAGE019
，
其中，
Figure 841249DEST_PATH_IMAGE020
为水膜温度，
Figure 133690DEST_PATH_IMAGE021
为气流温度；
Figure 803837DEST_PATH_IMAGE022
为水膜表面的水蒸气饱和蒸汽压，
Figure 498123DEST_PATH_IMAGE023
为气流的水蒸气饱和蒸汽压；
Figure 212001DEST_PATH_IMAGE024
为空气密度，
Figure 257318DEST_PATH_IMAGE025
为空气定压比热，
Figure 964112DEST_PATH_IMAGE026
为气体 常数；h为对流换热系数。
进一步地，超声探头组件探测到的结冰的总体积V 2 通过以下公式计算得到：
Figure 512905DEST_PATH_IMAGE027
。
本发明相对于现有技术至少具有如下技术效果：
1.本发明发现了造成水膜厚度测量不精确的原因在于：一是把水膜之下结冰层计算入水膜厚度；二是没有考虑试验平板的倾斜；
2.本发明在测量水膜厚度时，并没有将结冰层的厚度计入水膜厚度，且通过试验平板的倾斜角对水膜厚度进行了修正，从而使得水膜厚度的测量更加精确；
3.本发明在测量试验平板的倾斜角时，通过已有的设备即可完成，无需单独的测角仪，一方面，节省了测量设备和测量时间；另一方面，也避免了使用测角仪造成的累计误差。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面所描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种水膜厚度测量装置示意图；
图2是本发明实施例提供的试验平板的倾斜的示意图；
图3是本发明实施例提供的超声探头的分布；
图4是本发明实施例提供的一种水膜厚度测量方法的示意图；
图5是本发明实施例提供的试验平板的倾斜角度测量原理图；
图6是本发明实施例提供的倾斜角度的测量方法示意图。
具体实施方式
在下文中将参考附图对本发明的各方面进行更充分的描述。然而，本发明可以具体化成许多不同形式且不应解释为局限于贯穿本发明所呈现的任何特定结构或功能。相反地，提供这些方面将使得本发明周全且完整，并且本发明将给本领域技术人员充分地传达本发明的范围。基于本文所教导的内容，本领域的技术人员应意识到，无论是单独还是结合本发明的任何其它方面实现本文所公开的任何方面，本发明的范围旨在涵盖本文中所公开的任何方面。例如，可以使用本文所提出任意数量的装置或者执行方法来实现。另外，除了本文所提出本发明的多个方面之外，本发明的范围更旨在涵盖使用其它结构、功能或结构和功能来实现的装置或方法。应可理解，其可通过权利要求的一或多个元件具体化本文所公开的任何方面。
在此使用的术语仅仅是为了描述具体实施例，而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在，但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。
在此使用的所有术语(包括技术和科学术语)具有本领域技术人员通常所理解的含义，除非另外定义。应注意，这里使用的术语应解释为具有与本说明书的上下文相一致的含义，而不应以理想化或过于刻板的方式来解释。
对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。
本发明实施例中的水膜厚度测量方法，将采用如图1所示的水膜厚度测量装置，该测量装置布置于结冰风洞试验段1内，所述测量装置包括试验组件、供水组件、积水组件，所述供水组件包括储水容器2、活塞3、排水管4，所述活塞3可沿储水容器2滑动，所述排水管4与储水容器2的出口连接；所述试验组件包括试验平板7、超声探头组件8、冷却盘管9、保温材料10，所述超声探头组件8、冷却盘管9位于试验平板7和所述保温材料10之间；所述积水组件包括排水箱11和位于排水箱内的吸水材料12；
在实际进行测量时，结冰风洞试验段1为测量装置提供气流，储水容器2用于盛装水，活塞可以上下移动，当活塞上移到一定程度时，储水容器2内的水通过排水管4不断溢出，在气流的作用下，溢出的水将不断在试验平板7上铺开；另外，在正式试验前，通常需要通过冷却盘管9将试验平板7的温度降低到一定程度，因此，溢出的水在试验平板7上铺开时，将有结冰现象的存在；对于未结冰的水，将不断朝向位于试验平板7尾端的积水组件流动，吸水材料12将流入的水吸收，同时，积水组件中带有称重模块（未图示），该称重模块用于测量流入积水组件的水的质量。
进一步地，活塞3与曲柄连杆5连接，通过电机驱动曲柄连杆5运动，从而带动活塞上下移动，由于电机的转速可通过信号发生器进行控制，调节信号发生器的输出频率即可控制活塞和储水容器2内水的上移速度，实现精确供水控制。同时，在储水容器2上还设置有供水箱6，用于向储水容器2供水，优选地，在供水箱6和储水容器2之间还设置有供水阀15，以方便控制。
进一步地，在排水箱11的出口还设置有排水阀14，以方便排水。
现有技术中，往往直接测量水膜的厚度，但是将测量出的水膜厚度和计算机模拟出的水膜厚度相比，往往存在一定的差异；发明人通过观察发现，试验平板上除了形成水膜外，在水膜的下方还形成了结冰层，而现有技术中，将结冰层的厚度同样算入了水膜厚度内，因而造成了水膜厚度测量的不精确；值得说明的是，这一造成水膜厚度测量不精确的原因对于本领域技术人员而言并不是容易发现的；
另一方面，造成水膜厚度测量不精确的另一原因在于安装试验平板时，不可能达到绝对的水平，因此，对于水膜厚度的测量造成了误差；同样地，这一造成水膜厚度测量不精确的原因对于本领域技术人员而言也并不是容易发现的；
总之，发明人发现，造成水膜厚度测量不精确的原因在于：一是把水膜之下结冰层计算入水膜厚度；二是没有考虑试验平板的倾斜；
对于试验平板的倾斜造成的误差，以图2所示进行进一步说明：
aa 1 为超声探头测量的a点的冰层厚度值，a 1 a 2 为超声探头测量的a点的水膜厚度值，由于试验平板存在倾斜角α，实际上aa 3 才是a点的冰层的实际厚度值，a 3 a 4 才是a点的水膜的实际厚度值。其中，a 1 为过a点的竖直线与冰层表面的交点，a 2 为过a点的竖直线与水膜表面的交点，a 3 为过a点且垂直于试验平板7的辅助线与冰层表面的交点，a 4 为过a点且垂直于试验平板7的辅助线与水膜表面的交点。
如图3所示，超声探头组件8由多个超声探头组成，具体地，所述超声探头组件8由m ×n个超声探头组成，其中，m为超声探头的总行数，n为超声探头的总列数，超声探头表示为i x,y ，其中，x为超声探头的行号，y为超声探头的列号，1≤x≤m，1≤y≤n；其中，第一行的超声探头安装得尽量靠近试验平板7的第一边，试验平板7的第一边指的是靠近排水管4的一端；第一列的超声探头安装得尽量靠近与试验平板7的第一边垂直相交的第二边；最后一行的超声探头安装得尽量靠近试验平板7的第三边，第三边为与第一边相对的边；最后一列的超声探头安装得尽量靠近试验平板7的第四边，第四边为与第二边相对的边。
超声探头的行间距相等，超声探头的列间距相等；具体地，试验平板7的长度为l，宽度为w，则超声探头的行间距为l/(m-1)，超声探头的列间距为w/(n-1)。
如图4所示为本发明实施例中的一种水膜厚度测量方法的示意图，其包括如下步骤：
步骤S10:将测量装置布置于结冰风洞试验段1内，该测量装置包括试验组件、供水组件、积水组件，所述供水组件包括储水容器2、活塞3、排水管4，所述活塞3可沿储水容器2滑动，所述排水管4与储水容器2的出口连接；所述试验组件包括试验平板7、超声探头组件8、冷却盘管9、保温材料10，所述超声探头组件8、冷却盘管9位于试验平板7和所述保温材料10之间；所述积水组件包括排水箱11和位于排水箱内的吸水材料12；
步骤S20:将结冰风洞试验段内的气流速度和气体温度调节到预设值，通过冷却盘管9将试验平板7的温度调节到预设值；其中，所述试验平板7的长度为l，宽度为w，所述超声探头组件8由m×n个超声探头组成，其中，m为超声探头的总行数，n为超声探头的总列数，超声探头表示为i x,y ，其中，x为超声探头的行号，y为超声探头的列号，1≤x≤m，1≤y≤n；
步骤S30:活塞3不断上移使储水容器2内的水上升至排水管4的出口；
步骤S40：测量m×n个超声探头对应位置处的水膜厚度和结冰厚度，其中，超声探头i x,y 对应位置的水膜厚度表示为H water,x,y 和超声探头i x,y 对应位置的结冰厚度表示为H ice,x,y ；结合附图2，则aa 1 =H ice,x,y ，a 1 a 2 =H water,x,y ；
步骤S50：测量出试验平板的倾斜角度α；
步骤S60：根据试验平板的倾斜角度α，修正步骤S40中测量的水膜厚度和结冰厚 度，其中，水膜厚度H water,x,y 修正后表示为
Figure 663263DEST_PATH_IMAGE001
，结冰厚度H ice,x,y 修正后表示为
Figure 930296DEST_PATH_IMAGE028
，其中，
Figure 129197DEST_PATH_IMAGE029
。
结合附图2，则aa 3 =
Figure 345545DEST_PATH_IMAGE028
，a 3 a 4 =
Figure 604488DEST_PATH_IMAGE001
；
本发明实施例中，一方面，在测量水膜厚度时，考虑了水膜下的结冰，因此，测量的水膜厚度值相对于现有技术而言，更加精确；另一方面，本发明实施例中还测量了试验平板的倾斜角度α，并用以修正测量的水膜厚度和结冰厚度，本发明实施例的测量精度进一步得到提高。
以下将重点论述，本发明实施例中的有关倾斜角度α的测量。
本发明并没有采用直接测量的方式来测量试验平板的倾斜角度，而是通过了如图5所示的原理进行测量：
首先计算在试验平板7上结冰的总体积V 1 ，再计算超声探头组件8探测到的结冰的总体积V 2 ，再通过V 1 和V 2 计算如图5中的以ΔABC为底的柱形体的体积V 3 ，利用V 3 便可计算出图5中的ΔABC中的∠A,∠A即为倾斜角度α。其中，AB是试验平板7为倾斜状态下的垂直线，AC是试验平台7为水平状态的垂直线，B点和C点分别是冰层厚度与AB和AC的交点。
具体地，如图6所示，倾斜角度α的测量包括如下步骤：
步骤S51：计算试验平板7上结冰的总体积V 1 ；
步骤S52：计算所述超声探头组件8探测到的结冰的总体积V 2 ；
步骤S53：计算倾斜角度α，其中
Figure 421135DEST_PATH_IMAGE030
，其中，H AB 通过以下方程 联立计算得到：
Figure 158147DEST_PATH_IMAGE031
，
Figure 993116DEST_PATH_IMAGE032
表 示行号和列号均为1的超声探头测得的结冰厚度；
其中，AB是试验平板7为倾斜状态下的垂直线，AC是试验平板7为水平状态的垂直线，B点和C点分别是冰层厚度与AB和AC的交点，则H AB 为冰层与倾斜状态下试验平板垂直线相交的高度，为实际冰层厚度；H BC 为BC段的长度。
进一步地，所述步骤S51包括如下步骤：
步骤S511：计算流向试验平板的水的总质量
Figure 422961DEST_PATH_IMAGE007
，计算流出试验平板的水 的总质量
Figure 461324DEST_PATH_IMAGE008
，计算从试验平板表面蒸发的水的质量
Figure 2027DEST_PATH_IMAGE009
；
步骤S512：计算试验平板7上结冰的总质量
Figure 458547DEST_PATH_IMAGE010
，其中
Figure 793713DEST_PATH_IMAGE033
，根据
Figure 522635DEST_PATH_IMAGE012
计算V 1 ，
Figure 663766DEST_PATH_IMAGE034
，
Figure 161744DEST_PATH_IMAGE035
表示结冰密度。
进一步地，流向试验平板的水的总质量
Figure 979396DEST_PATH_IMAGE036
通过以下公式计算得到：
Figure 195613DEST_PATH_IMAGE037
，其中，p为活塞上移的行程，A为活塞的截面积，
Figure 140436DEST_PATH_IMAGE038
为水密度。
另外，当通过电机驱动曲柄连杆5运动，从而带动活塞上下移动时，p=N×S，其中，N为电机转动的次数；S为电机每转动一周活塞移动的距离。
进一步，从试验平板表面蒸发的水的质量
Figure 492920DEST_PATH_IMAGE039
通过以下公式计算得到：
Figure 248517DEST_PATH_IMAGE040
，
其中，
Figure 686452DEST_PATH_IMAGE020
为水膜温度，
Figure 107069DEST_PATH_IMAGE021
为气流温度；
Figure 641955DEST_PATH_IMAGE022
为水膜表面的水蒸气饱和蒸汽 压，
Figure 489826DEST_PATH_IMAGE023
为气流的水蒸气饱和蒸汽压；
Figure 726641DEST_PATH_IMAGE024
为空气密度，
Figure 950949DEST_PATH_IMAGE025
为空气定压比热，
Figure 278025DEST_PATH_IMAGE026
为气体常数；h为对流换热系数。
进一步地，关于超声探头组件8探测到的结冰的总体积V 2 则采用了离散化的思想，具体地，以每个超声探头为基准，将V 2 离散成(m-1)×(n-1)个长方体，首先计算每个长方体的体积，再累计求和即可，具体地，超声探头组件8探测到的结冰的总体积V 2 通过以下公式计算得到：
Figure 624693DEST_PATH_IMAGE027
。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种水膜厚度测量方法，其特征在于，包括如下步骤：
步骤S10:将测量装置布置于结冰风洞试验段（1）内，所述测量装置包括试验组件、供水组件、积水组件，所述供水组件包括储水容器（2）、活塞（3）、排水管（4），所述活塞（3）可沿储水容器（2）滑动，所述排水管（4）与储水容器（2）的出口连接；所述试验组件包括试验平板（7）、超声探头组件（8）、冷却盘管（9）、保温材料（10），所述超声探头组件（8）、冷却盘管（9）位于试验平板（7）和所述保温材料（10）之间；所述积水组件包括排水箱（11）和位于排水箱内的吸水材料（12）；
步骤S20:将结冰风洞试验段内的气流速度和气体温度调节到预设值，通过冷却盘管（9）将试验平板（7）的温度调节到预设值；其中，所述试验平板（7）的长度为l，宽度为w，所述超声探头组件（8）由m×n个超声探头组成，其中，m为超声探头的总行数，n为超声探头的总列数，超声探头表示为i x,y ，其中，x为超声探头的行号，y为超声探头的列号，1≤x≤m，1≤y ≤n；
步骤S30:活塞（3）不断上移使储水容器（2）内的水上升至排水管（4）的出口；
步骤S40：测量m×n个超声探头对应位置处的水膜厚度和结冰厚度，其中，超声探头i x,y 对应位置的水膜厚度表示为H water,x,y 和超声探头i x,y 对应位置的结冰厚度表示为H ice,x,y ；
步骤S50：测量出试验平板的倾斜角度α；
步骤S60：根据试验平板的倾斜角度α，修正步骤S40中测量的水膜厚度和结冰厚度，其 中，水膜厚度H water,x,y 修正后表示为
Figure 527759DEST_PATH_IMAGE001
，结冰厚度H ice,x,y 修正后表示为
Figure 290179DEST_PATH_IMAGE002
，其 中，
Figure 845663DEST_PATH_IMAGE003
。
2.如权利要求1所述的一种水膜厚度测量方法，其特征在于，所述步骤S50包括如下步骤：
步骤S51：计算试验平板（7）上结冰的总体积V 1 ；
步骤S52：计算所述超声探头组件（8）探测到的结冰的总体积V 2 ；
步骤S53：计算倾斜角度α，其中
Figure 300915DEST_PATH_IMAGE004
，其中，H AB 通过以下方程联立计算得 到：
Figure 200738DEST_PATH_IMAGE005
，
Figure 829165DEST_PATH_IMAGE006
表示位于试验平板的起始端下方的 超声探头测得的结冰厚度；
其中，AB是试验平板（7）为倾斜状态下的垂直线，AC是试验平板（7）为水平状态的垂直线，B点和C点分别是冰层厚度与AB和AC的交点，则H AB 为冰层与倾斜状态下试验平板垂直线相交的高度，为实际冰层厚度；H BC 为BC段的长度。
3.如权利要求2所述的一种水膜厚度测量方法，其特征在于，所述步骤S51包括如下步骤：
步骤S511：计算流向试验平板的水的总质量
Figure 865254DEST_PATH_IMAGE007
，计算流出试验平板的水的总质量
Figure 304457DEST_PATH_IMAGE008
，计算从试验平板表面蒸发的水的质量
Figure 691576DEST_PATH_IMAGE009
；
步骤S512：计算试验平板（7）上结冰的总质量
Figure 795798DEST_PATH_IMAGE010
，其中
Figure 14290DEST_PATH_IMAGE011
，并根据
Figure 811345DEST_PATH_IMAGE012
计算V 1 ，其中，
Figure 731765DEST_PATH_IMAGE013
，
Figure 905257DEST_PATH_IMAGE014
表 示结冰密度。
4.如权利要求3所述的一种水膜厚度测量方法，其特征在于，流向试验平板的水的总质 量
Figure 915939DEST_PATH_IMAGE015
通过以下公式计算得到：
Figure 946211DEST_PATH_IMAGE016
，其中，p为活塞上移的 行程，A为活塞的截面积，
Figure 42343DEST_PATH_IMAGE017
为水密度。
5.如权利要求3所述的一种水膜厚度测量方法，其特征在于，从试验平板表面蒸发的水 的质量
Figure 301418DEST_PATH_IMAGE018
通过以下公式计算得到：
Figure 432185DEST_PATH_IMAGE019
，
其中，
Figure 633359DEST_PATH_IMAGE020
为水膜温度，
Figure 216787DEST_PATH_IMAGE021
为气流温度；
Figure 778087DEST_PATH_IMAGE022
为水膜表面的水蒸气饱和蒸汽压，
Figure 763361DEST_PATH_IMAGE023
为气流的水蒸气饱和蒸汽压；
Figure 73119DEST_PATH_IMAGE024
为空气密度，
Figure 206160DEST_PATH_IMAGE025
为空气定压比热，
Figure 259567DEST_PATH_IMAGE026
为气体常 数；h为对流换热系数。
6.如权利要求3所述的一种水膜厚度测量方法，其特征在于，超声探头组件（8）探测到的结冰的总体积V 2 通过以下公式计算得到：
Figure 99347DEST_PATH_IMAGE027
。
CN202011413320.1A 2020-12-07 2020-12-07 一种水膜厚度测量方法 Active CN112197721B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011413320.1A CN112197721B (zh) 2020-12-07 2020-12-07 一种水膜厚度测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011413320.1A CN112197721B (zh) 2020-12-07 2020-12-07 一种水膜厚度测量方法

Publications (2)

Publication Number Publication Date
CN112197721A CN112197721A (zh) 2021-01-08
CN112197721B true CN112197721B (zh) 2021-02-19

Family

ID=74033755

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011413320.1A Active CN112197721B (zh) 2020-12-07 2020-12-07 一种水膜厚度测量方法

Country Status (1)

Country Link
CN (1) CN112197721B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113237796A (zh) * 2021-05-13 2021-08-10 北京建筑大学 基于蒸发介质润湿度的热湿传递特性测试系统及方法
CN114180072B (zh) * 2022-02-16 2022-04-12 中国空气动力研究与发展中心低速空气动力研究所 一种结冰厚度探测方法
CN114877837B (zh) * 2022-07-11 2022-09-09 中国空气动力研究与发展中心低速空气动力研究所 一种基于远场声压级变化的结冰厚度测量方法及系统
CN117302522B (zh) * 2023-11-28 2024-02-09 中国空气动力研究与发展中心低速空气动力研究所 一种用于飞行设备的低功耗超声波防除冰方法及装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2872897B1 (fr) * 2004-07-06 2006-10-13 Commissariat Energie Atomique Dispositif optique de mesure de l'epaisseur d'un milieu au moins partiellement transparent
RU2352493C1 (ru) * 2007-07-16 2009-04-20 Федеральное государственное унитарное предприятие "Центральный научно-исследовательский институт имени академика А.Н. Крылова" (ФГУП "ЦНИИ им. акад. А.Н. Крылова") Опытовый бассейн для испытаний моделей судов преимущественно во льдах
CN104897358B (zh) * 2015-06-17 2017-08-08 北京航空航天大学 一种应用于水膜发生与测量的实验装置
CN105222740B (zh) * 2015-09-24 2018-09-04 周志宏 一种多传感器联合测量冰厚的方法
CN106124562A (zh) * 2016-06-22 2016-11-16 中国航空工业集团公司沈阳发动机设计研究所 一种液态水含量测量方法
CN107021184B (zh) * 2017-03-22 2019-01-15 哈尔滨工程大学 一种船舶实水域破冰试验方法
CN110006620A (zh) * 2019-02-18 2019-07-12 中国空气动力研究与发展中心低速空气动力研究所 一种水膜测量系统
CN110455237B (zh) * 2019-08-09 2021-01-26 水利部南京水利水文自动化研究所 一种冰层厚度自动检测装置及方法
CN111307406B (zh) * 2020-05-06 2020-11-06 中国空气动力研究与发展中心低速空气动力研究所 一种结冰风洞液态水含量测量方法
CN111323200B (zh) * 2020-05-11 2020-08-07 中国空气动力研究与发展中心低速空气动力研究所 一种结冰风洞试验冰形面积计算方法

Also Published As

Publication number Publication date
CN112197721A (zh) 2021-01-08

Similar Documents

Publication Publication Date Title
CN112197721B (zh) 一种水膜厚度测量方法
CN103776763B (zh) 旋转圆柱结冰与冰层粘附力测量的实验装置及测量方法
Struk et al. An initial study of the fundamentals of ice crystal icing physics in the NASA propulsion systems laboratory
CN113899663B (zh) 一种大型结冰风洞中复杂结冰云雾mvd测量方法
CN114676560B (zh) 基于光纤结冰探测的云雾参数计算及冰形预测方法、介质
CN111307406B (zh) 一种结冰风洞液态水含量测量方法
US11781942B2 (en) Method for determining flow velocity distribution in roughness sublayers
Dehghani-sanij et al. Analysis of ice accretion on vertical surfaces of marine vessels and structures in arctic conditions
He et al. Wind tunnel test on the flow resistance of U-type water collecting devices for natural draft wet cooling towers
CN204613140U (zh) 用于云雾液态水含量的测量装置
Gao et al. Measurement of unsteady transition on a pitching airfoil using dynamic pressure sensors
CN113899784A (zh) 一种大型结冰风洞中复杂结冰云雾水含量测量方法
Weske Experimental investigation of velocity distributions downstream of single duct bends
Zhou et al. Experimental study on liquid drainage characteristics of metal foams under sloshing conditions
Davison et al. Naturally aspirating isokinetic total water content probe: Wind tunnel test results and design modifications
McClain et al. A climatic facility and apparatus for investigations of cold soaked fuel frost evolution
Dehkordi et al. The effect of droplet size and liquid water content on ice accretion and aerodynamic coefficients of tower legs
CN216524683U (zh) 一种用于结冰风洞试验中液态水含量的测量装置
Lilie et al. Description and results for a simple ice crystal detection system for airborne applications
Addy et al. A study of the effects of altitude on thermal ice protection system performance
CN109060295A (zh) 多功能水膜发生实验装置及方法
RU157695U1 (ru) Испытательный стенд по определению степени гидрофобности материалов для изготовления труб и ремонта трубопроводов
CN106645260A (zh) 一种用于土体有压冻融试验的自动供排水系统
CN108169819B (zh) 一种基于谐振式结冰探测器的过冷水含量测量方法
Esposito et al. Cloud calibration update of the CIRA icing 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