第4章-均勻反應(yīng)堆臨界理論2014

1、1均勻反應(yīng)堆的臨界理論均勻反應(yīng)堆的臨界理論主講：張競宇主講：張競宇2Contents前言前言均勻裸堆的單群擴散理論均勻裸堆的單群擴散理論有反射層反應(yīng)堆的單群擴散理論有反射層反應(yīng)堆的單群擴散理論功率分布展平概念功率分布展平概念3在上一章中我們討論了中子在在上一章中我們討論了中子在非增殖介質(zhì)非增殖介質(zhì)內(nèi)擴散內(nèi)擴散的規(guī)律和中子擴散方程的解法?，F(xiàn)在我們進一步的規(guī)律和中子擴散方程的解法?，F(xiàn)在我們進一步將其用于討論由核燃料和慢化劑等組成的將其用于討論由核燃料和慢化劑等組成的有限均有限均勻增殖介質(zhì)勻增殖介質(zhì)（反應(yīng)堆系統(tǒng)）內(nèi)的中子擴散問題。（反應(yīng)堆系統(tǒng)）內(nèi)的中子擴散問題。中心問題是討論反應(yīng)堆的臨界。中心問題是

2、討論反應(yīng)堆的臨界。4在反應(yīng)堆臨界理論中，主要研究兩方面的問題：在反應(yīng)堆臨界理論中，主要研究兩方面的問題：l 各種形狀的反應(yīng)堆達到臨界狀態(tài)的條件（臨界條件）各種形狀的反應(yīng)堆達到臨界狀態(tài)的條件（臨界條件）: e.g., 臨界時系統(tǒng)的體積大?。ㄅR界體積）和燃料成分（富臨界時系統(tǒng)的體積大小（臨界體積）和燃料成分（富集度）及其裝載量（臨界質(zhì)量）集度）及其裝載量（臨界質(zhì)量）。l 臨界狀態(tài)下系統(tǒng)內(nèi)中子通量密度（或功率）的空間分布。臨界狀態(tài)下系統(tǒng)內(nèi)中子通量密度（或功率）的空間分布。5確定輪方法計算實確定輪方法計算實際的反應(yīng)堆系統(tǒng)際的反應(yīng)堆系統(tǒng)l物理過程與中子能量的復雜物理過程與中子能量的復雜依賴關(guān)系依賴關(guān)系

3、“ “分群理論分群理論”l幾何的復雜性幾何的復雜性 “ “幾何近幾何近似似”處理處理l材料的復雜性材料的復雜性 “ “均勻化均勻化”處理處理 均勻反應(yīng)堆：均勻反應(yīng)堆：是指這樣一種堆，其堆芯的各種材料（燃料、是指這樣一種堆，其堆芯的各種材料（燃料、慢化劑、結(jié)構(gòu)材料等等）是均勻地混合在一起的。因此整慢化劑、結(jié)構(gòu)材料等等）是均勻地混合在一起的。因此整個堆芯的材料特性是一致的，核截面等數(shù)據(jù)都是一樣的。個堆芯的材料特性是一致的，核截面等數(shù)據(jù)都是一樣的。6均勻堆與非均勻堆均勻堆與非均勻堆l世界上數(shù)以千計的反應(yīng)堆中，只有一個名叫世界上數(shù)以千計的反應(yīng)堆中，只有一個名叫“水鍋爐水鍋爐”的實驗堆是均勻堆。其堆芯是

4、硫酸鈾酰的水溶液。的實驗堆是均勻堆。其堆芯是硫酸鈾酰的水溶液。l其他的都是非均勻堆，堆芯中的燃料和慢化劑是分開其他的都是非均勻堆，堆芯中的燃料和慢化劑是分開的，不混在一起。的，不混在一起。 既然如此，為何還要研究均勻反應(yīng)堆？既然如此，為何還要研究均勻反應(yīng)堆？Why?7研究思路：從容易的著手，逐步精確化研究思路：從容易的著手，逐步精確化均勻堆比較容易描述，建立的物理數(shù)學模型均勻堆比較容易描述，建立的物理數(shù)學模型比較簡單。但是，從中引出的基本概念有普遍比較簡單。但是，從中引出的基本概念有普遍應(yīng)用價值。應(yīng)用價值。工程設(shè)計中，對實際的非均勻反應(yīng)堆進行分析工程設(shè)計中，對實際的非均勻反應(yīng)堆進行分析時，也要

5、先進行時，也要先進行 “均勻化均勻化”，化為均勻堆。，化為均勻堆。8均勻裸堆的單群中子擴散方程的建立均勻裸堆的單群中子擴散方程的建立均勻裸堆的單群擴散方程的解均勻裸堆的單群擴散方程的解熱中子反應(yīng)堆的臨界條件熱中子反應(yīng)堆的臨界條件幾種幾何形狀裸堆的幾何曲率和中子通量密度幾種幾何形狀裸堆的幾何曲率和中子通量密度分布分布反應(yīng)堆曲率和臨界計算任務(wù)反應(yīng)堆曲率和臨界計算任務(wù)單群理論的修正單群理論的修正均勻裸堆的單群理論均勻裸堆的單群理論9裸堆裸堆：無反射層的反應(yīng)堆：無反射層的反應(yīng)堆單群單群：全部中子都在一個能群里。或者可以進一：全部中子都在一個能群里?；蛘呖梢赃M一步假設(shè)堆內(nèi)里所有中子都是熱中子。步假設(shè)堆內(nèi)

6、里所有中子都是熱中子。1 1. .單群理論的建立單群理論的建立10),(),(),(),(),(102trStrktrtrDttraa 對于由燃料與慢化劑組成的均勻增殖介質(zhì)反應(yīng)堆系統(tǒng)，單位時對于由燃料與慢化劑組成的均勻增殖介質(zhì)反應(yīng)堆系統(tǒng)，單位時間、單位體積內(nèi)的間、單位體積內(nèi)的裂變中子源強裂變中子源強為：為：),(),(trtrSfF 根據(jù)無限介質(zhì)增殖因子定義根據(jù)無限介質(zhì)增殖因子定義),(),(trktrSaF代入代入單群中子擴散方程單群中子擴散方程可得可得l D及及 a是對中子能譜平均后的數(shù)值；是對中子能譜平均后的數(shù)值；l 在反應(yīng)堆運行初期，須考慮外源中子，大多數(shù)情況下忽略外中子，認為在反應(yīng)堆

7、運行初期，須考慮外源中子，大多數(shù)情況下忽略外中子，認為裂變中子是反應(yīng)堆內(nèi)中子的唯一來源裂變中子是反應(yīng)堆內(nèi)中子的唯一來源?110),2(),2(tata221)()(1)()(LkdttdTtTDxx)()(),(tTxtx)()0 ,(0 xxa/2a/20 x無限平板反應(yīng)堆),(),(),(),(12txktxtxDttxaa(4-3)無外源無限平板反應(yīng)堆單群擴散方程無外源無限平板反應(yīng)堆單群擴散方程初始條件為初始條件為(4-4)邊界條件為邊界條件為(4-5),(1),(),(122txLktxttxD(4-6)由式(4-3)得利用分離變量法求解分離變量法求解，方程具有如下形式的解：(4-7)

8、將(4-7)式代入(4-6)式(4-8)2 2. .均勻裸堆的單群擴散方程的解均勻裸堆的單群擴散方程的解2aDL 12, 5 , 3 , 1nanBn, 3 , 2 , 1) 12(nanBnBxCBxAxsincos)(02cosBaA22)()(Bxx上式兩端必須等于某一常數(shù)，設(shè)為-B2，有0)()(22xBx或(4-9)波動方程(4-9)式的通解為由于初始通量密度分布0(x)關(guān)于x=0平面對稱，因此只能選擇滿足對稱條件的解，即BxAxcos)(由邊界條件(4-5)式可導出(x)滿足如下的邊界條件：(a/2)=0因此要求或(4-10)xanAxBAxnnnn) 12(coscos)(對于特

9、征值特征值Bn，波動方程的解 n(x) 稱為此問題的特征函數(shù)稱為此問題的特征函數(shù)。1322222/(1)(1)annnDLlDL BDL BnnltknCeT/) 1(nnnnlkdttdTtT1)()(1221)()(1nnnBLkdttdTtTDnnltknnexanAtx/)1(1) 12(cos),(對應(yīng)于每一個Bn值和n(x)，都有一個Tn(t)與之對應(yīng)該式可轉(zhuǎn)換為式中221nnBLkk(4-12)(4-13)(4-14)方程(4-12)解為其中C為待定常數(shù)。對于一維平板反應(yīng)堆，其中子通量密度的完全解中子通量密度的完全解就是對n=1到n=所有項的總和，即(4-15)14對于一定幾何形

10、狀和體積的反應(yīng)堆芯部，若對于一定幾何形狀和體積的反應(yīng)堆芯部，若B12對應(yīng)的對應(yīng)的k11，則，則(k1-1)為正值，中子通量密度為正值，中子通量密度 (x,t)將隨時間不斷增加，系將隨時間不斷增加，系統(tǒng)處于超臨界狀態(tài)。統(tǒng)處于超臨界狀態(tài)。若調(diào)整堆芯尺寸或改變材料成分，使若調(diào)整堆芯尺寸或改變材料成分，使k1 =1，則其余，則其余( kn-1)都將為負值。都將為負值。中子通量密度中子通量密度 (x,t)第一項將與時間無關(guān)，而其它各項將隨時間而衰減。第一項將與時間無關(guān)，而其它各項將隨時間而衰減。當時間足夠長時，當時間足夠長時，n1各項將衰減到零，系統(tǒng)處于穩(wěn)態(tài)，中子通量密度各項將衰減到零，系統(tǒng)處于穩(wěn)態(tài)，中

11、子通量密度按基波形式按基波形式(B=B1)分布，系統(tǒng)處于臨界狀態(tài)分布，系統(tǒng)處于臨界狀態(tài)。, 3 , 2 , 1) 12(nanBn221nnBLkknnltknnexanAtx/)1(1) 12(cos),(3 3. .熱中子反應(yīng)堆的臨界條件熱中子反應(yīng)堆的臨界條件三種情況：三種情況：150)()(22rBrg112121BLkk重要結(jié)論：重要結(jié)論：(1) 裸堆單群近似的臨界條件為：(4-17)B12為波動方程的最小特征值，記為Bg2，稱為特征曲率，稱為特征曲率；k1為有效增殖因子。(2) 反應(yīng)堆處于臨界狀態(tài)時，中子通量密度按最小特征值中子通量密度按最小特征值Bg2對應(yīng)的基波函對應(yīng)的基波函數(shù)分布

12、數(shù)分布，也就是說，此時反應(yīng)堆的中子通量密度空間分布滿足波動方程(4-18)16xaAxcos)(1)(1221aLkk無限平板反應(yīng)堆的臨界條件無限平板反應(yīng)堆的臨界條件為(4-19) 若系統(tǒng)材料組成給定若系統(tǒng)材料組成給定，則只有一個唯一的尺寸a0能使k1=1，即為臨界大??；當aa0時，則k11，為超臨界；當aa0時，k10，超臨界； 0，次臨界| |表示反應(yīng)堆偏離臨界狀態(tài)的程度。 反應(yīng)性單位：反應(yīng)性單位：PCM，1PCM=10-5；元，元，1元反應(yīng)性元反應(yīng)性=1 eff （反應(yīng)堆動力學）（反應(yīng)堆動力學）(4-49)33 單群是一種非常近似的方法。對于熱中子反應(yīng)堆，直接以上臨界條件有較大誤差單群是

13、一種非常近似的方法。對于熱中子反應(yīng)堆，直接以上臨界條件有較大誤差。 修正：修正：用用M2=L2+ 來替換上式中的來替換上式中的L211221gBLkk221LkBm11221gBMkk221MkBm熱中子反應(yīng)堆的 修正單群理論修正單群理論。例題：例題：P1006 6. .單群理論的修正單群理論的修正34反應(yīng)堆的最佳形狀反應(yīng)堆的最佳形狀這里所謂最佳最佳，是指用同樣材料，做成的反應(yīng)堆臨界體積最小。組成反應(yīng)堆的材料確定后，材料曲率就定組成反應(yīng)堆的材料確定后，材料曲率就定了，堆的臨界幾何曲率了，堆的臨界幾何曲率Bg也定了。但是，也定了。但是，對同一個幾何曲率對同一個幾何曲率Bg值，可以有不同的形值，可

14、以有不同的形狀和尺寸，對應(yīng)不同的體積狀和尺寸，對應(yīng)不同的體積。 35 2222222gxyzB B + B + Babcg22B()R222222.405()()grzBBBRH36 反應(yīng)堆的最佳形狀（續(xù)）幾何形狀最佳尺寸最小臨界體積長方體abc161/圓柱D1.08H148/球R130/ 3B3B3B37反射層的作用以及如何選擇反射層反射層的作用以及如何選擇反射層一側(cè)帶有反射層的反應(yīng)堆一側(cè)帶有反射層的反應(yīng)堆 帶有反射層的球形堆帶有反射層的球形堆 一側(cè)帶有反射層的柱形堆一側(cè)帶有反射層的柱形堆反射層節(jié)省反射層節(jié)省三、有反射層反應(yīng)堆的單群擴散理論三、有反射層反應(yīng)堆的單群擴散理論38 減少芯部中子的

15、泄漏，從而減小芯部的臨界體積，節(jié)省一減少芯部中子的泄漏，從而減小芯部的臨界體積，節(jié)省一部分核燃料；部分核燃料； 提高反應(yīng)堆的平均輸出功率，這是由于反射層的原因，其提高反應(yīng)堆的平均輸出功率，這是由于反射層的原因，其芯部中子通量密度分布比裸堆的中子通量密度分布更加平芯部中子通量密度分布比裸堆的中子通量密度分布更加平坦的緣故。坦的緣故。如何選擇反射層？如何選擇反射層？ 反射層材料散射截面要大反射層材料散射截面要大 反射層材料吸收截面要小反射層材料吸收截面要小 良好的慢化能力良好的慢化能力熱中子堆熱中子堆常用的反射層材料有：常用的反射層材料有：H2O, D2O, 石墨，鈹?shù)取Ｊ數(shù)?。快堆快堆常用?/p>

16、射層材料：常用反射層材料：重核重核 1 1. .反射層的作用以及如何選擇反射層反射層的作用以及如何選擇反射層391 裸堆裸堆 2 有反射層反應(yīng)堆有反射層反應(yīng)堆40221ccLkkB0)()()(2rkrrDcaccaccc芯部穩(wěn)態(tài)單群擴散方程芯部穩(wěn)態(tài)單群擴散方程 (角標角標 c)(4-52)u 該方程只有對于臨界系統(tǒng)臨界系統(tǒng)才成立。u 對于任意給定材料成分及幾何形狀與尺寸的反應(yīng)堆系統(tǒng)，它不一定處于穩(wěn)態(tài) 引入一個特征參數(shù)引入一個特征參數(shù)k來進行調(diào)整使其達到臨界。0)()()(2rkkrrDcaccaccc或者寫為0)()(22rBrcccLc 為芯部的擴散長度。(4-53)(4-54) 多區(qū)問題

17、多區(qū)問題: 分別寫出每一區(qū)的擴散方程分別寫出每一區(qū)的擴散方程證明證明: K為芯部的有效增殖因子為芯部的有效增殖因子2 2. .一側(cè)帶有反射層的反應(yīng)堆一側(cè)帶有反射層的反應(yīng)堆41221rrLk 0)()(22rkrrrr反射層穩(wěn)態(tài)單群擴散方程反射層穩(wěn)態(tài)單群擴散方程 （角標為（角標為 r）(4-56)式中(4-57)Lr為反射層的擴散長度。 邊界條件邊界條件為：(1) 在芯部和反射層的交界面上rcrrccDD(4-58)(2) 反射層的外推邊界上中子通量密度為零外推邊界上中子通量密度為零421. 帶有反射層的球形堆帶有反射層的球形堆（芯部半徑為（芯部半徑為R，反射層厚度為，反射層厚度為T）rrkAr

18、rkCrrrr)cosh()sinh()(rrTRkCrrr)(sinh)(rrBArcc)sin()(4-59)芯部方程式解芯部方程式解：反射層方程式解：反射層方程式解：(4-60)(tanhTRkCAr)sinh(1)cosh()sin()cos()sinh()sin(22TkRRTkkCDRRBRRBBADRTkCRRBArrrrccccrc根據(jù)：反射層外推邊界反射層外推邊界r=R+T處中子通量密度為零，有處中子通量密度為零，有上式代入(4-60)可得(4-61)方程(4-59)、(4-61)有兩個常數(shù)，由芯部與反射層交界面處邊界條件芯部與反射層交界面處邊界條件確定sinh x=e(x)

19、-e(-x)/2cosh x=ex+e(-x)/2tanh x=sinh x / cosh x431cot()1cothcccrrrRTDB RB RDLL帶反射層球形反應(yīng)堆單群臨界方程反射層球形反應(yīng)堆單群臨界方程給出了反應(yīng)堆臨界時幾何尺寸（給出了反應(yīng)堆臨界時幾何尺寸（R，T）和）和材料特性（材料特性（ Dc，Lr，Dr）之間的關(guān)系）之間的關(guān)系 當堆芯材料、反射層材料和尺寸已經(jīng)確定時，用此公式計算堆芯臨界尺寸的值。 當堆芯尺寸、反射層材料和尺寸已經(jīng)確定時，用它來計算達到臨界所需的堆芯材料成分 材料定，曲率知，算尺寸 尺寸知，算曲率，配材料。442. 側(cè)面帶有反射層的圓柱形反應(yīng)堆側(cè)面帶有反射層的

20、圓柱形反應(yīng)堆 (半徑半徑R，高度，高度H，側(cè)面反射層厚度，側(cè)面反射層厚度T)0)2,()2,(HrHrrc0),(zTRr;rrccrcDD0),(),(22zrBzrccc芯部反射層0),(),(22zrkzrrrr(4-63)(4-64)邊界條件邊界條件為在z=H/2處在r=R+T處(1) 在r=R處(4-65)(4-66)(4-67)45RR 0 芯部周圍有反射層芯部周圍有反射層 部分泄露出芯部的中子在反部分泄露出芯部的中子在反射層內(nèi)被散射而返回芯部，減少了中子損失，提高射層內(nèi)被散射而返回芯部，減少了中子損失，提高了中子的不泄露概率。了中子的不泄露概率。 因此在芯部材料性質(zhì)相同情況下，臨

21、界體積就要因此在芯部材料性質(zhì)相同情況下，臨界體積就要比裸堆的臨界體積小。比裸堆的臨界體積小。反射層節(jié)省反射層節(jié)省 ：反應(yīng)堆加上反射層所引起的臨界尺寸：反應(yīng)堆加上反射層所引起的臨界尺寸的減少。的減少。(4-81)2 2. .反射層節(jié)省反射層節(jié)省46反射層節(jié)省472222222405. 22405. 2effeffzrzrcHRHRBBB對于圓柱形反應(yīng)堆圓柱形反應(yīng)堆，反射層節(jié)省通常分別用徑向和軸向的反射層徑向和軸向的反射層節(jié)省來表示)22(,00HHRRzr(4-82) 可以把有可以把有反射層反應(yīng)堆的幾何曲率用芯部外形尺寸增大反射層反應(yīng)堆的幾何曲率用芯部外形尺寸增大 或或 2 來表示來表示。帶反射

22、層球形堆反射層球形堆22RBg圓柱形反應(yīng)堆圓柱形反應(yīng)堆R Reffeff、Heff稱為稱為等效半徑等效半徑、等效高度等效高度48)tanh(rrLTLrrrcccLTLRDRBRBDcoth1)cot(1rrccLTLBBtanh)tan(帶反射層球形反應(yīng)堆反射層節(jié)省反射層球形反應(yīng)堆反射層節(jié)省反射層厚度很小，即反射層厚度很小，即 T Lr = Lr。 過大的增加反射層厚度是沒有太大意義的！過大的增加反射層厚度是沒有太大意義的！rcrcLTRBLBcoth)cot(臨界方程臨界方程設(shè)Dc=Dr設(shè)反射層節(jié)省反射層節(jié)省 ，將R=R0- =/Bc- 代入上式得rrccLTLBBtanharctan1反

23、射層節(jié)省反射層節(jié)省 因為因為 很小，可得很小，可得tanhx = e(x)-e(-x) / ex+e(-x)49討論題討論題一個裸堆，加上反射層以后，幾何曲率與材料曲率還相等嗎？50討論題答案 不相等了。不相等了。 裸堆如果原先是臨界的，加上反射層以后，裸堆如果原先是臨界的，加上反射層以后，就超臨界了。就超臨界了。 如欲仍然維持臨界，只能如欲仍然維持臨界，只能 改變堆芯成分，降低材料曲率改變堆芯成分，降低材料曲率 減小堆芯尺寸，提高幾何曲率減小堆芯尺寸，提高幾何曲率51VHdVrVK)(1max 反應(yīng)堆內(nèi)的中子通量密度空間分布是不均勻的，而功率密度和中子通量反應(yīng)堆內(nèi)的中子通量密度空間分布是不均勻的