农业工程学报
農業工程學報
농업공정학보
2015年
12期
194-199
,共6页
刘妍华%曾志雄%郭嘉明%吕恩利%孟庆林
劉妍華%曾誌雄%郭嘉明%呂恩利%孟慶林
류연화%증지웅%곽가명%려은리%맹경림
二氧化碳%温室%数值方法%气肥增施
二氧化碳%溫室%數值方法%氣肥增施
이양화탄%온실%수치방법%기비증시
carbon dioxide%greenhouses%numerical methods%air fertilizer enrichment
CO2气体是温室作物光合作用的重要原料之一。为掌握温室CO2气肥增施性能,以温室CO2气肥增施为研究对象,采用计算流体力学(computational fluid dynamics,CFD)技术建立二维紊流数值计算模型。基于FLUENT软件,结合SIMPLE算法,采用有孔介质模型、k?ε模型、离心坐标(discrete ordinates,DO)模型,添加组分传输模型,对气肥喷射高度、气肥增施流量等因素对CO2增施性能进行数值模拟,得到温室内CO2浓度变化和分布规律。研究结果表明:气肥增施过程对作物区的温度场影响较小,温室的作物区域最大温度与最小温度差值不超过0.5℃,作物区域的气流流场以及温度场分布较为均匀;由于CO2的沉积效应,温室下部区域的CO2浓度相对较高;气肥喷射高度越高,CO2扩散的范围越大,沉积在作物区的CO2相对越少,CO2的浓度也相对较低;气肥增施流量越大,作物区域的CO2浓度上升越快。试验结果表明,CO2浓度模拟值与试验值差异不大于5%,模拟结果与试验结果较吻合,证明了模型的正确性。该研究对掌握温室CO2气肥增施性能的流场变化规律,开展温室气肥增施装备的优化设计具有一定的参考价值。
CO2氣體是溫室作物光閤作用的重要原料之一。為掌握溫室CO2氣肥增施性能,以溫室CO2氣肥增施為研究對象,採用計算流體力學(computational fluid dynamics,CFD)技術建立二維紊流數值計算模型。基于FLUENT軟件,結閤SIMPLE算法,採用有孔介質模型、k?ε模型、離心坐標(discrete ordinates,DO)模型,添加組分傳輸模型,對氣肥噴射高度、氣肥增施流量等因素對CO2增施性能進行數值模擬,得到溫室內CO2濃度變化和分佈規律。研究結果錶明:氣肥增施過程對作物區的溫度場影響較小,溫室的作物區域最大溫度與最小溫度差值不超過0.5℃,作物區域的氣流流場以及溫度場分佈較為均勻;由于CO2的沉積效應,溫室下部區域的CO2濃度相對較高;氣肥噴射高度越高,CO2擴散的範圍越大,沉積在作物區的CO2相對越少,CO2的濃度也相對較低;氣肥增施流量越大,作物區域的CO2濃度上升越快。試驗結果錶明,CO2濃度模擬值與試驗值差異不大于5%,模擬結果與試驗結果較吻閤,證明瞭模型的正確性。該研究對掌握溫室CO2氣肥增施性能的流場變化規律,開展溫室氣肥增施裝備的優化設計具有一定的參攷價值。
CO2기체시온실작물광합작용적중요원료지일。위장악온실CO2기비증시성능,이온실CO2기비증시위연구대상,채용계산류체역학(computational fluid dynamics,CFD)기술건립이유문류수치계산모형。기우FLUENT연건,결합SIMPLE산법,채용유공개질모형、k?ε모형、리심좌표(discrete ordinates,DO)모형,첨가조분전수모형,대기비분사고도、기비증시류량등인소대CO2증시성능진행수치모의,득도온실내CO2농도변화화분포규률。연구결과표명:기비증시과정대작물구적온도장영향교소,온실적작물구역최대온도여최소온도차치불초과0.5℃,작물구역적기류류장이급온도장분포교위균균;유우CO2적침적효응,온실하부구역적CO2농도상대교고;기비분사고도월고,CO2확산적범위월대,침적재작물구적CO2상대월소,CO2적농도야상대교저;기비증시류량월대,작물구역적CO2농도상승월쾌。시험결과표명,CO2농도모의치여시험치차이불대우5%,모의결과여시험결과교문합,증명료모형적정학성。해연구대장악온실CO2기비증시성능적류장변화규률,개전온실기비증시장비적우화설계구유일정적삼고개치。
Carbon dioxide (CO2) is one of the important raw materials for photosynthesis of crops in greenhouse, which can increase the harvest of the crops. In order to analyze the performance of CO2 fertilizer enrichment, this paper took indoor environment of greenhouse as the research object, and by means of unstructured mesh generation method, a two-dimensional turbulence computational model of greenhouse was built by ICEM CFD (the integrated computer engineering and manufacturing code for computational fluid dynamics) technique. After building 7532 triangle/quadrangle meshes, the skewness of mesh model was less than 0.75. To solve two-dimensional computational model of greenhouse, the CFD software FLUENT and the SIMPLE algorithm were used. Meanwhile, the porous model, the DO (discrete ordinates) model and thek-ε model were adopted. Governing equations of finite volume method were employed, including mass, momentum and energy conservation equations. In computational model, CO2 enriching inlet was set with the condition of mass flow inlet while outlet was set with the condition of pressure outlet through pre-calculation. And crop area was defined as porous material with the porosity of 0.8. Adopting steady-state solver to operate, environmental parameters such as indoor temperature, nitrogen (N2) concentration, oxygen (O2) concentration and CO2concentration should be initialized when flow field of greenhouse was stable. Then transient solver was needed to numerical simulation when acceleration of gravity was 9.8 m/s2 and step size was 0.01 s. By such methods, the effects of CO2jetting height, enriching flow speed and other factors on CO2 enrichment property were computed while the change of CO2concentration and distribution regularities were studied. After the simulations, some results were obtained. CO2 enrichment process had little effect on temperature of crop area, and both temperature distribution and air velocity distribution of crop area were uniform. As CO2 enrichment process lasted only for 5 min and solar radiation affected indoor environment of the greenhouse rarely, the temperature difference in the crop area was less than 0.5℃. On the other hand, cooler air was more likely to deposit to the bottom of the greenhouse and hotter air exhausted from wet curtain because of density difference of air flow. With the influence of buoyancy lift and newly added flow, air velocity of bottom side was higher. And the biggest velocity of air flow appeared in the area between jetting nozzle and bottom side of the greenhouse. Molar mass of CO2is greater than O2 and N2. And as a result of deposition effect of CO2, the bottom of the greenhouse had a higher CO2concentration while CO2 concentration in upper part was less, which helped the photosynthesis of crops. Moreover, CO2flowed to the exit of the greenhouse via bottom side. In crop area, CO2concentration was higher when closing to entrance. Although there was nearly no different in first 30 s, CO2jetting height had a certain effect on CO2enrichment time. CO2 concentration could be decreased when CO2jetting height was located at a high level. When jetting height was higher, spreading range of CO2was larger, CO2 deposited in crop area was less and CO2concentration was lower. However, CO2 concentration of crop area would be increased too fast when jetting height was too low. With the increase of CO2 enrichment flow, CO2concentration increased and CO2 enrichment time decreased. CO2 enrichment flow had a certain effect on CO2 enrichment time in crop area. Thus, transforming enrichment flow or time could control CO2concentration of the greenhouse which benefited photosynthesis of crops. A test was made to verify the accuracy of the numerical model. The simulation values and test values of CO2concentration at the same moment were in good agreement, and the biggest difference between them was less than 5%, which demonstrated that such models were correct. Result of this paper reveals the CO2 enrichment rules and the characteristics of flow field in greenhouse, which has a certain reference for the design and optimization of CO2 enriching equipments in greenhouse.