用零方程湍流模型模拟通风空调室内的空气流动

1、用零方程湍流模型模拟通风空调室内的空气流动simulationofindoorairflowinventilatedroombyzero- equa tionturbulencemodel摘要：利用带浮升力效应的湍流模型和一个新零 方程湍流模型对某房间内空气的混合对流流动进行了数值 模拟，通过比较发现，新零方程湍流模型与实验数据吻合 得更好，且可以很快获得收敛解。利用新零方程湍流模型 对房间内的等温流动、非等温流动进行了模拟，发现数值 计算结果和相应实验数据吻合得很好。由此考察和验证了 新模型对暖通空调领域中流动和传热问题的实用性和可靠 性，可以利用该模型快速、精确地设计和分析暖通空调领 域中

2、的室内空气流动问题。关键词：室内空气流动；数值模拟；湍流模型；零方 程abstractmi xedconvect ionflowina ventilated roomwassim ulatedusin gthek- £ mod elandanewz ero- equati onturbulen cemode1co mparisonof thenumeric alresultsa, ndthemeasu reddatasho wedthatthe newzero- eq uationturb ulencemode lgivesmore satisfacto ryresultsi nsh

3、orterti me. theisot hermalandn on- isother malairflow inanair-co nditionedr oomwerethensimulatedwiththezero-equationturbulencemode1.thesimulatedresuitsagreedwellwiththeexperiments.therefore,thenewzero-equationturbulencemodelisrecommendedforheating,ventilatingandairconditioningindustrybecauseitisquic

4、kandsoefficientlyaccurateforengineeringstudies .flow;numericalsimulakeyword sindoorairtion;turbu lencemodel :zeroequa tionmodell 974 年nielse n首先将计算流体动力学cfd(compu tationalf 1 uiddynamic s)技术应用于暖通空调工程领域，如今,可以利用cfd技术模拟预测空调房间内的空气流动,进行气流组织 设计与分析但是，在应用中也存在一些问题，如何快速、准 确地在模拟预测工程中需要优化比较的大量工况是其中最 为迫切的一个问题这主要取

5、决于湍流模型的选择.由于通 风空调房间室内空气流动多为湍流流动，而人们对湍流机理 的认识还不全面，故对于工程应用，目前多为半经验的、唯 象的模拟。限于目前的计算机能力，工程中最常采用的是 涡粘系数模型evm (eddyvi scositymod els)中的k- £两方 程模型或其变形。但是k- £模型对于等温流动情形能模拟 得很好，对于空调通风房间内的非等温，混合对流流动却 有较大的误差。nielsen于1 998年指出，对于热羽流， 贴壁射流，温度分层流动等需要不同的湍流模型进行模拟 方能取得满意的精度3,而这些流动形式往往并存于实际 上的空调房间，采用标准的k-

6、63;模型势必导致很大的误差。 而k- e模型的变形，如低雷诺数k- e模型可以取得较好的 精度，但是其计算量已相当于复杂二阶封闭模型或大涡模 拟的计算量，无法被工程界所接受。近年来，为满足暖通空调工程应用对数值模拟快速、 准确的需求，一些学者提出用零方程湍流模型对所关心的 问题进行模拟。由于湍流模型模拟是唯象的，半经验的， 故尽管零方程模型比较简单，但在专门的领域内却能获得 比复杂模型更符合实际的结果。本文介绍一个新的零方程 模型，并通过在室内空气流动数值模拟说明新模型在工程 中的应用。1两种湍流模型描述室内空气流动密度变化不大，通常采用boussin esq假 设4。文中所用两种湍流模型均

7、为evm,这是基于b oussinesq关于雷诺应力假设的湍流模型。基于以上假设， 可得室内空气湍流流动的控制方程，参见文5。evm的核 心是求解湍流动力粘度u t .本文讨论的重点为m it建筑技 术系开发的新零方程模型（下简称mit零方程模型）。由于 k- £模型是工程中常用的湍流模型，这里也利用它对本文研究的室内空气混合对流算 例进行计算。这两种湍流模型可大致描述如下： 带浮升力效应的模型5模型属两方程模型，它引入湍流动能k和湍流动 能耗散率e表示湍流粘性系数，u t =cd p 2k/ e (1 )其中：cd二，为常数，p为空气密度。考虑了浮升力影 响的k和e的微分方程可参考

8、文5。2) mit零方程模型这是在室内空气自然对流和混合对流的直接数值模拟 dns (dir ectlynumer icalsimula tion)结果的基础上提出 的湍流模型6 ，该模型针对房间内非等温流动的rayl eigh数范围(x1010),认为涡粘系数正比于流体密度、 当地速度和距壁面最近之距离，比例系数由直接数值模拟 的结果拟合而得u t=74 p vl (2)其中：v为当地时均速度，1为当地距壁面最近的距 离。该模型少求解2个微分方程，而仅求解关于质量、动 量和能量守恒的5个微分方程，故计算最省时间。2数值计算方法和工具简述本文采用清华大学建筑技术科学系建筑环境与设备研 究所开发

9、的sta ch-3进行模拟。它采用有限容积法将计算 区域离散为不均匀网格，差分格式采用混合模式，算法为 simpl e算法，动量方程釆用交错网格，边界条件采用壁面 函数法处理。3混合对流算例比较为比较mit零方程湍流模型和常用的模型对室内 空气流动数值模拟的性能，选择z hangj. s.等人做的室内 混合对流实验作为验证算例7。图1所示为实验工况，其 中w和ii分别为小室宽度和高度。从房间上方送风速度vd 为/s送入温度td为。c的空气，室内地面温度tf为。c,这 是一个强迫对流和自然对流并存的混合对流的典型例子， 有关参数见表lo其中：tc为出风温度，at fd=tf-td, 为送风与地面

10、温差，r ed为入口雷诺数，a rd为入口阿基 米得数，q为换气次数。图1实验工况示意图 表1实验条件参数表vd / (m s-1) td /°cte/°ctf/°c a tfd/°creda rdq/h- 1573 5气流组织实验在一个x x的小室中进行。实验数据 通过烟气流线方法获得。整个空气流型由一个贴附吊顶的 弯曲的自由空气射流，顶射流，和一个大的再循环漩涡体 现(见图2)图2实验的所得流型图采用两种湍流模型计算所得流型图如图3所示，与图2 的实验所得流型对比发现：两种模型都能将房间中心的再 回流漩涡旋很好的模拟出来。但是k-£模型不能

11、将房间右 下方的局部回流涡旋反映出来，而mit零方程模型比k-£ 模型更为准确地模拟出该混合对流的流型。进一步比较 x/w二和x/w二处温度、速度分别沿高度的变化发现（见图4 和图5） : mit零方程模型所得各点速度比模型所得结 果与实验数据更为接近。尤其是对于x/w二的位置（图4 b）, mit零方程模型模拟的速度变化趋势与实测值一致，在y/ h二处速度最大，从而正确模拟岀射流中心速度最大的射流 物理特性。总之，两个位置温度计算值都比实测值小（图 5）,这是因为计算中没有考虑辐射作用的缘故。尽管如此, 就两个位置的温度变化趋势而言，mit零方程模型与实测趋 势更吻合，且其计算值与

12、实验数据更接近。零方程模拟所 得房间中心涡旋位置偏差较大，说明模型也存在不足。图3计算所得流型图图4速度随高度分布 图5温度随高度分布由以上比较可见，对于室内空气混合对流流动这种比 较复杂的流动形式，采用模型并不能取得满意的结果, 甚至与实验数据不符。而mit零方程模型却能取得令人满 意的结果，尤其是采用该模型能更快地获得收敛结果，这 对于工程应用中大量的三维计算工况有着很大的实用价 值，因为在暖通空调房间气流组织设计中，往往需要对很 多三维工况进行模拟比较以获得最优化的设计。对于同样 的网格数，在同一台pii 1500, 128m内存微机上采用两种 模型计算达到收敛时间的比较结果见表2,可以

13、看出采用 mit零方程模型比模型快约4倍，且可取得满意的结 果。表2两种模型计算时间比较湍流模型网格划分计算至收敛时间/mink- £ 37 x32x 330mit 零方程 37x32x3 74mit零方程湍流模型的应用采用m it零方程模型可快速双不失工程上精度地对室 内空气流动进行数值模拟。为了考察该湍流模型在室内空 气流动中的实用性，此处再给出几个不同类型的算例，以 说明该模型在室内空气流动数值模拟中的优越之处。计算 工具为上述的清华大学建筑技术科学系自行开发的st ach-3o 4. 1室内等温流动采用的算例和实验资料是19 90年nielsen等人对图6 所示的房间进行的测

14、试数据8。该房间尺寸为：高度h二， 长高比l/h=,宽高比w/h二，送风口高为h, h/h=,回风口 高为t, t/h=,送风速度v0取为取为/s,水平入流。此处将 采有用x/h二和x/h二两个断面上，时均速度v沿高度方向y 的分布来进行对比和分析。图6实验房间示意图图7为上述两个断面上的实验结果和计算结果的对 比，由图可以看出，对于等温流动，mit零方程模型仍然能 取得和实验数据吻合得很好的模拟结果，并且，对于本次 计算所用网格数（37 x20x3）,在pi ii500, 128m内存微 机上只需3mi n即可获得收敛结果。图7速度沿高度方向分布4. 2室内非等温流动文9为了验证房间通风情况

15、数值计算的结果，进行了 如下实验：在一间长x宽x高为x x的小屋内非等温送风, 采用顶送风，送风口为宽度的条缝风口，回风口位于小屋 右下角，高为。斜向下4 5°送风，送风量为3 /s,送风温 度为。c,回风温度为22°c,送、回风温度为。c,室内由电 加热器模拟均匀分布热负荷19w/m3 o房间结构如图8所示。图8实验及计算用房间示意图通常人们关心的只是工作区的温度和速度值，故实验 所测点为房间正中(长度方向一半处，x=)高度依次为，, 处的温度和速度值。利用前述模型进行计算，所得结果如 图9所示。由图9可见，各实验点与模拟值的速度差值均 在/s以下；温度相差也很小，最大差

16、值仅为。c。说明计算 值和测量值吻合较好。在前述计算机上，对于本次计算采 用的网格数(24x16x14),大约计算10 min即可收敛。图9实验与计算值对比5讨论通过以上分析得知，采用mit零方程模型对室内空气 流动数值模拟可以获得和实验数据吻合得很好的结果。在 cfd发展初期，限于当时计算机技术的水平，对于湍流的模拟多 采用简单的湍流模型，如普朗特混合长度模型等。随着计 算机技术的不断发展以及对湍流认识的提高，人们开始釆 用一些复杂的湍流模型，如常用的k-£模型等。但是实践 表明，k- e模型对很多问题仍不能获得满意的模拟结果，故出现了更为复杂的湍流模型，如reynolds应力模型

17、（rsm）、代数应力模型（asm）等高阶封闭的模型。这些模 型需要求解的微分方程均在10个以上，对于工程中需要求 解的复杂，三维耗时太多，无法满足工程应用快速的需要, 而且，对于某些特定问题，这些复杂模型的模拟结果并不 比简单模型模拟得效果好。于是，近年来零方程等简单 模型又为众多学者所重视。需要指出的是，这些简单模型是建立在高级的湍流数 值模拟技术基础上的，如上述的mit零方程模型即是借助 直接数值模拟的结果提出的，这正表明了否定之否定的 哲学原理。当然，由于零方程模型本身的局限性，它在特 定领域内的适用性还需要在实践中接受进一步的检验。6结论1）对于室内空气自然对流和强迫对流共存的混合对流

18、 流动，釆用mit零方程模型能获得比带浮升力效应的k-e 模型更为准确的结果。2）mit零方程模型对室内等温和非等温流动能快速获 得模拟结果，并能保证一定准确度，可以用于暖通空调工 程室内空气流动的数值模拟，指导设计；3）针对性强的零方程湍流模型对特定问题能比复杂湍 流模型取得更好的结果，这可满足工程应用快速高效的要 求，但其适用性需要在实践中进一步检验。参考文献(referen ces)i c henqingyan , mosera, hu berapredi ctionofbuo yant, turbu lentf lowby alow-renol ds-numbe rk- £

19、model jashraetransaction.1990. 96(1) :564-57.2 chenqin gyan, xuwei ran. azero- equationtu rbulencemo delforindo orairflows imulationj. energ yandbuildi ng, 1998, 28 :137-144.3 nielse npv.thesel ectionoftu rbulencemo delsforpre dictionroo mairflowj . ashreatransactio ns, 1998, 10 4(1) ; 1119-11264 t rittondj.p hysicalflu iddynamics mclaren do