核燃料棒溫度場(chǎng)和系統(tǒng)熵增分布計(jì)算與研究

張功偉，張鈞波，張 敏

(1.南京理工大學(xué)能源與動(dòng)力工程學(xué)院, 江蘇南京210094;2.南京師范大學(xué)泰州學(xué)院，江蘇泰州225300)

?

核燃料棒溫度場(chǎng)和系統(tǒng)熵增分布計(jì)算與研究

張功偉1，張鈞波2，張 敏1

(1.南京理工大學(xué)能源與動(dòng)力工程學(xué)院, 江蘇南京210094;2.南京師范大學(xué)泰州學(xué)院，江蘇泰州225300)

基于熱力學(xué)第二定律以及壓水核反應(yīng)堆燃料棒穩(wěn)態(tài)傳熱偏微分方程的一般形式，通過(guò)熵增的數(shù)值計(jì)算，討論分析燃料棒內(nèi)熱量的傳遞過(guò)程和方向，以及能量品質(zhì)的得失。先對(duì)典型二維矩形域傳熱問(wèn)題進(jìn)行數(shù)值計(jì)算，并采用解析解對(duì)該數(shù)值方法進(jìn)行了驗(yàn)證。然后對(duì)含有內(nèi)熱源的單根燃料棒進(jìn)行傳熱計(jì)算，討論其溫度和熵增分布情況，通過(guò)熵增云圖分析可以展現(xiàn)燃料棒內(nèi)熱量的傳遞過(guò)程和品質(zhì)變化，為核反應(yīng)堆熱工設(shè)計(jì)提供有益參考。

壓水核反應(yīng)堆；燃料棒；熵增

保證核反應(yīng)堆安全運(yùn)行一直是近代核科學(xué)研究的重要課題，而研究燃料棒的熱量傳遞過(guò)程是其中一項(xiàng)重要內(nèi)容。對(duì)于壓水核反應(yīng)堆燃料棒傳熱系統(tǒng)的計(jì)算分析，一般都是從熱力學(xué)第一定律的角度出發(fā)[2]，通過(guò)數(shù)值傳熱學(xué)求解燃料棒內(nèi)溫度場(chǎng)分布，或是應(yīng)用流體力學(xué)計(jì)算軟件對(duì)燃料組件內(nèi)通道進(jìn)行三維數(shù)值模擬[3,4]，從而獲得燃料組件與外圍冷卻劑對(duì)流換熱的詳細(xì)狀況。然而對(duì)于燃料棒內(nèi)熱量傳遞的方向性，以及能量品質(zhì)得失的研究也十分必要的。

本文將從熱力學(xué)第二定律出發(fā)，通過(guò)計(jì)算系統(tǒng)熵增，重點(diǎn)討論燃料棒內(nèi)熱量傳遞過(guò)程的方向性。首先采用一個(gè)簡(jiǎn)單的二維矩形域傳熱問(wèn)題，對(duì)比分析求解出的數(shù)值解和解析解，并通過(guò)函數(shù)極小值求解，確定熱量的極值點(diǎn)，從而驗(yàn)證該數(shù)值方法的正確性。最后以單根燃料棒為模型，計(jì)算在核反應(yīng)堆正常工況下燃料棒的熵增分布情況，進(jìn)而可以從熱量傳遞和能量品質(zhì)變化的角度進(jìn)行更加合理的熱工設(shè)計(jì)，即更加快速有效的導(dǎo)出燃料棒內(nèi)反應(yīng)熱，保證反應(yīng)堆的安全運(yùn)行。

1 基本計(jì)算方程

在廣義正交曲線(xiàn)坐標(biāo)系內(nèi)一般形式的非穩(wěn)態(tài)導(dǎo)熱微分方程為：

(1)

式中：qv為單位體積的內(nèi)熱源，k為材料的熱導(dǎo)率。系數(shù)a=a1a2a3，根據(jù)所選坐標(biāo)系的不同，尺度系數(shù)a1、a2、a3既可以是常數(shù)，也可以是坐標(biāo)的函數(shù)。

假設(shè)系統(tǒng)內(nèi)部的過(guò)程是可逆的，則熵增函數(shù)可以表示為：

(2)

式中：δq為熱傳遞過(guò)程中微元熱量。

2 二維矩形域算例

2.1 方程的求解

在直角坐標(biāo)系內(nèi)，選取一長(zhǎng)為a，寬為b的矩形為計(jì)算域，假設(shè)該矩形常物性、無(wú)內(nèi)熱源，根據(jù)公式(1)，其穩(wěn)態(tài)導(dǎo)熱微分方程可簡(jiǎn)寫(xiě)為：

(3)

其邊界條件定為：

T(x,y)=x+y+xy

(4)

數(shù)值計(jì)算求得系統(tǒng)達(dá)到穩(wěn)態(tài)時(shí)的溫度Ta為：

(5)

(6)

系統(tǒng)總熵增為：

(7)

積分式(7)有如下兩種形式：

(8)

對(duì)式(8)分別進(jìn)行一次、二次求導(dǎo)，可求得系統(tǒng)總熵增隨x和y的變化率：

(9)

2.2 系統(tǒng)熵增的計(jì)算與分析

取系統(tǒng)ρcp=1，初始溫度T0=2℃，并取矩形幾何尺寸為a=2 m，b=1 m。通過(guò)公式(5)可以得到該矩形域每個(gè)點(diǎn)的溫度值，獲得溫度場(chǎng)，同時(shí)通過(guò)式(6)和式(7)，運(yùn)用數(shù)值方法積分得到矩形域熵增的分布情況。

如圖1矩形域溫度分布云圖所示，隨著坐標(biāo)x和y的不同，每一點(diǎn)的溫度值不斷變化，溫度最低點(diǎn)出現(xiàn)在原點(diǎn)，然后隨著x和y的增大，最終在頂點(diǎn)(2 , 1)處達(dá)到最大值。

圖1 矩形域的溫度分布圖Fig.1 Temperature distribution of rectangular region

根據(jù)式(6)，系統(tǒng)局部熵增變化情況可從沿x截面和y截面兩個(gè)方向分析，如圖2所示。從圖中可以看到系統(tǒng)局部熵增從x=0截面到x=2截面的變化過(guò)程，從開(kāi)始的TaT0，系統(tǒng)處于吸熱狀態(tài)，局部熵增沿y方向升高。同理，可以分析右圖局部熵增沿y截面的變化情況。

圖3為該矩形域總熵增的變化云圖。對(duì)比圖2，這里可以更清晰看出，在x=0及y=0截面處，TaT0，系統(tǒng)吸熱，總熵增隨著坐標(biāo)逐漸升高。同時(shí)可以發(fā)現(xiàn)隨著系統(tǒng)總吸熱量不斷增加，逐漸大于起初的放熱量，最終系統(tǒng)總熵增是正值。

根據(jù)二元函數(shù)極值準(zhǔn)則，可以求得系統(tǒng)熵增的極值點(diǎn)大約位于點(diǎn)(1.0 , 0.5)附近，此時(shí)總熵增值為-0.27。

3 燃料棒內(nèi)的熵增計(jì)算

前面通過(guò)一個(gè)簡(jiǎn)單的傳熱問(wèn)題，驗(yàn)證了該方法的合理性和準(zhǔn)確性，這里將對(duì)單根燃料棒內(nèi)的熵增問(wèn)題進(jìn)行分析。由于軸對(duì)稱(chēng)性，一根為常物性圓柱形燃料棒的溫度場(chǎng)與角度無(wú)關(guān)，則根據(jù)公式(1)，其穩(wěn)態(tài)導(dǎo)熱微分方程可簡(jiǎn)寫(xiě)為公式(10)：

(10)

在核反應(yīng)堆正常運(yùn)行工況下，燃料棒內(nèi)的熱源項(xiàng)通?？杀硎緸椋?/p>

(11)

建模時(shí)的坐標(biāo)原點(diǎn)取燃料棒中心，其上下壁面為定溫邊界條件，在側(cè)面r=R燃料棒包殼處與外圍冷卻劑進(jìn)行對(duì)流換熱，即模型的邊界條件可數(shù)學(xué)描述為：

圖2 系統(tǒng)分別沿x,y截面局部熵增Fig.2 The local entropy of system along the section x,y

圖3 系統(tǒng)沿xy平面的總熵增Fig.3 Total entropy of system along plane xy

(12)

式中：Tw(z)為r=R處的溫度，即Tw(z)=T(R,z)；Tf(z)為流體溫度，其線(xiàn)性化形式為：

(13)

從而可以求得此時(shí)燃料棒最終的穩(wěn)態(tài)溫度為，

(14)

同樣，這里假設(shè)燃料棒的初始溫度為T(mén)0，則在其達(dá)到穩(wěn)態(tài)時(shí)，系統(tǒng)的局部熵增可分成沿rφ和rz截面兩個(gè)方向，分別表示為式(15)和式(16)，總熵增為式(17)。

πρcpr2lnT0

(15)

(16)

πρcpr2zlnT0

(17)

這里選取單根燃料棒具有代表性的一段作為數(shù)值模型，其高度為L(zhǎng)=30 mm，半徑R=5 mm，初始溫度T0=(T1+T2)/2=568 K，其他各項(xiàng)熱工參數(shù)如表1所示。

表1 燃料棒熱工參數(shù)

由于對(duì)稱(chēng)性，取二分之一圓柱體燃料棒作為計(jì)算模型，圖4為其在rφ和rz截面的溫度場(chǎng)云圖，很容易看出燃料棒在這兩個(gè)平面的溫度場(chǎng)均以原點(diǎn)為中心向外環(huán)形擴(kuò)散，越靠近外圍溫度值越低，在中心附近溫度值最高。因?yàn)槠錅囟确植家?guī)律與角度無(wú)關(guān)，所以這里在rφ截面的熵增分析也只考慮隨半徑的變化情況。

圖4 燃料棒rφ,rz平面的溫度分布示意圖Fig.4 Temperature distribution of plane rφ,rz

圖5 底面沿半徑方向的溫度及局部熵增曲線(xiàn)圖Fig.5 Temperature distribution and local entropy of underside along radius direction

圖6 不同高度沿徑向的溫度及局部熵增曲線(xiàn)圖Fig.6 Temperature distribution and local entropy of different height along radius direction

從圖7的云圖可以更清晰看出上述的變化趨勢(shì)，由下往上，在底面有Ta

圖7 燃料棒xz平面的局部熵增示意圖Fig.7 Local entropy of plane xz

4 結(jié)論

對(duì)于壓水核反應(yīng)堆燃料棒傳熱系統(tǒng)的計(jì)算分析，溫度場(chǎng)的求解固然十分重要，但熱傳遞過(guò)程中熵增的計(jì)算也具有極大意義和實(shí)用價(jià)值。本文基于熱力學(xué)第二定律，從溫度分布和熵增兩個(gè)角度分析壓水堆燃料棒正常運(yùn)行工況下的傳熱過(guò)程，尤其在熱量傳遞的方向性上取得可信結(jié)果，從而可以對(duì)傳熱過(guò)程的優(yōu)劣程度進(jìn)行可靠性評(píng)估，也可以從能量品質(zhì)變化的角度優(yōu)化反應(yīng)堆熱工設(shè)計(jì)，對(duì)于進(jìn)一步提高反應(yīng)堆運(yùn)行的安全性有一定參考價(jià)值。

[1] 陳文振,于雷,郝建立.核動(dòng)力裝置熱工水力[M].北京:中國(guó)原子能出版社,2013.

[2] 劉會(huì)娟,張敏,彭文杰.圓柱形帶反射層反應(yīng)堆的數(shù)值傳熱計(jì)算[J].核科學(xué)與工程,2010,30(1):42-47.

[3] 羅磊,陳文振,陳志云,等.單個(gè)燃料元件熱工水力三維數(shù)值模擬[J].海軍工程大學(xué)學(xué)報(bào), 2011, 23(1): 63-66[4] 宋磊, 郭赟, 曾和義. 板狀燃料組件入口堵流事故下流場(chǎng)和溫度場(chǎng)的瞬態(tài)數(shù)值計(jì)算[J]. 核動(dòng)力工程, 2014, 35(3):6-10

[5] Zhang, M. Modeling of Radiative Heat Transfer and Diffusion Processes Using Unstructured Grid [D]. Cookeville: Tennessee Technological University, 2000.

[6] Patankar, S. V. Numerical Heat Transfer and Fluid Flow [M]. New York: Hemisphere Publishing, 1981.

Study on the Distribution of Temperature and Entropy for a Reactor Fuel Rod

ZHANG Gong-wei1, ZHANG Jun-bo2, ZHANG Min1

(1.School of Power & Energy Engineering, Nanjing University of Science & Technology, Nanjing, Jiangsu Prov. 210094；2. Taizhou College, Nanjing Normal University, Taizhou, Jiangsu Prov.225300)

According to the second law of thermodynamics and steady partial differential equations for heat transfer in a fuel rod of pressurized water reactor, entropy production of the rod is numerically calculated to discuss the process and direction of heat transfer, as well as the gain and loss of energy in the fuel rod. A typical heat transfer case for 2D rectangle region is conducted firstly and the analysis solution is used to validate the accuracy and robustness of the code. Finally, a single fuel rod with an internal heat source is numerically studied. The distribution for temperature and entropy in the rod is discussed. The contour map of entropy distribution indicates that the numerical method is able to reveal the transfer and variation of heat transfer in the rod clearly. It is somewhat useful for thermal performance design of nuclear reactor.

Pressurized water reactor; Fuel rods; Entropy production

2016-12-29

張功偉(1991—)，男，江蘇人，碩士研究生，從事計(jì)算流體力學(xué)及計(jì)算傳熱學(xué)研究方向

TL331

A

0258-0918(2017)02-0333-07