摘要
离散纵标法作为求解厚扩散极限中子输运问题的重要方法之一,其常用的空间离散格式,如有限差分格式极易在光学厚介质中产生数值扩散,并且粗网格精度不足和难以适配复杂几何等问题使得离散纵标法的应用具有一定的局限性。本研究采用伽辽金方法推导弱形式或变分形式的离散纵标方程,基于间断有限元思想构造高阶曲网格下的拉格朗日有限元基函数,建立中子输运方程的高阶间断有限元离散格式。选取了构造解算例、国际原子能机构(IAEA)EIR-2基准题和厚扩散极限算例进行建模与输运计算,对该空间离散格式的计算精度、收敛性进行了测试验证并分析了其厚扩散极限特性。数值结果表明,多介质问题中高阶间断有限元格式的计算结果与基准值相对误差小于1%,且在曲网格下也具有较高的计算精度和符合预期的收敛速度。该离散格式可有效解决曲网格下的厚扩散极限中子输运问题,具有较好的数值特性,且在光学厚扩散极限下具备扩散极限渐近保持性质。
Discrete ordinate method is one of the main numerical methods for solving the problem of neutron transport with thick diffusion limit.Its commonly used spatial discrete schemes,such as finite difference scheme,are easy to cause numerical diffusion in optical thick media,and the application of discrete ordinate method has certain limitations due to the lack of coarse mesh accuracy and difficulty in adapting to complex geometry.In this paper,Galerkin method is used to derive discrete ordinate equations in weak form or variational form.Based on the idea of discontinuous finite element,the Lagrangian finite element basis function in the higher-order curved meshes is constructed to establish the higher-order finite element discrete scheme of the transport equation.Two manufactured solution examples,the IAEA EIR-2 benchmark problem and the thick diffusion limit example,are selected for modeling and transport calculation,and the calculation accuracy and convergence of the spatial discrete scheme are tested and verified,and its thick diffusion limit characteristics are analyzed.The numerical results show that the relative error between the calculated results of high-order discontinuous finite element scheme and the reference value is less than 1%,and it also has high calculation accuracy and convergence rate in curved meshes.The discrete scheme can effectively solve the problem of neutron transport with thick diffusion limit in curved meshes.It has good numerical characteristics,and has the asymptotic preservation property of diffusion limit under optical thick diffusion limit.
作者
王新宇
张斌
陈义学
Wang Xinyu;Zhang Bin;Chen Yixue(School of Nuclear Science and Engineering,North China Electric Power University,Beijing,102206,China)
出处
《核动力工程》
EI
CAS
CSCD
北大核心
2023年第4期41-48,共8页
Nuclear Power Engineering
基金
国家自然科学基金(11975097)。
关键词
中子输运
离散纵标法
空间离散格式
间断有限元
厚扩散极限
Neutron transport
Discrete ordinate method
Spatial discretization
Discontinuous finite element
Thick diffusion limit
