【2019年整理】数据中心的热场模拟分析

题目：数据中心的热场模拟分析摘要：影响数据中心空调房间气流组织的主要因素是入口风速、进风口的位置、进回风口的相对位置等，本文采用 Fluent软件来模拟三维的较复杂的数据中心空调室内气流组织的流动形式，首先在Gambit里建立物理模型，然后对所建立的物理模型进行网格划分，选用空气流动的质量守恒方程、动量守恒方程、能量守恒方程以及 k-6方程，采用固定的入口风速、出口风速来对空调室内进行模拟分析，并在Fluent软件中以直观的方式表示出了各种不同气流组织方案下气流的温度场和速度场 (气流流型图)，分析讨论气流分布的规律、特点及其适用场合。结果证明，利用k-s湍流模型对空调房间的气流组织进行数值计算是可行的。这一数值模拟方法也能帮助设计人员预测其气流设计方案的效果。这对于保证良好的房间空调系统气流组织设计方案、提高室内空气品质、减少空调的能耗，对建筑物热环境的节能设计及室内空气质量的评价都有重要的意义。Influencetheairconditionroomofthedatacentercurrentofairorganizationofmainfactorisentrancewindvelocity,enterthepositionofplacewithadraught,enterreturntoplacewithadraughtofoppositepositionetc.,thisliterarygraceuseFluentsoftwareimitate3Dlyandmorecomplicationsofdatacenterairconditionindoorcurrentofairorganizationofflowform,firstintheGambitestablishmentphysicsmodel,thentoestablishmentofphysicsthemodelcarryonameshdividetheline,thequalitywhichchoosetouseaircurrentguardHengsquarethedistance,momentumguardHengsquarethedistance,energyguardHengsquaredistanceandk-'squaredistance,adoptionfixofentrancewindvelocity,exportwindvelocitycometoairconditionindoorcarryonimitateanalysis,andmeantbykeepingaviewintheFluentsoftwarevariousdissimilaritycurrentofairorganizationtheprojectdescendtemperaturefieldandspeedfield(thecurrentofairflowadiagram)ofcurrentofair,analysisregulation,characteristicsthatdiscusscurrentofair'sdistributeanditapplysituation.Resultcertificatethatmakeuseofswiftflowmodeltothecurrentofairoftheairconditionroomorganizationcarryonnumbercalculationisviable.Thisnumberimitateamethodalsoabilityhelpdesignpersonnelestimateitscurrentofairdesignprojectofeffect.Thisforassurancegoodsystemcurrentofairoftheroomairconditionorganizationdesignproject,exaltationindoorairquality,decreaseairconditionofabilityconsume,hottobuildingenvironmentofeconomyenergydesignandindoorairqualityofevaluationallhaveimportanceofmeaning.关键词：数据中心、气体流场、 CFDFluent软件问题重述：在一简化为长为10米，宽为9.5米，高为4米的数据中心机房内，在其机房正中间沿长度方向并排布置4排服务器机柜，每排由5个机柜组成。每个机柜为长1米，宽0.9米，高3米，每个机柜平均分成5层，每层称为一个机架，每个机架房一台服务器，每个服务器的功率为4千瓦（见图1）。由于服务器对工作环境要求比较高，特别是环境温度，通常数据中心机房需要对机房安装风冷系统及风排热系统。风冷系统将冷风（温度 13。）送到冷风道。风热系统将热风道的热风排到机房外。冷风道与热风道交错排列，冷风道宽 1.5米（靠墙边略小一些），地板是空隙率为50%勺空心砖，冷风从空心砖进入机房冷风道。热风道宽1米，热风通过热风道顶部的排风孔排出。机架的进风口面对冷风道， 出风口面对热风道，通过气流运动给服务器降温（见图2）。LQm你的任务如下：（1）建立该问题的机房空气流场模型，给出冷送风与热排风均为 0..400立方米/秒的室内流场。（2）建立该问题的热场分布模型，并在（ 1）条件下给出该问题的温度场。（3）如果外部环境温度为301机房的正常工作温度为不超过 30,计算该机房1天的制冷能耗。问题分析热通道冷通道热通道合理布置机柜对于确保机柜拥有适当温度和足够的空气非常重要。采取“面对面、背靠背”的机柜摆放及将空调设备与热通道对齐方式，这样在两排机柜的正面面对通道中间布置冷风出口，形成一个“冷通道”的冷空气区，冷空气流经设备后形成了热空气，再被排放到两排机柜背面中的“热通道”中，最后通过热通道上方布置的回风口回到空调系统，使整个机房气流、能量流流动通畅，不但提高了机房精密空调的利用率，而且还进一步提升了制冷效果。基本假设：在建立数学模型前对室内气体的流动先做以下假设1气体有无数无间隙的质点连续的组成2气体质点的物理参数是空间坐标和时间的连续函数模型建立求解我们采用k-8三维湍流模型方程，为了简化问题，做如下假设1室内气体满足牛顿内摩擦定律，为牛顿流体；2室内流体温度变化不大，密度可视为常数；3室内气体的流动形式为稳态紊流；4在紊流中心区，忽略能量方程中由于粘性作用而引起的能量耗散;5室内空气在房间内壁面上满足无滑移边界条件。6忽略人为地产热量，比如呼吸产热等气流流动的控制方程描述室内空气流动的三维N-S方程可由下列式表示(1)质量守恒方程：

(2)动量守恒方程:.t:Ui）•——（：UiUj）L 曲 -jij / — —：uiujj (2)动量守恒方程:.t:Ui）•——（：UiUj）L 曲 -jij / — —：uiujj 二xi3 二xt 一为(3)能量守恒方程:“i［、pr二i ；:Xi其中：P--密度(kg/m3)t--时间(s)xii--位置坐标(m)(i=1,2,3分别表示x,y,z三个方向)ui--xi方向的速度(m/s)邑--粘性系数(kg/(m.s))湍流模型k-8模型是个半经验公式，主要是基于湍流动能和扩散率。程是个由经验公式导出的方程。k-名模型假定流场完全是湍流，分子之间的粘性可以忽略。标准全是湍流的流场有效。K方程是个精确方程，名方k-8模型因而只对完标准k-⑤模型的方程湍流动能方程k和扩散方程s分别如下：C C C1/ Nt )乐； -一(%)+——(&Ui)=——尸+」!——+Gk+Gb-色-Ym+Ska Cxii %i艮 仃k j%■:、；22包Cdig7S式中：Gk一一由层流速度梯度而产生的湍流动能Gb一—由浮力产生的湍流动能；Ym一一在可压缩湍流中，过渡的扩散产生的波动C1,C2,C3,一—常量

Ok一—K方程的湍流普朗特数（无因次数）:外一一名方程的湍流普朗特数（无因次数）；Sk、Sr一一用户自定义的。湍流速度模型湍流速度从由下式确定k」t=：CLz其中C,是常量。模型常量如下：Cs=144,CZa=1.92,CV=0.09,Ck=1.0,6E=1.3。这些常量是从数据中得来的，包括空气、水的基本湍流。它们已经发现了怎样很好地处理墙壁束缚和自由剪切流。RNGk-g模型RNa-名模型是从暂态N-S方程中推出的，使用了一种叫 "RenormalizationGroup”的数学方法。解析性是由它直接从标准 k-s模型变来，还有其它的一些功能。RNGk-8模型的方程一（：k） 一（：k） ——「kUi）=——（：k」efft hi 二Xiik .『Gk Gb一一YM Sk；（：；）,—（：U）=—（：Jeff；（：；）,—（：U）=—（：Jefft :xi :xii ik二 二2京）C3Gbi2g『Rg SGk一一由层流速度梯度而产生的湍流动能Gb一—由浮力产生的湍流动能；YM一一在可压缩湍流中，过渡的扩散产生的波动C1,C2,C3,一—常量；^k一一K方程的湍流普朗特数（无因次数）:名方程的湍流普朗特数（无因次数）;Sk、Ss一一用户自定义的。(2)有效速度模型在RNGF^消除尺度的过程由以下方程:Cv：100方程4是一个完整的的方程，从中可以得到湍流变量怎样影响雷诺数，使得模型对低雷诺数和近壁流有更好的表现。在大雷诺数限制下方程得出：,、k2--t=)Cj—Cr=0.0845,来自RNG!论。这个彳1和标准k-名模型总的0.09很接近。在FLUENT^粘性的影响使用在方程的大雷诺数形式。当要计算低雷诺数是可以直接使用给出的方程。K-c模型中的模型湍流的产生在Gk项中，表现了湍流动能的产生，是按照标准 RNGI?旋流修正k-W模型而做的，从糟确的K方程，这项可以定义为：Gk=-口口一-k ij.Xi为了评估Gk和Boessinesq假设：Gk"S2S是系数，定义为：S=2S&k-8模型中受湍流7?力影响的k-名模型当重力和温度要出现在模拟中时，FLUENT^的k-鼻模型在K方程中考虑到了浮力的影响,相应的也在8方程中考虑了，浮力由下式给出：GbGbPrt：:xi这里：Prt一一湍流能量的普朗特数gi一一重力在i方向上的分量。对于标准和带旋流修正k-8模型，Prt的默认值是0.85.在RNG莫型里Prt=1/a，这里a对于标准和带旋流修正是由方程确定的，但是二.0是由方程确定的，但是二.0Prucp热膨胀系数P定义为:利用Fluent软件可以模拟出图像计算过程的残差图版I计尊较差限

Fig,4TCalculation Plct计算完成后，我们保存它的 case和data文件，然后可以通过后处理过程来显示和查看计算结果，我们截取x=1000的平面来显示出流场的温度、速度等数据的矢量图、等值线图，直观的表示出计算结果。19U*0G1期kMUhfOC1比1.加0C11・苏ECdoo1.口冲*的9321-26-JIE刘上川5丛21"Ml3hZM2加・。11O.OUt*DO船4-5工=l0QO平直的仝气速度分布的等flj线用值会JFiflk4-5Ai[V4locilyContourcfi=1000(n/s)212^001双“l»e@im，oo156e*OD14«6HX)1»e，001300121e*001UfOO10280092Se-OlB35e-0!741^01680OI556*014e5G-C139212»hC1187，C195222”23图4-6x=1000平面的空气速度分布的矢取用(B/s)Fig.4-6Airvelocityvectorofx=100C(n/s)(K)Lz3326023Xe<023383.27302326^023.2F023*83.2・83X0?Lz3M02330233023.1"0233023g0230^W303C2313023X02302^02303022W^Q22fi6e^22960*022432~22图4-7x=1000平面的空气温段分布的等值线图(K)Fig.4-7Airtenperaturccontourofx=1000(K)冷送风与热排风均为0..400立方米/秒的室内流场。可以模拟出该速度下的速度场和温度场：F油魂,3*，所用以巧期1靠，・陆「如fi*的।RRNS^1世配TEF3>C1e，蚓XI♦畸用.1iWfll"BMI2WH?1t跖Mlw蜘a也而叩冷送风与热排风均为 0..400立方米/秒的室内流场。32032＞更s.wct3tto