胡俊超 馬昊 程博

摘 要： 為探索解決露地軟籽石榴種植越冬問題的新途徑，設(shè)計、搭建了一種簡易式非固定防寒設(shè)施，針對該防寒設(shè)施冬季日間集熱特性，利用ANSYS Fluent 建立了穩(wěn)態(tài)模擬模型，對其進行了驗證試驗，并使用該模型對防寒設(shè)施參數(shù)（空氣間層和通風孔直徑）進行優(yōu)化試驗。驗證試驗表明，華氏度相對誤差ERE<1‰，防寒設(shè)施內(nèi)監(jiān)測點溫度模擬值的決定系數(shù)R2 和均方根誤差ERMSE 分別為0.962 5 與1.19 K，表明模型可信度較高。優(yōu)化試驗表明，防寒設(shè)施空氣間層厚度最佳3 cm、通風孔直徑最佳14 cm。

關(guān)鍵詞：軟籽石榴；防寒設(shè)施；CFD 模擬

中圖分類號：S628文獻標識碼：A文章編號：2095-1795（2023）12-0054-06

DOI：10.19998/j.cnki.2095-1795.2023.12.010

0 引言

軟籽石榴樹作為從國外引進的食用性狀優(yōu)良的果樹種類，豐富國內(nèi)水果市場品類的同時，為石榴種植戶帶來了相當可觀的經(jīng)濟收益。但軟籽石榴果樹普遍抗寒性較差，在–10 °C 溫度下，半日便會出現(xiàn)凍害[1]。目前，已有大量報道表明，低溫為突尼斯軟籽石榴引種推廣的主要限制因素[2-4]。國內(nèi)研究人員從軟籽石榴的防寒設(shè)施越冬、抗寒措施等方面，進行了一些探索研究，但研究進展較為緩慢，軟籽石榴抗寒越冬栽培技術(shù)尚未成熟。

覆土、覆蓋保溫材料是目前國內(nèi)主流的果樹防寒措施，據(jù)此，一些研究將防寒措施應(yīng)用于軟籽石榴。覆蓋防寒在防寒設(shè)施軟籽石榴種植中的應(yīng)用研究較少，馮一峰等[5] 對冷棚內(nèi)不同防寒措施下軟籽石榴的生長環(huán)境、生長狀況進行了相關(guān)研究，并得出秸稈＋黑膜的防寒設(shè)施內(nèi)防寒措施對軟籽石榴越冬保溫、防霉作用最佳。在露地軟籽石榴栽培上的應(yīng)用研究多一些，李銀芳等[6-9] 研究了不同覆蓋物下新疆果樹（包括軟籽石榴）的越冬效果及成本，并對防寒措施下的栽培管理提出了建議。宋娟等[10] 對不同防寒處理下匍匐栽培突尼斯軟籽石榴微域溫度進行了研究，綜合得出鋼筋架+草苫+厚覆土（30 cm）對突尼斯軟籽石榴微域溫度提升最明顯。國內(nèi)研究人員對軟籽石榴防寒越冬問題的解決途徑較為局限，保溫覆蓋較為簡陋，難以實現(xiàn)標準化、機械化種植，并且這些措施對于露地栽培的軟籽石榴防寒越冬局限性強，如在覆蓋物下的狹小空間內(nèi)，容易造成濕度過高，果樹枝條易發(fā)生霉變等問題。

為拓寬軟籽石榴防寒越冬問題的解決思路，本研究以露地傾斜栽培的軟籽石榴為研究對象，設(shè)計了一種非固定的防寒設(shè)施，以軟籽石榴防寒越冬防寒設(shè)施的模擬優(yōu)化為研究方向，在實物的基礎(chǔ)上，建立三維模型。利用CFD 技術(shù)，進行防寒設(shè)施內(nèi)日間穩(wěn)態(tài)溫度分析。將模型模擬結(jié)果與實測數(shù)據(jù)進行可靠性檢驗，在達到可靠性前提下，對不同尺寸參數(shù)配置下的防寒設(shè)施進行模擬、分析，優(yōu)選出最佳的參數(shù)配置，為防寒設(shè)施性能改良提供理論支持。

1 材料與方法

1.1 研究對象

1.1.1 防寒設(shè)施結(jié)構(gòu)設(shè)計

以露地傾斜平地栽培（東西向種植，植株垂直高100 cm、株距100 cm）的軟籽石榴苗為設(shè)計對象。設(shè)計相鄰3 棵軟籽石榴植株共用一個防寒設(shè)施。防寒設(shè)施整體呈三棱柱形，東西長400 cm、高100 cm、底寬約85 cm，防寒設(shè)施剖面結(jié)構(gòu)如圖1 所示。東西兩側(cè)為三角形彩鋼板，北側(cè)為傾斜角度60°的矩形彩鋼板（采光集熱板），彩鋼板厚度均為10 cm。采光集熱板上設(shè)有4 個通風孔，直徑10 cm，位置如圖2 所示。

采光集熱板仰面噴涂黑色油漆，其上覆蓋0.03 cm厚PVC 薄膜（集熱罩），膜與板間撐起1 cm 厚空氣間層；防寒設(shè)施南側(cè)懸掛厚度0.5 cm 的草苫和0.03 cm的PVC 薄膜，形成封閉圍護結(jié)構(gòu)。2022 年10 月15 日防寒設(shè)施實物在石河子大學(xué)試驗地完成搭建，如圖3所示。

1.1.2 防寒設(shè)施原理

以防寒設(shè)施及其內(nèi)部所占空氣為研究對象，日間、夜間能量平衡方程如下。

式中 Qin——室內(nèi)熱量

Qsolar——防寒設(shè)施內(nèi)凈太陽輻射量

Qf——各壁面與外界的凈輻射量

Qv——冷風滲透量

Qc——各壁面對流散熱量

Qs——地中傳熱量

考慮到冬季落葉果樹對光照要求低，防寒設(shè)施以犧牲光照為條件，增加了防寒結(jié)構(gòu)的面積，以降低外界溫度對防寒設(shè)施內(nèi)部溫度的影響。

1.2.4 邊界條件

（1）輻射模型。防寒設(shè)施模型日間集熱過程涉及壁面間的輻射熱傳導(dǎo)及透明介質(zhì)的輻射問題，需要使用輻射模型。選取Do 模型作為輻射模型，模型允許用來解決各類輻射問題，包括半透明壁面的輻射問題和不透明面與面之間的輻射問題，普適性強，應(yīng)用廣泛[14]。日間的太陽輻射量采用太陽輻射追蹤模型，輸入石河子市的地理位置數(shù)據(jù)、模擬日日期和時間，晴天氣象因子輸入0.75，陰天氣象因子輸入0.45。

（2）材料物性及邊界參數(shù)設(shè)置。采光集熱板結(jié)構(gòu)仰側(cè)至背側(cè)依次為0.2 cm 鐵板、10 cm 苯板和0.2 cm鐵板，考慮到上層鐵板涂有吸光材料，具有集熱作用，將其依次簡化為0.2 cm 材料為鐵平面、10 cm 材料為苯板固體域；兩側(cè)三角形擋板直接簡化為10 cm 材料為苯板平面。集熱罩簡化為0.03 cm 的PVC 材料平面。南側(cè)柔性材料由外至內(nèi)為0.03 cm 的PVC 膜＋0.5 cm的草苫，簡化為0.5 cm 草苫材料平面。地面簡化為10 cm 材料為土壤的平面。模擬所用材料物理參數(shù)如表1 所示。

對防寒設(shè)施日間太陽輻射下的穩(wěn)態(tài)溫度場進行模擬，邊界條件設(shè)置如表2 所示。

1.2.5 初始條件

以模擬日日間防寒設(shè)施內(nèi)各監(jiān)測點溫度值處于較穩(wěn)定狀態(tài)時的外界環(huán)境溫度和太陽輻射量作為模擬的初始條件。

1.3 試驗設(shè)計

1.3.1 試驗地點及設(shè)備

防寒設(shè)施仿真實物建于石河子大學(xué)北苑新區(qū)試驗地，地理位置為東經(jīng)85.94°、北緯44.27°。測溫設(shè)備使用彭云物聯(lián)溫濕度計，類型編號S10A-2 300，共4 臺。該設(shè)備溫度測量精度±0.3 K，記錄間隔5 min。

1.3.2 驗證試驗

試驗設(shè)置3 個溫度監(jiān)測點：防寒設(shè)施東側(cè)的上、下通風孔中心點和防寒設(shè)施內(nèi)近中心點。上通風孔中心點（T1）對應(yīng)模型坐標（?12.49，78，180），下通風孔中心點（T2）對應(yīng)模型坐標（?47.7，17.3，180），近中心點（T3）對應(yīng)模型坐標（0，40，0），坐標單位cm，x 正方向為正東，y 正方向為正上，z 正方向為正北。此外，在防寒設(shè)施外高50 cm 處設(shè)置一個測溫設(shè)備，監(jiān)測外環(huán)境溫度。測試時間2020 年11 月16—24 日，采用2020 年11 月20 日（晴天）及2020 年11月23 日（陰天） 14： 30—16： 30 防寒設(shè)施內(nèi)各監(jiān)測點平均溫度進行穩(wěn)態(tài)模擬驗證。

1.3.3 優(yōu)化試驗

集熱墻通風口直徑和空氣間層厚度對其集熱性能有較大的影響[15]。防寒設(shè)施防寒保溫作用直接或者間接來自于采光集熱板的集熱效應(yīng)，為達到最佳的保溫防寒效果，需對防寒設(shè)施中的采光集熱板相關(guān)尺寸參數(shù)進行優(yōu)化。

設(shè)置1、3、5 和7 cm 厚的空氣間層厚度，每種間層厚度下設(shè)置6、10、14 和19 cm 4 種通風孔直徑。通過更改模型參數(shù)，以晴天的輻射條件，對16 組模型進行晴天日間流場模擬，獲得不同參數(shù)下防寒設(shè)施同一時間段下的穩(wěn)態(tài)溫度場，以防寒設(shè)施內(nèi)3 個監(jiān)測點模擬溫度為依據(jù)，對不同配置下晴天日間的集熱能力進行判斷，得出最佳集熱采光板參數(shù)配置。

2 結(jié)果與分析

2.1 驗證試驗結(jié)果與穩(wěn)態(tài)場分析

2.1.1 驗證試驗結(jié)果

選用相對誤差（ERE）、均方根誤差（ERMSE）和決定系數(shù)（R2）對模擬值和實測值之間的符合度進行統(tǒng)計分析?？紤]到攝氏度相對誤差在0 °C 附近時，分子過小并不能合理地反映誤差值，本研究選擇華氏度千分相對誤差。ERMSE 計算公式見參考文獻[16]、ERE 計算公式見參考文獻[17]。

防寒設(shè)施內(nèi)監(jiān)測點溫度模擬值與實測值ERE 經(jīng)計算整理如表3 所示，華氏溫度的絕對誤差≤1.79 K，千分相對誤差均<1‰。防寒設(shè)施內(nèi)監(jiān)測點溫度模擬結(jié)果如圖4 所示， 模擬溫度與測試溫度的決定系數(shù)R2 和ERMSE 分別為0.962 5 和1.19 K，由圖4 可知，模擬值與實測值接近1∶1 線，表明模擬方案的可靠性較高。

2.1.2 穩(wěn)態(tài)場分析

截取計算區(qū)域內(nèi)上、下通風口中心處的溫度場剖面云圖，可以看出日間穩(wěn)態(tài)時防寒設(shè)施最高溫度出現(xiàn)在集熱采光板仰面中上部，溫度可達301.5 K，最低溫度出現(xiàn)在草苫近地面，為271.08 K。采光集熱板固體域溫度仰側(cè)向背側(cè)遞減，中部向四周遞減。不包括空氣間層在內(nèi)的流體域，整體呈現(xiàn)溫度隨高度增加而增加的趨勢，上通風及以上區(qū)域流體溫度較高，并且上下部溫差最大在14 K 左右。模擬誤差可能主要來源于地中傳熱。

2.2 優(yōu)化試驗結(jié)果分析

各監(jiān)測點空氣溫度隨空氣間層厚度和通風直徑變化的趨勢如圖5 所示。

2.2.1 空氣間層厚度對防寒特性的影響

由圖5 可知，空氣間層厚度對防寒設(shè)施的集熱性能有較大影響，隨著空氣間層厚度增加，各監(jiān)測點空氣溫度均呈現(xiàn)先升高后降低的趨勢。中心、下通風孔處空氣溫度在空氣間層厚度3 cm 時達到最大值，此后呈現(xiàn)下降趨勢。上通風孔處空氣溫度變化稍有不同，對于通風孔直徑6 和10 cm 的采光集熱板，上通風孔處最高溫度出現(xiàn)在空氣間層厚度5 cm 時；而對于通風孔直徑14 和18 cm 的集熱采光板，上通風孔空氣溫度在空氣間層厚度7 cm 時，仍然呈現(xiàn)上升趨勢，可能與空氣間層厚度試驗組設(shè)置不足有關(guān)?？紤]到上、下通風孔處空氣溫度對防寒設(shè)施外環(huán)境變化反應(yīng)較為敏感，對于防寒設(shè)施內(nèi)平均空氣溫度變化代表性較差，防寒設(shè)施集熱性能評定主要參考中心空氣溫度變化規(guī)律，因此，可以得出空氣間層厚度3 cm 時，防寒設(shè)施的集熱性能較好，中心空氣溫度最高達到285.91 K。

2.2.2 通風孔直徑尺寸對防寒特性的影響

由圖5 可知，通風口直徑變化對各監(jiān)測點空氣溫度的影響差異較大，隨著通風孔直徑的增加，上通風孔空氣溫度呈現(xiàn)下降趨勢，下通風孔空氣溫度大致呈現(xiàn)平緩上升趨勢，中心空氣溫度大致呈現(xiàn)先升高再降低的趨勢。上通風孔空氣溫度呈現(xiàn)與下通風孔和中心空氣溫度相反規(guī)律，可能與通風孔直徑增大利于間層熱空氣與防寒設(shè)施內(nèi)部冷空氣流通有關(guān)。在空氣間層厚度3 和5 cm 時，中心空氣溫度最大值出現(xiàn)在通風孔直徑14 cm 時，而在空氣間層厚度1 和7 cm 時，中心空氣溫度在通風孔直徑18 cm 時，仍有增長趨勢，未達到最大值。就中心空氣溫度變化趨勢來看，對于空氣間層厚度3 和5 cm 下，防寒設(shè)施最佳通風孔直徑在14 cm 左右，對于空氣間層厚度1 和7 cm 時，防寒設(shè)施最佳通風直徑在18 cm 及以上。

3 結(jié)束語

將集熱墻原理拓展應(yīng)用至農(nóng)業(yè)防寒設(shè)施中，提出一種簡易型的非固定軟籽石榴防寒設(shè)施，建立并驗證了防寒設(shè)施冬季日間穩(wěn)態(tài)模擬模型，并對防寒設(shè)施參數(shù)進行了模擬優(yōu)化，指導(dǎo)該防寒設(shè)施防寒性能的優(yōu)化改良。

對于軟籽石榴的非適宜種植地區(qū)，低溫越冬對已建成和即將建成的露地軟籽石榴生產(chǎn)種植有極大的阻礙，嚴重影響軟籽石榴花、果數(shù)量和質(zhì)量。非固定軟籽石榴防寒設(shè)施能為露地軟籽石榴越冬提供較好的越冬條件，利于軟籽石榴生產(chǎn)與推廣種植，并能為北方果樹種植防寒越冬提供新的途徑。

模擬試驗表明，建立的穩(wěn)態(tài)模型具有較高的可信度，為優(yōu)化模擬試驗提供理論依據(jù)。優(yōu)化試驗表明，空氣間層厚度3 cm、通風孔直徑14 cm 時，防寒設(shè)施的集熱效果最佳，試驗得到的防寒設(shè)施內(nèi)各監(jiān)測點溫度變化趨勢與黑賞罡等[17] 的研究結(jié)果相近。

本研究僅對防寒設(shè)施進行了特定時段的穩(wěn)態(tài)模擬，不能在瞬態(tài)模擬下獲得防寒設(shè)施內(nèi)日間溫度隨時間的變化規(guī)律，但由于防寒設(shè)施日間集熱穩(wěn)態(tài)溫度場直接影響防寒設(shè)施日間防寒效果、間接影響防寒設(shè)施夜間防寒效果，穩(wěn)態(tài)模擬的使用簡化了模擬試驗步驟，利于之后的優(yōu)化模擬試驗快速獲得結(jié)論。此外，未對防寒設(shè)施夜間防寒性能及其適宜的參數(shù)進行研究，保溫性能相關(guān)研究有待完善。

參考文獻

薛華柏，曹尚銀，郭俊英，等．突尼斯軟籽石榴氣候區(qū)劃北限及次適宜區(qū)的防寒栽培[J]．中國果樹，2010（2）：63-64．

李敏．突尼斯軟籽石榴凍旱的發(fā)生與預(yù)防[D]．泰安：山東農(nóng)業(yè)大學(xué)，2013．

LI Min． The happening and prevention of the cold and drought inTunisia soft seeds pomegranate[D]． Tai'an： Shandong AgriculturalUniversity，2013．

李明婉，唐琳，李宗圈，等．河南省丘陵石榴主產(chǎn)區(qū)2009 年凍害調(diào)查[J]．河南農(nóng)業(yè)科學(xué)，2010（11）：106-108．

王宇翔，高小峰，雷夢瑤，等．南陽地區(qū)發(fā)展軟籽石榴的凍害風險分析[J]．南陽師范學(xué)院學(xué)報，2022，21（3）：45-50．

WANG Yuxiang，GAO Xiaofeng，LEI Mengyao，et al．Risk analysisof freezing injury in developing soft seed pomegranate in Nanyangarea[J]． Journal of Nanyang Normal University， 2022， 21（ 3） ： 45-50．

馮一峰，王艷，楊植，等．不同防寒措施對防寒設(shè)施軟籽石榴越冬性的影響[J]．北方園藝，2021（10）：53-58．

FENG Yifeng， WANG Yan， YANG Zhi， et al． Effects of differentcold prevention measures on cold hardiness of soft seed pomegranate[J]．Northern Horticulture，2021（10）：53-58．

李銀芳，潘伯榮，孫永強，等．不同覆蓋物對果樹越冬的保溫作用[J]．北方園藝，2012（18）：19-22．

LI Yinfang，PAN Borong，SUN Yongqiang，et al．The heat preservationeffects of different coverings for fruit trees' overwintering[J]．Northern Horticulture，2012（18）：19-22．

李銀芳，潘伯榮，孫永強，等．果樹越冬不同覆蓋方法的成本計算[J]．北方園藝，2012（14）：197-200．

LI Yinfang，PAN Borong，SUN Yongqiang，et al．The cost calculationof different coverings for fruit trees' overwintering[J]． NorthernHorticulture，2012（14）：197-200．

李銀芳．利用簡易防寒設(shè)施防止果樹凍害的栽培模式[J]．北方園藝，2010（1）：92-94．

LI Yinfang．Cultivation mode of fruit tree for preventing freezing injuryusing minimum greenhouse[J]．Northern Horticulture，2010（ 1）：92-94．

李銀芳． 寒冷地區(qū)雙層覆蓋式果樹越冬技術(shù)[J]． 北方園藝，2008（12）：17-20．

LI Yinfang．The double coverage hibernalisation technology of fruit treein the cold regions[J]．Northern Horticulture，2008（12）：17-20．

宋娟，胡曉靜，唐誠，等．防寒處理提升南疆突尼斯軟籽石榴微域溫度效應(yīng)的研究[J/OL]．果樹學(xué)報：1-17. https：//doi.org/10.13925/j.cnki.gsxb.20220643．

SONG Juan， HU Xiaojing， TANG Cheng， et al． Temperature increasingeffect of cold-prevention treatments in micro-environment ofTunisian soft seed pomegranate trees in South Xinjiang[J/OL]． Journalof Fruit Science：1-17. https：//doi.org/10.13925/j.cnki.gsxb.20220643．

佟國紅，王鐵良，白義奎，等．日光溫室土壤溫度分布邊際效應(yīng)的數(shù)值模擬[J]．北方園藝，2010（15）：65-68．

TONG Guohong， WANG Tieliang， BAI Yikui， et al． Numericalmodeling of marginal effect on soil temperature distribution in a chinesesolar greenhouse[J]．North Horticulture，2010（15）：65-68．

PASUT W， DE CARLI M． Evaluation of various CFD modellingstrategies in predicting airflow and temperature in a naturally ventilateddouble skin fa?ade[J]．Applied Thermal Engineering，2012，37：267-274．

李宗翔，顧潤紅，張曉明，等．基于RNG k-ε 湍流模型的3D 采空區(qū)瓦斯上浮貯移[J]．煤炭學(xué)報，2014，39（5）：880-885．

LI Zongxiang， GU Runhong， ZHANG Xiaoming， et al． Simulationof gas migration in 3D goaf based on RNG k-? turbulence model[J]．Journal of China Coal Society，2014，39（5）：880-885．

BAJC T，TODOROVIC` M N，SVORCAN J．CFD analyses for passivehouse with Trombe wall and impact to energy demand[J]．Energy &Buildings，2014，98：39-44．

楊昭， 徐曉麗． 特朗勃壁溫度場分析[J]． 工程熱物理學(xué)報，2006（4）：568-570．

YANG Zhao， XU Xiaoli． Analysis on the temperature field of trombewall[J]． Journal of Engineering Thermophysics， 2006（ 4） ： 568-570．

程陳，馮利平，薛慶禹，等．日光溫室黃瓜生長發(fā)育模擬模型[J]．應(yīng)用生態(tài)學(xué)報，2019，30（10）：3 491-3 500．

CHENG Chen，F(xiàn)ENG Liping，XUE Qingyu，et al．Simulation modelfor cucumber growth and development in sunlight greenhouse[J]． Chinese Journal of Applied Ecology， 2019， 30（ 10） ： 3 491-3 500．

黑賞罡，姜曙光，楊駿，等．特朗伯墻體冬季集熱性能的CFD模擬分析[J]．四川建筑科學(xué)研究，2018，44（3）：116-121．

HEI Shanggang，JIANG Shuguang，YANG Jun，et al．Simulation ofthe thermal performance of a Trombe wall in winter using computationalfluid dynamics[J]．Sichuan Building Science Research，2018，44（3）：116-121．