李懿　武建新

摘 要：在非能動(dòng)安全殼冷卻系統(tǒng)CFD冷凝模型和蒸發(fā)模型的開發(fā)過程中，應(yīng)用了計(jì)算流體動(dòng)力學(xué)程序平臺(tái)CASTEM。采用了Buleev湍流模型和k-ε湍流模型作為冷凝模型，采用了G.Aiello發(fā)展液膜蒸發(fā)模型作為蒸發(fā)模型。經(jīng)過研究和試驗(yàn)證明，試驗(yàn)值和計(jì)算值之間，具有較為良好的吻合度，因此，在非能動(dòng)安全殼冷卻系統(tǒng)當(dāng)中，能夠采用建立的冷凝模型和蒸發(fā)模型對其穩(wěn)態(tài)下的傳熱傳質(zhì)特性進(jìn)行模擬。

關(guān)鍵詞：非能動(dòng)安全殼冷卻系統(tǒng)；CFD；冷凝模型；蒸發(fā)模型

前言

在先進(jìn)壓水堆的一些假想事故，例如失水事故當(dāng)中，利用非能動(dòng)安全殼冷卻系統(tǒng)對安全殼進(jìn)行排熱處理。設(shè)計(jì)了雙層的安全殼結(jié)構(gòu)，其中事故之后的傳熱面為內(nèi)層的鋼安全殼。在發(fā)生失水事故之后，破口處會(huì)有冷卻劑泄漏到安全殼當(dāng)中，在低壓狀態(tài)下閃蒸，利用冷凝機(jī)制、對流機(jī)制等，向鋼安全殼壁傳遞熱量，然后經(jīng)過壁面向外側(cè)流動(dòng)液膜傳遞。液膜受熱蒸發(fā)，會(huì)形成向上流動(dòng)的水蒸汽，從而從安全殼頂部開口出帶走熱量。

1 冷凝模型

在冷凝模型當(dāng)中，2D矩形通道是其計(jì)算區(qū)域，對低雷諾數(shù)下的NAVIER-STOKES方程進(jìn)行計(jì)算，能夠得到相應(yīng)的質(zhì)量守恒方程、動(dòng)量守恒方程、能量守恒方程、組分傳輸方程等。分別采用了Buleev湍流模型和Chien k-ε湍流模型[1]。其中，Buleev湍流模型具有良好的魯棒性，Chien k-ε湍流模型則具有良好的準(zhǔn)確性。通過COPAIN冷凝試驗(yàn)對冷凝模型進(jìn)行數(shù)據(jù)驗(yàn)證。結(jié)合實(shí)際情況，對79個(gè)驗(yàn)證工況進(jìn)行計(jì)算，其中，利用兩種模型分別計(jì)算了59個(gè)工況和20個(gè)工況。最終得出的總冷凝質(zhì)量流率計(jì)算值與試驗(yàn)值之間，存在著15%以內(nèi)的偏差，證明了冷凝模型具有良好的準(zhǔn)確性，能夠用于非能動(dòng)安全殼冷卻系統(tǒng)CFD模型的建立。

2 蒸發(fā)模型

在蒸發(fā)模型當(dāng)中，針對非能動(dòng)安全殼冷卻系統(tǒng)中，鋼安全殼外側(cè)的空氣環(huán)腔，利用G.Aiello進(jìn)行了相應(yīng)的1D模型的開發(fā)。在模型當(dāng)中，對液膜、空氣之間的質(zhì)量平衡、能量平衡等進(jìn)行了兼顧，具體分為蒸發(fā)傳熱、對流傳熱、傳導(dǎo)傳熱等不同的液膜主要傳熱機(jī)理。在模型當(dāng)中，主要是通過濕空氣能量平衡方程、液膜質(zhì)量平衡方程、濕空氣動(dòng)量平衡方程等進(jìn)行求解，同時(shí)將總熱流密度方程、對流熱流密度、液膜蒸發(fā)質(zhì)量流密度、液膜蒸發(fā)熱流密度等作為求解的輔助方程[2]。

如果工況中存在無液膜覆蓋的干壁狀態(tài)，在模型當(dāng)中，利用CF

X-4程度進(jìn)行計(jì)算，并且驗(yàn)證得到的結(jié)果，從而證明了模型與程序的計(jì)算結(jié)果比較相符。如果工況是有液膜覆蓋的濕壁狀態(tài)，可利用Park和Kang的蒸發(fā)試驗(yàn)進(jìn)行驗(yàn)證，計(jì)算數(shù)值采用干壁工況狀態(tài)下的標(biāo)定值。通過計(jì)算，能夠得出在軸向平均傳熱系數(shù)中，隨著入口空氣流速變化曲線和試驗(yàn)值的對比情況，從而得出在蒸發(fā)模型中，計(jì)算值和試驗(yàn)值比較相符。由此可見，蒸發(fā)模型在干壁工況和濕壁工況下，都具有較高的準(zhǔn)確度，因而能夠用于非動(dòng)能安全殼冷卻系統(tǒng)CFD模型的建立。

3 非能動(dòng)安全科冷凝系統(tǒng)CFD模型

根據(jù)前面開發(fā)完成的冷凝模型和蒸發(fā)模型，可以進(jìn)行非動(dòng)能安全殼冷卻系統(tǒng)CFD模型的建立。在該CFD模型當(dāng)中，根據(jù)邊界條件、幾何參數(shù)等相關(guān)的參數(shù)，例如蒸發(fā)通道的高度、寬度、壓力、入口溫度；冷凝板的高度、厚度；冷凝通道的高度、寬度、壓力、入口溫度、入口流速、質(zhì)量分?jǐn)?shù)等，就能夠?qū)Ψ莿?dòng)能安全殼冷卻系統(tǒng)的運(yùn)行特性進(jìn)行反映。在冷凝板處，空氣和蒸汽的混合物會(huì)發(fā)生冷凝，蒸汽的質(zhì)量分?jǐn)?shù)、近壁面的溫度等都會(huì)有所下降[3]。會(huì)朝著壁面方向，產(chǎn)生冷凝質(zhì)量流率。通過冷凝、對流等方式，熱量能夠傳遞到壁面，進(jìn)而繼續(xù)向壁面外側(cè)的液膜進(jìn)行傳遞。液膜在吸收熱量之后，會(huì)蒸發(fā)產(chǎn)生蒸發(fā)質(zhì)量流密度，進(jìn)而通過與空氣的對流和蒸發(fā)進(jìn)行熱量傳遞?？諝馐軣嵘仙?，會(huì)在向上流動(dòng)的過程中，蒸汽質(zhì)量分?jǐn)?shù)、空氣溫度等，也會(huì)逐漸升高。

在穩(wěn)態(tài)下，CFD模型中的參數(shù)分布情況，能夠反映出非能動(dòng)安全殼冷卻系統(tǒng)的變化趨勢，結(jié)果表明其與計(jì)算結(jié)構(gòu)基本吻合?？偀崃髅芏仍诶淠鍍蓚?cè)具有一定的不同，尤其是在兩側(cè)通道的出入口位置，具有更為明顯的差別。這主要是因?yàn)樵谕ǖ廊肟谔?，具有入口效?yīng)，同時(shí)受到第一個(gè)三維計(jì)算網(wǎng)格處理方式的影響造成的[4]。經(jīng)過驗(yàn)證，表明總傳熱量在冷凝板兩側(cè)是相同的。因此，非能動(dòng)安全殼冷卻系統(tǒng)CFD冷凝模型和蒸發(fā)模型能夠?qū)ζ浞€(wěn)態(tài)下的傳熱特性進(jìn)行有效的模擬，從而掌握其運(yùn)行狀態(tài)，確保其發(fā)揮出良好的作用和效果。

4 結(jié)束語

非動(dòng)能安全殼冷卻系統(tǒng)在實(shí)際應(yīng)用中具有十分重要的作用，為了確保其能夠發(fā)揮出更為良好的性能和效果，采用CFD來對冷卻系統(tǒng)的冷凝模型和蒸發(fā)模型進(jìn)行開發(fā)，并且以此為基礎(chǔ)建立非動(dòng)能安全殼冷卻系統(tǒng)CFD模型。經(jīng)過試驗(yàn)和研究證明，該模型具有較為良好的準(zhǔn)確性，因而能夠在實(shí)際應(yīng)用中發(fā)揮出良好的作用。

參考文獻(xiàn)

[1]蔣孝蔚，余紅星，孫玉發(fā)，等.非能動(dòng)安全殼冷卻系統(tǒng)傳熱傳質(zhì)模型研究[J].核動(dòng)力工程，2013，5：118-123.

[2]葉成，鄭明光，王勇，等.AP1000非能動(dòng)安全殼冷卻水WGOTHIC分析[J].原子能科學(xué)技術(shù)，2013，12：2225-2230.

[3]魏巍，魏來，陳艷芳，等.基于SimCont的AP1000非能動(dòng)安全殼冷卻系統(tǒng)仿真實(shí)現(xiàn)[J]. 原子能科學(xué)技術(shù)，2013，5：794-797.

[4]孫流莉，陽祥，趙瑞昌，等.安全殼冷卻試驗(yàn)蒸發(fā)特征數(shù)不確定度評(píng)定方法及其應(yīng)用[J].動(dòng)力工程學(xué)報(bào)，2015，8：687-692.

作者簡介：李懿（1991-），女，籍貫：內(nèi)蒙古，學(xué)歷：本科，職稱：研究生，研究方向：智能控制技術(shù)、多場耦合機(jī)電系統(tǒng)動(dòng)力學(xué)分析。

武建新（1971-），男，籍貫：內(nèi)蒙古，學(xué)位：博士，職稱：教授，研究方向：機(jī)器人技術(shù)、智能控制技術(shù)、多場耦合機(jī)電系統(tǒng)動(dòng)力學(xué)分析、機(jī)電系統(tǒng)優(yōu)化。