CN111241761B - 一种风力机叶片防冰涂料涂抹位置的确定方法 - Google Patents
一种风力机叶片防冰涂料涂抹位置的确定方法 Download PDFInfo
- Publication number
- CN111241761B CN111241761B CN202010080304.9A CN202010080304A CN111241761B CN 111241761 B CN111241761 B CN 111241761B CN 202010080304 A CN202010080304 A CN 202010080304A CN 111241761 B CN111241761 B CN 111241761B
- Authority
- CN
- China
- Prior art keywords
- temperature
- wall
- 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.)
- Active
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)组分输运模型:
湿空气为干空气和水蒸气的混合气体,组分输运模型将混合气体定义为一种二元混合气体,具有两者的物理属性,求解过程作为单相处理,湿空气密度按照不可压缩理想气体方程计算,湿空气的控制方程如下所述,质量方程:
式中:ρm为混合物密度,kg/m3;U为流速;
动量方程:
式中:P为混合压力,pa;τ为切应力张量;B为体受力向量,
能量方程:
组分质量:
式中:Sct为湍流施密特数;μt为湍流黏度,kg/(m·s);Di,m为质量扩散系数;
(2)膜凝结模型:
模型采用BELL蒸汽壁面冷凝模型,水蒸气扩散通过传质边界层,遇冷壁面释放潜热形成液膜,液体膜达到壁温,蒸汽扩散和凝结在气液截面持续发生,冷凝率取决于蒸汽的扩散速率,反过来又取决于边界层的水蒸汽浓度差,在边界层的一端是大量混合流体,浓度取决于混合物的参数,边界层的另一端是凝结液体水膜和气体混合物界面,假定在液体膜—水蒸汽混合物界面存在局部平衡,这意味着在界面处的水蒸汽浓度等于在膜温度下的饱和值,因此,膜界面温度决定了整个边界层的浓度差,并反过来控制冷凝率,假设如下所述:
(21)由于水蒸汽通过组分边界层扩散导致冷凝,这个假设用于导出冷凝水体积数学表达式,
(22)液膜和壁面存在热平衡,液膜温度Tfilm等于壁面温度Twall,
(23)气液界面是饱和状态。即交界面水蒸气分压力等于壁面温度下的饱和蒸汽压力,结合(22),壁面温度决定气液界面的蒸汽浓度,控制总的蒸汽扩散和通过边界层的浓度比,
(24)不考虑液膜运动带来的影响,由于流体剪切力和重力,液膜可能运动,但是运动不影响空气流动和凝结,
(25)分析仅限于膜状凝结,不考虑滴状冷凝,
(26)液膜的热阻忽略不计。
由于模型是基于以上假设建立的,凝结速率取决于水蒸气朝壁面的扩散速率,空气和水蒸气各组分的质量流量方程为:
用户自定义函数中,假定壁面温度等于或小于表面水蒸气分压力对应的饱和温度时发生凝结,如果温度高于饱和温度,水蒸气质量分数等于壁面邻近单元值,为了满足气液交界面局部热力学平衡的假设,当温度低于或等于饱和温度时,代表气液交界面的水蒸气质量分数被分配一个值,使得水蒸汽的分压等于水在局部壁面温度下的饱和压力,
混合物的体积质量源项为:
式中:Acell wall为壁面单元面积;Vcell为计算单体体积,
水蒸气的体积质量源项为:
(3)欧拉壁面模型:
欧拉壁面膜模型,实体表面膜不需划分网格,以表面虚拟膜代替,不影响核心流动动量场,多数情况下,表面薄膜的分离、脱落、蒸发不影响核心流动场的流动,通常情况下,模拟此种类型薄膜需要极大的计算资源,尤其在多相工况下精确计算相间通量,需要极细的网格来模拟水膜,而欧拉壁面膜模型不需要划分水膜网格,同时可以预测收集效率、耦合离散相模型,满足模拟壁面收集液滴形成水膜过程的需求,将水分考虑为小液滴,液滴在表面汇集成水膜,计算不考虑相变,
质量守恒方程:
膜动量守恒方程:
膜能量守恒方程:
式中:Ts为界面温度;Tf为膜平均温度;Vf为膜速度;kf为膜热导;Tw为壁面温度;假设膜的温度变化为线性分段变化,下半部分从壁面温度到膜平均温度Tf,上半部分由平均温度Tf到界面温度Ts;为液相撞击壁面的源项;为蒸发或凝结质量流率,
当离散粒子或液滴撞击壁面,会被液膜吸收,吸收后,其质量和动量增加到质量方程和动量方程的源项中,质量源项和动量源项分别为:
质量源项:
(4)水汽会在风力机叶片的前缘和尾缘位置凝结,故应在前缘和尾缘位置涂抹防冰涂料。
进一步的,使用的模型包括但不限于水蒸气凝结模型、欧拉壁面膜模型和颗粒磨损模型。
本发明的有益效果是:
通过本发明所公开的方法,确定不同条件下叶片表面上的结冰位置,在相应的结冰位置涂抹防冰涂料,既达到减轻结冰程度的目的,又节省涂料,降低成本。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本实施例的原理是通过CFD方法,采用数值模拟技术,模拟全尺寸叶片表面液滴撞击叶片表面的位置及水分汇集的位置,在该位置涂抹防冰涂料。
实施例1,本实施例公开了一种风力机叶片防冰涂料涂抹位置的确定方法,包括如下步骤:
(1)组分输运模型:
湿空气为干空气和水蒸气的混合气体,组分输运模型将混合气体定义为一种二元混合气体,具有两者的物理属性,求解过程作为单相处理,湿空气密度按照不可压缩理想气体方程计算,湿空气的控制方程如下所述,质量方程:
式中:ρm为混合物密度,kg/m3;U为流速;
动量方程:
式中:P为混合压力,pa;τ为切应力张量;B为体受力向量,
能量方程:
组分质量:
式中:Sct为湍流施密特数;μt为湍流黏度,kg/(m·s);Di,m为质量扩散系数;
(2)膜凝结模型:
模型采用BELL蒸汽壁面冷凝模型,水蒸气扩散通过传质边界层,遇冷壁面释放潜热形成液膜,液体膜达到壁温,蒸汽扩散和凝结在气液截面持续发生,冷凝率取决于蒸汽的扩散速率,反过来又取决于边界层的水蒸汽浓度差,在边界层的一端是大量混合流体,浓度取决于混合物的参数,边界层的另一端是凝结液体水膜和气体混合物界面,假定在液体膜—水蒸汽混合物界面存在局部平衡,这意味着在界面处的水蒸汽浓度等于在膜温度下的饱和值,因此,膜界面温度决定了整个边界层的浓度差,并反过来控制冷凝率,假设如下所述:
(21)由于水蒸汽通过组分边界层扩散导致冷凝,这个假设用于导出冷凝水体积数学表达式,
(22)液膜和壁面存在热平衡,液膜温度Tfilm等于壁面温度Twall,
(23)气液界面是饱和状态。即交界面水蒸气分压力等于壁面温度下的饱和蒸汽压力,结合(22),壁面温度决定气液界面的蒸汽浓度,控制总的蒸汽扩散和通过边界层的浓度比,
(24)不考虑液膜运动带来的影响,由于流体剪切力和重力,液膜可能运动,但是运动不影响空气流动和凝结,
(25)分析仅限于膜状凝结,不考虑滴状冷凝,
(26)液膜的热阻忽略不计。
由于模型是基于以上假设建立的,凝结速率取决于水蒸气朝壁面的扩散速率,空气和水蒸气各组分的质量流量方程为:
用户自定义函数中,假定壁面温度等于或小于表面水蒸气分压力对应的饱和温度时发生凝结,如果温度高于饱和温度,水蒸气质量分数等于壁面邻近单元值,为了满足气液交界面局部热力学平衡的假设,当温度低于或等于饱和温度时,代表气液交界面的水蒸气质量分数被分配一个值,使得水蒸汽的分压等于水在局部壁面温度下的饱和压力,
混合物的体积质量源项为:
式中:Acell wall为壁面单元面积;Vcell为计算单体体积,
水蒸气的体积质量源项为:
(3)欧拉壁面模型:
欧拉壁面膜模型,实体表面膜不需划分网格,以表面虚拟膜代替,不影响核心流动动量场,多数情况下,表面薄膜的分离、脱落、蒸发不影响核心流动场的流动,通常情况下,模拟此种类型薄膜需要极大的计算资源,尤其在多相工况下精确计算相间通量,需要极细的网格来模拟水膜,而欧拉壁面膜模型不需要划分水膜网格,同时可以预测收集效率、耦合离散相模型,满足模拟壁面收集液滴形成水膜过程的需求,将水分考虑为小液滴,液滴在表面汇集成水膜,计算不考虑相变,
质量守恒方程:
膜动量守恒方程:
膜能量守恒方程:
式中:Ts为界面温度;Tf为膜平均温度;Vf为膜速度;kf为膜热导;Tw为壁面温度;假设膜的温度变化为线性分段变化,下半部分从壁面温度到膜平均温度Tf,上半部分由平均温度Tf到界面温度Ts;为液相撞击壁面的源项;为蒸发或凝结质量流率,
当离散粒子或液滴撞击壁面,会被液膜吸收,吸收后,其质量和动量增加到质量方程和动量方程的源项中,质量源项和动量源项分别为:
质量源项:
(4)水汽会在风力机叶片的前缘和尾缘位置凝结,故应在前缘和尾缘位置涂抹防冰涂料。
本实施例使用的模型包括但不限于水蒸气凝结模型、欧拉壁面膜模型和颗粒磨损模型。
Claims (2)
1.一种风力机叶片防冰涂料涂抹位置的确定方法,其特征在于,包括如下步骤:
(1)组分输运模型:
湿空气为干空气和水蒸气的混合气体,组分输运模型将混合气体定义为一种二元混合气体,具有两者的物理属性,求解过程作为单相处理,湿空气密度按照不可压缩理想气体方程计算,湿空气的控制方程如下所述,质量方程:
动量方程:
能量方程:
式中:为净源项;q为净热通量;h为焓值;
组分质量:
(2)膜凝结模型:
模型采用BELL蒸汽壁面冷凝模型,水蒸气扩散通过传质边界层,遇冷壁面释放潜热形成液膜,液体膜达到壁温,蒸汽扩散和凝结在气液截面持续发生,冷凝率取决于蒸汽的扩散速率,反过来又取决于边界层的水蒸汽浓度差,在边界层的一端是大量混合流体,浓度取决于混合物的参数,边界层的另一端是凝结液体水膜和气体混合物界面,假定在液体膜—水蒸汽混合物界面存在局部平衡,这意味着在界面处的水蒸汽浓度等于在膜温度下的饱和值,因此,膜界面温度决定了整个边界层的浓度差,并反过来控制冷凝率,假设如下所述:
(21)由于水蒸汽通过组分边界层扩散导致冷凝,这个假设用于导出冷凝水体积数学表达式,
(22)液膜和壁面存在热平衡,液膜温度Tfilm等于壁面温度Twall,
(23)气液界面是饱和状态,交界面水蒸气分压力等于壁面温度下的饱和蒸汽压力,结合(22),壁面温度决定气液界面的蒸汽浓度,控制总的蒸汽扩散和通过边界层的浓度比,
(24)不考虑液膜运动带来的影响,由于流体剪切力和重力,液膜可能运动,但是运动不影响空气流动和凝结,
(25)分析仅限于膜状凝结,不考虑滴状冷凝,
(26)液膜的热阻忽略不计:
由于模型是基于以上假设建立的,凝结速率取决于水蒸气朝壁面的扩散速率,空气和水蒸气各组分的质量流量方程为:
用户自定义函数中,假定壁面温度等于或小于表面水蒸气分压力对应的饱和温度时发生凝结,如果温度高于饱和温度,水蒸气质量分数等于壁面邻近单元值,为了满足气液交界面局部热力学平衡的假设,当温度低于或等于饱和温度时,代表气液交界面的水蒸气质量分数被分配一个值,使得水蒸汽的分压等于水在局部壁面温度下的饱和压力,
混合物的体积质量源项为:
水蒸气的体积质量源项为:
(3)欧拉壁面模型:
欧拉壁面膜模型,实体表面膜不需划分网格,以表面虚拟膜代替,不影响核心流动动量场,多数情况下,表面薄膜的分离、脱落、蒸发不影响核心流动场的流动,通常情况下,模拟此种类型薄膜需要极大的计算资源,尤其在多相工况下精确计算相间通量,需要极细的网格来模拟水膜,而欧拉壁面膜模型不需要划分水膜网格,同时可以预测收集效率、耦合离散相模型,满足模拟壁面收集液滴形成水膜过程的需求,将水分考虑为小液滴,液滴在表面汇集成水膜,计算不考虑相变,
质量守恒方程:
膜动量守恒方程:
膜能量守恒方程:
式中:Ts为界面温度;Tf为膜平均温度;Vf为膜速度;kf为膜热导;Tw为壁面温度;假设膜的温度变化为线性分段变化,下半部分从壁面温度到膜平均温度Tf,上半部分由平均温度Tf到界面温度Ts;为液相撞击壁面的源项;为蒸发或凝结质量流率,
当离散粒子或液滴撞击壁面,会被液膜吸收,吸收后,其质量和动量增加到质量方程和动量方程的源项中,质量源项和动量源项分别为:
质量源项:
(4)水汽会在风力机叶片的前缘和尾缘位置凝结,故应在前缘和尾缘位置涂抹防冰涂料。
2.根据权利要求1所述风力机叶片防冰涂料涂抹位置的确定方法,其特征在于:使用的模型包括水蒸气凝结模型、欧拉壁面膜模型和颗粒磨损模型。
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 CN111241761A (zh) | 2020-06-05 |
CN111241761B true 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 (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 西安交通大学 | 一种压水反应堆安全壳两侧冷凝与蒸发耦合计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT506401B1 (de) * | 2008-06-24 | 2009-09-15 | Fronius Int Gmbh | Vorrichtung und verfahren zum reinigen von schweissbrennern |
-
2020
- 2020-02-05 CN CN202010080304.9A patent/CN111241761B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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)
Title |
---|
"基于欧拉壁面液膜模型的三维热气";李延 等;《北京航空航天大学学报》;20180531;正文第1.3节 * |
Also Published As
Publication number | Publication date |
---|---|
CN111241761A (zh) | 2020-06-05 |
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 | |
Birbarah et al. | External convective jumping-droplet condensation on a flat plate | |
Cui et al. | Study on the dynamic and thermal performances of a reversibly used cooling tower with upward spraying | |
Fortin et al. | New roughness computation method and geometric accretion model for airfoil icing | |
CN101561343B (zh) | 自然通风盐水模型实验装置 | |
Huang et al. | Laminar mixed convection heat and mass transfer in vertical rectangular ducts with film evaporation and condensation | |
CN114896906B (zh) | 一种考虑冰层和固体壁面中导热的积冰模拟方法 | |
CN111967201A (zh) | 一种基于数值仿真模型分析临界覆冰类型的方法 | |
Birbarah et al. | Internal convective jumping-droplet condensation in tubes | |
CN112985347A (zh) | 一种飞机结冰表面粗糙度计算方法 | |
Gu et al. | Effects of mass transfer on heat and mass transfer characteristics between water surface and airstream | |
CN111241761B (zh) | 一种风力机叶片防冰涂料涂抹位置的确定方法 | |
Shotlou et al. | Innovative method to reduce frost formation in liquid-to-air membrane energy exchangers (LAMEE) based on 3D CFD simulation | |
Qiu et al. | Development and validation of a numerical model for frost growth based on nucleation theory | |
Jia et al. | Research on the heat and mass transfer characteristics of fin-tube exchanger under low pressure environment | |
Ma et al. | Numerical simulation of frosting on wavy fin-and-tube heat exchanger surfaces | |
Enayatollahi et al. | Characterising the heat and mass transfer coefficients for a crossflow interaction of air and water | |
Sivasamy et al. | Numerical simulation of two‐dimensional laminar slot‐jet impingement flows confined by a parallel wall | |
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 | |
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的城市水体蒸发模拟方法 | |
Da Silveira et al. | Numerical simulation of ice accretion on the leading edge of aerodynamic profiles |
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 |