Solution of neutron transport equation using Daubechies’ wavelet expansion in the angular discretization

Daubechies’ wavelet expansion is introduced to discretize the angular variables of the neutron transport equation when the neutron angular flux varies very acutely with the angular directions. An improvement is made by coupling one-dimensional wavelet expansion and discrete ordinate method to make two-dimensional angular discretization efficient and stable. The angular domain is divided into several subdomains for treating the vacuum boundary condition exactly in the unstructured geometry. A set of wavelet equations coupled with each other is obtained in each subdomain. An iterative method is utilized to decouple the wavelet moments. The numerical results of several benchmark problems demonstrate that the wavelet expansion method can provide more accurate results by lower-order expansion than other angular discretization methods.