利用CFD方法研究室内湿度和结露分布_0.docVIP

利用CFD方法研究室内湿度和结露分布 摘要：本文首先建立了在气流场，温度场之外，能够 分析 室内湿度分布，结露分布， 计算 结露量动态形成过程的CFD计算模型，并给出了自由水面热湿同时移动的计算新 方法 ；通过与模型试验结果的比较，验证了计算模型的精确性；最后通过具体算例，发现通风方式对室内湿度分布和结露的形成 影响 很大，在室内湿度较大的场合，可以通过通风方式的优化设计更有效地解决除湿和结露 问题 。 关键词：CFD模拟 结露 湿度 通风 1.前言 湿度环境问题不同于其他室内污染问题，湿度过高或过低都将严重影响建筑物性能和居住者健康。湿度过低，人会产生干燥等不适感，引起墙体裂缝，木制板材变形，另外据北欧学者的 研究 ，流感病菌在低湿度下生存率明显增加；而湿度过高，一方面造成墙体表面及内部结露，降低墙体断热和耐久性，影响建筑物寿命[1]，另一方面，当湿度超过70%时，会带来霉菌（Fungi）的大量繁殖，引起过敏性皮炎，哮喘等疾病，影响居住者的身体健康[2-5]。我国地域辽阔，不同地区室内湿度环境呈现出不同的特点。因此，深入地研究室内湿度环境问题具有重要的意义。 随着计算机功能的飞速 发展 ，CFD仿真技术在建筑环境领域，如室内温度场分布，换气效率，人体周边微环境等研究中日益得到 应用 。但是到 目前 为止，国内外利用CFD技术分析研究室内湿度分布 规律 的还很少[6，7]，而以此研究室内结露的形成和发展的还没有。 2.考虑湿度和结露计算的CFD修正模型 本研究中采用标准k-ε紊流模型。但考虑到水蒸气含量较大时空气密度会出现变化，从而影响浮力的计算，依据近藤等提出的方法[8]，对模型进行了修正，在浮力项中导入βx。修正的湿结露CFD计算模型见表1。另外，本研究中还考虑了结露量的计算。由于结露的形成是一个动态的过程，提出了两个随时间变化的指标：至时间t，1）单位壁面积上的结露量CON(s,t)；2)壁面 总结 露量SUMCON(t)。计算方法及其与CFD模型的结合见图1。 另外，建筑内的湿源，如浴室和厨房等，由于自由水面面积较大，水温一般高于周围空气温度较多，水蒸气分子在扩散过程中还同时伴随有热的交换。如果不考虑这种热湿的同时传递将会给室内温湿度分布和气流计算结果带来较大的误差。由于这方面的研究较少[9]，我们通过试验发现自由水面的热湿传递量m和qm由以下公式确定，然后可以向上的热湿流束的形式作为内部边界条件代入CFD计算中： 图1 湿结露CFD计算模型中的结露计算流程 模型的计算公式 表1 连续方程： 运动方程： 输送方程： 输送方程： 熱输送方程： 水蒸气输送方程： 式中： 涡粘性系数等相关项： 　　　　（1） 　　　　（2） 式中 —室内换气次数，h-1； —对应于水温的饱和蒸汽压力，kpa； —周围空气的蒸汽分压力，kpa； —水温，； —周围空气温度，； —水蒸气的蒸发潜热，kJ/kg； —空气的定压比热，kJ/(kgK)。 —新提出的参数，它代表由水蒸气扩散引起的实际散热量与水蒸气全热交换时的最大散热量的比值。事实上，由于一部分水蒸气从水面蒸发的途中只进行了显热交换，没有发生相变，应是一个0-1之间的数值。利用简单的公式推导和试验拟合，它可以整理成下式。具体 内容 可 参考 文献 [10]。 　　　　（3） 3.模型试验 为验证所建湿结露CFD计算模型的正确性，在日本东京燃气公司技术研究所的人工气象室建立了模型小室，进行了模型试验，并比较了模拟结果与实验数据。 模型的概要见图2。模型小室由聚乙烯板制成，小室中通过地面上的加湿器和水的温控来调节和模拟各种加湿工况，加湿量由 电子 天平测量加湿器的重量变化求得。小室外壁上部和下部分别设开口，安装有小型轴流风机的通风短管可连接在开口上，可利用风机位置的变化来模拟不同的通风方式。小室内的风量由通风短管内安置的微风速仪测定。除了小室外部的温湿度外，在小室中心断面上布置了14个温湿度测定元件（THP-B4，日本神荣公司）进行温湿度分布的测定。 试验与模拟工况见表2。其中工况1为验证室内温湿度分布的稳态计算，工况2为验证结露形成与发展的动态计算。通风方式均为下送上排的机械排风方式。 试验工况 表2 工况编号 气象室条件 加湿器水温 （） 加湿量 （g/h） 通风量 （m3/h） 进风温度 温度（） 湿度（%） 1 14.3 47.0 44.6 27.5 17.0 14.3 2 19.4 45.0 68.1 101.3 7.9 25.2 图2 模型试验的概要与测点布置 4. 试验与计算结果的比较 4.1 温湿度分布的验证 图3给出了工况1的试验与CFD计算比较，其中试验数据为所有测点的测量值均达到稳定状态时的结果。模拟值中，qm=486W/m2为利用