A peer-reviewed journal published by K. N. Toosi University of Technology

Document Type : Research Article


MASNA engineering company‎, ‎P.O‎. ‎Box 1439951113‎, ‎Tehran‎, ‎Iran


‎Numerical solution of the multi-group static forward and adjoint neutron diffusion equation (NDE) using the Finite Elements Method (FEM) is investigated in detail‎. ‎A finite element approach based on the generalized least squares method is applied for the spatial discretization of the NDE in 3D-XYZ geometry‎. ‎A computer code called GELES was also developed based on the described methodology covering linear or quadratic tetrahedral elements generated via the mesh generator for an arbitrary shaped system‎. ‎A number of test cases are also studied to validate the proposed approach‎. ‎Moreover‎, ‎to assess the output dependency to the number of elements‎, ‎a sensitivity analysis is carried out at the end‎.


  • Multi-group static 3-D neutron di usion equation is solved using the nite element method.
  • Generalized least squares FEM, through a variational approach is applied for solving the NDE.
  • GELES code is developed based on the tetrahedral elements for an arbitrary shaped system.
  • To validate the approach, output of GELES were compared against the DONJON computer code.

 Acceptable accuracy for the neutron multiplication factor and the power distribution was achieved.


Ackroyd, R. (1986a). A finite element method for diffusion theory embracing nodal and difference methods. Progress in Nuclear Energy, 18(1-2):7–20.
Ackroyd, R. T. (1986b). Generalized least squares as a generator of variational principles and weighted residual methods for FEM transport methods. Progress in Nuclear Energy, 18(1-2):45–62.
Bell, G. I. and Glasstone, S. (1970). Nuclear Reactor Theory. Technical report, Division of Technical Information, US Atomic Energy Commission. Cavdar, S. and Ozgener, H. A. (2004). A finite element boundary element hybrid method for 2-D neutron diffusion calculations. Annals of Nuclear Energy, 31(14):1555–1582.
Center, A. C. (1977). ANL Benchmark Book-Report ANL-7416. Argonne National Laboratory, Argonne, IL.
Duderstadt, J. J. and Hamilton, L. J. (1976). Nuclear Reactor Analysis, volume 1. Wiley New York. H´ebert, A. (2008). A Raviart–Thomas–Schneider solution of the diffusion equation in hexagonal geometry. Annals of Nuclear Energy, 35(3):363–376.
H´ebert, A., Sekki, D., and Chambon, R. (2013). A User Guide for DONJON Version4. ´Ecole Polytechnique de Montr´eal Montr´eal QC, Canada, Tech. Rep. IGE-300.
Kang, C. M. and Hansen, K. F. (1973). Finite element methods for reactor analysis. Nuclear Science and Engineering, 51(4):456–495.
Kolev, N., Lenain, R., and Magnaud, C. (1999). AER Benchmark Specification Sheet. Test ID.: AER–FCM, 101. Lamarsh, J. R. (1975). Introduction to Nuclear Engineering. Lewis, E. E. (1981). Finite element approximation to the even-parity transport equation. In Advances in Nuclear Science and Technology, pages 155–225. Springer.
McConnell, A. J. (1951). The hypercircle method of approximation for a system of partial differential equations of the second order. In Proceedings of the Royal Irish Academy.
Section A: Mathematical and Physical Sciences, volume 54, pages 263–290. JSTOR.
Schulz, G. (1996). Solutions of a 3D VVER-1000 Benchmark. In Proc. 6-th Symposium of AER on VVER Reactor Physics and Safety, Kirkkonummi, Finland.
Wang, Y., Bangerth, W., and Ragusa, J. (2009). Threedimensional h-adaptivity for the multigroup neutron diffusion equations. Progress in Nuclear Energy, 51(3):543–555.
Zienkiewicz, O. C. and Taylor, R. L. (2005). The finite element method for solid and structural mechanics. Butterworth heinemann