【摘要】：核反應堆嚴重事故進程中,核燃料元件會發(fā)生熔化,甚至落入下封頭。熔融物的遷移及再定位等現(xiàn)象關(guān)系到嚴重事故的發(fā)展與緩解。為了給燃料棒熔化的局部現(xiàn)象研究提供模擬分析工具,本文基于原始的MPS(Moving Particle Semi-implicit)方法進行了二次開發(fā),以使其具備模擬熔化現(xiàn)象的功能。本文對原始MPS方法的改進包括:1)增加了能量方程求解模塊,使其具備了傳熱計算功能;2)增加了表面張力計算模型,并評價了不同表面張力模型之間差異;3)引入了基于焓值變化的粘性模型,使得固體與流體之間的粘性能平滑過渡;4)改進了粘性項的時間差分格式,將原MPS法的顯式差分格式改為隱式,這樣增大了計算時間步長,能夠節(jié)省計算資源;5)引入了PMS(Passive Moving Solid)的剛體計算模型,使得流固耦合計算成為可能。為考驗修改后的程序,本文通過一維非穩(wěn)態(tài)導熱和方形液滴振蕩分別驗證了傳熱模型和表面張力模型。為了驗證程序的綜合性能,本文采用改進后三維MPS程序模擬了伍德合金的熔化過程,并與實驗進行了對比驗證,結(jié)果表明改進后的MPS程序能夠正確預測熔化行為�；趬核讶剂显奈镄缘恼{(diào)研結(jié)果,本文利用改進后的MPS方法對四棒束燃料棒的熔化進程進行了初步模擬。
【學位授予單位】：上海交通大學
【學位級別】：碩士
【學位授予年份】：2015
【分類號】：TL364.4
【圖文】：

圖 1-1 MELCOR 中的燃料元件熔化模型[19]Fig.1-1 Melting model of fuel rods in MELCOR[19]在國內(nèi)，開展嚴重事故模擬仿真分析時仍主要依賴國外的軟件SCDAP/RELAP5 和 MECLOR 等。近幾年來，國內(nèi)許多研究單位正在加緊開展電關(guān)鍵設計軟件自主化技術(shù)相關(guān)的研究工作，以實現(xiàn)核電軟件的國產(chǎn)化，使我核電軟件開發(fā)技術(shù)達到國際水平。譬如國家核電軟件中心正在開發(fā)的 COSIN程序，集成了核物理、反應堆熱工水利和系統(tǒng)安全分析等多種核電關(guān)鍵設計分程序，可以用于核電站的燃料設計及性能分析、嚴重事故模擬、概率安全分析等。另一方面，基于拉格朗日的粒子法，由于使用自由移動的粒子，不需要劃和構(gòu)建網(wǎng)格，處理簡單且計算精度高，在計算流體力學領(lǐng)域迅速引起了關(guān)注。動粒子半隱式方法，即 MPS（Moving Particle Semi-implicit）方法作為拉格朗粒子法的一種，1997 年由日本東京大學 S.Koshizuka 教授提出，最初用于解決

圖 2-2 梯度模型示意圖[28]Fig. 2-2 Particle interaction for gradient model[28]斯算子模型 ，擴散方程如下：2dd t 擴散系數(shù)，在每個時間步長 t內(nèi)， 增加2d t，d 為子 i 與粒子 j 之間的遷移量為：02( )j ii j id tW r rn 入是為了使作用模型模擬結(jié)果與擴散方程的解析解相

界條件由表面面通常存在較大的變形，由于網(wǎng)格畸變以及數(shù)值耗散等問題傳出滿意的結(jié)果，而在 MPS 方法中，自由面的判斷則相對很簡單子外部無粒子存在，因此其粒子數(shù)密度較小，如圖 2-3 所示。在密度滿足以下條件：* 0in n是表面粒子，其中β為小于 1 的參數(shù)，對自由表面的判別有一定ka 指出在β取值范圍在 0.8 到 0.99 時，較為合理有效[28]。當求 方程時，自由表面粒子壓力設為 0。

【參考文獻】

相關(guān)期刊論文 前10條

1 段軍;周紅;陶書生;周濤;;福島核事故后嚴重事故研究重點的探討[J];環(huán)境保護與循環(huán)經(jīng)濟;2014年06期

2 段廣濤;陳斌;;移動粒子半隱式算法的并行化[J];工程熱物理學報;2013年11期

3 孫中國;李帝辰;陳嘯;梁楊楊;席光;;移動粒子半隱式法在流體機械數(shù)值模擬中的應用[J];排灌機械工程學報;2013年11期

4 左娟莉;田文喜;秋穗正;陳榮華;蘇光輝;;液態(tài)金屬內(nèi)單個氣泡上升行為的MPS法數(shù)值模擬[J];原子能科學技術(shù);2011年12期

5 潘徐杰;張懷新;;移動粒子半隱式法的大渦模擬研究與應用[J];船舶力學;2011年07期

6 潘徐杰;張懷新;孫學堯;;一種適用于移動粒子半隱式法的自由表面識別方法[J];上海交通大學學報;2010年01期

7 吳鵬;楊燕華;林萌;;穩(wěn)壓器安全閥卡開導致的堆芯熔化事故[J];核動力工程;2008年04期

8 楊燕華;徐濟捚;;國外核電廠嚴重事故研究規(guī)模與現(xiàn)狀[J];核安全;2003年02期

9 付霄華;SCDAP/RELAP5與MELCOR程序?qū)Χ研緭p傷過程預測的比較[J];核動力工程;2003年05期

10 張鎖春;光滑質(zhì)點流體動力學（SPH）方法（綜述）[J];計算物理;1996年04期

本文編號：2743670