蒙特卡罗方法课件(清华大学-林谦).pptVIP

计算体通量的模拟方法 在实际问题中，经常遇到要计算某一区域V0 的体通量。 在通量的定义部分已经介绍过，通量可以表示为粒子各次散射对通量的贡献和。因此，下面要介绍的各种估计方法，只叙述各次散射后的通量计算方法。 计算体通量的方法主要有以下几种。 解析(统计)估计方法 粒子 n 次散射（n＝0 时为源粒子）后的通量贡献为： 其中，s1和s2分别为粒子由点rn出发，沿Ωn方向到达 区 域V0的近端和远端的交点的距离。如果点rn在V0内， 则 s1＝0。如果粒子沿Ωn方向与V0有多段相交，则 为每段相交线段的通量贡献之和。如果粒子沿 Ωn方向与V0不相交，则 。 解析估计方法就是把体通量的贡献表达式直接计算出来。当系统为均匀介质时， 如果只是V0为均匀介质，则 如果V0由多层介质组成，则需分段计算积分。 在解析估计方法中，粒子每发生一次碰撞（包括零次散射），都要记录通量的贡献值。 径迹长度方法 设粒子从第 n 次散射到第 n＋1 次散射之间走过的径迹长度为 s ，则 n 次散射的通量贡献为： 径迹长度方法就是把粒子在V0内走过的径迹长度记录下来。 下面证明，径迹长度估计是无偏的。 碰撞密度方法 设粒子从第 n 次散射到第