Browse > Article
http://dx.doi.org/10.1016/j.net.2015.10.009

Development of Galerkin Finite Element Method Three-dimensional Computational Code for the Multigroup Neutron Diffusion Equation with Unstructured Tetrahedron Elements  

Hosseini, Seyed Abolfazl (Department of Energy Engineering, Sharif University of Technology)
Publication Information
Nuclear Engineering and Technology / v.48, no.1, 2016 , pp. 43-54 More about this Journal
Abstract
In the present paper, development of the three-dimensional (3D) computational code based on Galerkin finite element method (GFEM) for solving the multigroup forward/adjoint diffusion equation in both rectangular and hexagonal geometries is reported. Linear approximation of shape functions in the GFEM with unstructured tetrahedron elements is used in the calculation. Both criticality and fixed source calculations may be performed using the developed GFEM-3D computational code. An acceptable level of accuracy at a low computational cost is the main advantage of applying the unstructured tetrahedron elements. The unstructured tetrahedron elements generated with Gambit software are used in the GFEM-3D computational code through a developed interface. The forward/adjoint multiplication factor, forward/adjoint flux distribution, and power distribution in the reactor core are calculated using the power iteration method. Criticality calculations are benchmarked against the valid solution of the neutron diffusion equation for International Atomic Energy Agency (IAEA)-3D and Water-Water Energetic Reactor (VVER)-1000 reactor cores. In addition, validation of the calculations against the $P_1$ approximation of the transport theory is investigated in relation to the liquid metal fast breeder reactor benchmark problem. The neutron fixed source calculations are benchmarked through a comparison with the results obtained from similar computational codes. Finally, an analysis of the sensitivity of calculations to the number of elements is performed.
Keywords
Finite element method; Forward/adjoint fluxes; Galerkin; Gambit; Power distribution; Unstructured tetrahedron elements;
Citations & Related Records
연도 인용수 순위
  • Reference
1 L.F. Richardson, The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam, Philos. Trans. R. Soc. Lon. A 210 (1911) 307-357.   DOI
2 K. Ivanov, M. Manolova, T. Apostolov, An effective solution scheme of a three-dimensional reactor core model in hexagonal geometry, Comput. Phys. Commun. 82 (1994) 1-16.   DOI
3 V.G. Zimin, D.M. Baturin, Polynomial nodal method for solving neutron diffusion equations in hexagonal-z geometry, Ann. Nucl. Energy 29 (2002) 1105-1117.   DOI
4 A. Hebert, Development of the nodal collocation method for solving the neutron diffusion equation, Ann. Nucl. Energy 14 (1987) 527-541.   DOI
5 K.S. Smith, An Analytic Nodal Method for Solving the Two-Group, Multidimensional, Static and Transient Neutron Diffusion Equations, Massachusetts Institute of Technology, Cambridge, United States, 1979.
6 T.J. Barth, Aspects of unstructured grids and finite-volume solvers for the Euler and Navier-Stokes equations, AGARD, Special Course on Unstructured Grid Methods for Advection Dominated Flows 61, NASA Ames Research Center, Moffett Field, CA, United States, 1992, pp. 18-34. SEE N92-27671.
7 R. Li, Z. Chen, W. Wu, Generalized Difference Methods for Differential Equations: Numerical Analysis of Finite Volume Methods, CRC Press, United States, 2000.
8 S. Iqbal, An Adaptive Finite Element Formulation of the Boltzmann-Type Neutron Transport Equation, Ph.D thesis, Ghulam Ishaq Khan Institute of Engineering Sciences and Technology, Topi, Pakistan, 2007.
9 A. Hebert, A Raviart-Thomas-Schneider solution of the diffusion equation in hexagonal geometry, Ann. Nucl. Energy 35 (2008) 363-376.   DOI
10 S. Cavdar, H. Ozgener, A finite element/boundary element hybrid method for 2-D neutron diffusion calculations, Ann. Nucl. Energy 31 (2004) 1555-1582.   DOI
11 Y. Wang, W. Bangerth, J. Ragusa, Three-dimensional h-adaptivity for the multigroup neutron diffusion equations, Prog. Nucl. Energy 51 (2009) 543-555.   DOI
12 S.A. Hosseini, N. Vosoughi, Development of two-dimensional, multigroup neutron diffusion computer code based on GFEM with unstructured triangle elements, Ann. Nucl. Energy 51 (2013) 213-226.   DOI
13 M. Maiani, B. Montagnini, A Galerkin approach to the boundary element-response matrix method for the multigroup neutron diffusion equations, Ann. Nucl. Energy 31 (2004) 1447-1475.   DOI
14 J. Zhu, Z. Taylor, O. Zienkiewicz, The Finite Element Method: Its Basis and Fundamentals, Butterworth-Heinemann, Barcelona, Spain, 2005.
15 A.C. Center, Benchmark Problem Book, Report ANL-7416 (Suppl. 2), Argonne National Laboratory, Argonne, IL, 1977.
16 G. Schulz, Solutions of a 3D VVER-1000 Benchmark, in: Proceedings of the 6th Symposium of AER on VVER Reactor Physics and Safety, Kirkkonummi, Finland, 1996.
17 J. Wood, C. De Oliveira, A multigroup finite-element solution of the neutron transport equation-I: XY geometry, Ann. Nucl. Energy 11 (1984) 229-243.   DOI
18 J.J. Duderstadt, L.J. Hamilton, Nuclear Reactor Analysis, John Wiley & Sons, Inc, 1976.
19 J.R. Lamarsh, Introduction to Nuclear Reactor Theory, Addison-Wesley, Reading, MA, 1966.
20 G.L. Bell, S. Glasstone, Nuclear Reactor Theory, Van Nostrand Reinhold, New York, 1970.
21 T. Fowler, D. Vondy, G. Cunningham, Nuclear reactor core analysis code; CITATION, ORNL-TM-2496, Rev. 2, Oak Ridge National Laboratory, Oak Ridge, TN, July 1971.
22 D. Vondy, VENTURE-A Code Block for Solving Multigroup Neutronic Problems Applying Finite Difference Diffusion Approximation to Neutron Transport, 1977. ORNL-5062R1.
23 J. Lautard, S. Loubiere, C. Fedon-Magnaud, Cronos: A Modular Computational System for Neutronic Core Calculations, 1992.
24 S. Gonzalez-Pintor, D. Ginestar, G. Verdu, High order finite element method for the lambda modes problem on hexagonal geometry, Ann. Nucl. Energy 36 (2009) 1450-1462.   DOI