Development of Three-Dimensional Neutron Kinetics Code Based on High Order Nodal Expansion Method in Hexagonal-Z Geometry

To solve three-dimensional kinetics problems, a high order nodal expansion method for hexagonal-z geometry (HONEM) and a Runge-Kutta (RK) method are respectively adopted to deal with the spatial and temporal problem. In the HONEM, 1D partially-integrated flux are approximated by using four order polynomial. The two order polynomial is adopted to the approximation of partially-integrated leakages. The Runge-Kutta method is adopted as a tool for dispersing the time term of 3D kinetics equation. A flux weighting method (FWM) is used for obtaining homogenized cross sections of mix node. The three-dimensional hexagonal kinetics code has been developed based on this method and tested with two benchmark problems of VVER which are the control rod ejection without any feedback and with simple adiabatic Doppler feedback. The results calculated by this code agree well with the reference results and the code is validated.