CN111241761A - 一种风力机叶片防冰涂料涂抹位置的确定方法 - Google Patents

一种风力机叶片防冰涂料涂抹位置的确定方法 Download PDF

Info

Publication number
CN111241761A
CN111241761A CN202010080304.9A CN202010080304A CN111241761A CN 111241761 A CN111241761 A CN 111241761A CN 202010080304 A CN202010080304 A CN 202010080304A CN 111241761 A CN111241761 A CN 111241761A
Authority
CN
China
Prior art keywords
wall
temperature
mass
film
liquid
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
CN202010080304.9A
Other languages
English (en)
Other versions
CN111241761B (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.)
Ocean University of China
Original Assignee
Ocean University of China
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 Ocean University of China filed Critical Ocean University of China
Priority to CN202010080304.9A priority Critical patent/CN111241761B/zh
Publication of CN111241761A publication Critical patent/CN111241761A/zh
Application granted granted Critical
Publication of CN111241761B publication Critical patent/CN111241761B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种风力机叶片防冰涂料涂抹位置的确定方法，包括如下步骤：（1）组分输运模型：（2）膜凝结模型：（3）欧拉壁面模型：（4）水汽会在风力机叶片的前缘和尾缘位置凝结，故应在前缘和尾缘位置涂抹防冰涂料。通过本发明所公开的方法，确定不同条件下叶片表面上的结冰位置，在相应的结冰位置涂抹防冰涂料，既达到减轻结冰程度的目的，又节省涂料，降低成本。

Description

一种风力机叶片防冰涂料涂抹位置的确定方法
技术领域
本发明属于风力发电机组领域，特别涉及该领域中的一种风力机叶片防冰涂料涂抹位置的确定方法。
背景技术
现阶段风力机的防冰涂料价格昂贵，如果在风力机叶片全表面涂抹防防冰涂料，成本非常高。
发明内容
本发明所要解决的技术问题就是提供一种风力机叶片防冰涂料涂抹位置的确定方法。
本发明采用如下技术方案：
一种风力机叶片防冰涂料涂抹位置的确定方法，其改进之处在于，包括如下步骤：
(1)组分输运模型：
湿空气为干空气和水蒸气的混合气体，组分输运模型将混合气体定义为一种二元混合气体，具有两者的物理属性，求解过程作为单相处理，湿空气密度按照不可压缩理想气体方程计算，湿空气的控制方程如下所述，质量方程：
Figure BDA0002380059720000011
式中：ρm为混合物密度，kg/m3；U为流速；
动量方程：
Figure BDA0002380059720000012
式中：P为混合压力，pa；τ为切应力张量；B为体受力向量，
能量方程：
Figure BDA0002380059720000013
式中：
Figure BDA0002380059720000014
为净源项；q为净热通量；h为焓值；
组分质量：
Figure BDA0002380059720000015
式中：Yi为局部质量分数；
Figure BDA0002380059720000016
为扩散的质量；
Figure BDA0002380059720000017
式中：Sct为湍流施密特数；μt为湍流黏度，kg/(m·s)；Di,m为质量扩散系数；
(2)膜凝结模型：
模型采用BELL蒸汽壁面冷凝模型，水蒸气扩散通过传质边界层，遇冷壁面释放潜热形成液膜，液体膜达到壁温，蒸汽扩散和凝结在气液截面持续发生，冷凝率取决于蒸汽的扩散速率，反过来又取决于边界层的水蒸汽浓度差，在边界层的一端是大量混合流体，浓度取决于混合物的参数，边界层的另一端是凝结液体水膜和气体混合物界面，假定在液体膜—水蒸汽混合物界面存在局部平衡，这意味着在界面处的水蒸汽浓度等于在膜温度下的饱和值，因此，膜界面温度决定了整个边界层的浓度差，并反过来控制冷凝率，假设如下所述：
(21)由于水蒸汽通过组分边界层扩散导致冷凝，这个假设用于导出冷凝水体积数学表达式，
(22)液膜和壁面存在热平衡，液膜温度Tfilm等于壁面温度Twall，
(23)气液界面是饱和状态。即交界面水蒸气分压力等于壁面温度下的饱和蒸汽压力，结合(22)，壁面温度决定气液界面的蒸汽浓度，控制总的蒸汽扩散和通过边界层的浓度比，
(24)不考虑液膜运动带来的影响，由于流体剪切力和重力，液膜可能运动，但是运动不影响空气流动和凝结，
(25)分析仅限于膜状凝结，不考虑滴状冷凝，
(26)液膜的热阻忽略不计。
由于模型是基于以上假设建立的，凝结速率取决于水蒸气朝壁面的扩散速率，空气和水蒸气各组分的质量流量方程为：
Figure BDA0002380059720000021
Figure BDA0002380059720000022
式中：ω为质量分数；
Figure BDA0002380059720000023
为质量流量，kg/(m2·s)；υ为界面速度，m/s；D为质量扩散率，m2/s，
用户自定义函数中，假定壁面温度等于或小于表面水蒸气分压力对应的饱和温度时发生凝结，如果温度高于饱和温度，水蒸气质量分数等于壁面邻近单元值，为了满足气液交界面局部热力学平衡的假设，当温度低于或等于饱和温度时，代表气液交界面的水蒸气质量分数被分配一个值，使得水蒸汽的分压等于水在局部壁面温度下的饱和压力，
混合物的体积质量源项为：
Figure BDA0002380059720000024
式中：Acell wall为壁面单元面积；Vcell为计算单体体积，
水蒸气的体积质量源项为：
Figure BDA0002380059720000031
(3)欧拉壁面模型：
欧拉壁面膜模型，实体表面膜不需划分网格，以表面虚拟膜代替，不影响核心流动动量场，多数情况下，表面薄膜的分离、脱落、蒸发不影响核心流动场的流动，通常情况下，模拟此种类型薄膜需要极大的计算资源，尤其在多相工况下精确计算相间通量，需要极细的网格来模拟水膜，而欧拉壁面膜模型不需要划分水膜网格，同时可以预测收集效率、耦合离散相模型，满足模拟壁面收集液滴形成水膜过程的需求，将水分考虑为小液滴，液滴在表面汇集成水膜，计算不考虑相变，
质量守恒方程：
Figure BDA0002380059720000032
膜动量守恒方程：
Figure BDA0002380059720000033
式中：h为膜高度；ρ1液体密度；
Figure BDA0002380059720000034
为表面梯度算子；
Figure BDA0002380059720000035
平均膜速度；
Figure BDA0002380059720000036
为每单位壁面积质量源；g为重力加速度，
膜能量守恒方程：
Figure BDA0002380059720000037
式中：Ts为界面温度；Tf为膜平均温度；Vf为膜速度；kf为膜热导；Tw为壁面温度；假设膜的温度变化为线性分段变化，下半部分从壁面温度到膜平均温度Tf，上半部分由平均温度Tf到界面温度Ts；
Figure BDA0002380059720000038
为液相撞击壁面的源项；
Figure BDA0002380059720000039
为蒸发或凝结质量流率，
当离散粒子或液滴撞击壁面，会被液膜吸收，吸收后，其质量和动量增加到质量方程和动量方程的源项中，质量源项和动量源项分别为：
Figure BDA0002380059720000041
式中
Figure BDA0002380059720000042
为颗粒或液滴撞击表面的质量流率，
Figure BDA0002380059720000043
式中
Figure BDA0002380059720000044
分别为颗粒和液膜的速度，
质量源项：
Figure BDA0002380059720000045
式中：αd为第二相体积；ρd为第二相密度；vdn为垂直于壁面的相速度；A为壁面表面积,动量源项
Figure BDA0002380059720000046
Figure BDA0002380059720000047
为第二相速度向量；
(4)水汽会在风力机叶片的前缘和尾缘位置凝结，故应在前缘和尾缘位置涂抹防冰涂料。
进一步的，使用的模型包括但不限于水蒸气凝结模型、欧拉壁面膜模型和颗粒磨损模型。
本发明的有益效果是：
通过本发明所公开的方法，确定不同条件下叶片表面上的结冰位置，在相应的结冰位置涂抹防冰涂料，既达到减轻结冰程度的目的，又节省涂料，降低成本。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本实施例的原理是通过CFD方法，采用数值模拟技术，模拟全尺寸叶片表面液滴撞击叶片表面的位置及水分汇集的位置，在该位置涂抹防冰涂料。
实施例1，本实施例公开了一种风力机叶片防冰涂料涂抹位置的确定方法，包括如下步骤：
(1)组分输运模型：
湿空气为干空气和水蒸气的混合气体，组分输运模型将混合气体定义为一种二元混合气体，具有两者的物理属性，求解过程作为单相处理，湿空气密度按照不可压缩理想气体方程计算，湿空气的控制方程如下所述，质量方程：
Figure BDA0002380059720000048
式中：ρm为混合物密度，kg/m3；U为流速；
动量方程：
Figure BDA0002380059720000049
式中：P为混合压力，pa；τ为切应力张量；B为体受力向量，
能量方程：
Figure BDA0002380059720000051
式中：
Figure BDA0002380059720000052
为净源项；q为净热通量；h为焓值；
组分质量：
Figure BDA0002380059720000053
式中：Yi为局部质量分数；
Figure BDA0002380059720000054
为扩散的质量；
Figure BDA0002380059720000055
式中：Sct为湍流施密特数；μt为湍流黏度，kg/(m·s)；Di,m为质量扩散系数；
(2)膜凝结模型：
模型采用BELL蒸汽壁面冷凝模型，水蒸气扩散通过传质边界层，遇冷壁面释放潜热形成液膜，液体膜达到壁温，蒸汽扩散和凝结在气液截面持续发生，冷凝率取决于蒸汽的扩散速率，反过来又取决于边界层的水蒸汽浓度差，在边界层的一端是大量混合流体，浓度取决于混合物的参数，边界层的另一端是凝结液体水膜和气体混合物界面，假定在液体膜—水蒸汽混合物界面存在局部平衡，这意味着在界面处的水蒸汽浓度等于在膜温度下的饱和值，因此，膜界面温度决定了整个边界层的浓度差，并反过来控制冷凝率，假设如下所述：
(21)由于水蒸汽通过组分边界层扩散导致冷凝，这个假设用于导出冷凝水体积数学表达式，
(22)液膜和壁面存在热平衡，液膜温度Tfilm等于壁面温度Twall，
(23)气液界面是饱和状态。即交界面水蒸气分压力等于壁面温度下的饱和蒸汽压力，结合(22)，壁面温度决定气液界面的蒸汽浓度，控制总的蒸汽扩散和通过边界层的浓度比，
(24)不考虑液膜运动带来的影响，由于流体剪切力和重力，液膜可能运动，但是运动不影响空气流动和凝结，
(25)分析仅限于膜状凝结，不考虑滴状冷凝，
(26)液膜的热阻忽略不计。
由于模型是基于以上假设建立的，凝结速率取决于水蒸气朝壁面的扩散速率，空气和水蒸气各组分的质量流量方程为：
Figure BDA0002380059720000056
Figure BDA0002380059720000061
式中：ω为质量分数；
Figure BDA0002380059720000062
为质量流量，kg/(m2·s)；υ为界面速度，m/s；D为质量扩散率，m2/s，
用户自定义函数中，假定壁面温度等于或小于表面水蒸气分压力对应的饱和温度时发生凝结，如果温度高于饱和温度，水蒸气质量分数等于壁面邻近单元值，为了满足气液交界面局部热力学平衡的假设，当温度低于或等于饱和温度时，代表气液交界面的水蒸气质量分数被分配一个值，使得水蒸汽的分压等于水在局部壁面温度下的饱和压力，
混合物的体积质量源项为：
Figure BDA0002380059720000063
式中：Acell wall为壁面单元面积；Vcell为计算单体体积，
水蒸气的体积质量源项为：
Figure BDA0002380059720000064
(3)欧拉壁面模型：
欧拉壁面膜模型，实体表面膜不需划分网格，以表面虚拟膜代替，不影响核心流动动量场，多数情况下，表面薄膜的分离、脱落、蒸发不影响核心流动场的流动，通常情况下，模拟此种类型薄膜需要极大的计算资源，尤其在多相工况下精确计算相间通量，需要极细的网格来模拟水膜，而欧拉壁面膜模型不需要划分水膜网格，同时可以预测收集效率、耦合离散相模型，满足模拟壁面收集液滴形成水膜过程的需求，将水分考虑为小液滴，液滴在表面汇集成水膜，计算不考虑相变，
质量守恒方程：
Figure BDA0002380059720000065
膜动量守恒方程：
Figure BDA0002380059720000066
式中：h为膜高度；ρ1液体密度；
Figure BDA0002380059720000067
为表面梯度算子；
Figure BDA0002380059720000068
平均膜速度；
Figure BDA0002380059720000069
为每单位壁面积质量源；g为重力加速度，
膜能量守恒方程：
Figure BDA0002380059720000071
式中：Ts为界面温度；Tf为膜平均温度；Vf为膜速度；kf为膜热导；Tw为壁面温度；假设膜的温度变化为线性分段变化，下半部分从壁面温度到膜平均温度Tf，上半部分由平均温度Tf到界面温度Ts；
Figure BDA0002380059720000072
为液相撞击壁面的源项；
Figure BDA0002380059720000073
为蒸发或凝结质量流率，
当离散粒子或液滴撞击壁面，会被液膜吸收，吸收后，其质量和动量增加到质量方程和动量方程的源项中，质量源项和动量源项分别为：
Figure BDA0002380059720000074
式中
Figure BDA0002380059720000075
为颗粒或液滴撞击表面的质量流率，
Figure BDA0002380059720000076
式中
Figure BDA0002380059720000077
分别为颗粒和液膜的速度，
质量源项：
Figure BDA0002380059720000078
式中：αd为第二相体积；ρd为第二相密度；vdn为垂直于壁面的相速度；A为壁面表面积,动量源项
Figure BDA0002380059720000079
Figure BDA00023800597200000710
为第二相速度向量；
(4)水汽会在风力机叶片的前缘和尾缘位置凝结，故应在前缘和尾缘位置涂抹防冰涂料。
本实施例使用的模型包括但不限于水蒸气凝结模型、欧拉壁面膜模型和颗粒磨损模型。

Claims (2)

1.一种风力机叶片防冰涂料涂抹位置的确定方法，其特征在于，包括如下步骤：
(1)组分输运模型：
湿空气为干空气和水蒸气的混合气体，组分输运模型将混合气体定义为一种二元混合气体，具有两者的物理属性，求解过程作为单相处理，湿空气密度按照不可压缩理想气体方程计算，湿空气的控制方程如下所述，质量方程：
Figure FDA0002380059710000011
式中：ρm为混合物密度，kg/m3；U为流速；
动量方程：
Figure FDA0002380059710000012
式中：P为混合压力，pa；τ为切应力张量；B为体受力向量，
能量方程：
Figure FDA0002380059710000013
式中：
Figure FDA0002380059710000014
为净源项；q为净热通量；h为焓值；
组分质量：
Figure FDA0002380059710000015
式中：Yi为局部质量分数；
Figure FDA0002380059710000016
为扩散的质量；
Figure FDA0002380059710000017
式中：Sct为湍流施密特数；μt为湍流黏度，kg/(m·s)；Di,m为质量扩散系数；
(2)膜凝结模型：
模型采用BELL蒸汽壁面冷凝模型，水蒸气扩散通过传质边界层，遇冷壁面释放潜热形成液膜，液体膜达到壁温，蒸汽扩散和凝结在气液截面持续发生，冷凝率取决于蒸汽的扩散速率，反过来又取决于边界层的水蒸汽浓度差，在边界层的一端是大量混合流体，浓度取决于混合物的参数，边界层的另一端是凝结液体水膜和气体混合物界面，假定在液体膜—水蒸汽混合物界面存在局部平衡，这意味着在界面处的水蒸汽浓度等于在膜温度下的饱和值，因此，膜界面温度决定了整个边界层的浓度差，并反过来控制冷凝率，假设如下所述：
(21)由于水蒸汽通过组分边界层扩散导致冷凝，这个假设用于导出冷凝水体积数学表达式，
(22)液膜和壁面存在热平衡，液膜温度Tfilm等于壁面温度Twall，
(23)气液界面是饱和状态。即交界面水蒸气分压力等于壁面温度下的饱和蒸汽压力，结合(22)，壁面温度决定气液界面的蒸汽浓度，控制总的蒸汽扩散和通过边界层的浓度比，
(24)不考虑液膜运动带来的影响，由于流体剪切力和重力，液膜可能运动，但是运动不影响空气流动和凝结，
(25)分析仅限于膜状凝结，不考虑滴状冷凝，
(26)液膜的热阻忽略不计。
由于模型是基于以上假设建立的，凝结速率取决于水蒸气朝壁面的扩散速率，空气和水蒸气各组分的质量流量方程为：
Figure FDA0002380059710000021
Figure FDA0002380059710000022
式中：ω为质量分数；
Figure FDA0002380059710000023
为质量流量，kg/(m2·s)；υ为界面速度，m/s；D为质量扩散率，m2/s，
用户自定义函数中，假定壁面温度等于或小于表面水蒸气分压力对应的饱和温度时发生凝结，如果温度高于饱和温度，水蒸气质量分数等于壁面邻近单元值，为了满足气液交界面局部热力学平衡的假设，当温度低于或等于饱和温度时，代表气液交界面的水蒸气质量分数被分配一个值，使得水蒸汽的分压等于水在局部壁面温度下的饱和压力，
混合物的体积质量源项为：
Figure FDA0002380059710000024
式中：Acell wall为壁面单元面积；Vcell为计算单体体积，
水蒸气的体积质量源项为：
Figure FDA0002380059710000025
(3)欧拉壁面模型：
欧拉壁面膜模型，实体表面膜不需划分网格，以表面虚拟膜代替，不影响核心流动动量场，多数情况下，表面薄膜的分离、脱落、蒸发不影响核心流动场的流动，通常情况下，模拟此种类型薄膜需要极大的计算资源，尤其在多相工况下精确计算相间通量，需要极细的网格来模拟水膜，而欧拉壁面膜模型不需要划分水膜网格，同时可以预测收集效率、耦合离散相模型，满足模拟壁面收集液滴形成水膜过程的需求，将水分考虑为小液滴，液滴在表面汇集成水膜，计算不考虑相变，
质量守恒方程：
Figure FDA0002380059710000031
膜动量守恒方程：
Figure FDA0002380059710000032
式中：h为膜高度；ρ1液体密度；
Figure FDA0002380059710000033
为表面梯度算子；
Figure FDA0002380059710000034
平均膜速度；
Figure FDA0002380059710000035
为每单位壁面积质量源；g为重力加速度，
膜能量守恒方程：
Figure FDA0002380059710000036
式中：Ts为界面温度；Tf为膜平均温度；Vf为膜速度；kf为膜热导；Tw为壁面温度；假设膜的温度变化为线性分段变化，下半部分从壁面温度到膜平均温度Tf，上半部分由平均温度Tf到界面温度Ts；
Figure FDA0002380059710000037
为液相撞击壁面的源项；
Figure FDA0002380059710000038
为蒸发或凝结质量流率，
当离散粒子或液滴撞击壁面，会被液膜吸收，吸收后，其质量和动量增加到质量方程和动量方程的源项中，质量源项和动量源项分别为：
Figure FDA0002380059710000039
式中
Figure FDA00023800597100000310
为颗粒或液滴撞击表面的质量流率，
Figure FDA00023800597100000311
式中
Figure FDA00023800597100000312
分别为颗粒和液膜的速度，
质量源项：
Figure FDA00023800597100000313
式中：αd为第二相体积；ρd为第二相密度；vdn为垂直于壁面的相速度；A为壁面表面积,动量源项
Figure FDA00023800597100000314
Figure FDA00023800597100000315
为第二相速度向量；
(4)水汽会在风力机叶片的前缘和尾缘位置凝结，故应在前缘和尾缘位置涂抹防冰涂料。
2.根据权利要求1所述风力机叶片防冰涂料涂抹位置的确定方法，其特征在于：使用的模型包括但不限于水蒸气凝结模型、欧拉壁面膜模型和颗粒磨损模型。
CN202010080304.9A 2020-02-05 2020-02-05 一种风力机叶片防冰涂料涂抹位置的确定方法 Active CN111241761B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010080304.9A CN111241761B (zh) 2020-02-05 2020-02-05 一种风力机叶片防冰涂料涂抹位置的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010080304.9A CN111241761B (zh) 2020-02-05 2020-02-05 一种风力机叶片防冰涂料涂抹位置的确定方法

Publications (2)

Publication Number Publication Date
CN111241761A true CN111241761A (zh) 2020-06-05
CN111241761B CN111241761B (zh) 2022-10-28

Family

ID=70866320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010080304.9A Active CN111241761B (zh) 2020-02-05 2020-02-05 一种风力机叶片防冰涂料涂抹位置的确定方法

Country Status (1)

Country Link
CN (1) CN111241761B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110030723A1 (en) * 2008-06-24 2011-02-10 Fronius International Gmbh Apparatus and method for cleaning welding torches
CN102682145A (zh) * 2011-11-30 2012-09-19 天津空中代码工程应用软件开发有限公司 飞行结冰的数值模拟方法
CN107127190A (zh) * 2017-05-31 2017-09-05 中国空气动力研究与发展中心低速空气动力研究所 超疏水材料和振动除冰相耦合的主被动防冰装置及其方法
CN109751204A (zh) * 2019-02-18 2019-05-14 中国空气动力研究与发展中心低速空气动力研究所 一种风力机结冰数值模拟方法
CN110362918A (zh) * 2019-07-12 2019-10-22 西安交通大学 一种压水反应堆安全壳两侧冷凝与蒸发耦合计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110030723A1 (en) * 2008-06-24 2011-02-10 Fronius International Gmbh Apparatus and method for cleaning welding torches
CN102682145A (zh) * 2011-11-30 2012-09-19 天津空中代码工程应用软件开发有限公司 飞行结冰的数值模拟方法
CN107127190A (zh) * 2017-05-31 2017-09-05 中国空气动力研究与发展中心低速空气动力研究所 超疏水材料和振动除冰相耦合的主被动防冰装置及其方法
CN109751204A (zh) * 2019-02-18 2019-05-14 中国空气动力研究与发展中心低速空气动力研究所 一种风力机结冰数值模拟方法
CN110362918A (zh) * 2019-07-12 2019-10-22 西安交通大学 一种压水反应堆安全壳两侧冷凝与蒸发耦合计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李延 等: ""基于欧拉壁面液膜模型的三维热气"", 《北京航空航天大学学报》 *

Also Published As

Publication number Publication date
CN111241761B (zh) 2022-10-28

Similar Documents

Publication Publication Date Title
Fortin et al. Heat and mass transfer during ice accretion on aircraft wings with an improved roughness model
Dang et al. Numerical simulation of thermal performance for super large-scale wet cooling tower equipped with an axial fan
Cui et al. Study on the dynamic and thermal performances of a reversibly used cooling tower with upward spraying
Birbarah et al. External convective jumping-droplet condensation on a flat plate
Fortin et al. New roughness computation method and geometric accretion model for airfoil icing
Douvi et al. Aerodynamic performance investigation under the influence of heavy rain of a NACA 0012 airfoil for wind turbine applications
CN111967201A (zh) 一种基于数值仿真模型分析临界覆冰类型的方法
CN114896906B (zh) 一种考虑冰层和固体壁面中导热的积冰模拟方法
CN112985347A (zh) 一种飞机结冰表面粗糙度计算方法
Zhan et al. Numerical model of particle deposition on wet fin surfaces of heat exchanger under dehumidifying conditions
CN111241761B (zh) 一种风力机叶片防冰涂料涂抹位置的确定方法
Fortin et al. Prediction of 2D airfoil ice accretion by bisection method and by rivulets and beads modeling
CN114139393A (zh) 考虑水膜流动传热的部件电加热三维防冰数值模拟方法
Dimitrova et al. Influence of thermal effects on the wind field within the urban environment
CN110502780A (zh) 一种基于cfd-dem耦合模型的链篦机球团干燥速率仿真分析方法
Sivasamy et al. Numerical simulation of two‐dimensional laminar slot‐jet impingement flows confined by a parallel wall
Enayatollahi et al. Characterising the heat and mass transfer coefficients for a crossflow interaction of air and water
Douvi et al. Experimental and computational study of the effects of different liquid water content on the aerodynamic performance of a NACA 0012 airfoil at low Reynolds number
Qiu et al. Development and validation of a numerical model for frost growth based on nucleation theory
Comini et al. Latent and sensible heat transfer in air‐cooling applications
Yang et al. Numerical model of dusty-water bridge forming between neighbor fins of fin-and-tube heat exchangers
Javadi et al. Improving Power Plant Efficiency using Water Droplet Injection in Air Condensers
CN111709141A (zh) 一种基于Fluent的城市水体蒸发模拟方法
Gao et al. Calculation of natural ventilation in large enclosures
Yang et al. The analysis of the influence of packing and total pressure on cooling performance of the cooling tower

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