基于CFD的玻璃温室环境数值模拟及优化分析.pdf

第 44 卷第 1 期                                   应用科技                              Vol.44 No.1 2019 年 1 月                                Applied Science and Technology                                 Jan. 2019 DOI: 10.11991/yykj. 201901005 优先出版  基于 CFD 的玻璃温室环境数值模拟及优化分析  贾鹤鸣，韩骏骋，张森，孙康健，李瑶  东北林业大学  机电工程学院，黑龙江  哈尔滨  150040 摘   要： 为进一步研究和优化玻璃温室内部气流场及温度场的分布，基于计算流体力学 (computational fluid dynamics, CFD)对东北地区玻璃温室内部气流 场进行三维建模数值分析， 从而分析得出内部速度场、 温度场的分布模式，并针对温室内部气流分布情况提出 3 种结构优化策略。模拟结果表明：监测点处的气流速度及温度的平均相对误差分别为5.54%、 4%，数值模拟结果与试验测试结果吻合度较高。 3 种优化方案结果均优于初始方案，其中通过分析得出方案三将暖风扇置于离地 0.7 m 处的为最佳方案，植物区域处于最适风速的部分增加了 50%，气流均匀性最好且平均温度为23.49 ，有效改善温室内环境并使其更加适宜植物的培养。  关键词： CFD；玻璃温室；气流运动；温度场；速度场；多孔介质；数值模拟；均匀性；结构优化  中图分类号： TP273        文献标志码： A        文章编号 ：  Numerical analysis simulation and optimization of the environment in glass greenhouse based on CFD JIA Heming， HAN Juncheng， ZHANG Sen， Sun Kangjian， Li YaoCollege of Mechanical and Electrical Engineering, Northeast Forestry University, Harbin, Heilongjiang Province, 150040, China Abstract: In order to further study and optimize the distribution of air flow and temperature field in glass greenhouse, a three-dimensional numerical analysis of air flow field in glass greenhouse in Northeast China based on Computational Fluid Dynamics (CFD) is carried out in this paper. The distribution patterns of internal velocity field and temperature field are analyzed, and three structural optimization strategies are put forward according to the distribution of air flow in glass greenhouse. The simulation results show that the average relative errors of gas velocity and temperature at the monitoring points are 5.54% and 4%, respectively. The numerical simulation results are in good agreement with the experimental results. The results of the three optimization schemes are better than the initial ones. The best scheme is that the warm fan is placed at 0.7m away from the ground. The optimum wind speed is increased by 50% in the plant area. The air flow uniformity is the best and the average temperature is 23.49 C. It can effectively improve the greenhouse environment and make it more suitable for plant cultivation. Keywords: CFD; glass greenhouse; air movement; temperature field; velocity field; porous medium; numerical simulation; uniformity; structure optimization 1温室植物的产量取决于温室内环境的温度、空气湿度等小气候要素的适宜度1，尤其是较好的温度和气流分布，对于植物生长有着极其重要的作用。在所有环境因子中，适宜的温度是植物生长的首要条件，过低的温度甚至会造成温室内植物的死亡。因此，温室通风口的设置和风速都会对温室内气流的分布情况有较大的影响。程秀花等1提出室内气流的速度和分布情况较易受到植物的影响；汪郑邦2认为植物温室收稿日期： 2019-01-07.          基金项目： 国家自然科学基金项目 (51609048). 作者简介： 贾鹤鸣 , 男，副教授，博士；  通信作者： 贾鹤鸣， E-mail： jiaheminglucky99126.com 的气流分布对作物生长有较大的影响，不均匀的气流分布容易在作物周围形成换流，影响作物的蒸腾和光合作用； 刘文合等3的研究表明通过对温室进行升温，可以使得温室内气温显著升高，但对温室内土壤深处的温度影响不大；周伟等4通过试验证明了通风口的设置和风速都会对温室内气流的分布情况有较大的影响，同时温室的气流分布对植物蒸腾和光合作用有较大的影响。因此，如何合理优化温室结构，得到适宜的风速和温度并使两者均匀性提高是当前急需解决的重要问题。  网络出版时间：2019-03-25 14:24:21网络出版地址：http:/kns.cnki.net/kcms/detail/23.1191.U.20190325.1424.006.html 2                                          应用科技                                 第 44 卷  1  材料和方法  1.1  温室模型   试验在东北林业大学进行， 模拟温室建立在室内，受外界环境影响因素较小。温室的长宽高为 2 m 0.83 m 1.8 m,由于温室是密闭的，因此忽略温室的门，并对温室四周及顶部进行简化即将温室的顶部和四周均视为以浮法玻璃为材料。考虑到植物的冠层阻力，将植物模型简化为长 1.50 m、宽 0.6 m、高 0.5 m 的多孔介质， 置于温室正中央。 暖气扇规格为 0.4 m 0.4 m，置于温室西侧，离地面 1.3 m 高处。采用Gambit 进行以温室西北角为零点的温室模型的建立和网格化，得到的玻璃温室模型如图1 所示。 图 1  玻璃温室模型 1.2 监测点布置  在温室中选取 Z=0.415 m 处的截面，均匀分布 6个监测点，高度分别为 Y=0.6 m 和 1.2 m，将监测点依次记为16PP。在平面 X=1 m 处，共建立 4 个监测点，高度分别为 0.6 m 和 1.2 m，记为710P P 。通过将在监测点处得到的实际值与模拟值进行对比验证。采用风速仪和温度计对温室内风速进行测量，模拟数值室内温度由使用 Fluent 软件进行模拟，监测点的布置如图 2、 3 所示。 图 2  Z=0.415 m 截面监测点布置 图 3  Z=0.415 m 截面监测点布置 1.3 数学模型  1.3.1  CFD 模型 当 Fluent模拟温室环境时，使用标准 k-模型5将温室内的气流视为遵循质量守恒方程、动量守恒方程和能量守恒方程的湍流。 动量守恒方程指出，所有的流动在实际情况下都应该遵循质量守恒定律，定律用公式表示为： 0uvwxyz 式中的 u、 v、 w 分别为在 x、 y、 z 方向上温室内气流流速的分量  动量守恒用公式表现为：  222eff 222()()()uu uv uw Px yzxuuuxyz 222refeff ref222()()()-1( )vu vv vw PxyzyvvvgTTxyz 222eff 222()()()wu wv ww Px yzzwwwxyz 式中：eff代表有效粘度；代表空气热膨胀系数；g代表重力加速度； T 代表空气温度；refT代表空气的参考温度。  能量守恒用公式表示为  222eff222()()( )PpuT vT wT qxyzCTTTCxyz 式中：q代表热源；pC代表比热容；eff代表有效导热系数6 。  在处理封闭室内因温差引起的浮升力项时采用Boussinesq 假设7。 Boussinesq 假设是对流体流动现象的一种简化， 假设认为： 流体中的粘性耗散忽略不计；除密度外其他物性为常数；对密度仅考虑动量方程中与体积力有关的项，其余各项中的密度亦作为常数。重力项中的密度如若采用作参考温度，则在不同温度T 下温室内部的空气密度 可以用数学方程表示为：   00=1-TT式中：是温室内气流的热膨胀系数； T 为温室内的温度；0T为温室内部空气的参考温度。   流经多孔介质的空气可通过 Darcy-Forchheimer理论描述为  2=-FppCSuKK式中：S代表动量源项；pK代表多孔介质的渗透性第 1期                  贾鹤鸣， 韩骏骋等： 基于 CFD的玻璃温室环境数值模拟及优化分析                             3  系数；FC代表非线性动量损失系数；代表空气的动力粘度；代表空气密度； u 代表空气流速9。  将植物视为各向同性的多孔介质8，并建立其与温室内部风速的数学关系模型：  2=-LAV DSICu式中：LAVI 代表叶面积指数；DC 是植物冠层的阻力系数。 对于多孔介质动量源项，使用幂律模型，自定义经验系数0C 为0.395，1C 为0.2。植物通过呼吸作用产生的热量与温室中气流进行热交换，将植物通过呼吸作用产生的热量视为常值，在进行数值模拟时将植物区域看作“体积热源”。  1.3.2  计算域与网格划分 在使用 Gambit 划分网格时将整个玻璃温室作为计算域， 采用非结构化四面体网格 (混合 TGrid 网格类型以及 Tet/Hybrid混合单元 )将温室整体划分为四面体混合网格单元10，并局部加密暖风扇出风口和植物区域，以适应气流和温度可能出现的较大变化11。  1.3.3  边界条件 因为温室密闭性较好且置于室内，不考虑太阳辐射和漏风的影响，对模型进行简化；将暖风扇侧的边界条件设为速度入口条件，吹出水平方向的 25 暖风，平均速度为 1 m/s；覆盖层全部简化为玻璃，设为壁面条件，其与整个温室进行能量交换的方式主要为对流传热和传导10，初始温度设为 20 ；为简化能量交换过程，在进行模拟时不考虑其他的传热过程。将拥有大的热容量的土壤的边界条件设为壁面12，初始温度设为 20 。  1.3.4  求解运算 在进行模拟时，为修正速度，选用 SIMPLEC 算法；梯度项选用最小二乘法；压力项采用标准算法；能量和动量采用二阶迎风；为加快收敛速度，黏性和湍流选用一阶迎风。 2 FLUENT 仿真结果与验证 2.1 FLUENT 仿真结果 用 Fluent 进行模拟的结果如图 4、 5 所示。在吹入热风后，暖风口周围区域的气体流速变化明显，且存在较大梯度。气流从暖风扇水平送入温室后，略微向下方运动，在温室内部的东面、暖风扇正下方和温室底部形成涡流，且温度与风速有着直接的关系。  图 4  Y=0.5 m 和 Y=0.1 m 处速度云图 图 5  Y=0.5 m 和 Y=0.1 m 处温度云图 选取温室内植物冠层平面 Y=0.5 m 和植物内部平面 Y=0.1 m 的速度云图（图 4）进行分析，当风速为0.301.00 m/s 时，植物生长较好13。在使用 Fluent 进行模拟计算后，在距离地面 0.5 m 处的平均空气流速为 0.152 m/s，最低速度为 0.009 m/s，低于最佳风速。而植物冠层部分没有处于最适风速的区域出现。在平面 Y=0.1 m 处气流平均流速为 0.195 m/s，最小流速为0.013 m/s，低于适宜风速的最小值 0.3 m/s；最大风速为 0.552 m/s，虽略大于 0.3 m/s，但植物内部达到适宜植物生长的风速的部分仅占整体的 20%。综上所述，此时的气体流动不适于植物的生长。  选取温室内植物表面 Y=0.5 m 和植物内部 Y=0.1 m 的温度云图（图 5）进行分析。在暖风口附近温度最高，温度与其他区域相比有很大的梯度。在 2 个平面内最高温度可 24.45 ，平均温度 23.17 ，温差的最大值为 1.28 ，此温度较适于大部分植物的生长。  2.2 对比验证 在试验期间比较 10 个监测点的实际测量值和使用 Fluent 得出的计算值，其结果可从图 6、 7 中明显看出。  图 6  风速在监测点模拟值与实测值的比较 4                                          应用科技                                 第 44 卷  图 7  温度在监测点模拟值与实测值比较根据图 6，对比不同监测点的风速，所有观测点的绝对误差的平均值为 0.008 m/s，相对误差的平均值为 5.54%，绝对误差和相对误差的最大值出现在监测点6P 处，分别为 0.027 m/s 和 24.228%。虽然温室内风速的模拟值小于测量值，但大部分情况下，模拟结果与实际情况一致。因此可以认为模型的速度计算值可以较好地体现温室内部气体流速的情况。  根据图 7，对比不同监测点的温度， 10 个监测点绝对误差均值为 0.44 ，相对误差均值为 4%。10P 处出现最大绝对误差 1.66 和最大相对误差 2.78%， 模拟值与从监测点测出的值差别不大，出现误差的原因可能是温室模型密封不够好而造成漏风，导致室温略低。结合仪器位置和人员干扰对实验结果造成的影响14，从整体来看，温度的实际值和计算值误差很小，因此可以将温度的模拟值视作温室内温度的实际值。  3 温室内环境的优化设计 3.1 优化方案设计  温室内的速度场、温度场的均匀性受暖风扇的影响较大，因此通过在温室中改变进风口的速度和改变暖风口的位置来设计出 3 种改进方案，以满足植物生长的需要。方案一，升高暖风口速度，即暖风口速度从 1 m/s 升高到 2 m/s；方案二，下移暖风口的位置，将暖风口从 1.3 m 处下移到 1 m 处，如图 8 所示；方案三，下移暖风口的位置至 0.7 m 处，如图 9 所示。将原温室作为方案零来进行对比，温室、植物和暖风扇的尺寸视为定量不做改变，且方案一、方案二和方案三的网格划分与方案零相同，其余边界条件的设置也相同。通过对不同方案使用 Fluent 进行模拟计算，比较植物冠层表面和内部的风速及温度分布情况，依次来选择拥有适宜风速温度的最适宜植物生长的方案。  图 8  暖风扇离地 1 m 时的玻璃温室模型 图 9  暖风扇离地 0.7 m 时的玻璃温室模型  3.2 优化模拟分析 3.2.1 速度场模拟分析 经过改良后的方案气流分布模式基本相同：在温室上部和底部都有涡流的产生，暖风扇口附近的气流流速最大，并随着流动方向减弱，直到在东侧壁面附近开始沿温室顶部返回，温室西侧的植物表面部分和温室壁面附近风速较小，且会在植物内部的涡流处出现整个植物区域中风速的最大值。  在方案一中，从速度云图图 10 上看，选取的 2个平面空气流速平均值为 0.695 m/s，均高于方案零的0.152 m/s，且满足植物生长的最佳风速。风速在南北两侧墙处达到最大，为 1.562 m/s，远高于最适宜植物生长的风速 1 m/s，这样的速度容易对植物造成损伤，此时植物表面风速位于 0.31.0 m/s的部分约占整个植物冠层的 63%。在植物内部 Y=0.1 m 处最大风度为1.581 m/s，平均风速 0.928 m/s，处在最佳气流流速范围内的区域约占 27.9%，剩下的大部分区域风速均大于 1 m/s，不利于植物生长，由此可知不能仅增大风速来优化温室内部的气流分布情况。  图 10  方案一 速度场模拟 从方案二的云图图 11 可知， y=0.5 m 处气流的流速均值为 0.284 m/s，最大为 0.632 m/s； y=0.1 m 处平均风速为 0.423 m/s，最小风速为 0 m/s，最大风速为0.754 m/s，均小于方案一。处于适宜风速的占 25%，虽然相对于方案零适宜风速的占比有所增加，气流的均匀性与方案一相比较好，但适宜风速的占比仍小于方案一。  第 1期                  贾鹤鸣， 韩骏骋等： 基于 CFD的玻璃温室环境数值模拟及优化分析                             5  图 11  方案二 速度场模拟 分析方案三的云图图 12， 选取分析的平面平均气体流速为 0.383 m/s，最高流速可达 0.834 m/s，均符合最佳风速。在植物部分的最低风速为 0.07 m/s，大于方案零和方案二，但仍低于 0.3m/s。这时的植物冠层风速达到适宜风速的占 70%，高于其他方案，且均匀性最好。  图 12  方案三 速度场模拟 3.2.2 温度场模拟分析 分析方案一、方案二和方案三的温度情况（图1315） ，并将其与速度云图进行比对，发现气流速度较大的区域温度上升较多。 3 个方案的最高温度分别为 24.49、 24.45 和 24.72，平均温度为 23.93、 23.15和 23.49， 均高于方案零， 且植物区域的温度分布均匀性优于方案零。就温度而言，方案一的均值和最大值都为最高，其次为方案三。  图 13  方案一 温度场模拟 图 14  方案二 温度场模拟 图 15  方案三 温度场模拟 通过模拟分析实验中设计的 3 种温室的优化方案，可以发现： 1）方案一的气流均匀性稍差，且有较多的区域风速高于最佳风速，不利于植物生长； 2）方案二虽然略优于原始方案，但对于气流场和温度场并不能提供有效的优化； 3）方案三的气流均匀性最好，不存在流速大于 1 m/s 的区域，从而不会对植物的正常生长造成威胁，且风速大于方案零和方案二。总体来讲，方案三的暖风扇位置更适宜植物的生长。  4  结论  1）温度数值模拟结果与实验实测结果吻合良好，表明建立的 CFD 模型有效， 模拟玻璃温室内环境分布是可行的；  2）设计的种结构优化方案的风速值及均匀性均优于方案零。 综合考虑， 方案三气流循环模式最佳，气流均匀性最好，适宜植物生长。  本文仅基于风速和暖风扇位置进行温室的优化，在之后的研究中，可考虑将添加排气扇和增加暖风扇的数量作为温室的优化方式，考虑植物在阳光下的光合作用产生的热量和外部辐射等会影响温室内部气流和温度的条件， 探究温室内部气流与温度的分布情况。  参考文献 :  1 程秀花，毛罕平，伍德林等. 栽有番茄的玻璃温室内气流场分布CFD 数值模拟 J. 江苏 大学学报， 2010， 31（ 5） : 510-514. 2 汪郑邦 . 现代温室供热系统的设计及运行模拟研究 D. 哈尔滨工业大学， 2006：79-80. 3 刘文合 ,杨英英 ,徐占洋等 . 基于 CFD方法的温室散热系统结构优化模拟研究 J .灌溉排 水学报. 2016(04)：8-13. 4 周伟 ,汪小旵 ,李永博 . 可变边界条件下的 Venlo 温室温度场三维非稳态模拟J .农业机械学报， 2014， 45(11)： 6                                          应用科技                                 第 44 卷  304-310. 5 张伟建 ,王新忠 ,李亮亮等 . 夏季屋顶全开型温室遮阳网降温调控的CFD 分析J .江苏农业科学 , 2017(19)： 253-256. 6 刘焕 .基于 CFD 的人工光型植物工厂通风模拟与优化研究D.中国农业科学院 . 2018:10-29. 7 段明辉 . 冬季温室数值分析研究D，东北农业大学, 2014：7-9. 8 宿文等 . 自然通风对日光温室气温影响的模拟分析 J， 南京信息工程大学 , 2016： 1635-1642. 9 张艺萌 . 基于CFD的寒地水稻育秧大棚环境数值模拟分析与研究D.黑龙 江八一农垦大学, 2018， 28-54 10 赵杰强 , 赵云 . 机械通风连栋温室的温度场 CFD 模拟 J.中国农机化学报 , 2014,6: 76-79. 11 何国敏，现代化温室温度场数字化模拟研究D，南京农业大学， 2011：50-67. 12 沈敏， 食用菌工厂菇房内环境的模拟验证与优化研究 D, 江苏大学 ,2016:23-41. 13 吴飞青 .温室冬季热环境数值分析与实验研究D，浙江工业大学， 2010:41-71. 14 刘焕 .基于 CFD 的人工光植物工厂气流场和温度场的模拟及优化D，中 国农业科学院， 2017:4. 15 张芳，方慧，杨其长等基于模型的大跨度温室自然通风热环境模拟J中国农业气象，2017,38(4):221-229. 本文引用格式： 贾鹤鸣 , 韩骏骋 , 张森 , 等 . 基于 CFD 的玻璃温室环境数值模拟及优化分析 J. 应用科技 , DOI: 10.11991/yykj. 201901005. JIA Heming, HAN Juncheng, ZHANG Sen, et al. Numerical analysis simulation and optimization of the environment in glass greenhouse based on CFDJ. Applied science and technology, DOI: 10.11991/yykj. 201901005.