北京京创鑫业分散式空调室内流动及其传热过程的数值模拟运算.docVIP

分散式空调室内流动及传热过程的数值模拟 西安交通大学AirTS高大空间专用空调系列产品系统地解决了高大空间采暖、制冷、通风和热能回收等多种需求。适用于工业厂房、物流中心仓库、飞机及机车维修中心、体育馆、展览馆、汽车4S店、超市等4-33米高的建筑场所，技术成熟，安装简单，控制灵活，可靠性高，效能卓越。 1物理模型 1.1 几何结构 本文模拟的对象是一个高度为12米的车间，长和宽均取10米作为计算区域；采用的是空调，布置示意结构如图1所示。为了方便计算，在不影响计算结果的前提下对空调模型进行了合理的简化：在模型中只对空调的进出风口进行计算，将空调的出风口作为系统的入口边界条件，空调的入风口作为系统的出口边界条件，空调叶片角度从0-60°可调。本文关注的重点是车间内部的空气的流动及热交换过程，故而未考虑空调内部的流动及换热。计算中所采用的数据，均为，如表1所示。表1 参数 单位 数量 空调出口直径 mm 630 空调送风量 m3/h 970 空调风压压降 Pa 300 空调送风温度 ℃ 30 空调回风温度 ℃ 20 图1 计算区域示意图 1.2 边界条件 (1) 入口边界条件 入口采用速度入口边界条件，通过给定的空调送风量以及出口直径计算得出入口的速度为0.86m/s，入口空气温度为30℃。 (2) 出口边界条件 出口采用排气扇边界条件，将压力参考坐标原点设在入口面的中心处，此处的压力假定为1个标准大气压，依据给定的空调风压压降设计值300Pa，设定此处的压降为300Pa；回风温度设为20℃。 (3) 壁面边界条件 计算中考虑了壁面对流动的影响，对近壁面采用标准壁面；着重考虑各个不同壁面的换热条件：侧壁及屋顶既考虑其外壁面与外界空气之间的对流换热，又考虑其与天空之间的辐射换热，对流换热系数及辐射换热系数均根据一般气象条件依据手册查得Wall -9 24.36 0.9 -9 0.9 侧墙 Wall -9 24.36 0.9 -9 0.9 地面 Wall -9 0.9 图2 车间热物理模型 2数学模型 2.1 区域离散 数值计算的第一步是计算区域的离散，即对离散区域划分网格有限体积法（Finite Volume Method，简称FVM）是近年来发展最为迅速的一种离散化方法，目前在计算流体动力学领域得到了广泛应用本文采用就是有限体积法。有限体积法的区域离散实施过程为：把所计算的区域划分成多个互不重叠的子区域，即计算网格（grid），然后确定每个子区域中的节点位置及该节点所代表的控制体积。区域离散化过程结束后，可以得到节点、控制体积、界面和网格线四种几何要素。在划分网格的时候为了保证计算的精确度以及计算的速度对于三维区域应该尽量使用规则化的六面体网格。本文在进行网格的划分时使用的Cooper方法，将车间的上下两个平面作为源面划分六面体网格，提高了计算的准确度，减少了计算时间。在此基础上进行试算对网格的独立性的解进行分析从而确定最终采用的网格数量为12万。 图 网格划分示意图 2.2 控制方程 计算中考虑的是车间内部的三维稳态流动及换热过程，在区域离散的基础上对将每一个网格作为一个控制体，对进出控制体流体的流动及热交换进行分析并建立方程组：整个过程中需要建立遵循质量守恒的连续方程；遵循动量守恒的动量守恒方程以及遵循能量守恒定律的能量守恒方程；针对气体的湍流流动，使用工程上普遍适用的κ-ε双方程模型来推导动量方程中未知的湍流粘性，从而在数学上封闭了该方程组。 （1）连续方程 (公式1) 其中： ρ — 气相流体的密度，kg/m3； — 气相速度，m/s。 （2）动量方程 (公式2) 其中： — 气相流体的应力张量，Pa； — 气相流体的升力, m/s2； — 气相流体的质量力, m/s2。 （3）能量方程 (公式3) 其中： — 温度，K — 热源 — 比热容，j/kg-k； — 有效导热系数,w/m-k。 （4）湍流模型κ方程： (公式4)ε方程： (公式5) 其中： ， ， ，， 2.3 数值离散 前面已经给出了计算区域离散化以及气相流体的控制方程，在进行动力学计算之前，还需要将控制方程在网格上离散，即将偏微分格式的控制方程转化为各个节点上的离散方程组，再在计算机上求解这一系列的代数方程组，从而得到每一个网格节点上的解，进而得出整个区域上场的分布。 在建立离散方程时，很重要的