Development of a consistently formulated general order nodal method for solving the three-dimensional multi-group neutron kinetic equations

A new general high order consistent nodal method for solving the 3-D multigroup neutron kinetic equations in (x-y-z) geometry has been derived by expanding the flux in a multiple polynomial series for the space variables without the quadratic fit approximations of the transverse leakage and for the time variable and using a weighted-integral technique. The derived equation set is consistent mathematically, and therefore, we can expect very accurate solutions and less computing time since we can use coarse meshes in time variable as well as in spatial variables and the solution would converge exactly in fine mesh limit