A discontinuous finite element method based on B-spline wavelet on the interval for solving first-order neutron transport equation with discrete ordinate (SN) angular discretisation

2021 ◽  
Vol 162 ◽  
pp. 108511
Author(s):  
Hewei Guo ◽  
Wei Chen ◽  
Xinbiao Jiang ◽  
Duoyu Jiang ◽  
Xinyi Zhang
Author(s):  
Chao Fang ◽  
Hongchun Wu ◽  
Liangzhi Cao ◽  
Yunzhao Li

This paper presents a fast sub-grid scale (SGS) finite element method for the first order neutron transport equation. The spherical harmonics method is adopted for the angular discretization. The sub-grid scale discretization embeds discontinuous component in each element to provide a stabilization term for the continuous finite element formulation. Traditional SGS method uses Riemann decomposition and vacuum boundary assumption to decouple the discontinuous component. Here we propose a new method to perform the decoupling based on the assumption that the convection term of the discontinuous component is proportional to the residual of angular flux in each element. The computing costs for the establishment of the coefficient matrix of discontinuous component are reduced to O(1) from O(n3). Further more, the computing costs for the inversion of the coefficient matrix are reduced to O(n) from O(n3) by applying mass lumping technique. Numerical results show that the new method is not only more efficient but also yields more accurate solution than traditional sub-grid scale method.


Sign in / Sign up

Export Citation Format

Share Document