Document Type : Research Article

**Author**

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

**Abstract**

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.

**Highlights**

- Multi-group static 3-D neutron diusion 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.

**Keywords**

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