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

非均匀堆的均匀化计算步骤 栅元的均匀化计算。计算由燃料、包壳和慢化剂组成的一维圆柱栅元，求出多群中子通量密度空间—能群分布，归并求出少群或宽群的栅元均匀化截面 燃料组件的均匀化计算。由栅元均匀化计算结果，进行燃料组件的二维计算，得到组件均匀化少群截面常数 全堆芯的扩散计算。利用燃料组件少群均匀化常数进行全堆扩散计算，求出堆芯有效增值系数和中子通量密度及功率分布。 * 6.3 栅元均匀化群常数的计算 均匀化群常数计算的考虑 计算的主要问题是求介质内中子通量密度分布 栅元介质具有强吸收性 = 扩散理论不适用! 采用更精确的输运理论：SN方法、CPM方法、MC方法等 采用维格纳-赛兹（Wigner-Seitz）等效栅元近似 栅元均匀化群常数： 正方形栅元 六角形栅元 * 6.3.1 积分输运理论（碰撞概率法）的基本方程 假设在L系内中子与原子核的散射各向同性，r’处源Q(r’,E)所产生中子对r处的中子通量密度的贡献为 为连接r’与r点的直线路径的“光学距离”， 也就是以平均自由程为单位量度的距离。当?t为常数时，等于 其中 * 6.3.1 积分输运理论的基本方程 对于栅元计算，边界条件：等效栅元的边界为各向同性全反射且净中子流等于零，或者是入射流为零的真空边界。空间任意点r处的中子通量密度应等于单位时间内从所有其它点r’产生的源中子对r处中子通量密度贡献之和 积分中子输运方程 积分输运方程 vs. 扩散方程 * 输运方程的数值求解—碰撞概率法（CPM：Collision Probability Method） 1. 首先将系统划分为I个互不相交的均匀子区 当区域划分足够小时，可假设： 每一子区的截面参数为常数或可用该区的平均值表示 每一子区内的中子源强或中子通量密度等于常数（平源或平通量密度近似） 2. 对能量变量采用分群近似求解，把能量区间自E0到0分为E0E1…Eg…EG=0，G个能量区间，称为G群近似。 * 碰撞概率法（CPM：Collision Probability Method） 3. 在积分输运方程两端乘以?t，然后在每一子区体积内Vi及能量区间?Eg内对方程进行体积与能量积分，并按照分群近似方法处理，得 其中 分别表示第g群第i区平均中子通量密度和第g群第j区平均中子源强 * 源项： 忽略外源，得碰撞概率形式积分输运方程 群常数及首次碰撞概率可事先求得，上式为一线性方程组，可用迭代方法求解。CPM方法关键是首次碰撞概率计算，与几何及材料有关，专门程序计算。 Pij,g为第j区内产生的一个各向同性中子不经任何碰撞到达i区发生首次碰撞概率 * 6.3.2 碰撞概率方程的解及少群常数计算 迭代方程为： 其中： 用源迭代方法求解碰撞概率形式积分输运方程： 收敛准则： * 栅元的均匀化截面： 少群的群转移均匀化截面： 均匀化栅元的少群截面（n表示多群编号，g表示少群编号）： 6.4 燃料组件内均匀化通量密度分布及少群常数的计算 燃料组件计算的考虑： 燃料组件的计算可等效为二维问题，使问题简化 组件数很多，采用碰撞概率法计算量大，耗时 界面流方法(穿透概率法)：把组件分成若干个子区，各个子区之间通常用界面流必须连续的条件来耦合。能够使计算大大简化。 根据中子平衡原理，(i,j)子区Sn表面的出射流方程： 分别为(i,j)子区Sn表面上的出射和入射总中子流 为(i,j)子区平均中子通量密度 为入射到Sm表面上的g能群中子未经碰撞穿过子区从Sn表面逸出的首次穿透概率 根据中子守恒关系从输运方程可以推得(i,j)子区的中子平衡方程 方程右边第一项为子区内对Sn表面出射中子流的贡献，第二项表示入射到子区表面Sm上的中子未经碰撞直接穿出Sn的数目，两部分之和构成Sn表面总出射中子流 * 6.5 共振区群常数的计算 共振能区（几ev~0.01Mev）的考虑 共振能区，对某些核素存在强烈共振吸收现象 共振吸收截面与能量、栅元的几何结构、介质温度等有关 在多群库中，并不直接给出共振截面，而是给出一些共振参数数据，根据栅元的具体结构计算得出。 g群共振吸收截面 定义第i个共振峰有效共振积分 第g个能群有效共振积分 g群共振吸收截面转换为 * 6.5.1 非均匀栅元有效共振积分的计算 有效共振积分计算的考虑 非均匀栅格的空间-能量分布不均匀性，计算复杂 孤立棒栅元假设：假定燃料块间的距离大于中子在慢化剂内的平均自由程，如重水堆、石墨堆 只取一个栅元考虑：在孤立棒假设下，从一个燃料块飞出的共振中子不可能在穿过慢化剂时未经碰撞而仍以其原来能量进入相邻的另一个燃料块 假设栅元只是由燃料和慢化剂组成，且燃料由由一种元素组成，计算结果可扩展到多种元素 * 6.5.1 非均匀栅元有效共振积分的计算 在燃料快内能量高于E的中子