CN113252281B - 一种结冰云雾液滴尺寸分布的重构方法 - Google Patents

一种结冰云雾液滴尺寸分布的重构方法 Download PDF

Info

Publication number
CN113252281B
CN113252281B CN202110611640.6A CN202110611640A CN113252281B CN 113252281 B CN113252281 B CN 113252281B CN 202110611640 A CN202110611640 A CN 202110611640A CN 113252281 B CN113252281 B CN 113252281B
Authority
CN
China
Prior art keywords
droplet
diameter
measurement interval
drop
size
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
CN202110611640.6A
Other languages
English (en)
Other versions
CN113252281A (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 CN202110611640.6A priority Critical patent/CN113252281B/zh
Publication of CN113252281A publication Critical patent/CN113252281A/zh
Application granted granted Critical
Publication of CN113252281B publication Critical patent/CN113252281B/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/02—Wind tunnels
    • 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
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明适用于风洞试验技术领域，提供了一种结冰云雾液滴尺寸分布的重构方法，其包括如下步骤：将液滴尺寸测量区间划分为多个液滴尺寸通道，计算各个液滴尺寸通道内的液滴数密度、第一液态水含量、标准液态水含量；确定液滴直径测量区间内的截断液滴尺寸通道，确定搭接液滴尺寸通道，并计算搭接液滴尺寸通道内的液态水含量；重构结冰云雾液滴尺寸分布。本发明解决了现有技术中利用单一仪器只能进行较小的液滴尺寸范围的测量的技术问题，能够从不同仪器的测量结果中得到全液滴谱液滴尺寸分布及特征参数，也就是实现了结冰云雾的液滴的全尺寸测量和对结冰云雾的参数化表达。

Description

一种结冰云雾液滴尺寸分布的重构方法
技术领域
本发明属于风洞试验技术领域，尤其涉及一种结冰云雾液滴尺寸分布的重构方法。
背景技术
当飞机在云层中飞行时，云层中的过冷水滴（即温度低于冰点的液态水滴）会不断撞击飞机迎风面，导致其表面发生结冰现象。飞机结冰广泛存在于飞行实践中，并严重威胁飞行安全。
结冰风洞是开展飞机结冰研究及飞机部件防除冰系统验证的重要地面试验设备，其在飞机结冰适航审定中扮演着重要角色。结冰云雾环境模拟能力是结冰风洞性能的核心内容。结冰云雾微物理特征的准确测量与评估是结冰风洞准确模拟结冰云雾环境的关键，其中云雾液滴尺寸分布及特征参数是结冰云雾重要的微物理特征之一。液滴尺寸分布定义为云雾内各液滴直径对应的液滴浓度分布，液滴浓度通常包括液滴数量、数密度、液态水含量等参数。另外，体积中值直径是结冰云雾重要的微物理特征参数，体积中值直径定义为累积体积分数为0.5对应的液滴直径，即大于该直径的液滴体积与小于该直径的液滴体积相等，该参数基于液滴尺寸分布计算得到。
目前，国内结冰风洞通常采用单一原理的测量仪器测量云雾的液滴尺寸分布特征，然而，针对结冰风洞模拟的结冰云雾，尤其对于过冷大水滴结冰云雾，液滴直径范围覆盖微米尺度到毫米尺度，跨越三个量级，因此，单一仪器只能进行较小的液滴尺寸范围的测量，无法实现结冰云雾的液滴的全尺寸测量，也就无法对结冰云雾的参数化表达。
发明内容
本发明的目的在于提供一种结冰云雾液滴尺寸分布的重构方法，旨在实现结冰云雾的液滴的全尺寸测量。
本发明提供了一种结冰云雾液滴尺寸分布的重构方法，其包括如下步骤：
步骤S10：将液滴尺寸测量区间划分为多个液滴尺寸通道，计算各个液滴尺寸通道内的液滴数密度、第一液态水含量、标准液态水含量。
步骤S20：确定液滴直径测量区间内的截断液滴尺寸通道，确定搭接液滴尺寸通道，并计算搭接液滴尺寸通道内的液态水含量；
步骤S30：重构结冰云雾液滴尺寸分布。
进一步地，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道记为B i,j ，其中，B i,j =[D i,j-1 ，D i,j ]，i从1～L遍历，j从1～M i 遍历，i为液滴直径测量区间的序号，j为液滴尺寸通道的序号，D i,j-1 为第i个液滴直径测量区间内第j个液滴尺寸通道的最小液滴直径，D i,j 为第i个液滴直径测量区间内第j个液滴尺寸通道的最大液滴直径，L为液滴直径测量区间的个数，M i 为第i个液滴直径测量区间内的液滴尺寸通道的个数。
进一步地，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道内的液滴数密度记为C i,j ，其中，C i,j =N i,j /SV i,j ，i从1～L遍历，j从1～M i 遍历，N i,j 为第i个液滴直径测量区间内的第j个液滴尺寸通道内的液滴数量，SV i,j 为第i个液滴直径测量区间内的第j个液滴尺寸通道内的采样体积。
进一步地，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道内的第一液态水含量记为LWC i,j ，第i个液滴直径测量区间内的第j个液滴尺寸通道内的标准液态水含量记为NLWC i,j ，其中：
Figure 261868DEST_PATH_IMAGE001
，
Figure 788534DEST_PATH_IMAGE002
，
Figure 612133DEST_PATH_IMAGE003
，i从1～L遍历，j从1～M i 遍历，
Figure 265968DEST_PATH_IMAGE004
为第i个液滴直径测量区间内的第j个液滴尺寸通道内的平均液滴直径，ρ l 为液滴密度。
进一步地，所述步骤S20中，第i+1个液滴直径测量区间内的截断液滴尺寸通道记为
Figure 960255DEST_PATH_IMAGE005
，其中，
Figure 893707DEST_PATH_IMAGE006
，i从1～（L-1）遍历，k i 为第i个重叠区域对应的第i+1个液滴直径测量区间内的液滴尺寸通道的序号，k i 从1～M i+1 遍历。
进一步地，所述步骤S20中，第i个搭接液滴尺寸通道记为
Figure 939023DEST_PATH_IMAGE007
，第i个搭接液滴尺寸通道内的液态水含量记为
Figure 130970DEST_PATH_IMAGE008
，其中，
Figure 476501DEST_PATH_IMAGE009
，
Figure 564543DEST_PATH_IMAGE010
，i从1～（L-1）遍历，k i 从1～M i+1 遍历。
进一步地，所述步骤S30包括如下步骤：
步骤S31：组合重构全液滴谱液滴直径测量区间；
步骤S32：计算全液滴谱液滴直径测量区间内各通道内的液态水含量、液态水含量分数；
步骤S33：计算全液滴谱液滴直径测量区间的累积质量分数和体积中值直径。
进一步地，所述步骤S31中，重构后的全液滴谱液滴直径测量区间的液滴尺寸通道记为Bin s ，其中：
Figure 612002DEST_PATH_IMAGE011
s为重构后的全液滴谱液滴直径测量区间的液滴尺寸通道序号，D s-1为第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最小值，D s 为第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最大值。
进一步地，所述步骤S32中，第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道内的液态水含量记为LWC s ，其中：
Figure 342061DEST_PATH_IMAGE012
；
第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道内的液滴质量分数记为MF s ，其中：
Figure 542098DEST_PATH_IMAGE013
，
Figure 801041DEST_PATH_IMAGE014
，M为重构后的全液滴谱液滴直径测量区间的液滴尺寸通道总数。
进一步地，所述步骤S33中，重构后的全液滴谱液滴直径测量区间的第1～s液滴尺寸通道内的液滴质量分数之和记为CMF s ，其中，
Figure 837261DEST_PATH_IMAGE015
；体积中值直径记为MVD，其中，
Figure 839852DEST_PATH_IMAGE016
，s * 表示MVD所在的重构后的全液滴谱液滴直径测量区间的液滴尺寸通道序号。
本发明相对于现有技术的技术效果是：
本发明通过结冰云雾的重构，解决了现有技术中利用单一仪器只能进行较小的液滴尺寸范围的测量的技术问题，本发明的技术方案能够从不同仪器的测量结果中得到全液滴谱液滴尺寸分布及特征参数，也就是实现了结冰云雾的液滴的全尺寸测量和对结冰云雾的参数化表达。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面所描述的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种结冰云雾液滴尺寸分布的重构方法的流程示意图；
图2是本发明实施例的组合重构全液滴谱液滴直径测量区间的原理示意图；
图3是本发明实施例的组合重构后的液态水含量分数的分布图。
具体实施方式
在下文中将参考附图对本发明的各方面进行更充分的描述。然而，本发明可以具体化成许多不同形式且不应解释为局限于贯穿本发明所呈现的任何特定结构或功能。相反地，提供这些方面将使得本发明周全且完整，并且本发明将给本领域技术人员充分地传达本发明的范围。基于本文所教导的内容，本领域的技术人员应意识到，无论是单独还是结合本发明的任何其它方面实现本文所公开的任何方面，本发明的范围旨在涵盖本文中所公开的任何方面。例如，可以使用本文所提出任意数量的装置或者执行方法来实现。另外，除了本文所提出本发明的多个方面之外，本发明的范围更旨在涵盖使用其它结构、功能或结构和功能来实现的装置或方法。应可理解，其可通过权利要求的一或多个元件具体化本文所公开的任何方面。
在此使用的术语仅仅是为了描述具体实施例，而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在，但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。
在此使用的所有术语(包括技术和科学术语)具有本领域技术人员通常所理解的含义，除非另外定义。应注意，这里使用的术语应解释为具有与本说明书的上下文相一致的含义，而不应以理想化或过于刻板的方式来解释。
对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。
针对结冰风洞模拟的结冰云雾，尤其对于过冷大水滴结冰云雾，液滴直径范围覆盖微米尺度到毫米尺度，跨越三个量级，为了实现结冰云雾的液滴的全尺寸测量，本申请人根据液滴直径区间的量级选用多种测量仪器，根据仪器的测量范围，分别测量各个量级的液滴直径区间内的液滴尺寸分布，进而获得全液滴谱液滴尺寸分布及特征参数。
本申请人在通过测量仪器测量各个量级的液滴直径区间内的液滴直径和数量的过程中，需要将各个量级的液滴直径区间，根据仪器属性（分辨率等），分为多个紧密相接的液滴尺寸通道，通过测量各个液滴尺寸通道内的液滴数量、浓度等信息（包括液滴数密度、液态水含量等），以期获得全液滴谱液滴尺寸分布及特征参数。但是，本申请人发现，由于不同测量仪器之间的原理和性能不同，各个量级的液滴直径区间下的液滴尺寸通道宽度（液滴尺寸通道宽度指的该量级的液滴直径区间下的最大液滴直径与最小液滴直径之差）和液滴采样体积（液滴采样体积指的是采样周期内包含测量得到的云雾液滴的空气的总体积）存在差异，因此，无法直接从不同仪器的测量结果中得到全液滴谱液滴尺寸分布及特征参数。
如图1所示，本发明实施例提供了一种结冰云雾液滴尺寸分布的重构方法，其包括如下步骤：
步骤S10：将液滴尺寸测量区间划分为多个液滴尺寸通道，计算各个液滴尺寸通道内的液滴数密度、第一液态水含量、标准液态水含量；
由于不同仪器测量得到的液滴尺寸特征难以统一比较和分析，因此，需要首先对各仪器测量得到的液滴尺寸特征进行标准化处理，步骤S10的目的主要是标准化液滴尺寸测量区间内的液滴尺寸分布；
步骤S20：确定液滴直径测量区间内的截断液滴尺寸通道，确定搭接液滴尺寸通道，并计算搭接液滴尺寸通道内的液态水含量；
由于相邻的液滴尺寸测量区间之间，存在重叠区域，因此，在结冰云雾的重构过程中，需要考虑重叠区域的处理，步骤S20的目的主要是确定保留的重叠区域的具体区域、以及舍弃的重叠区域的具体区域；
步骤S30：重构结冰云雾液滴尺寸分布；
在重构结冰云雾之后，即获得了全液滴谱液滴尺寸分布及特征参数，实现结冰云雾的液滴的全尺寸测量，也就实现了结冰云雾的参数化表达。
进一步地，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道记为B i,j ，其中，B i,j =[D i,j-1 ，D i,j ]，i从1～L遍历，j从1～M i 遍历，i为液滴直径测量区间的序号，j为液滴尺寸通道的序号，D i,j-1 为第i个液滴直径测量区间内第j个液滴尺寸通道的最小液滴直径，D i,j 为第i个液滴直径测量区间内第j个液滴尺寸通道的最大液滴直径，L为液滴直径测量区间的个数，M i 为第i个液滴直径测量区间内的液滴尺寸通道的个数。
以下以包含三个典型液滴直径测量区间的结冰云雾为例，将液滴尺寸测量区间划分为多个液滴尺寸通道，如下表1所示：
表1
Figure 894396DEST_PATH_IMAGE017
从上表中可以看出：三个液滴直径测量区间均包括11个液滴尺寸通道；第1种液滴直径测量区的液滴尺寸通道为10~120μm，第2种液滴直径测量区的液滴尺寸通道为50~600μm，第3种液滴直径测量区的液滴尺寸通道为350~200μm；
进一步地，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道内的液滴数密度记为C i,j ，其中，C i,j =N i,j /SV i,j ，i从1～L遍历，j从1～M i 遍历，N i,j 为第i个液滴直径测量区间内的第j个液滴尺寸通道内的液滴数量，SV i,j 为第i个液滴直径测量区间内的第j个液滴尺寸通道内的采样体积。
在计算出液滴数密度记为C i,j 后，便可以计算第i个液滴直径测量区间内的第j个液滴尺寸通道内的液态水含量，具体地，在所述步骤S10中，将第i个液滴直径测量区间内的第j个液滴尺寸通道内的第一液态水含量记为LWC i,j ，其中，
Figure 855399DEST_PATH_IMAGE018
，
Figure 97024DEST_PATH_IMAGE019
，i从1～L遍历，j从1～M i 遍历，
Figure 168885DEST_PATH_IMAGE020
为第i个液滴直径测量区间内的第j个液滴尺寸通道内的平均液滴直径，ρ l 为液滴密度。
在计算出LWC i,j 后，便可以计算第i个液滴直径测量区间内的第j个液滴尺寸通道内的标准液态水含量，具体地，将第i个液滴直径测量区间内的第j个液滴尺寸通道内的标准液态水含量记为NLWC i,j ，其中，
Figure 546777DEST_PATH_IMAGE021
，i从1～L遍历，j从1～M i 遍历。
下表2展示了三个液滴直径测量区间各个液滴尺寸通道内的液滴数密度、液态水含量：
表2
Figure 199387DEST_PATH_IMAGE022
进一步地，所述步骤S20中，第i+1个液滴直径测量区间内的截断液滴尺寸通道记为
Figure 459467DEST_PATH_IMAGE023
，其中，
Figure 803861DEST_PATH_IMAGE024
，i从1～（L-1）遍历，k i 为第i个重叠区域对应的第i+1个液滴直径测量区间内的液滴尺寸通道的序号，k i 从1～M i+1 遍历；
进一步参见表1，可以得出，第2个液滴直径测量区间内的截断液滴尺寸通道为Bin 2,2 =[100,150]，第3个液滴直径测量区间内的截断液滴尺寸通道为Bin 3,2 =[500,650]；
进一步地，所述步骤S20中，第i个搭接液滴尺寸通道记为
Figure 832997DEST_PATH_IMAGE025
，第i个搭接液滴尺寸通道内的液态水含量记为
Figure 339064DEST_PATH_IMAGE026
，其中，
Figure 102752DEST_PATH_IMAGE028
，
Figure 516416DEST_PATH_IMAGE029
，i从1～（L-1）遍历，k i 从1～M i+1 遍历。
进一步参见表1，可以得出，第1个搭接液滴尺寸通道为Bin * 1 =[120, 150]，第2个搭接液滴尺寸通道为Bin * 2 =[600,650]；第1个搭接液滴尺寸通道内的液态水含量LWC * 1 =0.224 ×lg(150/120)=0.0217g/m 3 ，第2个搭接液滴尺寸通道内的液态水含量LWC * 2 =0.475×lg (650/600)=0.0165 g/m 3 。
在完成上述步骤之后，则可以正式开始结冰云雾的重构，具体地，所述步骤S30包括如下步骤：
步骤S31：组合重构全液滴谱液滴直径测量区间；
步骤S32：计算全液滴谱液滴直径测量区间内各通道内的液态水含量、液态水含量分数；
步骤S33：计算全液滴谱液滴直径测量区间的累积质量分数和体积中值直径。
进一步地，所述步骤S31中，重构后的全液滴谱液滴直径测量区间的液滴尺寸通道记为Bin s ，其中：
Figure 868900DEST_PATH_IMAGE030
s为重构后的全液滴谱液滴直径测量区间的液滴尺寸通道序号，D s-1为第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最小值，D s 为第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最大值。
如图2所示为组合重构全液滴谱液滴直径测量区间的原理示意图，首先是计算出第i+1个液滴直径测量区间内的截断液滴尺寸通道，再将第i+1个液滴直径测量区间内位于D i,Mi 左侧的部分舍弃，将第i+1个液滴直径测量区间内位于D i,Mi 右侧的部分保留，形成新的液滴直径测量区间。
进一步地，所述步骤S32中，第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道内的液态水含量记为LWC s ，其中：
Figure 342607DEST_PATH_IMAGE031
；
第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道内的液滴质量分数记为MF s ，其中：
Figure 780541DEST_PATH_IMAGE013
，
Figure 247164DEST_PATH_IMAGE014
，M为重构后的全液滴谱液滴直径测量区间的液滴尺寸通道总数。
如下表3所示为组合重构的全液滴谱液滴直径测量区间及全液滴谱液滴直径测量区间内各通道内的液态水含量、液态质量分数：
表3
Figure 250892DEST_PATH_IMAGE032
从上表可以看出，组合重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的个数为31个，第12个和22个液滴尺寸通道为新增的搭接液滴尺寸通道。
进一步地，所述步骤S33中，重构后的全液滴谱液滴直径测量区间的第1～s液滴尺寸通道内的液滴质量分数之和记为CMF s ，其中，
Figure 364341DEST_PATH_IMAGE015
；体积中值直径记为MVD，其中，
Figure 86309DEST_PATH_IMAGE033
，s * 表示MVD所在的重构后的全液滴谱液滴直径测量区间的液滴尺寸通道序号。
如图3所示为组合重构后的液滴尺寸分布图，其中，横坐标为全液滴谱液滴直径，纵坐标为累计质量分数，同时，计算出的MVD=350um。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种结冰云雾液滴尺寸分布的重构方法，其特征在于，包括如下步骤：
步骤S10：将液滴尺寸测量区间划分为多个液滴尺寸通道，计算各个液滴尺寸通道内的液滴数密度、第一液态水含量、标准液态水含量；
步骤S20：确定液滴直径测量区间内的截断液滴尺寸通道，确定搭接液滴尺寸通道，并计算搭接液滴尺寸通道内的液态水含量；
步骤S30：重构结冰云雾液滴尺寸分布；
所述步骤S30包括如下步骤：
步骤S31：组合重构全液滴谱液滴直径测量区间；
步骤S32：计算全液滴谱液滴直径测量区间内各通道内的液态水含量、液态水含量分数；
步骤S33：计算全液滴谱液滴直径测量区间的累积质量分数和体积中值直径。
2.如权利要求1所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道记为B i,j ，其中，B i,j =[D i,j-1 ，D i,j ]，i从1～L遍历，j从1～M i 遍历，i为液滴直径测量区间的序号，j为液滴尺寸通道的序号，D i,j-1 为第i个液滴直径测量区间内第j个液滴尺寸通道的最小液滴直径，D i,j 为第i个液滴直径测量区间内第j个液滴尺寸通道的最大液滴直径，L为液滴直径测量区间的个数，M i 为第i个液滴直径测量区间内的液滴尺寸通道的个数。
3.如权利要求2所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道内的液滴数密度记为C i,j ，其中，C i,j =N i,j /SV i,j ，i从1～L遍历，j从1～M i 遍历，N i,j 为第i个液滴直径测量区间内的第j个液滴尺寸通道内的液滴数量，SV i,j 为第i个液滴直径测量区间内的第j个液滴尺寸通道内的采样体积。
4.如权利要求3所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S10中，第i个液滴直径测量区间内的第j个液滴尺寸通道内的第一液态水含量记为LWC i,j ，第i个液滴直径测量区间内的第j个液滴尺寸通道内的标准液态水含量记为NLWC i,j ，其中：
Figure DEST_PATH_IMAGE001
，
Figure DEST_PATH_IMAGE002
，
Figure DEST_PATH_IMAGE003
，i从1～L遍历，j从1～M i 遍历，
Figure DEST_PATH_IMAGE004
为第i个液滴直径测量区间内的第j个液滴尺寸通道内的平均液滴直径，ρ l 为液滴密度。
5.如权利要求4所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S20中，第i+1个液滴直径测量区间内的截断液滴尺寸通道记为
Figure DEST_PATH_IMAGE005
，其中，
Figure DEST_PATH_IMAGE006
，
Figure DEST_PATH_IMAGE007
，
Figure DEST_PATH_IMAGE008
为第i+1个液滴直径测量区间内第k i 个液滴尺寸通道的最小液滴直径，
Figure DEST_PATH_IMAGE009
为第i个液滴直径测量区间内第M i 个液滴尺寸通道的最大液滴直径，
Figure DEST_PATH_IMAGE010
为第i+1个液滴直径测量区间内第k i 个液滴尺寸通道的最大液滴直径，i从1～（L-1）遍历，k i 为第i个重叠区域对应的第i+1个液滴直径测量区间内的液滴尺寸通道的序号，k i 从1～M i+1 遍历。
6.如权利要求5所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S20中，第i个搭接液滴尺寸通道记为
Figure DEST_PATH_IMAGE011
，第i个搭接液滴尺寸通道内的液态水含量记为
Figure DEST_PATH_IMAGE012
，其中，
Figure DEST_PATH_IMAGE013
，
Figure DEST_PATH_IMAGE014
，i从1～（L-1）遍历，k i 从1～M i+1 遍历。
7.如权利要求6所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S31中，重构后的全液滴谱液滴直径测量区间的液滴尺寸通道记为Bin s ，其中：
Figure DEST_PATH_IMAGE015
s为重构后的全液滴谱液滴直径测量区间的液滴尺寸通道序号，D s-1为第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最小值，D s 为第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最大值，x为重构序号，D i,x-1 为第i个液滴直径测量区间内第x个液滴尺寸通道的最小液滴直径，D i,x 为第i个液滴直径测量区间内第x个液滴尺寸通道的最大液滴直径，
Figure DEST_PATH_IMAGE016
为第i-1个液滴直径测量区间内第M i-1 个液滴尺寸通道的最大液滴直径，
Figure DEST_PATH_IMAGE017
为第i个液滴直径测量区间内第k i-1 个液滴尺寸通道的最大液滴直径，
Figure DEST_PATH_IMAGE018
为第i个液滴直径测量区间内第k i-1 +x个液滴尺寸通道的最小液滴直径，
Figure DEST_PATH_IMAGE019
为第i个液滴直径测量区间内第k i-1 +x个液滴尺寸通道的最大液滴直径。
8.如权利要求7所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S32中，第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道内的液态水含量记为LWC s ，其中：
Figure DEST_PATH_IMAGE020
；
第s个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道内的液滴质量分数记为MF s ，其中：
Figure DEST_PATH_IMAGE021
，
Figure DEST_PATH_IMAGE022
，M为重构后的全液滴谱液滴直径测量区间的液滴尺寸通道总数；
Figure DEST_PATH_IMAGE023
为第i个液滴直径测量区间内第x个液滴尺寸通道内的液态水含量，
Figure DEST_PATH_IMAGE024
为第i个液滴直径测量区间内第k i-1 +x个液滴尺寸通道内的液态水含量，
Figure DEST_PATH_IMAGE025
为第i-1个搭接液滴尺寸通道内的液态水含量。
9.如权利要求8所述的一种结冰云雾液滴尺寸分布的重构方法，其特征在于，所述步骤S33中，重构后的全液滴谱液滴直径测量区间的第1～s液滴尺寸通道内的液滴质量分数之和记为CMF s ，其中，
Figure DEST_PATH_IMAGE026
；体积中值直径记为MVD，其中，
Figure DEST_PATH_IMAGE027
，
Figure DEST_PATH_IMAGE029
，s * 表示MVD所在的重构后的全液滴谱液滴直径测量区间的液滴尺寸通道序号，
Figure DEST_PATH_IMAGE031
为第s * 个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最大值，
Figure DEST_PATH_IMAGE033
为第s * 个重构后的全液滴谱液滴直径测量区间的液滴尺寸通道的液滴直径最小值，
Figure DEST_PATH_IMAGE034
为第1～（s * -1）液滴尺寸通道内的液滴质量分数之和，
Figure DEST_PATH_IMAGE035
为第1～s * 液滴尺寸通道内的液滴质量分数之和。
CN202110611640.6A 2021-06-02 2021-06-02 一种结冰云雾液滴尺寸分布的重构方法 Active CN113252281B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110611640.6A CN113252281B (zh) 2021-06-02 2021-06-02 一种结冰云雾液滴尺寸分布的重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110611640.6A CN113252281B (zh) 2021-06-02 2021-06-02 一种结冰云雾液滴尺寸分布的重构方法

Publications (2)

Publication Number Publication Date
CN113252281A CN113252281A (zh) 2021-08-13
CN113252281B true CN113252281B (zh) 2021-09-21

Family

ID=77186030

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110611640.6A Active CN113252281B (zh) 2021-06-02 2021-06-02 一种结冰云雾液滴尺寸分布的重构方法

Country Status (1)

Country Link
CN (1) CN113252281B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113899663B (zh) * 2021-12-08 2022-02-18 中国空气动力研究与发展中心低速空气动力研究所 一种大型结冰风洞中复杂结冰云雾mvd测量方法
CN117890068B (zh) * 2024-03-14 2024-05-24 中国空气动力研究与发展中心低速空气动力研究所 一种结冰云雾颗粒尺寸重构方法及计算机可读存储介质

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5567865A (en) * 1994-09-14 1996-10-22 Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. Method for determining the size of airborne water droplets
CN101281046A (zh) * 2008-05-16 2008-10-08 重庆大学 大气结构物覆冰参数测量方法及多导体积冰器
CN103076290A (zh) * 2011-10-25 2013-05-01 罗斯蒙德宇航公司 确定云的液态水含量的方法
CN103072695A (zh) * 2011-10-25 2013-05-01 罗斯蒙德宇航公司 确定云的水滴大小分布的方法
JP2013144274A (ja) * 2012-01-13 2013-07-25 Sharp Corp 液滴量測定方法、粒度分布測定方法および液滴吐出装置
CN103308271A (zh) * 2012-03-13 2013-09-18 波音公司 过度冷却的大滴状物结冰条件的仿真系统
CN104117447A (zh) * 2013-04-24 2014-10-29 武汉航空仪表有限责任公司 一种用于结冰试验的喷雾参数控制方法
WO2016182992A1 (en) * 2015-05-12 2016-11-17 Croda, Inc. Method of analyzing spray particulates
CN109767518A (zh) * 2018-12-14 2019-05-17 南京航空航天大学 基于机翼结冰厚度确定气象参数mvd的反推方法
CN111199111A (zh) * 2020-01-15 2020-05-26 浙江大学 一种无人机双喷头雾滴粒径分布模拟方法
CN111291311A (zh) * 2020-05-06 2020-06-16 中国空气动力研究与发展中心低速空气动力研究所 一种积冰密度测量方法
CN111967201A (zh) * 2020-07-30 2020-11-20 华南理工大学 一种基于数值仿真模型分析临界覆冰类型的方法
CN112507572A (zh) * 2021-02-07 2021-03-16 中国空气动力研究与发展中心低速空气动力研究所 一种过冷大水滴结冰条件的优化匹配评估方法
CN112513611A (zh) * 2018-06-25 2021-03-16 巴斯夫涂料有限公司 油漆开发中的雾化期间液滴尺寸分布的测定方法和基于其的筛选方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120097543A1 (en) * 2010-10-25 2012-04-26 Conocophillips Company Method for predicting drop size distribution
CN104268399B (zh) * 2014-09-24 2017-03-15 空气动力学国家重点实验室 过冷大水滴条件下结冰风洞试验中模型参数的计算方法
CN105728229A (zh) * 2016-03-01 2016-07-06 中国飞机强度研究所 一种基于喷嘴动力参数的喷雾水滴平均直径快速估算方法
CN111310381B (zh) * 2020-05-14 2020-07-31 中国空气动力研究与发展中心低速空气动力研究所 一种三维水滴收集系数计算方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5567865A (en) * 1994-09-14 1996-10-22 Deutsche Forschungsanstalt Fur Luft-Und Raumfahrt E.V. Method for determining the size of airborne water droplets
CN101281046A (zh) * 2008-05-16 2008-10-08 重庆大学 大气结构物覆冰参数测量方法及多导体积冰器
CN103076290A (zh) * 2011-10-25 2013-05-01 罗斯蒙德宇航公司 确定云的液态水含量的方法
CN103072695A (zh) * 2011-10-25 2013-05-01 罗斯蒙德宇航公司 确定云的水滴大小分布的方法
JP2013144274A (ja) * 2012-01-13 2013-07-25 Sharp Corp 液滴量測定方法、粒度分布測定方法および液滴吐出装置
CN103308271A (zh) * 2012-03-13 2013-09-18 波音公司 过度冷却的大滴状物结冰条件的仿真系统
CN104117447A (zh) * 2013-04-24 2014-10-29 武汉航空仪表有限责任公司 一种用于结冰试验的喷雾参数控制方法
WO2016182992A1 (en) * 2015-05-12 2016-11-17 Croda, Inc. Method of analyzing spray particulates
CN112513611A (zh) * 2018-06-25 2021-03-16 巴斯夫涂料有限公司 油漆开发中的雾化期间液滴尺寸分布的测定方法和基于其的筛选方法
CN109767518A (zh) * 2018-12-14 2019-05-17 南京航空航天大学 基于机翼结冰厚度确定气象参数mvd的反推方法
CN111199111A (zh) * 2020-01-15 2020-05-26 浙江大学 一种无人机双喷头雾滴粒径分布模拟方法
CN111291311A (zh) * 2020-05-06 2020-06-16 中国空气动力研究与发展中心低速空气动力研究所 一种积冰密度测量方法
CN111967201A (zh) * 2020-07-30 2020-11-20 华南理工大学 一种基于数值仿真模型分析临界覆冰类型的方法
CN112507572A (zh) * 2021-02-07 2021-03-16 中国空气动力研究与发展中心低速空气动力研究所 一种过冷大水滴结冰条件的优化匹配评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"机翼结冰模拟中关键问题的研究";孟繁鑫;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20160131;全文 *
郭向东等." 结冰风洞中过冷大水滴云雾演化特性数值研究".《航空学报》.2020, *

Also Published As

Publication number Publication date
CN113252281A (zh) 2021-08-13

Similar Documents

Publication Publication Date Title
CN113252281B (zh) 一种结冰云雾液滴尺寸分布的重构方法
Raupach et al. Correction of raindrop size distributions measured by Parsivel disdrometers, using a two-dimensional video disdrometer as a reference
Gathman Optical properties of the marine aerosol as predicted by the Navy aerosol model
Strapp et al. Wind tunnel measurements of the response of hot-wire liquid water content instruments to large droplets
CN109901204A (zh) 一种基于伪距误差分布模型的gbas完好性性能评估方法
FR2960652A1 (fr) Procede de correction de mesures de reflectivite par detection d'isotherme et radar mettant en oeuvre ce procede
Stone et al. Atmospheric river reconnaissance observation impact in the Navy Global Forecast System
Kanagaratnam et al. A wideband radar for high-resolution mapping of near-surface internal layers in glacial ice
CN109492063A (zh) 一种高精度地图数据的质量评价方法
CN112556917B (zh) 一种利用测压装置进行测压的方法
CN111007541B (zh) 一种卫星导航地基增强系统仿真性能评估方法
Ma et al. Assimilation of all-sky radiance from the FY-3 MWHS-2 with the Yinhe 4D-Var system
Misra et al. SMAP RFI mitigation algorithm performance characterization using airborne high-rate direct-sampled SMAPVEX 2012 data
CN101192253A (zh) 用于燃气轮机性能诊断的状态初始化
CN206146647U (zh) 一种汽车尾气流量测量控制系统
Yano et al. A minimum bulk microphysics
King-Steen et al. NASA Glenn Icing Research Tunnel: 2018 change in drop-sizing equations due to change in cloud droplet probe sample area
Reitter et al. Evaluation of ice and snow content in the global numerical weather prediction model GME with CloudSat
CN114076727A (zh) 一种基于电阻率的冰的孔隙率测量方法
Anderson et al. Characterization of aerosols from simulated SAGE III measurements applying two retrieval techniques
CN108680934A (zh) 卫星导航信号的动态范围检测系统及检测方法
Rapone et al. Characterization of real world bus driving behavior for emission evaluation
Li et al. Influence of changes in surface layer properties on tire/pavement noise
Lean et al. ESA Contract Report
Duncan et al. EUMETSAT Contract Report

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