ntkFoam: An OpenFOAM based neutron transport kinetics solver for nuclear reactor simulation

Due to the complexity of detailed nuclear reactor numerical simulation, especially the complicated geometry and multi-physics coupling properties of the advanced reactor, a nuclear reactor kinetics solver for nuclear reactor engineering and design needs to be developed. Based on the open source C++ software OpenFOAM, this work establishes a neutron transport kinetics solver, namely ntkFoam, for the nuclear reactor kinetics simulation, from the governing equations to the detailed implementations. The coupling between multi-group neutron transport equations and the delayed neutron precursor balance equations are considered with using the finite volume method. By introducing the inverse power method to the OpenFOAM, the k -eigenvalue problems are calculated, and by coupling the neutron transport and delayed neutron precursor calculations using the Euler implicit scheme, the transient kinetics problems are simulated. Numerical results show that the proposed ntkFoam can simulate the multi-group neutron transport kinetics problem accurately and flexibility, and both the regular and irregular mesh configurations can be adopted. This work can provide some new perspectives and foundations to nuclear reactor coupling calculations.