管理類(lèi)《核反應(yīng)堆物理》第3部分(反應(yīng)堆臨界理論)

第3部分反應(yīng)堆臨界理論核電廠管理類(lèi)員工核電基礎(chǔ)理論培訓(xùn)——《核反應(yīng)堆物理》主要內(nèi)容3.1均勻裸堆的單群理論3.2有反射層的均勻堆3.3非均勻反應(yīng)堆3.1均勻裸堆的單群理論均勻：由于嚴(yán)格按實(shí)際的非均勻堆進(jìn)行中子擴(kuò)散或輸運(yùn)方程求解，非常復(fù)雜。實(shí)際計(jì)算都作“均勻化”處理，即認(rèn)為：堆內(nèi)的燃料、慢化劑、冷卻劑及結(jié)構(gòu)材料等是均勻混合的。裸堆：沒(méi)有反射層的反應(yīng)堆；中子源：有增殖介質(zhì)；單群理論：把熱中子反應(yīng)堆內(nèi)的所有中子都看成是熱中子，忽略中子能量的影響。雙群理論：比單群更精確的模型，把熱中子劃為一群，快中子為一群?；靖拍顔稳簲U(kuò)散方程在由燃料-慢化劑構(gòu)成的有限大小的均勻裸堆系統(tǒng)的芯部，單位時(shí)間、單位體積內(nèi)產(chǎn)生的中子數(shù)為：得到：

根據(jù)無(wú)限介質(zhì)增殖系數(shù)的定義：考慮啟動(dòng)過(guò)程的獨(dú)立的外中子源并利用斐克定律，得單群擴(kuò)散方程：根據(jù)前面章節(jié)所得單群中子擴(kuò)散方程：無(wú)限平板堆的單群擴(kuò)散方程解用D除上式各項(xiàng)，并注意到L2=D/∑a，得到：這是一個(gè)二階偏微分方程，通常用分離變量法求解：左端是x的函數(shù)，右端是t的函數(shù)，兩端必須均等于某一常數(shù)，令為-B2,得方程組：擴(kuò)散方程：無(wú)限平板堆的單群擴(kuò)散方程解所得方程可改寫(xiě)為：容易得出其通解為：由于初始通量密度關(guān)于x=0平面對(duì)稱，所以：由邊界條件:要求:A不能為零,所以：或Bn2稱為特征值，對(duì)應(yīng)一系列滿足方程的特征函數(shù)。n=1的稱為基波本征函數(shù)，n>1的統(tǒng)稱為諧波本征函數(shù)。所以：典型的波動(dòng)方程！無(wú)限平板堆的單群擴(kuò)散方程解討論時(shí)間相關(guān)項(xiàng)方程的解：對(duì)應(yīng)一確定Bn，有一確定的Tn(t)，用L2/(1+Bn2L2)乘以上式，有：其中：方程的解：對(duì)應(yīng)每個(gè)n：是滿足方程的解，其線性組合也是原問(wèn)題的解：無(wú)限平板堆的單群擴(kuò)散方程解其中Cn和An為待定系數(shù)，利用余旋函數(shù)系的正交關(guān)系可求得：代入上式，得：隨時(shí)間變化項(xiàng)知:n=1時(shí)，B1最小，k1最大。由：無(wú)限平板堆的單群擴(kuò)散方程解（1）當(dāng)k1<1,所有kn-1都小于1，通量密度按指數(shù)規(guī)律衰減，無(wú)法維持一個(gè)恒定中子通量密度分布.（3）k1等于1，這時(shí)只有對(duì)應(yīng)n=1的一項(xiàng)不隨時(shí)間變化，其余隨時(shí)間衰減.（2）當(dāng)k1>1,這時(shí)所有kn-1中至少有一項(xiàng)大于1，通量密度按指數(shù)規(guī)律增加，反應(yīng)堆也無(wú)法維持一個(gè)恒定中子通量密度分布.上面三種情況分別對(duì)應(yīng)次臨界、超臨界和臨界，如圖所示。無(wú)限平板臨界理論對(duì)于臨界的反應(yīng)堆，隨著時(shí)間變大，除去第一個(gè)模態(tài)(n=1)外的所有模態(tài)(n>1)都衰減了。漸進(jìn)的（或持久的）中子通量密度取決于n=1的模態(tài)。于是有：臨界系統(tǒng)中的中子通量密度分布為：幾何曲率：系數(shù)A由功率條件決定。無(wú)限長(zhǎng)平板堆單位面積所對(duì)應(yīng)的體積所發(fā)出的功率為：重要結(jié)果兩個(gè)重要結(jié)果：裸堆單群近似的臨界方程：當(dāng)反應(yīng)堆處于臨界狀態(tài)時(shí)，中子通量密度按最小特征值所對(duì)應(yīng)的基波特征函數(shù)分布，即穩(wěn)態(tài)反應(yīng)堆的中子通量密度空間分布滿足波動(dòng)方程：重要結(jié)果分析由臨界方程：滿足臨界方程。顯然，給定系統(tǒng)的材料組成，即給定：另一方面，若尺寸給定，必須調(diào)整堆芯燃料成分，使可以得到臨界尺寸，對(duì)無(wú)限平板堆，臨界方程為：則只有唯一的尺寸a0滿足臨界方程。重要結(jié)果分析來(lái)討論臨界方程中各項(xiàng)的物理意義：是熱中子在擴(kuò)散過(guò)程中的不泄漏幾率。不泄漏幾率PL定義為：中子泄漏率=根據(jù)于是：和前一章的臨界條件完全一樣，即k1就是前面定義的的反應(yīng)堆有效增殖系數(shù)k。

可發(fā)現(xiàn)：可以看出：從不泄漏幾率：反應(yīng)堆的中子泄漏與幾何曲率有關(guān)。從平板狀反應(yīng)堆的例子中可以看到，當(dāng)反應(yīng)堆體積增大時(shí)，就減小，因而正如所預(yù)期的那樣，不泄漏幾率也就增大。同樣，擴(kuò)散長(zhǎng)度L愈大，意味著中子自產(chǎn)生到被吸收所穿行的距離也愈大，因而從反應(yīng)堆中泄漏出去的幾率也就增大，不泄漏幾率PL就要變小。重要結(jié)果分析以及設(shè)有如圖所示一維石墨慢化反應(yīng)堆，k∞=1.06，L2=300cm2，試求：（1）臨界時(shí)反應(yīng)堆的厚度H和中子通量密度的分布(設(shè)外推距離為2cm)；（2）設(shè)取H=250cm，試求反應(yīng)堆的有效增殖系數(shù)k。解：（1）根據(jù)臨界條件:可得臨界反應(yīng)堆的幾何曲率為：所以:Bg=0.01414cm-1,根據(jù)Bg=π/a,有:由外推距離d=2cm,可求得臨界時(shí)反應(yīng)堆的厚度為:例題:H=a-2d=222.2-4=218.2cm根據(jù)無(wú)限平板型均勻裸堆單群擴(kuò)散方程的解:得臨界時(shí)中子通量密度分布為:（2）若H=250cm，則反應(yīng)堆的幾何曲率：反應(yīng)堆的不泄漏幾率PL和有效增殖系數(shù)k分別等于：例題:第一類(lèi)問(wèn)題：給定反應(yīng)堆材料成分(k∞、L2

給定)，確定臨界尺寸。與臨界方程等價(jià)的臨界條件根據(jù)臨界方程:等式左端的幾何曲率:只與反應(yīng)堆的形狀和大小有關(guān).等式右端k∞和L只與反應(yīng)堆芯部材料特性有關(guān).把右邊記作:稱為反應(yīng)堆的材料曲率.反應(yīng)堆臨界條件可表述為：材料曲率等于幾何曲率，即:例題總結(jié):上述例子演示了臨界問(wèn)題中要解決的兩類(lèi)問(wèn)題。第二類(lèi)問(wèn)題：給定反應(yīng)堆形狀尺寸，確定材料成分。由于型狀尺寸給定，為已知條件，例題總結(jié):實(shí)際計(jì)算中，在反應(yīng)堆材料成分和幾何尺寸均給定情況下，求有效增殖系數(shù)k和反應(yīng)性ρρ=0，反應(yīng)堆處于臨界ρ>0，反應(yīng)堆處于超臨界ρ<0，反應(yīng)堆處于次臨界這樣，材料成分（一般是燃料的濃縮度）便可確定。

若反應(yīng)堆的有效增殖系數(shù)k=0.90，計(jì)算反應(yīng)堆的反應(yīng)性。解：核電廠常用反應(yīng)性單位：PCM1PCM=10-5（△k/k）即:1＄=100￠（1元等于100分）1＄=700PCM討論反應(yīng)堆動(dòng)態(tài)問(wèn)題時(shí),反應(yīng)性常用“元”為單位：＄1元反應(yīng)性為1個(gè)βeff(有效緩發(fā)中子產(chǎn)額,若為0.007△k/k)例題總結(jié):有限高圓柱形均勻裸堆設(shè)一有限高圓柱形均勻裸堆,高為H,半徑為R,采用圓柱形坐標(biāo),坐標(biāo)原點(diǎn)位于軸線的半高度上.拉普拉氏算符的表達(dá)式為:中子通量密度分布是對(duì)稱的,與θ無(wú)關(guān),有限高圓柱均勻裸堆的波動(dòng)方程可以寫(xiě)為:邊界條件：(1)不計(jì)外推長(zhǎng)度時(shí),反應(yīng)堆外邊界上中子通量密度為零;(2)中子通量密度分布對(duì)稱.代入方程有限高圓柱形均勻裸堆采用分離變量法，令：上式等號(hào)左邊只與r有關(guān)，等號(hào)右邊只與z有關(guān)。因?yàn)閞和z是兩個(gè)彼此無(wú)關(guān)的獨(dú)立變量，要使改式保持相等關(guān)系，只有兩邊都是常數(shù)才行。該常數(shù)記為a2，則得到:得到：來(lái)看方程：有限高圓柱形均勻裸堆由邊界條件:方程的解為：其中，C1，C2是兩個(gè)待定常數(shù)。J0(x)，Y0(x)分別是零階第一類(lèi)及第二類(lèi)貝塞爾函數(shù)，它們隨x的變化見(jiàn)圖所示。有限高圓柱形均勻裸堆再來(lái)看方程：令：化為零階貝塞爾方程：該方程的通解為：零階貝塞爾函數(shù)圖有限高圓柱形均勻裸堆當(dāng)x＝0時(shí)，J0(0)＝1，而Y0(0)→-∞，而r＝0處堆內(nèi)中子通量密度有限，所以C2＝0。則上式可寫(xiě)成：有限高圓柱形均勻裸堆由邊界條件：正如前圖所示，在x＝x1，x2…等處，函數(shù)J0(x)為零，即J0(xn)＝0。故必須取如下的本征值：討論方程的解：與此同時(shí)，當(dāng)反應(yīng)堆達(dá)到臨界時(shí)，只有n＝1這個(gè)本征值才有物理意義。即：因?yàn)椋核裕号R界時(shí)，r向的中子通量密度分布為：綜上分析，有限高圓柱形均勻裸堆的中子通量密度分布公式為：有限高圓柱形均勻裸堆公式中常數(shù)A的確定：取環(huán)狀體積元

，在(r,z)處該體積元所發(fā)出來(lái)的功率為：從0到R對(duì)r積分，從-H/2到H/2對(duì)z積分可得到功率：代入中子通量密度的表達(dá)式，上式可化為：有限高圓柱形均勻裸堆裸堆由直徑為2.54cm的天然鈾金屬長(zhǎng)棒排列成柵距為0.152m的正方形柵格并懸掛在圓柱形容器內(nèi)所組成。容器內(nèi)盛有重水作為慢化劑。堆芯的高度、直徑之比為1.2。由材料的性質(zhì)已知＝8.6m-2。試估算使這一反應(yīng)堆剛好臨界的天然鈾質(zhì)量。例題：解：有限高圓柱堆的幾何曲率為：因?yàn)镠/(2R)＝1.2,即H＝2.4R代入：對(duì)于臨界反應(yīng)堆：所以：R2=0.8721m2，R=0.93m，H=2.4R=2.4x0.93=2.23m

燃料棒位于柵距為0.152m的正方形柵格節(jié)點(diǎn)上。因此每根棒分得的有效面積為0.0231m2。堆芯的橫截面積為：

R2＝

(0.93)2＝2.72m2。所以燃料棒數(shù)為：2.72/0.0231≈117.6。每根棒的直徑為2.54cm，其長(zhǎng)度和堆芯相同，即2.23m。燃料棒體積V為：例題(續(xù))：金屬鈾的密度＝19.0x103kg/m3，所以燃料總質(zhì)量為：m＝V＝0.132×19×103＝2.51×103kgm3一般說(shuō)來(lái)，對(duì)于無(wú)限長(zhǎng)有限厚的平板堆，臨界時(shí)，其厚度應(yīng)滿足：或?qū)τ谟邢薷邎A柱堆，臨界時(shí)，其高度H與半徑R須滿足：因?yàn)榕R界時(shí)：該式只給出R與H的關(guān)系而不能給出確定的值。但如果再加上一個(gè)條件，即要求最小臨界體積，便能從式中解得所謂最佳半徑或最佳高度。反應(yīng)堆的體積為V＝

R2H，臨界尺寸臨界尺寸將R2值代入體積表達(dá)式并求V的極小值，令：具有最小臨界體積的圓柱均勻堆，要求直徑與高大致有0.5的關(guān)系。

實(shí)際動(dòng)力堆堆芯尺寸投入運(yùn)行年電站或核船名稱輸出功率/MWe堆芯高×直徑/m1972施塔德6302.99×3.051972緬茵·楊基7933.6×3.51973哈欽森島18253.5×3.51974勇士11303.66×3.31976比布利斯11803.9×3.61962薩瓦娜號(hào)核商船1.7×1.61968奧托·哈恩號(hào)核商船1.12×1.151991秦山一期核電廠3002.90×2.4861995大亞灣核電廠9003.65×3.362002秦山三期7005.945×6.2862006田灣10003.53×3.16中子通量密度分布不均勻系數(shù)中子通量密度分布不均勻系數(shù)定義為堆芯最大熱中子通量密度與堆芯平均熱中子通量密度的比值，即：最大中子通量密度平均中子通量密度對(duì)于圓柱堆r=0處，其中子通量密度最大，即:圓柱體積:平均中子通量密度:不均勻系數(shù):有限高圓柱堆的中子通量密度分布不均勻系數(shù):中子通量密度分布不均勻系數(shù)中子通量密度分布徑向不均勻系數(shù)：中子通量密度分布軸向不均勻系數(shù)：不同堆芯幾何形狀的臨界均勻裸堆的曲率、中子通量密度分布與最小體積比較表均勻堆的中子通量密度表堆芯幾何形狀尺寸幾何曲率中子通量密度分布系數(shù)A最小體積無(wú)限平板厚1.57-長(zhǎng)方體3.88無(wú)限高圓柱半徑R2.32-有限高圓柱半徑R高H3.64球體半徑R3.29雙群擴(kuò)散理論

一群擴(kuò)散理論簡(jiǎn)單，只能給出近似結(jié)果，采用一群擴(kuò)散理論來(lái)分析將會(huì)帶來(lái)較大的偏差。能群數(shù)越多，計(jì)算結(jié)果越準(zhǔn)確，然而計(jì)算量是相當(dāng)大的。需要結(jié)合堆型綜合權(quán)衡。對(duì)于熱中子堆，利用雙群可得到較好的結(jié)果。堆型重水堆壓水堆高溫氣冷堆快中子反應(yīng)堆能群數(shù)1-24

46-18不同類(lèi)型反應(yīng)堆采用的離散化能群的數(shù)目能群序號(hào)中子能量范圍能區(qū)名稱中子譜1(10-0.821)MeV高能區(qū)裂變譜20.821MeV-5.53keV中能區(qū)1/E譜35.53keV-0.625eV共振能區(qū)1/E譜＋修正40.625eV-0eV熱能區(qū)麥克斯韋爾譜PWR少群能量劃分堆內(nèi)中子按能量劃分成兩群：熱中子歸為一群，簡(jiǎn)稱熱群；高于熱能的中子歸為一群，簡(jiǎn)稱快群；兩群分界能對(duì)于水堆約為0.6至1電子伏。在熱中子堆內(nèi)，快群中子主要是熱中子所引起的裂變產(chǎn)生，它又通過(guò)慢化和泄漏而消失；而熱中子則來(lái)源于快群中子的慢化，由于吸收和泄漏而消失。因此，根據(jù)中子平衡關(guān)系可列出反應(yīng)堆穩(wěn)態(tài)時(shí)芯部的快群及熱群的擴(kuò)散方程為：雙群擴(kuò)散理論下標(biāo)1,2代表快群和熱群?！苧為快群的移出截面，∑1→2為快群到熱群的散射截面。擴(kuò)散系數(shù)及所有的中子截面都是經(jīng)中子能譜平均后的平均值。

求解該方程，可得反應(yīng)堆內(nèi)雙群中子通量密度的典型分布曲線，如圖：

突起的原因：反射層的熱中子吸收較小，慢化能力較強(qiáng)。雙群擴(kuò)散理論用雙群擴(kuò)散理論，有效增殖系數(shù)k為：其中，Ps為快中子在慢化過(guò)程中的不泄漏幾率，Pd為熱中子擴(kuò)散過(guò)程的不泄漏幾率。兩者的定義為：雙群理論的臨界方程為：對(duì)于大型反應(yīng)堆，很小，若略去項(xiàng)，則得到：式中，M2=τ+L2，這樣便得到了修正單群理論的臨界方程。來(lái)看雙群理論的臨界方程為：雙群擴(kuò)散理論快群參數(shù)慢化劑DF/m

FT/m-1

/10-4m2H2O0.01134.1927D2O0.01290.985131Be0.005620.551102石墨0.010160.276368雙群擴(kuò)散理論熱群參數(shù)慢化劑密度/103kg/m3DT/m

aT/m-1/10-4m2H2O1.00.00161.978.1D2O1.100.00872.9x10-33.0x104Be1.850.0050.104480石墨1.600.00842.4x10-23500雙群擴(kuò)散理論例題設(shè)一輕水冷卻圓柱形反應(yīng)堆堆芯，其核參數(shù)為：求：（1）設(shè)堆芯高度H=3.55m，試求堆芯的臨界半徑。（2）如果給定堆芯半徑R=1.56m，求反應(yīng)堆的有效增殖因數(shù)為多少？（1）求出反應(yīng)堆堆芯的外推距離：m根據(jù)修正單群理論，材料曲率為：幾何曲率為：由例題解答m（2）反應(yīng)堆的有效增殖因數(shù)R=1.56m，幾何曲率為有效增殖因數(shù)為：例題解答3.2有反射層的均勻堆反射層的定義：包圍在反應(yīng)堆芯部外面用以反射從芯部泄漏出來(lái)的中子的物質(zhì)稱作反射層。反射層的作用：（1）減少芯部中子的泄漏，使芯部臨界尺寸比無(wú)反射層小，節(jié)省然料；（2）提高反應(yīng)堆的平均輸出功率；

反射層材料的要求：（1）散射截面∑s大；（2）吸收截面∑a要?。唬?）慢化能力要強(qiáng)。常用的反射層材料有：水、重水、石墨和鈹?shù)?。反射層的作用反射層?jié)省在芯部包有反射層以后，芯部臨界大小的減少量稱為反射層節(jié)省，用δ表示。對(duì)球形反應(yīng)堆：δ=R0-R對(duì)圓柱形反應(yīng)堆：徑向反射層節(jié)?。狠S向反射層節(jié)?。河玫刃愣训膸缀吻时硎荆河绊懛瓷鋵庸?jié)省的因素？以球形堆為例。球形反應(yīng)堆的反射層節(jié)省δ可以表示為：T為反射層厚度，Lr為反射層中的中子擴(kuò)散長(zhǎng)度。分兩種情況討論:（1）反射層厚度很小，即T<<Lr，反射層節(jié)省可以寫(xiě)成：δ≈T（2）反射層很厚，即T>>Lr，這時(shí)于是：δ≈Lr

即反射層厚度增加到一定值后，反射層節(jié)省δ就達(dá)到一個(gè)常數(shù)值。反射層節(jié)省反射層對(duì)中子通量密度分布的影響反射層節(jié)省反射層對(duì)徑向和軸向中子通量密度分布的影響：反射層節(jié)省反照率描述不同物理性質(zhì)介質(zhì)在交界面上粒子的入射和反射的情況。是描述出射該介質(zhì)的中子占入射該介質(zhì)的中子的份額。是一個(gè)中子通過(guò)介質(zhì)表面入射并從此表面返回的幾率。主要用來(lái)描述堆芯和反射層交界面處的邊界條件。根據(jù)斐克定律，反照率可以表示成：在進(jìn)行堆芯物理計(jì)算時(shí)，在堆芯和反射層