CN114858849B - 一种动态冰的导热系数获取方法 - Google Patents

一种动态冰的导热系数获取方法 Download PDF

Info

Publication number
CN114858849B
CN114858849B CN202210807281.6A CN202210807281A CN114858849B CN 114858849 B CN114858849 B CN 114858849B CN 202210807281 A CN202210807281 A CN 202210807281A CN 114858849 B CN114858849 B CN 114858849B
Authority
CN
China
Prior art keywords
hole
statistical
unit
ith
heat conduction
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
CN202210807281.6A
Other languages
English (en)
Other versions
CN114858849A (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 CN202210807281.6A priority Critical patent/CN114858849B/zh
Publication of CN114858849A publication Critical patent/CN114858849A/zh
Application granted granted Critical
Publication of CN114858849B publication Critical patent/CN114858849B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/18—Investigating or analyzing materials by the use of thermal means by investigating thermal conductivity

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本申请实施例公开了一种动态冰的导热系数获取方法，涉及动态冰特性研究技术领域。该方法包括：将动态冰划分为多个统计单元，将每个统计单元划分为多个网格单元，且划分方式不随孔的大小和分布而变化；建立统计单元中每个孔影响下的导热模型，并建立每个孔的影响范围在每个网格单元中的面积权重；根据每个孔影响下的导热模型和每个孔的影响范围在每个网格单元中的面积权重，得到统计单元中的每个网格单元的等效导热值；根据每个网格单元的等效导热值得到统计单元的平均导热值；根据多个统计单元的平均导热值得到动态冰的导热系数。通过上述方法，可以缓解获取动态冰导热系数的方法中，存在的孔隙大小相差较大时研究方法复杂等问题。

Description

一种动态冰的导热系数获取方法
技术领域
本申请涉及动态冰特性研究技术领域，更具体地，涉及一种动态冰的导热系数获取方法。
背景技术
过冷水撞击低温基底结冰，形成的带有非均匀分布的孔隙结构的冰即动态冰。导热系数是材料的固有属性之一，研究动态冰的导热系数是进一步研究动态冰特性的重要基础。
动态冰的孔隙结构是决定其物理特性的根本因素。由于动态冰中的孔隙密度和孔隙分布均受到结冰环境的影响，因此没有能够准确获取动态冰导热系数的方法。
动态冰的内部存在大量的非均匀分布的孔隙结构，因此动态冰可以被当做多孔异质材料。目前，获取异质材料导热系数的方法，主要包括有限元研究法。有限元研究法通过将异质材料划分成多个有限元网格，研究每个有限元网格的导热系数，最终获得整个异质材料的导热系数，但目前的有限元研究法在异质材料的孔隙大小相差较大时，存在计算复杂的问题。如专利CN111027236B，在划分有限元网格时，将不同的圆形孔隙结构划分成多个三角形有限元网格，通过得到每个有限元网格中的导热系数，最终获得整个材料的导热系数，该方案中，若材料的孔隙大小相差较小，在保证精度的情形下，可以使用较少的网格进行计算，但若材料的孔隙大小相差较大，为了保证精度，就需要将较大的孔隙划分更多的网格，因此计算过程会更复杂。
因此，现有技术在获取动态冰的导热系数时，存在孔隙大小相差较大时研究方法复杂等问题。
发明内容
本申请提出了一种动态冰的导热系数获取方法，将动态冰划分为多个统计单元，再将每个统计单元划分为多个网格单元，且划分方式不随动态冰中的孔的大小、形状和分布而变化，建立统计单元中的每个孔影响下的导热模型，以及每个孔的影响范围在每个网格单元中的面积权重，根据每个孔影响下的导热模型和每个孔的影响范围在每个网格单元中的面积权重，得到每个网格单元的等效导热值，在孔隙大小相差较大时，也可以根据导热模型和面积权重算出每个网格单元的等效导热值，无需划分更多的网格，从而简化网格单元的等效导热值的获取过程，然后根据每个网格单元的等效导热值得到统计单元的平均导热值，再根据多个统计单元的平均导热值得到动态冰的导热系数。如此，可以有效解决现有技术在获取动态冰的导热系数时，存在的孔隙大小相差较大时研究方法复杂等问题。
本申请实施例提供了一种动态冰的导热系数获取方法，该方法包括：S110.将所述动态冰划分为多个统计单元，将每个统计单元划分为多个网格单元，且划分方式不随孔的大小和分布而变化；S120.若所述统计单元中包含至少一个孔，则根据所述统计单元中每个孔的影响范围，建立每个孔影响下的导热模型，并建立每个孔的影响范围在每个网格单元中的面积权重，其中，每个孔影响下的导热模型为
Figure 820248DEST_PATH_IMAGE001
，i为所述统计单元中的孔的编号，
Figure 985650DEST_PATH_IMAGE002
为所述统计单元中的坐标点，且
Figure 20602DEST_PATH_IMAGE001
同时满足下述条件：对任意
Figure 650166DEST_PATH_IMAGE003
，
Figure 479582DEST_PATH_IMAGE004
，其中，
Figure 460176DEST_PATH_IMAGE005
为第i个孔所在的区域，
Figure 829978DEST_PATH_IMAGE006
为空气的导热系数，
Figure 454994DEST_PATH_IMAGE007
为所述第i个孔的面积，对任意
Figure 579945DEST_PATH_IMAGE008
，
Figure 657622DEST_PATH_IMAGE009
，其中，
Figure 96694DEST_PATH_IMAGE010
为所述第i个孔的影响范围所在的区域，
Figure 700851DEST_PATH_IMAGE011
为第i个孔的模拟导热值，
Figure 872069DEST_PATH_IMAGE012
为所述第i个孔的影响范围的面积，对任意
Figure 561676DEST_PATH_IMAGE013
，
Figure 476543DEST_PATH_IMAGE014
，
Figure 872889DEST_PATH_IMAGE015
为所述第i个孔的影响范围的边界，
Figure 616940DEST_PATH_IMAGE016
为静态冰的导热系数，对任意
Figure 403631DEST_PATH_IMAGE017
，
Figure 246822DEST_PATH_IMAGE018
；S130.根据所述每个孔影响下的导热模型和所述每个孔的影响范围在每个网格单元中的面积权重，得到所述统计单元中的每个网格单元的等效导热值；S140.根据所述每个网格单元的等效导热值得到所述统计单元的平均导热值；S150.根据多个所述统计单元的平均导热值得到所述动态冰的导热系数。
综上所述，本申请至少具有如下技术效果：
1.本申请通过将动态冰划分为多个统计单元，再将每个统计单元划分为多个网格单元，且划分方式不随孔的大小和分布而变化，为获取动态冰的导热系数建立较为简单的前提条件。
2.本申请通过建立每个孔影响下的导热模型，并建立每个孔的影响范围在每个网格单元中的面积权重，根据每个孔影响下的导热模型和每个孔的影响范围在每个网格单元中的面积权重，得到统计单元中的每个网格单元的等效导热值，在孔隙大小相差较大时，也可以根据导热模型和面积权重算出每个网格单元的等效导热值，无需划分更多的网格，从而使获取动态冰的导热系数的过程更加简单。
3.本申请通过使导热模型满足在孔的区域内以及在孔的影响范围内的积分条件，并使导热模型满足在孔的影响范围边界上以及孔的影响范围外的条件，从而建立每个孔影响下的导热模型，为获取动态冰的导热系数建立较为简单的前提条件。
4.本申请通过将每个孔影响下的导热模型与该孔在网格单元中的面积权重相乘，得到每个孔对该网格单元的等效导热值的影响值，再将每个孔的影响值累加起来，得到该网格单元的等效导热值，如此对每个网格单元进行计算，从而得到每个网格单元的等效导热值，为获取动态冰的导热系数建立较为简单的前提条件。
因此，本申请提供的方案可以有效解决现有技术在获取动态冰的导热系数时，存在的孔隙大小相差较大时研究方法复杂等问题。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出了本申请实施例1提供的一种动态冰的导热系数获取方法的流程示意图；
图2示出了本申请实施例1提供的动态冰的其中一个统计单元的示意图；
图3示出了本申请实施例1提供的孔及其影响范围所在区域的示意图；
图4示出了本申请实施例1提供的动态冰的微观结构显微照片的示意图；
图5示出了本申请实施例1提供的裂隙及其影响范围的面积的示意图；
图6示出了本申请实施例1提供的孔的影响范围在网格单元中的面积的示意图；
图7示出了本申请实施例1提供的等温和绝热边界条件的示意图；
图8示出了本申请实施例1提供的热流示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
以下，首先对本申请涉及的技术术语进行介绍。
动态冰：过冷水撞击低温基底结冰，形成的带有非均匀分布的孔隙结构的冰即动态冰，其中，孔隙结构可以包括气泡和裂隙。动态冰的生成可以是在飞机结冰的过程中。飞机结冰是飞行安全的主要威胁之一，飞机穿越包含过冷水滴的云层时，飞机的迎风表面会收集过冷水滴，过冷水滴受到扰动生成动态冰，进而破坏飞机的气动外形，因此，研究动态冰的特性，可以帮助理解飞机结冰和飞机防冰的物理过程，还可以帮助提升防冰系统设计精度，以及建立防冰系统和除冰系统的相似理论。动态冰的生成过程受到气流、物体外形、过冷水滴直径和当地温度等因素影响，不同条件下的动态冰的宏观外形和内部微观结构存在较大差异。通常将飞机生成的动态冰按照外观分为霜冰、明冰和混合冰，三种类型的根本差异在于其内部微观孔隙结构导致的密度差异。动态冰中的孔隙密度和孔隙分布均受到结冰环境的影响，而动态冰的孔隙结构是决定其物理特性的根本因素。在研究动态冰的特性时，由于动态冰的内部存在大量的非均匀分布的孔隙结构，因此动态冰可以被当做多孔异质材料。
静态冰：静态冰是指内部不含有孔隙结构的纯冰，静态冰是典型的多晶体，具有各向同性的物理性质及固定的熔点，且静态冰的导热系数约为2.2W/(m·K)。
动态冰的内部存在大量的非均匀分布的孔隙结构，因此动态冰可以被当做多孔异质材料。目前，获取多孔异质材料导热系数的方法，主要包括有限元研究法。有限元研究法通过将异质材料划分成多个有限元网格，研究每个有限元网格的导热系数，最终获得整个异质材料的导热系数，但目前的有限元研究法在对异质材料划分有限元网格时，有限元网格的大小、形状、边界、排列都随孔隙结构而变化，若材料的孔隙大小相差较小，在保证精度的情形下，可以使用较少的网格进行计算，但若材料的孔隙大小相差较大，如材料中存在微米量级和毫米量级的孔隙，为了保证精度，就需要将较大的孔隙划分更多的网格，计算过程比较复杂。
因此，为了解决上述缺陷，本申请实施例提供了一种动态冰的导热系数获取方法，该方法包括：将动态冰划分为多个统计单元，再将每个统计单元划分为多个网格单元，且划分方式不随动态冰中的孔的大小、形状和分布而变化，建立统计单元中的每个孔影响下的导热模型，以及每个孔的影响范围在每个网格单元中的面积权重，根据每个孔影响下的导热模型和每个孔的影响范围在每个网格单元中的面积权重，得到每个网格单元的等效导热值，从而简化网格单元的等效导热值的获取过程，然后根据每个网格单元的等效导热值得到统计单元的平均导热值，再根据多个统计单元的平均导热值得到动态冰的导热系数。从而有效解决现有技术在获取动态冰的导热系数时，存在的孔隙大小相差较大时研究方法复杂等问题。
下面对本申请所涉及到的动态冰的导热系数获取方法进行介绍。
实施例1
请参照图1，图1为本申请实施例1提供的一种动态冰的导热系数获取方法的流程示意图。本实施例中，该动态冰的导热系数获取方法可以包括以下步骤：
步骤S110：将所述动态冰划分为多个统计单元，将每个统计单元划分为多个网格单元，且划分方式不随孔的大小和分布而变化。
在本申请实施例中，可以根据动态冰材料的大小划分统计单元的个数，也可以根据所需精度划分统计单元个数，本申请在此不做限制，统计单元的个数为两个或两个以上。
统计单元的形状可以是正方形或长方形。
在本申请实施例中，每个统计单元中的网格单元的个数可以相等，也可以不相等，每个统计单元中网格的排列方式可以相同，也可以不相同。
网格单元的形状可以是任意类型的结构网格或非结构网格。
在本申请实施例中，网格单元的划分与动态冰中的孔的大小、形状、分布均没有关系。
作为一种可选实施方式，本申请以统计单元和网格单元是长方形为例进行说明，如图2所示，图2示出了动态冰的其中一个统计单元，该统计单元被划分为5行、6列网格单元，该统计单元包含7个孔。
在示例性实施例中，步骤S110可以包括子步骤S111。
子步骤S111：将所述动态冰划分为多个特征长度为
Figure 497675DEST_PATH_IMAGE019
的所述统计单元，且
Figure 276275DEST_PATH_IMAGE020
，其中，
Figure 674895DEST_PATH_IMAGE021
为所述动态冰的最小微观结构特征长度。
在该实施例中，根据最小微观结构特征长度划分统计单元的大小，根据经验可以得知，当
Figure 197143DEST_PATH_IMAGE020
时，实施本申请的方案更简单，也能保证结果精度。
本申请通过将动态冰划分为多个统计单元，再将每个统计单元划分为多个网格单元，且划分方式不随动态冰中的孔的大小和分布而变化，为获取动态冰的导热系数建立较为简单的前提条件。
步骤S120：若所述统计单元中包含至少一个孔，则根据所述统计单元中每个孔的影响范围，建立每个孔影响下的导热模型，并建立每个孔的影响范围在每个网格单元中的面积权重，其中，每个孔影响下的导热模型为
Figure 302503DEST_PATH_IMAGE001
，i为所述统计单元中的孔的编号，
Figure 376638DEST_PATH_IMAGE002
为所述统计单元中的坐标点，且
Figure 137920DEST_PATH_IMAGE001
同时满足下述条件：
对任意
Figure 322914DEST_PATH_IMAGE003
，
Figure 486042DEST_PATH_IMAGE004
，其中，
Figure 403183DEST_PATH_IMAGE005
为第i个孔所在的区域，
Figure 41974DEST_PATH_IMAGE006
为空气的导热系数，
Figure 640446DEST_PATH_IMAGE007
为所述第i个孔的面积，
对任意
Figure 782714DEST_PATH_IMAGE008
，
Figure 136335DEST_PATH_IMAGE009
，其中，
Figure 872210DEST_PATH_IMAGE010
为所述第i个孔的影响范围所在的区域，
Figure 133427DEST_PATH_IMAGE011
为第i个孔的模拟导热值，
Figure 271147DEST_PATH_IMAGE012
为所述第i个孔的影响范围的面积，
对任意
Figure 795670DEST_PATH_IMAGE013
，
Figure 877895DEST_PATH_IMAGE014
，
Figure 349328DEST_PATH_IMAGE015
为所述第i个孔的影响范围的边界，
Figure 200609DEST_PATH_IMAGE016
为静态冰的导热系数，
对任意
Figure 833716DEST_PATH_IMAGE017
，
Figure 340921DEST_PATH_IMAGE018
。
在本申请实施例中，每个统计单元中的孔的数量、形状和大小，可以通过对动态冰的微观结构显微照片进行图像识别得到，也可以通过随机模型生成，如通过高斯分布或对数正态分布等方式生成，还可以根据对材料微观结构的解析，使用马尔科夫链-蒙特卡洛方法（Markov Chain-Monte Carlo）或类似方法获得更加真实的多孔结构。
由于动态冰内部的孔隙结构分布不均匀，有的统计单元可能包含一个或多个孔，有的统计单元可能不包含孔。在示例性实施例中，若所述统计单元中不包含孔，则将静态冰的导热系数作为所述统计单元的平均导热值。
在本申请实施例中，统计单元中的每个孔对应一个影响范围。
在本申请实施例中，统计单元中的每个孔都会对该统计单元中的各个点产生影响，因此每个孔对该统计单元中的各个点都对应有一个导热模型，第i个孔对该统计单元中的各个点产生影响的导热模型为
Figure 475099DEST_PATH_IMAGE001
。
如图3所示，图3示出了其中一个统计单元中的第1个孔和第2个孔，
Figure 56253DEST_PATH_IMAGE022
为第1个孔所在的区域，
Figure 250474DEST_PATH_IMAGE023
为第1个孔的影响范围所在的区域，
Figure 979395DEST_PATH_IMAGE024
为第2个孔所在的区域，
Figure 792630DEST_PATH_IMAGE025
为第2个孔的影响范围所在的区域。
每个孔影响下的导热模型需要在四个定义域上满足以下四个条件。
在该孔内部，导热模型在该孔所在区域上的积分等于空气的导热系数与该孔的面积的乘积，即对任意
Figure 618504DEST_PATH_IMAGE026
，有：
Figure 593413DEST_PATH_IMAGE027
。
在该孔的影响范围内，导热模型在该孔的影响范围所在的区域上的积分，等于该孔的模拟导热值与该孔的影响范围的面积的乘积，即对任意
Figure 75210DEST_PATH_IMAGE028
，
Figure 551191DEST_PATH_IMAGE029
。其中，
Figure 372516DEST_PATH_IMAGE030
为第i个孔的模拟导热值，示例性的，对于圆形或可以等效为圆形的孔，其模拟导热值可以使用该公式计算：
Figure 377381DEST_PATH_IMAGE031
，
Figure 284158DEST_PATH_IMAGE032
，
Figure 235933DEST_PATH_IMAGE033
为第i个圆形或可以等效为圆形的孔的等效半径，
Figure 36399DEST_PATH_IMAGE034
为第i个圆形或可以等效为圆形的孔的影响范围的半径。
在该孔的影响范围的边界上，导热模型的取值无限趋近于静态冰的导热系数，即对任意
Figure 353111DEST_PATH_IMAGE035
，
Figure 871817DEST_PATH_IMAGE036
。
第四，在该孔的影响范围外，导热模型的取值等于静态冰的导热系数，即对任意
Figure 361704DEST_PATH_IMAGE037
，
Figure 485518DEST_PATH_IMAGE038
。
在本申请实施例中，满足上述四个条件的函数模型均可为导热模型。
在示例性实施例中，所述统计单元中的至少一个孔包括气泡和/或裂隙。
如图4所示，图4为动态冰的微观结构显微照片，动态冰中的孔可以是微米量级的气泡，也可以是毫米量级的裂隙。气泡及其影响范围的形状可以是圆形，也可以是椭圆，或者其他不规则形状。裂隙及其影响范围的形状可以是平行四边形，也可以是长方形，或者其他不规则形状。由于动态冰中的气泡和裂隙的形状并非全是规则的形状，因此可以将不规则形状简化或近似为规则的形状，再采用等效面积的方法获得气泡和裂隙的参数。
以将气泡的形状等效为圆形为例进行说明，如图3所示，若所述第i个孔为气泡，则所述
Figure 35448DEST_PATH_IMAGE007
为第i个气泡的等效面积，且
Figure 651237DEST_PATH_IMAGE039
，其中，
Figure 272711DEST_PATH_IMAGE040
为所述第i个气泡的等效半径，所述
Figure 657556DEST_PATH_IMAGE041
为所述第i个气泡的影响范围的面积，且
Figure 706283DEST_PATH_IMAGE042
，其中，
Figure 871686DEST_PATH_IMAGE043
为所述第i个气泡的影响范围的半径。
在本申请实施例中，还可以将气泡的形状等效为椭圆形或其他形状，其等效面积获取原理与将气泡的形状等效为圆形时的等效面积获取原理相同，在此不再赘述。
以将裂隙的形状等效为平行四边形为例进行说明，如图5所示，若所述第i个孔为裂隙，则所述
Figure 906638DEST_PATH_IMAGE007
为第i个裂隙的等效面积，且
Figure 536202DEST_PATH_IMAGE044
，其中，
Figure 365618DEST_PATH_IMAGE045
为所述第i个裂隙的等效底边，
Figure 80633DEST_PATH_IMAGE046
为所述第i个裂隙的等效高，所述
Figure 716014DEST_PATH_IMAGE047
为所述第i个裂隙的影响范围的面积，且
Figure 341030DEST_PATH_IMAGE048
，其中，
Figure 465981DEST_PATH_IMAGE049
为所述第i个裂隙的影响范围的底边，
Figure 543658DEST_PATH_IMAGE050
为所述第i个裂隙的影响范围的高。
在本申请实施例中，还可以将裂隙的形状等效为长方形或其他形状，其等效面积获取原理与将裂隙的形状等效为平行四边形时的等效面积获取原理相同，在此不再赘述。
本申请通过使导热模型满足在孔的区域内以及在孔的影响范围内的积分条件，并使导热模型满足在孔的影响范围边界上以及孔的影响范围外的条件，从而建立每个孔影响下的导热模型，为获取动态冰的导热系数建立较为简单的前提条件。
在示例性实施例中，步骤S120可以包括子步骤S122。
子步骤S122：建立每个孔的影响范围在每个网格单元中的面积权重
Figure 45047DEST_PATH_IMAGE051
，
Figure 321307DEST_PATH_IMAGE052
，其中，
Figure 758105DEST_PATH_IMAGE053
为所述统计单元中的网格单元的编号，i为所述统计单元中的孔的编号，
Figure 447712DEST_PATH_IMAGE054
为第i个孔的影响范围在第
Figure 362579DEST_PATH_IMAGE053
个网格单元内的面积，N为所述统计单元中包含的孔的数量。
在本申请实施例中，如图6所示，图6示出了动态冰的其中一个统计单元，该统计单元被划分为5行、6列网格单元，第1个孔和第2个孔均为圆形，作为一种可选实施方式，
Figure 758925DEST_PATH_IMAGE053
可以表示第e行、第f列的网格单元，本申请对此不做限制，图6中的用阴影标记的网格单元为第2行、第3列，其编号可以为（2,3）。
以图6中的第1个孔的影响范围在编号（2,3）的网格单元中的面积权重为例进行说明。
第1个孔的影响范围在编号（2,3）的网格单元中的面积
Figure 491258DEST_PATH_IMAGE055
，即图6中用阴影标记的网格单元的边与
Figure 277948DEST_PATH_IMAGE023
的虚线围成的四分之一圆的面积，将统计单元中包含的N个孔的影响范围在编号（2,3）的网格单元中的面积进行累加，若某个孔的影响范围与该网格单元没有重叠部分，则令该孔的影响范围在该网格单元内的面积为0，然后用
Figure 121139DEST_PATH_IMAGE055
除以累加得到的面积，从而得到第1个孔的影响范围在编号（2,3）的网格单元中的面积权重
Figure 575254DEST_PATH_IMAGE056
，并按照该步骤得到每个孔的影响范围在每个网格单元中的面积权重
Figure 150592DEST_PATH_IMAGE051
。
步骤S130：根据所述每个孔影响下的导热模型和所述每个孔的影响范围在每个网格单元中的面积权重，得到所述统计单元中的每个网格单元的等效导热值。
在示例性实施例中，步骤S130可以包括子步骤S131。
子步骤S131：根据所述每个孔影响下的导热模型
Figure 549212DEST_PATH_IMAGE001
和所述每个孔的影响范围在每个网格单元中的面积权重
Figure 805881DEST_PATH_IMAGE051
，得到所述统计单元中的每个网格单元的等效导热值
Figure 239137DEST_PATH_IMAGE057
,
Figure 985376DEST_PATH_IMAGE058
，其中，i为所述统计单元中的孔的编号，N为所述统计单元中包含的孔的数量。
在本申请实施例中，导热模型
Figure 12238DEST_PATH_IMAGE001
在一个网格单元中各点的数值处处相等，作为一种可选实施方式，在获得统计单元中的每个网格单元的等效导热值
Figure 197231DEST_PATH_IMAGE059
时，将
Figure 360359DEST_PATH_IMAGE002
取为每个网格单元格心的坐标。
本申请通过将每个孔影响下的导热模型与该孔在网格单元中的面积权重相乘，得到每个孔对该网格单元的等效导热值的影响值，再将每个孔的影响值累加起来，得到该网格单元的等效导热值，如此对每个网格单元进行计算，从而得到每个网格单元的等效导热值，为获取动态冰的导热系数建立较为简单的前提条件。
步骤S140：根据所述每个网格单元的等效导热值得到所述统计单元的平均导热值。
在示例性实施例中，步骤S140可以包括子步骤S141至子步骤S144。
子步骤S141：在x轴方向采用等温边界条件，在y轴方向采用绝热边界条件，根据所述每个网格单元的等效导热值
Figure 605396DEST_PATH_IMAGE059
和二维稳态导热控制方程
Figure 650712DEST_PATH_IMAGE060
，通过有限元方法计算得到所述每个网格单元的温度
Figure 514763DEST_PATH_IMAGE061
。
如图7所示，若统计单元的左下角为原点，下侧边界向右为x轴正方向，左侧边界向上为y轴正方向，则在x轴方向采用等温边界条件，在y轴方向采用绝热边界条件，可以是：统计单元的上下两侧完全绝热，统计单元的左侧温度为
Figure 922611DEST_PATH_IMAGE062
，右侧温度为
Figure 948336DEST_PATH_IMAGE063
，且
Figure 820563DEST_PATH_IMAGE064
。
在本申请实施例中，每个网格单元的温度
Figure 285042DEST_PATH_IMAGE061
在一个网格单元中各点的数值处处相等。
子步骤S142：根据所述每个网格单元的温度
Figure 422763DEST_PATH_IMAGE061
和傅里叶稳态导热定律
Figure 9602DEST_PATH_IMAGE065
得到所述每个网格单元边界上的热流
Figure 967193DEST_PATH_IMAGE066
。
在本申请实施例中，在x轴方向，每个网格单元的左边界上有热流流出，每个网格的右边界上有热流流入。
子步骤S143：将所述统计单元热流计算侧的所有网格单元在该侧边界上的热流相加，得到所述统计单元的总热流
Figure 563260DEST_PATH_IMAGE067
，其中，所述热流计算侧为所述统计单元在x轴方向的两侧中的任意一侧。
如图8所示，以统计单元被划分为5行、6列的网格单元，统计单元的上下两侧完全绝热，统计单元的左两侧温差恒定为例进行说明，热流计算侧为统计单元的左侧和右侧中的任意一侧。若热流计算侧为统计单元的左侧，则统计单元热流计算侧的所有网格单元在该侧边界上的热流包括：第1列的所有网格单元的左边界上的热流，或第6列所有网格单元的右边界上的热流。
子步骤S144：根据所述统计单元的总热流
Figure 352224DEST_PATH_IMAGE067
，得到所述统计单元的平均导热值
Figure 985331DEST_PATH_IMAGE068
，其中，
Figure 554853DEST_PATH_IMAGE069
为所述统计单元的在x轴方向的两侧的长度，
Figure 564397DEST_PATH_IMAGE070
为所述统计单元在x轴方向的两侧的温差。
步骤S150：根据多个所述统计单元的平均导热值得到所述动态冰的导热系数。
在示例性实施例中，步骤S150可以包括子步骤S151。
子步骤S151：将多个所述统计单元的平均导热值的收敛值或统计数确定为所述动态冰的导热系数。
具体地，可以在多个统计单元的平均导热值收敛时，将多个统计单元的平均导热值的收敛值确定为动态冰的导热系数，也可以在得到多个统计单元的平均导热值后，将多个统计单元的平均导热值的平均值、标准差等统计数中的任意一种确定为动态冰的导热系数。
本申请通过建立每个孔影响下的导热模型，并建立每个孔的影响范围在每个网格单元中的面积权重，根据每个孔影响下的导热模型和每个孔的影响范围在每个网格单元中的面积权重，得到统计单元中的每个网格单元的等效导热值，然后再根据每个网格单元的等效导热值得到统计单元的平均导热值，并根据多个统计单元的平均导热值得到动态冰的导热系数，即使在孔隙大小相差较大时，也可以根据导热模型和面积权重算出每个网格单元的等效导热值，无需划分更多的网格，从而使获取动态冰的导热系数的过程更加简单。
最后应说明的是：以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不驱使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。

Claims (9)

1.一种动态冰的导热系数获取方法，其特征在于，所述方法包括：
S110.将所述动态冰划分为多个统计单元，将每个统计单元划分为多个网格单元，且划分方式不随孔的大小和分布而变化；
S120.若所述统计单元中包含至少一个孔，则根据所述统计单元中每个孔的影响范围，建立每个孔影响下的导热模型，并建立每个孔的影响范围在每个网格单元中的面积权重，其中，每个孔影响下的导热模型为
Figure 320639DEST_PATH_IMAGE001
，i为所述统计单元中的孔的编号，
Figure 48424DEST_PATH_IMAGE002
为所述统计单元中的坐标点，且
Figure 333912DEST_PATH_IMAGE001
同时满足下述条件：
对任意
Figure 867661DEST_PATH_IMAGE003
，
Figure 922205DEST_PATH_IMAGE004
，其中，
Figure 86470DEST_PATH_IMAGE005
为第i个孔所在的区域，
Figure 859254DEST_PATH_IMAGE006
为空气的导热系数，
Figure 931115DEST_PATH_IMAGE007
为所述第i个孔的面积，
对任意
Figure 105745DEST_PATH_IMAGE008
，
Figure 175332DEST_PATH_IMAGE009
，其中，
Figure 435412DEST_PATH_IMAGE010
为所述第i个孔的影响范围所在的区域，
Figure 310964DEST_PATH_IMAGE011
为第i个孔的模拟导热值，
Figure 340100DEST_PATH_IMAGE012
为所述第i个孔的影响范围的面积，
对任意
Figure 846167DEST_PATH_IMAGE013
，
Figure 327964DEST_PATH_IMAGE014
，
Figure 7207DEST_PATH_IMAGE015
为所述第i个孔的影响范围的边界，
Figure 890850DEST_PATH_IMAGE016
为静态冰的导热系数，
对任意
Figure 364556DEST_PATH_IMAGE017
，
Figure 536912DEST_PATH_IMAGE018
；
S130.根据所述每个孔影响下的导热模型和所述每个孔的影响范围在每个网格单元中的面积权重，得到所述统计单元中的每个网格单元的等效导热值；
S140.根据所述每个网格单元的等效导热值得到所述统计单元的平均导热值；
S150.根据多个所述统计单元的平均导热值得到所述动态冰的导热系数。
2.根据权利要求1所述的动态冰的导热系数获取方法，其特征在于，所述步骤S120包括：
建立每个孔的影响范围在每个网格单元中的面积权重
Figure 488687DEST_PATH_IMAGE019
，
Figure 757995DEST_PATH_IMAGE020
，其中，
Figure 402603DEST_PATH_IMAGE021
为所述统计单元中的网格单元的编号，i为所述统计单元中的孔的编号，
Figure 62254DEST_PATH_IMAGE022
为第i个孔的影响范围在第
Figure 817720DEST_PATH_IMAGE021
个网格单元内的面积，N为所述统计单元中包含的孔的数量。
3.根据权利要求1至2任一所述的动态冰的导热系数获取方法，其特征在于，所述步骤S130包括：
根据所述每个孔影响下的导热模型
Figure 675955DEST_PATH_IMAGE001
和所述每个孔的影响范围在每个网格单元中的面积权重
Figure 491464DEST_PATH_IMAGE019
，得到所述统计单元中的每个网格单元的等效导热值
Figure 638412DEST_PATH_IMAGE023
,
Figure 197569DEST_PATH_IMAGE024
，其中，i为所述统计单元中的孔的编号，N为所述统计单元中包含的孔的数量。
4.根据权利要求3所述的动态冰的导热系数获取方法，其特征在于，所述统计单元中的至少一个孔包括气泡和/或裂隙。
5.根据权利要求4所述的动态冰的导热系数获取方法，其特征在于，若第i个孔为气泡，则所述
Figure 175889DEST_PATH_IMAGE007
为第i个气泡的等效面积，且
Figure 896721DEST_PATH_IMAGE025
，其中，
Figure 530964DEST_PATH_IMAGE026
为所述第i个气泡的等效半径，所述
Figure 628233DEST_PATH_IMAGE027
为所述第i个气泡的影响范围的面积，且
Figure 461060DEST_PATH_IMAGE028
，其中，
Figure 618372DEST_PATH_IMAGE029
为所述第i个气泡的影响范围的半径；
若所述第i个孔为裂隙，则所述
Figure 536650DEST_PATH_IMAGE007
为第i个裂隙的等效面积，且
Figure 640872DEST_PATH_IMAGE030
，其中，
Figure 62626DEST_PATH_IMAGE031
为所述第i个裂隙的等效底边，
Figure 656418DEST_PATH_IMAGE032
为所述第i个裂隙的等效高，所述
Figure 796413DEST_PATH_IMAGE027
为所述第i个裂隙的影响范围的面积，且
Figure 704326DEST_PATH_IMAGE033
，其中，
Figure 246165DEST_PATH_IMAGE034
为所述第i个裂隙的影响范围的底边，
Figure 10859DEST_PATH_IMAGE035
为所述第i个裂隙的影响范围的高。
6.根据权利要求3所述的动态冰的导热系数获取方法，其特征在于，则所述步骤S140包括：
在x轴方向采用等温边界条件，在y轴方向采用绝热边界条件，根据所述每个网格单元的等效导热值
Figure 638150DEST_PATH_IMAGE036
和二维稳态导热控制方程
Figure 84174DEST_PATH_IMAGE037
，通过有限元方法计算得到所述每个网格单元的温度
Figure 480521DEST_PATH_IMAGE038
；
根据所述每个网格单元的温度
Figure 416116DEST_PATH_IMAGE038
和傅里叶稳态导热定律
Figure 530702DEST_PATH_IMAGE039
得到所述每个网格单元边界上的热流
Figure 514839DEST_PATH_IMAGE040
；
将所述统计单元热流计算侧的所有网格单元在该侧边界上的热流相加，得到所述统计单元的总热流
Figure 31271DEST_PATH_IMAGE041
，其中，所述热流计算侧为所述统计单元在x轴方向的两侧中的任意一侧；
根据所述统计单元的总热流
Figure 137767DEST_PATH_IMAGE042
，得到所述统计单元的平均导热值
Figure 739650DEST_PATH_IMAGE043
，其中，
Figure 527477DEST_PATH_IMAGE044
为所述统计单元的在x轴方向的两侧的长度，
Figure 898416DEST_PATH_IMAGE045
为所述统计单元在x轴方向的两侧的温差。
7.根据权利要求1所述的动态冰的导热系数获取方法，其特征在于，则所述步骤S150包括：
将多个所述统计单元的平均导热值的收敛值或统计数确定为所述动态冰的导热系数。
8.根据权利要求1所述的动态冰的导热系数获取方法，其特征在于，则所述步骤S120还包括：
若所述统计单元中不包含孔，则将静态冰的导热系数作为所述统计单元的平均导热值。
9.根据权利要求1所述的动态冰的导热系数获取方法，其特征在于，则所述步骤S110包括：
将所述动态冰划分为多个特征长度为
Figure 175813DEST_PATH_IMAGE046
的所述统计单元，且
Figure 999413DEST_PATH_IMAGE047
，其中，
Figure 590931DEST_PATH_IMAGE048
为所述动态冰的最小微观结构特征长度。
CN202210807281.6A 2022-07-11 2022-07-11 一种动态冰的导热系数获取方法 Active CN114858849B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210807281.6A CN114858849B (zh) 2022-07-11 2022-07-11 一种动态冰的导热系数获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210807281.6A CN114858849B (zh) 2022-07-11 2022-07-11 一种动态冰的导热系数获取方法

Publications (2)

Publication Number Publication Date
CN114858849A CN114858849A (zh) 2022-08-05
CN114858849B true CN114858849B (zh) 2022-10-04

Family

ID=82625663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210807281.6A Active CN114858849B (zh) 2022-07-11 2022-07-11 一种动态冰的导热系数获取方法

Country Status (1)

Country Link
CN (1) CN114858849B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115544818B (zh) * 2022-12-05 2023-02-28 中国空气动力研究与发展中心低速空气动力研究所 多层异质薄壁固体导热计算的网格划分及导热计算方法
CN116187140B (zh) * 2023-03-09 2023-10-13 成都流体动力创新中心 动态冰弹性模量预测方法
CN116306127B (zh) * 2023-03-09 2024-02-02 成都流体动力创新中心 动态冰微结构建模方法
CN115901841B (zh) * 2023-03-13 2023-05-30 中国电建集团昆明勘测设计研究院有限公司 水库拱桥形冰盖形成模拟及整体形态演变测量方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6238086B1 (en) * 1998-02-02 2001-05-29 Nec Corporation Method of calculating thermal resistance in semiconductor package accommodating semiconductor chip within a case which can be applied to calculation for semiconductor package with radiation fins
CN102162802A (zh) * 2011-01-11 2011-08-24 河海大学 全级配混凝土导热系数预测模型的构建方法
CN102313759A (zh) * 2011-07-26 2012-01-11 河海大学 一种硬化砂浆导热系数的预测模型的构建方法
CN106596630A (zh) * 2016-12-20 2017-04-26 东北电力大学 一种颗粒污垢等效导热系数的计算方法
CN111027236A (zh) * 2019-10-30 2020-04-17 武汉大学 多相复合材料的等效导热系数的细观尺度研究方法
CN113899784A (zh) * 2021-12-08 2022-01-07 中国空气动力研究与发展中心低速空气动力研究所 一种大型结冰风洞中复杂结冰云雾水含量测量方法
CN113901591A (zh) * 2021-12-08 2022-01-07 中国空气动力研究与发展中心低速空气动力研究所 热线单元干功率计算方法及基于此的液态水含量计算方法
CN114441590A (zh) * 2021-12-22 2022-05-06 中国航天空气动力技术研究院 一种针对梯度防热材料传热及力学特性的确定方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101699275B (zh) * 2009-10-30 2011-06-22 东华大学 一种高蓬松易变形材料导热性能的检测装置及检测方法
CN105181734B (zh) * 2015-09-25 2017-10-24 西南交通大学 一种形状记忆合金热机械疲劳实验装置
US10557810B2 (en) * 2016-10-06 2020-02-11 The Regents Of The University Of California Ice nucleating particle spectrometer
CN111443106B (zh) * 2020-05-15 2023-02-21 中南大学 非均质材料等效导热系数测试方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6238086B1 (en) * 1998-02-02 2001-05-29 Nec Corporation Method of calculating thermal resistance in semiconductor package accommodating semiconductor chip within a case which can be applied to calculation for semiconductor package with radiation fins
CN102162802A (zh) * 2011-01-11 2011-08-24 河海大学 全级配混凝土导热系数预测模型的构建方法
CN102313759A (zh) * 2011-07-26 2012-01-11 河海大学 一种硬化砂浆导热系数的预测模型的构建方法
CN106596630A (zh) * 2016-12-20 2017-04-26 东北电力大学 一种颗粒污垢等效导热系数的计算方法
CN111027236A (zh) * 2019-10-30 2020-04-17 武汉大学 多相复合材料的等效导热系数的细观尺度研究方法
CN113899784A (zh) * 2021-12-08 2022-01-07 中国空气动力研究与发展中心低速空气动力研究所 一种大型结冰风洞中复杂结冰云雾水含量测量方法
CN113901591A (zh) * 2021-12-08 2022-01-07 中国空气动力研究与发展中心低速空气动力研究所 热线单元干功率计算方法及基于此的液态水含量计算方法
CN114441590A (zh) * 2021-12-22 2022-05-06 中国航天空气动力技术研究院 一种针对梯度防热材料传热及力学特性的确定方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Experimental Thermal and Fluid;Yi Xian;《experimental thermal and fluid science》;20211231;第124卷;1-3 *
土热传导系数及模型的研究现状和展望;张楠;《岩土力学》;20160630;第37卷(第6期);1550-1562 *

Also Published As

Publication number Publication date
CN114858849A (zh) 2022-08-05

Similar Documents

Publication Publication Date Title
CN114858849B (zh) 一种动态冰的导热系数获取方法
Barros et al. Forcing symmetry exchanges and flow reversals in turbulent wakes
Lee et al. Observing and simulating wind-turbine wakes during the evening transition
Liu et al. Experimental investigation of ventilated partial cavitating flows with special emphasis on flow pattern regime and unsteady shedding behavior around an axisymmetric body at different angles of attack
CN110816885A (zh) 一种浮空器结冰特性数值仿真与试验验证系统
Ozcer et al. FENSAP-ICE: Numerical prediction of ice roughness evolution, and its effects on ice shapes
CN111695203B (zh) 一种反蜂群无人机气动布局设计和性能评估方法
Jiang et al. Model for ice wet growth on composite insulator and its experimental validation
US20240135057A1 (en) Method for numerical reconstruction and heat transfer characteristics evaluation of microstructure of thermal barrier coatings containing microcracks
KR20140115835A (ko) 공기 흡입구의 표면 결빙 시뮬레이션 방법
Wang et al. A numerical study on the ventilation coefficients of falling lobed hailstones
CN111310381A (zh) 一种三维水滴收集系数计算方法
CN113468655A (zh) 基于飞行试验数据的高速边界层转捩判据标定方法和系统
Kim et al. Estimation of roughness parameters within sparse urban-like obstacle arrays
Sirui et al. A review of icing prediction techniques for four typical surfaces in low-temperature natural environments
Liu et al. An experimental investigation on the convective heat transfer process over an ice roughened airfoil
Nath et al. Parametric investigation of aerodynamic performance degradation due to icing on a symmetrical airfoil
Steiner et al. Ice roughness and its impact on the ice accretion process
Cao et al. Insight into rime ice accretion on an aircraft wing and corresponding effects on aerodynamic performance
CN109492322A (zh) 一种储煤筒仓内部煤体自燃位置预测方法
Tabatabaei et al. Numerical study of aerodynamic characteristics of a symmetric NACA section with simulated ice shapes
Zhu et al. An Adaptive Cartesian Method for Prediction of the Unsteady Process of Aircraft Ice Accretion
Tabatabaei et al. Time‐Dependent Effects of Glaze Ice on the Aerodynamic Characteristics of an Airfoil
Szilder et al. Numerical Simulation of In-Flight Icing via a Particle-Based Morphogenetic Method
Ferdous et al. Prediction of ice accretion and aerodynamic performance analysis of NACA 2412 aerofoil

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