蔬菜氣霧栽培箱不同氣流循環方式的流場和溫度場CFD模擬

蔣 蘋，楊希文，羅亞輝，龍莉霞，石毅新，胡文武

蔬菜氣霧栽培箱不同氣流循環方式的流場和溫度場CFD模擬

蔣 蘋1,2，楊希文1，羅亞輝1,2，龍莉霞1，石毅新1,2，胡文武1,2

（1. 湖南農業大學工學院，長沙 410128；2. 南方糧油作物協同創新中心，長沙 410128）

蔬菜氣霧栽培箱內的空氣流動和溫度對箱體內的環境調節及農作物生長具有重要作用，農作物周圍空氣的均勻性流動能促進農作物的生長速率。為探究蔬菜氣霧栽培箱內流場及溫度場分布規律，基于計算機流體力學（CFD，computational fluid dynamics）方法，利用FLUENT軟件，結合標準湍流模型、有孔介質模型、作物冠層質熱交換模型等，建立了蔬菜氣霧栽培箱不同氣流循環方式下的CFD模型。并對氣霧栽培箱內的環境進行優化設計，設計了3種氣流循環方案：頂面進側面出，側面進頂面出，側面進側面出。對送回風口的不同位置布局進行了研究，并對3種氣流循環方案進行了數值模擬。模擬結果可知：氣流為頂面進側面出方案中，風速位于生菜生長適宜風速值區域占58.1%，適宜溫度值區域占93.6%，通風死角區域占比0.844%；氣流為側面進頂面出方案中，生菜生長適宜風速值區域占59.6%，適宜溫度值區域占99.98% ，通風死角區域占比0.069%；氣流為側面進側面出方案中，風速位于生菜生長適宜風速值區域占54.3%，適宜溫度值區域占92.4%，通風死角區域占比16.7%。分析對比后得到側面進頂面出為最佳氣流循環方案。并對此進行了試驗測試，結果表明：氣霧栽培箱內溫度、風速模擬值和實測值進行對比，溫度平均相對誤差為3.9%，均方根誤差為0.86℃。風速平均相對誤差為3.5%，均方根誤差為0.26 m/s，模擬值和實測值誤差較小，模擬效果良好，驗證了CFD模型的準確性。該研究為蔬菜氣霧栽培箱內的流場及溫度變化規律，內部環境調節，裝置優化設計提供了參考依據。

風速；溫度；計算機流體力學；優化設計；蔬菜氣霧栽培箱

0 引 言

農作物氣霧栽培是無土栽培領域一種新型栽培方式，它不僅能提高作物水肥利用率，還能改善作物根系缺氧，提高土地利用率，減輕環境污染等。氣霧栽培模式主要有立柱式、金字塔式、梯形式[1]。目前部分學者對氣霧栽培設施也有相關研究，如高建民等[2]、劉義飛[3]、劉巖[4]通過設計新型的氣霧栽培設施進行農作物的氣霧栽培生產，以此提高農作物產量和品質。本課題組設計了一種蔬菜氣霧栽培箱裝置，蔬菜氣霧栽培箱是一種密閉性農作物栽培生產系統，適宜在溫室大棚、植物工廠、家庭陽臺等場所對農作物進行栽培生產。影響蔬菜氣霧栽培箱內農作物生長的主要因素有氣流速度、溫度、濕度、光照等。而氣流速度和溫度是影響蔬菜氣霧栽培箱內農作物冠層生長的重要因素，氣流速度和溫度對農作物的蒸騰作用及光合作用有著重要作用，影響農作物的品質和產量。研究氣流速度和溫度的變化規律對推廣氣霧栽培技術有重要意義。目前CFD（CFD，computational fluid dynamics）主要在溫室、塑料大棚、植物工廠等相關領域有研究。在氣霧栽培領域采用CFD方法進行流場和溫度場的數值模擬缺乏理論研究。采用CFD數值模擬方法對提高氣霧栽培箱內部農作物栽培的生長效率有著重要作用。目前相關研究有：王新忠等[5]利用CFD模型對全開型玻璃溫室自然通風進行了流場分析，結果表明側窗為45°，天窗調整至60°時，溫室整體溫度從38.4 ℃降至36.9 ℃，能達到較好的降溫效果。翁衛兵等[6]采用CFD方法對冷藏車廂內流場的流動特性進行了模擬，結果表明根據特定夾角開啟隔熱氣簾風機同時調節風向，能有效改善溫度場均勻性分布程度。林加勇等[7]對公豬舍內風速場和溫度場環境分布規律進行模擬研究，結果表明風速場模擬值與實測值誤差較小，模擬裝豬時的豬舍，溫度分布結構與裝豬前相似。高云等[8]對多環境參數控制的豬養殖箱內氣流場進行了數值模擬，結果表明養殖箱內部的氣流場形成循環，適宜通風且無死角，養殖箱能夠在較長時間穩定運行。大量文獻研究都證明了流場和溫度場對農作物的生長有著極大影響。目前研究中尚未使用CFD方法用于氣霧栽培種植技術中，本文在前人研究基礎上使用CFD方法確定農作物在氣霧栽培中的適宜風速大小及溫度范圍。建立的CFD模型能為氣霧栽培技術的環境建模和控制提供理論參考，彌補物理試驗的不足。

為了調節蔬菜氣霧栽培箱體氣流及溫度分布的均勻性，合理布置氣霧栽培箱進風口和出風口，避免氣流在氣霧栽培箱內部及農作物葉片周圍產生繞流，影響農作物正常的光合作用和蒸騰作用，該文采用計算機流體力學方法對氣霧栽培箱內部流場和溫度場進行數值模擬，并通過試驗驗證氣霧栽培箱內部流場及溫度場分布規律。同時通過合理設置進風口速度大小及進出風口位置，使氣霧栽培箱內的氣流速度與溫度兩者達到適宜及均勻的狀態，讓生菜在氣霧栽培箱內的生長環境處于最優狀態，達到高質高產的目的，以期為氣霧栽培箱的設計及優化提供理論參考。

1 材料與方法

1.1 試驗裝置

本課題組設計的蔬菜氣霧栽培箱裝置結構如圖1所示。外部形狀為一個正六棱，尺寸邊長0.35 m，高度1.6 m。氣霧栽培箱框架由不銹鋼方鋼焊接而成，頂面、底面及四周用透明亞克力板（PMMA，polymeric methyl methacrylate）進行封裝。內部栽培架則由6塊擠塑板組成，每一塊擠塑板的尺寸長1.6 m，寬0.35 m。農作物種植在擠塑板（XPS，extruded polystyrene）的定植孔內。氣霧栽培箱底部營養液槽里面的營養液通過增壓泵增壓，把營養液通過頂部的高壓霧化噴頭噴射到農作物根系，以此提供給農作物生長所需的養分和水分。氣霧栽培箱內的送風系統采用頂面送回風或側面送回風的方式，頂面送回風口呈120°角均勻布置3個，尺寸長為6 cm，寬為6 cm，對應為圖1中頂板上3個正方形開口；側面送回風口呈120°角均勻布置在3塊亞克力板上，每塊亞克力板上均勻布置3個送回風口，尺寸長6 cm，寬6 cm，對應為圖1中亞克力板上9個正方形開口。

1. 亞克力底板 2.過濾網 3.XPS定植板 4.定植孔 5.防水遮陽膜 6.不銹鋼支撐架 7.亞克力板 8.頂板 9.排風扇 10.不銹鋼環形噴霧圓環 11.霧化噴頭 12.進回風口 13.PVC管 14.增壓泵 15.營養液槽

1.2 試驗條件

試驗在湖南農業大學工學院智能氣候室內進行，氣霧栽培箱位于智能氣候室內。以智能氣候室為環境背景倉，對農作物生長環境參數進行設置，模擬出外界環境真實值。本文以蔬菜氣霧栽培箱為研究對象，以生菜為試材。對于生菜生長所需溫度，根據張成波等[9]研究表明生菜發芽期適宜空氣溫度為18～20 ℃，幼苗期光期空氣溫度為15～20 ℃，采收期光期空氣溫度18～20 ℃。根據Kitaya，Thongbai等[10-12]研究可知，氣流速度在0.3～1.3 m/s范圍內適宜植物生長，植物蒸騰作用及光合作用隨著氣流速度的增大而增強，植物周圍空氣流動能促進植物生長。以智能氣候室為環境背景倉，設置適合生菜生長所需的溫濕度，氣霧栽培箱外部環境溫度設置為18 ℃，相對濕度控制在85%，光照時間間隔設置為12 h。

1.3 試驗方法

本文以蔬菜氣霧栽培箱為研究對象，以生菜為試材。對蔬菜氣霧栽培箱送回風口的不同位置布局設計出了3種氣流循環方案，利用FLUENT軟件，對蔬菜氣霧栽培箱內的流場和溫度場進行數值模擬，得到最優方案。然后通過試驗對最優的方案進行驗證。試材生菜種植在在氣霧栽培箱內部擠塑板上的定植孔內，在擠塑板上的生菜冠層布置風速及溫度監測點，每個監測點距離擠塑板水平距離約5 cm，每塊擠塑板上均勻布置3個監測點，從上往下間隔依次為40 cm，左右等距，6個面總共布置18個測量風速值及溫度值監測點，如圖2所示。使用長為8 cm，寬為8 cm的軸流風機調控風速，可布置在12個進回風口處，如圖2所示，風速可控范圍在0～3.8 m/s。溫度采用溫度傳感器測量（深圳燦輝電子有限公司生產），布置在溫度監測點處，溫度測量范圍為-50～110℃，測試精度：±0.5 ℃，采樣周期2 s。風速值測量采用手持數字風速儀進行測量，在風速監測點處測量，分辨率為0.001 m/s，風速測量誤差±3%。

注：l表示18個風速和溫度值監測點。￡表示軸流風機安裝位置。

2 CFD模型

2.1 數值模型

2.2 生菜作物冠層質熱交換模型

氣霧栽培箱內部生長的生菜把截獲的光照輻射以顯熱和潛熱的方式與周圍空氣進行質熱交換，根據作物冠層的能量平衡原理建立數學模型[14]：

2.3 生菜作物多孔介質模型

2.4 定義材料和邊界條件

對氣霧栽培箱模型進行簡化，該型氣霧栽培箱為封閉式，因在智能氣候室內對真實環境進行模擬設置，所以不存在太陽光直射，光照來源于智能氣候室內的植物生長燈，所以忽略太陽輻射，將氣霧栽培箱6個亞克力板壁面簡化為溫度面（wall-temperature）；將氣霧栽培箱XPS栽培板設置為絕熱壁面（wall）；進風口設置為速度入口（velocity-inlet）；出風口設置為壓力出口，壓力設置為0；進出風溫度設置為實測值18 ℃[18-20]。模擬的主要參數和取值如表1所示。

表1 模擬參數設置和取值

2.5 計算域和網格劃分

對氣霧栽培箱模型進行簡化，建模比例為1∶1，選取氣霧栽培箱葉片生長區作為計算域。在Ansys Workben-ch的ICEM模塊中進行網格劃分，由于六面體網格節點易于計算，結構化網格生成的網格質量高，采用結構化六面體方式對氣霧栽培箱進行網格劃分，利用Fluent軟件對網格文件進行后處理。最后計算出網格節點數為531 628，網格數量為3 070 692，網格劃分質量良好。圖3為蔬菜氣霧栽培箱網格劃分圖。

圖3 氣霧栽培箱網格劃分

3 CFD模擬和氣流循環優化設計

3.1 模擬優化方案設計

為探究氣霧栽培箱內部生菜葉片生長區域的風速及溫度分布，選取氣霧栽培箱內部生菜冠層平面進行分析，通過相關文獻[10-12]可知生菜作物冠層平面適宜風速處于0.3～1.30 m/s。

為了使氣霧栽培箱氣流和溫度滿足生菜作物生長需要，根據鄧書輝等[21]、姚家君等[22]、胥芳等[23]、呂曉蘭等[24]學者對不同通風循環模式下流場和溫度場的研究，本文設計了3種氣流循環方式，即方案1，頂面進側面出；方案2，側面進頂面出；方案3，側面進側面出。如圖4所示。利用CFD模擬并比較3種方案中生菜生長冠層平面風速及溫度分布規律，得出最優氣流循環方案。

注：紅色平面表示截面A、B和C。

方案1、方案2、方案3三維模型建立都相同，網格劃分也都一致。將回風口設置為壓力出口邊界條件，將進風口設置為速度入口邊界條件。為保證生菜生長區域在所需的適宜風速范圍內，對不同方案的通風量，入口風速依據下列公式進行計算[25]：

3.2 方案模擬仿真分析

3.2.1 頂面進側面出

方案1是頂面3個送風口進風，3個側面9個回風口出風。對溫度分布圖5a和5b分析可知，從送風口處向下呈兩邊回流的趨勢分布。回流溫度的分布與氣流的分布存在一定耦合關系，風速值較大區域溫度值較小，風速值較小區域溫度值較大。高溫區域主要集中在進口截面的兩側，這是因為空氣從進口進去后會在兩側形成兩個渦，與主流區域的能量交換較少，所以溫度相對較高。作物冠層平面最大溫度為24.55 ℃；溫度平均值為19.05 ℃；適宜溫度所占百分比為93.6%。

對進口截面氣流速度分布圖5c和5d分析可知，氣流經由氣霧栽培箱頂面3個送風口進入箱體后，主流區域直接沖向底部，沖擊底部后再向上流動。送風口下方氣流變化梯度較大，氣流方向垂直向下，并且主要集中在送風口的下部，然后沿著豎直方向向下移動，在氣霧栽培箱內形成漩渦流。部分氣流從氣霧栽培箱側面回風口流出，導致氣霧栽培箱內部風速減小。對速度分布流線圖5e分析可知，速度從進口進入后，主流區域會在觸底后向上流動，并從兩側的出口流出，整體呈現出對稱流動，3個出口中靠近底部區域的出口流量比較大，而靠近上側的出口流量比較小。作物冠層平面最大風速為1.43 m/s；內部風速平均值為0.19 m/s；適宜風速所占百分比為58.1%。

3.2.2 側面進頂面出

方案2是3個側面9個送風口進風，頂面3個回風口出風。對溫度分布圖6a和6b分析可知，送風口處溫度較低，回風口處溫度稍高，底部區域的溫度略高于中間區域，這是因為從3個側面進風以后，空氣主要向上流動，造成了底部區域的溫度略高于中間區域，風速小的地方溫度較高，風速大的地方溫度較低，但整體內部溫度分布較為均勻，溫度與氣流的分布存在一定的耦合關系。作物冠層平面最大溫度為21.68 ℃；最小溫度為17.85 ℃；溫度平均值為18.38 ℃；適宜溫度所占百分比為99.98%。

圖5 方案1溫度和速度分布云圖及速度分布流線圖

對氣流速度分布圖6c和6d分析可知，氣流經由氣霧栽培箱3個側面9個送風口進入箱體后，進口速度較高，沖擊內壁面后向兩側流動，速度分布相對比較均勻。從進口截面上的溫度分布云圖可以看出，溫度分布比較均勻，沒有明顯的高溫與低溫區域。兩個進口面中間的區域溫度略高于進口面進風口上下方向。氣流方向朝遠離送風口的方向速度逐漸變小。氣流從氣霧栽培箱頂面回風口流出時氣流速度有變大的趨勢。對速度分布流線圖6e分析可知，流動整體分布比較均勻，沒有明顯的高速與低速。底部區域的流動比較小，這是因為出口在上側，流體進去后，主要向上流動。作物冠層平面最大風速為1.62 m/s；內部風速平均值為0.19 m/s；適宜風速所占百分比為59.6%。

圖6 方案2溫度和速度分布云圖及速度分布流線圖

3.2.3 側面進側面出

方案3是1個側面3個送風口進風，2個側面6個回風口出風。對溫度分布圖7a和7b分析可知，送風口處溫度向前逐漸變高，兩邊逐漸降低。回風口處溫度稍高，靠近進口端面的溫度比較低，出口端面的溫度比較高。靠近出口截面上的溫度比較高。溫度整體呈現出明顯的左右不均勻性。作物冠層平面最大溫度為25.95 ℃；最小溫度為17.65 ℃；溫度平均值為19.2 ℃；適宜溫度所占百分比為92.4%。

對氣流速度分布圖7c和7d分析可知，進風口上下方向氣流變化梯度較大，進口面上的速度比較高，其余5個面域的速度比較低，速度分布的均勻性比較差。送風口出氣流速度較大，氣流方向朝遠離送風口的方向速度向四周逐漸變小，3個送風口的氣流分布上下不均勻，存在回風口的2個側面氣流分布梯度較小，氣流從回風口流出時氣流速度有變大的趨勢，靠近出口的截面上整體速度比較低。對速度分布流線圖7e分析可知，進口段流線分布比較密集，該段的流速較高，進去后從兩側出口流出，進口端正對面的區域中，流動情況比較差。作物冠層平面最大風速為1.41 m/s；內部風速平均值為0.08 m/s；適宜風速所占百分比為54.3%。

圖7 方案3溫度和速度分布云圖及速度分布流線圖

3.2.4 方案優化討論

對3個方案溫度和風速分布云圖分析及表2、3可知。風速送風口處風速一般比較大，溫度較低，所以溫度的分布與氣流的分布存在一定耦合關系，風速較大區域溫度較小。3個方案風速最小值都為0，方案1風速最大值為1.43 m/s，風速平均值為0.19 m/s。風速位于生菜生長適宜風速值0.3～1.3 m/s的區域占58.1%，溫度平均值為19.05 ℃，溫度位于生菜生長適宜溫度值15～20 ℃的區域占93.6%。方案2風速最大值為1.62 m/s，風速平均值為0.19 m/s。風速位于生菜生長適宜風速值0.3～1.3 m/s的區域占59.6%，溫度平均值為18.38 ℃，溫度位于生菜生長適宜溫度值15～20 ℃的區域占99.98%。方案3風速最大值為1.41 m/s，風速平均值為0.08 m/s。風速位于生菜生長適宜風速值0.3~1.3 m/s的區域占54.3%，溫度平均值為19.2℃，溫度位于生菜生長適宜溫度值15~20℃的區域占92.4%。綜合比較3種方案，方案2為最優方案。

表2 不同氣流循環方案作物冠層平面風速

表3 不同氣流循環方案作物冠層平面溫度

3.2.5 方案通風死角討論

氣霧栽培箱內部進行農作物栽培時，由于部分區域通風不均勻會產生死角，導致沒有氣流速度分布，影響溫度的分布均勻性，會影響農作物的正常生長。所以對進出風口的設計布局顯得尤為重要，包括進出風口大小、數量、布局等都會影響氣流的均勻流動。設計合理的進出風口可避免農作物生長過程中引起的一系列生產問題。對文中3種方案的死角問題進行討論，通過數值模型截面分析可知（如表4和圖8所示）。

注：紅框表示通風死角區域

方案1中，通風死角區域占比0.844%，主要分布在通風進口面和通風出口面之間的交界處靠近上側的部分，這是因為風從上側進入后，到達底部區域后方向旋轉，向上流動，從出口面的3個出口流出，越向上風流量越小，所以在上側區域會出現部分死角。方案2中，通風死角區域占比僅0.069%，沒有明顯的通風死角，主要原因是進出風口布局較為合理。方案3中，通風死角區域較大，占16.7%，主要分布在2個出口所在面的中間區域，這是因為風從進口進去后，直接從兩側流向出口，只有極少空氣流向兩個出口的中間區域。通過對3中方案進行對比可知，方案2的通風性較好，通風死角區域占比最小。

表4 不同方案通風死角所占百分比及溫度平均值

3.3 試驗結果分析與CFD模型驗證

為驗證仿真模型的可靠性和正確性，在蔬菜氣霧栽培箱內種植生菜的情況下，對最優方案2各點溫度和風速的實測值與模擬值進行試驗驗證，如圖9所示。對圖9a的18個溫度監測點的模擬值與實測值進行對比分析可知，測點11的誤差值最大，最大誤差為1.35 ℃，相對誤差最大6.9%，平均相對誤差為3.9%，均方根誤差為0.86 ℃。對圖9b的18個風速監測點的模擬值與實測值進行對比分析可知，測點6的絕對誤差最大，最大誤差為0.54 m/s，相對誤差最大9.5%，平均相對誤差為3.5%，均方根誤差為0.26 m/s。

從圖中對數據比較可以發現，實測值與模擬值存在偏差，但整體的溫度和風速分布較為均勻，所建立的CFD模型是準確的，可用此模型對其他方案進行驗證分析。

4 結 論

1）本文中設計了3種氣流循環方式：方案1為頂面進側面、方案2為側面進頂面出、方案3為側面進側面出。對3個方案溫度和風速分布云圖分析后得到風速送風口處風速一般比較大，溫度較低，風速較大區域溫度較小，溫度的分布和氣流的分布存在一定耦合關系。方案1風速位于生菜生長適宜風速值區域占58.1%，適宜溫度值區域占93.6%，通風死角區域占比0.844%。方案2風速位于生菜生長適宜風速值區域占59.6%，適宜溫度值區域占99.98%，通風死角區域占比0.069%。方案3風速位于生菜生長適宜風速值區域占54.3%,，適宜溫度值區域占92.4%，通風死角區域占比16.7%。綜合比較3種方案并結合生菜生長的溫度和風速，方案2為最優方案。

2）將18個點的模擬值與實測值進行對比，溫度平均相對誤差為3.9%，均方根誤差為0.86 ℃。風速平均相對誤差為3.5%，均方根誤差為0.26 m/s。模擬值和實測值誤差較小，模擬效果良好。建立的CFD模型可準確的描述蔬菜氣霧栽培箱內部的流場和溫度場變化規律和空間分布特性。

圖9 氣霧栽培箱溫度及風速監測點模擬值與實測值比較

3）本文中模擬了3種氣流循環方案，得到側面進頂面為最佳氣流循環模式，并對送回風口的不同位置布局進行了研究。如果繼續進行優化研究，可結合正交試驗方法，將送回風口大小、風速、數量及位置等設置為不同因素，對該類蔬菜氣霧栽培箱的送回風口的通風進行優化設計，避免通風死角區域，模擬并設計出最佳通風模式。

4）所建立的CFD模型，模擬預測了不同氣流循環方式下蔬菜氣霧栽培箱內部溫度及風速的變化規律。本文對單個的蔬菜氣霧栽培箱流場和溫度場進行了模擬，可對大規模商業化種植，包括在溫室大棚內氣霧栽培設備布局、大型霧化工廠種植、家庭陽臺霧化裝置等，氣霧栽培技術的環境建模和控制提供更多理論依據和信息，彌補物理試驗的不足。

[1] 張建國，何春梅，凌敏，等. 氣霧栽培的應用與研究綜述[J]. 林業與環境科學，2017，33(4)：130－134.

Zhang Jianguo, He Chunmei, Ling Min, et al. Review on the application and research progress of aeroponic cultivation[J]. Forestry and Environmental Science, 2017, 33(4): 130－134. (in Chinese with English abstract)

[2] 高建民，黃桂珍，尹文楚，等. 桁架式超聲霧化栽培器的霧滴沉降和根際溫濕度變化規律[J]. 農業工程學報，2013，29(6): 185－192. Gao Jianmin, Huang Guizhen, Yin Wenchu, et al. Variation of root-zone droplet deposition, temperatureand humidity of truss ultrasonic atomization cultivator[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(6): 185－192. (in Chinese with English abstract)

[3] 劉義飛. 基于LabVIEW的溫室番茄霧培控制系統研究[D]. 北京：中國農業科學院，2014. Liu Yifei. Aeroponics Control System for Tomato Cultivation in Greenhouse Based on LabVIEW[D]. Beijing: Chinese Academy of Agricultural Sciences, 2014. (in Chinese with English abstract)

[4] 劉巖. 植物工廠中霧培成套裝置的設計與應用[D]. 杭州：浙江大學，2014. Liu Yan. Design and Application of Aeroponics Package Unit in Plant Factory[D]. Hangzhou: Zhejiang University, 2014. (in Chinese with English abstract)

[5] 王新忠，張偉建，張良，等. 基于CFD的夏季屋頂全開型玻璃溫室自然通風流場分析[J]. 農業機械學報，2016，47(10): 332－337. Wang Xinzhong, Zhang Weijian, Zhang Liang, et al. 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. (in Chinese with English abstract)

[6] 翁衛兵，房殿軍，李強，等. 冷藏運輸廂內流場和溫度場協同控制[J].農業機械學報，2014, 45(6): 260－265. Weng Weibing, Fang Dianjun, Li Qiang, et al. Cooperative control of flow field and temperature field in refrigerated transport carriage[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(6): 260－265. (in Chinese with English abstract)

[7] 林加勇，劉繼軍，孟慶利，等. 公豬舍夏季溫度和流場數值CFD模擬及驗證[J].農業工程學報，2016，32(23): 207－212. Lin Jiayong, Liu Jijun, Meng Qingli, et al. Numerical CFD simulation and verification of summer indoor temperature and airflow field in boar building[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(23): 207－212. (in Chinese with English abstract)

[8] 高云，陳震撼，王瑜，等. 多環境參數控制的豬養殖箱設計及箱內氣流場分析[J]. 農業工程學報，2019，35(2)：203－212. Gao Yun, Chen Zhenhan, Wang Yu, et al. Design for pig breeding chamber under multiple environment variable control and analysis of internal flow field[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(2): 203－212. (in Chinese with English abstract)

[9] 張成波，楊其長. 植物工廠概論[M]．北京：中國農業科學技術出版社，2005：56－156.

[10] Kitaya Y, Tsuruyama J, Shibuya T, et al. Effects of air current speed on gas exchange in plant leaves and plant canopies[J]. Advances in Space Research the Official Journal of the Committee on Space Research, 2003, 31(1): 177－182.

[11] Thongbai P, Kozai T, Ohyama K. CO2and air circulation effects on photosynthesis and transpiration of tomato seedlings[J].Scientia Horticulturae, 2010, 126(3): 338－344.

[12] Kitaya Y. Effects of air velocity on photosynthesis of plant canopies under elevated CO2levels in a plant culture system[J]. Advances in Space Research the Official Journal of the Committee on Space Research, 2004, 34(7): 1466－1469.

[13] Versteeg H K, Malalasekera W. An Introduction to Computational Fluid Dynamics[M]. London: Longman Group Ltd. , 1995.

[14] Boulard T, Wang S. Experimental and numerical studies on the heterogeneity of crop transpiration in a plastic tunnel[J]. Computers and Electronics in Agriculture, 2002, 34(1/2/3): 173－190.

[15] Bruse M. Development of a Numerical Model for the Simulation of Exchange Processes Between Small Scale Environmental Design and Microclimate in Urban Areas[D]. Bochum: University of Bochum, 1998.

[16] Bruse M．Development of a Micro-scale Model for the Calculation of Surface Temperature in Structured Terrain[D]. Bochum: University of Bochum，1995．

[17] Kichah A, Bournet P E, Migeon C, et al. Measurement and CFD simulation of microclimate characteristics and transpiration of an impatiens pot plant crop in a greenhouse[J]. Biosystems Engineering, 2012, 112(3): 252 －252.

[18] Wilson J D. Numerical studies of flow through a windbreak[J]. Journal of Wind Engineering & Industrial Aerodynamics, 1985, 21(2): 119－154.

[19] 王福軍. 計算流體動力學分析[M]．北京：清華大學出版社，2004: 1－57.

[20] 張師帥. 計算流體動力學及其應用[M]．武漢：華中科技大學出版社，2011：1－110.

[21] 鄧書輝，施正香，李保明，等. 擋風板對低屋面橫向通風牛舍內空氣流場影響的PIV測試[J]. 農業工程學報，2019，35(1)：188－194. Deng Shuhui, Shi Zhengxiang, Li Baoming, et al. Influence of baffle on air flow field in low profile cross ventilated cattle barn analysis based on PIV test[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(1): 188－194. (in Chinese with English abstract)

[22] 姚家君，郭彬彬，丁為民，等. 基于鵝舍氣流場CFD模擬的通風系統結構優化與驗證[J]. 農業工程學報，2017，33(3)：214－220. Yao Jiajun, Guo Binbin, Ding Weimin, et al. Structure optimization and validation of goose house ventilation system based on airflow field simulation by CFD[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(3): 214－ 220. (in Chinese with English abstract)

[23] 胥芳，蔡彥文，陳教料，等. 濕簾-風機降溫下的溫室熱/流場模擬及降溫系統參數優化[J]. 農業工程學報，2015，31(9)：201－208. Xu Fang, Cai Yanwen, Chen Jiaoliao, et al. Temperature/ flow field simulation and parameter optimal design for greenhouses with fan-pad evaporative cooling system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(9): 201－208. (in Chinese with English abstract)

[24] 呂曉蘭，張美娜，常有宏，等. 果園風送噴霧機導流板角度對氣流場三維分布的影響[J]. 農業工程學報，2017，33(15)：81－87. Lü Xiaolan, Zhang Meina, Chang Youhong, et al. Influence of deflector angles for orchard air-assisted sprayer on 3D airflow distribution[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(15): 81－87. (in Chinese with English abstract)

[25] 蔡增基，龍天渝. 流體力學泵與風機[M]．北京：中國建筑工業出版社，2009：156－179.

CFD simulation of flow field and temperature field of different air flow circulation modes in vegetable aeroponics cultivation box

Jiang Pin1,2, Yang Xiwen1, Luo Yahui1,2, Long Lixia1, Shi Yixin1,2, Hu Wenwu1,2

(1.410128,; 2.410128,)

The air flow and temperature in the vegetable aeroponics cultivation box play an important role in the internal environment regulation and crop growth. The uniform air flow around crops can promote the growth rate of the crop. In order to study the distribution of flow field and temperature field in vegetable aeroponics cultivation box, CFD models of vegetable aeroponics cultivation box under different air flow circulation modes were established based on computational fluid dynamic (CFD) method, FLUENT software, standard turbulence model, porous medium model, crop canopy heat exchange model. The environment in the aeroponics cultivation box was optimized, and three air flow circulation schemes were designed(air flows from top to side, side to top and side to side). The layout of the different positions of the return air outlets was studied, and the numerical simulations of the three airflow circulation schemes were carried out. The simulation results showed that in the scheme of air flow entering from the top and exiting from the side, the suitable wind speed value area of lettuce accounts for 58.1%, the suitable temperature value area accounts for 93.6%, and the ventilation dead zone area accounts for 0.844%.In the scheme of air flow side entering and exiting the top surface, the suitable wind speed value area of lettuce accounts for 59.6%, the suitable temperature value area accounts for 99.98%, and the ventilation dead zone area accounts for 0.069%. In the scheme of air flow side entering side exit, the suitable wind speed value area for lettuce growth accounts for 54.3%, the suitable temperature value area accounts for 92.4%, and the ventilation dead zone area accounts for 16.7%. After comparison, the side entering the top surface was obtained as the optimal air circulation scheme. The temperature and wind speed simulation values and measured values in the aeroponics cultivation box were compared. The results showed that the average temperature relative error was 3.9% and the root mean square error was 0.86 ℃. The average relative error of wind speed was 3.5%, and the root mean square error was 0.26 m/s. The error between simulated and measured values was small, and the simulation effect was good. The accuracy of the CFD model was verified. In crop cultivation, aeroponics was a new cultivation method in the field of soilless cultivation, which could be applied to greenhouses, plant factories, family balconies and other places through microclimate control technology. In this paper, by regulating the wind speed in the aeroponics box, the climate growth environment and temperature uniformity of the crops can be ensured, and the problems of the roots and canopies of the crops can be effectively solved, and the crop yield and quality can be improved. Airflow velocity and temperature are important factors for crop growth in aeroponics boxes, which affect the transpiration and photosynthesis of crops. Therefore, it is of great significance to study the changes of air velocity and temperature to promote air aeroponics technology. This study provides a reference for the flow field and temperature change law in vegetable aeroponics cultivation box, internal environment adjustment and device optimization design.

wind speed; temperature; computational fluid dynamics; optimization design; vegetable aeroponics cultivation box

2019-02-20

2019-07-16

國家重點研發計劃項目（2017YFD0700903-2）；湖南省重點研發計劃項目（2018NK2063）

蔣 蘋，教授，博士，博士生導師，主要從事農業機械智能控制研究。Email：teacher_jp@163.com

10.11975/j.issn.1002-6819.2019.16.026

S317

A

1002-6819(2019)-16-0233-09

蔣 蘋，楊希文，羅亞輝，龍莉霞，石毅新，胡文武.蔬菜氣霧栽培箱不同氣流循環方式的流場和溫度場CFD模擬[J]. 農業工程學報，2019，35(16)：233－241. doi：10.11975/j.issn.1002-6819.2019.16.026 http://www.tcsae.org

Jiang Pin, Yang Xiwen, Luo Yahui, Long Lixia, Shi Yixin, Hu Wenwu. CFD simulation of flow field and temperature field of different air flow circulation modes in vegetable aeroponics cultivation box[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(16): 233－241. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.16.026 http://www.tcsae.org