設(shè)計(jì)-行波堆物理特性分析

1、摘要行波堆是一種新概念反應(yīng)堆，它依靠初始裝載的易裂變核（如 U235）產(chǎn)生的中子將可裂變核形成增殖波先行，燃耗（如 U238）轉(zhuǎn)化為易裂變核（Pu239）的特性，行的長的反應(yīng)堆堆型。行波堆運(yùn)行時(shí)不需要控制棒等反應(yīng)性控制裝置，僅依靠反應(yīng)堆自身調(diào)節(jié)其反應(yīng)性，達(dá)到臨界的目的，也不需要冷卻劑流量的控制，在反應(yīng)堆壽期內(nèi)不需換料并能達(dá)到很高的燃耗深度，從而減少對核廢料的處理，并能更進(jìn)一步地防止核擴(kuò)散，減少廢料處理，達(dá)到了核能發(fā)電的綠色可持續(xù)需求、高固有安全需求、高經(jīng)濟(jì)性需求和防核擴(kuò)散需求，符合中國的能源需求和環(huán)保策略，有較好的預(yù)期前景。本文首先對行波堆的發(fā)展歷史進(jìn)行了系統(tǒng)調(diào)研，了解幾種主流的行波堆的概念設(shè)

2、計(jì)和運(yùn)行原理，以此為基礎(chǔ)建立了一個(gè)簡單的行波堆模型（該模型為圓柱形堆芯，直徑 65cm，高度 200cm，堆芯分為點(diǎn)火區(qū)和增殖區(qū)，堆芯外面 著冷卻劑層，厚 10cm，冷卻劑層外面為反射層，厚 15cm）。其次，學(xué)習(xí)了采用程序進(jìn)行反應(yīng)堆分析的方法，并借助開發(fā)的先進(jìn)程序 RMC 對提出的行波堆模型進(jìn)行了臨界計(jì)算和燃耗計(jì)算，分析了 Keff、中子通量和功率分布、主要錒系核素等物理參數(shù)隨燃耗的變化趨勢。結(jié)果表明所設(shè)計(jì)的行波堆方案在燃耗過程中，呈現(xiàn)出了標(biāo)準(zhǔn)的通量波、功率波和密度波，符合通常意義上行波堆的燃耗特性。本文的工作為將來行波堆的概念優(yōu)化、物理分析和工程應(yīng)用做了初步準(zhǔn)備和技術(shù)積累，在基礎(chǔ)層面上具有

3、一定價(jià)值。：行波堆，蒙卡方法，RMC，物理特性，燃耗計(jì)算IABSTRACTTraveling wave reactor heap is a new concept, it depends on the initial loading ofneutrons of fissile nuclear fuel (such as U235) to be generated by the fisofnuclear fuel (such as U238)o fissile fuel (Pu239) characteristic,wave formingproliferation, reaction burn

4、up wave underwent piles of long life type. Traveling wavereactor control rods, etc. do not need to run the reactivity control devi, rely solelyon self-regulation of its reactivity reactor reached critical pures, does not need tocontrol the flow of coolanthe reactor lifetime without refueling and can

5、 reach highburnup, thereby reducing the nuclear waste treatment and prevention of nuclearproliferation can further reduce waste disal, reaching the green sustainable demandfor nuclearer generation, high inherent safety requirements demand higheconomic and anti-nuclear Diffudemand, in line with China

6、s energy needs andenvironmental strategies, expected bettrospects.ly, the development history of the traveling wave reactor systematically research to understand the concepts of design and operating principles of severalmainstream traveling wave reactor, as a basis for the establishment of a simplet

7、raveling wave reactor m(This mis a cylindrical core, the diameter of 65cm,height of 200cm, and the proliferation of the seedo the core region, the coolant outof the core layer is wrapped withickness of 10cm, outside the coolant layer is areflective layer, a thickness of 15cm). Secondly, the study of

8、 the Monte Carloprocedure reactorysis method, and with Tsinghua University to develop advancedRMC Monte Carlo program for the traveling wave reactor mproed by thecritical computing and burnup calculations,er distribution, mainly actinides and othyzes Keff, neutron flux volume andhysical parameters o

9、f the trend withburnup. The results showt the design of the traveling wave reactor programburnup pros, showing a standard flux wave, waveer and density waves, in linewith the usual sense of upward wave reactor burnup characteristics. Future work ofthis pr is the concept of the traveling wave reactor

10、 optimization, physicalysisand preliminary engineering done to prepare and technology accumulation, with acertain valuehe basic research level.KEYWORDS:TWR,Monte Carlo Methods,RMC,Physical Characteristics,BurnupII目錄摘要IABSTRACTII第 1 章1.11.21.31.4第 2 章2.12.22.3緒論1選題背景和意義1國內(nèi)外現(xiàn)狀3行波堆的應(yīng)用前景4的主要內(nèi)容及預(yù)期目標(biāo)6蒙卡方法

11、與 RMC 程序7蒙卡方法介紹7RMC 簡介9RMC 輸入概述102.3.12.3.22.3.3輸入模塊10輸入卡10輸入格式112.4RMC 特殊設(shè)置122.4.12.4.22.4.32.4.42.4.5幾何12材料12臨界計(jì)算12計(jì)數(shù)器12燃耗計(jì)算12本章小結(jié)13第 3 章 行波堆設(shè)計(jì)方案14行波堆模型及參數(shù)14RMC 建模15本章小結(jié)16第 4 章 計(jì)算結(jié)果分析174.14.24.3keff 分析17通量、功率分布分析17核素密度變化分析19本章小結(jié)22第 5 章 結(jié)論及展望23參考文獻(xiàn)25附致錄26謝371第 1 章緒論1.1 選題背景和意義地球上的全部的能源總的來說可以分為兩類，即核

12、能和引力能兩種物理機(jī)制，核能分為裂變能與聚變能，聚變能來自，在地球上出現(xiàn)的形式豐富多彩，所有的生物生長都依賴于光，生物通過光合作用，把能轉(zhuǎn)化為生物質(zhì)能，經(jīng)過數(shù)百萬年在的變化，經(jīng)過化學(xué)和物理、熱化學(xué)、生物化學(xué)等過程，將生物質(zhì)能轉(zhuǎn)變?yōu)榛渲邪禾?、石油、天然氣等。如果是由能光伏發(fā)電形式收集，就是現(xiàn)在發(fā)展的光點(diǎn)，如果以集熱的方式所使用，就是現(xiàn)在的太陽能發(fā)電。這些都是人為直接轉(zhuǎn)化運(yùn)用的形式。以上所述這些基于來自能核聚變在地球上的表現(xiàn)，因?yàn)槠溲h(huán)周期比較短，很快被被大家稱為可再生能源。能所恢復(fù)，所以真正基于地球自身的物理機(jī)制是核裂變能，核裂變能產(chǎn)生熱量，是自然的體現(xiàn)；而在人工應(yīng)用方面的體現(xiàn)，則表現(xiàn)

13、為通常所說的核能，是通過人工核反應(yīng)堆體現(xiàn)出來的，具體表現(xiàn)為通過反應(yīng)堆核裂變產(chǎn)生熱量，利用這些熱量產(chǎn)生蒸汽，進(jìn)而蒸汽用來發(fā)電。核裂變能的特殊性，使它直接來源于這些核裂變能的能量密度是化學(xué)能的百萬倍，以此特性打造的核裂變能的應(yīng)用體系，顯然是與其他能源系特和不同的。由于核裂變能的利用需要濃縮鈾或天然鈾作為，通過產(chǎn)生蒸汽來發(fā)電，因此核裂變能的大規(guī)模發(fā)展受到了濃縮鈾產(chǎn)能和循環(huán)發(fā)展的制約。也就是說，核能源規(guī)模的擴(kuò)張將增加對技術(shù)創(chuàng)新的相應(yīng)要求。現(xiàn)在開放和循環(huán)技術(shù)，可以短期和中長期核電規(guī)模不太大的滿足核能應(yīng)用的需求，大規(guī)模的擴(kuò)張是不現(xiàn)實(shí)的，而且從資源利用、廢物處理這些角度，都需要封閉式的循環(huán)系統(tǒng)以及相應(yīng)的核反

14、應(yīng)堆堆型，比如使用混合快中子增值反應(yīng)堆等。系統(tǒng)的先進(jìn)反應(yīng)堆、高溫氣冷堆、基于以上原因行波式核反應(yīng)堆的概念被提出。行波反應(yīng)堆不同于當(dāng)前正在運(yùn)行的商業(yè)核反應(yīng)堆堆芯，通過操作合理分布，核就可以從高濃縮鈾的點(diǎn)火區(qū)點(diǎn)燃，U235 裂變產(chǎn)生的多余中子將周圍增殖區(qū)不易裂變的 U238 轉(zhuǎn)化成易裂變的 Pu239，當(dāng)達(dá)到一定的濃度之后，keff 穩(wěn)定在 1 左右，形成自持裂變反應(yīng)，同時(shí)開始在原位生成的，形成行波。行波堆以增殖波先行續(xù)得形式運(yùn)行，裝料足夠的話可以連續(xù)運(yùn)行數(shù)十年甚至上百年。為了維持核反應(yīng)堆的持續(xù)運(yùn)行，堆芯的區(qū)保持常規(guī)核反應(yīng)堆的大小和質(zhì)量，按照常規(guī)的方1圖 1-1 行波堆燃燒策略式利用冷卻劑將核裂變

15、產(chǎn)生的熱量帶出堆芯，利用蒸汽發(fā)生器產(chǎn)生的蒸汽用來發(fā)電，其余的部分為燒盡或待增殖的。除了啟動(dòng)時(shí)點(diǎn)火區(qū)需要用到濃縮鈾，其他所有的都可以使用貧鈾或天然的鈾礦，因此不需要鈾的分離和濃縮。形象的說，行波堆就像蠟燭，用火柴點(diǎn)燃后就會(huì)逐漸燒盡，而且可以點(diǎn)燃其他的蠟燭。在實(shí)際工程化應(yīng)用過程中，行波堆可以選擇使用定期移動(dòng)，增殖波的空間的固定處理方式。如圖 1-1 所示，這就是行波堆的基本概念。行波堆工業(yè)系統(tǒng)不需要大規(guī)模的，可以逐步減少濃縮鈾的產(chǎn)品，大大縮減了昂貴、復(fù)雜的處理體系。系統(tǒng)包括多個(gè)步驟，從鈾的開采、鈾的濃縮、轉(zhuǎn)化、終端處理等。利用行波反應(yīng)堆，可以大大簡化了步驟的數(shù)目，可以逐步地減少負(fù)極的產(chǎn)能以及復(fù)雜昂

16、貴的系統(tǒng)。因此行波堆的發(fā)展不但可以提早進(jìn)入核能的大規(guī)?？沙掷m(xù)發(fā)展，同時(shí)還極大地減小了核擴(kuò)散的風(fēng)險(xiǎn)。比如鈾的濃縮和是核擴(kuò)散風(fēng)險(xiǎn)最大的環(huán)節(jié)。的后續(xù)處理，因此，與其他核系統(tǒng)相比，行波反應(yīng)堆更可能贏得大多數(shù)國家的支持，在世界范圍廣泛開展真正體現(xiàn)核能與煤、油、氣等化石相比，具有無碳排放、清潔的優(yōu)點(diǎn)。因此，行波堆的發(fā)展能滿足核能發(fā)電的綠色可持續(xù)需求、高固有安全需求、高經(jīng)濟(jì)性需求和防核擴(kuò)散需求。從長久使用核能資源的角度，行波堆和快堆的作用是類似的。他們有很高的利用資源，如果開發(fā)成功，能讓人類用上千百年。從核的利用形式上來看，行波堆使用貧鈾作為運(yùn)行期間不換料，沒有乏。從核不擴(kuò)散角度來看，行波堆由于使用貧鈾后處

17、理或者只是簡單后處理。，行波堆一次裝料能夠持續(xù)運(yùn)行幾十年，除最初的點(diǎn)火區(qū)必須要少量濃縮鈾外，2其余所有都可以從天然鈾、貧鈾或乏中選材。行波燃燒理論從根本上改革了快堆的轉(zhuǎn)化和模式，它使用了易裂變核素的生成與消耗之間的自平衡設(shè)計(jì)，能夠讓反應(yīng)堆數(shù)十年穩(wěn)態(tài)運(yùn)行而無需燃燒控制，安全高效地轉(zhuǎn)化和利用鈾-钚資源，處理輕水堆的乏。在資源利用、防核擴(kuò)散、環(huán)境保護(hù)、經(jīng)濟(jì)和安全方面都具有突出的優(yōu)點(diǎn)，具有全面第四代反應(yīng)堆相關(guān)性能水平的性能行波堆的這種特性，能夠?qū)F(xiàn)有的沉重廢物負(fù)擔(dān)轉(zhuǎn)化成巨大的經(jīng)濟(jì)價(jià)值。據(jù)估算現(xiàn)存的 70 萬噸核廢料可產(chǎn)生 100的電供給本國使用。中國的鈾、數(shù)十年復(fù)雜昂貴的處理體系和核天然鈾在行波堆的發(fā)

18、展構(gòu)建上，可能避免乏擴(kuò)散的漫長歷史，主動(dòng)進(jìn)入核能利用更為先進(jìn)的時(shí)代，最大程度保護(hù)人類、國家、環(huán)境利益式創(chuàng)新?？偨Y(jié)起來，行波堆技術(shù)可以簡單地概括為核增殖，是一個(gè)先，行進(jìn)的理想狀態(tài)的東西，滿足可持續(xù)性、防核擴(kuò)散、高經(jīng)濟(jì)性和安全波堆可以讓鈾資源的利用率提高近百倍，廢物量減少數(shù)十倍，是一種可以把一個(gè)百年的資源為數(shù)千年的技術(shù)。但是就目前的發(fā)展情況來看，雖然行波堆作為一種新概念的快堆，在原理上是可行的，符合中國的環(huán)保策略和能源需求，有較好的預(yù)期前景，但因?yàn)槟透咻椛鋼p傷材料和高燃耗等相關(guān)技術(shù)問題，其在工程可實(shí)施性方面仍1.2 國內(nèi)外巨大。現(xiàn)狀行波堆的物理概念源于上世紀(jì)四十年代，學(xué)術(shù)界有一些相關(guān)性的學(xué)術(shù)，但由

19、于缺乏計(jì)算工具，材料和要求高于已知范疇，致使工作無法深入到工程化應(yīng)用開波快堆是一種原位增殖的快堆。它能夠通過一次通過循環(huán)方式對天然鈾資源實(shí)現(xiàn)高效利用，是核能系統(tǒng)上的性的創(chuàng)新。行波反應(yīng)堆的尚處于起步階段，不僅缺乏工程實(shí)踐，在設(shè)計(jì)上也僅處于堆芯概念設(shè)計(jì)階段1-3。九十年代氫彈之父 Edwardler 與其Lowell Wood 開始把行波堆作為未來大規(guī)模利用的能源系統(tǒng)進(jìn)行，逐步在國際上影響帶動(dòng)了一些機(jī)構(gòu)的超前學(xué)術(shù)，譬如 2000 年 Huan Dam 對“自穩(wěn)臨界波”進(jìn)行了基本分析和計(jì)算 ，Xue-Nong Chen 等對“孤立燃耗波”的物理求解等4-6。H.Sekimoto 教授他們的團(tuán)隊(duì)也做了

20、大量的基礎(chǔ)CANDLE 堆的概念和行波堆的概念非常相似，工作。H.Sekimoto 教授主要是研發(fā)第四代先進(jìn)反應(yīng)堆，主攻鉛鉍合金冷卻快堆，包括模塊式小型堆，作為能源使用和核廢物嬗變應(yīng)用解決 Po 污染和清除技術(shù)；2003 年-2008 年提出并開發(fā) CANDLE（蠟燭堆）行波式增殖、核反應(yīng)堆概念和技術(shù)。他CANDLE 是一種新的核of Neutron Flux, Nuclide反應(yīng)堆的燃耗策略，縮寫代表 Constant Axial Sh3Densities ander ShDuring Life of Energy Production（恒定的軸向的中子通量形狀，核子密度和功率形狀在能量產(chǎn)生

21、的時(shí)間內(nèi)），而且也代表著蠟燭般的燃耗。他表示，當(dāng)此燃耗策略被采用時(shí)，雖然被固定在反應(yīng)堆堆芯，但是隨著功率的輸出，燃燒區(qū)域按照成比例的速度移動(dòng)，沿著堆芯軸向，不改變核子數(shù)密度的空間分布，中子通量和功率密度。反應(yīng)性和反應(yīng)堆特性不會(huì)改變。最顯著，使用這種策略時(shí)，它是沒有必要使用控制棒進(jìn)行燃耗的控制。因此，蠟燭核反應(yīng)堆是安全的，而同樣重要的是，讓感到更安全。蠟燭燃耗有各種其他突破性的優(yōu)點(diǎn)。當(dāng)此燃耗作用在快堆，具有優(yōu)異的中子經(jīng)濟(jì)效益，可以得到優(yōu)異的性能。它可以使用天然鈾或貧作為并且大約的 40會(huì)被燃燒。就可以繼續(xù)使用大量貧鈾已經(jīng)可用，因此，如果能夠使用它作為，核能幾乎是千年而無需進(jìn)一步的鈾礦開采，濃縮，

22、和后處理。另外，使用的的量被大大降低。同時(shí)他還，雖然在使用蠟燭燃耗很大的優(yōu)勢，但是可以使用之前眾多的技術(shù)發(fā)展是必要的。然而，對于塊高溫氣冷反應(yīng)堆，目前幾個(gè)國家正在開發(fā)，蠟燭燃耗可無需額外的技術(shù)開發(fā)應(yīng)用。他還描述了蠟燭燃耗應(yīng)用到高溫氣冷反應(yīng)堆和快堆會(huì)具有優(yōu)良的中子經(jīng)濟(jì)7-14。除了行波堆的學(xué)術(shù)進(jìn)展之外，2006 年“發(fā)明”創(chuàng)立“泰拉能源”公司對行波堆進(jìn)行工程化研發(fā)，該公司蘭在美核學(xué)會(huì) 2009 年“核”九的封面中較為詳細(xì)地介紹了“行波堆”和公司發(fā)展情況。相比之下，國內(nèi)對行波堆的并不是很充分，在物理分析、數(shù)值模擬、工程設(shè)計(jì)等各方面的都很少，對這種新型堆的可行性和物理特性的了解并不是很深。國家核電技

23、術(shù)公司與泰拉能源公司簽署了技術(shù)交流合作諒解備忘錄。國家核電核工程究。設(shè)計(jì)院和泰拉能源合作，共同推出了中國項(xiàng)目初步可行性研1.3 行波堆的應(yīng)用前景行波堆(Traveling Wave Reactor)是完全自動(dòng)化的長壽期的反應(yīng)堆模型。它的基本是用一個(gè)相對較小的臨界系統(tǒng)來形成一個(gè)臨界波，點(diǎn)燃相對較大的次臨界系統(tǒng)。假設(shè)有一個(gè)易裂變核素核密度低，但可裂變核素核子密度相對較豐富的圓柱形次臨界系統(tǒng)，通過在其一端添加外（這個(gè)可以是器產(chǎn)生的中子也一個(gè)臨界的系統(tǒng)產(chǎn)生的中子），可以讓這一端的可裂變核素轉(zhuǎn)化為易裂變核素，keff 穩(wěn)定在 1 后去掉外，反應(yīng)堆仍可以自行增殖，形成一個(gè)穩(wěn)定向前推進(jìn)的臨界波。當(dāng)兩種效應(yīng)相

24、互抵消，臨界波能繼續(xù)增殖，形成增殖波先行，燃耗行的現(xiàn)象，這就是行波堆的基本概念。這個(gè)臨界波有一些特性，在波的前端，可裂變核素在中子的作用下轉(zhuǎn)化為易裂變核素，keff 上升，在波的后面，轉(zhuǎn)化的易裂變核素由于消耗而減少，且生成了大量的裂變產(chǎn)物，4導(dǎo)致 Keff 下降。行波反應(yīng)堆燃耗策略具有以下優(yōu)點(diǎn)：反應(yīng)堆運(yùn)行時(shí)無需反應(yīng)性控制裝置：反應(yīng)堆控制更方便；無過剩反應(yīng)性；由反應(yīng)性引起的事故可以避免；沒有控制棒吸收中子，因此中子利用率更高。隨著燃耗反應(yīng)堆特性不會(huì)發(fā)生變化：如功率峰因子、keff 等不隨燃耗發(fā)生變化；不同燃耗期反應(yīng)堆操作策略不變。不需要控制冷卻劑的流量：徑向各點(diǎn)的功率特性不隨燃耗的變化而變化，冷

25、卻劑的流量也不需要改變。徑向功率分布的優(yōu)化更簡單：由于徑向每一點(diǎn)功率不隨燃耗發(fā)生變化，只需對其中一個(gè)燃耗點(diǎn)的徑向功率分布進(jìn)行優(yōu)化，而不需要考慮燃耗的影響。所以，優(yōu)化方案更簡單、更徹底。能夠有很大的效果來防止核的擴(kuò)散，削減廢料處理方面有很好的效果：增加堆芯的高度還可延長反應(yīng)堆的壽期，因此可以設(shè)計(jì)長的反應(yīng)堆。在反應(yīng)堆運(yùn)行期間內(nèi)不需要更換而且能達(dá)到很高的燃耗深度，不僅減少了核廢料的處理，而且進(jìn)一步地防止了核擴(kuò)散。6） 除了燃耗區(qū)以外，的 K都低于 1，發(fā)生反應(yīng)性事故的幾率很低，的和更加簡便、安全。缺點(diǎn)包括以下三點(diǎn)：1）材料問題：包殼、結(jié)構(gòu)材料的壽期問題。2）沿軸向的冷卻劑壓降增大：由于冷卻劑壓降方面

26、的限制，堆芯軸向高度不能太高。3）堆芯軸向的功率進(jìn)行優(yōu)化的空間很?。憾研镜妮S向的功率波以及功率峰值都是行波堆運(yùn)行期間不可避免的問題。要求設(shè)計(jì)時(shí)需要針對功率峰因子進(jìn)行優(yōu)化設(shè)計(jì)。行波堆還可以用壓水堆運(yùn)行產(chǎn)生的的廢料作為總結(jié)如下：進(jìn)行燃燒，它的優(yōu)點(diǎn)可以1）第一個(gè)循環(huán)后不再需要富集的 U235。U238 在中子的作用下可以轉(zhuǎn)化出大量的易裂變的 Pu239，可以維持反應(yīng)堆的臨界，而且其反應(yīng)性幾乎等于零，行波堆停堆時(shí)只要引入少量的負(fù)反應(yīng)性，反應(yīng)堆就會(huì)停止運(yùn)行；如果需要再次進(jìn)行啟動(dòng)，引入即可啟動(dòng)。燃耗深度可達(dá) 40%（不同2）行波堆壽期比較長，燃燒后卸載出來的堆型卸載出來的燃耗深度也會(huì)不相同）。中國核工業(yè)已

27、經(jīng)發(fā)展了很多年，并主要堆型壓水堆核電站運(yùn)營多年，積累了大量的核廢料，如果這些廢物在行波反應(yīng)器的應(yīng)用，不僅變廢為寶，產(chǎn)生大量的5電能，而且還可以大大減少核廢料的數(shù)量，降低了加工難度。1.4的主要內(nèi)容及預(yù)期目標(biāo)本文首先對行波堆的發(fā)展歷史進(jìn)行了系統(tǒng)調(diào)研，了解幾種主流的行波堆的概念設(shè)計(jì)和運(yùn)行原理，以此為基礎(chǔ)建立了一個(gè)簡單的行波堆模型。其次，學(xué)習(xí)了采用序 RMC程序進(jìn)行反應(yīng)堆分析的方法，并借助開發(fā)的先進(jìn)程對行波堆模型進(jìn)行了臨界計(jì)算和燃耗計(jì)算，分析了 Keff、中子通量和功率分布、主要錒系核素等物理參數(shù)隨燃耗的變化趨勢。第一章主要通過參考一些國內(nèi)外的相關(guān)文獻(xiàn)對行波堆的選題背景意義和國內(nèi)外現(xiàn)狀以及應(yīng)用前景進(jìn)

28、行系統(tǒng)的了調(diào)研，并總結(jié)了調(diào)研的資料進(jìn)行具體的描述，最后介紹本文的主要內(nèi)容及預(yù)期目標(biāo)。第二章主要介紹了蒙卡方法的相關(guān)內(nèi)容，詳細(xì)介紹了堆用蒙卡程序 RMC 程序，包括 RMC 程序的輸入概述和一些特殊的設(shè)置。第三章對行波堆的設(shè)計(jì)方案進(jìn)行了完整的描述，參照文獻(xiàn)一的方案，確定行波堆的模型，作出基本圖形，并確定行波堆的物理參數(shù)和相關(guān)材料，然后用 RMC程序建立行波堆模型。第四章主要是利用堆用蒙卡 RMC 程序建立的行波堆模型進(jìn)行了行波堆臨界計(jì)算和燃耗計(jì)算，對 RMC 程序計(jì)算輸出的數(shù)據(jù)進(jìn)行分析，作出 keff、通量波、功率波以及相關(guān)核素密度變化的圖形，分析 Keff、中子通量和功率分布、主要錒系核素等物

29、理參數(shù)隨燃耗的變化趨勢。第五章主要是對全文的總結(jié)和對行波堆未來發(fā)展的展望，提出了本文波堆方案的不足，并為方案優(yōu)化給出一些建議。行6第 2 章蒙卡方法與 RMC 程序方法，也被稱為統(tǒng)計(jì)模擬方法，是 21 世紀(jì) 40 年代中期，由于科學(xué)技術(shù)和計(jì)算機(jī)的發(fā)明的開發(fā)，而們一種以概率統(tǒng)計(jì)理論為指導(dǎo)的一類非常重要的數(shù)值計(jì)算方法。由于能解粒子輸運(yùn)問題， 方法被認(rèn)為是反應(yīng) 堆堆芯設(shè)計(jì)的先進(jìn)方法，因此本文進(jìn)行行波堆相關(guān)計(jì)算就是利用了堆用蒙卡方法 RMC 程序，本章主要簡單介紹了蒙卡方法，然后對 RMC 程序的進(jìn)行了詳細(xì)的介紹，包括 RMC 輸入概述以及一些特殊設(shè)置15。2.1 蒙卡方法介紹中子輸運(yùn)理論的基本原理是

30、中子在介質(zhì)中運(yùn)動(dòng)的和的進(jìn)程。隨著核能科學(xué)技術(shù)發(fā)展，它已經(jīng)成為了核科學(xué)技術(shù)領(lǐng)域中一個(gè)獨(dú)立的重要基礎(chǔ)理論科學(xué)，在核反應(yīng)堆物理與、核技術(shù)工程應(yīng)用以及軍事等領(lǐng)域中得到了廣泛應(yīng)用。蒙特卡羅方法便是目前中子輸運(yùn)的重要方法之一。方法也稱為統(tǒng)計(jì)模擬法、隨機(jī)模擬法。二十世紀(jì)四十年代中期隨著科學(xué)技術(shù)的發(fā)展和電子計(jì)算機(jī)的發(fā)明，人們提出了一種以概率統(tǒng)計(jì)理論為基礎(chǔ)的一類非常重要的數(shù)值計(jì)算方法。由在第二次中研制的“曼計(jì)劃”計(jì)劃的成員 S.M.烏和 J.馮首先提出。 他們用這種概率論方法在計(jì)算機(jī)上模擬粒子鏈?zhǔn)椒磻?yīng)。通過對中子運(yùn)動(dòng)歷史的統(tǒng)計(jì)，得到期望參數(shù)的估計(jì)值。由于這個(gè)方法使用了隨機(jī)數(shù)和抽樣，J.馮等人將其以著名通常人們所

31、用的方法可以大致地分成兩類：一類是所求解的問題本身具有內(nèi)在的隨機(jī)性，借助計(jì)算機(jī)的運(yùn)算能力可以直接模擬這種隨機(jī)的過程。例如在核物理中，分析中子在反應(yīng)堆中的傳輸過程。中子與原子核作用受到量子力學(xué)規(guī)律的制約，人們只能知道它們相互作用發(fā)生的概率，卻無法準(zhǔn)確獲得中 子與原子核作用時(shí)的位置以及裂變產(chǎn)生的新中子的行進(jìn)速率和方向?？茖W(xué)家依據(jù)其概率進(jìn)行隨機(jī)抽樣得到裂變位置、速度和方向，這樣模擬大量中子的行為后，經(jīng)過統(tǒng)計(jì)就能獲得中子傳輸?shù)姆秶?，作為反?yīng)堆設(shè)計(jì)的依據(jù)。另一種類型是解決問題可以轉(zhuǎn)化為一個(gè)隨機(jī)分布的數(shù)的特征，如發(fā)生隨機(jī)概率，或隨量的期望值。通過隨機(jī)抽樣的方法，隨機(jī)事件估計(jì)概率，或取樣頻率估計(jì)的隨量的功能

32、特性的數(shù)量，而作為這個(gè)問題的一個(gè)解決方案。該方法用于解決復(fù)雜問題的積分。如下圖 3-1 所示為蒙卡方法解決中子與原子核相互作用的隨機(jī)抽樣框圖：7賭城（Monte Carlo）來命名，從此稱為了隨機(jī)模擬的代名詞。圖 2-1 隨機(jī)抽樣框圖粒子輸運(yùn)的問題具有顯著隨機(jī)性質(zhì)，粒子輸運(yùn)過程是一種 過程。一 個(gè)從源出發(fā)的粒子，在其運(yùn)動(dòng)方向上哪一點(diǎn)碰撞是偶然的，但是有一定的概率分布；和原子核（或原子）碰撞，還有幾種概率不全相同的類型；碰撞之后的能量和方向，也遵從了一定的概率分布。在這個(gè)過程中，粒子可能被介質(zhì)吸收，也可能從系統(tǒng)逃脫，運(yùn)動(dòng)過程結(jié)束；否則，粒子將繼續(xù)下一次類似的運(yùn)動(dòng)過程。但是下一次碰撞的位置、碰撞后

33、的能力和方向只與本次碰撞有關(guān)，而與粒子之前的碰撞無關(guān)。在求解粒子輸運(yùn)問題方面，首先要確定粒子輸運(yùn)的所有物理過程，包括介質(zhì)的原子核組成、粒子發(fā)射的分布、粒子的源分布、與原子核發(fā)生各種反應(yīng)的還要考慮粒子 t 和權(quán)重 W。假設(shè)（不帶電）粒子在兩次碰撞之中按直線運(yùn)動(dòng)，運(yùn)動(dòng)方向與能量均沒有改變，那么粒子在介質(zhì)中的運(yùn)動(dòng)過程可用下面狀態(tài)序列描述：S0 , S1, S2 ,.,SM 1, SM（2-1）這個(gè)序列稱為粒子隨機(jī)運(yùn)動(dòng)歷史，模擬一個(gè)粒子的運(yùn)動(dòng)過程轉(zhuǎn)變成確定狀態(tài)序列的過程。模擬粒子的運(yùn)動(dòng)過程分為兩個(gè)步驟：第一步確定粒子的初始狀態(tài) S0 ，第二步確認(rèn)由 Sm1 到 Sm （m=1,2,M）。第二步由可以分

34、為兩個(gè)過程：第一個(gè)過程是確定碰撞點(diǎn)的位置 Zm ，稱為輸運(yùn)過程；第二個(gè)過程是確定碰撞后粒子的運(yùn)動(dòng)方向和能量，稱為碰撞過程。粒子輸運(yùn)計(jì)算的主要問題有:(1)臨界問題：臨界問題的是中子增殖系統(tǒng)的中子平衡分布，實(shí)質(zhì)上有效增殖因數(shù)的求解是積分方程的特征值問題。用8概率以及粒子運(yùn)動(dòng)歷史結(jié)束的條件等，之后確定粒子的狀態(tài)序列和狀態(tài)參數(shù)，狀態(tài)參數(shù)包括粒子的能力 E，空間位置 r 和運(yùn)動(dòng)方向 ，用 S r( E, , ) 表示，有時(shí)方法求解臨界問題，和問題模擬中子運(yùn)動(dòng)歷史基本相同，不同點(diǎn)在于需要通過源迭代得到穩(wěn)定的分布；（2）問題：在輻射源周圍，放置具有一定形狀和厚度的物，求被后射線的能譜、強(qiáng)度和角分布，蒙特卡

35、羅方法求解問題就是解具有穩(wěn)定外源項(xiàng)的非積分方程；（3）微擾問題：微擾問題計(jì)算反應(yīng)堆的微小變化引起的對有效增殖因數(shù)的改變量；（4）通量問題： 通量是粒子輸運(yùn)問題中的基本物理量。用蒙卡方法計(jì)算通量，是粒子輸運(yùn)問題的基本問題；（5）伴隨問題：從數(shù)學(xué)方面看，粒子輸運(yùn)問題的解與伴隨問題的解有一定的關(guān)系，可以通過解粒子輸運(yùn)方程的伴隨問題來求原來問題的解15-19。由于能解決粒子輸運(yùn)問題方法被認(rèn)為是反應(yīng)堆堆芯設(shè)計(jì)的先進(jìn)方法。隨著計(jì)算機(jī)的飛速發(fā)展，全堆(蒙卡)計(jì)算已經(jīng)成為一個(gè)熱點(diǎn)。建立完善的通用程序，以避免大量的重復(fù)性工作，可以是程序，基于蒙特卡羅方法開展技術(shù)和成果的改進(jìn)和修正的計(jì)算的基礎(chǔ)，這些結(jié)果可以進(jìn)一步

36、提高程序。通用程序通常具有以下特點(diǎn)：具有靈活的幾何處理能力參數(shù)通用化，使用方便元素和介質(zhì)材料數(shù)據(jù)齊全能量范圍廣，功能強(qiáng)，輸出量靈活全面含有簡單可靠又能普遍適用的抽樣技巧具有較強(qiáng)的繪圖功能14RMC 簡介2.29堆用蒙卡分析程序 RMC(Reactor Monte Carlo code)是由工程物理系核能科學(xué)與工程管理反應(yīng)堆工程計(jì)算分析（REAL 團(tuán)隊(duì)）研發(fā)的、用于反應(yīng)堆堆芯計(jì)算分析的三維輸運(yùn)蒙卡程序。RMC 的研發(fā)始于 2001 年，至今經(jīng)歷了三個(gè)階段，已投入人力大于 40 人年。RMC 程序針對反應(yīng)堆計(jì)算分析中的基本需求，同時(shí)結(jié)合新概念反應(yīng)堆系統(tǒng)設(shè)計(jì)時(shí)幾何結(jié)構(gòu)靈活、中子能譜復(fù)雜及材料組分多樣

37、、各向異性及泄漏強(qiáng)（某些特定情況）等特點(diǎn)進(jìn)行研發(fā)，是多物理多尺寸耦合核能系統(tǒng)數(shù)值分析 的物理計(jì)算 。RMC 能夠處理復(fù)雜幾何結(jié)構(gòu)、采用連續(xù)能量點(diǎn)截面對復(fù)雜能譜和材料進(jìn)行描述，并能夠根據(jù)實(shí)際問題的需要對臨界問題本征值、本征函數(shù)計(jì)算，系統(tǒng)燃耗模擬，瞬態(tài)過程分析（含 截面處理）等進(jìn)行計(jì)算。并針對蒙卡方法的特點(diǎn)， RMC 中研發(fā)并應(yīng)用了幾何處理 、核截面處理優(yōu)化、輸運(yùn)過程模擬新方法、源收斂判斷與 、計(jì)數(shù)器優(yōu)化、并行算法優(yōu)化等提高計(jì)算效率的方法和技巧。RMC 可以處理包括重復(fù)結(jié)構(gòu)在內(nèi)的任意幾何結(jié)構(gòu),采用了 能量格架以及哈希表等一些新的算法對能量和截面搜索進(jìn)行改善。此外,RMC 具有豐富的2.3 RMC

38、輸入概述2.3.1 輸入模塊RMC 的輸入文件按照模塊劃分，各模塊的名稱及相應(yīng)功能如下所示：SURFACE 模塊：定義曲面類型和曲面方程。UNIVERSE 模塊：描述某個(gè)完整的幾何空間。RMC 采用基于層級空間的幾何描述，輸入文件中可能存在多個(gè) UNIVERSE 模塊。MATERIAL 模塊：定義材料組成。CRITICALITY 模塊：定義臨界計(jì)算參數(shù)，包括粒子數(shù)、初始源等。TALLY 模塊：定義計(jì)數(shù)器，包括通量、功率、反應(yīng)率等。CONVERGENCE 模塊：定義源收斂和參數(shù)。BURNUP 模塊：定義燃耗計(jì)算參數(shù)，包括燃耗柵元、功率、時(shí)間步長等。PR模塊：定義輸出內(nèi)容。PLOT 模塊：定義畫圖

39、參數(shù)。2.3.2 輸入卡各個(gè)輸入模塊具有特定的輸入卡：表 2-1 輸入模塊10輸入模塊輸入卡功能描述定義一個(gè)曲面。包括曲面類型，曲面方程參數(shù)，曲面邊界條件，等SurfaurfUniverseCell定義空間內(nèi)的某個(gè)柵元。包括填充材料，幾何形狀，體積，溫度，等MaterialMat定義一種材料。包括材料密度，核素組成，等CriticalityerIter指定源迭代的初始 keff和粒子數(shù)InitSrc指定源迭代計(jì)算的初始裂變源分布收斂技巧,很好地解決了通用蒙卡程序(MCNP)收斂慢的問題,并增加了溫度相關(guān)截面處理、燃耗計(jì)算、中子動(dòng)力學(xué)計(jì)算、混合方法、熱工耦合等功能。續(xù)表 2-1輸入模塊2.3.3

40、 輸入格式RMC 輸入文件的格式應(yīng)注意以下幾點(diǎn)：（1） 每個(gè)模塊以相應(yīng)的標(biāo)識(shí)，模塊之間以空行隔開。形如表 2-2 輸入格式（2）輸入卡頂格寫，輸入卡中的選項(xiàng)卡以空格間隔。如果輸入卡一行未寫完，可換行后空格續(xù)寫。例如：注釋符使用“/”。RMC 輸入文件不區(qū)分大小寫。在 windows 下，不建議使用 txt 格式的文本文件作為輸入文件。建議使用 UltraEdit 轉(zhuǎn)換為 Dos 格式。11CellTally 2type = 1filter = 101energy = 6.25E-720.0cell = 2 0 1:289Universe0 Universe1 Surface Material

41、CriticalityTallyCellTally定義柵元計(jì)數(shù)器。統(tǒng)計(jì)某個(gè)或多個(gè)柵元內(nèi)的積分通量、功率、吸收反應(yīng)率或裂變反應(yīng)率BurnupBurnCell指定燃耗區(qū)er指定總功率TimeStep指定時(shí)間步長AceLib指定重要核素所匹配的 ACE 截面數(shù)據(jù)庫Outpell指定需要輸出核素密度的柵元2.4 RMC 特殊設(shè)置2.4.1 幾何復(fù)雜幾何描述是蒙卡程序相對于確定論程序的重要優(yōu)勢之一。與世界上其它絕大多數(shù)蒙卡程序一樣，RMC 采用基于層級空間的幾何描述系（universe-based geometry system） 。RMC 的幾何描述系統(tǒng)包括三類基本的幾何描述單元，即，曲面（surfa

42、ce），柵元（cell）和空間（universe）。一般地，物理系統(tǒng)由多個(gè)或單個(gè)層級空間組成，每個(gè)空間由一定數(shù)量的柵元算來定義。，柵元通過曲面方向 （sense） 的交并運(yùn)2.4.2 材料材料輸入模塊描述材料組成，包括材料密度、核素截面數(shù)據(jù)庫、核素份額，以及定義與連續(xù)能量或多群 ACE 截面數(shù)據(jù)庫相關(guān)的參數(shù)選項(xiàng)。2.4.3 臨界計(jì)算蒙卡臨界計(jì)算采用裂變源迭代法，即，用當(dāng)前代產(chǎn)生的裂變源作為下一代的初始源。用戶需要在輸入文件中指定用于源迭代的基本參數(shù)，包括初始有效增殖因子（keff）、每代中子數(shù)和中子代數(shù)，同時(shí)還需要指定初始裂變源分布。在臨界計(jì)算模塊中，用戶還可以選擇隨機(jī)數(shù)發(fā)生器的類型和參數(shù)，以

43、及選擇并行臨界計(jì)算的并行模式。2.4.4 計(jì)數(shù)器RMC 包含三類計(jì)數(shù)器，分別是柵元計(jì)數(shù)器（cell tally）、網(wǎng)格計(jì)數(shù)器（mesh tally）和截面計(jì)數(shù)器（cross-section tally）。柵元計(jì)數(shù)器用來統(tǒng)計(jì)柵元內(nèi)的宏觀物理量，包括積分通量、功率、裂變反應(yīng)率和吸收反應(yīng)率等。RMC 采用 Cell-Map方法，能夠高效地處理大規(guī)模柵元計(jì)數(shù)器。網(wǎng)格計(jì)數(shù)器與柵元計(jì)數(shù)器類似，但它統(tǒng)計(jì)的宏觀物理量不是基于柵元，而是基于預(yù)先劃定的網(wǎng)格。截面計(jì)數(shù)器統(tǒng)計(jì)柵元內(nèi)指定核素、指定反應(yīng)類型的反應(yīng)截面（微觀反應(yīng)率） 。以上三種計(jì)數(shù)器都支持分能群統(tǒng)計(jì)，即，按照指定的能量區(qū)間分別計(jì)數(shù)。2.4.5 燃耗計(jì)算蒙卡

44、燃耗計(jì)算是蒙卡臨界計(jì)算和點(diǎn)燃耗計(jì)算的相互耦合。傳統(tǒng)的蒙卡燃耗程序（如 MCBurn、MCODE）一般采用第接口，通過外耦合的方式，循環(huán)調(diào)用蒙卡輸運(yùn)程序（如 MCNP）和點(diǎn)燃耗程序（如 ORIGEN-2）。12在 RMC 程序當(dāng)中，內(nèi)嵌有開發(fā)的點(diǎn)燃耗計(jì)算模塊 DEPTH。DEPTH 采用矩陣指數(shù)法，能夠精確、高效地處理含約 1500 種核素的精細(xì)燃耗鏈。RMC 燃耗計(jì)算的基本流程是，首先通過臨界計(jì)算（連續(xù)能量）模塊，得到中子通量、單群反應(yīng)截面等數(shù)據(jù)，傳遞給點(diǎn)燃耗模塊 DEPTH；DEPTH 完成點(diǎn)燃耗計(jì)算，得到新的核素密度，傳遞給 RMC 臨界計(jì)算模塊。通過數(shù)據(jù)的往返傳遞，從而完成燃耗計(jì)算的全過

45、程。與傳統(tǒng)的蒙卡燃耗程序（如 MCBurn，MCODE）相比，RMC 燃耗計(jì)算功能的基本特點(diǎn)包括：含內(nèi)耦合燃耗計(jì)算模塊，能夠處理含約 1500 種核素的精細(xì)燃耗鏈；整合了 ORIGEN-S 和 ORIGEN-2 的最新點(diǎn)燃耗數(shù)據(jù)庫。支持含重復(fù)結(jié)構(gòu)的燃耗計(jì)算，無需用戶為每個(gè)燃耗區(qū)指定不同的初始材料，大大減少了用戶輸入。支持對大規(guī)模燃耗區(qū)的“并行臨界+并行點(diǎn)燃耗”計(jì)算。在并行點(diǎn)燃耗計(jì)算模式中，燃耗區(qū)被平均分配到各個(gè)進(jìn)程，各自獨(dú)立地完成點(diǎn)燃耗計(jì)算20。本章小結(jié)本章主要介紹了蒙卡方法與 RMC 程序。對蒙卡方法的介紹中首先提到隨著核能科學(xué)技術(shù)發(fā)展，中子輸運(yùn)理論已經(jīng)成為了核科學(xué)技術(shù)領(lǐng)域中一個(gè)獨(dú)立的重要基

46、礎(chǔ)理論科學(xué)，方法便是目前中子輸運(yùn)的重要方法之一，隨后介紹了蒙卡方法的背景以及分類，最后詳細(xì)講述了蒙卡方法解決粒子輸運(yùn)問題。對 RMC程序的介紹首先講述 RMC 的來源是由研發(fā)的，然后簡單介紹了 RMC 相比蒙卡程序在解決反應(yīng)堆問題方面的各種優(yōu)點(diǎn)，最后詳細(xì)講解了這些優(yōu)點(diǎn)。13第 3 章行波堆設(shè)計(jì)方案本文主要行波堆物理特性，首先要做的就是建立行波堆模型，而建立行波堆模型要設(shè)計(jì)行波堆的參數(shù)與材料，本文借鑒文獻(xiàn)一（，. 行波堆初步概念設(shè)計(jì)（2010）的行波堆設(shè)計(jì)方案，利用 RMC程序建立行波堆模型，作出行波堆模型圖，列出行波堆相關(guān)參數(shù)，并介紹了輸入卡的書寫。3.1行波堆模型及參數(shù)圖 3-1 行波堆模型

47、表 3-1 堆芯物理參數(shù)14堆芯物理參數(shù)值堆芯直徑65cm堆芯高度200cm冷卻劑厚度10cm反射層厚度15cm點(diǎn)火區(qū)高度20cm增殖區(qū)高度180cm圖 3-1 為堆芯分布示意圖，表 3-1 為堆芯各項(xiàng)物理參數(shù)，圖中上端的紅域表示 U235 的較高富集度區(qū)，是行波堆的點(diǎn)火區(qū)，長度為 20cm；下邊黃域?yàn)樾胁ǘ训脑鲋硡^(qū)，其中 U235 富集度較低，但 U238 富集度較高，長度為 180cm；第二層的區(qū)域?yàn)殂U鉍冷卻劑層，第三層區(qū)域?yàn)槭瓷鋵訁^(qū)域。行波堆運(yùn)行時(shí)從上邊點(diǎn)火區(qū)啟動(dòng)，點(diǎn)燃下邊增殖區(qū)，當(dāng) keff 變?yōu)?1 的時(shí)候，燃燒區(qū)域沿著堆芯軸向按照成比例的速度移動(dòng)，這就是行波堆的燃燒模式。下表

48、3-2 為行波堆堆芯所用材料：表 3-2 堆芯材料3.2 RMC 建模參照文獻(xiàn)一（，. 行波堆初步概念設(shè)計(jì)（2010）的方案，利用 RMC 程序建立行波堆模型，分析 Keff、功率分布、增值比、主要核素密度等關(guān)鍵物理參數(shù)的變化規(guī)律，進(jìn)行燃耗運(yùn)算，最后分析行波堆相關(guān)物理特性。編輯行波堆輸入卡（見附錄），本文利用 RMC 程序設(shè)計(jì)的行波堆堆芯為圓柱形，垂直放置，堆芯高度微 200cm，直徑為 65cm，堆芯有 3 層：最內(nèi)層為區(qū)，區(qū)分為點(diǎn)火區(qū)和增殖區(qū)，上部 20cm 為點(diǎn)火區(qū)，剩下為增殖區(qū)，點(diǎn)火區(qū)為 U235百分比為 15%的金屬鈾混合物，增殖區(qū)為天然鈾；中間層為冷卻劑，冷卻劑材料為液態(tài)金屬鉛鉍合

49、金，冷卻劑厚度為 10cm；最外層為反射層，材料為石墨，厚度為 15cm，堆芯最外層表面設(shè)定為全反射邊界，并且中子沖出最外層表面后停止追蹤。為了便于統(tǒng)計(jì)堆芯通量和功率分布，堆芯區(qū)共分為 40 個(gè)柵元，點(diǎn)火區(qū)分為 4 個(gè)柵元， 增殖區(qū)分為 36 個(gè)柵元。首先用臨界計(jì)算輸入模塊（CRITICALITY）統(tǒng)計(jì) keff 變化規(guī)律，初始 keff 值設(shè)置為 1.0，初始裂變源設(shè)置在點(diǎn)火區(qū)的中心，每一代中子數(shù)為 50000，共計(jì)算 250 代，前 50 代為非活躍代；然后是用計(jì)數(shù)器輸入模塊(Tally)統(tǒng)計(jì)堆芯通量和功率的分布，本文中計(jì)數(shù)器輸入模塊用到的為柵元計(jì)數(shù)器；最后用燃耗計(jì)算模塊（BURNUP）進(jìn)

50、行燃耗計(jì)算,燃耗柵元為堆芯 40 個(gè)柵元，共 47 個(gè)時(shí)間步，步長Ace 截面數(shù)據(jù)庫為.30c。為天，堆芯功率為 600MW，15堆芯材料名稱冷卻劑材料Pb&Bi反射層材料C反應(yīng)堆功率600MW點(diǎn)火區(qū)高富集度鈾增殖區(qū)低富集度鈾輸入卡中 UNIVERSE 0 為行波堆堆芯柵元?jiǎng)澐帜K，CELL1-4 為點(diǎn)火區(qū)，CELL5-40 為增殖區(qū)，CELL41-43 為冷卻劑，CELL44-46 為反射層，CELL47 為堆芯外部空間；mat 表示柵元材料；void=1 表示中子進(jìn)入該區(qū)域停止。SURFACE 為曲面模塊，surf 為曲面輸入卡，定義柵元曲面；cz 表示軸心在 Z 軸的曲面，pz 表示垂直

51、于 Z 軸的曲面；數(shù)字為曲面半徑；bc=1 表示全反射邊界條件。MATERIAL 為材料輸入模塊，定義各個(gè)柵元材料，數(shù)據(jù)庫為.30c。CRITICALITY 為臨界計(jì)算輸入模塊，erIter 為源迭代輸入卡；Keff0 選項(xiàng)卡設(shè)置初始有效增殖系數(shù)，缺省值為 keff0 = 1；Population 選項(xiàng)卡設(shè)置每代中子數(shù)（5000），非活躍代代數(shù)（50）和總代數(shù)（250）；InitSrc 為初始裂變源輸入卡的。Tally 為計(jì)數(shù)器輸入模塊，CellTally 為柵元計(jì)數(shù)器，Type=1 表示柵元通量， Type=2 表示柵元功率；Filter 表示計(jì)數(shù)柵元的層級；Cell 表示需要計(jì)數(shù)的柵元（1

52、 到 40）。BURNUP 表示燃耗計(jì)算模塊，BurnCell 為燃耗柵元；TimeStep 輸入卡指定各個(gè)燃耗步的時(shí)間步長，為天；er 輸入卡指定各個(gè)燃耗步的功率密度 ，為 W/gHM；SubStep 輸入卡指定點(diǎn)燃耗計(jì)算的步數(shù) （范圍 1 - 9999） ， 缺省值為 10 步；Inherent 輸入卡指定用于繼承重要核素的吸收份額和質(zhì)量份額，缺省值分別為 0.99 和 0.999；AceLib 輸入卡指定繼承的核素所采用的數(shù)據(jù)庫（.30c）；Outp的文件中。本章小結(jié)ell 卡用于輸出指定柵元的核素密度，存于后綴名為“.den”本章主要內(nèi)容是借鑒文獻(xiàn)一（，. 行波堆初步概念設(shè)計(jì)（2010

53、）的行波堆設(shè)計(jì)方案，首先作出了行波堆模型圖，設(shè)計(jì)行波堆相關(guān)物理參數(shù)和材料，根據(jù)這些參數(shù)建立 RMC 程序行波堆模型，對模型進(jìn)行了簡單的介紹，并詳細(xì)介紹輸入卡的書寫。利用 RMC 程序計(jì)算編寫的輸入卡，為分析行波堆相關(guān)物理特性進(jìn)行準(zhǔn)備。16第 4 章計(jì)算結(jié)果分析本章的主要內(nèi)容是利用 RMC 程序，計(jì)算編輯的行波堆輸入卡，對該堆芯進(jìn)行數(shù)值計(jì)算，根據(jù)計(jì)算出的數(shù)據(jù)分析 keff 變化、中子通量分布、堆芯功率分布和 U235、U238、Pu239、Pu240、Pu241 的核子密度變化，并根據(jù)計(jì)算數(shù)據(jù)畫出相應(yīng)圖形，對行波堆能否實(shí)現(xiàn)自動(dòng)化的長壽期運(yùn)行進(jìn)行驗(yàn)證。4.1 keff 分析利用 RMC 程序，對該

54、堆芯進(jìn)行數(shù)值計(jì)算，計(jì)算結(jié)果：1.41.31.21.11.00.90.80.70.60.501234567891011燃耗步（步數(shù)）圖 4-1 堆芯 Keff 隨時(shí)間變化圖 4-1 表示堆芯 keff 隨時(shí)間的變化。從圖中可以看出來，在開始的時(shí)間內(nèi)（前幾個(gè)燃耗步），keff 的值因?yàn)?U235 的大量裂變而比較大，當(dāng) Pu239 的產(chǎn)生率和消耗率相等的時(shí)候，keff 逐漸下降達(dá)到穩(wěn)定值，直到功率波移動(dòng)到行波堆的尾部時(shí)，由于耗盡，keff 下降。4.2 通量、功率分布分析利用 RMC 程序，對該堆芯進(jìn)行數(shù)值計(jì)算，計(jì)算結(jié)果：17keff4035302520151050-501020柵元3040圖 4

55、-2 堆芯功率波隨時(shí)間變化108642001020柵元3040圖 4-3 堆芯通量波隨時(shí)間變化率分布在經(jīng)過初始的功率波動(dòng)下降之后，形成了沿堆芯軸在圖 4-2向的穩(wěn)定功率波。由于上端點(diǎn)火區(qū)初始的 U235 濃度很高，所以圖 5 中左端功率出現(xiàn)峰值。而緊接著功率下降是由于點(diǎn)火區(qū) U235 的消耗導(dǎo)致的，與此同時(shí)在堆芯功率波的右側(cè) U238 在中子的作用下轉(zhuǎn)化為 Pu239，導(dǎo)致功率向右偏移。這樣在功率波的右端 U238 轉(zhuǎn)化為 Pu239，功率上升，在功率波的左端由于 Pu239 的18通量（n/cm2s）功率（MW）2天12天36天60天70天180天360天360天360天360天2天12天3

56、6天60天70天180天360天360天360天360天消耗和裂變產(chǎn)物的生成，功率下降，從而形成了自左向右的功率波，也正是由于U238 轉(zhuǎn)化為 Pu239，維持了功率的穩(wěn)定。4.3 核素密度變化分析行波堆是完全自動(dòng)化的長壽期反應(yīng)堆堆型，因此堆芯內(nèi)不需要控制棒、可溶性毒物等反應(yīng)性控制裝置，需要靠堆芯自身的特性實(shí)現(xiàn)其完全自動(dòng)化。圖 4-4、圖 4-5、圖 4-6、圖 4-7、圖 4-8 分別是行波堆設(shè)計(jì)方案中 U235、U238、Pu239、Pu240 和 Pu241 在堆內(nèi)隨時(shí)間的分布曲線。圖 4-4 U235 核子密度分布隨時(shí)間變化圖 4-4 表示位于堆芯點(diǎn)火區(qū)的 U235 隨時(shí)間的變化。從圖

57、中可以看出，在初始時(shí)刻（前幾個(gè)燃耗步），U235 富集度比較大，當(dāng)反應(yīng)堆開始運(yùn)行的時(shí)候，U235裂變，產(chǎn)生大量中子，裂變產(chǎn)生的多余中子將周圍不易裂變的 U238 轉(zhuǎn)化成易裂19變的 Pu239，經(jīng)過幾個(gè)燃耗步之后，U235 迅速減少，隨著反應(yīng)堆的運(yùn)行 U235 的密度逐漸達(dá)到穩(wěn)定的分布，直到功率波移動(dòng)到行波堆的尾部時(shí)，隨著的消耗耗，U235 的密度開始下降，耗盡時(shí)，U235 密度下降到幾乎為 0。圖 4-5U238 核子密度分布隨時(shí)間變化圖 4-5 表示堆芯的 U238 隨時(shí)間的變化。從圖中可以看出，在初始時(shí)刻（前幾個(gè)燃耗步），U238 的變化非常小，當(dāng)反應(yīng)堆開始運(yùn)行的時(shí)候，U235 裂變，產(chǎn)

58、生大量中子，裂變產(chǎn)生的多余中子將周圍不易裂變的 U238 轉(zhuǎn)化成易裂變的 Pu239，經(jīng)過前幾個(gè)燃耗步之后，臨界區(qū)的 U238 開始減少，而臨界區(qū)后面的 U238 密度保持不變，隨著反應(yīng)堆的運(yùn)行 U238 逐漸消耗，直到功率波移動(dòng)到行波堆的尾部時(shí)，U238 的含量迅速減少直到 耗盡。圖 4-6 Pu239 核子密度分布隨時(shí)間變化20圖 4-6 表示堆芯 Pu239 隨時(shí)間的變化。從圖中可以看出，在初始時(shí)刻（前幾個(gè)燃耗步），當(dāng)反應(yīng)堆開始運(yùn)行的時(shí)候，U235 裂變產(chǎn)生的大量中子將周圍不易裂變的U238 轉(zhuǎn)化成易裂變的Pu239，Pu239 含量迅速上升，隨著反應(yīng)堆的運(yùn)行，Pu239的含量逐漸下降，

59、隨著功率的波移動(dòng)，Pu239 的變化形成穩(wěn)定分布，直到行波堆的尾部時(shí)，耗盡 Pu239 含量下降。0.00100.00080.00060.00040.00020.000001020柵元3040圖 4-7Pu240 核子密度分布隨時(shí)間變化0.000250.000200.000150.000100.000050.0000001020柵元3040圖 4-8Pu241 核子密度分布隨時(shí)間變化圖 4-7、4-8 中可以看出，Pu240 和 Pu241 的核子密度在行波堆內(nèi)隨著反應(yīng)21數(shù)密度（n/cm3）數(shù)密度（n/cm3）2天12天36天60天70天180天360天360天360天360天2天12天36

60、天60天70天180天360天360天360天360天堆的運(yùn)行都會(huì)形成各自穩(wěn)定的分布。從本文以上的計(jì)算結(jié)果可以看出，目前該堆的設(shè)計(jì)存在幾個(gè)方面的問題：第一是行波堆在啟動(dòng)時(shí)，堆芯軸向功率的波動(dòng)太大，具體可見圖 5；第二是行波堆在啟動(dòng)時(shí) keff 低于 1.0，這種情況下，如果想讓反應(yīng)堆啟動(dòng)，必須加入外，考慮到提供外的時(shí)間不能太長（一般從幾天到幾周是可以接受的），最好的辦法是直接調(diào)整到 1 的 Keff 值啟動(dòng)；第三是行波堆啟動(dòng)時(shí) keff 的波動(dòng)太大。因?yàn)樾胁ǘ褯]有反應(yīng)性控制裝置，所以只能允許 keff 在 1.0 附近的一個(gè)可接受的范圍內(nèi)出現(xiàn)小幅的波動(dòng)。以上結(jié)果可以看出，行波堆在概念上是存在的