液態(tài)金屬湍流換熱的直接數(shù)值模擬研究.pdf

1、液態(tài)金屬在新一代核能技術(shù)領(lǐng)域中具有非常重要的應(yīng)用，但低普朗特數(shù)Pr的物性使得它們湍流換熱的表現(xiàn)與空氣、水大不相同，因此需要深入研究其特性。液態(tài)金屬湍流換熱的實驗研究困難，提供的數(shù)據(jù)主要是換熱系數(shù)、壓降等宏觀量，無法給出精確的湍流場信息，因此需要借助數(shù)值手段。而數(shù)值方法中工程應(yīng)用最廣泛的雷諾平均模擬RANS方法使用的湍流熱通量模型是針對空氣或水的，研究表明這不適用于低Pr數(shù)流體研究，無法提供準(zhǔn)確溫度場信息。而直接數(shù)值模擬DNS方法，不采用

2、任何模型并且可以提供精確的瞬時三維場數(shù)據(jù)，所以成為本文研究液態(tài)金屬湍流換熱特性的主要方法。
　　本文主要采用DNS方法開展了以鉛鉍合金LBE為代表的液態(tài)金屬的湍流換熱機理研究。在圓管中重點分析了湍流混合對流研究中浮力輔助流下的換熱機理;在后臺階中研究了不同Pr數(shù)對強迫對流換熱的影響，以及比較了浮力效應(yīng)在后臺階浮力輔助流和抵抗流中的差異。
　　圓管湍流混合對流研究中發(fā)現(xiàn)輔助流下近壁面各量及換熱隨浮力增大呈現(xiàn)先減小后增強的變化，

3、這主要是由于浮力對平均場的改變導(dǎo)致湍流發(fā)生變化。浮力較弱時，壁面熱流體加速抑制了湍流邊界層中的噴射-掃掠事件，湍流被削弱，甚至出現(xiàn)層流化;隨著浮力的進(jìn)一步增強，流向速度峰值逐漸接近加熱壁面，徑向分布呈現(xiàn)M型，此時剪切應(yīng)力反向增大，湍動能產(chǎn)生項增大，湍流增強。相比于空氣湍流混合對流，浮力對液態(tài)金屬流動的影響更顯著。
　　后臺階強迫對流顯示低Pr數(shù)流體的努賽爾數(shù)Nu極大值點與再附著點重合;Pr數(shù)增大，溫度脈動和換熱都會增強，它們的極大

4、值位置都會向上游移動;剪切層帶來的高強度渦旋沖擊是再附區(qū)換熱增強的主要原因;回流區(qū)內(nèi)湍流熱輸運的作用比較顯著，即使是對于低Pr數(shù)流體?；謴?fù)區(qū)內(nèi)低Pr數(shù)流體的湍流普朗特數(shù)Prt分布不均勻，這給液態(tài)金屬后臺階湍流換熱的RANS模擬帶來挑戰(zhàn)。
　　浮力對液態(tài)金屬后臺階湍流換熱有重要影響。輔助流下浮力加速熱流體，使得二次渦增大并超過主渦，壁面摩擦系數(shù)增大;回流區(qū)縮小，剪切層的縮短使湍流強度整體減弱;主渦向上游攜帶的熱量減少因而Nu數(shù)增強;

5、理查森數(shù)Ri=0.2時可見浮力的作用開始使近壁面湍流增強。抵抗流下很小的Ri就能使回流區(qū)迅速擴大，主渦占據(jù)整個回流區(qū)，二次渦幾乎消失，壁面摩擦系數(shù)減小;剪切層拉長使得湍流強度增強;主渦向上游攜帶熱量增加因而回流區(qū)內(nèi)Nu數(shù)減弱。
　　此外，基于本文后臺階DNS數(shù)據(jù)和棒束子通道換熱實驗關(guān)系式，本文還對RANS方法湍流熱通量封閉模型進(jìn)行改進(jìn)，采用Prt模型取代常數(shù)0.85。對不同Prt模型進(jìn)行比較分析后發(fā)現(xiàn)，基于當(dāng)?shù)亓鲌鰠?shù)的Kays模