高等流體力學課程論文

1、研究生課程考核試卷（適用于課程論文、提交報告）科 目： 高等流體力學 教 師： 葉建 姓 名： 楊偉，陸祺 ，張鵬青 學 號： 20131002051，20131002021t ，20131002026t 專 業(yè)： 動力工程及工程熱物理，核科學與技術, 動力工程及工程熱物理 類 別： 學術，學術 ,學術 上課時間： 2013 年9月至2013年 11 月 考 生 成 績：卷面成績平時成績課程綜合成績閱卷評語： 閱卷教師 (簽名) 氣冷堆燃料元件表面流場及溫度場模擬 1.本文主要研究內容國內外已經有許多學者利用計算流體力學（CFD）對模塊式高溫氣冷堆球床內的流場、溫度場做過大量研究，并取得了不少

2、成果。這方面的研究涉及了燃料元件與壁面或燃料元件之間接觸和間隙情況下的傳熱和流動的特征，孔隙率對溫度分布均勻性的影響等等。本文的主要研究內容是數(shù)值計算獲得堆芯元件表面冷卻劑氦氣的速度及溫度分布規(guī)律，并研究不同堆積方式對分布規(guī)律的影響。2.球床內流動與換熱數(shù)學模型實際上，在球床對內，燃料元件的分布是隨機不規(guī)則分布的，有些部分可能是正四棱錐堆積，有些部分是立方體形式堆積；當然，還有可能是其它的堆積方式。本文主要研究的是正四棱錐堆積于立方體堆積時對溫度、速度分布的影響。因此，建立的物理模型由立方體和正四棱錐兩種，采用的計算軟件為Fluent。2.1 物理模型在球床對內，球形燃料元件的不同堆積方式對氦

3、氣的流動及其溫度分布會產生的影響不同。因此，我們需要建立模型模擬以下情況：（1）球床內溫度、壓力和速度分布規(guī)律；（2）堆積方式變化時溫度分布和流動的改變?；谶@兩方面，我們需要建立不同堆積方式的模型。圖2.1 正四棱錐堆積圖2.2和圖2.3分別為立方體堆積模型和正四棱錐堆積模型，其中可以明顯看到，左邊為燃料元件的分布，右邊為流體區(qū)域。在兩個模型當中，燃料球的尺寸均為60mm。圖2.1中為正四棱錐堆積方式的平面尺寸圖。第一層燃料球球心之間的距離為7.002cm,第一層與第二層燃料球球心的距離為3.3912cm。通過對以上兩個簡單的模型進行數(shù)值計算，可以獲得球床內冷卻劑氦氣的速度及溫度分布規(guī)律。通

4、過對以上兩個模型的計算結果比較和分析不同堆積方式對球床內溫度分布及氦氣流動的影響。圖2.2 立方體堆積模型圖2.3 正四棱錐堆積模型圖2.4非結構網格2.2 k-模型為了使方程封閉，特別是在高雷諾數(shù)下的計算，應采用k-標準模型，該模型是在湍動能k方程基礎之上引入湍流脈動動能的耗散率方程形成的，于1972年提出。其中k為湍動能，為湍流耗散率，分別來自湍流動能方程式和耗散方程式。它們的計算公式如下： k=3/2*(u*I)2(2.4)式中u為流體流速，I為湍流強度。=vui'xkui'xk=(ui'xkui'xk)(2.5)由量綱分析可以得粘性系數(shù)i與湍動能k和耗散

5、率的關系式如下：i=Cuk2(2.6)其中有：=Cu34*k32l(2.7)式中Cu為無量綱常數(shù)，通常取0.09，k為湍動能，l為湍流尺度。2.3物理模型及求解器設置在求解控制方程時采用Fluent軟件,選用基于壓力的求解器，物理模型采用非穩(wěn)態(tài)模型，由于本文三維模型網格單元基本為多面體結構，因此選least-quares cell based這種壓力梯度來計算控制方程中的導數(shù)項。對于近壁面的處理方式采用標準壁面函數(shù)，其中壓力-速度解耦一般使用SIMPLE和SIMPLEC，對非定常流動問題或者包含比平均網傾斜度更高的網格則使用PISO算法，基于本模型網格不規(guī)則且是非定常流動模型，因此選用PISO

6、算法，其計算結果會更精確一些；但是這種算法相比較起SIMPLEC算法，可能收斂較困難，因此需要根據實際需要調整松弛因子，該因子減小有利于收斂，但是會增加收斂時間，所以需謹慎使用。由于球床內的燃料元件堆積形成的空隙通道不規(guī)則，流體壓力變化劇烈，PRESTO！主要用于高旋流，壓力急劇變化流。因此壓力基分離求解器的壓力插值算法建議使用PRESTO！。這里可以先采用一階方程進行初始迭代，再進行二階迭代計算，若MUSCL算法能夠滿足收斂條件則亦用MUSCL算法進行更精確的計算。2.4邊界條件設定及其實現(xiàn)在本模型當中，由于只取局部進行分析，設定氦氣為在高壓下的局部模型當中為不可壓縮流體。根據HTR-10的

7、運行參數(shù)知，運行壓力為3MP，由于在堆芯取的局部模型壓降相對于運行壓力來說比較小，相對壓力取0Pa時有利于分析。進出口溫度分別為250 和700 ，取其之間的任意溫度均可，為方便計算，取模型入口溫度為600K。整個堆芯入口質量流量為4.3kg/s,根據模型入口截面的大小可以算出正四棱錐模型入口流量為0.03312kg/s,立方體模型的入口流量為0.02454kg/s。對于進出口的湍動能k和湍流耗散率，可由上述公式計算得出。在流體流動過程當中，球體壁面溫度分布不均勻，不能設定球體表面溫度，第一類邊界條件不適用，這里采用第二類邊界條件。同時選用無滑移光滑表面。從HTR-10的運

8、行參數(shù)可知，其熱功率為10MW，燃料球總數(shù)量為27000個，尺寸為60mm，有以下公式可得出：q=反應堆總熱功率燃料球總表面積=32747.93W/m2所以,設定壁面的熱流量為32747.93W/m2.3 計算結果與分析本章根據第三章所建立的數(shù)學模型進行計算，得到固定床反應堆堆芯內冷卻劑氦氣的速度及其溫度分布規(guī)律,并討論了燃料球尺寸、孔隙率等對上述規(guī)律的影響。圖3.1速度、溫度、壓力及擾動分布3.1正四棱錐堆積方式從整個模型的分析可以看出，我們已得到溫度、速度甚至壓力變化的大概規(guī)律，但是要想更好更具體地分析冷卻劑氦氣的溫度、速度以及壓力分布，我們還需要從單一燃料元件著手，研究其表面壓力以及溫度

9、等的分布情況。在整個球床堆芯內，除球床周圍邊界的燃料元件外，幾乎所有燃料元件的溫度、速度和壓力的分布規(guī)律基本相同。因此，取模型當中的單一燃料元件來分析即可，在此模型當中，分別取兩種情況下其中心的單一燃料元件對比分析。如圖3.2所示，在元件間接觸點附近出現(xiàn)溫度的極大值點。在燃料元件末端也出現(xiàn)了較高的溫度，這是由于在這個地方冷卻劑流速很低，且有渦流形成，使該區(qū)域內冷卻劑出現(xiàn)滯留甚至倒流，導致對流換熱很差，從而使元件的局部溫度很高。 圖3.2 元件表面溫度分布如圖3.3所示，在Z軸-1.6cm和1.6cm的地方都出現(xiàn)了壓力的極小值點，這些地方出現(xiàn)這種現(xiàn)象的原因是由于這些位置間隙通道較小，且形成類似擴

10、壓管結構，導致壓力迅速降低而流速迅速升高。圖3.3 元件表面壓力分布如圖3.4所示，該圖為元件表面沿冷卻劑流動方向的速度分布圖，可發(fā)現(xiàn)速度最高點與壓力的極小值點相對應。原因如上所述，均是由于元件間狹小縫隙通道處形成擴壓管結構，在此處壓力迅速降低，同時速度升高。圖3.4 元件表面速度分布因此可以得出結論，在燃料元件堆積形式不變，改變尺寸或者增大流速、熱流密度等的情況下，冷卻劑在球床內的流動的規(guī)律基本一致，分布形式大致相同，只是在形式或者程度上發(fā)生了變化。3. 2立方體堆積方式當燃料元件堆積形式發(fā)生改變時，必然對對球床燃料元件溫度、速度和壓力分布和變化的影響。在不改變燃料元件尺寸的前提下，對其堆積

11、形式做一些改變，分析這一單一變量對上述分布的影響。把上述正四棱錐堆積方式變?yōu)榱⒎襟w堆積形式，其分布如下：如圖3.5所示，該圖是模型中心的單一元件溫度分布圖。在坐標-3cm和3cm處出現(xiàn)了溫度的最高值，而在0點附近也出現(xiàn)了溫度分布的較高點。其中坐標-3cm和3cm處出現(xiàn)溫度極值的原因有兩個：一是此處位于上下兩層燃料元件接觸點；二是此處流速較低，與正四棱錐模型一樣，有渦流形成，該區(qū)域內冷卻劑出現(xiàn)滯留甚至倒流，對流換熱很差。0點附近出現(xiàn)溫度較高點的原因單一，就是燃料元件之間在此處相互接觸，接觸點附近流體流動較慢，對流換熱較差。圖3.5 元件表面溫度分布如圖3.6所示，沿流動方向，壓力在元件表面的變化

12、與正四棱錐堆積形式的變化有很大區(qū)別。在立方體排布形式當中，元件上半部分壓降變化很慢，但到下半部分時，壓力會迅速降低，這與上面所述的原理一樣，在燃料元件相互形成的較小空隙處形成了類似擴壓管的通道，導致這個位置壓力突然降低。兩者有一個共同點就是這個壓力迅速降低的位置都是在與周圍元件形成的較小空隙處。但是元件堆積方式不同，這個較小間隙位置就處于元件不同位置。圖3.6 元件表面壓力分布如圖3.7所示，速度在0點附近出現(xiàn)了最大值點，在-3cm和3cm處最小，幾乎為零。與圖3.13相對應分析可知，在-3cm和3cm處為燃料元件接觸點，也是燃料元件在流動方向的最低和最高位置，在這些位置會出現(xiàn)渦流甚至倒流，速度最低。而在0點附近周圍元件形成了類似于擴壓管的結構，壓降最快，速度也增加最快。圖3.7 元件表面速度分布4.小結本文不僅計算得到了燃料元件以正四棱錐方式堆積的球床內溫度、速度和壓力的分布規(guī)律，而且得到不同尺寸和不同元件堆積方式模型下的分布規(guī)律；討論了孔隙率對其