中子迁移方程有限无迭代算法

摘要

本文讨论中子迁移方程Ａｕ＝Ｆ的数值方法，这里Ａ是迁移算子，Ａ＝Ｂ＋Ｊ，而Ｂ＝Ω．Ｖ＋σＩ，Ω是中子运动方向单位矢量，△是梯度算子，Ｉ是恒等算子，σ是横载面积，Ｊ是积分算子，Ｆ是外中子源。本文提出上述方程的迭代算法：Ｂ－ｕ＾（ｎ＋１）＝－Ｊｊ＾（ｎ）＋Ｆ，结合不连续Ｇａｌｅｒｋｉｎ方法给出迭代有限元算法，讨论了算法的收敛性和误差估计，给出计算实例，理论分析和计算结果表明此方法结构简单，计算方便和精度