Abstract

It is difficult to control discretization errors and reduce computational cost at the same time for multi-scale neutron transport problems, and the adaptive mesh refinement technique is one of the powerful and effective methods for high-resolution transport calculation. We propose a multilevel mesh adaptivity algorithm and develop a discrete ordinates calculation framework on 3-D Cartesian meshes. Data management and transport sweep are optimized based on the multilevel pyramid data structure containing pointer-type and array-type variables. A new spatial-moment-ratio error indicator is derived to measure the leading term of distribution functions and to drive the local mesh refinement. The numerical properties of spatial moment factors are better than the normalized gradient of angular flux. Compared with uniform refinement, our adaptivity algorithm reaches the same accuracy level with a computational cost saving of 60%-80% for heterogeneous problems.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call