点探测器通量新方位角偏倚抽样计算方法.pdfVIP

点探测器通量的新方位角偏倚抽样计算方法 范佳锦1 裴鹿成2 (1．清华大学工程物理系2．中国原子能科学研究院) 摘要MonteCarlo计算点探测器通量的方法很多，但真正达到可应用目标的仅有指向概率 方法和环探测器方法两种。环探测器方法能解决指向概率方法的方差无限问题，但它在计算时要 求系统必须满足环对称条件，应用范围极窄。将方位角进行偏倚抽样计算，虽然估计量仍然是无 界的，但无界的阶与环探测器方法相同，并且方差有限，更为重要的是，它适合所有情况的点探测 器通量计算问题。本文给出了一种新的方位角偏倚抽样方法，相对于原先方法，能有效减小计算 量，从而获得更高的计算效率。 一、引 言 Carlo方法的应用问题中占有非常重要的地位，先后发展出很多方法， 计算点通量在Monte 但真正达到可应用目标的仅有指向概率方法，或称“下次事件估计方法”(简记为NEE方法)和 Carlo程序MCNP中，除提供这两种方法外，还曾提供过 环探测器方法两种【l。]。在超级Monte “多一次碰撞通量估计方法”(简记为OMCFE方法)，实际应用过程中，因效果不佳，在MCNT3B 及其以后的版本中被删除掉了，目前国外还在工作希望能对其进行改进H。，但是效果仍不甚理 想。指向概率方法与环探测器方法的估计量都是无界的，但是，前者无界的阶为O(1／R2)(R 为碰撞点到探测点的距离)，后者无界的阶为0(1／尺幽)(R舶为碰撞点到探测器环的最小距 Carlo估计的收敛速度，后者的方差有限；前者适合 离)；前者的方差无限，直接影响相应Monte 所有情况下的点探测器通量计算问题，后者却只适合满足环对称情况下(环探测器上每一点的 通量都等于点探测器的通量)的点探测器通量计算问题，使其应用范围受到了极大限制。方位 角偏倚抽样方法的基本原理是，利用粒子迁移距离已抽样确定和散射角在(0，23)上的均匀性， 确定这样一个环，粒子具体是在环上哪一点发生碰撞，采用偏倚抽样方法，即对方位角采用偏 倚抽样积分法，使其距离点探测器越近者越多抽样，从而提高计算结果的收敛速度。本文利用 空间立体角问的相互关系，给出了一种新的方位角偏倚抽样方法，相对于原先方法，能有效减 小计算量，从而提高了计算结果的收敛速度。 二、点探测器通量的数学期望表达式 不影响本文对计算点探测器通量方位角偏倚方法的讨论，在这将只考虑由源出发经过m 次(m1)散射后发射的一个粒子(中子或光子)，再经过一次散射对点探测器产生的通量贡 献。Fm(r。)(图1)。 根据L(r。)的定义，可以将其表示成如下数学期望形式： 一20一 图1粒子的随机游动历史 獬 ；。+。)d‰×-『P(c。s臼州I；川)dcos‰， ×艿(；m+z一；a)南，r(tm+-I；m+，，二m+·)dtm+， 对于随机变量9。 于 (2) “+0|+0 瓮赢； 射一一一 籽时一∑。 的数学期望。 式中 rm——由源出发经m次散射后发射粒子的位置； 叫。——由源出发经过m次散射后发射粒子的方向矢量； 9。——由源出发经过m次散射后发射粒子的方位角，即平面(叫。一。，叫。)与平面(∞。一，， rd—rm)问的二面角； cos口。——由源出发经过m次散射的粒子，散射前后间夹角余弦； ￡。——由源出发经过m次散射后发射的粒子，自由飞行的距离。 三、计算点探测器通量的方位角偏倚抽样方法 定义+rm处的坐标为(茹。，YoZO)，4rm+1处的坐标为(z，,Yl，z，)。为了方便起见，在绝对坐标 系。彬的基础上建立相对坐标系O魍z，使其原