第五部分蒙特卡罗方法在计算机上的实现教学课件教学文稿.pptVIP

中子非弹性散射后能量的确定 中子非弹性散射后，能量与质心系散射角θC的关系是： 其中， 为第 K 个能级的阈能， 为第 K 个能级的激发态能量。 如果确定了实验室系散射角θL，则根据下式 确定 后，再计算出 Em+1。 光子康普顿（Compton）散射后能量的确定 光子发生康普顿散射后，其能量分布密度函数为： 其中， K(α) 为归一因子。 ， 和 分别为光子散射前后的能量，以 m0c2 为单位，m0为电子静止质量，c 为光速。 光子康普顿散射能量分布的抽样方法为： x 的抽样确定后，散射后的能量为： ＞ ＞ ＞ ≤ ≤ ≤ 碰撞后散射角的随机抽样 粒子碰撞后运动方向Ωm+1的确定，一般与散射角有关。由已知分布抽样确定散射角后，再确定Ωm+1。常见的散射角分布有如下几种： 质心系各向同性分布 散射角在质心系服从各向同性分布时，其抽样方 法为 。质心系散射角θC抽样确定后， 需转换成实验室系散射角θL。 在中子弹性散射情况下，转换公式为： 其中 A 为碰撞核质量， 。 在中子非弹性散射情况下，转换公式为： 其中， 为第 K 个能级的阈能。 中子弹性散射勒让德 (Legendre) 多项式分布 中子弹性散射角分布常以勒让德多项式的展开形式给定。散射角余弦 x＝cosθ的分布密 度函数为： 其中 Pl(x) 为 l 阶勒让德多项式。 该分布即为 n 阶勒让德近似展开。 勒让德多项式由以下递推公式确定： 考虑新的分布： 当选取 x0，x1，… xn 为 Pn+1(x)＝0 的根，且 时，fa(x) 依照勒让德多项式展开的前 n 项与 f (x) 的展开形式相同。因此，可以用 fa(x) 作为 f (x) 的近似分布。 在实际问题中，由于勒让德多项式展开项数不够，可能出现某个 为负值的现象。此时可以采用如下近似分布： 其中： 对于该近似分布，可用加抽样方法进行抽样： 此时，由于偏倚抽样而引起的纠偏因子为 wK ，也就是说，粒子的权重要乘上wK。 光子康普顿散射角分布 光子的康普顿散射角与其散射前后的能量有关 , 它的分布密度函数为： 抽样方法为： 碰撞后运动方向Ωm+1的确定 实验室系散射角θL确定后，依据不同的坐标系的表现形式，有不同的确定方法。 确定方向余弦 um+1，vm+1，wm+1 其中， 方位角 在 [0, 2π] 上均匀分布。 当 时，不能使用上述公式，可用下面的简单公式： 确定Ωm+1的球坐标 (θm+1，φm+1) 设Ωm的球坐标分别为 (θm,φm)，其中，θ为粒子运动方向与 z 轴的夹角， φ为粒子运动方向在 x y 平面上投影的方位角。则Ωm+1的球坐标 (θm+1,φm+1) 分别由下式确定： 球形几何的随机游动公式 一般几何的随机游动公式可以应用到球形几何，而对球对称问题，使用特殊形式更为方便。 下次碰撞点的径向位置 rm+1的确定 两次碰撞点间的距离 L 确定之后，下次碰撞点的径向位置 rm+1的计算公式为： 设系统的外半径为R，如 rm+1≥R，则粒子逃出系统。 粒子碰撞后瞬时运动方向的确定 在球对称系统中，粒子运动方向用其与径向夹角余弦来描述。使用球面三角公式，粒子碰撞后瞬时运动方向与径向夹角余弦 cosθm+1的计算公式为： 其中， 为在 [0, 2π] 上均匀分布的方位角， 为在 rm+1 点进入碰撞前瞬时运动方向与 rm+1 径向之间的夹角。 点到给定边界面的距离 在抽样确定输运距离、判断粒子是否穿透系统时， 常遇到求由 rm 出发，沿Ωm 方向到达某个区域表面的距离问题。在记录对结果的贡献时，也常使用类似的量。区域表面通常是平面或二次曲面。 求到达区域表面的距离问题，实际上是求直线（或半射线）与平面或二次曲面的交点问题。这是 蒙特卡罗方法解粒子输运的各种实际问题时 , 所遇到的基本几何问题。 点到平面的距离 点 沿方向 的直线方程为： 该直线到达方程为 的平面的距离为： 当 与平面平行时，即 直线与平面无交点。 如果 l 为负值，直线与平面也无交点。这时，粒子的