Finite element based composite solution for neutron transport problems

A finite element treatment for solving neutron transport problems is presented. The employs region-wise discontinuous finite elements for the spatial representation of the neutron angular flux, while spherical harmonics are used for directional dependence. Composite solutions has been obtained by using different orders of angular approximations in different parts of a system. The method has been successfully implemented for one dimensional slab and two dimensional rectangular geometry problems. An overall reduction in the number of nodal coefficients (more than 60% in some cases as compared to conventional schemes) has been achieved without loss of accuracy with better utilization of computational resources. The method also provides an efficient way of handling physically difficult situations such as treatment of voids in duct problems and sharply changing angular flux. It is observed that a great wealth of information about the spatial and directional dependence of the angular flux is obtained much more quickly as compared to Monte Carlo method, where most of the information in restricted to the locality of immediate interest. (author)