摘要：基于計算流體動力學原理，建立了換熱器流體流動和傳熱過程的控制方程組，應用算子分裂法，將換熱器控制方程分裂為基本的對流、擴散、Stokes 1\'7題；采用8～20節點三維等參單元，推導出有限元離散化方程組。編制了流體流動和傳熱的數值仿真系統程序，應用仿真系統進行傳熱數值模擬計算和實驗驗證，通過對比兩者的Nu與Re的關系，可以看出仿真系統的計算值和實驗測量值非常接近。實驗證明，采用編制的折流桿換熱器數值仿真系統進行傳熱模擬計算是切實可行的。

　　關鍵詞：折流桿換熱器CFD 傳熱 數值模擬 仿真系統

　　引 言

　　折流桿換熱器以桿式支撐替代原弓形擋板，具有抗振、高效、低壓降等優點。與傳統的折流板管殼式換熱器相比，折流桿換熱器內部結構有較大變化，其殼程內部采用折流桿組成的折流柵取代折流板作管間支撐物，使殼程流體由橫向流動變為平行流動，不僅大大減少了傳熱死區，而且大幅度減少了流體因多次反復折流而損失的殼程壓降。與現行管殼式換熱器相比，折流桿換熱器的總傳熱系數提高了35％ 以上，殼側阻力降減小了50％以上。計算流體動力學(CFD)是進行傳熱、傳質、動量傳遞及燃燒、多相流和化學反應等研究的核心和重要技術，具有成本低和能模擬較復雜過程等優點¨0 。在給定的參數下用計算機對現象進行一次數值模擬相當于進行一次數值實驗，通過CFD計算結果就可以進行評估，并且可以比較快捷、準確及直觀地反映出流體在換熱器中流動的過程，如速度場、壓力場、溫度場等的分布，可以得到用實驗方法無法得到的流場和溫度場分布的細觀信息，易于分析各種結構對傳熱過程的影響。因此，對新型換熱器的研制具有深遠的影響。

    筆者在此闡述了采用CFD技術的折流桿換熱器數值仿真系統的編制思路和功能，并對其正確性進行了考核，對新型換熱器及其強化傳熱元件的研究與開發方面具有重要意義。

    換熱器CFD數值算法建模目前在計算流體力學中，采用有限元法計算三維流體問題，都是采用8節點三維單元，該單元形式簡單，但精度低，誤差大，而且對復雜曲邊邊界適應性差。筆者針對一般Galerkin有限元法的不足，為改進有限元法數值計算的收斂性和提高計算精度，建立有限元法改進的離散格式，采取了如下一些改進措施：為改進收斂性，應用算子分裂法，將換熱器控制方程分裂為基本的對流、擴散、Stokes問題；為提高計算精度，提出并推導出二階時間精度離散格式；為提高8節點 維等參單元的精度和復雜邊界的適應性，采用8～20節點三維等參單元，并推導出其空間離散格式；為進一步提高數值仿真的收斂性，利用最小二乘法原理，推導出一個較規范統一的耗散權函數。

    1．控制方程

    設換熱器橫截面為X一Y平面，沿換熱器長度方向為z軸。流體流動處于層流狀態，殼程流體控制方程有：

    連續方程

                     

                     

                     

                     

                     

                     

                     

                      

    數值仿真系統的編制

    數值仿真系統主要由以下3個部分組成：基本前處理模塊(FEMPREP)、有限元流動與傳熱計算模塊(FEMFLOW)、計算結果后處理模塊(FEM—POST)。仿真專家系統的程序流程圖如圖1所示。

                   

    程序采用Visual C++、Fo~ran90混合語言編寫，并編譯成可執行文件，其中Fo~ran90采用Fo~ran Power Station 4．0編譯。在程序中，采用了模塊化設計、動態內存管理技術、內外存數據交換等技術。仿真專家系統各模塊的功能概述如下。

    1．基本前處理模塊(FEMPREP)

    該模塊首先需要定義單元類型和材料屬性。單元類型是8—2O節點三維等參單元。材料屬性主要有管程、殼程流體密度、比熱、導熱系數、粘度等各個參數；換熱管材料的導熱系數。其次生成模型的幾何體，如生成換熱器橫截面( —Y平面)的幾何結構。第三，劃分有限元網格。對于換熱器，根據圖2所示的網格劃分原則，計算各節點的位置及進行單元網格劃分(二維)，生成二維有限元網格，沿 方向延伸(拉伸)，即可得到三維有限元網格，并形成網格數據文件。第四，施加載荷或邊界條件。例如初始條件、邊界條件等。第五，定義分析類型和分析選項。例如根據流體流動狀況，人工選擇是層流計算還是湍流計算；是否需要進行能量(溫度)計算；迭代方法選擇(牛頓一拉斐遜迭代法和預處理的共軛梯度法)；收斂準則；最大迭代次數等多個選項。最后形成用于計算需要的數據文件。基本前處理模塊的應用，可有效降低有限元分析的數據輸入量，提高效率。

                    

    2．有限元流動與傳熱計算模塊(FEMFLOW)

    有限元流動與傳熱計算模塊是數值仿真系統的核心和關鍵，計算成敗、計算效率、計算結果的精確性以及物理真實性、可靠性主要體現在該模塊。該模塊以式(5)、(6)、(10)、(1 1)等為基本有限元離散化方程組，通過建立單元剛度矩陣，組建總剛度矩陣，采用增量牛頓一拉斐遜迭代法、增量預處理的共軛梯度法進行求解。該模塊的程序流程圖如圖3所示。INPUT過程模塊打開在基本前處理模塊中所形成的數據文件，獲取計算所需的全部數據；ADDRES過程模塊根據模型單元和節點數據計算總剛度矩陣對角元的存儲地址以及剛度矩陣的元素總長；FLOW 3D過程模塊調用8—2O節點三維等參單元，以建立單元剛度矩陣，并組集總體剛度矩陣；ALGOR過程模塊是牛頓一拉斐遜迭代法、預處理的共軛梯度法的求解算法；CALSOL過程模塊采用Lu分解法求解總剛度矩陣；EQUIT過程模塊進行平衡迭代計算。最后可計算得到各點的速度、壓力、溫度值。

                

    3．計算結果后處理模塊(FEMPOST)

    該模塊采用粒子跟蹤技術、等值面(等值線)技術等將有限元流動與傳熱計算模塊得到的計算結果轉化為等值線圖、矢量圖或流線圖表示，以較直觀地顯示出流體流動和傳熱的數值模擬計算結果。數值仿真系統的驗證

   為考核仿真系統計算折流桿換熱器的能力，采用如下的折流桿換熱器：筒體內徑Di為145mm，管束長徑比(管束長度與殼體直徑之比)為Li／Di= 13，換熱管為Φl4mm×2mm，37根，管間距為19mm；換熱管采用正方形排列；折流柵間距Lb=

50mm，折流圈規格為Φ143 mm×Φ134 mm，折流桿布置為單排管間布桿，折流桿直徑為5mm，殼程進口接管為Φ76 mm×4mmm，殼程進口接管與管板的距離為100 mm。管程通入熱流體為105℃飽和水蒸氣，殼程通入20℃空氣。應用數值仿真系統采用層流選項計算出Re≤2 300時Nu值與值Re的關系已繪于圖4，實驗獲得的Nu與Re的關系也繪于圖4中。

            

    從圖4中可以看出，仿真系統的計算值和實驗測量值非常接近，兩者的lgNu與lgRe關系呈良好線性，且斜率均為0．6。仿真系統計算的實驗考核證明，采用仿真系統進行數值模擬分析是適用的，用于折流桿換熱器的計算是可行的。

    結 論

    筆者探討了基于CFD技術的折流桿換熱器數值模擬實現方式，研究了其有限元離散化方法，推導出有限元離散化方程組。闡述了折流桿換熱器流動和傳熱數值仿真系統的編制思路以及關鍵模塊的程序流程圖。應用仿真系統進行傳熱數值模擬計算，通過實驗進行了驗證，證明應用折流桿換熱器數值仿真系統進行模擬計算是切實可行的。折流桿換熱器仿真系統的研究具有很重要的工程應用價值，為新型換熱器的研究、設計開發提供了有效的手段和工具。