Nuclear Engineering and Technology

Korean Nuclear Society (한국원자력학회)
권호 표지
DOI QR코드

DOI QR Code

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

  • Received : 2015.02.28
  • Accepted : 2015.10.19
  • Published : 2016.02.25
Download PDF
⟨ Previous Next ⟩

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

References

  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. https://doi.org/10.1098/rsta.1911.0009
  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. https://doi.org/10.1016/0010-4655(94)90126-0
  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. https://doi.org/10.1016/S0306-4549(01)00083-4
  4. A. Hebert, Development of the nodal collocation method for solving the neutron diffusion equation, Ann. Nucl. Energy 14 (1987) 527-541. https://doi.org/10.1016/0306-4549(87)90074-0
  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. https://doi.org/10.1016/j.anucene.2007.07.016
  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. https://doi.org/10.1016/j.anucene.2004.04.006
  11. Y. Wang, W. Bangerth, J. Ragusa, Three-dimensional h-adaptivity for the multigroup neutron diffusion equations, Prog. Nucl. Energy 51 (2009) 543-555. https://doi.org/10.1016/j.pnucene.2008.11.005
  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. https://doi.org/10.1016/j.anucene.2012.07.032
  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. https://doi.org/10.1016/j.anucene.2004.04.003
  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. https://doi.org/10.1016/0306-4549(84)90054-9
  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. D. Vondy, VENTURE-A Code Block for Solving Multigroup Neutronic Problems Applying Finite Difference Diffusion Approximation to Neutron Transport, 1977. ORNL-5062R1.
  22. J. Lautard, S. Loubiere, C. Fedon-Magnaud, Cronos: A Modular Computational System for Neutronic Core Calculations, 1992.
  23. 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. https://doi.org/10.1016/j.anucene.2009.07.003
  24. 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.

Cited by

  1. Response Surface Methodology Control Rod Position Optimization of a Pressurized Water Reactor Core Considering Both High Safety and Low Energy Dissipation vol.19, pp.2, 2016, https://doi.org/10.3390/e19020063
  2. Extraction of theoretical equation for the gamma ray buildup factor of the three-layered spherical shield vol.14, pp.4, 2016, https://doi.org/10.1088/1748-0221/14/04/p04011
  3. High accurate three-dimensional neutron noise simulator based on GFEM with unstructured hexahedral elements vol.51, pp.6, 2016, https://doi.org/10.1016/j.net.2019.04.006
  4. Multiblock Adaptive Mesh Refinement for the SN Transport Equation Based on Lattice Boltzmann Method vol.193, pp.11, 2019, https://doi.org/10.1080/00295639.2019.1620052
  5. Comparative study between the developed simulators based on the direct and indirect methods for calculation of the energy spectrum of the neutron source vol.14, pp.11, 2019, https://doi.org/10.1088/1748-0221/14/11/p11019
  6. Calculation of effective point reactor kinetic parameters combined with the Galerkin finite element method vol.31, pp.7, 2016, https://doi.org/10.1007/s41365-020-00781-y
  7. A Comparative Analysis of Neutron Transport Calculations Based on Variational Formulation and Finite Element Approaches vol.13, pp.20, 2016, https://doi.org/10.3390/en13205424
  8. Solution of direct and inverse conduction heat transfer problems using the method of fundamental solutions and differential evolution vol.37, pp.9, 2016, https://doi.org/10.1108/ec-01-2020-0017