温室墙体中覆铝箔封闭空气腔热工性能模拟分析

p第 33 卷 nbsp; nbsp;第 2 期 nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; 农 nbsp;业 nbsp;工 nbsp;程 nbsp;学 nbsp;报 nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp;Vol.33 nbsp;No.2 2017 年 nbsp; nbsp; nbsp;1月 nbsp; nbsp; nbsp; nbsp; nbsp; Transactions of the Chinese Society of Agricultural Engineering nbsp; nbsp; nbsp; nbsp; Jan. 2017 nbsp; nbsp; 227 nbsp; 温室墙体中覆铝箔封闭空气腔热工性能模拟分析张潇丹，颉建明※，郁继华，吕 nbsp; 剑 nbsp;（甘肃农业大学园艺学院，兰州 nbsp;730070） nbsp;摘 nbsp; 要 通过建立封闭空气腔二维稳态流动传热模型和温室墙体一维非稳态导热模型，模拟计算封闭腔内空气温度分布，研究了日光温室墙体中覆铝箔封闭空气腔的热工性能。结果表明壁面覆铝箔可有效减少封闭空气腔的辐射换热量；封闭空气腔的热阻随封闭腔高度的增加而增大，高度达 1.5 m 后，热阻趋于不变；封闭空气腔的厚度小于 0.03 m 时，其热阻随厚度增加而增大，厚度超过 0.03 m 后，热阻逐渐减小；覆铝箔封闭空气腔高度为 1.5 m、厚度为 0.03 m、内外壁面温差为 2～ 20 K 时，热阻为 0.70～ 0.55 Km2/W，保温隔热效果相当于 0.81～ 0.64 m 厚夯实黏土结构、 0.55～ 0.43 m 厚红砖砌体结构墙体或 0.20～ 0.16 m 厚煤渣、 0.06～ 0.05 m 厚珍珠岩、 0.03～ 0.02 m 厚聚苯板隔热材料。 3 组 30 mm 厚覆铝箔封闭空气腔加 480 mm 红砖复合墙体 （ 360 mm 红砖墙 3 组 30 mm 封闭空气腔 120 mm 红砖墙， 240 mm 红砖墙 3 组 30 mm封闭空气腔 240 mm 红砖墙） ，其夜间向室内放热量较单一 480 mm 红砖墙体提高 99.5～ 104.2，与相同结构聚苯板红砖复合墙体无明显差距。 nbsp;关键词 墙体；温度；对流换热；日光温室；封闭空气腔；平均导热系数；热阻 nbsp;doi 10.11975/j.issn.1002-6819.2017.02.031 中图分类号 S625.1 nbsp; nbsp; nbsp; nbsp; nbsp; 文献标志码 A nbsp; nbsp; nbsp; nbsp; nbsp; 文章编号 1002-68192017-02-0227-07 张潇丹，颉建明，郁继华，吕 nbsp;剑. 温室墙体中覆铝箔封闭空气腔热工性能模拟分析[J]. 农业工程学报，2017，332227－233. nbsp; nbsp;doi 10.11975/j.issn.1002-6819.2017.02.031 nbsp; nbsp;http//www.tcsae.org Zhang Xiaodan, Xie Jianming, Yu Jihua, L Jian. Simulation analysis of thermal properties of air enclosure covered with aluminum foil in wall of solar greenhouse[J]. Transactions of the Chinese Society of Agricultural Engineering Transactions of the CSAE, 2017, 332 227 － 233. in Chinese with English abstract nbsp; nbsp;doi 10.11975/j.issn.1002-6819.2017.02.031 nbsp; nbsp;http//www.tcsae.org 0 nbsp;引 nbsp;言日光温室由前屋面、后屋面、后墙、山墙构成，后墙作为日光温室的承重结构，在保温设计上也具有重要的意义[1-3]。从 20 世纪 30 年代发展至今，后墙结构和材料不断更新换代，从单一材料墙体到多层异质复合墙体，其蓄热、保温、隔热性能都有了大幅提升[4-7]。理想的后墙结构由蓄热层、保温层、隔热层 3 部分构成。蓄热层日间储存热量，夜间放出热量，是后墙内部温度随时间波动较大的部分；保温层为热稳定层，温度几乎没有波动，除了阻止热量向外流出，还起了缓冲作用，保证了蓄放热的均匀性； 隔热层一般采用导热系数小于 0.05 W/mK的材料，用于截断经由保温层流出的热量[8-15]。由于静止的空气导热系数很小，是良好的保温隔热材料，近年来封闭空气腔结构在各领域的应用日趋广泛。在日光温室的改造升级中，也出现了由内外层砖墙和中间空气腔构成的中空墙体[16]。亢树华等[17]通过对不同材料和结构的日光温室后墙进行温度观测，得出隔热性能最好的是珍珠岩，其次是煤渣，再次是锯末，最差是中空。白义奎收稿日期 2016-08-21 nbsp; nbsp;修订日期 2016-10-14 基金项目国家公益性行业（农业）专项资助（ 201203002） nbsp;作者简介张潇丹，女，从事设施园艺研究。兰州 nbsp; 甘肃农业大学园艺学院，730070。 Email ※通信作者颉建明，男，教授，博导，主要从事设施农业教学与研究。兰州 nbsp; 甘肃农业大学园艺学院， 730070。 Email 等[18]基于前人对封闭腔空气自然对流传热的研究，对缀铝箔聚苯板空心墙体保温性能进行了理论分析，结果表明空气层厚度在 0.02～ 0.10 m 间的缀铝箔聚苯板空心墙体具有较好的保温隔热效果。李凯等[19]、凌浩恕等[20]研究了带有竖向空气层的相变蓄热墙体的热工性能。 在已有的研究中， 一般是通过试验测量得到数据再归纳分析，或凭借经验公式直接进行计算分析，但是，通过试验测量得到的数据有限， 已有的经验公式也往往不能适用于日光温室墙体热环境， 要探索封闭空气腔在日光温室墙体保温设计中的合理结构， 这些方法具有一定的局限性。 nbsp;封闭腔的空气自然对流换热是计算流体力学和计算传热学中的一个经典课题[21]。 Davis[22]最早给出了正方形封闭腔中，无量纲数 Ra 在 103～ 106内时空气自然对流换热的 Nue 基准解。 Catton[23]给出了封闭腔高厚比在 1～ 10内的努赛尔数关联式。李光正等[24-25]、李世武等[26]、董韶峰等[27]、黄建春等[28]运用多种数值解法模拟计算了封闭腔内的空气对流传热问题。 这些研究均是将封闭空气腔作为一个独立系统进行数值计算，并没有结合具体情况进行研究。本文根据日光温室墙体的建筑结构、内部温度变化及热量传递规律，建立封闭腔二维稳态流动传热模型和墙体一维非稳态导热模型， 采用 MATLAB 软件编制程序，模拟计算封闭腔内空气温度分布，分析温室墙体中不同模式封闭空气腔的热工特性，以期为封闭空气腔在日光温室墙体中的应用提供理论依据。 nbsp;农业生物环境与能源工程农业工程学报（ http//www.tcsae.org） nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; 2017 年 nbsp; nbsp;228 1 nbsp;封闭空气腔热工计算 封闭空气腔内存在多种传热形式，其中热传导占10， 对流换热占 20～ 30， 辐射换热占 60～ 70[29]。根据传热学理论，封闭空气腔的表面对流换热系数 hc由下式计算 nbsp;/ceAhNu Lλ⋅ nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; （ 1） nbsp;式中 λA为空气的导热系数，在 273 K（ 0 ℃）时， λA 0.025 W/mK； L 为封闭腔的厚度， m；努赛尔数 Nu为表征流体中对流换热与导热传热之比的无量纲参数，数值上等于封闭腔表面处的无量纲温度梯度， 封闭空气腔内外壁面高度方向上各点的 Nu 的均值为平均努赛尔数 Nue。 nbsp;封闭空气腔内辐射换热发生在两个面积相等、存在温差的壁面之间，辐射换热系数 hr由下式计算 nbsp;14 412 12121 / 1 / 1 / rbhTTTσεε−⋅ − ⋅ − − nbsp; （ 2） nbsp;式中 σb为 Stefan-Boltzmann 常数，取值 5.6710-8； ε1和 ε2分别为 2 个壁面在一定温度下的发射率 ； T1和 T2为 2 个壁面的热力学温度， K。辐射换热占封闭空气腔内传热的很大比例，严重影响其保温隔热性能。两表面间的辐射换热量和两表面温差、 角系数及表面发射率相关。室温条件下，红砖、混凝土等建筑材料的表面发射率均在 0.9 以上，铝箔的表面发射率约为 0.2，若用其将封闭腔沿温差方向隔开， 封闭腔内的辐射换热量将大幅降低[30]。 nbsp;热阻 R 作为封闭空气腔保温隔热性能的评价指标，由下式计算 nbsp;1/ crRhh nbsp; nbsp; nbsp; nbsp; （ 3） nbsp;式中 hchr为封闭空气腔的总传热系数， W/m2K。封闭空气腔内对流换热系数 hc与辐射换热系数 hr的当量导热系数 λe由下式计算 nbsp;ecrhhLλ ⋅ nbsp; nbsp; nbsp; nbsp; nbsp;（ 4） 2 nbsp;封闭空气腔内部流动传热数值模拟计算 nbsp;2.1 nbsp;物理模型 1）封闭空气腔结构设计 如图 1a 所示，根据日光温室墙体建筑结构，封闭空气腔的高度 H 分别设为 0.3、 0.6、 0.9、 1.2、 1.5、 1.8、2.1、 2.4 m， 厚度 L 分别设为 0.01、 0.02、 0.03、 0.04、 0.05、0.08、 0.10、 0.12、 0.24 m。 2）封闭空气腔温度设计 nbsp; nbsp;假设封闭空气腔的上下壁面均为绝热壁面，参考多层异质复合墙体中温度的变化和冬季室外平均温度，在温室墙体保温层模拟中，墙体温度偏高且沿厚度方向的变化较小，将封闭腔内外壁面的温差 ΔT1分别设为 2、 4、6、 8、 10 K，内壁面温度 T1依次取值 275、 277、 279、281、 283 K，外壁面温度 T2取值 273 K；在墙体隔热层模拟中，墙体温度偏低且沿厚度方向的变化较大，将内外壁面的温差 ΔT2分别设为 10、 20 K， T1依次取值 273、283 K， T2取值 263 K。 nbsp;注 T1、 T2为封闭空气腔内外壁面温度， L、 H 为封闭空气腔厚度和高度；Δx、 Δy 为水平方向和垂直方向的空间步长， W、 E、 N、 S、 P 为 5 个相邻的控制容积节点。 nbsp;Note T1, T2are temperatures of inside and outside surface of the air enclosure; L, H are thickness and height of the air enclosure; Δx, Δy are space steps in horizontal direction and vertical direction; W, E, N, S, P are five adjacent control volume nodes. 图 1 nbsp;封闭空气腔的物理模型及网格划分图 nbsp;Fig.1 nbsp;Physical model and meshing figure of air enclosure 2.2 nbsp;数学模型及计算模拟方法 建立封闭空气腔二维稳态流动传热模型，模拟计算腔内温度分布，获得其内外壁面处无量纲温度梯度（ Nue），得到封闭空气腔内对流换热系数 hc。以下是描述二维不可压缩流体稳态流动传热的 控制方程 nbsp;流函数方程 nbsp;2222xyϕϕω∂∂∂∂（ 5） nbsp;涡量方程 nbsp;uv Tgx yxxyy xρω ρω ω ωηηρβ∂ ∂∂∂∂∂ ∂ −∂ ∂∂∂∂∂ ∂（ 6） nbsp;能量方程 nbsp;ppuT vT T Tx yxCxyCyρρ λ λ∂ ∂∂∂∂∂ ∂ ∂∂∂∂∂（ 7） nbsp;式中 φ、 ω、 T 为所求变量，分别代表流函数（ m2/s）、涡量（ s-1）和温度（ K）； ρ 为气体密度， kg/m3； η 为动力黏度， Ns/m2； g 为重力加速度， m/s2； β 为气体膨胀系数， K-1； λ 为气体导热系数， W/mK； Cp为气体定压比热容， J/kgK； u、 v 分别为气体流速在 x、 y 方向上的分量， m/s。 nbsp;采用有限容积法，把封闭腔沿 x、 y 方向划分为均匀网格，生成网格后，将描写流动传热的偏微分方程离散为各节点上的代数方程。如图 1b，阴影部分为控制容积，x 方向的空间步长 Δx 为 0.005 m， y 方向的空间步长 Δy最大为 0.015 m。以能量方程为例稳态流动中，控制容积内热力学能的增加率为 0， 通过空气对流进入控制容积的净热量等于该控制容积内空气向四周扩散的净热量。对于内节点 P 所代表的控制容积，与其热平衡密切相关的控制容积分别为 W、 E、 N、 S。对流作用带有强烈的方向性，采用一阶迎风格式离散对流项；扩散作用具有均匀性，采用中心差分格式离散扩散项，流函数、涡量、能量方程将离散为 3 组代数方程。 nbsp;第 2 期 nbsp;张潇丹等温室墙体中覆铝箔封闭空气腔热工性能模拟分析 nbsp;229 采用迭代法联解方程组，设置两部分迭代外迭代和内迭代。在一个外迭代层次上，按温度涡量流函数的内迭代路线运算，内迭代采用 Gauss-Seidel 迭代法，以前后两次待求量值的差值小于允许值作为收敛判据；内迭代收敛完成后，再进入下一层外迭代计算，外迭代以封闭腔内外壁面处无量纲温度梯度（ Nue）在前后两层外迭代中的偏差小于允许值作为收敛判据。编写MATLAB 计算程序， 涡量迭代收敛的允许值未超过 10-3，温度、流函数收敛的允许值未超过 10-4，外迭代收敛的允许值为 10-5，内外迭代均在 5 000 步以内完成收敛。 nbsp;边界条件为能量方程中，上下壁面为绝热壁面，热流量恒为 0，内外壁面温度为确定值；涡量方程中，采用 Thom 公式[21]持续调整壁面涡量值；流函数方程中，在壁面上对速度取无滑移边界条件，流函数恒为 0。 nbsp;2.3 nbsp;自然对流验证 封闭空气腔内外壁面温度为确定值，腔内空气自然对流与壁面热辐射无耦合关系，可单独验证。当封闭空气腔厚度为 0.10～ 0.12 m，高度为 0.6～ 1.2 m，内外壁面温差为 10 K（保温层）时，封闭空气腔物理模型与 Catton[23]所研究的物理模型基本一致，表 1 为利用文献经验公式计算得到的 Nue结果与本文模拟计算得到的 Nue结果， 经对比，本文计算结果是正确的，说明通过模拟计算研究封闭空气腔热工性能是可行的。 nbsp;表 1 nbsp;封闭空气腔的平均努塞尔数 Nue计算结果 nbsp;Table 1 nbsp;Average Nusselt number of air enclosure 封闭空气腔高度 H Height of air enclosure /m 封闭空气腔厚度 L Thickness of air enclosure /m 方法 nbsp; 0.6 0.9 1.2 经验公式法[23]6.192 5.959 5.207 0.10 数值模拟法（本文） nbsp;6.115 5.766 5.497 经验公式法[23]8.046 7.271 6.766 0.12 数值模拟法（本文） nbsp;7.572 7.144 6.850 注 L 为封闭空气腔厚度， Ra 为瑞利数，封闭腔高度为 0.3 m。 Note L is thickness of air enclosure. Ra is Raleigh number, height of air enclosure is 0.3 m. 图 2 nbsp;不同 Ra 下封闭空气腔内等温线分布 nbsp;Fig.2 nbsp;Isotherm distributions in air enclosure under different Raleigh number Ra 为表示流体流动状态发生过渡的条件的无量纲数。封闭空气腔高度为 0.3 m，厚度为 0.02～ 0.10 m，内外壁面温差为 10 K（保温层）时， Ra 在 103～ 106间，且随厚度的增大而增大。该条件下封闭腔内等温线的分布见图 2，温度从高温壁向低温壁逐渐递减， Ra 较小时，等温线为一簇几乎平行的竖直线，封闭腔内空气几乎没有流动；随着 Ra 的增大，等温线逐渐发生波动，高温壁面等温线向上凸起，低温壁面等温线向下凸起，出现了明显的温度边界层，腔内自然对流开始发展；当 Ra 达到106后，封闭腔中部等温线几乎为水平线，腔内空气自然对流加剧，开始由层流向湍流发展。模拟结果呈现的封闭腔内空气流动规律与理论描述[31]基本一致，说明了本研究所采用的数值模拟方法的可靠性。 nbsp;3 nbsp;结果与分析 3.1 nbsp;影响封闭空气腔保温隔热性能的因素 3.1.1 nbsp;封闭空气腔的不同壁面处理及高度对热阻的影响 封闭空气腔厚度为 0.03 m，外壁面温度为 273 K，内外壁面温差为 5 和 10 K，壁面覆铝箔（双面，下同）和未覆铝箔时，不同高度对封闭空气腔热阻的影响见图 3。 nbsp;注封闭空气腔厚度为 0.03 m，外壁面温度为 273 K，下同。 nbsp; Note Thickness of air enclosure is 0.03 m, temperature of outside surface is 273 K, the same as below. 图 3 nbsp;不同壁面处理下封闭空气腔随其高度、温差的变化 nbsp;Fig.3 nbsp;Thermal resistance under different surfaces, temperature difference, height of air enclosure 由图 3 可知，壁面覆铝箔的封闭空气腔，厚度为0.03 m、高度为 0.3～ 2.4 m、内外壁面温差为 5 K 时，与壁面未覆铝箔的封闭空气腔相比，热阻增加约 210.8；内外壁面温差为 10 K 时，热阻增加约 196.9。封闭空气腔的内外壁面覆铝箔，减小了壁面发射率，降低了腔内的辐射换热量，提高了封闭空气腔的保温隔热性能。 nbsp;高度增加，封闭空气腔的热阻增大，且趋势逐渐变缓，达到一定高度后，热阻趋于不变。图 3a 中，封闭空气腔壁面未覆铝箔，当封闭空气腔厚度为 0.03 m，内外壁面温差为 5 和 10 K 时，高度从 0.3 m 增加到 1.5 m，热阻分别增加 6.8、 7.1；高度从 1.5 m 增加到 2.4 m，热阻均增加 1.5。图 3b 中，封闭空气腔壁面覆铝箔，内外壁面温差为 5 和 10 K 时，高度从 0.3 m 增加到 1.5 m，热阻分别增加 22.0、 20.6；高度从 1.5 m 增加到 2.4 m，热阻分别增加 4.2、 5.1。考虑到日光温室墙体建筑结构和实际施工情况[1]， 温室墙体保温设计中封闭空气腔的适宜高度为 1.5 m。 nbsp;3.1.2 nbsp;封闭空气腔的不同壁面处理及厚度对热阻的影响 封闭空气腔高度为 1.5 m，外壁面温度为 273 K，内外壁面温差为 5 和 10 K，壁面覆铝箔和未覆铝箔时，不农业工程学报（ http//www.tcsae.org） nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; 2017 年 nbsp; nbsp;230 同厚度对封闭空气腔热阻的影响见图 4。 nbsp;由图 4 可知，壁面覆铝箔的封闭空气腔，高度为1.5 m、厚度为 0.03～ 0.12 m、内外壁面温差为 5 K 时，与壁面未覆铝箔的封闭空气腔相比， 热阻增加约 209.9；内外壁面温差为 10 K 时，热阻增加约 196.4，再次表明封闭空气腔壁面覆铝箔可大幅提升其保温隔热性能。 nbsp;注封闭空气腔高度为 1.5 m。 nbsp;Note Height of air enclosure is 1.5 m. 图 4 不同壁面处理下封闭空气腔热阻随厚度、温差的变化 nbsp;Fig.4 nbsp;Thermal resistance under different surfaces, temperature difference, thickness of air enclosure; nbsp;图 4a、 4b 中，随着厚度增大，封闭空气腔热阻的变化趋势基本一致。封闭空气腔高度为 1.5 m 时，其热阻在封闭腔厚度小于 0.03 m 时随厚度的增加而增大，在厚度为 0.03 m 时存在最大值；厚度在 0.03～ 0.12 m 间，内外壁面温差为 5 K 时，热阻随厚度增加而减小，内外壁面温差为 10 K 时，厚度对热阻无明显影响；当厚度超过0.12 m 后，热阻逐渐减小。封闭空气腔在 0.03～ 0.12 m厚度范围内时，其热阻较大、保温隔热性能较优。 3.2 nbsp;覆铝箔封闭空气腔在温室墙体保温设计中的热工性能分析 日光温室墙体保温设计中，不同厚度、不同内外壁面温差下，高度为 1.5 m 的覆铝箔封闭空气腔在保温层、隔热层模拟中的热阻见图 5，平均导热系数（同一厚度封闭空气腔在不同内外壁面温差下的当量导热系数的均值）见表 2。 nbsp;注 ΔT1、 ΔT2分别为外壁面温度 T2为 273、 263 K 时封闭空气腔内外壁面的温差。封闭空气腔高度为 1.5 m。 nbsp;Note ΔT1, ΔT2 are temperature differences between inside and outside surface of the air enclosure when T2 is 273K, 263 K, respectively. Height of air enclosure is 1.5 m. 图 5 nbsp;不同厚度、温差下覆铝箔封闭空气腔的热阻 nbsp;Fig.5 nbsp;Thermal resistance under different temperature difference and thickness of air enclosure covered with aluminum foil 表 2 nbsp;不同厚度覆铝箔封闭空气腔的平均导热系数 Table 2 nbsp;Average thermal conductivity under different thickness of air enclosure covered with aluminum foil 平均导热系数 nbsp;Average thermal conductivity/WmK-1 封闭空气腔的厚度 nbsp;Thickness of air enclosure/m 封闭腔在保温层 nbsp;Air enclosure in heat preservation layer 封闭腔在隔热层 nbsp;Air enclosure in thermal insulation layer 0.02 0.036 0.037 0.03 0.047 0.052 0.04 0.063 0.071 0.05 0.080 0.090 0.08 0.130 0.145 0.10 0.165 0.181 0.12 0.198 0.217 注 平均导热系数指同一厚度封闭空气腔在不同内外壁面温差下的当量导热系数的均值，保温层模拟中的内外壁面温差为 2～ 10 K，隔热层模拟中的内外壁面温差为 10～ 20 K。 nbsp;Note Average thermal conductivity is the mean of equivalent thermal conductivity of the same thickness air enclosure in different temperature difference between inner and outer surface; In heat preservation layer and thermal insulation layer, temperature difference is 2-10 K, 10-20 K, respectively. 3.2.1 nbsp;覆铝箔封闭空气腔在温室墙体保温层中的热工参数 在日光温室墙体保温层模拟中， 封闭空气腔高度外壁面温度为 273 K，内外壁面温差为 210 K。由图 5 可知，高度为 1.5 m 的覆铝箔封闭空气腔，内外壁面温差为 2 K时，热阻在封闭空气腔厚度为 0.04 m 时存在最大值；内外壁面温差为 4～ 10 K 时，热阻在封闭空气腔厚度为0.03 m 时存在最大值，此时封闭空气腔的保温隔热性能最优。表 2 中，封闭空气腔厚度为 0.020.12 m 时，其平均导热系数随厚度增加而增大、为 0.0360.198 W/m⋅K。厚度为 0.03 m 的封闭空气腔， 内外壁面温差为 2～ 10 K 时，其平均导热系数为 0.047 W/mK ，热阻为 0.70 ～0.58 Km2/W（图 5），保温隔热效果相当于 0.81～ 0.67 m厚夯实黏土结构、 0.55～ 0.45 nbsp;m 厚红砖砌体结构墙体或0.20～ 0.17 m 厚煤渣、 0.06～ 0.05 m 厚珍珠岩、 0.03～0.025 m 厚聚苯板隔热材料[31]。 3.2.2 nbsp;覆铝箔封闭空气腔在温室墙体隔热层中的热工参数 在日光温室墙体隔热层模拟中，封闭空气腔的外壁面温度为 263 K，内外壁面温差为 10～ 20 K。由图 5 可知，日光温室墙体中封闭空气腔高度为 1.5 m，内外壁面温差为 10、 20 K 时，热阻在封闭空气腔厚度为 0.03 m 时存在最大值，保温隔热效果最佳。表 2 中，封闭空气腔厚度为 0.02～ 0.12 m 时，其平均导热系数随厚度增加而增大、为 0.037～ 0.217 W/mK。 厚度为 0.03m 的覆铝箔封闭空气腔，内外壁面温差为 10～ 20 K 时，其平均导热系数为0.052 W/mK，热阻为 0.60～ 0.55 Km2/W，保温隔热效果相当于 0.70～ 0.64 m 厚夯实黏土结构、 0.47～ 0.43 m 厚红砖砌体结构墙体或 0.17～ 0.16 m 厚煤渣、 0.06～ 0.05 m厚珍珠岩、 0.025～ 0.02 m 厚聚苯板隔热材料[31]。 3.3 nbsp;带有覆铝箔封闭空气腔的温室墙体的传热分析 以甘肃酒泉典型晴天（ 2014 年 12 月 7 日）下日光温室内外气温及太阳辐射量作为边界条件（表 3），选择480 mm 红砖墙作为计算模型， 将组 30 mm 厚覆铝箔封闭第 2 期 nbsp;张潇丹等温室墙体中覆铝箔封闭空气腔热工性能模拟分析 nbsp;231 空气腔组合起来作为保温隔热材料，放置在红砖墙厚度方向上的不同位置，形成 5 个墙体结构方案（表 4）。建立一维非稳态传热模型，编写 MATLAB 程序，计算温室墙体的内表面日间（ 0900－ 1700）总蓄热量和夜间（ 1700−次日 0900）总放热量。根据李明等[13]的方法确定了在表 3 温室内外环境因素扰动下 480 mm红砖墙的蓄热层约为 160 mm，温度随厚度波动较小的保温层约为200 mm，墙体靠近室外一侧、温度波动较大的厚度约为120 mm。依据表 2 结果，当封闭空气腔在红砖墙内部偏外侧（方案 2）及正中间（方案 4）时，其平均导热系数均按照保温层模拟取值，为 0.047 W/mK。 nbsp;表 3 nbsp;酒泉地区 2014 年 12 月 7 日日光温室内、外气温和墙体内表面太阳辐射量的逐时数据 Table 3 nbsp;Hourly temperature inside and outside greenhouse and solar radiation in JiuQuan Dec 7, 2014 时刻 nbsp;Time 室外气温 nbsp;Outdoor temperature/℃ nbsp;室内气温 nbsp;Indoor temperature/℃ nbsp;墙体内表面太阳辐射量Solar radiation on inner surface of wall/Wm-20000 -9.3 9.6 0 0100 -9.3 9 0 0200 -9.8 8.7 0 0300 -12.1 8.2 0 0400 -11 7.7 0 0500 -11 7.3 0 0600 -13.3 7.1 0 0700 -14 6.7 0 0800 -13.3 6.5 0 0900 -12.1 6.3 0 1000 -7.1 8 125 1100 -4.6 14.2 302 1200 -2.2 24 382 1300 -0.8 20.1 473 1400 1 20.7 434 1500 1.9 19.8 409 1600 1.9 19.4 285 1700 1 14.5 53 1800 -1.7 15.6 0 1900 -3.1 13.9 0 2000 -4.6 12.5 0 2100 -7.1 11.5 0 2200 -8.3 10.8 0 2300 -8.8 10.3 0 表 4 nbsp;日光温室墙体内表面全天蓄放热量 Table 4 nbsp;Heat accumulate and release of inner surface of greenhouse wall nbsp;方案 nbsp;编号 nbsp;Scheme number 墙体结构 nbsp;Wall structure 日间总蓄热量Heat accumulate at daytime/ MJm-2夜间总放热量Heat release at nighttime/ MJm-21 480 mm 红砖墙 nbsp;5.49 2.14 2 360 mm 红砖墙 30 mm 30 mm30 mm封闭空气腔 120 mm 红砖墙 nbsp;4.75 4.27 3 360 mm 红砖墙 90 mm 聚苯板120 mm 红砖墙 nbsp;4.71 4.30 4 240 mm 红砖墙 30 mm 30 mm30 mm封闭空气腔 240 mm 红砖墙 nbsp;4.69 4.37 5 240 mm 红砖墙 90 mm 聚苯板240 mm 红砖墙 nbsp;4.67 4.38 注聚苯板的密度为 18 gm-3。 nbsp;Note Destiny of EPS is 18 gm-3. 晴天条件下，方案 2、 4 较方案 1 的墙体内表面日间总蓄热量分别降低 13.5、 14.6，夜间总放热量分别提高 99.5、 104.2，覆铝箔封闭空气腔覆铝箔封闭空气腔大幅提高了温室墙体保温性能；方案 2 与方案 3、方案 4与方案 5 的墙体日间总蓄热量及夜间向温室内的总放热量无明显差距，说明覆铝箔封闭空气腔可以替代聚苯板在日光温室墙体建造中应用。 nbsp;4 nbsp;结 nbsp;论 本文基于传热学理论，根据日光温室墙体的建筑结构、内部温度变化及热量传递规律建立了二维稳态流动传热模型，模拟分析了日光温室墙体中不同模式封闭空气腔的热工性能，得出以下结论 nbsp;1）中空墙体壁面覆铝箔可有效减少辐射换热量，大幅提高其保温隔热性能；封闭空气腔的热阻随其高度的增加而减小，高度达 1.5 m 后，热阻趋于不变；封闭空气腔的厚度小于 0.03 m 时，其热阻随厚度增加而增大，厚度超过 0.03 m 后，热阻逐渐减小。 nbsp;2）覆铝箔封闭空气腔的适宜高度为 1.5 m，厚度为0.03 m。当封闭空气腔内外壁面温差为 2～ 20 K 时，该结构的热阻为 0.70～ 0.55 Km2/W，保温隔热效果相当于0.81～ 0.64 m 厚夯实黏土结构、 0.55～ 0.43 m 厚红砖砌体结构墙体或 0.20～ 0.16 m厚煤渣、 0.06～ 0.05 m厚珍珠岩、0.03～ 0.02 m 厚聚苯板隔热材料。 nbsp;3） 3 组 30 mm 厚覆铝箔封闭空气腔加 480 mm 红砖复合墙体 （ 360 mm 红砖墙 3 组 30 mm 封闭空气腔120 mm 红砖墙， 240 mm 红砖墙 3 组 30 mm 封闭空气腔 240 mm红砖墙） ， 其夜间向室内放热量较单一 480 mm红砖墙体提高 99.5～ 104.2， 与相同结构聚苯板红砖复合墙体无明显差距。覆铝箔封闭空气腔可以替代聚苯板在日光温室墙体建造中应用。 nbsp;[参 nbsp;考 nbsp;文 nbsp;献] [1] 邹志荣，周长吉 . 温室建筑与结构 [M]中国农业出版社，2012. nbsp;[2] 刘志杰，郑文刚，胡清华，等 . 中国日光温室结构优化研究现状及发展趋势 [J]. 中国农学通报， 2007， 232 449－ 453. nbsp;Liu Zhijie, Zheng Wengang, Hu Qinghua, et al. Current situation and development on structure optimization of solar greenhouse in China[J]. Chinese Agricultural Science Bulletin, 2007, 232 449－ 453. in Chinese with English abstract [3] 陈青云 . 日光温室的实践与理论 [J]. 上海交通大学学报农业科学版， 2008， 265 343－ 350. nbsp;Chen Qingyun. Progress of practice and theory in sunlight greenhouse[J]. Journal of Shanghai Jiaotong University Agriculture Sciences, 2008, 265 343－ 350. in Chinese with English abstract [4] 魏晓明，周长吉，曹楠，等 . 中国日光温室结构及性能的演变 [J]. 江苏农业学报， 2012， 284 855－ 860. Wei Xiaoming, Zhou Changji, Cao Nan, et al. Evolution of structure and perance of Chinese solar greenhouse[J]. Jiangsu Journal of Agriculture Sciences, 2012, 284 855－860. in Chinese with English abstract 农业工程学报（ http//www.tcsae.org） nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; nbsp; 2017 年 nbsp; nbsp;232 [5] 柴立龙，马承伟，籍秀红，等 . 北京地区日光温室节能材料使用现状及性能分析 [J]. 农机化研究， 20078 17－ 21. Chai Lilong, Ma Chenwei, Ji Xiuhong, et al. Present situation and the perance analysis of the energy-saving material in solar greenhouse of Beijing region[J]. Journal of Agricultural Mechanization Research, 20078 17－ 21. in Chinese with English abstract [6] 佟国红，白义奎，赵荣飞，等 . 日光温室复合墙与土墙热性能对比分析 [J]. 沈阳农业大学学报， 2011， 426 718－722. nbsp;Tong Guohong, Bai Yikui, Zhao Rongfei, et al. Thermal property investigations for layered walls and earth walls in Chindes solar greenhouse[J]. Journal of shenyang Agricultural University, 2011, 426/p