第6 章 栅格的非均匀效应与均匀化群常数的计算(反应堆物理分析)

1、目 录p6.1 栅格的非均匀效应p6.2 栅格的均匀化处理p6.3 栅元均匀化群常数的计算p6.4 燃料组件内均匀化通量密度分布及少 群常数的计算p6.5 共振区群常数的计算p6.6 栅格几何参数的选择6.1 栅格的非均匀效应p堆的分类(按照堆芯内燃料和慢化剂的分布形式)均匀堆：均匀堆：燃料和馒化剂均匀混合在一起，例如把铀和慢化剂制成铀盐溶液非均匀堆：非均匀堆：在非均匀堆中，把燃料制成块状，如圆柱状(棒状)、环状、球状、片状(平板状)等，按一定的几何形式放置在慢化剂中构成所谓栅格结构的堆芯。目前的动力反应堆几乎都是非均匀的。非均匀反应堆非均匀反应堆世界上第一座反应堆是世界上第一座反应堆是非均匀

2、非均匀反应堆反应堆? ?当时还没有铀浓缩技术铀浓缩技术，可以得到的核燃料只有天然铀。天然铀与任何慢化剂（即使用吸收中子极少的石墨或重水）组成的均匀混合物都不可能达到临界均匀混合物都不可能达到临界。费米和西拉德等人敏锐地看出，把燃料做成块状块状，与慢化剂分开，可以得到更好的效果。把天然铀块插入石墨块，把这样的石墨块一块块堆起来，堆到一定大小，就达到了临界。他们把此临界装置称为Nuclear Pile天然铀与慢化剂的均匀混合物为什么达不到临天然铀与慢化剂的均匀混合物为什么达不到临界？界？p欲达临界，欲达临界，k k 一定要大于一定要大于 1 1。p天然铀石墨的均匀混合物天然铀石墨的均匀混合物, ,

3、各因子可达到的最大值是各因子可达到的最大值是1.050.591.330.855pfk有关概念p 反应堆栅格反应堆栅格(Lattice):(Lattice):在非均匀反应堆中，按照某种有规则的图形布置的燃料和其它材料的阵列p 栅元栅元(Cell)(Cell)：反应堆各栅格中具有相同材料组成和几何形状的单元p 燃料棒尺寸与堆芯内中子平均自由程相比： 轻水堆：约大于1cm,与棒径相当，必须细致考虑非均匀效应； 快堆：约为数十cm，可粗略处理非均匀效应VENUS-1# ADS sub-critical reactorSCWR-thermal spectrum reactorAP1000-PWR1.1.

4、 热中子分布：热中子分布：p 主要在慢化剂中产生：燃料慢化能力比慢化剂小p 主要被燃料核吸收：燃料吸收热中子能力比慢化剂大p 形成从慢化剂流向燃料快的热中子流p 空间自屏效应空间自屏效应：外层燃料核对内层燃料核的屏蔽作用p 缺点缺点：燃料吸收热中子能力下降，热中子吸收系数 减小非均匀栅格内中子的空间分布非均匀栅格内中子的空间分布1-热中子热中子; 2-共振中子共振中子; 3-裂变中子裂变中子2. 共振中子分布：共振中子分布：p共振中子主要在慢化剂中产生p燃料具有很大共振吸收截面造成强烈空间自屏效应空间自屏效应p共振中子与燃料核碰撞概率减小，与慢化剂碰撞概率增大p优点：优点：燃料核共振吸收能力减

5、小, 逃脱共振俘获概率增大非均匀栅格内中子的空间分布非均匀栅格内中子的空间分布1-热中子热中子; 2-共振中子共振中子; 3-裂变中子裂变中子3.3.裂变中子分布：裂变中子分布：p在燃料块内产生p约60%具有1.1MeV以上的能量，可引起238U裂变p优点：优点：燃料制成块状，增加了快中子增殖效应, 增大非均匀栅格内中子的空间分布非均匀栅格内中子的空间分布1-热中子热中子; 2-共振中子共振中子; 3-裂变中子裂变中子栅格的非均匀效应p非均匀栅格内的中子通量密度分布是不均匀的；p空间自屏效应空间自屏效应 ：p变大，f 变小（小于p的效应）， 略有上升；p合理选择燃料快的直径或厚度、栅距等，可以

6、提高无限介质增殖因子K（四因子公式）6.2 栅格的均匀化处理p非均匀反应堆有成千上万根燃料棒 , 严格计算非常复杂p非均匀反应堆中子通量密度分布：宏观 + 微观p实际计算中，非均匀反应堆等效为均匀反应堆非均匀反应堆等效为均匀反应堆均匀化处理均匀化处理p均匀化思想：均匀化思想：用一个等效的均匀介质来代替非均匀栅格，使得计算结果（特征物理量，如中子反应率）计算结果（特征物理量，如中子反应率）与非均匀栅格相等或近似。p核心问题核心问题：如何确定等效均匀化介质的各种中子截面参数或有效群常数p首选物理量首选物理量：栅元内各能群的各种中子反应率中子反应率GgsfaxdVdEErErdVdEErggEVxE

7、Vgx, 1;,),(),(),(, 因而：,( ,)( ,)( ,)ggEgVx ggEVx r Er E dVdEr E dVdEffDDSS=蝌蝌,( ,)( ,) ( ,), ,;1,gx gEVEgVr E dVdEx r Er E dVdExa f sgGffDDS=S=蝌蝌( ,)( ,)ggEVEVr E dVdEr E dVdEffDD=蝌蝌近似认为：非均匀反应堆的计算p把栅格均匀化把栅格均匀化，考虑非均匀效应计算出等效均匀化系统的均匀化群常数；p把非均匀系统等效成均匀系统把非均匀系统等效成均匀系统。利用上面计算的均匀化常数，采用均匀反应堆理论计算临界大小、中子通量密度分布、

8、功率分布等这种处理方法叫做非均匀反应堆的均匀化处理非均匀反应堆的均匀化处理堆芯的均匀化截面的计算非均匀堆的均匀化计算步骤 栅元的均匀化计算。计算由燃料、包壳和慢化剂组成的一维圆柱栅元，求出多群中子通量密度空间能群分布，归并求出少群或宽群的栅元均匀化截面燃料组件的均匀化计算。由栅元均匀化计算结果，进行燃料组件的二维计算，得到组件均匀化少群截面常数全堆芯的扩散计算。利用燃料组件少群均匀化常数进行全堆扩散计算，求出堆芯有效增值系数和中子通量密度及功率分布。216.3 栅元均匀化群常数的计算p 均匀化群常数计算的考虑 计算的主要问题是求介质内中子通量密度分布主要问题是求介质内中子通量密度分布 栅元介质

9、具有强吸收性 = 扩散理论不适用扩散理论不适用! ! 采用更精确的输运理论：SN方法、CPM方法、MC方法等 采用维格纳采用维格纳-赛兹（赛兹（Wigner-Seitz）等效栅元近似）等效栅元近似栅元均匀化群常数：R/0.56419R0.52504aaa正方形栅元六角形栅元,),(),(),(,dVdEErdVdEErErggEVEVxgx 226.3.1 6.3.1 积分输运理论（碰撞概率法）的基本方程积分输运理论（碰撞概率法）的基本方程trr/ 假设在L系内中子与原子核的散射各向同性，r处源Q(r,E)所产生中子对r处的中子通量密度的贡献为),(rrE为连接r与r点的直线路径的“光学距离光

10、学距离”，也就是以平均自由程为单位量度的距离以平均自由程为单位量度的距离。当t为常数时，等于2exp ( ,)( ,)4E rrQ r Err0),(),(rrtdllErrE其中236.3.1 6.3.1 积分输运理论的基本方程积分输运理论的基本方程24),(exp),(),(dVrrrrEErQErV24),(exp),(rrrrEErQ 对于栅元计算，边界条件：等效栅元的边界为各向同性全反射且净中子流等于零，或者是入射流为零的真空边界。空间任意点r处的中子通量密度应等于单位时间内从所有其它点r产生的源中子对r处中子通量密度贡献之和积分中子输运方程积分中子输运方程积分输运方程积分输运方程

11、vs. 扩散方程扩散方程24输运方程的数值求解输运方程的数值求解碰撞概率法碰撞概率法（CPM：Collision Probability Method）jiVVVVIiVjiii,0,;, 1,12iI1. 首先将系统划分为I个互不相交的均匀子区均匀子区当区域划分足够小时，可假设： 每一子区的截面参数为常数或可用该区每一子区的截面参数为常数或可用该区的的平均值表示平均值表示 每一子区内的中子源强或中子通量密度每一子区内的中子源强或中子通量密度等于常数等于常数（平源或平通量密度近似）（平源或平通量密度近似）2. 对能量变量采用分群近似求解分群近似求解，把能量区间自E0到0分为E0E1EgEG=0

12、，G个能量区间，称为G群近似。25碰撞概率法碰撞概率法（CPM：Collision Probability Method）IjjgijjgiigigtVPQV1, gjgiEVjjgEViigdVdEErQVQdVdEErV),(1),(1,12iI3. 在积分输运方程两端乘以t，然后在每一子区体积内Vi及能量区间Eg内对方程进行体积与能量积分，并按照分群近似方法处理，得其中24),(exp),(),(dVrrrrEErQErV分别表示第g群第i区平均中子通量密度和第g群第j区平均中子源强26 IjjgijGgjggfgjggiigigtVPkV1,1, ,GggjggfGggjgjggjgS

13、kQ1, ,1, ,源项源项：忽略外源，得碰撞概率形式积分输运方程碰撞概率形式积分输运方程群常数及首次碰撞概率可事先求得，上式为一线性方程组，可用迭代方法求迭代方法求解。解。CPMCPM方法方法关键是首次碰撞概率计算，与几何及材料有关关键是首次碰撞概率计算，与几何及材料有关，专门程序计算。IjjgijjgiigigtVPQV1, ,20exp();()4ijrrgt g iij gjigt gVVjrrPdV dVrrdlVrr Pij,g为第j区内产生的一个各向同性中子不经任何碰撞到达i区发生首次碰撞概率首次碰撞概率276.3.2 6.3.2 碰撞概率方程的解及少群常数计算碰撞概率方程的解及

14、少群常数计算IjjnjfnjjnjfnVQkVQk1)2(,)2(I1)1(,)1(1GgnjgfnjgjfQ1)1(, ,)1(,IjjgijnjfIjGgnggijjnjgjgginigigtVPQkPVV1,)1(,11)1(,)(, ,)(,迭代方程为迭代方程为：其中：用源迭代方法源迭代方法求解碰撞概率形式积分输运方程：收敛准则收敛准则：2,1)()1()()()1()(max;nnnigigigginnnkkk IjjgijGgjggfgjggiigigtVPkV1,1, ,28GgtfaxVVIiiigIiiigigxgx, 2 , 1;,栅元的均匀化截面栅元的均匀化截面：igni

15、ingngniiinnnggVV, , GgtfaxVVgnIiiingnIiiininxgx, 2 , 1;,少群的群转移均匀化截面少群的群转移均匀化截面：均匀化栅元的少群截面均匀化栅元的少群截面（n n表示多群编号，表示多群编号，g g表示少群编号）表示少群编号）：6.4 燃料组件内均匀化通量密度分布及少群常数的计算p燃料组件计算的考虑：燃料组件计算的考虑：燃料组件的计算可等效为二维问题，使问题简化组件数很多，采用碰撞概率法计算量大，耗时界面流方法(穿透概率法)：把组件分成若干个子区，各个子区之间通常用界面流必须连续的条件来耦合。能够使计算大大简化。J1J -1J =21i =2 31i-

16、i( , )i jmS3sJ-3sJ+1sJ-1sJ+2sJ-2sJ+4sJ-4sJ+2s3s4s1s, i j根据中子平衡原理，(i,j)子区Sn表面的出射流方程：,41( , )( , )( , )( , )( , )1, ;1, ;1,;1,4n gnm gnmggSijgS VSS SmJi jV Q i j Pi jJi j Pi jiIjJgGn+-=+=1( , )( , )() ( , )( , )( , )GggfggggggQ i ji ji ji jSi jkcnf =S+S+,n gSJ分别为(i,j)子区Sn表面上的出射和入射总中子流( , )gi jf为(i,j)子

17、区平均中子通量密度( , )ngS VPi j为入射到Sm表面上的g能群中子未经碰撞穿过子区从Sn表面逸出的首次穿透概率3sJ-3sJ+1sJ-1sJ+2sJ-2sJ+4sJ-4sJ+2s3s4s1s, i j,41( , )( , )( , )( , )( , )1, ;1, ;1,;1,4n gnm gnmggSijgS VSS SmJi jV Q i j Pi jJi j Pi jiIjJgGn+-=+=1( , )( , )() ( , )( , )( , )GggfggggggQ i ji ji ji jSi jkcnf =S+S+,1()()( , )( , )lu gru glu

18、 gru gt gggux yi jJJJJi jQ i jVf+-=+-+ S=根据中子守恒关系从输运方程可以推得(i,j)子区的中子平衡方程 方程右边第一项为子区内对Sn表面出射中子流的贡献，第二项表示入射到子区表面Sm上的中子未经碰撞直接穿出Sn的数目，两部分之和构成Sn表面总出射中子流326.5 6.5 共振区群常数的计算共振区群常数的计算gEFgiigadEEI)(,FVFFdVErVE),(1)(iEFaidEEEI)()(ggEFEFagadEEdEEE)()()(,p 共振能区（几共振能区（几ev0.01Mevev0.01Mev）的考虑）的考虑 共振能区，对某些核素存在强烈共振

19、吸收现象 共振吸收截面与能量、栅元的几何结构、介质温度等有关 在多群库中，并不直接给出共振截面，而是给出一些共振参数数据，根据栅元的具体结构计算得出。pg g群共振吸收截面群共振吸收截面定义第i个共振峰有效共振积分第g个能群有效共振积分giiEFagIdEEEIg)()(g g群共振吸收截面转换为群共振吸收截面转换为336.5.1 6.5.1 非均匀栅元有效共振积分的计算非均匀栅元有效共振积分的计算p 有效共振积分计算的考虑有效共振积分计算的考虑非均匀栅格的空间-能量分布不均匀性，计算复杂孤立棒栅元假设孤立棒栅元假设：假定燃料块间的距离大于中子在慢化剂内的平均自由程，如重水堆、石墨堆只取一个栅

20、元考虑：在孤立棒假设下，从一个燃料块飞出的共振中子不可能在穿过慢化剂时未经碰撞而仍以其原来能量进入相邻的另一个燃料块假设栅元只是由燃料和慢化剂组成，且燃料由由一种元素组成，计算结果可扩展到多种元素346.5.1 6.5.1 非均匀栅元有效共振积分的计算非均匀栅元有效共振积分的计算/,)1 ()()(dEdEEEEVFEEFFFsF在燃料快内能量高于E的中子与燃料核弹性碰撞后进入E与E+dE能量范围内的中子数这些中子在燃料快内发生首次碰撞的数目燃料快内发生首次碰撞的数目为/,0)1 ()()()(1dEdEEEEEPVFEEFFFsFFP PF0F0(E):(E): 在燃料块内产生的均匀和各向同

21、性分布、能量为E的中子未经碰撞逸出块外在慢化剂内发生首次碰撞的概率，即首次飞行逃脱概率首次飞行逃脱概率。F和F为燃料块内和慢化剂内共振中子平均通量密度35EEEPdEEEEEPEEFtFEEFFFsFFFtF)()()1 ()()()(1)()(,0/,0,EE1)(0)()(0,0,EPVEPVMMMtFFFt根据中子平衡原理，得燃料快内的中子慢化方程燃料快内的中子慢化方程假设共振峰很窄，慢化剂内能谱分布慢化剂内能谱分布为：互易关系式互易关系式在慢化剂内慢化到慢化剂内慢化到E E与与E+dEE+dE能量范围内的中子在燃料快中发生首次碰撞的数目能量范围内的中子在燃料快中发生首次碰撞的数目/,0

22、)1 ()()()(dEdEEEEEPVMEEMMMsMM(6-35)假设 t,Mt,M s,Ms,M, , (6-33) (6-33)式可转换为式可转换为dEEPVFFFt0,(6-36)(6-37)P PM0M0(E):(E): 慢化剂内能量为E的中子在燃料块内发生首次碰撞的概率。含F (E)的积分方程，若PF0已知，求解方程可得F (E) 。(6-33)366.5.2 6.5.2 等价原理等价原理首次飞行逃脱概率首次飞行逃脱概率P PF0F0非常复杂，非常复杂，近似方法求解近似方法求解。SVl4定义燃料棒平均弦长燃料棒平均弦长(6-38)le1假想的假想的“逃脱逃脱”宏观截面宏观截面(6

23、-39)VNSlNNFFFee41假想的假想的“逃脱逃脱”微观截面微观截面(6-40)lPFtFteeF,011维格纳有理近似公式维格纳有理近似公式(6-41),()()(I)(0)(, 1,1)()(,eEFaeFsFaeFpEFaFeFsFaeFpFTFdEEdEEENRIMNREEEii近似宽共振，近似窄共振(6-42)对比(6-42)与(2-75): 若非均匀反应堆的e等于均匀堆的FMsMNN/,，则两者的有效共振积分相等，即为等价原理等价原理。376.5.4 6.5.4 温度对共振吸收的影响温度对共振吸收的影响p多普勒效应多普勒效应：= 非均匀堆共振吸收随温度而改变能量自屏能量自屏空间自屏空间自屏p“能量自屏能量自屏”效应效应)1(11,01101,eFpiEeFpFaFaiEeFsFaeFpFaiEEdEEEdEIii386.5.4 6.5.4 温度对共振吸收的影响温度对共振吸收的影响)1 (,01101eFpiiEEI)1 (,02202eFpiDEEI02IDIi =共振吸收随温度升高而升高共振吸收随温度升高而升高p“能量自屏能量自屏”效应效应396.5.4 6.5.4 温度对共振吸收的影响温度对共振吸收的影响=共振吸收随温度升高而升高共振吸收随温度升高而升高 燃料温度升高，由于多普勒效应，能