蒙特卡罗方法在核技术中的应用.ppt

常用的通用蒙特卡罗程序简介 （1）MORSE程序 较早开发的通用蒙特卡罗程序，可以解决中子、 光子、中子－光子的联合输运问题。采用组合几何 结构，使用群截面数据，程序中包括了几种重要抽 样技巧，如俄国轮盘赌和分裂技巧，指数变换技 巧，统计估计技巧和能量偏移抽样等。程序提供用 户程序，用户可根据需要编写源分布以及记录程序。 常用的通用蒙特卡罗程序简介 （2）EGS程序 EGS是Electron-Gamma Shower 的缩写，它是 一个用蒙特卡罗方法模拟在任意几何中，能量从几个 KeV到几个TeV的电子-光子簇射过程的通用程序包。 由美国Stanford Linear Accelerator Center提供。 EGS于1979年第一次公开发表，提供使用。EGS4是 1986年发表的EGS程序的最新版本。 常用的通用蒙特卡罗程序简介 （3）MCNP程序 MCNP是美国Los Alamos国家实验室开发的大型 多功能通用蒙特卡罗程序，可以计算中子、光子和 电子的联合输运问题以及临界问题，中子能量范 围从10-11MeV至20MeV，光子和电子的能量范围从 1KeV至1000MeV 。程序采用独特的曲面组合几何 结构，使用点截面数据，程序通用性较强，与其 它程序相比，MCNP程序中的减方差技巧是比较多 而全的。 * * * * * * * * * * * * * * * * * * * * * * 4.1　加抽样方法 对如下加分布给出一种抽样方法： 其中Pn≥0， ，fn(x)为与参数n有的分布密 度函数，n=1，2，…。 加分布的抽样方法为:首先抽样确定n’，然后由 fn’(x)中抽样x，即： 4. 随机抽样的一般方法 如下形式的分布称为乘分布： 　　其中H(x)为非负函数， f1(x)为任意分布密度 函数。 　　令M为H(x)的上界，乘抽样方法如下： 4.2　乘抽样方法 ≤ ＞ 　　在实际问题中，经常会遇到如下形式的分布： 其中Hn(x)为非负函数，fn(x) 为任意分布密度函 数，n=1，2，…。不失一般性，只考虑n=2的情况： 将 f(x) 改写成下式的加分布形式： 4.３　乘加抽样方法 其中： 乘加抽样方法为： ≤ 其中： 例7.光子散射后能量分布的抽样 　　令光子散射前后的能量分别为 和 （以 m0c2 为单位，m0为电子静止质量，c 为光速）， 则 x 的分布密度函数为： 　　该分布即为光子散射能量分布，它是Klin－ Nishina 公式确定的。其中 K(α) 为归一因子： 把光子散射能量分布改写成如下形式： 在［1,1+2α］上定义如下函数: 则有 使用乘加抽样方法： ax ＞ ＞ ? 4 27 ＞ ≤ ≤ ≤ 2 2 2 2 1 1 1 3 2 3 2 1 2 3 2 ) 1 ( 2 1 29 4 27 1 1 2 1 2 1 2 1 f f f f f f f f f X X X X X X X X X = - ￡ + = + ￡ = ú ú ? ù ê ê ? é + ? ? ? è ? - + ￡ + + = x ax a x a a x a 光子散射能量分布的抽样方法为： 例8. 散射方位角余弦分布的抽样 散射方位角φ在［0,2π］上均匀分布，则其正 弦和余弦sinφ和cosφ服从如下分布： 直接抽样方法为： 令φ=2θ，则θ在［0,π］上均匀分布，作变换 其中0≤ρ≤1，0≤θ≤π，则 (x,y) 表示上半个单位圆内的点。 如果 (x,y) 在上半个单位圆内均匀分布，则 θ在［0π］上均匀分布，由于 因此抽样sinφ和cosφ的问题就变成在上半个单位圆内均匀抽样 (x,y) 的问题。 为获得上半个单位圆内的均匀点，采用挑选 法，在上半个单位圆的外切矩 形内均匀投点（如图）。 舍弃圆外的点，余下的就 是所要求的点。 抽样方法为： 为了确定各向同性散射 ， 根据公式： 对于各向同性散射，θ在［0,π］上均匀分 布，φ在［0,2π］上均匀分布。 则: 例9. 各向同性散射方向余弦的抽样 因此，抽样确定 （另类）抽样方法： 定义两个随机变量： 　　可以证明，当 　时，随机变量 x 和 y 服从如下分布: 定义区域为： 则w＝cosθ的分布可以用上述分布表示成积分分 布的形式： 令 ，则属于上述