球几何中子输运保正线性间断有限元格式-强激光与粒子束.PDF

第 卷第 期 强 激 光 与 粒 子 束 , 29 7 Vol.29 No.7 年 月 , 2017 7 HIGH POWERLASERANDPARTICLEBEAMS Jul. 2017 球几何中子输运保正线性间断有限元格式* , , 洪振英 袁光伟 魏军侠 ( , ) 北京应用物理与计算数学研究所 北京 100094 : , 摘 要 针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题 构造了保正线性间断有 , 。 , 限元格式 该格式保持中子角通量 阶矩和 阶矩 现有方法计算中子角通量非负时 采用传统的线性间断有 0 1 , ; , , 。 限元方法 求解线性方程组 原方法计算出现负通量 则采用构造的保正格式 求解非线性方程组 编制了球几 , 。 , 何中子输运问题保正格式程序模块 并集成到应用程序 数值算例表明构造的保正格式计算的中子通量非负 , 。 有效降低数值误差 提高数值计算的精度 : ; ; ; 关键词 球几何 输运方程 保正格式 间断有限元方法 中图分类号: 文献标志码: : / O241.82 A doi10.11884HPLPB201729.160320 [ ] 1-2 、 、 、 、 中子输运方程广泛应用于核反应堆设计与分析 惯性约束聚变 临界性安全 医学理疗和图像 核武器 。 , , 安全以及油井探测等领域 多物理耦合计算中 输运问题通常需要考虑多介质计算 间断有限元格式适用于多 。 , ( , ) 介质问题的数值模拟 当相邻的轻重介质性质差别较大时 现有的间断有限元格式 DG LD计算的中子通量 , 。 不能确保为正 这违背中子通量必须非负的物理要求 通量计算出现负值会导致辐射流体与中子输运耦合问