熔鹽堆用冷凍閥的熱-結(jié)構(gòu)特性研究

李啟明 唐忠鋒 傅 遠(yuǎn) 王納秀

（中國(guó)科學(xué)院上海應(yīng)用物理研究所 上海 201800）

熔鹽堆用冷凍閥的熱-結(jié)構(gòu)特性研究

李啟明 唐忠鋒 傅 遠(yuǎn) 王納秀

（中國(guó)科學(xué)院上海應(yīng)用物理研究所 上海 201800）

反應(yīng)堆安全是任何核能系統(tǒng)開發(fā)與設(shè)計(jì)的重要組成。冷凍閥是熔鹽堆的核心技術(shù)之一，其熱力學(xué)特性直接關(guān)系到熔鹽堆核能系統(tǒng)的固有安全。本文旨在通過(guò)研究冷凍閥的熱-結(jié)構(gòu)特性并進(jìn)行優(yōu)化，探索可靠的安全設(shè)計(jì)，提高第四代核能系統(tǒng)的固有安全性，利用ANSYS軟件建立了冷凍閥的三維有限元模型，對(duì)其關(guān)閉狀態(tài)和開啟過(guò)程進(jìn)行了分析研究。結(jié)果表明：1) 冷凍閥的扁平部位由于內(nèi)部空間熔鹽少，熔鹽熔融和凍結(jié)容易控制，是實(shí)現(xiàn)開關(guān)功能和熔鹽堆過(guò)熱安全泄放的核心部位；2) 獲得了不同因素（如換熱系數(shù)、加熱功率、保溫尺寸等）對(duì)冷凍閥在關(guān)閉狀態(tài)與開啟過(guò)程中的溫度場(chǎng)及應(yīng)力大小的影響規(guī)律，為冷凍閥的優(yōu)化設(shè)計(jì)及安全運(yùn)行管理提供了依據(jù)；3) 由于初始模型中部分設(shè)計(jì)不合理（如保溫尺寸和加熱功率等），導(dǎo)致冷凍閥熱應(yīng)力過(guò)高，容易引起疲勞損傷，經(jīng)分析優(yōu)化和改進(jìn)后，最大熱應(yīng)力明顯減小，綜合性能得到很大提高。

熔鹽堆，冷凍閥，熱-結(jié)構(gòu)耦合，有限元

中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)專項(xiàng)TMSR (Thoriumbased Molten Salt Reactor，釷基熔鹽堆)屬于第四代先進(jìn)反應(yīng)堆，具有獨(dú)特的經(jīng)濟(jì)性和固有安全性等優(yōu)點(diǎn)[1,2]。冷凍閥是實(shí)現(xiàn)熔鹽堆固有安全的一個(gè)重要設(shè)備，其采用“非能動(dòng)”設(shè)計(jì)理念，即利用熔鹽自身的“冷卻-凝固/加熱-熔融”原理來(lái)實(shí)現(xiàn)冷凍閥的開啟與關(guān)閉。冷凍閥在熔鹽回路試驗(yàn)系統(tǒng)中的布置如圖1所示，閥體為水平安裝，當(dāng)反應(yīng)堆穩(wěn)定運(yùn)行時(shí)，閥內(nèi)熔鹽處于冷凝狀態(tài)，可有效隔斷管路內(nèi)的熔鹽通過(guò)該閥的流動(dòng)，實(shí)現(xiàn)其關(guān)閉功能（即關(guān)閉狀態(tài)）；借助加熱裝置（如短路加熱[3]、高頻感應(yīng)加熱等）使熔鹽熔化來(lái)實(shí)現(xiàn)冷凍閥的開啟（即開啟過(guò)程）。與常規(guī)機(jī)械閥相比，具有無(wú)運(yùn)動(dòng)部件、零泄漏和兼有過(guò)熱安全排放等優(yōu)點(diǎn)。

圖1 熔鹽回路試驗(yàn)系統(tǒng)中的冷凍閥布置Fig.1 Freeze-Valve arrangement on engineering test loop.

在冷凍閥研究歷史上，美國(guó)ORNL (Oak Ridge National Laboratory，橡樹嶺國(guó)家實(shí)驗(yàn)室)的相關(guān)工作已經(jīng)證明了冷凍閥的可行性，并且通過(guò)測(cè)試，積累了大量實(shí)驗(yàn)數(shù)據(jù)[4,5]，但沒(méi)有對(duì)冷凍閥的熱力學(xué)性能及其結(jié)構(gòu)特性做深入研究。目前國(guó)內(nèi)未見相關(guān)研究的文獻(xiàn)報(bào)道。鑒于冷凍閥自身不規(guī)則的幾何結(jié)構(gòu)、較大的溫度梯度及工作在高溫和循環(huán)熱載荷等惡劣條件下，容易產(chǎn)生疲勞破壞，為較全面地掌握冷凍閥的綜合性能并進(jìn)行優(yōu)化，提高閥體的可靠性，進(jìn)而確保熔鹽堆的安全性能，有必要對(duì)冷凍閥的熱-結(jié)構(gòu)特性進(jìn)行深入研究。

本文以FLiNaK (LiF-NaF-KF，46.5-11.5-42.0 mol%)為熔鹽材料，Inc600 (Inconel600)為閥體材料，分析了冷凍閥在關(guān)閉狀態(tài)（穩(wěn)態(tài)）和開啟過(guò)程（瞬態(tài)）中的溫度場(chǎng)及相應(yīng)的熱應(yīng)力，研究了相關(guān)因素對(duì)冷凍閥綜合性能的影響規(guī)律，討論了降低熱應(yīng)力的途徑，對(duì)熔鹽回路試驗(yàn)系統(tǒng)的冷凍閥設(shè)計(jì)提出了建議。

1 分析模型

中國(guó)科學(xué)院上海應(yīng)用物理研究所研制的冷凍閥結(jié)構(gòu)如圖2所示，其總長(zhǎng)LL為250 mm，扁平部位長(zhǎng)度FL為50 mm，內(nèi)部空間高度12.5 mm，端部圓管部分長(zhǎng)度PL為50 mm，內(nèi)徑40 mm，壁厚4mm。其中，IL為閥體端部伴熱保溫尺寸，HL為加熱尺寸。閥體材料Inc600在不同溫度下的物性參數(shù)如表1，泊松比為0.3（鎳鉻鋼材泊松比一般為0.25?0.30[6]）。熔鹽材料FLiNaK的密度、熱導(dǎo)率及比熱與溫度有關(guān)，關(guān)系式為[7]：

式中，ρ為密度，kg/m3，λ為熱導(dǎo)率，W/m·oC，Cp為比熱，J/kg·oC，T為溫度，oC。相變溫度為454oC，相變焓約為400 kJ/kg[8]。

根據(jù)相關(guān)文獻(xiàn)[9]，封閉空間內(nèi)流體的自然對(duì)流傳熱效果可以用有效熱導(dǎo)率來(lái)取代分子熱導(dǎo)率，同時(shí)狹小空間內(nèi)流體的有效熱導(dǎo)率與分子熱導(dǎo)率相差很小，考慮到冷凍閥內(nèi)空間很小且熔鹽的粘度較大，自然對(duì)流對(duì)傳熱的貢獻(xiàn)可忽略不計(jì)，因此熔鹽無(wú)論是固態(tài)還是液態(tài)，均視為實(shí)體，只考慮其熱傳導(dǎo)。冷凍閥的ANSYS 13.0 workbench[10]有限元模型如圖3所示，全部為六面體單元，其中，熱分析模型(圖3(a))中的閥體及熔鹽實(shí)體的兩端面均施加溫度載荷，閥體外表面為自然對(duì)流換熱及熱輻射散熱。將熔鹽實(shí)體抑制后的結(jié)構(gòu)分析模型如圖3(b)，其載荷為熱分析結(jié)果的溫度場(chǎng)，兩端無(wú)約束。

圖2 冷凍閥結(jié)構(gòu)示意圖Fig.2 Sketch map of a Freeze-Valve.

圖3 冷凍閥有限元模型 (a) 熱分析模型；(b) 結(jié)構(gòu)分析模型Fig.3 Finite element model of Freeze-Valve. (a) thermal analysis model; (b) structure analysis model

表1 Inconel600在不同溫度下的材料物性參數(shù)Table 1 Material parameters at different temperatures of Inconel600.

2 關(guān)閉狀態(tài)分析

2.1溫度場(chǎng)分析

在熔鹽回路試驗(yàn)系統(tǒng)中，與冷凍閥兩端相接的管道進(jìn)行了伴熱保溫（即保持管道及其內(nèi)部熔鹽溫度恒定）。有限元模型中將端面施加恒定溫度載荷700℃（當(dāng)前系統(tǒng)設(shè)計(jì)的最高運(yùn)行溫度）。冷凍閥暴露在空氣中，初步取外表面的自然對(duì)流換熱系數(shù)為7 W/m·oC、平均輻射率為0.2，外界環(huán)境溫度為22oC。初始設(shè)計(jì)的試驗(yàn)?zāi)Ｐ椭?，閥體端部圓管部分未加伴熱保溫，即IL為0 mm。熱平衡時(shí)的溫度分布云圖正向剖視如圖4(a)所示，圖4(b)為縱向剖視圖。從圖4看出，不管是正向還是縱向，穩(wěn)態(tài)下冷凍閥及其內(nèi)部熔鹽的溫度，兩端最高700 oC、中間最低304 oC，沿軸線由兩端向中間遞減，且等溫線向中心略有凹陷，這是由于兩端溫度恒定，閥體表面散熱的緣故?？傮w上，中間部位的熔鹽溫度在熔點(diǎn)以下，即冷凍閥阻塞關(guān)閉。

圖4 冷凍閥原型的溫度分布云圖 (a) 正向剖視圖；(b) 縱向剖視圖Fig.4 The typical temperature contours of Freeze-Valve. (a) forward cutaway view; (b) lateral cutaway view

為研究表面對(duì)流換熱系數(shù)、表面輻射率、伴熱保溫尺寸及保溫溫度對(duì)冷凍閥溫度分布的影響效果，圖5給出了閥中心軸線溫度分布與表面對(duì)流換熱系數(shù)、表面輻射率、伴熱保溫尺寸及保溫溫度的關(guān)系曲線。由于結(jié)構(gòu)溫度分布對(duì)稱，取其一半做圖，以閥中心點(diǎn)為0點(diǎn)。

圖5(a)顯示了對(duì)流換熱系數(shù)分別為2、7、10、15 W/m·oC時(shí)的溫度分布曲線。從圖5看出，隨著換熱系數(shù)的增大，其溫度梯度升高，整體溫度降低，中部處于熔點(diǎn)以下的熔鹽量增加。

圖5(b)示出了表面輻射率分別為0、0.05、0.1、0.3、0.5條件下的溫度變化曲線。由圖可知，表面輻射率對(duì)溫度分布的影響與對(duì)流換熱系數(shù)相似，即表面輻射率越大，溫度梯度越大，中部處于熔點(diǎn)以下的熔鹽量越多。

冷凍閥的溫度分布與其端部伴熱保溫尺寸（IL分別為20、30、40和45 mm）的關(guān)系如圖5(c)所示?？梢钥闯觯?) 兩端伴熱保溫尺寸增加，整體溫度水平明顯升高，處于熔點(diǎn)以下的熔鹽量減少；2) 當(dāng)伴熱保溫尺寸增加到45 mm時(shí)，熔鹽溫度全部處于熔點(diǎn)以上；3) 為使冷凍閥在關(guān)閉狀態(tài)下保持凍結(jié)阻塞，其內(nèi)部一定體積的熔鹽溫度應(yīng)在熔點(diǎn)以下（即凍結(jié)），因此，由圖可知當(dāng)前模型兩端的伴熱保溫尺寸不宜超過(guò)40 mm。需要說(shuō)明的是，由于設(shè)置了伴熱保溫，其溫度恒定，因此，端部部分尺寸的溫度保持一定（曲線中的橫線部分）。

圖5(d)顯示了冷凍閥兩端伴熱保溫溫度分別為500oC、700oC、800oC、1000oC時(shí)的溫度分布曲線?？梢钥闯觯S著伴熱保溫溫度升高，溫度梯度略有增加，整體溫度水平上升，處于熔點(diǎn)以下的熔鹽量降低，即內(nèi)部?jī)鼋Y(jié)熔鹽量減少。

通過(guò)上述分析得出：1) 冷凍閥在關(guān)閉狀態(tài)下的溫度分布及其內(nèi)部?jī)鼋Y(jié)的熔鹽量受表面換熱系數(shù)、表面輻射率、兩端伴熱保溫尺寸及保溫溫度影響明顯；2) 當(dāng)前模型在自然對(duì)流換熱及表面輻射率為0.2的條件下，冷凍閥兩端圓管部分的伴熱保溫尺寸不宜超過(guò)40 mm，若系統(tǒng)設(shè)計(jì)的最高運(yùn)行溫度高于700oC，伴熱保溫尺寸應(yīng)相應(yīng)減小。

圖5 中心軸的溫度分布與冷凍閥的表面對(duì)流換熱系數(shù)(a)、表面輻射率(b)、伴熱保溫尺寸(c)及兩端保溫溫度(d)的關(guān)系MPMS：熔鹽熔點(diǎn)Fig.5 Axis temperature distribution v.s. heat transfer coefficient(a), emissivity(b), insulation length(c) and the holding temperature(d). MPMS: the melting point of the molten salt

2.2應(yīng)力分析

為研究冷凍閥在關(guān)閉狀態(tài)下的溫度分布及其自身的幾何結(jié)構(gòu)對(duì)應(yīng)力分布的影響，將表面自然對(duì)流換熱系數(shù)為7 W/m·oC、表面輻射率為0.2、端部伴熱保溫溫度為700 oC、保溫尺寸為10 mm的熱分析結(jié)果（溫度場(chǎng)）作為熱載荷施加到閥體上，其典型應(yīng)力分布如圖6，其中圖6(a)為整體應(yīng)力分布云圖，圖6(b)為剖視圖。從圖6看出以下幾點(diǎn)：

1) 冷凍閥的最大熱應(yīng)力產(chǎn)生于扁平部位周圍（圖6(a)中深色區(qū)域）；

2) 外表面與內(nèi)表面的應(yīng)力分布略有差別：a)在外表面，扁平部位旁邊出現(xiàn)了應(yīng)力集中，幾何過(guò)渡區(qū)域應(yīng)力很小（圖6(a)中淺色區(qū)域）；b) 在內(nèi)表面，扁平部位中間和幾何過(guò)渡區(qū)域都出現(xiàn)了應(yīng)力集中（圖6(b)中深色區(qū)域）。

圖6 冷凍閥典型應(yīng)力分布云圖 (a) 整體應(yīng)力分布云圖；(b) 剖視圖Fig.6 Typical stress contours of Freeze-Valve. (a) global view; (b) cutaway view

為進(jìn)一步研究閥體兩端圓管部位的伴熱保溫尺寸對(duì)應(yīng)力集中的影響，改變伴熱保溫尺寸，最大應(yīng)力的變化規(guī)律如圖7所示?？梢钥闯觯?/p>

1) 增加閥體兩端圓管部位的伴熱保溫尺寸，最大應(yīng)力上升，當(dāng)保溫尺寸為40 mm時(shí)，最大應(yīng)力約50MPa；

2) 結(jié)合圖5(c)的分析結(jié)果得出：冷凍閥在關(guān)閉狀態(tài)下產(chǎn)生的應(yīng)力集中應(yīng)主要是由于溫度升高導(dǎo)致的熱膨脹受自身不規(guī)則的幾何結(jié)構(gòu)約束造成，溫度越高，熱膨脹越大，應(yīng)力越大。若要進(jìn)一步降低冷凍閥在關(guān)閉狀態(tài)下的應(yīng)力水平，應(yīng)主要從改進(jìn)結(jié)構(gòu)設(shè)計(jì)方面入手。

圖7 冷凍閥最大應(yīng)力與保溫尺寸的關(guān)系Fig.7 Maximum stress v.s. insulation length of Freeze-Valve.

3 開啟過(guò)程分析

3.1溫度場(chǎng)分析

利用ANSYS軟件分析冷凍閥的開啟過(guò)程時(shí)，為簡(jiǎn)便起見，將閥體作為內(nèi)熱源模擬短路加熱或高頻感應(yīng)加熱，熱生成率初始值設(shè)置為7.0×107W/m3，熱生成的閥體長(zhǎng)度HL為150 mm。初步取外表面自然對(duì)流換熱系數(shù)7 W/m·oC，輻射率為0.2，外界環(huán)境溫度為22 oC，伴熱保溫尺寸IL為0 mm。

冷凍閥開啟過(guò)程的初始溫度為關(guān)閉狀態(tài)時(shí)的溫度場(chǎng)，加熱35 s（設(shè)計(jì)時(shí)間）后的溫度分布如圖8(a)所示。由圖可知，只有貼近閥壁的薄層熔鹽溫度在熔點(diǎn)以上，內(nèi)部很大一部分熔鹽溫度處于熔點(diǎn)以下，即尚未熔化，因此當(dāng)前設(shè)計(jì)的加熱長(zhǎng)度為150 mm，其兩端50 mm圓管未伴熱保溫情況下，按設(shè)計(jì)要求35 s時(shí)間內(nèi)不能使熔鹽完全熔融導(dǎo)通。

為研究延長(zhǎng)加熱時(shí)間對(duì)冷凍閥開啟過(guò)程的影響效果，加熱60 s后的溫度分布如圖8(b)，可以看出：1)扁平部位內(nèi)部的熔鹽溫度已經(jīng)達(dá)到熔點(diǎn)以上，但過(guò)渡區(qū)域內(nèi)還有大塊熔鹽的溫度尚在熔點(diǎn)以下（即未熔）；2) 閥體局部溫度過(guò)高(954 oC)，已超過(guò)材料的許用溫度[11]。

通過(guò)上述分析可知，延長(zhǎng)加熱時(shí)間不能使冷凍閥有效導(dǎo)通。應(yīng)提高關(guān)閉狀態(tài)下冷凍閥中部的溫度，減少閥內(nèi)凍結(jié)的熔鹽量，結(jié)合2.1的分析結(jié)果，通過(guò)降低閥體表面對(duì)流換熱系數(shù)或輻射率，增加兩端的伴熱保溫尺寸或保溫溫度等措施來(lái)實(shí)現(xiàn)。

圖8 各種不同工況下的冷凍閥溫度分布云圖 (a) 加熱長(zhǎng)度150 mm，未增加保溫；(b) 加熱長(zhǎng)度150 mm，增加加熱時(shí)間Fig.8 Typical temperature contours of Freeze-Valve in different operating conditions. (a) heating length=150 mm, no insulation; (b) heating length=150 mm, with improved heating time

為進(jìn)一步研究加熱功率密度及閥體端部的伴熱保溫尺寸對(duì)冷凍閥開啟過(guò)程的影響效果。保溫尺寸分別為40和25 mm的溫度分布與不同加熱功率(1.0×107、3.0×107、5.0×107W/m3)的關(guān)系如圖9所示。從圖9看出：

(1) 冷凍閥兩端伴熱保溫40 mm時(shí)的整體溫度明顯比25 mm時(shí)的高。

(2) 加熱功率密度增加，中間扁平部位內(nèi)的熔鹽溫度上升，其它部位的溫度幾乎沒(méi)有變化。

(3) 當(dāng)伴熱保溫尺寸為40 mm，加熱功率密度達(dá)5.0×107W/m3后，中心軸熔鹽溫度全部處于熔點(diǎn)以上，即完全熔化，冷凍閥完全導(dǎo)通。

(4) 當(dāng)伴熱保溫尺寸為40 mm，加熱功率密度為3.0×107W/m3時(shí)，閥體扁平部位內(nèi)還有少許熔鹽塊未熔化，如圖10，其長(zhǎng)度及厚度均小于扁平部位內(nèi)部空間尺寸，將被管路內(nèi)的液態(tài)熔鹽帶出，使冷凍閥導(dǎo)通。此時(shí)，為防止緊急泄料時(shí)管道堵塞，與閥體出口端相接管道的管徑不宜小于閥體端部直徑。

(5) 當(dāng)伴熱保溫尺寸為25 mm時(shí)，中間扁平部位內(nèi)貼近閥壁的熔鹽熔化，但閥內(nèi)尚有大塊熔鹽溫度都小于熔點(diǎn)，即未熔，如圖11，不能使冷凍閥有效導(dǎo)通。因此，伴熱保溫尺寸不應(yīng)小于25 mm。

圖9 中心軸溫度分布與加熱功率的關(guān)系Fig.9 The axis temperature distribution v.s. heating power of Freeze-Valve.

圖10 伴熱保溫尺寸40 mm，加熱功率密度2.0×107W/m3時(shí)的閥內(nèi)未熔熔鹽體(a) 前視圖；(b) 側(cè)視圖Fig.10 Unmelted molten salt block in Freeze-Valve with 40 mm insulation length and 2.0×107W/m3power density. (a) front view; (b) side view

圖11 伴熱保溫尺寸25 mm，加熱功率密度2.0×107W/m3時(shí)的閥內(nèi)未熔熔鹽體(a) 前視圖；(b) 側(cè)視圖Fig.11 Unmelted molten salt block in Freeze-Valve with 25 mm insulation length and 2.0×107W/m3power density. (a) front view; (b) side view

為進(jìn)一步研究冷凍閥兩端伴熱保溫40 mm時(shí)的閥體溫度分布與加熱功率的關(guān)系，加熱功率分別為1.0×107、3.0×107和5.0×107W/m3時(shí)閥體的溫度分布如圖12所示。從圖看出：1) 當(dāng)功率密度小于7.0×107W/m3時(shí)，冷凍閥的溫度梯度較小；2) 當(dāng)功率密度超過(guò)7.0×107W/m3后，閥體的溫度梯度及局部溫度隨加熱功率密度的增加明顯升高，加熱功率密度為1.0×108W/m3時(shí)的閥體最高溫度已經(jīng)超過(guò)900 oC，超過(guò)了材料的許用溫度。因此，為防止材料失效，應(yīng)嚴(yán)格控制冷凍閥開啟時(shí)的加熱功率大小。

圖12 冷凍閥的溫度與加熱功率的關(guān)系Fig.12 Temperature distribution of Freeze-Valve v.s. heating power.

3.2應(yīng)力分析

為研究冷凍閥在開啟過(guò)程中的溫度分布及其自身的幾何結(jié)構(gòu)對(duì)應(yīng)力分布的影響，將3.1分析的溫度場(chǎng)作為熱載荷施加到閥體上，冷凍閥在不種工況下的最大應(yīng)力如表2。從表2看出：

1) 在其它條件相同的情況下，加熱功率密度越高，應(yīng)力越大，應(yīng)變也越大；

2) 加熱長(zhǎng)度為60 mm時(shí)的最大應(yīng)力明顯比加熱長(zhǎng)度為150 mm時(shí)的最大應(yīng)力高；

3) 相同工況下增加伴熱保溫尺寸，最大應(yīng)力減小。

通過(guò)上述分析可得：在允許的范圍內(nèi)，增加加熱長(zhǎng)度和伴熱保溫尺寸，減小加熱功率密度，對(duì)降低熱應(yīng)力有利。

表2 冷凍閥在各種工況下的最大應(yīng)力Table 2 The maximum stress of Freeze-Valve in different operating conditions.

另外，由表2可知，當(dāng)加熱尺寸為150 mm、伴熱保溫25 mm時(shí)，熱應(yīng)力較小，但此時(shí)冷凍閥不能完全導(dǎo)通，因此，雖然此時(shí)應(yīng)力較小，但功能上不能滿足要求。當(dāng)冷凍閥的加熱尺寸為150 mm、伴熱保溫尺寸為40 mm、加熱功率密度在2.0×107?7.0×107W/m3范圍內(nèi)時(shí)，最大應(yīng)力較小，最大應(yīng)變也未超過(guò)0.5‰。

為進(jìn)一步研究閥體在不同加熱功率密度下的應(yīng)力分布情況，圖13顯示了冷凍閥在加熱長(zhǎng)度150mm、保溫尺寸40 mm時(shí)，加熱功率密度分別為2.0×107和7.0×107W/m3時(shí)的應(yīng)力分布。從圖可看出：

1) 當(dāng)加熱功率密度為2.0×107W/m3時(shí)，應(yīng)力集中于冷凍閥扁平部位兩側(cè)(圖13(a))；

2) 當(dāng)加熱功率密度為7.0×107W/m3時(shí)，最大應(yīng)力還產(chǎn)生于冷凍閥過(guò)度區(qū)域(圖13(b))。

根據(jù)上述分析結(jié)果，結(jié)合冷凍閥溫度分布(圖12)得到：

1) 當(dāng)加熱功率密度較小時(shí)，冷凍閥溫度梯度較小，應(yīng)力集中主要應(yīng)是由于閥體的熱膨脹受其自身不規(guī)則的幾何結(jié)構(gòu)約束所致；

2) 當(dāng)加熱功率密度較大時(shí)，最大應(yīng)力還產(chǎn)生于溫度梯度最大的部位。

通過(guò)上述分析可知，在冷凍閥開啟過(guò)程中，有效控制加熱功率的大小對(duì)降低熱應(yīng)力非常重要。

圖13 冷凍閥在不同加熱功率密度下的應(yīng)力分布 (a) 2.0×107W/m3；(b) 7.0×107W/m3Fig.13 The stress distribution of Freeze-Valve with different power density. (a) 2.0×107W/m3; (b) 7.0×107W/m3

4 結(jié)論

利用ANSYS軟件，采用熱-結(jié)構(gòu)耦合的方法，研究了冷凍閥在關(guān)閉狀態(tài)及開啟過(guò)程中的性能特點(diǎn)，并分析了環(huán)境溫度、換熱系數(shù)、輻射率、保溫溫度及加熱功率、加熱長(zhǎng)度、伴熱保溫尺寸對(duì)冷凍閥溫度場(chǎng)及應(yīng)力分布的影響規(guī)律，得出如下結(jié)論：

(1) 冷凍閥的熱力學(xué)性能受表面換熱系數(shù)、輻射率、伴熱保溫尺寸及保溫溫度影響明顯，當(dāng)前冷凍閥兩端的伴熱保溫尺寸不能超過(guò)40 mm，且不宜低于25 mm，即冷凍閥中間部位暴露在空氣中的長(zhǎng)度不能小于170 mm且不宜大于200 mm。

(2) 溫度梯度越大，熱應(yīng)力越大。在允許范圍內(nèi)，增加冷凍閥兩端圓管部位的伴熱保溫尺寸或（和）加熱尺寸，對(duì)提高冷凍閥的熱力學(xué)性能和降低熱應(yīng)力有利。若增加加熱功率，有利于冷凍閥實(shí)現(xiàn)開啟功能，但其應(yīng)力升高。

(3) 在當(dāng)前試驗(yàn)系統(tǒng)最高運(yùn)行溫度700 oC前提下，為保證冷凍閥在穩(wěn)態(tài)下實(shí)現(xiàn)關(guān)閉功能，35 s內(nèi)實(shí)現(xiàn)有效開啟，并使冷凍閥的應(yīng)力水平最小，閥體的加熱尺寸為150 mm時(shí)，建議其兩端圓管部位增加40 mm伴熱保溫，加熱功率密度控制在2.0×107?7.0×107W/m3范圍內(nèi)。

基于本文的研究結(jié)果，下一步研究工作將結(jié)合實(shí)測(cè)數(shù)據(jù)檢驗(yàn)并修正邊界條件，提高分析的可靠性。并根據(jù)冷凍閥的材料特性曲線（如S-N曲線等）及熔鹽腐蝕性能[12]，分析冷凍閥在疲勞和腐蝕作用下的使用壽命，為今后熔鹽堆的冷凍閥研究積累數(shù)據(jù)和經(jīng)驗(yàn)。

需要說(shuō)明的是，本文的研究?jī)H從熱-結(jié)構(gòu)的角度分析影響冷凍閥熱應(yīng)力集中的因素，沒(méi)有考慮與冷凍閥相連的管路約束及其熱膨脹對(duì)冷凍閥應(yīng)力分布的影響。

1 JAMES A L. The fourth generation of nuclear power[J]. Progress in Nuclear Energy, 2002, 40(3?4): 301?307

2 Forsberg Charles W. Molten-Salt-Reactor Technology Gaps[C]. Proceedings of ICAPP 06, 2006: 6295

3 朱正罡. 短路加熱[J]. 醫(yī)藥工程設(shè)計(jì), 1992, (06): 20?25 ZHU Zhenggang. Short-circuit heating[J]. Pharmaceutical Engineering Design, 1992, (06): 20?25

4 Bettis E S. The aircraft reactor experiment design and construction[J]. Nuclear science and engineering: II, 1957: 804?825

5 ORNL-TM-0128[R], ORNL, 1962

6 張絲雨. 最新金屬材料牌號(hào)、性能、用途及中外牌號(hào)對(duì)照速用速查實(shí)用手冊(cè)[M]. 中國(guó)科技文化出版社, 2005: 146 ZHANG Siyu. The practical handbook of the latest metal on the grades, performance, use and antitheses of Chinese and foreign brands[M]. China Science and Technology Culture Press, 2005: 146

7 Benes O, Konings R J M. Thermodynamic properties and phase diagrams of fluoride salts for nuclear applications[J]. Journal of Fluorine Chemistry, 2009, 130 (1): 22?29

8 Khokhlov V, Korzun I, Dokutovich V, et al. Heat capacity and themal conductivity of molton ternary lithium sodium and zirconium fluorides mixtures[J]. Journal of Nuclear Materials, 2011, 410(03): 32?38

9 Bergman Theodore L, Lavine Adrienne S, Incropera Frank P, et al. Fundamental of heat and mass transfer[M]. Seventh Edition. John Wiley & SONS, Inc, New York, 2011: 621?627

10 ANSYS Inc., ANSYS 13.0 Workbench User's Guide[CP], 2012

11 ASME boiler and pressure vessel code (II), Part D: Properties (Metric) Materials[S]. 2010

12 鄭俊義, 余笑寒, 劉敏, 等. 鎳基合金受熔融氟化鹽腐蝕的同步輻射XRF分析[J]. 核技術(shù), 2011, 34(5): 336?440 ZHENG Junyi, YU Xiaohan, LIU Min, et al. SRX-ray analysis of Ni-based alloys corroded in molten fluoride salt[J]. Nuclear Techniques, 2011, 34(5): 336?440

Research on thermal and mechanical behaviour of a Freeze-Valve for molten salt reactor

LI Qiming TANG Zhongfeng FU Yuan WANG Naxiu
(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China)

Background: Reactor safety is an important component of developing and designing any nuclear energy systems. The Freeze-Valve is one of the core technologies of the molten salt reactor, which thermal-structural property is directly related to the inherent safety of the molten salt reactor nuclear system. Purpose: The purpose of this paper is to improve the inherent safety of the fourth-generation nuclear energy systems, by researching and optimizing the thermodynamic properties of Freeze-Valve, and exploring reliable safety design. Methods: A 3D Finite Element model to simulate the thermal-structural coupling behaviors was established by the ANSYS software to analyze the performances and the properties of a Freeze-Valve. Results: The analysis result shows that the flat part of the Freeze-Valve plays a very important role in engineering applications. The effects of different factors (heat transfer coefficient, heating power, and insulation size and so on) on the temperature and the stress field of the Freeze-Valve during operation are obtained, which provide some basis for optimization of design and safe operation. The thermal stress of the initial model of the Freeze-Valve is so large that it is easy to cause fatigue failure, owing to the unreasonable initial design (insulation size and heating power, for example). Conclusions: After the relative analysis and structural optimization, the maximum stress of the Freeze-Valve is significantly reduced, and the performance has been greatly improved.

Molten salt reactor, Freeze-Valve, Thermal-structural coupled, Finite element method

TL35

10.11889/j.0253-3219.2013.hjs.36.040649

中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)資助(No. XDA01020304)

李啟明，男，1984年出生，2010年于中國(guó)科學(xué)院上海應(yīng)用物理研究所獲碩士學(xué)位，核技術(shù)及應(yīng)用專業(yè)，助理研究員

2012-10-31，

2013-01-05

CLC TL35