王金翔,謝婭婭

基于有限元的蔬菜大棚通風(fēng)模擬研究

王金翔,謝婭婭*

荊楚理工學(xué)院, 湖北 荊門 448000

本文基于有限元仿真模型模擬了不同通風(fēng)方式下大棚內(nèi)的空氣流場分布規(guī)律，實測數(shù)據(jù)驗證了模型的準(zhǔn)確性。研究發(fā)現(xiàn)出風(fēng)口與進風(fēng)口相對時，靠近進風(fēng)口附近可種植的蔬菜高度高于出風(fēng)口附近。出風(fēng)口與進風(fēng)口不相對時，靠近出風(fēng)口一定范圍內(nèi)的空氣速度較為平穩(wěn)，有利于蔬菜種植。頂部通風(fēng)有助于改變大棚頂部的空氣循環(huán)模式。在遠(yuǎn)離進風(fēng)口的頂部設(shè)置出風(fēng)口可促進空氣流動速度增加。本文的研究結(jié)果可為蔬菜大棚通風(fēng)口設(shè)計提供理論指導(dǎo)。

有限元模型; 蔬菜大棚; 通風(fēng)模擬

隨著生活水平的提高，人們對各種各樣的蔬菜的需求量增加，餐桌上時常可見反季節(jié)的蔬菜，這都?xì)w功于蔬菜大棚大面積推廣應(yīng)用。由于大棚能夠改變蔬菜生長的微氣候，為其營造適宜的生長環(huán)境，無需根據(jù)天氣選擇種植的蔬菜品種。通風(fēng)在大棚內(nèi)部微氣候調(diào)控中起到非常重要的作用。通風(fēng)的作用是交換熱量、降低濕度，為大棚提供新鮮空氣保證其內(nèi)部的二氧化碳濃度。但如果大棚通風(fēng)不暢，會出現(xiàn)溫度過高、濕度過大、二氧化碳濃度偏低和有害氣體積累等問題，影響蔬菜的正常生長，造成病蟲害嚴(yán)重、作物減產(chǎn)、品質(zhì)低劣。因此，通風(fēng)方式對大棚環(huán)境調(diào)控規(guī)律研究尤為重要[1-3]。

隨著計算機技術(shù)的發(fā)展，基于有限元的計算流體力學(xué)模擬（CFD）被廣泛用于流體流動規(guī)律研究。大棚由于尺寸較大，實驗研究存在一定的難度和誤差，因此眾多研究者利用CFD研究大棚在應(yīng)用過程中的空氣流動和溫度分布。何科奭等[4]利用計算流體力學(xué)模擬的方式研究不同通風(fēng)方式下溫室內(nèi)氣流和溫度場，探究通風(fēng)方式對通風(fēng)率、室內(nèi)外溫差和室內(nèi)氣候均勻性的影響規(guī)律，揭示溫室微氣候形成機理?，F(xiàn)有的通風(fēng)方式主要有機械通風(fēng)和自然通風(fēng)，研究者研究了不同通風(fēng)方式的特征。趙杰強等[5]、黃全豐等[6]利用Fluent模擬了機械通風(fēng)形式下，連棟大棚內(nèi)的溫度分布規(guī)律。部分研究者[7-11]通過Fluent模擬了不同通風(fēng)方式下大棚內(nèi)的空氣流動和溫度場分布，研究大棚的通風(fēng)操作方式。目前的研究成果[12-20]均是針對不同類型的大棚進行模擬，未充分討論研究不同通風(fēng)模式下大棚內(nèi)的流動規(guī)律，以及根據(jù)該流動規(guī)律選擇種植事宜的蔬菜品種。

本文建立有限元仿真模型，利用Fluent軟件模擬了不同通風(fēng)方式下大棚內(nèi)的空氣流場分布規(guī)律，對比分析了通風(fēng)口與進風(fēng)口相對位置不同的影響，探究了不同頂部通風(fēng)口布置方式對大棚內(nèi)的流場影響規(guī)律。本文的研究結(jié)果可為大棚通風(fēng)口的設(shè)計和大棚內(nèi)蔬菜種植方式提供指導(dǎo)。

1 模擬模型及對象

1.1 模型方程

在進行通風(fēng)操作時，溫室內(nèi)空氣流動過程中產(chǎn)生的雷諾數(shù)Re一般大于4 000，具有較高的湍流特性，因此本文選用-湍流模型來描述空氣流動。同時本文假設(shè)大棚內(nèi)氣體為連續(xù)、穩(wěn)定、不可壓縮牛頓流體，因此用有限元的方法求解空氣流動問題時，控制方程滿足連續(xù)性方程、動能方程、湍流能量耗散方程和湍流能量耗散率方程如下。

（1）連續(xù)性方程：

（2）動量方程：

（3）湍流能量耗散方程和能量耗散率方程

1.2 模擬對象及方法

本文選擇華中地區(qū)典型的金屬骨架塑料溫室大棚為研究對象，寬度8 m，長度20 m，脊高3 m，頂高5.0 m，主進風(fēng)口寬0.8 m。采用非結(jié)構(gòu)化網(wǎng)格對大棚進行劃分，進出風(fēng)口進行了網(wǎng)格加密。進口風(fēng)速為1.2 m?s-1。模擬計算時采用的邊界條件和初始條件如表1所示。

表 1 模擬邊界條件和初始條件

2 結(jié)果與分析

2.1 網(wǎng)格無關(guān)性檢驗

為了避免網(wǎng)格大小影響模擬計算結(jié)果準(zhǔn)確性，本文采用3種網(wǎng)格尺寸對進風(fēng)口和出風(fēng)口相對的大棚進行模擬研究。各種網(wǎng)格尺寸大小下，通風(fēng)30 min之后沿著大棚長度方向距地面2 m高度處的空氣速度如圖1所示。

圖 1 網(wǎng)格無關(guān)性檢驗

圖1中，曲線為模擬值，星號為實際測量值。由圖可見，不網(wǎng)格大小模擬獲得的空氣速度沿長度方向變化趨勢基本一致，而且與實測值之間誤差較小。因此，模擬過程采用的模型準(zhǔn)確，可很好地反應(yīng)實際通風(fēng)過程，同時網(wǎng)格尺寸對模擬結(jié)果影響較小。后文的模擬中采用的網(wǎng)格尺寸為0.05 m。

2.2 進風(fēng)口位置及數(shù)量對蔬菜大棚內(nèi)空氣流動的影響

進風(fēng)口和出風(fēng)口的相對位置會影響通風(fēng)時空氣在蔬菜大棚內(nèi)的流動軌跡，進而影響不同高度處蔬菜和空氣接觸的量，因此本節(jié)考察了三種不同出風(fēng)口位置對大棚內(nèi)的空氣流動影響。圖2所示的是不同出風(fēng)方式下大棚沿長度方向的中軸面上的空氣速度分布，其中middle代表出風(fēng)口與進風(fēng)口相對，lower代表出風(fēng)口在進風(fēng)口相對位置的下方，upper代表出風(fēng)口在進風(fēng)口相對位置的上方，出風(fēng)口的尺寸均為0.8 m。由圖2（a）可見，出風(fēng)口與進風(fēng)口相對時，空氣達到大棚中心位置后的速度可以保持在0.5 m?s-1左右，之后速度逐漸減小。在出口附近高度方向上空氣的速度梯度較大。因此，圖2（a）的通風(fēng)方式下，靠近進風(fēng)口的地方比較適合種植高度與進風(fēng)口高度差不多的蔬菜品種，如黃瓜、西紅柿等。而在出風(fēng)口附近比較適合種植低于出風(fēng)口的蔬菜品種，如芹菜、青菜、豌豆等，這樣增加不同種類的蔬菜與空氣接觸的機會。由圖2（b）和（c）可見，當(dāng)出風(fēng)口在進風(fēng)口相對位置的下方和上方時，另外一個在距地面1 m的位置時，空氣在未到達大棚中心位置前，其速度已降至0.4 m?s-1左右，之后風(fēng)速繼續(xù)降低且在高度方向上速度梯度較小?？拷M風(fēng)口的地方比較適合種植高度與進風(fēng)口高度差不多，且需要通風(fēng)量較大的蔬菜品種。而在出風(fēng)口附近比較適合種植低于出風(fēng)口的蔬菜品種，或與進風(fēng)口高度差不多但需要的通風(fēng)量較小的蔬菜品種，這樣增加不同種類的蔬菜與空氣接觸的機會。

圖 2 不同出風(fēng)方式下大棚沿長度方向的中軸面上的空氣速度分布

圖3所示的是不同出風(fēng)方式下大棚沿長度方向的中軸面上的空氣速度矢量分布。由圖可見，在大棚的頂部形成一個空氣流動大循環(huán)，不利于空氣的交換，因此進風(fēng)口以上的位置不適合種植過高的蔬菜品種?？拷M風(fēng)口2 m左右的位置存在一個空氣流動小循環(huán)，不利于空氣的交換，因此靠近進風(fēng)口不適宜種植蔬菜。

圖 3 不同出風(fēng)方式下大棚沿長度方向的中軸面上的空氣速度矢量分布

圖4所示的是三種不同出風(fēng)口設(shè)置方式下0.8 m和2 m高度處空氣速度沿長度方向分布。由圖可見，出風(fēng)口與進風(fēng)口不相對時，大棚內(nèi)的空氣速度比兩個風(fēng)口相對時較小，主要是由于空氣運動的軌跡有一定的增加。但是出風(fēng)口與進風(fēng)口不相對時，靠近出風(fēng)口一定范圍內(nèi)的空氣速度較為平穩(wěn)。

圖 4 三種不同出風(fēng)口設(shè)置方式下不同高度處空氣速度沿長度方向分布

圖5所示的是三種不同出風(fēng)口設(shè)置方式下長度5 m、10 m、15 m處空氣速度沿高度方向分布。由圖可見，三種不同出風(fēng)口下，大棚長度L=5 m、高度H=3 m左右存在一個空氣速度非常小的區(qū)域，這是由于大棚上不氣流大循環(huán)引起的，因此該高度范圍內(nèi)不適宜種植蔬菜。大棚長度L＞10 m后空氣的速度沿軸向分布梯度較小。

圖 5 三種不同出風(fēng)口設(shè)置方式下不同長度位置處空氣速度沿高度方向分布

2.3 頂部出風(fēng)口打開方式對大棚內(nèi)空氣流動的影響

圖6和圖7分別是不同頂部出風(fēng)方式下大棚沿長度方向的中軸面上的空氣速度和空氣速度矢量分布，四種通風(fēng)方式分別為：1個頂部中間出風(fēng)口、1個頂部中間出風(fēng)口和1個左側(cè)出風(fēng)口、1個頂部中間出風(fēng)口和1個右側(cè)出風(fēng)口、1個頂部中間出風(fēng)口和1個左側(cè)出風(fēng)口和1個右側(cè)出風(fēng)口，出風(fēng)口的尺寸均為0.8 m。由圖可見，頂部出風(fēng)方式對大棚2 m高度以下范圍內(nèi)的空氣流速影響較小，這主要是由于空氣沿著水平方向經(jīng)過大棚且進風(fēng)口中心距地面高度為1.4 m，大棚頂部的空氣流速較小。由圖6可見，在遠(yuǎn)離進風(fēng)口的頂部設(shè)置出風(fēng)口可促進遠(yuǎn)離進風(fēng)口端的位置處的空氣流動，提高空氣流速，有助于蔬菜的種植。由圖7可見，頂部通風(fēng)有助于改變大棚頂部的空氣循環(huán)模式，有助于頂部空氣置換出大棚。

圖 6 不同頂部出風(fēng)方式下大棚沿長度方向的中軸面上的空氣速度分布

圖 7 不同頂部出風(fēng)方式下大棚沿大棚長度方向的中軸面上的空氣速度矢量分布

3 結(jié) 論

本文基于有限元仿真模型模擬了不同通風(fēng)方式下蔬菜大棚內(nèi)的空氣流場分布規(guī)律。模擬獲得的空氣速度分布結(jié)果和實際測量值對比說明模型的準(zhǔn)確性。分別模擬了通風(fēng)口與進風(fēng)口不同相對位置下大棚內(nèi)的流場分布，結(jié)果說明出風(fēng)口與進風(fēng)口相對的通風(fēng)方式下，靠近進風(fēng)口的地方比較適合種植高度與進風(fēng)口高度差不多的蔬菜品種，如黃瓜、西紅柿等。而在出風(fēng)口附近比較適合種植低于出風(fēng)口的蔬菜品種，如芹菜、青菜、豌豆等，這樣增加不同種類的蔬菜與空氣接觸的機會。出風(fēng)口與進風(fēng)口不相對時，靠近出風(fēng)口一定范圍內(nèi)的空氣速度較為平穩(wěn)。同時探究了不同頂部通風(fēng)口布置方式對大棚內(nèi)的流場影響規(guī)律。頂部通風(fēng)有助于改變大棚頂部的空氣循環(huán)模式，有助于頂部空氣置換出大棚。在遠(yuǎn)離進風(fēng)口的頂部設(shè)置出風(fēng)口可促進遠(yuǎn)離進風(fēng)口端的位置處的空氣流動，提高空氣流速，有助于蔬菜的種植。本文的研究結(jié)果說明需要根據(jù)不同通風(fēng)方式下大棚內(nèi)的流場分布確定種植的蔬菜類型，同時頂部設(shè)置通風(fēng)口有助于通風(fēng)。

[1] 張同鋒.蔬菜大棚微氣候環(huán)境的實驗與數(shù)值模擬研究[D].西安:西安建筑科技大學(xué),2006

[2] 孟凡英,劉美.自然通風(fēng)單棟溫室內(nèi)流場的CFD模擬[J].黑龍江農(nóng)業(yè)科學(xué),2012(2):128-131

[3] Stefano B, Patrizia T, Filippo B,. Efficacy of greenhouse natural ventilation: Environmental monitoring and CFD simulations of a study case [J]. Energy and Buildings, 2016,125:276-286

[4] 何科奭,陳大躍,孫麗娟,等.不同風(fēng)況和開窗配置對夏季單棟塑料溫室微氣候的影響[J].農(nóng)業(yè)機械學(xué)報,2017,48(12):311-318,339

[5] 趙杰強,趙云.機械通風(fēng)連棟溫室的溫度場CFD模擬[J].中國農(nóng)機化學(xué)報,2014,35(6):76-79

[6] 黃全豐,趙云.機械通風(fēng)溫室夏季內(nèi)部溫度場及動態(tài)特性CFD模擬[J].農(nóng)機化研究,2013,35(1):30-33

[7] 劉煥.基于CFD的人工光型植物工廠通風(fēng)模擬與優(yōu)化研究[D].北京:中國農(nóng)業(yè)科學(xué)院,2018

[8] Eduardo, Rosales-Vicelis J, Efrén, Fitz-Rodríguez, López-Diaz, J. CFD modeling of the environment of a naturally ventilated greenhouse for central Mexico conditions [C]. Humberto Source: ASABE 2018 Annual International Meeting, 2018

[9] Lukas S, Madjid M. Geometric optimization of thermally activated small scale solar chimneys for natural ventilation applying CFD simulations and measurements [C]. Thomas Source: 30th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems, ECOS 2017, 2017, 30th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems, ECOS 2017.

[10] Wang XZ, Zhang WJ, Zhang L,. Analysis of flow field for full open-roof glass greenhouse with nature ventilation in summer based on CFD [J].Transactions of the Chinese Society for Agricultural Machinery, 2016,47(10):332-337

[11] 張芳.基于CFD技術(shù)的大跨度溫室噴霧降溫?zé)岘h(huán)境模擬[D].北京:中國農(nóng)業(yè)科學(xué)院,2017

[12] 曹陽.蔬菜大棚內(nèi)部流場三維紊流數(shù)值模擬及大棚斷面優(yōu)化研究[D].揚州:揚州大學(xué),2004

[13] Kim RW, Hong SW, Norton T,. Computational fluid dynamics for non-experts: Development of a user-friendly CFD simulator (HNVR-SYS) for natural ventilation design applications [J]. Biosystems Engineering, 2020,193:232-246

[14] Villagrán EA, Romero EJB, Bojacá CR. Transient CFD analysis of the natural ventilation of three types of greenhouses used for agricultural production in a tropical mountain climate [J]. Biosystems Engineering, 2019,188:288-304

[15] 薛曉萍,宿文.基于CFD的自然通風(fēng)對日光溫室濕度分布模擬分析[J].海洋氣象學(xué)報,2019,39(4):90-96

[16] Saberian A, Sajadiye SM. The effect of dynamic solar heat load on the greenhouse microclimate using CFD simulation [J]. Renewable Energy, 2019,138:722-737

[17] 尹奎超.某溫室大棚自然通風(fēng)的CFD模擬分析[J].建筑熱能通風(fēng)空調(diào),2020,39(4):88-91

[18] 何國敏,汪小旵,孫國祥.現(xiàn)代化溫室自然通風(fēng)時濕熱環(huán)境CFD模擬研究[J].西南大學(xué)學(xué)報(自然科學(xué)版),2011,33(9):136-141

[19] Sun XY, Yang H, Liu QF,. Event-triggered control for greenhouse temperature under natural ventilation based on computational fluid dynamics [J]. Systems Science and Control Engineering, 2021,9(S2):93-101

[20] Lebbal C, Bougoul S, Zeroual S. Simulation of natural ventilation inside tunnel greenhouse [J]. International Journal of Heat and Technology, 2020,38(3):752-757

Simulation for the Ventilation Process in Vegetable Greenhouses Based on Finite Element Method

WANG Jin-xiang, XIE Ya-ya*

448000,

This paper simulates the distribution of air flow field in greenhouses under different ventilation modes based on finite element simulation model. The accuracy of the model was verified by experimental data. It was found that when the air outlet was opposite to the air inlet, the height of vegetables planted near the air inlet was higher than that near the air outlet. When the air outlet is not relative to the air inlet, the air velocity within a certain range near the air outlet is relatively stable. Top ventilation helps to change the air circulation pattern at the top of the greenhouse. The air outlet at the top far from the air inlet can promote the increase of air flow speed. The results of this study can provide theoretical guidance for the design of greenhouse vents.

Finite element method; vegetable greenhouse; simulation ventilation

S625.51

A

1000-2324(2022)02-0314-06

10.3969/j.issn.1000-2324.2022.02.020

2021-11-05

2021-12-06

湖北省教育廳科學(xué)研究計劃項目(B2021263);荊楚理工學(xué)院校級科研項目(YD202101)

王金翔(1982-),男,碩士,講師,主要研究方向:通信與信息系統(tǒng)、機器人工程. E-mail:wjx1982@126.com

Author for correspondence. E-mail:33353302@qq.com