1、基于基于Airpak的绿色机房模型仿真评价与控制的绿色机房模型仿真评价与控制吉林大学吉林大学概述概述机房模型建立及求解变量说明3问题重述与分析3 12模型假设模型假设模型评价与推广3 54一一.问题重述与分析问题重述与分析?随着人们对高密度计算,多任务计算的需求,越来越多的高性能数据中心或互联网中心（DC,IDC）建成，为了保证机房内设备健康运行,数据中心制冷系统必须根据机房内热点的温度（室内最高温度）向机房送配冷气.合理地给服务器分配工作任务,能够降低机房内热点的温度,达到节能目的.绘制机房冷,热通道的热分布及流场分布,求出室内热点的位置;机房热分布数学模型的建立与比较;机房实际任务量下的最
2、优任务分配方案求解;按照电子信息系统机房设计规范 C级要求控制机房温度,服务器设计任务量一定条件下,对空调的送风速度或送风温度的合理控制.?机房内流体为不可压缩牛顿流体,且符合Boussinesq假设;?室内流体为粘性流体,流体作定常流动,忽略由流体粘性力做功所引起的耗散热;?假定流场具有高的紊流Re数,流体的紊流粘性各向同性;?不考虑机房外部传热的影响;?将机房近似看作封闭系统;?不考虑机房内其他辅助设备(如PDU等)对整个模型影响;?流动为紊态流动;?忽略固体壁面和室内物体表面的热辐射.二.模型假设时间维数扩散通量W(t)空调系统的工作负荷三.变量说明变量变量说明各通道的横坐标值,纵坐标值
3、,高度值通用变量速度矢量源项动量热源,污染源等引起的源项为湍流脉动动能的耗散率各通道主机的横坐标值,纵坐标值,高度值湍流速度编号为的机柜任务分配因子编号为的机柜附近处热点温度编号为的机柜的发热功率变量变量说明空气动力粘度,取机房中机柜编号,R机房内热点处的温度机柜满负荷工作时发热功率标定温度法 中的标定温度距离变量湍流脉动功能空气密度,取1.237三.变量说明四四.机房模型建立及求解机房模型建立及求解问题一:绘制机房冷,热通道的热分布及流场分布,求出室内热点的位置1.冷,热通道的热分布热通道冷通道问题一:绘制机房冷,热通道的热分布及流场分布,求出室内热点的位置2.冷,热通道的速度分布热通道冷通
4、道四四.机房模型建立及求解机房模型建立及求解问题一:绘制机房冷,热通道的热分布及流场分布,求出室内热点的位置3.机房内最高温度位置四四.机房模型建立及求解机房模型建立及求解利用MATLAB编程可求解室内温度最高位置在通道3（热通道）,处,最高温度为56.1725问题二:机房热分布数学模型的建立与比较1.CFD分析方法四四.机房模型建立及求解机房模型建立及求解建立模型（Build Model）生成网格（Generate Mesh）检测结果（Examine Results）计算求解（Calculate Solution）AirpakFLUENT问题二:机房热分布数学模型的建立与比较2.研究对象的物
5、理模型边界条件边界条件整个计算空间大约 283m3,考虑到机房的复杂性,采用混合网格对计算区域进行离散,即采用变步长的方法对网格进行划分,对气流速度变化较快的区域进行不同程度的网格加密,根据题目中所给数据插值来施加周期边界条件:a)送风口和回风口分别设为风速进口边界和出口边界;b)机架设为热通量边界,数值上取其总发热量与散热表面积之比;c)机房的边界设为恒温边界.四.机房模型建立及求解问题二:机房热分布数学模型的建立与比较3.网格划分由于机房内结构复杂,为了更准确的反映房内的气流组织,在划分网格时对有物体的地方,如空调,机柜群等处网格相对其他地方更细密.右图为不同截面的网格分布.四.机房模型建
6、立及求解问题二:机房热分布数学模型的建立与比较4.计算求解建立数学模型是对所研究的流动问题进行数学描述,连续性方程,动量方程,能量方程都可以写成下列通用控制微分方程的形式:连续性方程:K方程:方程:由于控制方程是非线性的,所以在达到收敛精度之前,需进行迭代求解.我们采用的Airpak利用有限体积法,把计算区域划分为离散的控制体积网格,在每个控制体积上积分控制方程,形成计算变量的代数方程,然后进行迭代求解.四.机房模型建立及求解问题二:机房热分布数学模型的建立与比较5.检查结果I.热分布图按附件数据给定相关参数后得到机房的热分布图.冷通道（通道2）热分布热通道（通道3）热分布1）下图为冷,热通道
7、的热分布图.四.机房模型建立及求解问题二:机房热分布数学模型的建立与比较5.检查结果I.热分布图2）下图分别为距离空调位置 2.4m,5m和7.2m的截面热分布图.距空调2.4m截面:距空调位置5m:距空调位置7.2m截面:四.机房模型建立及求解问题二:机房热分布数学模型的建立与比较5.检查结果I.热分布图3）下图分别为模型的整体热分布立体图和气流图.热分布立体图热分布气流图四.机房模型建立及求解问题二:机房热分布数学模型的建立与比较5.检查结果II.速度分布图按附件数据给定相关参数后得到机房的速度分布图.冷通道（通道2）速度分布热通道（通道3）速度分布1）下图为冷,热通道的速度分布图.四.机
8、房模型建立及求解问题二:机房热分布数学模型的建立与比较5.检查结果II.速度分布图2）图19,图20,图21分别为距离空调位置2.4m,5m和7.2m的截面速度分布图.距空调2.4m截面:距空调位置5m:距空调位置7.2m截面:四.机房模型建立及求解问题二:机房热分布数学模型的建立与比较5.检查结果II.速度分布图3）图17,图18分别为模型的整体速度分布立体图和气流图.速度分布立体图速度分布气流图四.机房模型建立及求解问题二:机房热分布数学模型的建立与比较5.检查结果III.结果分析与比较根据附件1绘制的热通道热分布云图可知在高度 Z一定时,随着距离空调位置Y的增加,温度呈现先增大后减小,然
9、后又增大又减小的双峰分布;在距离空调位置Y一定时,随着高度Z的增加,温度逐渐增大,然后逐渐下降的趋势,但减小趋势不明显.根据Airpak绘制的热通道（也取通道3）热分布可知,在高度Z一定时,随着距离空调位置Y的增加,温度先增大后减小,然后又增大又减小的双峰分布.在距离空调位置Y一定时,随着高度Z的增加,温度总体呈现先增大后减小的趋势,减小趋势也不明显.同样比较可知,对于冷通道附件1和Airpak绘制的图形也吻合.通过Airpak仿真得出的机房热点温度为56.1833,出现在 Y=6.4617m,Z=1.8261m 的位置.利用MATLAB编程可求解室内温度最高位置在通道3（热通道）Y=6.83
10、82m,Z=2.145m 处,最高温度为56.1725.四四.机房模型建立及求解机房模型建立及求解问题三:机房实际任务量下的最优任务分配方案求解1.基本思路通过查阅文献可以得知,合理地给服务器分配工作任务,能够降低机房内热点的温度,达到节能目的.其基本思路是对于热点（温度高）区域附近的服务器需减少工作量,以降低产热;热点区域服务器减少的工作量由冷点（温度低）区域附近服务器承担,使冷点区域的温度上升,最终使整个机房内温度不会出现明显差异（热点与冷点的温度差减小）,达到降低机房内热点温度值的目的.采用方法1.普通方法2.智能方法.倒数法.倒序法.标定温度法平均分配四.机房模型建立及求解问题三:机房
11、实际任务量下的最优任务分配方案求解2.服务器实际任务量为0.8的普通平均分配方案为使服务器的实际任务量为0.8,可将每个机柜的功率改为原来的80%（由第二问Airpak仿真建模可知任务量为1时的发热功率为6700W）,即每个机柜的发热功率为5360W.采用此方案时,通过Airpak仿真建模可以求出机房的整体热分布立体图和热分布气流图如下.普通平均分配方案热分布立体图普通平均分配方案热分布气流图此方案下的热点温度为此方案下的热点温度为47.5795,47.5795,机房内平均温度为机房内平均温度为26.7.26.7.四四.机房模型建立及求解机房模型建立及求解问题三:机房实际任务量下的最优任务分配
12、方案求解3.服务器实际任务量为0.8的智能分配方案I.方案1 倒数法 智能分配方案此方案采用“倒数法”对各机柜的功率重新分配任务量,即编号为的机柜任务分配因子的关系式为:重新分配后编号为的机柜的发热功率),(nmnm,?8,41,11nmnmnmnmnmRR?nmnmPP,?),(nm四.机房模型建立及求解问题三:机房实际任务量下的最优任务分配方案求解3.服务器实际任务量为0.8的智能分配方案I.方案1 倒数法 智能分配方案具体做法如下.2.将（1）中计算的数据取倒数,得到;3.计算各机柜任务分配因子,即用每个机柜的温度倒数除以温度倒数总和;4.按上面的各机柜任务分配因子重新分配各机柜的发热功
13、率;5.将发热功率输入到用Airpak建立的仿真模型中,可以求得采用此方案时,机房的整体热分布立体图和热分布气流图.1.利用Airpak求出各机柜热道附近的热点温度;此方案的热点温度为44.8839,机房平均温度为27.2.四四.机房模型建立及求解机房模型建立及求解问题三:机房实际任务量下的最优任务分配方案求解3.服务器实际任务量为0.8的智能分配方案II.方案2 倒序法 智能分配方案此方案采用“倒序法”对各机柜的功率重新分配任务量,即将各机柜附近处热点温度进行排序,将热点温度最大的机柜赋予最低的热点温度,而热点温度最低的机柜赋予最高的温度,并将赋予温度与比作为各机柜任务分配因子,即则重新分配
14、后编号为的机柜的发热功率的计算公式为nmR,nmR?8,41,nmnmnmRnm,?8,41,nmnmnmnmnmRR?),(nmnmP,nmnmPP,?四.机房模型建立及求解问题三:机房实际任务量下的最优任务分配方案求解3.服务器实际任务量为0.8的智能分配方案II.方案2 倒序法 智能分配方案参照方案1的计算步骤我们可以得出各机柜任务分配因子和各机柜发热功率将发热功率输入到用Airpak建立的仿真模型中,可以求得采用此方案时,机房的整体热分布立体图和热分布气流图(如右图所示).四四.机房模型建立及求解机房模型建立及求解此方案下的热点温度为46.2828,机房内平均温度为27.1.问题三:机
15、房实际任务量下的最优任务分配方案求解3.服务器实际任务量为0.8的智能分配方案III.方案3 标定温度法 智能分配方案此方案采用 标定温度法 对各机柜的功率重新分配任务量,即考虑标定温度 T与各机柜附近处热点温度差值,具体来讲各机柜任务分配因子公式为则重新分配后编号为的机柜的发热功率的计算公式为对于标定温度,我们给定四.机房模型建立及求解问题三:机房实际任务量下的最优任务分配方案求解3.服务器实际任务量为0.8的智能分配方案III.方案3 标定温度法 智能分配方案四四.机房模型建立及求解机房模型建立及求解参照方案1的计算步骤我们可以得出各机柜任务分配因子和各机柜发热功率将发热功率输入到用Air
16、pak建立的仿真模型中,可以求得采用此方案时,机房的整体热分布立体图和热分布气流图(如右图所示).此方案下的热点温度为41.0293,机房内平均温度为27.7.问题三:机房实际任务量下的最优任务分配方案求解4.服务器实际任务量为0.8的各智能分配方案评价与优化I.服务器实际任务量为0.8的各智能分配方案评价通过使用Airpak仿真分析能够得到上述三种智能分配方案均能降低机房内热点温度,其中 倒序法 智能分配方案对降低机房内热点温度效果最不明显,而 标定温度法 智能分配方案能明显降低机房内热点温度尬倒数法 智能分配方案介于两者之间.四四.机房模型建立及求解机房模型建立及求解倒数法 智能方案倒序法
17、 智能方案标定温度法 智能方案热点温度:44.8839热点温度:46.2828热点温度:41.0293问题三:机房实际任务量下的最优任务分配方案求解4.服务器实际任务量为0.8的各智能分配方案评价与优化?标定温度法 智能分配方案优化我们通过改变 标定温度法 智能分配方案中的标定温度来找到最优的标定温度T,各方案中T 值见下表.方案1方案2方案3方案4方案16T（）5051525365通过不同方案中的标定温度 T 求出不同方案各机柜的发热功率,并将其输入到用Airpak建立的仿真模型中,通过比较各方案的热点值的不同,我们得出进而得到服务器实际任务量为 0.8时最优任务分配方案是方案最优任务分配方
18、案是方案7,7,即标定温度即标定温度T=56T=56时.四四.机房模型建立及求解机房模型建立及求解问题三:机房实际任务量下的最优任务分配方案求解4.服务器实际任务量为0.8的各智能分配方案评价与优化?标定温度法 智能分配方案优化以下是该最优任务分配方案（T=56的 标定温度法 智能分配方案）机房的整体热分布立体图和热分布气流图.从热分布立体图和热分布气流图中可以明显观查到此时各机柜温度分布相当,温度均匀.四四.机房模型建立及求解机房模型建立及求解此方案下的热点温度为此方案下的热点温度为40.4291,40.4291,机房内平均温度为机房内平均温度为28,28,问题三:机房实际任务量下的最优任务
19、分配方案求解5.服务器实际任务量为0.5时最优任务分配方案采取与服务器实际任务量为 0.8时同样的思路,即首先分别采取倒数法 智能分配方案尬倒序法 智能分配方案以及 标定温度法 智能分配方案可以得出不同方案各机柜的发热功率,并将其输入到用Airpak建立的仿真模型中,通过比较各方案的热点值的不同,我们同样可以得到 标定温度法 智能分配方案相对于其他方法而言是最好的,然后,通过改变 标定温度法 智能分配方案中标定温度T,得到不同的优化方案（如下表）.方案1方案2方案3方案4方案16T（）3536373844四四.机房模型建立及求解机房模型建立及求解问题三:机房实际任务量下的最优任务分配方案求解5
20、.服务器实际任务量为0.5时最优任务分配方案将不同方案中的标定温度T 求出不同方案各机柜的发热功率,并将其输入到仿真模型中,通过比较各方案的热点值的不同,我们得出进而得到服务器实际任务量为0.5时最优分配方案是方案3,即标定温度T=37时.T=37的 标定温度法 智能分配方案中机房的整体热分布立体图和热分布气流图.从热分布立体图和热分布气流图中可以明显观查到此时各机柜温度分布相当,温度均匀.此方案下的热点温度为此方案下的热点温度为29.3265,29.3265,机房内平均温度为机房内平均温度为22.122.1四四.机房模型建立及求解机房模型建立及求解问题四:按照规范中C级要求控制机房温度,在任
21、务量一定条件下,对空调的送风速度或送风温度进行合理控制从电子信息系统机房设计规范中我们可以了解到,各级电子信息系统机房技术要求中对于机房温度具有如下规定:项目技术要求备注A级B级C级主机房温度（开机时）2311828不得结算合理控制机房温度（开机时）机房制冷功耗出风口速度V送风槽温度t四.机房模型建立及求解问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制1.符合电子信息系统机房设计规范C级规定的送风槽出口速度V,送风槽出口温度t范围由题目中附件二提供的数据可知,空调的送风速度在0.6-1.2之间,首先通过不断的改变送风速度和改变送风温度,可以得
22、到不同的机房温度模型,并将其数值输入到Airpak中进行仿真分析.并从中读取出不同方案的平均温度,进而得到送风速度,送风温度和机房平均的温度值PT三者间的若干个散点的关系.V（m/s）t（）0.310.50.710.951.101.111025.223.323.122.421.5211227.225.425.724.523.5231732.230.730.729.528.628四四.机房模型建立及求解机房模型建立及求解问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制1.符合电子信息系统机房设计规范C级规定的送风槽出口速度V,送风槽出口温度t范围
23、将得到的散点通过MATLAB中的Curve Fitting Toolbox工具箱拟合成一个曲面并求出该曲面.根据电子信息系统机房设计规范C级要求,找到机房平均的温度为18-28时对应的曲面即满足要求的送风速度和送风温度之间的关系面.四.机房模型建立及求解通过MATLAB中的Curve Fitting Toolbox工具箱拟合成的曲面方程为:f(x,y)=p00+p10*x+p01*y+p20*x2+p11*x*y+p02*y2+p30*x3+p21*x2*y+p12*x*y2其中参数为:p00=22.43 p02=-0.01091p30=-26.26 p10=-45.36 p20=59.98
24、p21=-0.3609p01=1.185 p11=0.3865 p12=0.004225 拟合函数的检验拟合函数的检验:通过对曲面的拟合进行检验（如表13）可以看出拟合出的曲面是正确的,是可信的,用拟合后的曲面进行求解是可行的.问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制1.符合电子信息系统机房设计规范C级规定的送风槽出口速度V,送风槽出口温度t范围四四.机房模型建立及求解机房模型建立及求解判别因子平方和误差判定系数均方根误差数值0.76540.99620.29162.出口速度 V,送风槽出口温度 t 与机房内热点值R的关系同样地,由题目中
25、附件二提供的数据可知,空调的送风速度在0.6-1.2之间,首先通过不断的改变送风速度和改变送风温度,可以得到不同的机房温度模型,并将其数值输入到Airpak中进行仿真分析.并从中读取出不同方案的热点温度,进而得到送风速度,送风温度和机房热点的温度值R三者间的若干个散点关系.0.310.50.650.710.830.951.041.11.1051.11027.37 27.10 27.07 26.89 26.11 25.42 25.12 24.52 24.1224.02431.41 31.42 31.22 31.00 30.15 29.50 29.17 28.58 28.1428.00835.48
26、 35.70 35.36 35.09 34.18 33.57 33.24 32.61 32.1532.0012.847.57 40.78 40.34 40.00 39.01 38.45 38.10 37.48 36.9536.791749.76 45.14 44.70 44.30 43.23 42.72 42.35 41.71 41.1640.98VtVt问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制四四.机房模型建立及求解机房模型建立及求解同样地,由题目中附件二提供的数据可知,空调的送风速度在0.6-1.2之间,首先通过不断的改变送风速度和
27、改变送风温度,可以得到不同的机房温度模型,进而得到出口速度 V,送风槽出口温度 t 与机房内热点值R的关系.问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制2.出口速度 V,送风槽出口温度 t 与机房内热点值R的关系四.机房模型建立及求解拟合函数的检验:通过对曲面的拟合进行检验可以看出拟合出的曲面是正确的,用拟合出的曲面进行求解是可行的.f(x,y)=p00+p10*x+p01*y+p20*x2+p11*x*y+p02*y2+p30*x3+p21*x2*y+p12*x*y2+p03*y3+p40*x4+p31*x3*y+p22*x2*y2+p1
28、3*x*y3+p04*y4其中参数为:p00=26.17 p12=-0.1416p21=4.147p31=-0.9991p22=-0.0896p10=2.615 p03=-0.0005p13=0.01177p22=-0.0896p11=-3.16p01=1.403 p40=2.52p04=-0.0003p02=0.1106p30=-10.32p20=4.185 p31=-0.9991问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制2.出口速度 V,送风槽出口温度 t 与机房内热点值R的关系四四.机房模型建立及求解机房模型建立及求解判别因子平方和
29、误差判定系数均方根误差数值4.45750.98950.8010通过MATLAB中的Curve Fitting Toolbox工具箱拟合成的曲面方程为:3.出口速度 V,送风槽出口温度t 与空调功耗的关系由资料可以得出空调系统消耗电能产生冷/热空气输入建筑物室内,用以改变室内热环境.因此,致冷或加热空气的能耗和空气总量就是k空调系统的工作负荷 W(t),设制冷或加热空气的能耗函数为W1(t),它由热/冷量Qu(t)和空调设备的参数k决定;空调风机的能耗函数为W2(t),由风量Vk和风机参数k决定.W(t)的函数关系式如下:),(),()(21?tVkWtQWtWu?问题四:按照规范中C级要求控制
30、机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制四.机房模型建立及求解4.实际任务量为设计任务量的0.8时,送风速度V和送风温度t控制的最优方案由上述分析可以得出,出风槽出口速度 V,送风槽出口温度 t的取值与平均温度PT,房间内热点值R以及空调功耗W(t)密切相关,为了充分的分析评价,我们通过改变出风槽出口速度,送风槽出口温度 的取值建立多个控制方案.0.600.650.701.151.2010.0(10.0,0.60)(10.0,0.65)(10.0,0.70)(10.0,1.15)(10.0,1.20)10.5(10.5,0.60)(10.5,0.65)(10.5,0
31、.70)(10.5,1.15)(10.5,1.20)11.0(11.0,0.60)(11.0,0.65)(11.0,0.70)(11.0,1.15)(11.0,1.20)15.5(15.5,0.60)(15.5,0.65)(15.5,0.70)(15.5,1.15)(15.5,1.20)16.0(16.0,0.600(16.0,0.65)(16.0,0.70)(16.0,1.15)(16.0,1.20)t问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制四四.机房模型建立及求解机房模型建立及求解V4.实际任务量为设计任务量的0.8时,送风速度V和
32、送风温度t控制的最优方案综合考虑每个控制方案分别评价房间内热点温度R和功耗W(t),和电机的发热,求得最优的控制方案是送风速度V=0.70m/s,送风温度t=13,此时机房的整体热分布立体图和热分布气流图如下.问题四:按照规范中C级要求控制机房温度,在任务量一定条件下,对空调的送风速度或送风温度进行合理控制四四.机房模型建立及求解机房模型建立及求解此方案下的热点温度为29.3265,机房内平均温度为22.1.五.模型评价与推广模型优缺点分析1.优点模型假设合理 我们根据实际情况对模型进行了合理的简化,在模型的求解过程中既保证了精度,又提高了效率;模型设置合理 我们充分考虑机房热点和机柜温度的影
33、响因素,并对房间内的流体进行了合理的假设,选择了在此条件下求解精度较高的k-方程,设置合理的边界条件并运用Airpak软件建模,进行仿真求解.不断进行查找资料和仿真分析,充分考虑数据中心的节能等问题,进一步优化我们的模型,再与实际相结合得出我们的最优模型;求解结果合理 通过上述模型求解的数据与附录中所给数据差距很小,得出的结论与查阅的文献资料相符,因而我们得出的结果是合理的;模型假设,建立,求解思路清晰 可开发潜力大,对同类模型有重要的参考意义.2.缺点将机房内流体为不可压缩牛顿流体,而实际情况下气体是可以压缩的,因此,模型可能会存在误差;模型优缺点分析五.模型评价与推广将机房看成封闭的系统,
34、并忽略了其他辅助设备（如PDU等）的影响;由于机房空调采用上部进风方式,导致靠近空调上部空间的气流聚集,不易散热,温度较高,易形成 热点,不利于绿色机房的设计;机房是全封闭式,进风口设在空调上部,机房内内气流的回风只能通过机房房顶空间回风进风口,造成气流组织不是很好.可以将地域与季节因素考虑进去,对于不同的地域可以采取一些相应的改进措施,如可以加装一些空气干燥设备或加装新风换气系统,以减少空调制冷时的额外功耗;机房模型推广五.模型评价与推广在实际应用中可以先考虑不同服务器,空调等的功耗和发热量的不同,对模型进行修正和微调,以增大结果的可靠性和可信性.在实际的推广应用中,可将人的因素和室内外热交换考虑进去;对于采取其他结构布置的机房,可以在Airpak模型搭建的过程中加入相应的仪器设备模型（blocks,vents,fans等）使模型能够适应更多条件;讲解者:李鑫组员:熊钊 于壮 李鑫(吉林大学)基于基于Airpak的绿色机房模型仿真评价与控制的绿色机房模型仿真评价与控制
