方 慧，伍 綱，程瑞鋒，仝宇欣

基于CFD的植物工廠冠層微環(huán)境管道通風(fēng)效果模擬

方 慧，伍 綱，程瑞鋒，仝宇欣※

（1. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所，北京 100081；2. 農(nóng)業(yè)農(nóng)村部設(shè)施農(nóng)業(yè)節(jié)能與廢棄物處理重點(diǎn)實(shí)驗(yàn)室，北京 100081）

為增加植物工廠多層栽培模式中作物冠層內(nèi)部氣流擾動(dòng)，該研究設(shè)計(jì)了一種冠層微環(huán)境管道通風(fēng)裝置，利用計(jì)算流體力學(xué)軟件（Computational Fluid Dynamics，CFD）構(gòu)建三維栽培模型。在模型中，將植物冠層區(qū)域考慮為多孔介質(zhì)，多孔介質(zhì)的黏滯阻力系數(shù)為25，慣性阻力系數(shù)為1.3；LED燈管設(shè)置為熱源模型，熱源模型散熱量為297 525 kW/m3。利用模型模擬并實(shí)測入口速度為8 m/s時(shí)植物冠層表面和上部空間氣流速度，發(fā)現(xiàn)空氣流域的平均誤差為16%，模擬值與實(shí)測值吻合較好。利用驗(yàn)證的模型模擬不同入口速度下植物冠層的氣流分布，結(jié)果表明：進(jìn)氣速度設(shè)置為8 m/s，該速度下植物冠層內(nèi)部適宜速度區(qū)域體積占比為67.58%，冠層表面適宜速度區(qū)域面積占比為55.23%，冠層內(nèi)部氣流平均速度為0.14 m/s。研究表明植物工廠冠層微環(huán)境管道通風(fēng)模式能有效增加冠層氣流擾動(dòng)。

通風(fēng)；植物工廠；氣流速度；多孔介質(zhì)；CFD

0 引 言

植物工廠是一種通過設(shè)施內(nèi)環(huán)境高精度控制實(shí)現(xiàn)作物周年連續(xù)生產(chǎn)的高效農(nóng)業(yè)生產(chǎn)方式，單位面積蔬菜產(chǎn)能可達(dá)露地40倍以上，被認(rèn)為是保障食物安全的重要途徑之一，受到世界各國的廣泛重視[1-4]。規(guī)模化、集約化和智能化的植物工廠生產(chǎn)方式已成為未來農(nóng)業(yè)的發(fā)展趨勢。現(xiàn)階段，植物工廠環(huán)境調(diào)控主要通過氣流循環(huán)解決，送風(fēng)的方式對植物工廠微環(huán)境有著決定性作用，直接影響系統(tǒng)內(nèi)部的溫度、氣流速度、區(qū)域溫差以及通風(fēng)調(diào)溫系統(tǒng)能耗等方面。

在植物工廠生產(chǎn)中為節(jié)省室內(nèi)空間，植物工廠以多層式栽培方式為主。一般栽培床層間距為30～40 cm，去除作物的生長高度，作物冠層以上空間非常有限，盡管植物工廠設(shè)計(jì)的氣流較大，但由于作物和栽培床的阻擋，導(dǎo)致氣流在作物周圍形成繞流，在栽培架的中間位置，氣流速度幾乎為零，使得植物生長環(huán)境中氣流不均勻，影響作物蒸騰速率和光合作用，容易誘導(dǎo)植物生理性病害發(fā)生，尤其以葉燒現(xiàn)象最為普遍[5-8]。Hesham等[9]研究認(rèn)為植物工廠葉燒的主要原因是內(nèi)部氣流分布不合理引起。由于植物工廠內(nèi)部氣流分布不均勻，過道風(fēng)速較大，植物冠層風(fēng)速較低，從而導(dǎo)致植物運(yùn)輸鈣離子的能力降低，加之植物工廠內(nèi)部作物生長速率較快，使得作物新葉片鈣的濃度降低，從而使植物葉片缺鈣形成葉燒，該缺鈣現(xiàn)象并非是由于植物根部吸收鈣不足導(dǎo)致的，而是由于植物對鈣離子的轉(zhuǎn)運(yùn)能力不足，導(dǎo)致沒有足夠的鈣來供給植物生長從而發(fā)生葉燒現(xiàn)象[9]。為減少葉燒現(xiàn)象，最有效的辦法是增加作物冠層氣流速度。在作物氣流速度研究方面，Kitaya等[10-12]通過試驗(yàn)表明氣流速度在0.02～1.30 m/s范圍內(nèi)，隨著氣流速度的提升，黃瓜葉片凈光合速率和蒸騰速率作用分別增加了1.2和2.8倍，該研究結(jié)果表明作物冠層氣流速度的增加，對作物光合和蒸騰速率有積極的影響。作物冠層的低速區(qū)域會(huì)影響作物蒸騰速率及其與環(huán)境間的對流換熱，進(jìn)而影響作物微環(huán)境溫度[13-14]。

為優(yōu)化植物冠層區(qū)域氣流速度分布，設(shè)計(jì)合理的植物工廠通風(fēng)模式，本文從通風(fēng)管道氣流分流原理出發(fā)，理論上分析氣孔分布對氣流均勻性的影響，以確定氣孔數(shù)量及分布；以植物工廠種植生菜為例（L...）通過CFD模擬試驗(yàn)，探明在植物和LED燈管影響下，最優(yōu)的進(jìn)氣速度，以期為植物工廠通風(fēng)設(shè)計(jì)提供參考。

1 計(jì)算模型與驗(yàn)證試驗(yàn)

1.1 模型構(gòu)建與網(wǎng)格劃分

植物工廠單層栽培架如圖1所示。栽培空間高度為400 mm，栽培板長1 500 mm，寬700 mm，共32個(gè)定植孔（4行8列），孔橫向間距為130 mm，縱向間距為170 mm；LED燈管長1 200 mm，7根燈管均勻布置于頂部，燈管間距為80 mm。通氣管道主管布置于栽培槽一側(cè)，支管布置于栽培孔行間，共3行，氣流通過支管上的通氣孔進(jìn)入植物冠層。根據(jù)Colin等[15]研究表明多孔通氣管的管徑比（氣孔總面積/管橫截面積）為1.0～1.5時(shí)，每個(gè)通氣孔的氣流分配較均勻。通氣裝置主管道選用內(nèi)徑為36 mm的PVC軟管，根據(jù)孔徑比要求，支管選用內(nèi)徑為20 mm的PVC管。為保證氣流能均勻流入到植物冠層，根據(jù)支管孔徑比要求和常規(guī)鉆頭尺寸，每根支管上設(shè)置個(gè)38個(gè)通氣孔（2行19列），孔徑為4 mm，孔間距為70 mm，計(jì)算得到通氣孔的管徑比為1.5。

利用CFD模型中的Design Modeler模塊構(gòu)建單層栽培空間模型，將構(gòu)建的模型輸入到ANSYS Meshing模塊中，利用Proximity and Curvature方法對幾何體進(jìn)行網(wǎng)格劃分，并對支管出氣孔和LED燈管進(jìn)行網(wǎng)格進(jìn)行加密處理，共生成4 126 129個(gè)網(wǎng)格，820 592個(gè)節(jié)點(diǎn)，利用skewness計(jì)算出網(wǎng)格最大偏斜度為0.7，最小偏斜度為6.0×10-6，平均值為0.23，網(wǎng)格質(zhì)量優(yōu)，可用于案例模擬分析[16]。

1.2 模擬方程

1.3 材料屬性與邊界條件

根據(jù)前期測試，通風(fēng)管道布置于栽培層一側(cè)，管徑為50 mm時(shí)，管道入口處的最優(yōu)氣流速度應(yīng)為3 m/s[19]，即最優(yōu)氣體流量為6.0×10-3m3/s。本研究為增加植物內(nèi)部氣流擾動(dòng)，將通風(fēng)管道設(shè)置于植物冠層內(nèi)部，入口處通風(fēng)管內(nèi)徑為36 mm，以氣體總流量為定值，則入口風(fēng)速需達(dá)到5.9 m/s。本設(shè)計(jì)將通風(fēng)管道直接布置于植物冠層內(nèi)部，勢必會(huì)增加氣流阻力，因此需增加氣流速度。為優(yōu)化植物冠層內(nèi)氣流速度分布，本模擬研究設(shè)置4個(gè)氣流入口速度，分別為6、7、8和9 m/s，入口氣流溫度設(shè)定為24 ℃。栽培裝置模型上端和下端設(shè)置為壁面（Wall）；四周設(shè)置為壓力出口（Pressure outlet），出口回流溫度為操作環(huán)境溫度24 ℃。模擬計(jì)算中涉及的材料包括空氣、植物、鋁和玻璃，各材料熱物理屬性參數(shù)如表1所示。

表1 材料熱物理參數(shù)

模型中為考慮植物阻擋和燈管散熱對氣流的影響，將植物考慮為多孔介質(zhì)，燈管設(shè)置為熱源邊界。

1.3.1 植物多孔介質(zhì)參數(shù)

植物對氣流有一定的阻力，會(huì)導(dǎo)致氣流動(dòng)量損失，進(jìn)而對植物工廠內(nèi)部通風(fēng)環(huán)境產(chǎn)生負(fù)面影響。在CFD模型中將生菜假設(shè)為多孔介質(zhì)，根據(jù)Darcy-Forchheimer定律，生菜對氣流的阻擋壓降梯度為

1.3.2 LED燈管熱源參數(shù)

LED燈管電光轉(zhuǎn)化效率與芯片材料和制作過程有關(guān)。Nelson等[21]測試不同類型LED燈珠發(fā)現(xiàn)，紅色燈珠和藍(lán)色燈珠的電光轉(zhuǎn)化效率分別為32%和49%。本試驗(yàn)以植物工廠中常用紅藍(lán)組合LED燈管為例計(jì)算燈管熱功率。LED燈管由220顆紅燈珠和20顆藍(lán)組成（型號：IGL-T8-36B1R11，廣州市智卉光田農(nóng)業(yè)科技有限責(zé)任公司），紅燈珠電功率為0.15 W，藍(lán)燈珠電功率為0.21 W。燈管散熱功率計(jì)算式為

LED燈管由燈罩、燈珠芯片和散熱鋁板組成。在CFD模型中可將LED燈管簡化為玻璃燈罩和鋁基板散熱器[20]，LED燈管模型如圖2所示，燈管直徑為20 mm，鋁基板截面寬度為17 mm，高度為4 mm，散熱器熱流密度計(jì)算式為

1.4 驗(yàn)證試驗(yàn)

為驗(yàn)證模型，于2021年6月20日在試驗(yàn)型植物工廠栽培架中布置通風(fēng)管道進(jìn)行通風(fēng)試驗(yàn)。植物工廠位于中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所內(nèi)，植物工廠規(guī)格（長×寬×高）為4.1 m×2.3 m×2.6 m。管道風(fēng)機(jī)轉(zhuǎn)速通過調(diào)壓器調(diào)節(jié)（TDGC2-0.5KVA，德碩智能電器有限公司）。通風(fēng)管入口處氣流速度為8 m/s，氣流溫度為24 ℃，植物工廠環(huán)境溫度為24 ℃。對植物工廠單層栽培裝置中成熟期生菜冠層內(nèi)部和上部空氣流域的氣流速度進(jìn)行測試，測點(diǎn)布置如圖3所示。在中間位置的LED燈管正下方布置30個(gè)測點(diǎn)，分布于栽培層上部、中部和生菜冠層表面，用于記錄氣流速度。氣流速度測量選用紅外熱線風(fēng)速儀（Climomaster6501-BG，日本加野麥克斯公司），測量范圍為0.01～5.00 m/s，測量精度為標(biāo)準(zhǔn)值的±2%。

上部空氣流域測點(diǎn)Measured points in the up air interior；中部空氣流域測點(diǎn)Measured points in the middle air interior；冠層表面測點(diǎn)Measured points on canopy Surface

2 結(jié)果與分析

2.1 模型驗(yàn)證

為驗(yàn)證模型結(jié)果的準(zhǔn)確性，在試驗(yàn)型植物工廠中進(jìn)行管道通風(fēng)試驗(yàn)，試驗(yàn)期生菜種植密度為22棵/m2。風(fēng)管進(jìn)氣口氣流速度為8 m/s，試驗(yàn)裝置運(yùn)行穩(wěn)定后取溫度值和氣流速度值。

氣流速度模擬結(jié)果與實(shí)測結(jié)果對比如圖4所示。從圖4可以看出冠層上部空氣流域的氣流速度模擬值與實(shí)測值較接近。空氣流域20個(gè)氣流速度測點(diǎn)的平均誤差為16%；植物冠層10個(gè)氣流速度測點(diǎn)的平均誤差較大為23%?？諝饬饔虻哪M值與實(shí)測值比較接近，但在植物冠層表面模擬值與實(shí)測值誤差較大，其主要原因是在模型中將植物設(shè)置為均勻的多孔介質(zhì)，而在實(shí)際生產(chǎn)中植物并非均勻物質(zhì)，植株內(nèi)部孔隙度較小，而植株與植株間的孔隙度較大，且由于氣流的干擾，植物呈現(xiàn)波動(dòng)狀態(tài)，進(jìn)而影響氣流速度。

2.2 進(jìn)氣速度對冠層表面氣流分布影響

植物工廠常規(guī)通風(fēng)模式中，由于作物冠層外部葉片對氣流的阻擋，導(dǎo)致冠層內(nèi)部風(fēng)速小，甚至出現(xiàn)停滯區(qū)，通風(fēng)效率較低[22]。為了突破冠層邊界層阻力及葉片封閉效應(yīng)，本試驗(yàn)將通風(fēng)管道置于栽培板上。氣流通過通風(fēng)管道氣孔進(jìn)入作物冠層，由于葉片的阻擋，氣流在作物冠層內(nèi)部發(fā)生繞流，從葉片空隙流入到周圍環(huán)境中。圖5為不同進(jìn)氣速度下生菜冠層表面氣流分布云圖。從圖5可以看出延氣流流向，冠層表面氣流速度較大。隨著入口氣流速度的增加，高速區(qū)域面積逐漸向外延展擴(kuò)大。通過計(jì)算冠層表面面積加權(quán)的平均速度得到，當(dāng)進(jìn)氣速度為6、7、8和9 m/s時(shí)，冠層表面的平均速度分別為0.11、0.13、0.15和0.17 m/s。

Kitaya等[23]的研究表明植物冠層內(nèi)部氣流速度從0.01增加到0.2 m/s時(shí)，能有效增加植物蒸騰速率和凈光合速率，氣流速度從0.2增加到1.0 m/s，蒸騰速率持續(xù)增加，但當(dāng)氣流速度從0.5增加到1.0 m/s時(shí)，凈光合速度穩(wěn)定不變。Goto等[24]的研究表明植物冠層表面氣流速度在0.3～1.0 m/s時(shí)，能有效降低葉燒病發(fā)生率。為比較不同入口速度下植物冠層氣流分布特征，將冠層氣流速度分為停滯區(qū)（<0.1 m/s）、適宜區(qū)（0.1 m/s≤≤1.0 m/s）和高速區(qū)（>1.0 m/s），利用CFD軟件計(jì)算植物冠層氣流速度分別在停滯區(qū)、適宜區(qū)和高速區(qū)的面積占冠層總面積的比例，如表2所示。當(dāng)進(jìn)氣流速為6 m/s時(shí)，冠層表面停滯區(qū)域的面積較大，占總表面積的66.67%。為減少停滯區(qū)域面積，將進(jìn)氣速度逐步提高到7、8和9 m/s，并利用CFD進(jìn)行模擬，發(fā)現(xiàn)隨著進(jìn)氣速度的增加，停滯區(qū)面積所占比例由66.67%降低到31.43%，適宜區(qū)面積所占比例由33.33%增加到68.57%，且冠層表面未出現(xiàn)高速區(qū)域。

表2 不同進(jìn)氣速度下冠層表面氣流速度分布比較

2.3 進(jìn)氣速度對冠層內(nèi)部氣流分布影響

氣流經(jīng)管道小孔進(jìn)入作物冠層內(nèi)部，增加冠層內(nèi)部氣流擾動(dòng)，在一定程度上能改善植物內(nèi)部環(huán)境。在冠層內(nèi)部軸方向選取=224 、715和1 274 mm截面如圖6所示，截面通過延氣流方向的第3排、10排和第18排通氣孔。速度分布云圖如圖7所示，由圖7可以看出，不同進(jìn)氣速度下植物冠層內(nèi)部氣流分布輪廓一致。氣流在管孔附近速度較大，進(jìn)入植物冠層區(qū)域后逐漸衰減，在氣孔夾角間存在停滯區(qū)域。通過計(jì)算冠層內(nèi)部體積加權(quán)平均值發(fā)現(xiàn)，進(jìn)口處氣流速度為6、7、8和9 m/s時(shí)，植物冠層內(nèi)部的平均速度分別為0.10、0.12、0.14和0.16 m/s。

隨著入口速度的增加，停滯區(qū)域體積所占比例減少，但高速區(qū)域體積所占比例增加。進(jìn)口處氣流速度分別為6、7、8和9 m/s時(shí)，速度適宜區(qū)域體積占比分別為42.84%、52.36%、67.58%和77.08%（表3），高速區(qū)域體積比分別為0.01%、0.02%、0.04%和0.06%。雖然高速區(qū)域體積比較小，但通過測量發(fā)現(xiàn)，當(dāng)氣流速度由6增加到9 m/s時(shí)，氣孔高速區(qū)域與管道外壁最大垂直距離由40增加到60 mm，通風(fēng)管道半徑為10 mm，即高速區(qū)域距離管中心的距離由50增加到了70 mm。在實(shí)際生產(chǎn)中，栽培孔間距為130 mm，通風(fēng)管道布置于植物株間，管孔附近的高速區(qū)域距離管中心距離大于65 mm時(shí)對植物生長有影響。為避免植物冠層內(nèi)部氣流速度過大，同時(shí)滿足冠層內(nèi)部適宜的氣流速度，選擇進(jìn)氣口最佳氣流速度為8 m/s。

表3 冠層內(nèi)部氣流速度分布比較

3 討 論

本研究的目的是探討植物工廠采用冠層微環(huán)境管道通風(fēng)裝置后不同進(jìn)氣速度對生菜冠層氣流分布及均勻性的影響。對于單個(gè)栽培層，植物和補(bǔ)光燈對氣流的影響密不可分。本文提出將生菜簡化為多孔介質(zhì)，將補(bǔ)光燈設(shè)置為熱源，模擬植物冠層內(nèi)部及上部空間的氣流速度。通過模擬和實(shí)測的數(shù)據(jù)對比發(fā)現(xiàn)，在植物冠層內(nèi)部和上部空間模擬值與實(shí)測值吻合較好，但在植物冠層上表面，兩者差異較大，主要原因是模擬中將植物考慮為均勻的多孔介質(zhì)，而實(shí)際生產(chǎn)中作物冠層之間空隙分布并不均勻，且氣流流動(dòng)過程中作物也會(huì)隨之?dāng)[動(dòng)，勢必會(huì)影響氣流走向。

目前在畜禽舍通風(fēng)研究方面，為更能真實(shí)反映畜禽對氣流的影響，探索性的將畜禽考慮為實(shí)體模型，如程瓊儀[25]在測試雞群對氣流阻力時(shí)，將雞身簡化為橢球型、輕微簡化型和雞身型，該方法更能精確模擬畜禽舍中氣流分布。植物不同于雞群，不僅植株與植株間存在空隙，植株內(nèi)部亦存在空隙，且植物葉片厚度小，在網(wǎng)格劃分中若將植物考慮為實(shí)體模型，勢必會(huì)增加網(wǎng)格數(shù)量，降低計(jì)算速度。因此，需對植物外形進(jìn)行輕微簡化，以期在不增加網(wǎng)格數(shù)量下實(shí)現(xiàn)對冠層氣流分布精確模擬。

4 結(jié) 論

本研究針對現(xiàn)有植物工廠中氣流主要集中在過道，栽培區(qū)域出現(xiàn)停滯區(qū)的問題，設(shè)計(jì)了一種局部管道通風(fēng)裝置，并利用CFD軟件模擬了不同進(jìn)氣速度下植物冠層內(nèi)部和冠層表面氣流分布，為植物工廠通風(fēng)設(shè)計(jì)提供參考。研究主要結(jié)論如下：

1）利用CFD軟件構(gòu)建了植物工廠冠層微環(huán)境管道通風(fēng)模型，在模型中將植物冠層考慮為多孔介質(zhì)，多孔介質(zhì)的黏滯阻力系數(shù)設(shè)置為25，慣性阻力系數(shù)設(shè)置為1.3，將LED燈管考慮為熱源，熱流密度設(shè)置為297 525 W/m3。通過模擬和實(shí)測進(jìn)氣速度為8 m/s，進(jìn)氣溫度為24 ℃條件下植物冠層表面氣流速度和冠層上部空間氣流速度，得到空間流域模擬值與試驗(yàn)值的平均誤差為16%，構(gòu)建的植物工廠冠層微環(huán)境管道通風(fēng)模型能準(zhǔn)確描述進(jìn)氣速度對植物冠層上部空氣間氣流速度的影響。

2）利用驗(yàn)證的模型模擬不同入口氣流下植物冠層的氣流分布，結(jié)果表明：植物工廠冠層微環(huán)境管道通風(fēng)裝置能有效增加作物冠層內(nèi)部氣流擾動(dòng)，當(dāng)入口氣流速度從6增加到9 m/s時(shí)，植物冠層表面的平均速度由0.11 m/s增加到0.17 m/s，冠層表面適宜速度區(qū)域面積占比由33.33%增加到68.57%，植物冠層內(nèi)部的平均速度由0.11增加到0.17 m/s，冠層內(nèi)部適宜速度區(qū)域體積占比由42.84%增加到77.08%。為避免冠層內(nèi)部高速氣流對植物的影響，可將進(jìn)氣速度設(shè)置為8 m/s。

[1] 楊其長. 植物工廠[M]. 北京：清華大學(xué)出版社，2019.

[2] Kozai T. Smart Plant Factory: The Next Generation Indoor Vertical Farms[M]. Singapore: Springer, 2018.

[3] Kozai T. Plant factory in Japan current situation and perspectives[J]. Chronica Horticulturae, 2013, 53(2): 8-11.

[4] Kozai T, Chun C, Ohyama K. Closed systems with lamps for commercial production of transplants using minimal resources[J]. Acta Horticulturae, 2004, 630: 239-254.

[5] Lee I B. Trend on CFD application in agriculture[J]. Rural Resources, 2013, 55(1): 51-57.

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

Liu Huan. Research on Ventilation Simulation and Optimization of Plant Factory with Artificial Light Based on CFD[D]. Beijing: Chinese Academy of Agricultural Science, 2019. (in Chinese with English abstract)

[7] 王君，楊其長，魏靈玲，等. 人工光植物工廠風(fēng)機(jī)和空調(diào)協(xié)同降溫節(jié)能效果[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(3)：177-183.

Wang Ju, Yang Qichang, Wei Lingling, et a1. Energy saving effect on cooperating cooling of conditioner and air exchanger in plant factory with artificial light[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(3): 177-183. (in Chinese with English abstract)

[8] 張宇，肖玉蘭. 人工光植物工廠生產(chǎn)成本構(gòu)成和經(jīng)濟(jì)效益分析[J]. 長江蔬菜，2017，(4)：34-40.

Zhang Yu, Xiao Yulan. Analysis on production costs and economic benefits for production of plant factory with artificial light[J]. Journal of Changjiang Vegetables, 2017(4): 34-40. (in Chinese with English abstract)

[9] Hesham A A, Tong Y, Yang Q. Lettuce plant growth and tipburn occurrence as affected by airflow using a multi-fan system in a plant factory with artificial light[J]. Journal of Thermal Biology, 2020, 88: 1-12.

[10] Kitaya Y, Tsuruyama J, Shibuya T, et al. Effects of air current speed, light intensity and CO2concentration on photosynthesis and transpiration of plant leaves[C]. In 35th COSPAR Scientific Assembly. Paris: 35th COSPAR Scientific Assembly, 2004: 16478-16485.

[11] Omasa K, Nouchi I, De Kok L J. Plant Responses to Air Pollution and Global Change[M]. Tokyo: Spring, 2005.

[12] 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.

[13] Baek M S, Kwon S Y, Lim J H. Improvement of the crop growth rate in plant factory by promoting airflow inside the cultivation[J]. International Journal of Smart Home, 2016, 10(2): 63-74.

[14] 李琨，鄒志榮. 植物工廠生菜根際通風(fēng)對冠層及根際環(huán)境影響. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(7)：178-187.

Li Kun, Zou Zhirong. Environmental effects of root zone ventilation on canopy and rhizosphere of lettuce in plant factory[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(7): 178-187. (in Chinese with English abstract)

[15] Colin M W, Nevin D A. Design of air distribution systems for closed greenhouses[J]. Acta Horticulturae(ISHS), 1994, 361: 93-104.

[16] 李海峰，吳冀川，劉建波，等. 有限元網(wǎng)格剖分與網(wǎng)格質(zhì)量判定指標(biāo)[J]. 中國機(jī)械工程，2012，23(3)：368-377

Li Haifeng, Wu Jichuan, Liu Jianbo, et al. Finite element mesh generation and mesh quality evaluation index[J]. China Mechanical Engineering, 2012, 23(3): 368-377. (in Chinese with English abstract)

[17] Shih T H, Liou W W, Shabbir A, et al. A new k-ε eddy viscosity model for high Reynolds number turbulent ?ow[J]. Computers & Fluids, 1995, 24(3): 227-238.

[18] 程秀花. 溫室環(huán)境因子時(shí)空分布CFD模型構(gòu)建及預(yù)測分析研究[D]. 鎮(zhèn)江：江蘇大學(xué)，2011.

Cheng Xiuhua. CFD Model Construction and Prediction Analysis of Spatial and Temporal Distribution of Greenhouse Environmental Factors[D]. Zhenjiang: Jiangsu University, 2011. (in Chinese with English abstract)

[19] 吳晨溶，程瑞鋒，方慧，等. 基于CFD的植物工廠管道通風(fēng)模擬及優(yōu)化[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2021，26(1)：78-87.

Wu Chenrong, Cheng Ruifeng, Fang Hui, et al. Simulation and optimization of air tube ventilation in plant factory based on CFD[J]. Journal of China Agricultural University, 2021, 26(1): 78-87. (in Chinese with English abstract)

[20] 張晨，方慧，程瑞鋒，等. 基于風(fēng)洞系統(tǒng)的生菜空氣動(dòng)力學(xué)研究[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2019，24(12)：96-103.

Zhang Chen, Fang Hui, Cheng Ruifeng, et al. Wind tunnel measurement of aerodynamic properties of a lettuce canopy[J]. Journal of China Agricultural University, 2019, 24(12): 96-103. (in Chinese with English abstract)

[21] Nelson J A, Bugbee B. Supplemental greenhouse lighting: Return on investment for LED and HPS fixtures[J]. Controlled Environments, 2013, 7: 1-9.

[22] Fang H, Li K, Wu G, et al. A CFD analysis on improving lettuce canopy airflow distribution in a plant factory considering the crop resistance and LEDs heat dissipation[J]. Biosystems Engineering, 2020, 200(12): 1-12.

[23] 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-82.

[24] Goto E, Takakura T. Prevention of lettuce tipburn by supplying airto inner leaves[J]. Transactions of the ASABE, 1992, 35: 641-645.

[25] 程瓊儀. 疊層籠養(yǎng)蛋雞舍夏季通風(fēng)氣流CFD模擬與優(yōu)化[D]. 北京：中國農(nóng)業(yè)大學(xué)，2018.

Cheng Qiongyi. The CFD Simulation and Improvements of Ventilation in Stacked-Cage Laying Hen Housing[D]. Beijing: China Agricutural University, 2018. (in Chinese with English abstract)

Numerical simulation of ventilation efficiency with perforated air tubes for controlling canopy environment in a plant factory based on CFD

Fang Hui, Wu Gang, Cheng Ruifeng, Tong Yuxin※

(1.,,100081,; 2.,,100081,)

Much more cultivation layers and high canopy density have posed great challenges to the even distribution of airflow in plant factories. It is highly required for the optimal design of airflow distribution systems, thereby improving the physiological orders and post-harvest quality for a high commercial value of the products. In this study, a ventilation system was designed with the three parallel lines of perforated air tubes, further to generate a constant airflow inside the plant canopy. The air tube was fabricated in a diameter of 20 mm with a length of 1 500 mm using the layout and size of a cultivation shelf. Each air tube was equipped with two rows of jets. A three-dimensional computational fluid dynamics (CFD) model was developed to analyze the distribution of airflow in the plant’s canopy under different inlet velocities. Furthermore, the lettuce leaves were assumed as the porous medium with the drag coefficient of 0.02, resulting in a viscous resistance factor and the inertial resistance factor of 25 and 1.3, respectively. A LED lamp (?20 mm) was assumed as an aluminum substrate and a glass part. An energy term with a constant heat generation of 297 525 W/m3was also defined as the heat source lost by the volumetric convection from the LED lamp aluminum substrate part for the heat dissipation. The average error of simulated air velocity was 16%, compared with the measured. Three groups of air velocity (lower than 0.1 m/s, between 0.1 and 1.0 m/s, and higher than 1.0 m/s) were selected to calculate the proportion of the areas of the plant’s canopy surface and the volume proportion inside the plant’s canopy. The results showed that the percentage of air velocity at the lettuce canopy surface between 0.1 and 1 m/s increased from 33.33% to 68.57%, while, the percentage of volumes inside the plant’s canopy increased from 42.84% to 77.08%, respectively, when the air velocity increased from 6 to 9 m/s. By contrast, the percentage of volumes inside the plant’s canopy increased from 0.01% to 0.06%, with the air velocity higher than 1.0 m/s. The maximum distance between the ventilation pipe and high airflow increased from 50 to 70 mm. Therefore, the air velocity was recommended as 8 m/s, further to avoid the air velocity higher than 1 m/s plants canopy. At the inlet air flow velocity, the percentage of volumes inside the plant’s canopy with the air velocity between 0.1 and 1 m/s was 67.58%, the percentage of areas on the plant’s canopy surface with the air velocity between 0.1 and 1.0 m/s was 55.23%, and the volume-weighted air velocity in the canopy was 0.14 m/s. This finding can provide a new solution to improve the canopy environment in a plant factory.

ventilation; plant factory; airflow velocity; porous media; CFD

方慧，伍綱，程瑞鋒，等. 基于CFD的植物工廠冠層微環(huán)境管道通風(fēng)效果模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(22)：188-193.doi：10.11975/j.issn.1002-6819.2021.22.021 http://www.tcsae.org

Fang Hui, Wu Gang, Cheng Ruifeng, et al. Numerical simulation of ventilation efficiency with perforated air tubes for controlling canopy environment in a plant factory based on CFD[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(22): 188-193. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2021.22.021 http://www.tcsae.org

2021-06-30

2021-11-06

國家重點(diǎn)研發(fā)計(jì)劃戰(zhàn)略性科技創(chuàng)新合作專項(xiàng)（2020YFE0203600）

方慧，副研究員，研究方向?yàn)樵O(shè)施環(huán)境模擬方面。Email：fanghui@caas.cn

仝宇欣，博士，副研究員，研究方向?yàn)樵O(shè)施園藝方面。Email：tongyuxin@caas.cn

10.11975/j.issn.1002-6819.2021.22.021

S625.4

A

1002-6819(2021)-22-0188-06