馬棟梁 周 濤 齊 實(shí) 陳 杰 夏榜樣 肖澤軍1（華北電力大學(xué) 核科學(xué)與工程學(xué)院 北京 102206）

2（華北電力大學(xué) 核熱工安全與標(biāo)準(zhǔn)化研究所 北京 102206）

3（非能動(dòng)核能安全技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室 北京 102206）

4（中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室 成都 610041）

5（中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆熱工水力技術(shù)重點(diǎn)實(shí)驗(yàn)室 成都 610041）

超臨界水?dāng)M臨界點(diǎn)的熱膨脹系數(shù)計(jì)算研究

馬棟梁1,2,3周 濤1,2,3齊 實(shí)1,2,3陳 杰1,2,3夏榜樣4,5肖澤軍4,51（華北電力大學(xué) 核科學(xué)與工程學(xué)院 北京 102206）

2（華北電力大學(xué) 核熱工安全與標(biāo)準(zhǔn)化研究所 北京 102206）

3（非能動(dòng)核能安全技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室 北京 102206）

4（中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室 成都 610041）

5（中國(guó)核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆熱工水力技術(shù)重點(diǎn)實(shí)驗(yàn)室 成都 610041）

在臨界點(diǎn)附近時(shí)，超臨界水的熱膨脹系數(shù)很大，隨著壓力的升高，擬臨界點(diǎn)的熱膨脹系數(shù)迅速下降，擬臨界點(diǎn)的熱膨脹系數(shù)是計(jì)算臨界相變轉(zhuǎn)換數(shù)的基礎(chǔ)，對(duì)于分析超臨界壓力下的流動(dòng)不穩(wěn)定性非常重要。因此，計(jì)算超臨界水在擬臨界點(diǎn)的熱膨脹系數(shù)，對(duì)于了解和掌握超臨界水堆中能量的轉(zhuǎn)換或熱量傳遞非常重要。運(yùn)用MATLAB曲線擬合工具箱，對(duì)超臨界水的擬臨界點(diǎn)的膨脹系數(shù)進(jìn)行了擬合回歸分析。擬合得出了超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)的計(jì)算公式。該公式具有結(jié)構(gòu)簡(jiǎn)單易于計(jì)算的特點(diǎn)，最大絕對(duì)誤差為0.20 K-1，最大的相對(duì)誤差為0.19%，計(jì)算精度滿足工業(yè)研究與分析的要求。

超臨界水，擬臨界點(diǎn)，膨脹系數(shù)，擬合公式

在超臨界狀態(tài)下，液態(tài)水和汽態(tài)水是沒(méi)有明顯分界點(diǎn)的。在超臨界區(qū)，用擬臨界階段狀態(tài)點(diǎn)（擬臨界點(diǎn)）作為區(qū)分?jǐn)M液態(tài)和擬汽態(tài)的分界點(diǎn)。在擬臨界點(diǎn)時(shí)，對(duì)應(yīng)的溫度值為擬臨界點(diǎn)溫度值，對(duì)應(yīng)此時(shí)的定壓比熱容達(dá)到最大值。在超臨界水從液態(tài)向汽態(tài)逐漸轉(zhuǎn)變的過(guò)程中，其熱力學(xué)的各種物性參數(shù)會(huì)有著各種較為奇特的變化趨勢(shì)。所以在這些物性參數(shù)中，熱膨脹系數(shù)作為一個(gè)導(dǎo)出變量，是一個(gè)重要的參數(shù)。該參數(shù)在超臨界核反應(yīng)堆的流動(dòng)不穩(wěn)定性及換熱分析計(jì)算中具有重要的應(yīng)用。Ambrosini等[1]在分析超臨界壓力下的流動(dòng)不穩(wěn)定的時(shí)候，提出了臨界相變轉(zhuǎn)換數(shù)(Trans-pseudocritical Number, NTPC)的概念。在臨界相變轉(zhuǎn)換數(shù)的計(jì)算公式中，擬臨界點(diǎn)熱膨脹系數(shù)是其中需要計(jì)算的一個(gè)重要變量。擬臨界點(diǎn)的膨脹系數(shù)，即為在一定超臨界壓力下，對(duì)應(yīng)在擬臨界溫度點(diǎn)時(shí)的熱膨脹系數(shù)。值得一提的是，“超臨界熱膨脹系數(shù)”跟“擬臨界點(diǎn)的熱膨脹系數(shù)”兩者概念是不同的?！俺R界熱膨脹系數(shù)”是在一定壓力下，隨著溫度變化時(shí)由連續(xù)變化的一組數(shù)值組成。而“擬臨界點(diǎn)的熱膨脹系數(shù)”是“超臨界熱膨脹系數(shù)”中的峰值點(diǎn)[2]。計(jì)算超臨界水在擬臨界點(diǎn)的熱膨脹系數(shù)，對(duì)于理解和掌握超臨界水堆中能量的轉(zhuǎn)換或熱量的傳遞非常重要。在能源動(dòng)力工程中的其他各個(gè)領(lǐng)域分析研究中，超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)的計(jì)算亦有著重要而廣泛的應(yīng)用。但是在大量的能源與動(dòng)力類的文獻(xiàn)[2-7]中，超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)的具體計(jì)算數(shù)值卻鮮有報(bào)道。依據(jù)國(guó)際水和水蒸汽性質(zhì)協(xié)會(huì)(International Association for the Properties of Water and Steam, IAPWS)提供的1997年工業(yè)用計(jì)算標(biāo)準(zhǔn)(IAPWS-IF97)[8]，利用液態(tài)水和汽態(tài)水的熱膨脹系數(shù)差值計(jì)算的方法，可以直接計(jì)算得到在不同超臨界壓力下的液態(tài)水和汽態(tài)水的熱膨脹系數(shù)。確定擬臨界點(diǎn)溫度之后，從而可以相應(yīng)得到擬臨界點(diǎn)的熱膨脹系數(shù)。但是直接計(jì)算時(shí)，計(jì)算工作量往往過(guò)大，且計(jì)算過(guò)程較為繁瑣。超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)計(jì)算模型的提出，能夠簡(jiǎn)單迅速地計(jì)算得到在各種不同壓力下的擬臨界點(diǎn)的熱膨脹系數(shù)值，且精度范圍滿足超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)分析的要求。該公式的提出，為進(jìn)一步編程分析計(jì)算超臨界水的相變轉(zhuǎn)換數(shù)等參數(shù)具有重要意義，為便捷快速地量化分析超臨界水的流動(dòng)換熱特性奠定了基礎(chǔ)，有利于更好地對(duì)超臨界水進(jìn)行流動(dòng)不穩(wěn)定性分析。

1 超臨界水熱膨脹系數(shù)計(jì)算模型

1.1 計(jì)算范圍

IAPSW-IF97公式將液態(tài)水和汽態(tài)水的整個(gè)有效區(qū)域劃分為5個(gè)子區(qū)域，如圖1所示。

圖1 IAPWS-IF97的分區(qū)Fig.1 IAPWS-IF97 subregion.

在圖1中的分區(qū)計(jì)算公式中，第三區(qū)為超臨界區(qū)，為需要計(jì)算熱膨脹系數(shù)的區(qū)域。其中劃分第三區(qū)和第二區(qū)的邊界計(jì)算公式為：

超臨界區(qū)的計(jì)算公式適用范圍如下：623.15K≤T≤T(P)Eq.(1)，Ps(T)≤P≤100MPa。

1.2 定義公式

根據(jù)工程熱力學(xué)的定義，液態(tài)水和汽態(tài)水的熱膨脹系數(shù)是在等壓情況下比體積隨著溫度的變化率。計(jì)算公式如下：

式中：pβ為熱膨脹系數(shù)，K-1；υ為比體積，m3·kg-1；T為熱力學(xué)溫度，K。

在超臨界區(qū)域范圍內(nèi)，第三區(qū)給出的基本方程形式是Helmholtz比自由能方程(,T)fρ（對(duì)應(yīng)的無(wú)量綱形式的Helmholtz比自由能方程為(,)φ δ τ）。在超臨界區(qū)域中，(,)φ δ τ以折算密度δ和折算溫度τ為獨(dú)立的變量。壓力的關(guān)系表達(dá)式如下：

式中：ρ為密度，kg·m-3；φδ為φ( δ ,τ)對(duì)折算密度的偏導(dǎo)數(shù)；基準(zhǔn)密度ρc=322 kg·m-3。

1.3 計(jì)算公式

因?yàn)樵诔R界區(qū)域中給出的基本方程不是直接以P和T為獨(dú)立變量的，所以，要使用隱含導(dǎo)數(shù)求導(dǎo)的方法，將P看成是相對(duì)于T獨(dú)立的變量，然后將壓力的關(guān)系方程兩邊同時(shí)對(duì)T求偏導(dǎo)數(shù)，即：

解上述方程可得：

將式(6)帶入熱膨脹系數(shù)的定義式(3)中，得到熱膨脹系數(shù)的計(jì)算公式如下：

式中：δτφ為(,)φ δ τ對(duì)折算密度和折算溫度的二階偏導(dǎo)數(shù)；δδφ為(,)φ δ τ對(duì)折算密度的二階偏導(dǎo)數(shù)。

2 超臨界水熱膨脹系數(shù)參數(shù)分析

2.1 不同壓力的熱膨脹系數(shù)

當(dāng)水在超臨界狀態(tài)下，液態(tài)水和汽態(tài)水是沒(méi)有明顯分界點(diǎn)的。當(dāng)水處于溫度較低時(shí)，可以認(rèn)為此時(shí)水的狀態(tài)是和亞臨界以下的液態(tài)水是一致的。當(dāng)水處于溫度較高時(shí)，可以認(rèn)為此時(shí)水的狀態(tài)是和亞臨界以下的汽態(tài)水是一致的。但是在從液態(tài)水向汽態(tài)水過(guò)渡的過(guò)程中，此時(shí)處于一種液態(tài)水和汽態(tài)水的混合狀態(tài)。在一定的壓力條件下，在某一臨界溫度點(diǎn)附近時(shí)，此時(shí)液態(tài)水和汽態(tài)水的熱膨脹系數(shù)達(dá)到最大值，之后迅速減小。超臨界水在不同壓力條件下的熱膨脹系數(shù)變化趨勢(shì)如圖2所示。

圖2 不同壓力下超臨界水的熱膨脹系數(shù)變化趨勢(shì)Fig.2 Thermal expansion coefficient of supercritical water under different pressure.

由圖2可見(jiàn)，在壓力一定的條件下，超臨界水在從液態(tài)水向汽態(tài)水狀態(tài)轉(zhuǎn)變的過(guò)程中，熱膨脹系數(shù)首先迅速升高，達(dá)到峰值之后又迅速下降。超臨界水的熱膨脹系數(shù)所對(duì)應(yīng)的溫度在擬臨界溫度點(diǎn)附近，此時(shí)的熱膨脹系數(shù)為擬臨界點(diǎn)的熱膨脹系數(shù)。隨著壓力的升高，擬臨界點(diǎn)的熱膨脹系數(shù)的峰值逐漸降低，同時(shí)對(duì)應(yīng)的擬臨界溫度值卻在逐漸升高。

通過(guò)對(duì)式(3)進(jìn)行差值計(jì)算，選擇其中熱膨脹系數(shù)的峰值，從而得出不同壓力下的擬臨界點(diǎn)的溫度值和擬臨界點(diǎn)的熱膨脹系數(shù)的數(shù)值，如表1所示。

通過(guò)表1計(jì)算得出超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)，可以直觀地看出在各種壓力下的擬臨界點(diǎn)的熱膨脹系數(shù)數(shù)值的大小，可以供相關(guān)研究人員對(duì)超臨界壓力下的水的擬臨界點(diǎn)的熱膨脹系數(shù)進(jìn)行查詢參考使用。

表1 不同壓力對(duì)應(yīng)的擬臨界點(diǎn)的熱膨脹系數(shù)Table1 Pseudo-critical thermal expansion coefficient under different pressure.

2.2 擬臨界點(diǎn)的熱膨脹系數(shù)變化趨勢(shì)

在不同壓力下的超臨界水?dāng)M臨界溫度值和擬臨界點(diǎn)的熱膨脹系數(shù)的變化趨勢(shì)情況如圖3所示。

圖3 擬臨界溫度和擬臨界點(diǎn)的熱膨脹系數(shù)變化趨勢(shì)Fig.3 Trend of pseudo-critical temperature and pseudo-critical thermal expansion coefficient.

通過(guò)圖3可以看出，隨著壓力升高，超臨界水的擬臨界溫度值逐漸升高，而超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)是迅速下降的，下降速度呈現(xiàn)為指數(shù)級(jí)別的速度下降。當(dāng)下降到一定程度之后，擬臨界點(diǎn)的熱膨脹系數(shù)變化趨勢(shì)趨于平緩，說(shuō)明達(dá)到一定壓力之后，隨著溫度的變化，超臨界水的熱膨脹變化是比較小的。

3 擬合公式的確定及驗(yàn)證

3.1 公式的確定

由于在表1中計(jì)算所得的各壓力下的擬臨界點(diǎn)的熱膨脹系數(shù)是通過(guò)逐點(diǎn)溫度差值計(jì)算熱膨脹系數(shù)而得到的，計(jì)算工作量較大，且工作較為繁瑣。所以根據(jù)表1中計(jì)算得到的擬臨界點(diǎn)的熱膨脹系數(shù)的原始數(shù)據(jù)，利用MATLAB擬合數(shù)據(jù)工具箱對(duì)超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)值進(jìn)行了擬合分析。通過(guò)擬合分析得到了如下兩種公式結(jié)構(gòu)的擬合公式，計(jì)算結(jié)果的回歸情況較為理想。擬合公式分別如下：

式中：c1=0.327(0.317,0.337)；c2=-7.112(-7.342, -6.882)；d1=-44(-44.01,-43.98)；d2=484(483.7, 484.3)。其中括弧內(nèi)的系數(shù)參數(shù)范圍為回歸系數(shù)為95%置信區(qū)間的界限參數(shù)。

3.2 公式的驗(yàn)證

擬合公式一和擬合公式二計(jì)算擬臨界膨脹系數(shù)的計(jì)算結(jié)果及相對(duì)誤差和絕對(duì)誤差如表2、3所示。

為了將表2、3中的擬合公式計(jì)算結(jié)果進(jìn)行直觀對(duì)比，表2、3中直接計(jì)算結(jié)果和擬合公式一和擬合公式二得出的結(jié)果對(duì)比情況如圖4所示。

通過(guò)圖4中擬合公式中的計(jì)算結(jié)果和直接計(jì)算得出的擬臨界點(diǎn)的熱膨脹系數(shù)的對(duì)比可以看出，擬合公式一在壓力較低時(shí)與實(shí)際計(jì)算得出的數(shù)值誤差較小，但是壓力超過(guò)23MPa之后擬合公式一得出的計(jì)算結(jié)果普遍比實(shí)際值偏高。擬合公式二的誤差平方和(Sum of the Squared Errors, SSE)跟公式一的誤差平方和相比較小，均方根誤差(Root Mean Square Error, RMSE)也比公式一的均方根誤差要小很多。擬合公式二的回歸系數(shù)為1，表明擬合結(jié)果的實(shí)際趨勢(shì)和實(shí)際計(jì)算結(jié)果的趨勢(shì)是完全線性相關(guān)的。但是在壓力低于23MPa的壓力范圍內(nèi)，通過(guò)擬合公式二計(jì)算得出的結(jié)果跟擬合公式一相比偏高較多。

表2 擬合公式一和擬臨界點(diǎn)的熱膨脹系數(shù)的結(jié)果對(duì)比Table2 Contrast of pseudo-critical thermal expansion coefficient with the result of fitting Formula 1.

表3 擬合公式二和擬臨界點(diǎn)的熱膨脹系數(shù)的結(jié)果對(duì)比Table3 Contrast of pseudo-critical thermal expansion coefficient with the result of fitting Formula 2.

圖4 直接計(jì)算結(jié)果和擬合公式一(a)、擬合公式二(b)計(jì)算結(jié)果的對(duì)比Fig.4 Contrast of directly calculation result with the result of fitting Formula 1 (a) and fitting Formula 2 (b).

3.3 誤差分析

擬合公式一和擬合公式二的計(jì)算結(jié)果的絕對(duì)誤差和相對(duì)誤差分布情況如圖5所示。通過(guò)圖5(b)在公式一和公式二的相對(duì)誤差分布情況也可以看出，在壓力較低時(shí)公式一的相對(duì)誤差較小，在壓力較高時(shí)公式二的相對(duì)誤差較小。所以綜合考慮后確定，當(dāng)超臨界壓力小于等于22.5MPa時(shí)，建議選用公式一；當(dāng)超臨界壓力大于等于23.0MPa時(shí)，建議選用公式二。當(dāng)超臨界壓力介于22.5MPa和23.0MPa之間時(shí)，則采用擬合公式一和擬合公式二組合的方式計(jì)算。具體組合過(guò)渡形式見(jiàn)§3.4的最終公式。組合計(jì)算形式的采用，保證了擬臨界點(diǎn)的熱膨脹系數(shù)的計(jì)算在不同壓力區(qū)間變化時(shí)的計(jì)算結(jié)果具有連續(xù)光滑性。

圖5 擬合公式一與擬合公式二的絕對(duì)誤差(a)和相對(duì)誤差(b)對(duì)比Fig.5 Absolute error contrast (a) and relative error contrast (b) between fitting Formula 1 and fitting Formula 2.

3.4 公式最終形式

經(jīng)過(guò)§3.3的誤差分析，最終確定的公式形式如式(13)所示：

用IAPWS-IF97標(biāo)準(zhǔn)計(jì)算擬臨界點(diǎn)的熱膨脹系數(shù)時(shí)，需要首先計(jì)算已知壓力下對(duì)應(yīng)的擬臨界點(diǎn)溫度，進(jìn)而通過(guò)壓力和溫度，計(jì)算出該壓力和溫度下的熱膨脹系數(shù)，即擬臨界點(diǎn)的熱膨脹系數(shù)。而采用新擬合的式(13)計(jì)算擬臨界點(diǎn)的熱膨脹系數(shù)時(shí)，由于擬臨界點(diǎn)的熱膨脹系數(shù)是壓力的函數(shù)，所以僅知道壓力，即可通過(guò)擬合公式求出該超臨界壓力下對(duì)應(yīng)擬臨界點(diǎn)的熱膨脹系數(shù)。對(duì)比該兩種不同的計(jì)算方法，運(yùn)用本文公式可以加快反應(yīng)堆程序的計(jì)算速率，進(jìn)而減少計(jì)算時(shí)間。

4 結(jié)語(yǔ)

根據(jù)液態(tài)水和汽態(tài)水的熱膨脹系數(shù)計(jì)算定義，直接計(jì)算得出了超臨界水在各個(gè)壓力下的擬臨界點(diǎn)的熱膨脹系數(shù)。然后根據(jù)直接計(jì)算得出的各個(gè)壓力下的擬臨界點(diǎn)的熱膨脹系數(shù)，運(yùn)用MATLAB曲線擬合工具箱，對(duì)超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)進(jìn)行了擬合回歸分析，并得出了兩組計(jì)算超臨界水的擬臨界點(diǎn)的熱膨脹系數(shù)的計(jì)算公式。

1) 超臨界水?dāng)M臨界點(diǎn)的熱膨脹系數(shù)，在臨界點(diǎn)附近時(shí)達(dá)到最大值，隨著壓力的上升，擬臨界點(diǎn)的熱膨脹系數(shù)迅速下降，當(dāng)壓力超過(guò)25.0MPa之后，擬臨界點(diǎn)的熱膨脹系數(shù)的變化趨于平緩。

2) 當(dāng)超臨界壓力低于22.5MPa時(shí)，建議選用擬合公式一計(jì)算擬臨界膨脹系數(shù)；當(dāng)超臨界壓力大于等于23.0MPa時(shí)，建議選用擬合公式二計(jì)算擬臨界膨脹系數(shù)。當(dāng)壓力介于22.5MPa和23.0MPa之間時(shí)，采用上述兩公式組合的方式進(jìn)行計(jì)算。最終推薦計(jì)算公式形式如式(13)所示。最終計(jì)算公式最大的絕對(duì)誤差為0.201K-1，最大的相對(duì)誤差為0.19%，滿足計(jì)算分析擬臨界點(diǎn)的熱膨脹系數(shù)的要求，且具有結(jié)構(gòu)簡(jiǎn)單易于編程計(jì)算的特點(diǎn)。

1 Ambrosini W, Sharabi M. Dimensionless parameters in stability analysis of heated channels with fluids at supercritical pressures[J]. Nuclear Engineering and Design, 2008, 238: 1917-1929. DOI: 10.1016/j. nucengdes.2007.09.008

2 Pioro I L, Khartabil H F, Duffey R B. Heat transfer to supercritical fluids flowing in channels-empirical correlations (survey)[J]. Nuclear Engineering and Design, 2004, 230: 69-91. DOI: 10.1016/j.nucengdes.2003.10.010

3 李佳, 史秀敏. 基于IAPWS-IF97水和蒸汽體膨脹系數(shù)和壓縮系數(shù)[J]. 燃?xì)馀c熱力, 2012, 4(4): A01-A04. DOI: 10.3969/j.issn.1000-4416.2012.04.001 LI Jia, SHI Xiumin. Volumetric expansion coefficient and compressibility coefficient of water and steam based on IAPWS-IF97[J]. Gas and Heat, 2012, 4(4): A01-A04. DOI: 10.3969/j.issn.1000-4416.2012.04.001

4 傅晟威, 周翀, 楊燕華, 等. 跨臨界瞬態(tài)數(shù)值模擬中的高精度擬臨界溫度計(jì)算[J]. 上海交通大學(xué)學(xué)報(bào), 2011, 3(3): 413-417 FU Shengwei, ZHOU Chong, YANG Yanhua, et al. Calculation of high-precision pseudo-critical temperatures for trans-critical transient analysis[J]. Journal of Shanghai Jiaotong University, 2011, 3(3): 413-417

5 李精精，周濤，段軍, 等. 基于遺傳神經(jīng)網(wǎng)絡(luò)方法的流動(dòng)不穩(wěn)定起始點(diǎn)研究[J]. 核動(dòng)力工程, 2014, 4(2): 63-66 LI Jingjing, ZHOU Tao, DUAN Jun, et al. Study on onset of flow instability by genetic neural network[J]. Nuclear Power Engineering, 2014, 4(2): 63-66

6 陳瑋瑋, 方賢德, 商輝, 等. 超臨界壓力下豎直管內(nèi)水的傳熱關(guān)系式研究[J]. 工程熱物理學(xué)報(bào), 2016, 1(1): 104-110 CHEN Weiwei, FANG Xiande, SHANG Hui, et al. Investigation of heat transfer correlations of water under supercritical pressure in vertical tubes[J]. Journal of Engineering Thermo Physics, 2016, 1(1): 104-110

7 Liao Y X, Lucas D, Krepper E, et al. Flashing evaporation under different pressure levels[J]. Nuclear Engineering and Design, 2013, 265: 801-813. DOI: 10.1016/j.nucengdes.2013.09.027

8 Wagner W, Kruse A. 水和蒸汽的性質(zhì)[M]. 項(xiàng)紅衛(wèi), 譯.北京: 科學(xué)出版社, 2003 Wagner W, Kruse A. Properties of water and steam[M]. XIANG Hongwei, Tr. Beijing: Science Press, 2003

Calculation research on pseudo-critical thermal expansion coefficient of supercritical water

MA Dongliang1,2,3ZHOU Tao1,2,3QI Shi1,2,3CHEN Jie1,2,3XIA Bangyang4,5XIAO Zejun4,51(School of Nuclear Science and Engineering, North China Electric Power University, Beijing 102206, China)
2(Institute of Nuclear Thermal-hydraulic Safety and Standardization, North China Electric Power University, Beijing 102206, China)
3(Beijing Key Laboratory of Passive Safety Technology for Nuclear Energy, Beijing 102206, China)
4(Key Laboratory of Reactor System Design Technology, Nuclear Power Institute of China, Chengdu 610041, China)
5(Science and Technology on Reactor System Design Technology Laboratory, Nuclear Power Institute of China, Chengdu 610041, China)

Background: The supercritical water thermal expansion coefficient is very large when the reactor is approaching the critical point, but drops rapidly with the increase of pressure. The pseudo-critical thermal expansion coefficient is the basis for calculation of the critical phase transition transformation number which is of great importance in the analysis of flow instability of supercritical pressure water, therefore its calculation is very important to understand and master the conversion of energy or heat transfer in supercritical water reactor. Purpose: This study aims to derive a fitting formula for direct calculation of the pseudo-critical expansion coefficient. It is useful for thermal-hydraulic program analysis. Methods: First of all, the pseudo-critical temperature need to be calculated, then a regression analysis for the pseudo-critical thermal expansion coefficient of supercritical water is fitted by using

Supercritical water, Pseudo-critical, Thermal expansion coefficient, Fitting formula

MA Dongliang, male, born in 1982, graduated from North China Electric Power University with a master’s degree in 2008, doctoral student, focusing on nuclear thermal hydraulic and safety

ZHOU Tao, E-mail: zhoutao@ncepu.edu.cn

TL99

10.11889/j.0253-3219.2016.hjs.39.120601

中核核反應(yīng)堆熱工水力技術(shù)重點(diǎn)實(shí)驗(yàn)室基金(No.20130901)、核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室基金(No.2015BJ0151)資助

馬棟梁，男，1982年出生，2008年于華北電力大學(xué)獲碩士學(xué)位，博士研究生，研究方向?yàn)楹藷峁にεc安全

周濤，E-mail: zhoutao@ncepu.edu.cn

Supported by Nuclear Reactor Thermal Hydraulic Technology Key Laboratory Fund (No.20130901), Key Laboratory of Nuclear Reactor System Design Technology (No.2015BJ0151)

2016-10-12，

2016-11-10

MATLAB curve fitting toolbox. The fitting calculation formula about the pseudo-critical thermal expansion coefficient is derived from those fitting data in the form of combination formula. Results: The calculation result is continuous smoothness with maximum absolute error no more than 0.20 K-1, and the biggest relative error is 0.19%. Conclusion: The formula is simple in structure and easy to calculate with high efficiency, and the calculation precision meets the requirements of industry research and analysis.