二维柱坐标系中子输运问题的边界型算法

传统的确定论方法在研究高维度的中子输运方程时需要求解大型矩阵的逆矩阵,该过程将增加计算量、降低计算效率,因此本文采用了一种新型的边界型方法——半边界法(the Half Boundary Method)来解决此类问题。半边界方法利用已有的中子输运方程推导得出相邻两点之间的中子通量密度关系通式,然后迭代地应用此通式求出模型内部任一点与边界点之间的中子通量密度关系,即可得到整个模型的中子通量密度分布。半边界方法在已知一半边界条件处的中子通量密度即可求得模型内部任一点的中子通量密度,提高了计算效率。本文将半边界法在真空边界条件模型中得出的结果与NECP-MCX软件得到的数值结果进行比较,通过分析发现不同空间变量和角度变量离散数量对半边界法数值结果的影响,当空间变量M = 20,N = 40和角度变量H = 50,K = 32时,半边界法的平均误差为0.76%,具有较高精确度。