蒙特卡罗方法在ADS屏蔽计算中的应用.pdfVIP

第 25 卷 第 2 期 核 动 力 工 程 Vol. 25. No.2 2 0 0 4 年 4 月 Nuclear Power Engineering Apr. 2 0 0 4 文章编号：0258-0926(2004)02-0106-04 蒙特卡罗方法在ADS屏蔽计算中的应用 廖义香 （南华大学核科学技术学院，湖南，衡阳，421001 ） 摘要：利用蒙特卡罗方法计算了新一代核能系统加速器驱动系统(ADS)中质子束管内的中子归一化注量 率分布以及通过质子束管入口和其它外表面逸出的归一化中子注量率，得出了一些对ADS 系统的设计有重 要意义的结论。 关键词：蒙特卡罗方法；加速器驱动系统(ADS)；屏蔽计算 + 中图分类号：TL329 . 2 文献标识码：A 1 前 言 后通过对事件的模拟“试验”得到该事件出现的 随着近代加速器技术的进展，借助强流中能 频率，并用它近似代替事件出现的概率，从而得 质子加速器驱动的次临界核反应堆系统(ADS) 到问题的解。用统计学的术语来表述就是把所求 愈来愈引起学术界的关注。ADS 是外源中子驱 问题转换成某个随机变量的期望值，通过模拟“试 动次临界反应堆，为一种新型的核能系统，其中 验”方法，获得这个随机变量的平均值。 子余额比同条件下的临界反应堆要高。这种由外 可以把蒙特卡罗方法解题步骤归结为： 源中子导致的中子余额增量有利于嬗变核废料 (1)构造或描述概率过程：对于本身就具有随 或增加核燃料的转换，因此，ADS 又称放射性 机性质的问题，如粒子输运问题，主要是正确描 洁净系统，能够克服核能系统固有的环境和资源 述和模拟这个概率过程；对于本来不是随机性质 [1] 的确定性问题，比如计算定积分、解线性方程组 问题，从而使核能成为可持续发展的能源 。 蒙特卡罗方法(Monte Carlo)又称随机抽样技 和偏微分方程边值问题，要用蒙特卡罗方法求解， 巧或统计实验方法，与基于确定论的一般数值方 就必须事先构造一个人为的概率过程，它的某些 法有着本质的区别。它能够较逼真地描述事物的 参量正好是所要求问题的解，即要将不具有随机 特点以及物理实验的真实过程，故解决问题与实 性质的问题转化为随机性质的问题。 际非常符合，可以得到很圆满的结果。具有受几 (2)实现从已知概率分布抽样：构造了概率模 何条件限制少等诸多优点，因此能够解决一些一 型以后，由于各种概率模型都可以看作是由各种 般数值方法难以解决的三维复杂问题。 各样的概率分布构成的，因此产生已知概率分布 本文利用蒙特卡罗方法计算了 ADS 系统中 的随机变量(或随机向量) ，就成为实现蒙特卡罗 质子束管内的归一化中子注量率分布，以及通过 方法模拟实验的基本手段，这也是蒙特卡罗方法 质子束管入口和其他外表面溢出的归一化中子注 被称为随机抽样的原因。随机数是实现由已知分 量率，得出了一些对ADS 系统的设计有意义的结 布抽样的基本工具。在由已知分布抽样过程中， 论。 将随机数作为已知量，用适当的数学方法可以由 它产生具有任意已知分布的简单子样。 2 蒙特卡罗方法的基本思想