農業溫室系統的多相流傳遞過程研究進展

張建文，賀亞如，王慶貴

(1.北京化工大學 流體力學與傳熱研究室，北京 100029；2.西藏山南雅礱綠色畜禽有限公司，西藏 山南 856100)

?

農業溫室系統的多相流傳遞過程研究進展

張建文1，賀亞如1，王慶貴2

(1.北京化工大學 流體力學與傳熱研究室，北京 100029；2.西藏山南雅礱綠色畜禽有限公司，西藏 山南 856100)

摘要：分析了研究溫室多相傳遞過程的特點，對溫室內復雜的多相流動與傳熱問題進行了闡述，進而對目前關于溫室傳熱、傳質過程的CFD(計算流體力學)研究進行了綜述，總結了CFD在溫室多相流動與傳熱傳質研究中的優勢，并對CFD研究溫室多相流系統的前景進行了展望。

關鍵詞：精細農業；溫室系統；流體流動；傳熱傳質；CFD

0引言

多相流本身是個復雜的流動過程，而溫室中的多相流過程又同時涉及傳熱傳質、植物作用、相間交換等問題。近年來，溫室農業正逐步向精細化、集約化、科技化發展。因而研究溫室中的多相流過程對于指導溫室農業生產具有非常重要的意義。

精細農業是一種現代化的農業管理理念，旨在實現定點的精確化，定量的精確化，定時的精確化[1]?，F代溫室農業生產適時引入精細農業思想，以期通過實時獲取農田中每個小區域的作物信息，診斷作物長勢和產量在空間上的差異，區別對待，按需實施定位調控，以實現作物優質高產[2-5]。

溫室環境精細化控制的首要任務是實現溫室環境因子的定點監測。溫室環境定點監測主要是通過在溫室內布控傳感器來檢測溫室內環境因子，但是由于溫室是一個復雜的多相傳遞系統，溫室內環境因子非常復雜，完全依賴傳感器檢測實現起來較為困難，且將使成本大大增加，因此完全依靠傳感器實現溫室環境控制具有一定的難度[6]。

近年來，越來越多的學者開始研究溫室多相流系統，以實現對溫室環境的控制。溫室內定點的作物信息主要包括水分、養分、光照、溫度、濕度等，對這些環境因素的控制構成了溫室精細控制的基礎，因此，研究溫室內的多相流傳遞過程具有重要的意義。

1溫室中的多相流問題

1.1溫室中的多相流特點

溫室農業生產系統中，包括溫室空間中的氣體環境、水環境、土壤以及植物冠層等微環境；期間涉及到水分、養分在土壤、植物、大氣間的傳遞過程(圖1)，構成了溫室農業生產過程中的復雜多相流傳遞問題。

溫室內的多相流傳遞特點主要表現為：(1)傳遞主體復雜。以水分傳遞為例，包括植物與土壤間的水分傳遞，植物與空氣間的水分傳遞，土壤與空氣間的水分傳遞，以及由灌溉引起的溫室內與外界環境間的水分傳遞。(2)傳遞涉及多種作用：植物的生理作用，如蒸騰作用、呼吸作用、光合作用等。(3)傳遞過程涉及多種溫室調控設施，如灌溉設施、通風設施、加溫降溫設施及防蟲網等，這些調控設施的存在使得溫室多相流問題呈現出一定的特殊性。(4)傳熱過程兼具對流傳熱與輻射傳熱，包括空氣對流傳熱、太陽輻射及太陽輻射導致的二次(局部)對流傳熱等。(5)傳質傳熱介質包括植物，構成了特殊的傳質傳熱阻力。

圖1　溫室中的多相流傳遞

溫室生產過程包括了諸多同時發生的過程，并且這些過程具有不同的反應時間和反應模式，而且存在多種相互作用。概括起來，溫室中的多相流傳遞主要包括多相流傳熱及傳質過程，其中多相流組分傳遞又包括空氣傳遞和養分傳遞。

1.2溫室中的多相流傳質過程

1.2.1溫室空氣環境的傳遞過程溫室中的空氣環境傳遞主要指溫室內外的空氣交換。一般溫室生產過程中，為了保證溫室的保溫效果，通常需要減少通風換氣率。但在生產過程中又有通風換氣的需求[7]：(1)必要的通風換氣可防止溫室晝夜溫度逆差，促進作物正常生長；(2)密閉溫室容易導致溫室內CO2缺乏，則植物光合作用不能順利進行，通風能夠補充CO2；(3)排除有毒(或有害)氣體，并調節溫室內空氣的相對濕度；(4)夜間通風也可減少溫室內作物的呼吸作用，減少同化物消耗，提高作物產量；(5)夏季通風能夠降低溫室內溫度。

由此可見，溫室內外的空氣交換十分必要，而溫室通風系統的設計又直接影響著溫室內的空氣流動，故溫室通風系統的設計十分關鍵。

1.2.2溫室中的養分傳遞溫室中的養分傳遞包括養分在土壤、植物、空氣間的傳遞與流動，其中不僅涉及到相變過程，還涉及到溫室與外界環境在溫室調控設施作用下的養分交換。以溫室中的水分傳遞為例，水分傳遞的同時又伴隨著水的相變過程：土壤中的水分通過蒸發作用散失到空氣中；夜間溫室內水蒸氣的冷凝過程；植物通過蒸騰作用從土壤表面吸收水分，又通過葉面上的氣孔散發到空氣中；灌溉設施引起的溫室內外的水分傳遞；另外溫室的通風系統也促進溫室內外的水汽傳遞。

1.3溫室中的多相流傳熱

根據換熱情形的差異，多相流傳熱可分為與外界無加熱或冷卻等熱量交換過程的絕熱多相流和有熱量交換的多相流傳熱過程。在有熱交換的多相流流動過程中通常伴隨著單組分的相變過程，如液體汽化成蒸汽或蒸汽凝結成液體。

由于太陽的輻射作用，及溫室中配置的加熱降溫設備，溫室生產中常常伴隨著熱交換及輻射傳輸過程。

由通風引起的能量交換(通風引起溫室內外能量交換)；加熱系統與室內空氣之間的熱交換(冬季為了維持溫室內溫度，需通過加熱)；土壤與周圍環境的熱交換(土壤與溫室內空氣進行對流換熱，與植物進行輻射換熱)；植物冠層與空氣間的能量交換。

由以上對溫室系統多相流問題的分析可知，實現溫室內環境的精細控制具有一定的難度。

2CFD對溫室農業生產多相傳遞過程的控制

常規意義上溫室的精細控制主要是利用遙感技術獲取溫室內作物生長環境、生長狀況和空間變異的大量時空變化信息，并針對性地提供處方信息。由于溫室內多相流問題十分復雜，溫室環境是個多變量耦合、時變、非線性的復雜系統，需要在溫室內大量布設傳感器，實施時具有一定難度且成本較高。

CFD研究能夠對各種結構、環境參數下的溫室內的各項參數分布進行一個準確直觀的預測，與傳感器相比能夠實現對溫室內各個點的環境信息的預測，更加精細、詳盡，且能夠大大降低成本。另外，CFD研究建立在數值計算的基礎上，能夠指導人們從機理出發，實現對溫室的精細調控，從而指導溫室結構的選擇及環境參數的調控。因此，CFD在國內外溫室研究中的運用越來越普遍[8]。

2.1溫室傳質過程的CFD研究

2.1.1溫室通風(空氣傳遞)控制的CFD研究良好的通風系統對于保持溫室內氣候的穩定性分布至關重要。溫室內的通風情況嚴重影響著溫室內空氣的交換及溫室的內部氣候條件，而這些氣候條件又關系到溫室內作物的生長[9]。CFD研究能夠模擬不同通風裝置下的溫室內部的氣候分布，減少實驗次數，降低成本，并且能夠針對自然通風口的配置提出指導意見。

Mistriotis等[10]研究了四連棟地中海型溫室。他們先后考慮了3種情形：(1)2個側窗(開度為45°)以及4個水平開的頂窗；(2)開度為45°的側窗，第3、4個頂窗開度為85°；(3)只開頂窗，開度為85°。由于頂窗關于中間的屋脊對稱，故模擬時只研究2個連棟(圖2)。在第一種情形下，新鮮空氣從側窗進入，熱空氣主要從中間的頂窗R2處流出。通風口R2扮演的角色表明更大的頂窗開度將會提高溫室內的通風速率(第二種情形，頂窗開度85°)。這種情形下，第二股冷空氣流從頂窗R1較低處流入，R1下半部分相當于一個入口，但上半部分是暖空氣出口。當側窗關閉后(第三種情形)，R1是冷流入口，新鮮空氣徑直流向地面，將熱空氣推向溫室中心，熱空氣主要從R2流出，只有一小部分從R1流出。不難看出，第三種情形下的通風效率不如前兩種情形。該研究證實了開通側窗的必要性，即側窗主要是冷空氣入口，而頂窗主要是熱空氣出口。

圖2　Mistriotis的通風口配置研究

沈明衛等[11]為了比較華北地區連棟塑料溫室3種典型的不同頂窗配置形式對其自然通風效果的影響，運用CFD對配置不同頂窗的溫室內自然通風氣流場進行了三維穩態模擬。通過對不同頂窗配置的溫室內氣流場分布情況的比較，對于溫室頂窗配置的選擇具有一定的指導意義。

郭萬東等[12]利用CFD軟件CFX對華南型單棟塑料溫室的表面風壓進行了數值模擬。將模擬結果與前人所做的風洞實驗結果比較后發現，數值模擬和風洞實驗的風壓分布趨勢相近，表明CFD方法能夠定性和定量地預測溫室表面風壓的分布規律。

Francisco Domingo Molina-Aiz等[13]用CFD數值模擬方法對Almeria型溫室內的自然通風的風速效應進行了分析。用多孔介質模擬防蟲網和植物間隙中的空氣流動。結果表明，實驗值與模擬值的吻合度較高。對只配置側窗和同時配置有頂窗的Almeria型溫室內的溫度及空氣速率分布進行了預測，證明了配置頂部通風口的Almeria型溫室內通風的必要性，對Almeria型溫室通風口的設計具有重要的指導意義。

2.1.2溫室濕度(水分)控制的CFD研究濕度控制是冬季溫室最重要的問題之一；過高的相對濕度(RH)和植物表面自由水的存在很容易引起真菌性病害的發展。

Davide Piscia等[14]對四連棟溫室內夜間條件下的冷凝過程進行了數值分析。通過UDF導入了薄膜冷凝模型。用CFD模擬了頂棚溫度對溫室濕度的影響，并用實驗進行了比對。模擬結果與實驗數據的比對驗證了CFD模型在預測溫室氣候變量方面的準確性。

模擬表明，夜間條件下，頂棚是溫室內溫度最低的表面；頂棚溫度和溫室濕度比(RH)存在著很強的關聯。這就意味著通過控制頂棚溫度來控制溫室內的濕度及冷凝是有可能的。這一CFD冷凝模型可用于溫室濕度控制的設計策略中，尤其是無供暖設備的溫室中。

譚勝男等[15]采用CFD方法對Venlo型玻璃溫室在夏天自然通風狀況下的噴霧降溫過程進行了數值模擬。在充分考慮太陽輻射影響和室內水蒸氣傳輸過程基礎上，結合離散相模型，構建了求解室內系統的3D數學模型，并對邊界條件的設置進行了探討。對Venlo型溫室在噴霧降溫措施下室內空氣溫度、濕度的變化以及在整個溫室空間的溫度分布進行了數值模擬與預測。結果顯示溫室內測點溫、濕度的模擬值與實測值的平均相對誤差分別為4.9%和5.7%，模擬結果與測試結果吻合良好。

Bartzanas等[16]在未使用輻射模型的情況下用UDF來描述植物水和熱交換。用包含有氣孔阻力的方程來模擬植物的顯熱和潛熱交換。模擬值與試驗值間的誤差介于12%～15%。但是Wang等認為影響蒸騰作用的主要因素是太陽輻射。另外一些研究者Sapounas等[17]用帶有UDF的植物多孔介質模型來模擬作物蒸散。

程秀花等[18]對栽有番茄作物的溫室內部的溫濕度進行了數值模擬，將溫室內水分蒸發與作物蒸騰以體積熱源的方式添加，對比實驗結果發現此模擬的模型處理具有有效性。

Abderzak Kichah等[19]通過實驗驗證了計算流體力學(CFD)在模擬溫室內的溫濕度及植物蒸騰速率等方面的準確性。在100 m2的溫室隔間內，對長凳上的盆栽鳳仙花屬植物進行了實驗觀察。同時對植物及其周圍環境間的熱能及水汽交換進行了二維模擬。植物子模型將植物作為多孔介質并且針對潛熱和顯熱傳遞加入了特定的源相。通過實驗與模擬結果的比較發現，CFD模型不僅能夠模擬植物冠層內的溫度分布，同時還能模擬一天內植物葉面溫度的垂直分布。實驗和模擬的比較表明，CFD模擬在預測從頂部到底部的植物的葉面溫度上有很好的準確性，同時對一天內植物蒸騰速率的變化趨勢提供了一個好的預估(r>0.92)。24 h內對整體蒸騰速率的測量值和模擬預測值間的差值<15%。最終證明該模型可通過預測植物水需求，及時調整灌溉策略。

2.2溫室傳熱過程的CFD研究

溫室內微氣候并不只受對流傳遞的影響，輻射也是很重要的一個影響。但是，在早期的很多應用CFD研究中，都沒有同時解決輻射傳遞和對流傳遞的耦合問題。

Fodaros等[20]對日光照射下弧形隧道溫室的氣候進行了數值模擬，結果表明輻射對溫室系統的影響是呈動態變化的?？梢姕厥覂鹊妮椛溥^程相當復雜，故在CFD模擬研究時通常用簡化的方式來處理[21]。對太陽及熱輻射效應的處理，多數情況下都是在物理邊界上設定具體的壁面溫度[22-23]或熱通量[24-25]來代替直接在計算域的邊界上設置輻射條件。Kim等[26]建立了一個用于計算單棟溫室內表面太陽輻射通量的簡單模型。截至目前，僅有少數研究包含有求解光強的輻射傳輸方程(RTE)的輻射機理。

雖然Lee等[27]用P1模型來求解RTE。但是，到目前為止，應用最廣泛的仍是DO輻射模型。由于DO輻射模型兼具準確性、計算經濟性、靈活性的優點，使得模擬半透明壁面輻射成為可能。

Chen等[28]基于低雷諾數k-ε湍流動能原理，應用DO輻射模型來模擬溫室內輻射熱傳遞，用多孔介質來模擬植物并考慮植物的顯熱和潛熱傳導，建立了一個加熱溫室的CFD模型(圖3)。并開發了一種新穎的基于CFD和EPM的溫室熱控制策略，這種熱控制策略能夠節省能量消耗量，并能維持適當的室內溫度和更好的控制系統性能。最后用實驗證明了以CFD-EPM為基礎的溫室熱控制策略能夠更準確地預測溫度(圖4)，從而最大限度地降低溫室內的能量消耗，提高溫室系統性能。

圖3　溫室供暖系統的熱平衡

圖4　溫室內空氣溫度的模擬值與實測值的對比圖

譚勝男[29]以Venlo型現代化溫室為研究對象，利用CFD數值模擬軟件和現場試驗，對溫室內微環境進行了模擬、驗證和預測。模擬過程中選擇DO輻射模型，用多孔介質模型模擬溫室內種植的黃瓜作物。將模擬結果與試驗測試結果進行對比驗證分析，發現結果吻合良好，證明了所建立的CFD模型的有效可行性。

Wang等[30]基于中國中部典型的扇形垂直截面的塑料大棚溫室建立的CFD模型中的微氣候進行了研究。為了模擬棚頂及覆蓋層的對流轉移和輻射交流的耦合作用，在原來的CFD模型中加入了輻射模型，代替了通常所用的基于能量守恒的耦合方法。并且在作物冠層建模中獨創性地采用了分形滲透模型。對數值模擬結果與測量實驗數據的分析比較表明，該模型是成功的。隨后該模型被用于探究溫室內微氣候變量分布，表明此類型溫室內的氣流形式、溫度、濕度都與鋸齒狀溫室中不同。該研究表明這種新穎的CFD模型可作為溫室設計及微氣候研究的有效工具。

此外，還有部分學者將CFD模擬與優化控制相結合來實現對溫室環境的控制。如李永博等[31]在CFD模擬的基礎上引入均勻性評價指標，通過數值求解得出溫室內溫度控制的最優化方案，指導溫室內溫度的精細控制。此方法對于溫室內其他環境因子的研究同樣適用。

3結論

CFD研究對溫室內定點信息的預測更加精細、詳盡，不僅能預測即時信息，還能夠預測未來一段時間內溫室環境因子的分布情況，能夠從機理出發指導溫室的通風結構設計、灌溉、施肥等操作。此外，CFD研究與常規的實驗研究相比，能大大減少傳感設備的布控，降低成本。因而CFD研究是實現溫室精細化控制的強有力工具。

CFD研究在處理有植物生長的溫室時仍具有一定的局限性。目前已有的研究中，大多將植物看作多孔介質，植物的蒸騰作用、光合作用、呼吸作用等生理作用通過源相來導入。而源相中關于植物的作用往往是依據現有的作物生長模型來處理，如對植物蒸騰作用用Penman-Monteith方程來描述。這就涉及到各種參數的確定，且對于不同作物或同一作物的不同生長階段，這些參數均有差別，因此利用CFD對溫室內植物冠層的處理有一定的難度。對植物與溫室環境的作用機制及各種相關聯參數進行進一步的研究，將對溫室生產的精細化控制產生更深遠的意義。

參考文獻：

[1] 楊聯安.我國精確農業發展策略初探[J].測繪與空間地理信息,2004,27(2):3-5.

[2] 趙國強,王君,穆曉濤.精細農業研究進展[J].氣象與環境科學,2007,30(1):84-88.

[3] Blackmore S, Godwin R, Fountas S. The analysis of spatial and temporal trends in yield map data over six years[J]. Biosystems Engineering, 2003, 84(4): 455-466.

[4] 李民贊.精細農業特征、現狀與發展[J].農機市場,2008(3):38-41.

[5] 劉剛.精細農業與GPS農業應用[J].全球定位系統,2000,25(2):12-16.

[6] 許安祥.我國溫室農業的發展現狀與趨勢[J].排灌機械,2000，18(1):4-8.

[7] 曾飚,張凱,黃丹楓.空氣熱交換系統在溫室環境控制中的應用[J].上海交通大學學報,2012,30(5):40-44.

[8] 劉珊,張林華,張峰.日光溫室熱濕環境影響機理研究進展[J].山東建筑大學學報:2009,24(6):587-593.

[9] Pierre-Emmanuel Bounet, Thietty Boulard. Effect of ventilator configuration on the distributed climate of greenhouses: A review of experimental and CFD studies[J]. Computers and Electronics in Agriculture, 2010, 74(2): 195-217.

[10] Mistriotis A, Arcidiacono C, Picuno P, et al. Computational analysis of ventilation in greenhouse at zero and low-wind-speeds[J]. Agricultural & Forest Meteorology, 1997, 88(1/4): 121-135.

[11] 沈明衛,郝飛麟.連棟塑料溫室自然通風流場的穩態模擬——頂窗配置形式對溫室內流場的影響[J].浙江大學學報,2007,41(4):668-672.

[12] 郭萬東,徐泳,張天柱.溫室結構風壓的CFD模擬[J].中國農業大學學報,2005,10(1):56-59.

[13] Francisco Domingo Molina-Aiz, Diego Luis Valera, Antoio Jesus Alvarez. Measurement and simulation of climate inside Almeria-type greenhouse using computational fluid dynamics[J]. Agricultural and Forest Meteorology, 2004，125(1/2): 33-51.

[14] Davide Piscia, Juan I Montero, Esteban Baeza, et al. A CFD greenhouse night-time condensation model[J]. Biosystems Engineering, 2012,111(2): 141-154.

[15] 譚勝男,汪小旵,保智敏,等.溫室內噴霧降溫系統的CFD模擬[J].江蘇農業學報,2013,29(2):283-287.

[16] Bartzanas T, Boulard T, Kittas C. Effect of vent arrangement on windward ventilation of a tunnel greenhouse[J]. Biosystems Engineering, 2004, 88(4): 479-490.

[17] Sapounas A, Nikita C, Bartzanas T, et al. Aspects of CFD modeling of fan and pad evaporative cooling system in greenhouses[M]. 2nd PALENC Conf. and 28th AIVC Conf. on Building Low Energy Cooling and Advanced Ventilation Technologies in the 21st Century, Crete Island, Greece, 2007.

[18] 程秀花,毛罕平,倪軍.溫室環境—作物濕熱系統CFD模型構建與預測[J].農業機械學報,2011,42(2):172-179.

[19] Abderzak Kichah, Pierre-Emmanuel Bournet, Christophe Migeon, et al. Measurement and CFD simulation of microclimate charecteristics and transpiration of an Impatiens pot plant crop in a greenhouse [J]. Biosystems engineering, 2012, 112(1): 22-34.

[20] Fodaros D K, Baxevaou C A, Bartzanas T, et al. Numerical simulation of thermal behavior of a ventilated arc greenhouse during a solar day[J]. Renewable Energy, 2010, 35(7): 1380-1386.

[21] Bartzanas T, Kacira M, Zhu H, et al. Computational fluid dynamics applications to improve crop production systems[J].Computers and Electronics in Agriculture, 2013, 93(93): 151-167.

[22] Campen J B, Bot G P A. Determination of greenhouse-specific aspects of ventilation using three-dimensional computational fluid dynamics[J]. Biosystems Engineering, 2003, 84(1): 69-77.

[23] Wang S, Boulard T. Measurement and prediction of solar radiation distribution in full scale greenhouse tunnel[J]. Agronomie , EDP Sciences, 2000, 20(1): 41-50.

[24] Montero J I, Munoz P, Anton A, et al. Computational fluid dynamic modelling of night-time energy fluxes in unheated greenhouses[J]. Acta Horticulturae , 2005, 693: 403-410.

[25] Bournet P E, Ould Khaoua S A, Boulard T. Numerical prediction of the effect of vent arrangements on the ventilation and energy transfer in a multi-span glasshouse using a bi-band radiation model[J]. Biosystems Engineering, 2007, 98(2): 224-234.

[26] Kim K, Jeong-Yeol Yoon, Hyuck-Jin Kwon, et al. 3-D CFD analysis of relative humidity distribution in greenhouse with a fog cooling system and refrigerative dehumidifiers[J]. Biosystems Engineering, 2008, 100(2): 245-255.

[27] Lee I B, T H Short. Predicted effects of internal horizontal screens on natural ventilation of a multi-span greenhouse//[C]. In: 91st Annual International Meeting of ASAE. Orlando ( USA), 1998.

[28] Chen J L, Xu F, Tan D P, et al. A control method for agricultural greenhouses heating based on computational fluid dynamics and energy prediction model[J]. Applied Energy, 2015, 141(1): 106-118.

[29] 譚勝男.基于CFD的現代化溫室環境數值模擬與優化研究[D].南京:南京農業大學,2013.

[30] Wang X W, Luo J Y, Li X P. CFD based study of heterogeneous microclimate in a typical chinese greenhouse in central China[J]. Journal of Integrative Agriculture, 2013, 12(5): 914-923.

[31] 李永博,周偉,李鵬飛,等.基于CFD模擬模型的溫室溫度場均勻性控制[J].農業機械學報, 2012, 43(4):156-161.

(責任編輯：管珊紅)

收稿日期：2016-01-14

基金項目：工信部2012年物聯網專項資金資助。

作者簡介：張建文(1969—)，男，山西夏縣人，教授，博士，研究方向：流體力學與傳熱。

中圖分類號：S625

文獻標志碼：A

文章編號：1001-8581(2016)07-0126-05

Research Advance in Multiphase Flow Transfer Process of Agricultural Greenhouse System

ZHANG Jian-wen1, HE Ya-ru1, WANG Qing-gui2

(1. Lab of Fluid Flow and Heat Transfer, Beijing University of Chemical Technology, Beijing 100029, China;2. Shannanyalong Green Animal & Plant Limited Company in Tibet, Shannan 856100, China)

Abstract：The author analyzed the characteristics of multiphase flow transfer process of agricultural greenhouse system, expounded several problems in complex multiphase flow and heat transfer inside the greenhouse, further reviewed the current researches on computational fluid dynamics (CFD) of heat and mass transfer process in the greenhouse, summarized the advantages of CFD in the study of multiphase flow and heat-mass transfer in greenhouse system, and discussed the prospects of CFD in the study of greenhouse multiphase flow system.

Key words：Precision agriculture; Greenhouse system; Multiphase fluid flow; Heat and mass transfer; CFD