2012年“深圳杯”全国大学生数学建模夏令营C题北京大学.pptxVIP

绿色数据中心温度分布及流场的分析邱铭达周沛劼吕一鼎北京大学数学科学学院问题的分析与概括如何在给定了机房中一些位置处的温度信息的情况下,描绘机房整体的温度分布情况?如何对于机房内温度和流场进行数学建模?如何在模型的基础上,寻找在一定任务负载的情况下的优化机房温度的方法?第一部分:描绘温度分布基本思路:进行三维插值工具:GNU OctaveHermite Interpolation第二部分：机房的建模提要：抽象与模拟对机房的物理抽象房间内的物体室内空气流动的流体力学问题和框架使用Computational Fluid Dynamic(CFD)软件 Ansys Fluent? 模拟流场和热传递的过程机房的物理抽象房间本身（封闭空间）空凋送风口（冷空气入口）Gauge PressureTemperature空调出风口（排气口）机器（冷风进入机器，被加热后再被风扇排出）进一步细化：4组机柜，每组8个机柜，每个机柜5个机架，共160台服务器。服务器的模型单台服务器的管道模型单台服务器被理想化为一个管道，管道中有一个风扇(fan)和一个散热器(radiator)冷空气借助风扇从冷通道进入管道，通过散热器并被加热，再排入热通道散热器的物理属性：Loss Coefficient （与本问题关系不大）Heat Flux表征服务器的发热能力为什么用Heat Flux表征服务器发热？服务器向房间中的排热是依靠radiator来完成的热流密度仅与两个因素相关——热传导系数、温差温差大，散热快；温差小，散热慢，最终达到平衡态，服务器放热与散热相同单台服务器的管道模型Gambit: 机房的构建机房内空气流动的物理分析框架流体力学的基本控制方程考虑三维直角坐标系，社流体的速度矢量 在三个坐标上的分量分别是u，v，w，压力为p，流体密度为 。这里的四个量都是空间坐标及时间的函数。根据质量守恒、动量守恒以及能量守恒定律，我们得到如下的流体力学控制方程：连续型方程：动量方程：能量方程：可压缩流的控制方程由于我们考虑的问题中有比较明显的压强变化（数百帕斯卡），所以我们决定选用可压缩流而非不可压缩流来解决问题有以上基本控制方程我们可以推导出可压缩流体的湍流方程：其中，连续性方程：雷诺方程：标量输运方程：简化的湍流模型：双方程k-? 模型理论方程在计算上难以实现考虑k-? 模型原因：使用最广泛、最有效；不需要就不同情况修正模型参数；对于限制流（即有壁面约束）特别有效。不使用其他简化模型原因： ——混合长度模型不适合在湍流量的对流和扩散不可忽略（如有回流）的场合使用； ——雷诺应力模型过于复杂，收敛过慢双方程k-? 模型概述：首先定义如下涡粘性(是无量纲的常数)：之后将标准模型应用于k、 ? 的输运方程有以下两式：k-? 模型的边界条件k-? 模型方程由于有梯度扩散而成为椭圆型方程，其特性与其他椭圆型方程相同，需要下列边界条件：入口：给定 k、? 的分布出口或者对称轴满足：自由流： k=0；?=0在探索性设计计算中，如果没有详细的边界信息，CFD工程用户一般不需要自己测量k和? 这两个量；利用文献资料中的k和? 就可以使得工作进行下去。在具体的计算中，由于其是非常复杂的偏微分方程，所以我们使用了计算流体力学（CFD）软件fluent来实现。模拟工具：fluent 6.3.26（1）生成计算网格；（2）输入并检查网格；（3）选择求解器；（4）选择求解的模型；（5）确定流体的材料物性；（6）确定边界类型及其边界条件；（7）调节计算控制参数；（8）流场初始化；（9）求解计算模拟效果模拟效果基本结论：图像定性分析冷通道中，温度随着高度的增大而升高热通道中，有热点的存在直观物理解释流场模拟的结果分析和评价从图像上看基本满足附件一的定性要求： ——冷通道，越靠后面温度越高且高处温度更高。 ——热通道，温度数值符合较好，且存在热点，位置基本符合模型的特点：用管道模型简化了服务期内部的散热过程计算速度相对较快，结果定性正确，定量基本吻合第三部分：机房温度的优化问题的细化已知机房的任务分配和送风口的冷风温度，怎样选取合理的冷风的出口速度？已知机房的任务分配和送风口的冷风速度，怎样选取合理的冷风的出口温度？已知机房的任务分配,怎样选取合理的冷风出口速度与出口温度的组合？总任务量一定的时候，散热方式不变，怎样合理分配任务量？问题初探：参数条件的选取T\P300pa600pa900pa1200pa1500pa1800pa2100pa279K299.9216293.3067291.1075290.379289.453288.9865288.6153282K303.0908296.3723294.089293.0737292.4359291.971291.5789285K306.1105299.5