风力机机舱内自然对流数值模拟 戚中浩.docVIP

中国工程热物理学会 学科类别： 热机气动热力学 学术会议论文 编号：092103 风力机机舱内自然对流数值模拟 戚中浩 张晓东 康顺 华北电力大学电站设备状态监测与控制教育部重点实验室 北京102206 联系电话：010 E-mail：kangs@ncepu.edu.cn 摘要:本文对不同环境温度下，风力机机舱内自然对流换热现象进行了三维数值模拟。研究了在确定的热源壁面温度条件下，机舱内对流换热强度和流动形态随环境温度的变化关系，得到了换热面平均尔数与表征换热温差的瑞利数Ra的关系曲线。分析了热源体位置变化对换热和机舱内流态的影响，进行了热源体布置优化，本文的研究基于商业软件，Spalart-Allmaras湍流模型。 关键词： 自然对流 二维封闭方腔 风力机机舱 0引言 一次能源的枯竭，让风力发电越来越受到世界各国的重视并得到了广泛的开发和利用。然而大气湍流、阵风、雷电、结冰等环境气候条件都会影响到风力机安全运行和正常出力。在以往的风力机事故中，机舱，及机舱内部零件发生故障占了很大的比例。那么，如何保证风力机机舱、机舱内部零部件等在不同的环境条件，尤其是在低温、结冰等恶劣环境下安全正常运行也就十分必要。如文献[1]对风力机机舱在高温和低温两种情况下进行了数值分析，研究了机舱外的风速、风温对机舱内部温度的影响。文献[2]对空气流过风力机机舱内发电机的热分析。文献[3]对转子和机舱间相互影响做了气动分析。随着风电装机容量的快速增长，将有越来越多的风电场在气候条件复杂的地区建设， 本文结合东方电气集团东方汽轮机有限公司2.5兆瓦风力机研发项目，对风力机机舱内部环境进行了CFD数值模拟研究，分析了环境温度和热原体位置对机舱内部自然对流换热和流动形态的影响。所建立的数值计算模型通过二维封闭方腔的自然对流试验数据进行了模型验证。 机舱的几何结构见图1，模型简化的几何尺寸在下面分析中详细给出。 本文假定空气视为不可压缩流体浮升力计算 1.2数值方法文章采用商业软件NUMECA，Fine/hexa软件包中的HEXA_STREAM求解器，求解时间相关的雷诺时均方程，使用二阶中心格式离散，时间推进采用四阶Runge-Kutta方法。湍流的封闭模型为S-A模型。采用多重网格加速。 迭代计算收敛、完成的条件为。其中表示U、V、W、P和T，i、j、k表示某一个节点坐标，m表示当前迭代步数。 2算例分析 2.1二维方腔内自然对流换热 二维封闭方腔如图2，垂直高为2.5m，水平宽为0.5m。方腔竖直壁面为恒温、水平壁面绝热，左侧，温度为侧，温度为 数值计算使用非结构网格，网格布置进入层流底层。 2.1.2 y=1.25m处，速度V、温度T水平方向的分布 图3 a是方腔内自然对流换热流动状态分布极 图3 b为温度分布。处与实验值吻合较好， 2.1.3 温度等值线图、速度矢量图、速度云图 图4是二维方腔内温度场和速度场分布。图4 a是温度等值线图，从图上知道，温度的极大值出现在左上角，极小值出现在右下角。在Thhh、Tc壁面附近形成热边界层，温度变化剧烈，温度梯度大，而在方腔的中心区，温度层状分布。这与图3 b温度曲线变化是一致的。图4 b、c是方腔内部被加热的空气被冷却的空气图5是风力机机舱的简化模型，机舱沿x、y、z轴是L(m)、W(m)、H(m)，在机舱底面上有一个凸出的长方体热源，长宽高是L ’ (m)、W ’ (m)、H ’ (m)，热源在机舱底面上的位置用热源底面中心表示，中心点坐标是（M，N），其中M= m+ L ’/，N= n+ W ’/2。所有机舱壁面设为冷壁面，温度为，热源所有壁面设为热壁面，温度为5个 数值壁面附近区域热源平均努尔数 2.2.1机舱壁温度对机舱内自然对流换热影响 热源位置、壁面温度（40℃）保持不变，机舱壁面温度从30℃变到-20℃，热源面i的平均努尔数瑞利数如图6所示通过方程（1）定义 （i=a、b、c、d、e） （1） Ai第i个热源面面积，是温度无量纲 瑞利数Ra为 （2） 其中温差=Th-Tc，是空气热扩散系数，ν是空气运动粘滞系数，是空气体膨胀系数，H是特征长度，取机舱内热源的高度。 从图6上可以知道，每个热源面随Ra增大均逐渐增大。热源面a的a大于其它4个面的平均努塞尔数，热源面e（顶面）的e比其它4个热源面小。是由于热源面a比其它4个热源面接触冷空气多，热交换量要大，