质子在重核内级联过程的蒙特卡罗计算

用更逼真的核模型(费米型的核子密度分布)进行了模拟计算,计算中用折线近似粒子的曲线运动轨迹,并考虑了短子在核表面处的折射和反射。用自由核子截面确定核内的碰撞,并且自始至终地使用了统计抽样技术。入射粒子的能量被限于350Mev 以下,即π介子不大可能产生的能量范围。初步计算结果与实验数据进行了对比。     

用S_N方法求解二维动态多速中子迁移方程

本文简要地介绍了二维动态多速中子迁移方程的S_N解法,并给出一些计算结果。为了验证方法的可靠性,我们计算了七个临界装置。计算结果表明:在无反射层的均匀介质中,对角度的网点采取4×7,空间网点采取7×7以上,就可以得到满意的结果,不必作通量校正。但在加上反射层时,空间网点需增加到12×8。在YH1机上计算,对无反射层均匀介质,计算每个时间分点,作通量校正时需2.7-3.6秒,不作通量校正需1.5-2.1秒。程序向量化后计算每个时间分点需0.48秒以下,提高速度5倍左右。本计算程序可以给出瞬时中子通量的空间分布、能谱分布及角分布,同时给出装笠参数,时间常数λ,有效倍增因子K_6等。     

空间辐射环境中的辐射效应

天然空间辐射环境与辐射效应基本机制是空间飞行器抗辐射加固研究中的两个关键问题, 本文总结了天然空间辐射环境的总体性质, 并对辐射效应基本机制进行了简要的分析, 指出了目前空间飞行器抗辐射加固研究的重点。     

相空间有限元方法在解二维中子输运方程中的应用

本文采用相空间有限元方法求解了柱形临界多群中子输运问题。其中对于方程中的坐标变量用分片连续线性多项式作为试探函数,对于方程中的角度变量用分片连续双线性多项式作为试探函数。整个求解空间区域和角度区域分别采用三角形和矩形单元划分,然后利用迦辽金方法得到一个以网格点处角通量为未知数的线性联立代数方程组,方程组中的系数矩阵的存储采用了压缩存储技术。最后用高斯消元法解此有限元方程组,表明相空间有限元方法计算收敛性较好、计算精度高。     

电子在材料中输运的蒙特卡罗模拟

介绍了电子输运问题的蒙特卡罗方法的基本原理和步骤。电子输运过程的模拟采用了“浓缩历史”蒙特卡罗方法 ,电子随机游动步长根据能量步划分 ,在每个电子游动步内 ,电子的能量损失为非辐射 (碰撞 )能量损失和辐射 (韧致辐射 )能量损失之和 ,其中非辐射能量损失通过对Landau公式的随机抽样得到 ,辐射能量损失以及韧致辐射光子的产生通过利用Bethe Heitler的随机抽样处理 ,在每个电子游动步结束时电子的出射角的抽样利用Gouddsmit Saunderson多重散射公式 ,电子 -电子的散射利用Moller截面公式进行抽样。最后对电子在航天用合金铝平板材料中的输运进行了模拟计算并给出了计算结果     

二维非定常中子输运的并行算法

编制了求解时间相关二维中子输运问题的并行计算程序,并在YH-1计算机上进行过许多验证计算。对六个临界装置的计算表明,临界参数的计算值与实验值之间的偏差在参数值的实验误差范围之内,这证明计算程序是正确的。不同并行算法的计算结果完全相同也证实了程序编制的正确性。计算表明,并行程序的计算效率比相同参数下的串行程序提高约30倍以上。为了尽可能提高并行计算的效率,除了采取通常的向量化技巧外,在算法上提出了并维技巧和斜对角扫描的方法。     

均匀平面γ放射源上空γ射线的能谱和角分布

本文用蒙特·卡罗方法计算了在光滑地表面上均匀分布的γ放射源上空1米、300米、500米三个高度处的γ辐射场。对于1.25Mev和0.67Mev两种源能量(分别相应于Co~(60)和Cs~(137)的γ射线能量),得到了辐射场上空γ光子通量的能量分布和角分布。对于每种源能量,计算了0.04～1.5Mev之间8个能量间隔和0°～180°之间9个角度间隔的散射通量;单独计算了这种单能源上空各高度处的直射束通量,并与测量的和国外计算的数据进行了对比。     

带电粒子束传输中发散范围的计算

本文讨论了带电粒子束在真空环境中的定常传输。对束的径向发散提出了一个计算方法。根据数值计算结果,提出一个计算束半径的近似公式。