首次碰撞源方法减弱射线效应研究.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 首次碰撞源方法减弱射线效应 刘永康 胡永明 清华大学核能与新能源技术研究院，北京，100084 摘要：通过构造梭长是与距点源的距离为二阶小量的立方体单元，改善了在与源区相邻位置的首次碰撞源 计算效果，利用点源近似修正进一步提高了整个系统的计算精度，编制了适用于三维离散纵标计算程序NTS 的首次碰撞源计算程序FCS，通过对包含真空区域的Kobayashi三维基准问题的校核计算结果，本方法能够 有效减弱离散纵标计算射线效应． 关键词：射线效应；首次碰撞源；离散纵标法；真空区域 1引言 离散纵标方法是中子输运方程数值计算方法中最常采用的确定性方法之一，得到了广泛 的应用，但是由于在离散纵标方法中。只允许中子沿着有限的特定方向运动，会产生所谓的“射 线效应”，在散射截面很小和系统中只有少量分布源的情况下，例如在研究堆辐照孔道附近的 中子通量计算中，堆射线效应会造成很大的误差，甚至会严重歪曲计算结果，为了消除或减 弱射线效应的影响，已经提出了很多实用方法，如果要完全消除射线效应的影响，可以采用 某种方法使离散纵标方程等价于某种形式的PIl方程，即“虚源方法”，但是虚源方法存在收敛 慢的缺点Ill。在大多数应用中。将射线效应减少到某种限度之下就可以满足要求了，在这种 情况下，需要采取某种方法减轻而非完全消除射线效应的影响，例如本文采用的首次碰撞源 方法，该方法虽然不能完全消除射线效应，但在绝大多数情况下，能给出足够精确的结果。 首次碰撞源方法的原理非常简单，即通过半解析方法求出由于分布源产生的未碰撞中子注量 率，利用该注量率计算出空间每个网格的自散射反应率，这些反应率相当于在该网格的一个 分布源，这样，就将少量的孤立源分布到更广泛的空阅范围中，从而减轻了射线效应的影响。 本文利用首次碰撞源减缓射线效应方法编制了三维首次碰撞源(FCS)计算程序，与自 行开发的三维离散纵标中子输运计算程序NTS配合，可以用较短的计算时问，将Kobayashi 基准问题12J的计算误差控制在20％以内，有效的减弱了射线效应的影响。 2首次碰撞源减缓射线效应的方法 为方便起见，在本文的讨论中，假定都是单群输运问题，公式变量都不带能群下标。 在空间某点毛处存在点源吼，则系统中任意位置声处的未碰撞中子注量率为 ∥。(I；，盎)=谚1印制 16 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 式中户=l乒一无I是点源所在点乞到系统中任意点j}的距离， 口(户)=fq(乒)z是从点源所在点磊到任意点声的光学平均自由程。 如果能将(1)式对角度和空间积分，则可以得到空间任意网格由于点源go而产生的未碰撞 中子注量率，在X-Y-Z几何下，即 ’蛇旷e e ef．暖删肚纰 (2) 41r 格边界指标 蛇，，j是单元平均的未碰撞中子注量率 ∥，国是角空间独立变量 万代表DiracDelta函数 但是，除了极特殊的情况外，(2)式是无法解析积分的，数值积分要达到满意的精度要求 所需要的计算量也非常大[31，一般情况下只能通过近似方法来求解。 2．1穿透几率形式的首次碰撞源方法 在文献[41中，Alcouffe等提出了一种能满足离散纵标输运平衡的首次碰撞源方法，其基 本原理是：(1)首先假定在点源和目标单元(可以理解成探测器)之间是真空，即中子流没 有衰减，这样(2)式中的指数衰减项可以去掉：(2)在(1)求出的真空中子流上乘上从点源到探 测器的指数衰减即得到探测器位置的实际中子流：(3)利用前面求出的中子流，通过中子守 恒方程，计算出探测器位置的未碰撞中子注量率．在文献14J中，Alcouffe等通过射线追踪方 法求出从点源到网格顶点的指数衰减，并用网格每个顶点指数衰减的算术平均值作为该面的 平均指数衰减。这种方法是二阶精度的，能够较好的解决那些在一个较大范围内存在少量孤 立源，而计算只需要考虑离源较远的位置的问题，对于需要考虑离源较近位置的问题，该