事實上，日光溫室的圍護結構內部的傳熱過程并非穩態，為了更為客觀的掌握溫室內部的熱環境，就要進行墻體非穩態傳熱分析，進而了解墻體的熱工性能及其熱傳導過程，進而改善其隔熱保溫措施，這對日光溫室能耗分析來說是關鍵的一部分，有助于節能型日光溫室設計和運用的發展。
一、墻體非穩態傳熱
墻體非穩態傳熱過程、墻體內部的溫度分布以及熱流量隨著時間的變化規律的求解，實質上就是求解以下兩個偏微分方程：上述兩個偏微分方程的解法主要有以下兩種：積分變換法和有限差分法。本文所關心的問題是墻體表面發生的溫度和熱流量變化，至于墻體內每個部位的具體溫度變化則不必細究。日光溫室圍護結構的整個傳熱過程溫度變化范圍較小，可以將墻體導熱系數又和導溫系數Q當做是常數。因此，選擇使用積分變換法來計算該常系數線性偏微分方程組定解問題。按照不同的邊界條件分解方式，積分變換法可分為諧波反應法和瞬態反應法。由于太陽輻射等外部條件時刻發生變化，引起改變的不確定的因素時刻存在，難以用規律的正弦函數來描述，所以諧波反應法不可能作為持續性的能耗分析工具0瞬態反應法又包括Z傳遞函數法和反應系數法，兩種方法都是把連續擾量離散成按時間序列進行分布的單元擾量。板壁的非穩態傳熱中，常見等腰三角波作為單元擾量，即要求獲得板壁熱力系統的Z傳遞系數或者反應系數，而兩者具有相似的含義和求解方法。以下就分別對兩種方法進行簡單的介紹。
(1)反應系數法
反應系數法在求解問題時基本思路是，把隨著時間連續變化的擾量曲線離散成為按時間序列進行分布的單元擾量；然后再求出墻體的熱力系統對于單位單元擾量的反應系數；最后，運用已求得的反應系數疊加積分后算出單位面積板壁的圍護結構逐時傳熱得熱量。Mitalas和Stephenson提出并且詳細介紹了反應系數法在負荷計算中起到的應用。墻體的反應系數序列YG)的定義是墻體對于其單位等腰溫度的三角波輸入熱流所具有的輸出值的等時間序列。它能夠描述墻體受到室外溫度擾量的動態響應過程，通過三角波的疊加來逼近室外空氣綜合溫度變化情況，從而獲得墻體熱力系統對于任意室外擾量所產生的響應。反應系數法適用于任何一種擾量，它有別于諧波法的一個重要的特點就是，反應系數法并不要求邊界條件具有周期性，因此它可以用于全年性的動態模擬。為了獲得反應系數，需要對墻體傳熱或者吸熱傳遞函數的分母B(s)求根，也就是求出s傳遞函數的極值點，然后運用Heaviside展開式求解得到最終的結果，通常稱為直接求根法。但B(s)有無窮多個負實根(-s.)，所以求出所有的根是不現實的。根小到一定程度時，其對反應系數產生的影響非常微小，一般只需求得一定數量的根值，就可以保證反應系數所具有的精度負荷要求。通過計算和比較，我們發現，計算Y(0}和Y(1)時，取傳遞函數法的理論基礎是由當前的輸入值和前一時刻的輸入值和輸出值共同確定。
和反應系數相同，傳遞函數也能夠描述墻體動態熱特性，但是所需求的系數項比起反應系數來說少得多，這就使得計算所耗費的時間和計算機存貯空間大大的減小。系統的z傳遞函數，也即脈沖傳遞函數，其實是系統單位脈沖響應的z變換。但這里所說到的板壁圍護結構z傳遞函數，是指板壁熱力系統的單位等腰三角脈沖產生的響應的z變換，相當于系統中設置了單位等腰三角波保持器。所以，板壁圍護結構傳熱和吸熱的z傳遞函數，指的就是板壁吸熱和傳熱反應系數序列的z變換。相比于其他的方法，傳遞函數法最為重要的性質就是其輸入輸出的值是等時間間隔序列值。因此，作為輸入值，太陽輻射強度和室外溫度的逐時記錄數據，只需要做很少的預處理。此方法主要的局限性是它只能夠用于線性時不變的系統；而傳遞函數則可以通過直接求根法來求解墻體的吸熱或傳熱z傳遞函數來獲得，盡管直接求根法不可避免的失根現象在這里仍然存在；z傳遞函數作為反應系數的近似形式，必然會導致一定的誤差。綜合上述三種常用墻體非穩態傳熱的計算方法各自的優缺點，本文決定運用反應系數法來進行墻體非穩態傳熱的計算。
二、反應系數計算法
(1)擾量的離散
用等腰三角波來進行擾量的離散，能夠將室外氣候的變化規律較為精準的表達。
三、定義非穩態邊界條件選取8月30日早9：00^-16：30作為非穩態計算的時間段。計算域的確定、網格的劃分、控制方程及模擬模型的相關設置都與穩態模擬相同，主要邊界條件的設置為：以上邊界中除了需要用自定義函數(UDF)來定義的非穩態邊界條件而外，其余邊界條件的確定均與穩態數值模擬時相同，這里不再重復了。下面著重的介紹非穩態邊界條件中入口溫度、膜以及土壤邊界條件的設置情況。因為早9：00—16：30這個時段入口溫度、膜、土壤溫度和后墻風口里來自室外的新風溫度波動度都比較大，膜溫隨著太陽輻射強度的變化也一直在改變，所以定義為定溫或者絕熱邊界都與實際的情況具有較大的差距。并且這些溫度都具有一定的變化規律，早晚較低而中午較高，即與時間t存在密切的關系。上述五個回歸程式的相關系數分別為0.9833，0.9633，0.9843，0.9799，0.9044，說明所測入口溫度、膜、土壤溫度、后墻溫度、后墻風口溫度與時間呈顯著相關，因此可以用上述五式表達入口溫度、膜、土壤溫度、后墻溫度、后墻風口溫度隨時間的變化規律。將土述五個回歸方程式使用c語言源代碼表達出來，建立一個c語言源代碼文件，編譯UDF函數并連接到FLUENT，在FLUENT中執行UDF，進行計算，詳見附錄。
三、三維穩態計算結果與分析
設置完非穩態邊界條件以后，仍然取1e-3作為計算動量方程、連續性方程、k及￡方程時的收斂精度標準，輻射和能量方程的計算收斂精度標準為1e-6。迭代計算的步長值設為60s，每個步長內最大迭代次數取80次，每問隔30min自動保存當前時刻的計算結果。圖4-7給出了9：30，13：30與16：30溫度場的計算結果，分別選Y=1800mm(對準風口)、Y=6300mm(兩風口之間)兩個截的溫度值來進行比較分析。圖4-7反映出了溫室兩個具有代表性的剖面在不同時刻呈現出的溫度分布云圖。
本次非穩態模擬從初始時刻9：30開始到16：30結束，以上分別選擇了上午9：30(t=1.800e+03)、中午13：30(t=1.6200e+03)、下午16：30(t=2.7000e+03)三個時刻的溫度分布云圖來進行分析。由圖中能看出：從9：30到13：30溫室內溫度逐步升高，薄膜附近的室氣溫度變化幅度較大；到16：30時，室內空氣較13：30已經有了明顯的降低，這與實測結果吻合，在第二章中我們看到，16：00過后室外溫度開始呈下降趨勢。模擬的結果進一步驗證了太陽輻射是導致溫室內空氣溫度上升的主要原因。比較不同時刻的溫度分布云圖能夠看出，不同時刻室內空氣溫度沿高度Z方向的布有所不同，但基本規律相似，最高溫都出現在溫室頂部，這是由于空氣受到熱浮力往高處上升，并且通過前后風口進入室內的新風在室內吸取熱量后經風壓作用由頂部窗口排出。與測試數據相比較可知，整個計算空間的模擬得到的溫度分布趨勢和實測數據的規律基本相同。
四、實驗數據驗證非穩態模擬結果
在上午9：30的模擬結果中取出試驗測點處的溫度值與實測值做比較，如圖4-8所示，從圖4-8中可以看出，三維非穩態的模擬結果能夠準確的反應出溫室內不同點的溫度變化趨勢。模擬的結果整體比試驗數據偏小，最大誤差為3.04℃，位于下端靠近土壤表面處，最小誤差僅為0.001℃。造成模擬值偏小于實測值的原因，主要是室外氣候的隨機以及多變性在時刻改變自然通風口的空氣流動方向、速度大小、風量和溫度等，使得在本次模擬計算中，邊界條件很難精確的符合實際情況。
從圖中我們可以看到，雖然模擬結果顯示的溫室內空氣溫度分布與實測值在數值上存在一定差異，但除了邊界點外，最大誤差僅在2.2℃以內；并且空間分布的總體趨勢是一致的，所以本次對日光溫室熱環境的三維非穩態數值模擬是成功的。相比于穩態數值計算結果，非穩態數值模擬結果與實測值的吻合程度更高。所以，得出結論：本章節建立的三維非穩態CFD模型及包括的邊界條件有效，利用fluent來進行夏季自然通風條件下日光溫室的熱環境三維非穩態數值計算是可行的。http://m.gzsfwl.com
