2012深圳杯数学建模夏令营C题解答Word下载.docx

1、全国评阅编号（由全国组委会评阅前进行编号）：3仿真机房建模摘要：本文主要研究根据机房的基础设施状态，按照行业规范要求合理地26布置机柜，分布任务，尽量避免局部地区过热问题，我们根据题目所给定的一些数据和信息分析建立了扩散传播模型、采用了线性回归方法、并基于k e两方程紊流模型建立数学模型。首先，本文将绿色机房设计聚焦于机房内热环境分析，然后我们根据附件1 给出的数据，通过采用 MATLAB 软件绘制出机房温度和流场分布，之后求出温度的极值，从而求出温度最高点。其次，在确定该机房热分布时，考虑到机柜的对称性，本文将每个机柜作为一个热源，建立扩散传播模型，将热量扩散模拟为不同浓度空气在机房进行流通

2、，通过冷通道散发的冷气流与机柜散发的热气流接触，使得机房不同方位的温度不同，根据最后推导的空气浓度表达式，基于任务量，得出整体机房的热分布情况。再次，经过对题目的分析，基于上一部分所得模型和附件 2 中的数据，在充分考虑机房对称性的基础上，选取任务量对称的温度，采用 matlab 软件计算出温度极值点，再整合机柜任务量对温度极值点的权向量，得出机柜任务量与温度的相关度，再通过线性回归，最后得到最低温度的分配方案。最后，模型四以机房空气流组织为研究对象，采用k e两方程紊流模型建立数学模型，对机房的冷通道出口风速与机房室内的温度的关系进行了详细的动态描述，使得气流动态化，从而总结出通过控制冷通道

3、出风口风速来调节室内的温度，达到电子信息系统机房设计规范C 级要求。在本次模型建立中，我们得出以下结论：1. 热量在温度极值点会发生集中现象2. 在机房对称的情况下，温度分布还是会收到各种干扰，从而导致温度分布不均3. 当机房开机时应该加大冷通道的通风速度，关机的时候降低冷通道的通风速度关键词：热分布机房气流扩散传播一、问题重述1.1 问题背景由于高密度计算、多任务计算的需要，越来越多的高性能数据中心或互联网中心（DC、IDC）正逐渐建成。在现代的数据中心内，由于刀片服务器成本与性价比高，体积小而被广泛使用。由于自身能源与冷却条件限制，这类大规模的数据中心或许每年需要花费数百万美元，主要用于计

4、算设备及系统冷却所需的能源费用。因此有必要提高数据中心设备的能效，极大化数据中心的能源利用率及计算能力。大约在上世纪 90 年代后期，IBM、HP 等公司首先提出绿色数据中心的概念，并受到世界各国的广泛重视。1.2 提出问题（1） 根据附件 1 的数据，绘出冷、热通道的热分布及流场分布及室内最高温度位置。（2） 建立描述该问题热分布的数学模型及算法，并与测试案例进行比较。（3） 如果定义该机房的总体任务量为 1，根据你的模型及附件 1 的流场数据，确定服务器实际任务量为 0.8 及 0.5 的最优任务分配方案，并给出室内最高温度。（4） 如果按照电子信息系统机房设计规范（附件 3）C 级要求控

5、制机房温度，讨论服务器设计任务量一定条件下，如何控制空调的送风速度或送风温度（可以通过送风槽的出口风速与温度来描述）。二、问题分析（1） 由于整个机房内不同方位的温度是有差异的，根据附件 1 的数据，采用 MATLAB 软件绘制出冷热通道的热分布图与流场分布图，并根据热分布图得出了机房室内的最高温度。（2） 由于机柜散发热量及冷通道的气体传播的途径有很多，本题忽略其他因素，只考虑在空气中传播，建立大气扩散模型，通过与测试案例得出机房室内的热分布情况。（4） 根据电子信息系统机房设计规范C 级要求，在一定任务量的基础上，以机房内的气流组织为研究对象，采用k e两方程紊流模型建立数学模型， 并得出

6、了空调如何控制风速来调节机房室内的温度。三、问题假设（1） 假设每个机柜散发的热量是相通的，且热量分布稳定。（2） 假设机柜散发热量是以相同溶度的热空气在机房传播。（3） 假设机柜散发热与冷通道的冷气在机房的分布是不受机房风速影响。（4） 假设机房室内气体为不可压缩流体，且满足 Boussinesq 假设：认为流体密度的变化仅对浮升力产生影响。（5） 假定流场具有高的紊流 Re 数，流体的紊流粘性具有各向同性。（6） 气流为低速流动，可忽略由流体粘性力做功所引起的耗散热。（7） 假设一点的温度只与一个机柜相关。（8） 假设温度与机柜任务量满足指数函数关系。（9） 假设温度最高点只发生在温度极值

7、点上四、符号说明ki （i = x，y，z）：机房内任意一点的扩散系数C（x，y，z，t）：机房内任意一点热空气的浓度Ux ，Uy，Uz：为 x，y，z 方向风速Kx，Ky，Kz：为 x，y，z 方向上的扩散系数p0 ：单位时间内机柜散发热空气的速率Q0 ：从机柜散发出的热流总量Q1 ：在（t，t + Dt ）内通过某空间域的热流量Q2 ：在某空间域热流量的增量x ， y ， z ：分别为用浓度标准差表示的 x，y，z 轴上热量的扩散参数vs ：热空气与冷空气接触的系数；S ：机柜散发热的扩散的系数K 、t ：分别为风速与任意扩散时刻H ：机柜散发出的热源距地面的高度G：为重力加速度；K：为流

8、体紊流动能；P：为时均压力；C：为热源强度；Pr：为充分紊流时的普朗特数；T：为流体温度；T0：为参考温度；Ui：为速度分量，当 i=1、2、3 时分别代表 X、Y、Z 方向；：为紊流能量耗散率；m、s：分别为层流和紊流的动力粘性系数；：为流体密度；：为流体体积膨胀系数。五、模型建立及求解5.1 模型一经过对附件 1 的数据和题目的分析，直接使用 MATLAB 软件将该机房的二维热分布及二维的流场分布图绘制如下图 1 由附件 1 得到的通道二、三热分布图和流场分布图：图 1.1 通道二热分布图图 1.2 通道二流场分布图图 1.4 通道三热分布图图 1.4 通道三流场分布图图 1 由附件 1

9、得到的通道二、三热分布图和流场分布图同时，通过 MATLAB 求出极值（MATLAB 计算程序见附录）： 通道二：Imregionalmax = 41.3707Points =3.36972.7000通道三：Imregionalmax = 47.85713.36972.40915.2 模型二机房的热分布模型记热量开始扩散的时刻t = 0 ，并且以热源作为扩散中心，则可将 t 时刻机房内任意一点（x，y，z）的热空气浓度记为C（x，y，z，t）。假设单位时间内通过单位法向量面积的热流量与浓度梯度成正比，则有q = -ki gradC其中ki （i = x，y，z）是扩散系数，grad 表示浓度梯

10、度。（1）假设机房空间W 的体积为 V，包括机房空间W 的曲面为一规则的球面，设其表面面积 S，外法线向量为n = - x ，- y ， ，表1zz则在（t，t + Dt）内流通过空间域W 的热流量可以表示为：= t +DtQ1t q ndsdts（2）在空间W 所包围的区域内空气中的热流量增量可表示为：Q2 = C（x，y，z，t + Dt ）- C（x，y，z，t）dVV（3）t+Dt由机柜散发的的总热流量为： Q0 = t p0 dVdtW（4）由质量守恒可得出： Q0 = Q1 + Q2 ，即，C（x, y, z, t + Dt ）- C（x, y, z, t ）dV +t +Dt t

11、 ndsdt =p0 dVdt（5）VsW根据曲面积分 Gauss 公式得： nds= divqdVsV（6）则式子（5）可以转换成C（x，y，z，t + Dt ）- C（x，y，z，t） D Dtt dV+ 由于：t +Dt tdivqdVdt = t（7） （+ D ）（）t +Dt kdiv（gradC ）dtC = lim Cx，y，z，tt - C= lim t（8）tDt 0DtDt 0Dt故式子（7）即可转换成： t C dV Dt +divqdV Dt =p0 dV Dt（9）VV化简得：dV +WdivqdV = p（10）t 0根据 A.Fick 扩散微分方程式中：C +CC 2VtU x X+ U y Y+ U z Z= KXX 2 + KY Y 2 + KZ Z 2（11）其中：C 为常温气体浓度；t 为时间；Ux，Uy，Uz 为 x，y，z 方向风速；Kx，Ky，Kz 为 x，y，z 方向上的扩散系数。假设机房热量是在无风条件下扩散