蒙卡与SN耦合方法在压水堆堆腔漏束计算中应用.pdfVIP

第十一届全国蒙特卡罗方法及其应用学术交流会论文 蒙卡与 S 耦合方法在压水堆堆腔漏束计算中的应 N 用 郑 征，吴宏春，曹良志，郑友琦，张宏博，王梦琪 （西安交通大学 核科学与技术学院，陕西 西安 710049） 摘要：摘要：堆腔漏束计算属于深穿透问题，而且几何结构复杂。离散纵标方法(S )计算速度较快，但是耗费内存，存在射线 摘要摘要：： N 效应，蒙特卡罗方法几何描述精确，但是抽样效率太低，计数方差较大。为克服以上困难，并充分利用两种算法各自的 优点，本文研究了基于离散纵标方法和蒙特卡罗方法的耦合方法。主要思想是利用接口程序 DORT2MCNP 根据离散纵标方 法提供的中子角通量，计算出面源的分布概率，然后由修改后的 SOURCE 子程序生成包括上下圆面源和圆柱面源的组合 源，提供给 MCNP 程序进行计算分析。初步计算结果表明，该耦合方法是正确的，可用于堆腔漏束的计算分析。 关键词：关键词：组合源，耦合方法，DORT,MCNP,堆腔漏束，SOURCE子程序 关键词关键词：： 中图分类号：中图分类号：TL328 文献标志码文献标志码：：A 中图分类号中图分类号：： 文献标志码文献标志码：： 反应堆压力容器和一次混凝土屏蔽墙之间构成的空间一般称为堆腔。反应堆运行过程中，堆芯中的一 部分中子和光子从压力容器顺着堆腔泄漏出去，形成堆腔漏束。堆腔漏束的计算分析对于控制和降低反 应堆主回路设备的活化水平具有重要的指导意义[1-2]。 堆腔漏束计算属于深穿透问题，通常采用离散纵标方法(S )和蒙特卡罗方法进行计算分析。两种方法 N 各有优缺点。离散纵标方法计算速度较快，但是比较耗费内存，而且存在射线效应；蒙特卡罗方法几何 描述精确，采用连续点截面数据库，计算精度高，但是比较耗时。为了综合两种方法的优点，国内外的 科研人员开展了基于离散纵标方法和蒙特卡罗方法的耦合方法的研究。比如美国橡树岭国家实验室的 [3] [4] M. B. Emmett 等人开发的DOMINO 程序，可以将DOT程序（DORT 程序的早期版本）和 MORSE程序（蒙 特卡罗程序）耦合起来进行计算。DOMINO 程序的主要思想是将 DOT 程序提供的角通量转化为出射流， 然后根据出射流计算出中子随能量、空间、极角和辐角分布的归一化概率。此后美国西屋公司对 DOMINO [5-6] 程序进行了升级，并且将其应用于压水堆的屏蔽计算和中子漏束分析 。 本文基于 DOMINO 程序的主要思想，编制了 DORT程序和 MCNP[7]程序的接口程序DORT2MCNP，根据DORT 程序提供的中子角通量，计算出面源的分布概率。修改后的 SOURCE 子程序可以生成包括上下圆面源和 圆柱面源的组合源，提供给 MCNP 程序计算分析。初步数值结果表明，该耦合方法是正确的，可用于堆 腔漏束的计算分析。 1 耦合方法计算流程 耦合方法的计算流程主要分为下面三步： （1） 建立简化的二维 R-Z 几何模型，使用 DORT 程序计算得到整个区域的中子角通量分布。这一步 计算不为得到精确的计算结果，而是为后续的概率分布提供信息，所需的信息全部包括在 DORT 程序 的输出文件中。 （2 ） 确定组合面源的位置，分段数目，能群数目等参数。利用接口程序 DORT2MCNP 读取角通量等 信息，将其转化为穿过上下圆面和圆柱面的出射流，然后计算出各面上的中子随能量、空间、极角和辐 角分布的归一化概率。组合源