风能设备仿真软件：Fluent二次开发_（24）.风力发电机组结冰与除冰仿真.docx

PAGE1

PAGE1

风力发电机组结冰与除冰仿真

1.结冰仿真原理

1.1结冰的物理过程

风力发电机组在寒冷气候条件下工作时，叶片表面可能会结冰。结冰对风力发电机组的性能有显著影响，包括降低风能利用率、增加叶片重量和不平衡问题、增加机械应力等。因此，对结冰过程的仿真至关重要。

结冰过程主要包括以下几个物理阶段：

水滴撞击叶片表面：当含有水滴的空气流过叶片时，水滴会撞击叶片表面。

水滴沉积：部分水滴会在叶片表面沉积，形成液态水层。

水层冻结：液态水层在低温下逐渐冻结，形成冰层。

冰层增长：随着更多水滴的撞击和沉积，冰层会逐渐增长。

1.2结冰模型

在Fluent中，结冰过程可以通过多种模型进行仿真，常见的模型包括：

Messinger模型：用于模拟水滴撞击和沉积过程。

Leeds模型：用于模拟水层冻结过程。

ICEPAK模型：用于模拟冰层增长过程。

1.2.1Messinger模型

Messinger模型主要用于计算水滴在叶片表面的撞击和沉积。该模型基于水滴的动能、动量和热能交换，计算水滴撞击叶片后的沉积效率。模型的数学表达式如下：

η

其中：

η为水滴的沉积效率。

L为水滴的撞击动能。

Lc

ΔL

1.2.2Leeds模型

Leeds模型用于模拟水层的冻结过程。该模型基于水层的温度、厚度和热传导特性，计算水层的冻结速率。模型的数学表达式如下：

?

其中：

h为冰层厚度。

t为时间。

qf

ρi

Lf

1.2.3ICEPAK模型

ICEPAK模型用于模拟冰层的增长过程。该模型考虑了水滴的持续撞击、水层的冻结和冰层的重力滑移。模型的数学表达式如下：

?

其中：

hw

hs

1.3结冰仿真设置

在Fluent中设置结冰仿真需要以下几个步骤：

定义流体和固体材料：包括空气、水滴和叶片材料。

设置初始和边界条件：包括温度、湿度、风速等。

选择和配置结冰模型：选择适合的结冰模型并配置相关参数。

网格划分：对计算域进行适当的网格划分。

求解设置：选择求解器和求解方法，设置迭代参数。

1.3.1定义流体和固体材料

在Fluent中，首先需要定义流体和固体材料的属性。例如，定义空气、水滴和叶片材料的密度、比热、导热系数等。

#定义空气属性

air_properties={

density:1.225,#kg/m^3

specific_heat:1005,#J/kg·K

thermal_conductivity:0.0257#W/m·K

}

#定义水滴属性

water_drop_properties={

density:997,#kg/m^3

specific_heat:4181,#J/kg·K

thermal_conductivity:0.598#W/m·K

}

#定义叶片材料属性

blade_material_properties={

density:2500,#kg/m^3

specific_heat:880,#J/kg·K

thermal_conductivity:0.7#W/m·K

}

1.3.2设置初始和边界条件

初始和边界条件的设置对于结冰仿真至关重要。例如，设置温度、湿度、风速等。

#设置初始条件

initial_conditions={

temperature:273.15,#K

humidity:0.5,#相对湿度

wind_speed:10#m/s

}

#设置边界条件

boundary_conditions={

inlet:{

temperature:268.15,#K

humidity:0.8,#相对湿度

wind_speed:10#m/s

},

outlet:{

pressure:101325#Pa

},

blade_surface:{

temperature:273.15,#K

roughness:0.001#m

}

}

1.3.3选择和配置结冰模型

选择合适的结冰模型并配置相关参数。例如，选择Messinger模型并设置临界动能和动能扩散范围。

#选择Messinger模型

model=Messinger

#设置模型参数

model_parameters=