CN107833275A - 一种动态结冰微观结构三维建模方法 - Google Patents

一种动态结冰微观结构三维建模方法 Download PDF

Info

Publication number
CN107833275A
CN107833275A CN201711214269.XA CN201711214269A CN107833275A CN 107833275 A CN107833275 A CN 107833275A CN 201711214269 A CN201711214269 A CN 201711214269A CN 107833275 A CN107833275 A CN 107833275A
Authority
CN
China
Prior art keywords
dimensional
dynamic ice
hole
microstructure
modeling method
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
CN201711214269.XA
Other languages
English (en)
Other versions
CN107833275B (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.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational 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 Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN201711214269.XA priority Critical patent/CN107833275B/zh
Publication of CN107833275A publication Critical patent/CN107833275A/zh
Application granted granted Critical
Publication of CN107833275B publication Critical patent/CN107833275B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种动态结冰微观结构三维建模方法，提出以0‑1三维矩阵的形式表达结冰内部微观结构，其中0表示的是结冰，1表示的是结冰内的气泡孔隙。基于孔隙为球形的假设，以均匀分布的方式随机生成孔隙的球心坐标，并以特定分布的方式随机生成孔隙的直径。最终结合球心位置和相应的直径信息确定三维矩阵内各坐标点处的取值，从而获得结冰三维微观结构。这一结冰三维微观结构建模方法中的相关依据源于结冰实验数据，具有较高的可信度。同时，该方法为结冰的定量分析提供必要的支撑。

Description

一种动态结冰微观结构三维建模方法
技术领域
本发明涉及动态结冰微观特征分析技术领域，具体涉及一种动态结冰微观结构的三维建模方法。
背景技术
当过冷水滴撞击低温基底，在满足一定条件后会发生冻结，随着夹杂过冷水气流的不断撞击，基底表面形成越来越厚的结冰。与传统结冰较为不同的是，这种结冰具有典型的动态过程，在微观上，表现为水滴的不断冻结累积，并且相互间形成孔隙。这一微观结构的孔隙特性与密度、波的传播速度、导热系数等密切相关，直接影响结冰冰形计算、结冰探测和防除冰等。然而，目前缺少对动态结冰微观结构的定量刻画，使得对结冰相关问题的精细化研究缺少必要的理论依据。针对这一问题，基于统计的相关知识，提出了结冰微观结构的三维建模方法。
发明内容
本发明所要解决的技术问题是提供一种动态结冰微观结构三维建模的方法。该方法基于结冰内部气泡孔隙为球形的假设，将动态结冰微观结构抽象为三维矩阵，并结合孔隙直径和数量的分布规律，以及孔隙的随机性，确定其在结冰内的分布位置，最终达到以0-1（其中0表示结冰，1表示结冰内部的孔隙）矩阵的形式表征动态结冰三维微观结构的目的。
本发明的动态结冰微观结构三维建模方法，其特征在于包括以下步骤：
a. 给定动态结冰的三维区域。给定一个大小为的三维0-矩阵T，并确定矩阵大小与实际结冰尺寸的转换比例r，其中r的单位为个/长度单位，例如个/微米，个/毫米等。
b. 给定动态结冰内部孔隙数量N。
c. 以随机的形式确定每个孔隙的球心坐标。记球心为，。
d. 以随机的形式确定每个孔隙对应球的直径。记直径为，。
e. 生成动态结冰三维微观结构对应的0-1三维矩阵。对矩阵中的元素P进行判断取值，若存在整数，使得点P的与球心间的距离小于或者等于，则矩阵中的元素P的取值。
本发明的动态结冰三维微观结构建模方法中的孔隙数量、位置分布、半径分布等相关信息的确定规则均源于结冰实验数据，具有较高的可信度。该方法定量化提出了结冰的三维微观结构，为结冰的相关定量分析提供必要的支撑。
附图说明
图1为依据均匀分布，在区间[1, 19683000]上随机生成的2000个整数值随序号的分布图；
图2为依据均匀分布，在区间[0, 1]上随机生成的2000个数值随序号的分布图；
图3为生成的空隙直径分布图；
图4为所生成的动态结冰的三维微观结构图；
图5为截面x=135的孔隙分布示意图；
图6为截面y=135的孔隙分布示意图；
图7为截面z=135的孔隙分布示意图。
具体实施方式
下面结合附图和实施例进一步对本发明详细说明，以下实施例只是描述性的，非限定性的，不能以此限定本发明的保护范围。
实施例1
采用本发明方法，建立10×10×10 (mm3)大小的结冰微观结构三维模型。具体如下：
(1) 给定动态结冰的三维区域。给定一个大小为270×270×270的三维0-矩阵T，令矩阵大小与实际结冰尺寸的转换比例(个/毫米)，即该矩阵代表的真实区域大小为10×10×10(mm3)。
(2) 给定动态结冰内部孔隙数量。
(3) 以随机的形式确定每个孔隙的球心坐标。按照均匀分布方式，在区间[1,19683000]上生成2000个随机整数，其取值分布如图1所示，其中x轴表示2000个整数的序号，y轴表示其对应的整数值；以运算确定孔隙球心的坐标，其中函数int()表示取整运算。
(4) 确定每个孔隙对应球的直径。孔隙直径的生成基于的分布函数，其中。具体为：按照均匀分布规则，在区间[0,1]上生成2000个随机数，其图像如图2所示，并以此作为分布函数逆函数的自变量值，求出因变量值，并以作为孔隙的直径，其中。所生成孔隙的直径分布如图3所示，其中x轴表示直径，y轴表示直径小于x的孔隙数量占比。
(5) 生成动态结冰三维微观结构对应的0-1三维矩阵。按照的顺序，依次判断球心周围矩形区域(区域内点的坐标满足，y和z的取值类似)内的点，若其与球心间的距离小于或者等于，则矩阵在该坐标点处的取值，否则不进行任何操作。图4中给出了所生成的动态结冰的三维微观结构图，图5，图6，图7分别是截面x=135，y=135，z=135的示意图。

Claims (4)

1.一种动态结冰微观结构三维建模方法，其特征在于包括以下步骤：
(1) 给定动态结冰的三维区域；
给定一个大小为的三维0-矩阵T，并确定矩阵大小与实际结冰尺寸的转换比例r；
(2) 给定动态结冰内部孔隙数量N；
(3) 以随机的形式确定每个孔隙的球心，其中；
(4) 基于某分布确定每个孔隙对应球的直径，其中；
(5) 生成动态结冰三维微观结构对应的0-1三维矩阵；
对矩阵中的元素P进行判断取值，若存在整数，使得点P与球心间的距离小于或者等于，则矩阵中的元素P的取值。
2.根据权利要求1所述的动态结冰微观结构三维建模方法，其特征在于：步骤(3)中所述每个孔隙的球心坐标的确定方式为按照均匀分布的规则在 ]区间上生成N个随机整数；以运算确定孔隙的球心坐标，其中函数int()表示取整运算。
3.根据权利要求1所述的动态结冰微观结构三维建模方法，其特征在于：步骤(4)中所述分布为，其中直径为任意正数。
4.根据权利要求3所述的动态结冰微观结构三维建模方法，其特征在于：所述。
CN201711214269.XA 2017-11-28 2017-11-28 一种动态结冰微观结构三维建模方法 Active CN107833275B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711214269.XA CN107833275B (zh) 2017-11-28 2017-11-28 一种动态结冰微观结构三维建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711214269.XA CN107833275B (zh) 2017-11-28 2017-11-28 一种动态结冰微观结构三维建模方法

Publications (2)

Publication Number Publication Date
CN107833275A true CN107833275A (zh) 2018-03-23
CN107833275B CN107833275B (zh) 2020-12-08

Family

ID=61646043

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711214269.XA Active CN107833275B (zh) 2017-11-28 2017-11-28 一种动态结冰微观结构三维建模方法

Country Status (1)

Country Link
CN (1) CN107833275B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109118534A (zh) * 2018-07-13 2019-01-01 同济大学 基于机器视觉图像实时获取冰风洞内模型表面上结冰冰型详细信息的方法
CN114076727A (zh) * 2022-01-10 2022-02-22 中国空气动力研究与发展中心低速空气动力研究所 一种基于电阻率的冰的孔隙率测量方法
CN115374599A (zh) * 2022-07-13 2022-11-22 河海大学 一种基于随机分布的3d纤维状孔隙冰数值建模方法
CN116187140A (zh) * 2023-03-09 2023-05-30 成都流体动力创新中心 动态冰弹性模量预测方法
CN116306127A (zh) * 2023-03-09 2023-06-23 成都流体动力创新中心 动态冰微结构建模方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070110873A1 (en) * 2003-07-07 2007-05-17 Nestec, S.A. Aerated frozen suspensions with adjusted creaminess and scoop ability based on stress-controlled generation of superfine microstructures
CN102682144A (zh) * 2011-11-30 2012-09-19 天津空中代码工程应用软件开发有限公司 直升机旋翼飞行结冰的数值模拟方法
WO2014120961A1 (en) * 2013-01-30 2014-08-07 University Of Pittsburgh - Of The Commonwealth System Of Higher Education Compositions for prevention of ice build-up
CN104298886A (zh) * 2014-10-20 2015-01-21 上海电机学院 航空发动机旋转部件的结冰3-d数值模拟方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070110873A1 (en) * 2003-07-07 2007-05-17 Nestec, S.A. Aerated frozen suspensions with adjusted creaminess and scoop ability based on stress-controlled generation of superfine microstructures
CN102682144A (zh) * 2011-11-30 2012-09-19 天津空中代码工程应用软件开发有限公司 直升机旋翼飞行结冰的数值模拟方法
WO2014120961A1 (en) * 2013-01-30 2014-08-07 University Of Pittsburgh - Of The Commonwealth System Of Higher Education Compositions for prevention of ice build-up
CN104298886A (zh) * 2014-10-20 2015-01-21 上海电机学院 航空发动机旋转部件的结冰3-d数值模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JOSHUA BLAKE .ECT: "Simulating the Freezing of Supercooled Water Droplets Impacting a Cooled Substrate", 《52ND AEROSPACE SCIENCES MEETING》 *
李伟斌 等: "动态结冰微观孔隙结构定量分析", 《航空学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109118534A (zh) * 2018-07-13 2019-01-01 同济大学 基于机器视觉图像实时获取冰风洞内模型表面上结冰冰型详细信息的方法
CN114076727A (zh) * 2022-01-10 2022-02-22 中国空气动力研究与发展中心低速空气动力研究所 一种基于电阻率的冰的孔隙率测量方法
CN114076727B (zh) * 2022-01-10 2022-05-13 中国空气动力研究与发展中心低速空气动力研究所 一种基于电阻率的冰的孔隙率测量方法
CN115374599A (zh) * 2022-07-13 2022-11-22 河海大学 一种基于随机分布的3d纤维状孔隙冰数值建模方法
CN115374599B (zh) * 2022-07-13 2023-06-16 河海大学 一种基于随机分布的3d纤维状孔隙冰数值建模方法
CN116187140A (zh) * 2023-03-09 2023-05-30 成都流体动力创新中心 动态冰弹性模量预测方法
CN116306127A (zh) * 2023-03-09 2023-06-23 成都流体动力创新中心 动态冰微结构建模方法
CN116187140B (zh) * 2023-03-09 2023-10-13 成都流体动力创新中心 动态冰弹性模量预测方法
CN116306127B (zh) * 2023-03-09 2024-02-02 成都流体动力创新中心 动态冰微结构建模方法

Also Published As

Publication number Publication date
CN107833275B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN107833275A (zh) 一种动态结冰微观结构三维建模方法
Carrier et al. Evaporation of water: evaporation rate and collective effects
Padoan et al. Supersonic turbulence in the Perseus molecular cloud
CN101706967B (zh) 一种真实感云层的综合模拟方法
CN107423498A (zh) 一种高致密度离散颗粒多相体系的建模方法
CN106021828A (zh) 一种基于格子-玻尔兹曼模型的流体模拟方法
CN112613243B (zh) 一种流体力学模拟的方法、装置及计算机可读存储介质
CN109117861A (zh) 一种顾及空间位置的点集多层次聚集性分析方法
CN110457771A (zh) 一种基于高程偏差传递的dem水流方向计算方法
Wang et al. User recruitment system for efficient photo collection in mobile crowdsensing
EP3179390A1 (en) Method and apparatus for modeling movement of air bubble based on fluid particles
CN107590327A (zh) 一种水泥基材料颗粒三维随机投放的分层分区方法
CN112364435B (zh) 多次抛撒箔条云气动融合轨道再入实时跟踪统计模拟方法
CN107206691A (zh) 用于设置增材制造工艺的三维物体的打印属性的方法
CN112233231B (zh) 一种基于云计算的城市三维实景漫游方法及系统
CN109117519B (zh) 一种油箱结构瞬态温度场的模拟方法
CN108920872A (zh) 针对dsmc方法的bcp粒子定位实现方法及系统
CN115098827A (zh) 一种积雪环境下的红外场景数值计算方法
CN114564901B (zh) 一种结合随机函数的汽车涂层抗石击性仿真评价方法
CN113158520A (zh) 一种用于冷冻靶系统中的燃料冰层界面追踪模拟方法
CN107704653B (zh) 一种水泥基材料三维颗粒分组均匀化投放方法
CN107239559A (zh) 一种基于矢量网格进行空间动目标位置计算的方法
CN106021643B (zh) 一种随机缩松缺陷对铸件力学性能影响的数值分析方法
CN114330080A (zh) 一种飞行器面对称跨流域流场的预测方法
CN104168579B (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