《利用CFD方法研究室内湿度和结露分布》由会员分享，可在线阅读，更多相关《利用CFD方法研究室内湿度和结露分布（8页珍藏版）》请在金锄头文库上搜索。
1、1利用 CFD 方法研究室内湿度和结露分布摘要：本文首先建立了在气流场，温度场之外，能够分析室内湿度分布，结露分布，计算结露量动态形成过程的 CFD 计算模型，并给出了自由水面热湿同时移动的计算新方法；通过与模型试验结果的比较，验证了计算模型的精确性；最后通过具体算例，发现通风方式对室内湿度分布和结露的形成影响很大，在室内湿度较大的场合，可以通过通风方式的优化设计更有效地解决除湿和结露问题。 关键词：CFD 模拟 结露 湿度 通风 1.前言湿度环境问题不同于其他室内污染问题，湿度过高或过低都将严重影响建筑物性能和居住者健康。湿度过低，人会产生干燥等不适感，引起墙体裂缝，木制板材变形，另外据北欧
2、学者的研究，流感病菌在低湿度下生存率明显增加；而湿度过高，一方面造成墙体表面及内部结露，降低墙体断热和耐久性，影响建筑物寿命1，另一方面，当湿度超过 70%时，会带来霉菌（Fungi）的大量繁殖，引起过敏性皮炎，哮喘等疾病，影响居住者的身体健康2-5。我国地域辽阔，不同地区室内湿度环境呈现出不同的特点。因此，深入地研究室内湿度环境问题具有重要的意义。随着计算机功能的飞速发展，CFD 仿真技术在建筑环境领域，如室内温度场分布，换气效率，人体周边微环境等研究中日益得到应用。但是到目前为止，国内外利用 CFD 技术分析研究室内湿度分布规律的还很少6，7，而以此研究室内结露的形成和发展的还没有。2.考
3、虑湿度和结露计算的 CFD 修正模型2本研究中采用标准 k- 紊流模型。但考虑到水蒸气含量较大时空气密度会出现变化，从而影响浮力的计算，依据近藤等提出的方法8，对模型进行了修正，在浮力项中导入 x。修正的湿结露 CFD 计算模型见表 1。另外，本研究中还考虑了结露量的计算。由于结露的形成是一个动态的过程，提出了两个随时间变化的指标：至时间 t，1）单位壁面积上的结露量 CON(s,t)；2)壁面总结露量SUMCON(t)。计算方法及其与 CFD 模型的结合见图 1。另外，建筑内的湿源，如浴室和厨房等，由于自由水面面积较大，水温一般高于周围空气温度较多，水蒸气分子在扩散过程中还同时伴随有热的交换
4、。如果不考虑这种热湿的同时传递将会给室内温湿度分布和气流计算结果带来较大的误差。由于这方面的研究较少9，我们通过试验发现自由水面的热湿传递量 m 和 qm 由以下公式确定，然后可以向上的热湿流束的形式作为内部边界条件代入 CFD 计算中：图 1 湿结露 CFD 计算模型中的结露计算流程模型的计算公式 表 1 连续方程： 运动方程： 输送方程： 3输送方程： 熱输送方程： 水蒸气输送方程： 式中： 涡粘性系数等相关项： （1）（2）式中 室内换气次数，h-1；对应于水温的饱和蒸汽压力，kpa；周围空气的蒸汽分压力，kpa；水温，；周围空气温度，；水蒸气的蒸发潜热，kJ/kg；空气的定压比热，kJ
5、/(kgK)。新提出的参数，它代表由水蒸气扩散引起的实际散热量与水蒸气全热交换时的4最大散热量的比值。事实上，由于一部分水蒸气从水面蒸发的途中只进行了显热交换，没有发生相变，应是一个 0-1 之间的数值。利用简单的公式推导和试验拟合，它可以整理成下式。具体内容可参考文献10。（3）3.模型试验为验证所建湿结露 CFD 计算模型的正确性，在日本东京燃气公司技术研究所的人工气象室建立了模型小室，进行了模型试验，并比较了模拟结果与实验数据。模型的概要见图 2。模型小室由聚乙烯板制成，小室中通过地面上的加湿器和水的温控来调节和模拟各种加湿工况，加湿量由电子天平测量加湿器的重量变化求得。小室外壁上部和下
6、部分别设开口，安装有小型轴流风机的通风短管可连接在开口上，可利用风机位置的变化来模拟不同的通风方式。小室内的风量由通风短管内安置的微风速仪测定。除了小室外部的温湿度外，在小室中心断面上布置了14 个温湿度测定元件（THP-B4，日本神荣公司）进行温湿度分布的测定。试验与模拟工况见表 2。其中工况 1 为验证室内温湿度分布的稳态计算，工况 2为验证结露形成与发展的动态计算。通风方式均为下送上排的机械排风方式。试验工况 表 2 工况编号5气象室条件加湿器水温（）加湿量（g/h）通风量（m3/h）进风温度温度（）湿度（%）114.347.044.6627.517.014.3219.445.068.1
7、101.37.925.2图 2 模型试验的概要与测点布置4. 试验与计算结果的比较4.1 温湿度分布的验证图 3 给出了工况 1 的试验与 CFD 计算比较，其中试验数据为所有测点的测量值均达到稳定状态时的结果。模拟值中，qm=486W/m2 为利用式 2）和 3）计算的实际水蒸气扩散散热量， （此时 Fm 约为 0.52）。为了对比，我们又假设了水蒸气蒸发时7都保持气相，没有发生相变化，扩散散热全部由显热交换构成的情况（qm=184W/m2）以及水蒸气蒸发时为全热交换的情况（q=1026W/m2）。由图可知，水蒸气扩散散热量对室内温度分布影响很大，如采用 qm=1026W/m2 顶棚附近的温
8、度比测定值高 2 度左右。某些研究11在计算浴室热湿负荷时，主张以全热交换来概算水蒸气扩散散热量势必造成很大的计算误差。相比之下，因为小室内湿度主要由水蒸气质量平衡决定，扩散散热值对小室内湿度的平均水平影响不大，但由于热流束的浮力效果不同，水面附近的气流方式导致湿度分布发生微妙的变化。综合地看，采用 q=486W/m2 的计算结果，无论温湿度，与实测都最为吻合。图 4 为测定断面上的流场,温湿度场的 CFD 模拟结果.由图可见,在此断面上从水面处形成的热湿羽流几乎没有受到小室内通风的影响,温湿度成层现象非常明显。图 3 实测与模拟的温湿度分布比较（工况 1）图 4 CFD 计算结果（左：气流场
9、；中：温度场；右：湿度场）4.2 结露形成与发展过程的验证图 5 给出了工况 2 的试验与 CFD 计算结果的比较。试验和模拟时间均为 30 分钟左右。为了更好地形成表面结露并防止出现小室内空气湿度达到 100%的情况，本试验工况进行时，进口处通过预加热装置对进风加热至 25.2 度。由图可见，实测和模拟都显示经过 20 分钟左右顶棚处的测点（P1-7 和 P2-7）湿度达到饱和，表明出现了结露。这说明虽然到目前为止结露问题还没有好的直接测定方法，通过比较小室内的温湿度动态分布，计算与模拟值随时间的变化规律基本一致，可以认为利用此计算模型来分析结露问题是可行的。8我们针对工况 2，利用湿结露
10、CFD 计算模型对结露进行了动态模拟，模拟时间为 1 小时。图 6 给出了 4 个时间点的结露分布。试验开始 20 分钟左右，结露首先在小室后上部角落出现，然后以较快的速度沿顶棚和侧壁发展，在 45 分钟以后结露面积基本达到稳定，但结露量继续增加。从结露量上看，小室后部两个侧壁和顶棚的结露量较多。这是因为前部侧壁和地面离进风口较近，热风使这些壁面温度升高的缘故。图 5 实测与模拟的温湿度动态变化比较（工况 2）图 6 CFD 模拟的结露分布随时间变化图 5.通风对结露的影响通风是解决结露问题的重要手段之一，但是到目前为止很少有这方面的量化分析。我们利用 CFD 湿结露计算模型通过 3 个算例针
11、对不同的通风量和通风方式对结露的影响进行了初步的探讨。算例 1：即工况 2；算例 2：通风量由 7.9 增至 9.4m3/h，其他条件不变；算例 3：通风方式改为上送下排的机械排风方式，其他条件不变。图 7 为算例 2 和 3 的结露分布模拟结果(t=60min)。与图 6 相比较，由于通风量增加，算例 2 的结露面积相应地减少，特别是侧壁靠下部的区域结露基本消失。算例 3 的通风量与算例 1 完全相同，但结露分布形状完全不同。因为变成上送下排的方式，抑制了水面附近的浮力效果造成的水蒸气上升现象，顶棚与侧壁上部的结露减少。图 8 给出了三个算例的各壁面及总的结露量的计算结果(t=60min)。
12、由9图可知，算例 2 和算例 3 的总结露量分别只有算例 1 的 26%和 20%。图 7 不同通风量与通风方式对结露分布的影响（左：算例 2；右：算例 3）图 8 不同通风量与通风方式下结露量的变化 6.结论为利用 CFD 技术研究室内的湿度分布与结露问题，本研究开发了湿结露 CFD 计算模型，通过对模型的验证和算例计算，可得以下结论：1） 无论是不考虑结露，只考虑湿度分布的稳态计算，还是考虑结露量的非稳态计算，CFD 模型的计算结果都与试验结果较为吻合。这说明可以利用此模型进行室内湿度场和结露的详细分析。特别是结露问题，在到目前为止还没有有效的测定方法的情况下，CFD 的应用提供了非常重要
13、的研究手段； 2） 本研究中还涉及到热湿同时传递问题，给出了新的计算方法并应用于 CFD 的计算。3） 即使是同样的风量，不同的通风方式对室内湿度分布和整体的湿度水平影响很大，在室内湿度较大的场合，可以通过通风方式的优化设计更有效地解决除湿和结露问题；7.参考文献1 Straube JF. Moisture in buildings. ASHRAE Journal 2002; (1):15-9.2 Sterling EM, Arundel A, Sterling TD. Criteria for human exposure to humidity in occupied buildings.
14、 ASHRAE Transactions 1985;(91):103 Murray AB, Ferguson AC, Morrison BJ. Ensitizaiton to house dust mites in a North American cit y. Journal of Allergy and Clinical Immunology 1985;76(1):108-12.4 Bates JM, Rorek DA, Ballantye MH. Dust mite counts and mite allergens in family homes before and after dr
15、y extraction carpet cleaning, Proceedings of the 6th International Conference on Air Quality and Climate, 1993, Vol.3, p.33-8.5 Harving H, Korsgaard J, Dahl J, Beck HI, Bjerring P. House dust mites and atopic dermatitis. Annals of Allergy 1990;65:25-31.6 Chao NT, Wang WA, Chiang CM. Study of control
16、 strategy using outdoor air to reduce winter indoor humidity in Taiwanese apartments demonstrated by ventilation design for a bathroom. ASHRAE Transactions. 1996; 102(1): 182-91.7 Kolokotroni M, Saiz N, Littler J. Moisture movement: A study using tracer gas techniques and CFD modeling. Building Services Engineering and Technology. 1992;13(2): 113-17.8 近藤靖史、長澤康弘、藤村淳一.湿度浮力影響考慮室内温熱環境 予測（）室内空気中水蒸気空間温度分布与影響.日本建築学会計画 系論文集. No.5
