兩種通風模式下的畜禽養(yǎng)殖場所內(nèi)環(huán)境場模擬研究

摘 要：農(nóng)村禽畜的生存環(huán)境對畜禽健康和生產(chǎn)性能的影響很大。其生存環(huán)境一般由空氣的溫度、濕度、氣流速度和太陽輻射等溫熱因素綜合而成。畜禽養(yǎng)殖場所內(nèi)的空氣流動狀況對畜禽生存環(huán)境有重大影響，并且直接影響其他因素作用的發(fā)揮。該文利用計算流體力學的方法分析了在風機進口走廊出口條件下畜禽養(yǎng)殖場所內(nèi)的速度場和溫度場的分布情況，并著重對比了兩種通風模式下的環(huán)境場內(nèi)的狀況。研究結(jié)果將為禽畜養(yǎng)殖場所的通風設計提供理論依據(jù)和指導。

關鍵詞：計算流體動力學 畜禽 養(yǎng)殖場所 環(huán)境模擬 兩種通風模式

中圖分類號：U463.65+1 文獻標識碼：A 文章編號：1674-098X（2012）12（c）-00-03

影響農(nóng)村畜禽健康和生產(chǎn)性能的因素很多，既有品種、飼料和防疫問題等因素，也有生存環(huán)境問題等因素。當品種、飼料和防疫問題基本解決以后，禽畜的生存環(huán)境對畜禽健康和生產(chǎn)性能將起到?jīng)Q定作用。禽畜的生存環(huán)境一般由空氣的溫度、濕度、氣流（風）速度和太陽輻射等溫熱因素綜合而成。文獻[1]指出當環(huán)境溫度達到32 ℃以上時，雞會發(fā)育緩慢，產(chǎn)蛋率下降，蛋殼變薄易碎，甚至死亡。文獻[2]指出氣流對家禽健康的影響，主要出現(xiàn)在寒冷的環(huán)境下，低溫高速氣流進入畜禽養(yǎng)殖場所，易于使畜禽發(fā)生關節(jié)炎、神經(jīng)炎、肌肉炎等，甚至引起凍傷。在以上眾多因素中，通風換氣是第一要素。文獻[3]指出通風換氣不僅可以在高溫的夏季通過加大氣流促使畜禽的散熱使其感到舒適，以緩解高溫對家禽的不良影響，還可以防止養(yǎng)殖場所內(nèi)潮濕，保障場所內(nèi)空氣清新，改善畜禽養(yǎng)殖場所內(nèi)的環(huán)境質(zhì)量。畜禽養(yǎng)殖場所內(nèi)通常采用機械通風方式進行場所內(nèi)的通風換氣。因此，如何合理的為畜禽養(yǎng)殖場所設計機械通風系統(tǒng)，如何采用合理的方法預測畜禽養(yǎng)殖場所內(nèi)的環(huán)境狀態(tài)，保證畜禽養(yǎng)殖場所內(nèi)的環(huán)境狀況滿足畜禽生產(chǎn)和健康的要求，就顯得非常重要。

計算流體動力學是通過計算機數(shù)值計算和圖像顯示，對包含有流體流動和熱傳導等相關物理現(xiàn)象的系統(tǒng)所做的分析，它可以得到在某種條件下的流場各個位置上的物理量，以及這些物理量隨時間的變化。該文采用計算流體動力學的方法對兩種通風模式下養(yǎng)豬場所內(nèi)的流場進行模擬，利用得到的氣流場以及溫度場等數(shù)據(jù)評價養(yǎng)殖場所內(nèi)通風條件是否符合要求，以期望進一步對通風系統(tǒng)進行優(yōu)化。

1 畜禽養(yǎng)殖場所物理模型

選擇桐廬農(nóng)村某種豬場的養(yǎng)殖場所為研究對象，其物理模型尺寸及邊界條件均為實際測量參數(shù)。種豬養(yǎng)殖場所內(nèi)空間尺寸為：5 m（長）×4 m（寬）×4.6 m（高），其中走廊空間尺寸為1 m（長）×4 m（寬）×4.6 m（高）。由于豬欄為柵欄形式，欄高0.8 m，柵欄之間縫隙很大。在一側(cè)墻上裝有1臺風機，風機直徑為0.85 m，5葉片，輪轂直徑為0.61 m，風機位于墻面的正中心。豬舍底部設有1條排尿溝，截面尺寸為0.1 m×0.1 m。豬舍兩側(cè)墻上的窗戶尺寸均為0.1 m×0.1 m。兩種通風模式下養(yǎng)殖場所內(nèi)的立體布置圖如圖1所示。

圖1 兩種通風模式下養(yǎng)殖場所內(nèi)立體圖

在對養(yǎng)殖場所進行建模時，不可能將場所內(nèi)的每一件物件都進行非常準確的建模，有的應該省略，有點應該適當簡化，其目的是為了降低建模和網(wǎng)格劃分難度，節(jié)省計算時間。故對該養(yǎng)殖場所做出如下省略及簡化：（1）由于豬舍內(nèi)柵欄較低且縫隙很大，故在建模時省略。（2）豬體是豬舍內(nèi)的主要熱源，對豬體進行合理的處理關系到養(yǎng)殖場所內(nèi)氣流場和溫度場模擬的準確性。由于豬體是不規(guī)則形狀，且在場所內(nèi)為運動狀態(tài)，故該文直接假定豬舍內(nèi)的兩面墻為發(fā)熱源。（3）由于養(yǎng)殖場所內(nèi)的窗戶與走廊相比可以近似忽略不計，故在建模時忽略了窗戶的

影響。

2 模型的網(wǎng)格劃分

物理模型建立完成以后，需要對模型進行網(wǎng)格劃分。整個流體區(qū)域被劃分為5個區(qū)域，進口區(qū)、出口區(qū)、風扇區(qū)、豬舍區(qū)、走廊區(qū)。該文假定風機入口為進口區(qū)，走廊一端出口為出口區(qū)。在風扇區(qū)采用四面體Tgrid網(wǎng)格，在其他區(qū)域采用六面體Cooper網(wǎng)格，在場所內(nèi)四壁建立Inflated Boundary，這些表面上的Maximum Thickness設置為0.01 m。整個模型的網(wǎng)格總數(shù)在220萬以上。

3 計算方法及邊界條件

考慮到養(yǎng)殖場所內(nèi)機械通風時，場所內(nèi)氣流狀態(tài)保持穩(wěn)定狀態(tài)，故該文選擇定常計算進行模擬。定常模擬主要方程選用三維雷諾平均守恒型Navier-Stokes 方程。

模擬計算采用不可壓分離式求解器（Segregated）隱式方案，采用RNG k-ε模型進行湍流模擬。在豬舍墻壁周圍，粘性流體的速度滿足無滑移條件，即相對壁面速度為零。壁面設置為豬的體溫，豬體壁面溫度為38 ℃，發(fā)熱方式為輻射發(fā)熱。計算域的入口假定為風機入口，設定為流量入口。該文選定流量為0.03 kg/s。出口假定為走廊一端口，設定為壓力出口。該文選定出口壓力為大氣壓。由于豬舍的流體均為空氣，故設置為理想氣體。

4 結(jié)果分析

4.1 速度場分析

圖3顯示了兩種通風模式下養(yǎng)殖場所流場區(qū)域子午面的流線圖，從圖中可以看出，空氣由于風機的負壓作用進入豬舍，風機橫向模式下在豬舍靠近入口墻壁的兩個較大區(qū)域形成旋渦，而后再流向走廊一端出口。風機縱向模式下在豬舍整個流場區(qū)域有很多漩渦，氣流脈動更加劇烈，通風效果更加明顯。有利于種豬的散熱，從而有利于健康及生產(chǎn)。圖4顯示了兩種通風模式下養(yǎng)殖場所流場區(qū)域子午面的速度云圖，從圖中可以看出，高速氣流在風機縱向通風模式下比在橫向通風模式下分布更廣，橫向通風模式下的高速氣流主要分布在風機出口附近，而縱向通風模式下的高速氣流主要分布在風機出口以及豬舍的中心區(qū)域，最大風速可達到6.09 m/s。

4.2 溫度場分析

圖5顯示了兩種通風模式下養(yǎng)殖場所流場區(qū)域子午面的溫度場云圖，從圖中可以看出，除了兩側(cè)墻壁設定溫度為38 ℃以外，其他區(qū)域溫度都要低于38 ℃，且靠近風機的流場區(qū)域的溫度場明顯低于其他區(qū)域，最低溫度為30 ℃，這是由于風機的散熱所導致的。從靠近墻壁截面的溫度場也能得出相似的結(jié)論。橫向通風模式下風機進口附近的溫度明顯低于其他區(qū)域，且分布有明顯的梯度，而縱向通風模式下養(yǎng)殖場所內(nèi)的溫度明顯低于橫向通風模式。

5 結(jié)語

畜禽養(yǎng)殖場所內(nèi)的環(huán)境對畜禽的健康以及生產(chǎn)的影響很大，而通風換氣對于養(yǎng)殖場所內(nèi)的環(huán)境最為關鍵。該文通過對某種豬養(yǎng)殖場所內(nèi)的環(huán)境進行數(shù)值模擬，得到以下結(jié)論：

（1）在養(yǎng)殖場所內(nèi)，橫向通風模式下的氣流在靠近風機兩側(cè)的墻壁附近有明顯的旋渦，而縱向通風模式下的氣流在整個場所內(nèi)都可以看到很多明顯的漩渦。

（2）高速氣流在風機縱向通風模式下比在橫向通風模式下分布更廣，橫向通風模式下的高速氣流主要分布在風機出口附近，而縱向通風模式下的高速氣流主要分布在風機出口以及豬舍的中心區(qū)域，最大風速可達到6.09m/s。

（3）除了兩側(cè)墻壁設定溫度為38 ℃以外，其他區(qū)域溫度都要低于38 ℃，且靠近風機的流場區(qū)域的溫度場明顯低于其他區(qū)域，最低溫度為30 ℃。橫向通風模式下風機進口附近的溫度明顯低于其他區(qū)域，且分布有明顯的梯度，而縱向通風模式下養(yǎng)殖場所內(nèi)的溫度明顯低于橫向通風

模式。

綜合以上結(jié)論，由此可判斷縱向通風模式更適合該養(yǎng)殖場所的通風條件。

參考文獻

[1] 韓華，劉繼軍，馬宗虎.計算流體力學在畜舍中的研究現(xiàn)狀與進展[J].中國畜牧雜志，2007，43（11）：38-40.

[2] 佟國紅，李保明.用CFD方法模擬日光溫室溫度環(huán)境初探[J].農(nóng)業(yè)工程學報，2007，23（7）：178-185.

[3] 佟國紅，張國強.豬舍內(nèi)氣流變化的模擬研究[J].沈陽農(nóng)業(yè)大學學報，2007，38（3）：379-382.