Calculation of effective point reactor kinetic parameters combined with the Galerkin finite element method

Calculation of effective point reactor kinetic parameters combined with the Galerkin finite element method 1 2 corresp first-author Corresponding author, bxyuan@caep.cn Corresponding author, bxyuan@caep.cn bxyuan@caep.cn https://orcid.org/0000-0002-4134-5799 1 1 2 2 2 China Academy of Engineering Physics, Institute of Nuclear Physics and Chemistry, Mianyang 621900, China China Academy of Engineering Physics, Institute of Nuclear Physics and Chemistry, Mianyang 621900, China Institute of Applied Physics and Computational Mathematics, Beijing 100094, China Institute of Applied Physics and Computational Mathematics, Beijing 100094, China Precise point reactor kinetic parameters are essential in the study of reactor dynamics. Point reactor kinetic parameters include the effective delayed neutron fraction and prompt neutron lifetime. In this work, effective point reactor kinetic parameters, which can be applied to unstructured grids, were calculated based on the Galerkin finite element method. First, two-dimensional and three-dimensional benchmarks were used to verify the calculation of steady-state neutronics parameters. Then, the Tehran research reactor core was divided into hexahedral meshes, and the forward flux and adjoint flux were calculated. Finally, the effective point reactor kinetic parameters of the Tehran research reactor were obtained by comprehensively processing the steady-state neutronics parameters. Key words: Finite element method Reactor Kinetic parameters 1 Introduction 1 The finite element method can handle unstructured meshes and can thus be adapted to complex geometries. Since the 1970s, this method has been gradually applied and promoted in the field of neutronics computing [ 1 1 - 4 4 1-4 ]. Zhang et al. [ 5 5 ] proposed a strategy to accelerate solving the neutron diffusion equation via the finite element method, and the acceleration ratio reached 10 3 . Hosseini et al. [ 6 6 - 9 9 6-9 ] proposed a finite element method for solving steady-state neutronics and neutron noise problems. Avvakumov et al. [ 10 10 ] also solved a typical neutron dynamics problem via the finite element method. The finite element platform FEniCS was used to develop the neutron dynamics code, and the extensible tool package SLEPc was used to solve the eigenvalue of the coefficient matrix. González-Pintor et al. [ 11 11 ] proposed a high-order finite element method to solve the lambda eigenvalue problem of hexagonal geometry. For analyzing reactor dynamics, effective point reactor kinetic parameters need to be obtained in advance, and finite element methods can be used. Saadi et al.[ 12 12 ] reviewed the calculation methods of effective point reactor kinetic parameters and compared deterministic and stochastic methods by using the same benchmark. Arkani et al. [ 13 13 ] calculated the effective point reactor kinetic parameters of the Tehran research reactor (TRR) by using the six-group weighting method. Lashkari et al. [ 14 14 ] calculated the effective point reactor kinetic parameters of the TRR mixed core through the core analysis tool MTR_PC. The effective delayed neutron fraction was found to decrease with burnup, while the prompt neutron lifetime increased. With increasing highly enriched uranium fuel assemblies, the prompt neutron lifetime increases, while the effective delayed neutron fraction does not change significantly. Hosseini et al. [ 15 15 ] calculated and measured the effective point reactor kinetic parameters of the TRR. MTR_PC was used to calculate the neutronics parameters of the TRR in cold and hot states. By using WIMS and CITATION, Zaker et al. [ 16 16 ] calculated the effective point reactor kinetic parameters of two TRR cores with both highly enriched uranium and low enriched uranium. To improve the steady-state calculation in the current deterministic method of effective point reactor kinetic parameters, the calculation of these parameters and the Galerkin finite element method are combined in this work based on the existing framework of FEMN [ 17 17 ]. This paper is organized as follows: the basic concepts of the Galerkin finite element method are outlined, and the 2D and 3D steady-state benchmarks based on the Galerkin finite element method are verified in Sect. 2. In Sect. 3, the TRR benchmarks used worldwide are described, and the effective point reactor kinetic parameters of the TRR are obtained by the improved strategy proposed above. In the final section, the calculation results are discussed, and a brief summary is provided. 2 Galerkin Finite Element Method 1 2.1 Galerkin finite element formulation of the neutron diffusion problem 2 The effective delayed neutron fraction is defined as follows: <math display="block" id="M1"><mrow><msub><mi>β</mi><mi>i</mi></msub><msup><mrow/><mrow><mtext>eff</mtext></mrow></msup><mo>=</mo><mfrac><mrow><mo>∫</mo><mo>∫</mo><mo>∫</mo><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>g</mi><mo>=</mo><mn>1</mn></mrow><mi>G</mi></msubsup><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup><msub><mi>χ</mi><mrow><mi>d</mi><mi>i</mi></mrow></msub><msup><mrow/><mi>q</mi></msup><msub><mi>β</mi><mi>i</mi></msub><msup><mrow/><mi>q</mi></msup><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></msubsup><msub><mrow><mrow><mo>(</mo><mrow><mtext>ν</mtext><msub><mo>∑</mo><mtext>f</mtext></msub></mrow><mo>)</mo></mrow></mrow><msup><mi>g</mi><mo>′</mo></msup></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mtext>d</mtext><mi>V</mi></mrow><mrow><mo>∫</mo><mo>∫</mo><mo>∫</mo><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>g</mi><mo>=</mo><mn>1</mn></mrow><mi>G</mi></msubsup><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup><msub><mi>χ</mi><mi>p</mi></msub><msup><mrow/><mi>q</mi></msup><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></msubsup><mi>ν</mi><msub><mtext>Σ</mtext><mtext>f</mtext></msub><msup><mrow/><msup><mi>g</mi><mo>′</mo></msup></msup><mrow><mo>(</mo><mtext>r</mtext><mo>)</mo></mrow><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mtext>d</mtext><mi>V</mi></mrow></mfrac><mo>.</mo></mrow></math> For group g , <math id="M2"><mrow><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup></mrow></math> is adjoint flux, <math id="M3"><mrow><msub><mi>χ</mi><mrow><mtext>d</mtext><mi>i</mi></mrow></msub><msup><mrow/><mi>q</mi></msup></mrow></math> is the delayed neutron spectrum of fissionable isotope <math id="M4"><mi>q</mi></math> in group i , <math id="M5"><mrow><msub><mi>β</mi><mi>i</mi></msub><msup><mrow/><mi>q</mi></msup></mrow></math> is the delayed neutron yield for fissionable isotope <math id="M6"><mi>q</mi></math> in group i , <math id="M7"><mtext>ν</mtext></math> is the number of neutrons emitted per fission, <math id="M8"><mrow><msub><mrow><mrow><mo>(</mo><mrow><msub><mo>∑</mo><mtext>f</mtext></msub></mrow><mo>)</mo></mrow></mrow><msup><mi>g</mi><mo>′</mo></msup></msub></mrow></math> is the fission cross section, <math id="M9"><mrow><msub><mo>∅</mo><mi>g</mi></msub></mrow></math> is forward flux, and <math id="M10"><mrow><msub><mi>χ</mi><mi>p</mi></msub><msup><mrow/><mi>q</mi></msup></mrow></math> is the prompt neutron spectrum of fissionable isotope <math id="M11"><mi>q</mi></math> . Similarly, the definition of prompt neutron lifetime can be given as follows: <math display="block" id="M12"><mrow><mi>l</mi><mo>=</mo><mfrac><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub><mo>∫</mo><mo>∫</mo><mo>∫</mo><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>g</mi><mo>=</mo><mn>1</mn></mrow><mi>G</mi></msubsup><msub><mo>∅</mo><mtext>g</mtext></msub><msup><mrow/><mtext>*</mtext></msup><mfrac><mn>1</mn><mrow><msub><mi>v</mi><mi>g</mi></msub></mrow></mfrac><msub><mo>∅</mo><mtext>g</mtext></msub><mtext>d</mtext><mi>V</mi></mrow><mrow><mo>∫</mo><mo>∫</mo><mo>∫</mo><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>g</mi><mo>=</mo><mn>1</mn></mrow><mi>G</mi></msubsup><msub><mo>∅</mo><mtext>g</mtext></msub><msup><mrow/><mtext>*</mtext></msup><msub><mi>χ</mi><mi>i</mi></msub><msup><mrow/><mi>g</mi></msup><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></msubsup><mi>ν</mi><msub><mtext>Σ</mtext><mtext>f</mtext></msub><msup><mrow/><msup><mi>g</mi><mo>′</mo></msup></msup><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mtext>d</mtext><mi>V</mi></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M13"><mrow><mfrac><mn>1</mn><mrow><msub><mi>v</mi><mi>g</mi></msub></mrow></mfrac></mrow></math> is the inverse velocity of group <math id="M14"><mrow><mtext> </mtext><mi>g</mi></mrow></math> . The 3D steady-state equation of neutronics is as follows: <math display="block" id="M15"><mrow><mo>−</mo><mo>∇</mo><mo>⋅</mo><msub><mi>D</mi><mi>g</mi></msub><mo>∇</mo><msub><mo>∅</mo><mi>g</mi></msub><mo>+</mo><msub><mo>∑</mo><mrow><mi>r</mi><mo>,</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><mi>g</mi></msub><mo>=</mo><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></munderover><msub><mo>∑</mo><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>→</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mo>+</mo><mfrac><mrow><msub><mtext>χ</mtext><mi>g</mi></msub></mrow><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></munderover><msub><mrow><mrow><mo>(</mo><mrow><mtext>ν</mtext><msub><mo>∑</mo><mtext>f</mtext></msub></mrow><mo>)</mo></mrow></mrow><msup><mi>g</mi><mo>′</mo></msup></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mo>.</mo></mrow></math> For group g , <math id="M16"><mrow><msub><mi>D</mi><mtext>g</mtext></msub></mrow></math> is the diffusion coefficient, <math id="M17"><mrow><msub><mo>∑</mo><mrow><mi>r</mi><mo>,</mo><mi>g</mi></mrow></msub></mrow></math> is the removal cross section, <math id="M18"><mrow><msub><mo>∑</mo><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>→</mo><mi>g</mi></mrow></msub></mrow></math> is the scattering cross section from group <math id="M19"><msup><mi>g</mi><mo>′</mo></msup></math> to <math id="M20"><mi>g</mi></math> , <math id="M21"><mrow><msub><mtext>χ</mtext><mi>g</mi></msub></mrow></math> is the fission spectrum, and <math id="M22"><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></math> is the effective multiplication factor. The boundary conditions are described as follows: <math display="block" id="M23"><mrow><mrow><mo>(</mo><mrow><mtext>α</mtext><msub><mi>D</mi><mi>g</mi></msub><mfrac><mrow><mo>∂</mo><msub><mo>∅</mo><mi>g</mi></msub></mrow><mrow><mo>∂</mo><mi>n</mi></mrow></mfrac><mo>+</mo><mtext>β</mtext><msub><mo>∅</mo><mi>g</mi></msub></mrow><mo>)</mo></mrow><mo>|</mo><msub><mtext>Г</mtext><mn>1</mn></msub><mo>　</mo><mo>=</mo><mn>0</mn><mo>　</mo><mo>　</mo><mtext>α</mtext><mo>≠</mo><mtext>0</mtext><mo>,</mo></mrow></math> where <math id="M24"><mrow><msub><mtext>Г</mtext><mn>1</mn></msub></mrow></math> is the boundary, n is the normal direction of <math id="M25"><mrow><msub><mtext>Г</mtext><mn>1</mn></msub></mrow></math> , and <math id="M26"><mtext>α</mtext></math> and <math id="M27"><mtext>β</mtext></math> are constants. By following the Galerkin finite element method, the weak form of the diffusion equation can be written as follows: <math display="block" id="M28"><mrow><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mo>−</mo><mi>v</mi><mo>∇</mo><mo>⋅</mo><msub><mi>D</mi><mi>g</mi></msub><mo>∇</mo><msub><mo>∅</mo><mi>g</mi></msub><mtext>d</mtext><mi>V</mi><mo>+</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><msub><mo>∑</mo><mrow><mtext>r</mtext><mo>,</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><mi>g</mi></msub><mtext>d</mtext><mi>V</mi><mo>=</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></munderover><msub><mo>∑</mo><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>→</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mtext>d</mtext><mi>V</mi><mo>+</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><mfrac><mrow><msub><mtext>χ</mtext><mi>g</mi></msub></mrow><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mtext>G</mtext></munderover><msub><mrow><mrow><mo>(</mo><mrow><mtext>ν</mtext><msub><mo>∑</mo><mtext>f</mtext></msub></mrow><mo>)</mo></mrow></mrow><msup><mi>g</mi><mo>′</mo></msup></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mtext>d</mtext><mi>V</mi><mo>.</mo></mrow></math> After integration by parts, the finite element form of the 3D steady-state neutronics equation can be obtained: <math display="block" id="M29"><mrow><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><msub><mi>D</mi><mi>g</mi></msub><mo>∇</mo><msub><mo>∅</mo><mi>g</mi></msub><mo>⋅</mo><mo>∇</mo><mi>v</mi><mtext>d</mtext><mi>V</mi><mo>+</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><msub><mo>∑</mo><mrow><mtext>r</mtext><mo>,</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><mi>g</mi></msub><mtext>d</mtext><mi>V</mi><mo>+</mo><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∬</mo></mstyle><mrow><msub><mtext>Г</mtext><mn>1</mn></msub></mrow><mtext> </mtext></munderover><mfrac><mtext>β</mtext><mtext>α</mtext></mfrac><msub><mo>∅</mo><mi>g</mi></msub><mi>v</mi><mtext>d</mtext><mi>S</mi><mo>=</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></munderover><msub><mo>∑</mo><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>→</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mtext>d</mtext><mi>V</mi><mo>+</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><mfrac><mrow><msub><mtext>χ</mtext><mi>g</mi></msub></mrow><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></munderover><msub><mrow><mrow><mo>(</mo><mrow><mtext>ν</mtext><msub><mo>∑</mo><mtext>f</mtext></msub></mrow><mo>)</mo></mrow></mrow><msup><mi>g</mi><mo>′</mo></msup></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><mtext>d</mtext><mi>V</mi><mo>.</mo></mrow></math> The 3D steady-state adjoint equation of neutronics is as follows: <math display="block" id="M30"><mrow><mo>−</mo><mo>∇</mo><mo>⋅</mo><msub><mi>D</mi><mi>g</mi></msub><mo>∇</mo><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup><mo>+</mo><msub><mo>∑</mo><mrow><mtext>r</mtext><mo>,</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup><mo>=</mo><munder><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>≠</mo><mi>g</mi></mrow></munder><msub><mo>∑</mo><mrow><mi>g</mi><mo>→</mo><msup><mi>g</mi><mo>′</mo></msup></mrow></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><msup><mrow/><mtext>*</mtext></msup><mo>+</mo><mfrac><mrow><msub><mrow><mrow><mo>(</mo><mrow><mtext>ν</mtext><msub><mo>∑</mo><mtext>f</mtext></msub></mrow><mo>)</mo></mrow></mrow><mi>g</mi></msub></mrow><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></munderover><msub><mtext>χ</mtext><msup><mi>g</mi><mo>′</mo></msup></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><msup><mrow/><mtext>*</mtext></msup><mo>.</mo></mrow></math> Similarly, the finite element form of the 3D steady-state neutronics adjoint equation can be obtained as follows: <math display="block" id="M31"><mrow><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><msub><mi>D</mi><mi>g</mi></msub><mo>∇</mo><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup><mo>⋅</mo><mo>∇</mo><mi>v</mi><mtext>d</mtext><mi>V</mi><mo>+</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><msub><mo>∑</mo><mrow><mtext>r</mtext><mo>,</mo><mi>g</mi></mrow></msub><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup><mtext>d</mtext><mi>V</mi><mo>+</mo><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∬</mo></mstyle><mrow><msub><mtext>Г</mtext><mn>1</mn></msub></mrow><mtext> </mtext></munderover><mfrac><mtext>β</mtext><mtext>α</mtext></mfrac><msub><mo>∅</mo><mi>g</mi></msub><msup><mrow/><mtext>*</mtext></msup><mi>v</mi><mtext>d</mtext><mi>S</mi><mo>=</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>≠</mo><mi>g</mi></mrow><mi>G</mi></munderover><msub><mo>∑</mo><mrow><mi>g</mi><mo>→</mo><msup><mi>g</mi><mo>′</mo></msup></mrow></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><msup><mrow/><mtext>*</mtext></msup><mtext>d</mtext><mi>V</mi><mo>+</mo><mo>∫</mo><mo>∫</mo><msub><mo>∫</mo><mi>Ω</mi></msub><mi>v</mi><mfrac><mrow><msub><mrow><mrow><mo>(</mo><mrow><mtext>ν</mtext><msub><mo>∑</mo><mtext>f</mtext></msub></mrow><mo>)</mo></mrow></mrow><mi>g</mi></msub></mrow><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>G</mi></munderover><msub><mtext>χ</mtext><msup><mi>g</mi><mo>′</mo></msup></msub><msub><mo>∅</mo><msup><mi>g</mi><mo>′</mo></msup></msub><msup><mrow/><mtext>*</mtext></msup><mtext>d</mtext><mi>V</mi><mo>.</mo></mrow></math> 2.2 Benchmark test 2 1) 2D steady-state neutronics benchmark The 2D benchmark TWIGL is used for the code test [ 18 18 ]. The geometric layout of the 1/4 core is shown in Fig. 1 Fig. 1 , and the macroscopic cross sections used are listed in Table 1 Table 1 . Both triangular and quadrilateral elements can be used for 2D geometry meshing, and Fig. 2 Fig. 2 shows the 2D mesh adopted by the present study. The steady-state neutronics parameters are calculated by FEMN. <math id="M32"><mrow><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mo>*</mo></msubsup></mrow></math> and adjoint flux are shown in Table 2 Table 2 and 3 3 and 4 4 , respectively. 1 1 Material types 1 1 Energy group 1 1 <math id="M33"><mrow><msub><mi>D</mi><mi>g</mi></msub><mo> </mo><mrow><mo>(</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>c</mi></mstyle><msup><mstyle mathvariant="bold" mathsize="normal"><mi>m</mi></mstyle><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 1 <math id="M34"><mrow><msub><mo>∑</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>t</mi></mstyle><mo>,</mo><mi>g</mi></mrow></msub><mrow><mo>(</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>c</mi></mstyle><msup><mstyle mathvariant="bold" mathsize="normal"><mi>m</mi></mstyle><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 1 <math id="M35"><mrow><msub><mrow><mrow><mo>(</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>ν</mi></mstyle><msub><mo>∑</mo><mstyle mathvariant="bold" mathsize="normal"><mi>f</mi></mstyle></msub></mrow><mo>)</mo></mrow></mrow><mi>g</mi></msub><mo> </mo><mrow><mo>(</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>c</mi></mstyle><msup><mstyle mathvariant="bold" mathsize="normal"><mi>m</mi></mstyle><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 1 <math id="M36"><mrow><msub><mo>∑</mo><mrow><msup><mi>g</mi><mo>′</mo></msup><mo>→</mo><mi>g</mi></mrow></msub><mrow><mo>(</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>c</mi></mstyle><msup><mstyle mathvariant="bold" mathsize="normal"><mi>m</mi></mstyle><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 1 1 1 1 1 1 1 1.4 1 1 0.01 1 1 0.007657 1 1 0.01 1 1 1 1 2 1 1 0.4 1 1 0.15 1 1 0.218772 1 1 0 1 1 2 1 1 1 1 1 1.3 1 1 0.008 1 1 0.00328158 1 1 0.01 1 1 1 1 2 1 1 0.5 1 1 0.05 1 1 0.06563 1 1 0 1 1 Item 1 1 Calculated value/reference value (deviation) 1 1 <math id="M38"><mrow><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mo>*</mo></msubsup></mrow></math>   1 1 1.001312/1.001310 (0.000002) —reference value —calculated value —deviation % 1 1 0.2138 0.2103 1.62% 1 1 0.2095 0.2060 1.65% 1 1 0.2006 0.1972 1.68% 1 1 0.1866 0.1835 1.68% 1 1 0.1674 0.1645 1.75% 1 1 0.1431 0.1405 1.85% 1 1 0.1164 0.1125 3.39% 1 1 0.0839 0.0823 1.96% 1 1 0.0553 0.0520 5.96% 1 1 0.0238 0.0229 3.92% 1 1 0.4910 0.4892 0.37% 1 1 0.4814 0.4796 0.37% 1 1 0.4615 0.4596 0.41% 1 1 0.4301 0.4282 0.44% 1 1 0.3862 0.3843 0.48% 1 1 0.3298 0.3281 0.51% 1 1 0.2631 0.2614 0.63% 1 1 0.1909 0.1895 0.71% 1 1 0.1200 0.1186 1.14% 1 1 0.0533 0.0520 2.43% 1 1 0.8132 0.8129 0.04% 1 1 0.7985 0.7981 0.05% 1 1 0.7674 0.7669 0.07% 1 1 0.7173 0.7165 0.11% 1 1 0.6455 0.6446 0.14% 1 1 0.5513 0.5503 0.18% 1 1 0.4365 0.4354 0.26% 1 1 0.3099 0.3088 0.36% 1 1 0.1909 0.1895 0.75% 1 1 0.0839 0.0822 2.02% 1 1 1.1600 1.1597 0.02% 1 1 1.1415 1.1412 0.03% 1 1 1.1014 1.1009 0.05% 1 1 1.0344 1.0337 0.07% 1 1 0.9350 0.9341 0.10% 1 1 0.8001 0.7990 0.14% 1 1 0.6305 0.6293 0.19% 1 1 0.4363 0.4355 0.19% 1 1 0.2629 0.2615 0.55% 1 1 0.1145 0.1125 1.78% 1 1 1.4239 1.4243 0.03% 1 1 1.4063 1.4067 0.03% 1 1 1.3660 1.3663 0.02% 1 1 1.2935 1.2936 0.01% 1 1 1.1782 1.1780 0.02% 1 1 1.0136 1.0131 0.05% 1 1 0.7998 0.7990 0.10% 1 1 0.5509 0.5503 0.11% 1 1 0.3295 0.3280 0.44% 1 1 0.1429 0.1405 1.68% 1 1 1.5655 1.5667 0.08% 1 1 1.5556 1.5569 0.09% 1 1 1.5289 1.5302 0.09% 1 1 1.4693 1.4705 0.08% 1 1 1.3563 1.3571 0.06% 1 1 1.1778 1.1781 0.02% 1 1 0.9344 0.9341 0.03% 1 1 0.6449 0.6447 0.04% 1 1 0.3857 0.3843 0.36% 1 1 0.1672 0.1645 1.62% 1 1 1.5634 1.5654 0.12% 1 1 1.5683 1.5705 0.14% 1 1 1.5739 1.5762 0.15% 1 1 1.5566 1.5588 0.14% 1 1 1.4687 1.4705 0.12% 1 1 1.2925 1.2936 0.08% 1 1 1.0334 1.0337 0.03% 1 1 0.7164 0.7165 0.01% 1 1 0.4295 0.4282 0.31% 1 1 0.1864 0.1834 1.59% 1 1 1.4381 1.4419 0.27% 1 1 1.4617 1.4660 0.29% 1 1 1.5145 1.5191 0.31% 1 1 1.5731 1.5762 0.19% 1 1 1.5276 1.5302 0.17% 1 1 1.3645 1.3663 0.13% 1 1 1.1000 1.1009 0.08% 1 1 0.7664 0.7669 0.06% 1 1 0.4609 0.4596 0.27% 1 1 0.2003 0.1972 1.54% 1 1 1.3344 1.3388 0.29% 1 1 1.3734 1.3778 0.32% 1 1 1.4610 1.4660 0.34% 1 1 1.5668 1.5706 0.24% 1 1 1.5536 1.5569 0.21% 1 1 1.4043 1.4067 0.17% 1 1 1.1398 1.1411 0.12% 1 1 0.7973 0.7981 0.10% 1 1 0.4807 0.4796 0.23% 1 1 0.2092 0.2061 1.50% 1 1 1.2871 1.2910 0.30% 1 1 1.3339 1.3382 0.33% 1 1 1.4367 1.4419 0.36% 1 1 1.5613 1.5653 0.26% 1 1 1.5632 1.5668 0.23% 1 1 1.4216 1.4243 0.19% 1 1 1.1582 1.1597 0.13% 1 1 0.8120 0.8129 0.11% 1 1 0.4903 0.4892 0.22% 1 1 0.2135 0.2103 1.18% —reference value —calculated value —deviation % 1 1 0.2656 0.2612 1.64% 1 1 0.2603 0.2559 1.70% 1 1 0.2492 0.2449 1.72% 1 1 0.2319 0.2278 1.75% 1 1 0.2080 0.2042 1.81% 1 1 0.1777 0.1744 1.83% 1 1 0.1424 0.1397 1.92% 1 1 0.1042 0.1021 2.02% 1 1 0.0660 0.0644 2.40% 1 1 0.0277 0.0266 3.97% 1 1 0.6510 0.6486 0.37% 1 1 0.6384 0.6359 0.40% 1 1 0.6120 0.6094 0.43% 1 1 0.5703 0.5677 0.45% 1 1 0.5120 0.5096 0.48% 1 1 0.4373 0.4350 0.53% 1 1 0.3487 0.3465 0.64% 1 1 0.2528 0.2509 0.74% 1 1 0.1584 0.1567 1.10% 1 1 0.0660 0.0644 2.39% 1 1 1.0984 1.0933 0.08% 1 1 1.0785 1.0791 0.06% 1 1 1.0365 1.0371 0.06% 1 1 0.9687 0.9689 0.02% 1 1 0.8718 0.8716 0.03% 1 1 0.7445 0.7440 0.07% 1 1 0.5880 0.5872 0.14% 1 1 0.4130 0.4115 0.37% 1 1 0.2528 0.2508 0.78% 1 1 0.1042 0.1020 2.09% 1 1 1.6698 1.6683 0.09% 1 1 1.6432 1.6416 0.10% 1 1 1.5855 1.5836 0.12% 1 1 1.4891 1.4870 0.14% 1 1 1.3460 1.3437 0.17% 1 1 1.1517 1.1494 0.20% 1 1 0.9013 0.8985 0.31% 1 1 0.5878 0.5874 0.07% 1 1 0.3485 0.3465 0.58% 1 1 0.1422 0.1396 1.79% 1 1 2.0651 2.0650 0.00% 1 1 2.0395 2.0395 0.00% 1 1 1.9811 1.9809 0.01% 1 1 1.8760 1.8755 0.02% 1 1 1.7089 1.7080 0.05% 1 1 1.4700 1.4688 0.08% 1 1 1.1514 1.1494 0.18% 1 1 0.7440 0.7440 0.00% 1 1 0.4369 0.4349 0.45% 1 1 0.1775 0.1745 1.69% 1 1 2.2705 2.2715 0.04% 1 1 2.2561 2.2574 0.06% 1 1 2.2174 2.2186 0.06% 1 1 2.1311 2.1321 0.05% 1 1 1.9672 1.9677 0.03% 1 1 1.7082 1.7081 0.01% 1 1 1.3451 1.3438 0.01% 1 1 0.8710 0.8716 0.07% 1 1 0.5115 0.5095 0.39% 1 1 0.2077 0.2043 1.63% 1 1 2.2515 2.2527 0.05% 1 1 2.2588 2.2604 0.07% 1 1 2.2698 2.2716 0.08% 1 1 2.2559 2.2586 0.12% 1 1 2.1302 2.1322 0.09% 1 1 1.8746 1.8756 0.05% 1 1 1.4876 1.4870 0.04% 1 1 0.9676 0.9688 0.12% 1 1 0.5695 0.5677 0.32% 1 1 0.2316 0.2278 1.63% 1 1 1.9432 1.9501 0.36% 1 1 1.9767 1.9840 0.37% 1 1 2.0742 2.0840 0.47% 1 1 2.2687 2.2716 0.13% 1 1 2.2154 2.2186 0.15% 1 1 1.9789 1.9809 0.10% 1 1 1.5835 1.5837 0.01% 1 1 1.0351 1.0370 0.18% 1 1 0.6111 0.6094 0.28% 1 1 0.2488 0.2449 1.57% 1 1 1.7717 1.7762 0.26% 1 1 1.8258 1.8310 0.28% 1 1 1.9756 1.9841 0.43% 1 1 2.2566 2.2604 0.17% 1 1 2.2532 2.2573 0.18% 1 1 2.0366 2.0396 0.14% 1 1 1.6408 1.6415 0.04% 1 1 1.0769 1.0793 0.23% 1 1 0.6374 0.6359 0.24% 1 1 0.2599 0.2559 1.55% 1 1 1.7066 1.7112 0.27% 1 1 1.7710 1.7762 0.29% 1 1 1.9414 1.9500 0.44% 1 1 2.2484 2.2527 0.19% 1 1 2.2670 2.2716 0.20% 1 1 2.0617 2.0651 0.17% 1 1 1.6672 1.6682 0.06% 1 1 1.0967 1.0993 0.24% 1 1 0.6501 0.6486 0.23% 1 1 0.2652 0.2612 1.51% Compared with the reference value [ 18 18 ], the deviation of <math id="M39"><mrow><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mtext>*</mtext></msubsup></mrow></math> is 0.000002, while the deviation of adjoint flux in the outermost assemblies is less than 4%. The deviation is mainly affected by the mesh size, and it can be further reduced by mesh refinement. 2) 3D steady -state neutronics benchmark The 3D benchmark IAEA-3D PWR is used for the code test [ 19 19 ]. A schematic of the 3D benchmark is shown in Fig. 3 Fig. 3 , and the macroscopic cross sections used are listed in Table 5 Table 5 . Both triangular prism and hexahedral elements can be used for 3D geometric meshing. Figure 4 Figure 4 shows the 3D mesh adopted by the present study. The steady-state neutronics parameters are calculated by FEMN. k eff , <math id="M40"><mrow><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mo>*</mo></msubsup></mrow></math> , and power distribution are shown in Table 6 Table 6 and 7 7 . 2 1 Material types 1 1 Energy group 1 1 <math id="M41"><mrow><msub><mi>D</mi><mi>g</mi></msub><mo> </mo><mrow><mo>(</mo><mrow><msup><mrow><mtext>cm</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 1 <math id="M42"><mrow><msub><mo>∑</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>t</mi></mstyle><mo>,</mo><mi>g</mi></mrow></msub><mo> </mo><mrow><mo>(</mo><mrow><msup><mrow><mtext>cm</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 1 <math id="M43"><mrow><msub><mrow><mrow><mo>(</mo><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>ν</mi></mstyle><msub><mo>∑</mo><mstyle mathvariant="bold" mathsize="normal"><mi>f</mi></mstyle></msub></mrow><mo>)</mo></mrow></mrow><mi>g</mi></msub><mrow><mo>(</mo><mrow><msup><mrow><mtext>cm</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 1 <math id="M44"><mrow><msub><mo>∑</mo><mrow><mi>g</mi><mo>→</mo><msup><mi>g</mi><mo>'</mo></msup></mrow></msub><mrow><mo>(</mo><mrow><msup><mrow><mtext>cm</mtext></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>)</mo></mrow></mrow></math>   1 2 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310009&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310018&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310013&type= 1 1 Fuel1 1 1 1 1 1 1.500 1 1 0.010 1 1 0.000 1 1 0.020 1 1 1 1 2 1 1 0.400 1 1 0.085 1 1 0.135 1 1 — 1 2 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310022&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310030&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310025&type= 1 1 Fuel1+rod 1 1 1 1 1 1.500 1 1 0.010 1 1 0.000 1 1 0.020 1 1 1 1 2 1 1 0.400 1 1 0.130 1 1 0.135 1 1 — 1 2 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310033&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310041&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310037&type= 1 1 Fuel2 1 1 1 1 1 1.500 1 1 0.010 1 1 0.000 1 1 0.020 1 1 1 1 2 1 1 0.400 1 1 0.080 1 1 0.135 1 1 — 1 2 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310046&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310054&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310050&type= 1 1 Reflector 1 1 1 1 1 2.000 1 1 0.000 1 1 0.000 1 1 0.040 1 1 1 1 2 1 1 0.300 1 1 0.010 1 1 0.000 1 1 — 1 2 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310058&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310068&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10310063&type= 1 1 Reflector+Rod 1 1 1 1 1 2.000 1 1 0.000 1 1 0.000 1 1 0.040 1 1 1 1 2 1 1 0.300 1 1 0.055 1 1 0.000 1 1 — 1 1 Item 1 1 Value 1 1 Reference k eff 1 1 1.02900 1 1 Calculated k eff 1 1 1.02904 1 1 Calculated <math id="M47"><mrow><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mtext>*</mtext></msubsup></mrow></math>   1 1 1.02910 —reference value —calculated value —deviation % 1 1 0.7290 1 1 1 1 1 1 1 1 1 1 1 1 0.7237 1 1 0.73% 1 1 1.2810 1.2643 1.31% 1 1 1.3970 1.3802 1.21% 1 1 1 1 1 1 1 1 1 1 1.4220 1.4068 1.07% 1 1 1.4320 1.4160 1.12% 1 1 1.3680 1.3547 0.98% 1 1 1 1 1 1 1 1 1.1930 1.1815 0.96% 1 1 1.2910 1.2798 0.87% 1 1 1.3110 1.3019 0.70% 1 1 1.1780 1.1725 0.47% 1 1 1 1 1 1 0.6100 0.6113 0.21% 1 1 1.0720 1.0663 0.54% 1 1 1.1810 1.1779 0.26% 1 1 0.9720 0.9701 0.20% 1 1 0.4760 0.4794 0.72% 1 1 1 1 0.9530 0.9536 0.07% 1 1 1.0550 1.0568 0.17% 1 1 1.0890 1.0923 0.31% 1 1 0.9230 0.9275 0.49% 1 1 0.7000 0.7046 0.65% 1 1 0.5970 0.6093 2.06% 1 1 0.9590 0.9664 0.78% 1 1 0.9760 0.9829 0.71% 1 1 1.0000 1.0077 0.77% 1 1 0.8660 0.8789 1.50% 1 1 0.6100 0.6233 2.01% 1 1 1 1 0.77700.78961.62% 1 1 0.75700.76901.58% 1 1 0.71100.72622.14% 1 1 1 1 1 1 Compared with the reference value [ 19 19 ], the deviation of <math id="M48"><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></math> is 0.00004. The deviation of heat deposition of the inner fuel assembly is less than 1%. Although the deviation of the outer fuel assembly close to the reflector is high, the overall deviation is less than 2.2%. The deviation is mainly affected by the mesh size, and it can be further reduced by mesh refinement. The forward flux and adjoint flux can also be adequately determined by FEMN. The distribution of forward flux and adjoint flux for the 3D benchmark are shown in Fig. 5 Fig. 5 and Fig. 6 Fig. 6 , respectively. 3 Results and Discussion 1 3.1 Description of the TRR 2 The TRR is a 5-MW pool-type light-water moderated, heterogeneous solid fuel reactor, in which the water is also used for cooling and shielding. The TRR core is immersed in either section of a two-section, concrete pool filled with water. Utilization of the reactor is essential for research, training, and production of radioisotopes [ 12 12 ]. Substantial research has been performed on the point reactor kinetic parameters of the TRR; thus, it was chosen as a benchmark to facilitate the calculation comparison. The first core configuration of the TRR is shown in Fig. 7 Fig. 7 [ 13 13 ]. The standard fuel element (SFE) consists of 19 fuel plates with 20% enrichment of uranium 235. Each fuel plate consists of two unit areas of fuel cladding and one fuel meat; the meat is made of U3O8 powder dispersed in a pure aluminum matrix, and the cladding and other structural materials are Al6061. The coolant channel is located between the fuel plates. The control fuel element (CFE) consists of 14 fuel plates with 20% enrichment of uranium 235, and each fuel plate consists of two unit areas of fuel cladding and one fuel meat. The shim safety control rod is composed of two absorbing plates and a knuckle subassembly. The absorber plate is an alloy of silver, indium, and cadmium (80, 15, and 5% wt., respectively), while the regulating rod is made of stainless steel. All absorbing rods are fork type. The main geometrical data of the CFE and SFE can be found in reference [ 12 12 ], and the configuration is shown in Fig. 8 Fig. 8 [ 13 13 ]. The calculation of TRR's effective point reactor kinematic parameters involves three steps: 1) Preparation of the six-group cross section and six-group neutron velocity by WIMS and 2DSN; 2) Calculation of the six-group forward flux and adjoint flux with the Galerkin finite element method; 3) Calculation of the effective point reactor kinetic parameters. Because the TRR core is asymmetric, the use of triangular prism elements for geometric meshing and six-group calculation will cause substantial memory consumption, which is impractical for personal computers. Therefore, this work uses a memory-saving hexahedron meshing method, and Fig. 9 Fig. 9 shows the meshes used in the calculation. To distinguish the energy spectrum between prompt and delayed neutrons, the total number of energy groups must be greater than two groups, because in a typical two-group structure, neutrons can only be produced in the fast group, so the effect of delayed neutrons is not reflected. The neutron group structure, transient neutron energy spectrum, delayed neutron energy spectrum, and delayed neutron parameters of 235 U and 238 U can be found in reference [ 12 12 ]. 3.2 Results and discussion 2 The steady-state neutronics parameters of the TRR were calculated using the FEMN, and the excess reactivity in this state is 6344 pcm. The forward flux and adjoint flux in the axial plane of 51.5 cm are shown in Figs. 10 Figs. 10 and 11 11 , respectively. The forward flux and adjoint flux distribution of the TRR in this state is not given in relevant literature, so only the calculation results are given here for peer reference. For FEMN, the relative error of the steady state calculation has been discussed, and interested readers can refer to Sect. 2. From the perspective of excess reactivity, the deviation between the calculated value (6344 pcm) and the reference value (6481 pcm) of reference [ 12 12 ] is 2%, which proves that the steady-state calculation is reliable from the side. Equations (1) and (2) proposed above are used to process the effective point reactor kinetic parameters of the TRR, and the calculated results are shown in Table 8 Table 8 . Compared with the recommended value of reference [ 12 12 ], the relative error of the effective delayed neutron fraction is 2%, while the relative error of the prompt neutron lifetime is 10%. The main reason for the relative error is that the cross sections and neutron energy spectra used by different researchers are inconsistent, which eventually leads to some slight differences in forward flux, adjoint flux, and group neutron velocity. 1 1 Item 1 1 1 1 Value (calculated/recommended [ 12 12 ]) 1 1 Effective delayed neutron fraction 1 1 Group 1 1 1 0.00026/0.00024 1 1 1 1 Group 2 1 1 0.00148/0.00158 1 1 1 1 Group 3 1 1 0.00133/0.00142 1 1 1 1 Group 4 1 1 0.00288/0.00286 1 1 1 1 Group 5 1 1 0.00094/0.00085 1 1 1 1 Group 6 1 1 0.00020/0.00031 1 1 1 1 Total 1 1 0.00709/0.00726 1 1 Prompt neutron lifetime (µ sec) 1 1 1 1 43/48 4 Conclusion 1 1) Based on the internally developed code FEMN, a method for determining effective point reactor kinetic parameters was developed and validated. 2) Compared with the reference values, the relative errors of the effective delayed neutron fraction and prompt neutron lifetime are 2% and 10%, respectively. The main factor that led to the errors was the inconsistency between the cross sections and the neutron energy spectra. References A. Hébert , Application of a dual variational formulation to finite element reactor calculations . Ann. Nucl. Energy . 20 ( 12 ), 823 - 845 ( 1993 ). https://doi.org/10.1016/0306-4549(93)90076-2 https://doi.org/10.1016/0306-4549(93)90076-2 Application of a dual variational formulation to finite element reactor calculations L.A. Semenza , E.E. Lewis , E.C. Rossow , The application of the finite element method to the multigroup neutron diffusion equation . Nucl. Sci. Eng . 47 ( 3 ), 302 - 310 ( 1972 ). https://doi.org/10.13182/NSE72-A22416 https://doi.org/10.13182/NSE72-A22416 The application of the finite element method to the multigroup neutron diffusion equation Y. Naito , R. Kawamura , Y. Tokuno et al ., A computer program for solving the three-dimensional multigroup diffusion equation by finite element method with 20-node isoparametric element (FEM3DJAR) . Prog. Nucl. Energ . 18 ( 1-2 ), 207 - 214 ( 1986 ). https://doi.org/10.1016/0149-1970(86)90027-2 https://doi.org/10.1016/0149-1970(86)90027-2 A computer program for solving the three-dimensional multigroup diffusion equation by finite element method with 20-node isoparametric element (FEM3DJAR) K. AZEKURA , New finite element solution technique for neutron diffusion equations . J. Nucl. Sci. Technol . 17 ( 2 ), 89 - 97 ( 1980 ). https://doi.org/10.1080/18811248.1980.9732552 https://doi.org/10.1080/18811248.1980.9732552 New finite element solution technique for neutron diffusion equations C.Y Zhang , G. Chen , Fast solution of neutron diffusion problem by reduced basis finite element method . Ann. Nucl. Energy . 111 , 702 - 708 ( 2018 ). https://doi.org/10.1016/j.anucene.2017.09.044 https://doi.org/10.1016/j.anucene.2017.09.044 Fast solution of neutron diffusion problem by reduced basis finite element method S.A. Hosseini , N. Vosoughi , Neutron noise simulation by GFEM and unstructured triangle elements . Nucl. Eng. Des . 253 , 238 - 258 ( 2012 ). https://doi.org/10.1016/j.nucengdes.2012.08.023 https://doi.org/10.1016/j.nucengdes.2012.08.023 Neutron noise simulation by GFEM and unstructured triangle elements S.A. Hosseini , F. Saadatian-Derakhshandeh , Galerkin and generalized least squares finite element: a comparative study for multi-group diffusion solvers . Prog. Nucl. Energ . 85 , 473 - 490 ( 2015 ). https://doi.org/10.1016/j.pnucene.2015.07.009 https://doi.org/10.1016/j.pnucene.2015.07.009 Galerkin and generalized least squares finite element: a comparative study for multi-group diffusion solvers S.A. Hosseini , Development of galerkin finite element method three-dimensional computational code for the multigroup neutron diffusion equation with unstructured tetrahedron elements . Nucl. Eng. Technol . 48 , 43 - 54 ( 2016 ). https://doi.org/10.1016/j.net.2015.10.009 https://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 S.A. Hosseini , Sensitivity analysis of the galerkin finite element method neutron diffusion solver to the shape of the elements . Nucl. Eng. Technol . 49 , 29 - 42 ( 2017 ). https://doi.org/10.1016/j.net.2016.08.006 https://doi.org/10.1016/j.net.2016.08.006 Sensitivity analysis of the galerkin finite element method neutron diffusion solver to the shape of the elements A.V. Avvakumov , P.N. Vabishchevich , A.O. Vasilev et al ., Solution of the neutronics code dynamic benchmark by finite element method . AIP Conference Proceedings . 1773 , 110003 , ( 2016 ). https://doi.org/10.1063/1.4965007 https://doi.org/10.1063/1.4965007 Solution of the neutronics code dynamic benchmark by finite element method S. González-Pintor , D. Ginestar , G. Verdú , High order finite element method for the lambda modes problem on hexagonal geometry . Ann. Nucl. Energy . 36 , 1450 - 1462 ( 2009 ). https://doi.org/10.1016/j.anucene.2009.07.003 https://doi.org/10.1016/j.anucene.2009.07.003 High order finite element method for the lambda modes problem on hexagonal geometry M.K. Saadi , A. Abbaspour , Effective point kinetic parameters calculation in Tehran research reactor using deterministic and probabilistic methods . Nucl. Sci. Tech . 28 , 171 ( 2017 ). https://doi.org/10.1007/s41365-017-0323-7 https://doi.org/10.1007/s41365-017-0323-7 Effective point kinetic parameters calculation in Tehran research reactor using deterministic and probabilistic methods M. Arkani , M. Hassanzadeh , S. Khakshournia , Calculation of six-group importance weighted delayed neutron fractions and prompt neutron lifetime of MTR research reactors based on Monte Carlo method . Prog. Nucl. Energ . 88 , 352 - 363 ( 2016 ). https://doi.org/10.1016/j.pnucene.2015.12.005 https://doi.org/10.1016/j.pnucene.2015.12.005 Calculation of six-group importance weighted delayed neutron fractions and prompt neutron lifetime of MTR research reactors based on Monte Carlo method A. Lashkari , H. Khalafi , H. Kazeminejad , Effective delayed neutron fraction and prompt neutron lifetime of Tehran research reactor mixed-core . Ann. Nucl. Energy . 55 , 265 - 271 ( 2013 ). https://doi.org/10.1016/j.anucene.2012.12.001 https://doi.org/10.1016/j.anucene.2012.12.001 Effective delayed neutron fraction and prompt neutron lifetime of Tehran research reactor mixed-core S.A. Hosseini , N. Vosoughi , M.B. Ghofrani et al ., Calculation, measurement and sensitivity analysis of kinetic parameters of Tehran Research Reactor . Ann. Nucl. Energy . 37 , 463 - 470 ( 2010 ). https://doi.org/10.1016/j.anucene.2010.01.018 https://doi.org/10.1016/j.anucene.2010.01.018 Calculation, measurement and sensitivity analysis of kinetic parameters of Tehran Research Reactor M. Zaker , Effective delayed neutron fraction and prompt neutron lifetime of Tehran research reactor . Ann. Nucl. Energy . 30 , 1591 - 1596 ( 2003 ). https://doi.org/10.1016/S0306-4549(03)00081-1 https://doi.org/10.1016/S0306-4549(03)00081-1 Effective delayed neutron fraction and prompt neutron lifetime of Tehran research reactor B.X Yuan , H. Zeng , W.K Yang et al , Advances in FEMN: the code for nuclear noise analysis based on ICEM-CFD . Proc. ASME. ICONE26, Volume 3: Nuclear Fuel and Material, Reactor Physics, and Transport Theory, V003T02A006 , July 22-26, 2018 . https://doi.org/10.1115/ICONE26-81169 https://doi.org/10.1115/ICONE26-81169 . Advances in FEMN: the code for nuclear noise analysis based on ICEM-CFD K.P Wang , H.C Wu , L.Z Cao , The application of the analytic basis function expansion triangular-z nodal method for neutron diffusion calculation . China. Sciencepap . 10 ( 23 ), 2774 - 2778 ( 2015 ). https://doi.org/10.3969/j.issn.2095-2783.2015.23.017 https://doi.org/10.3969/j.issn.2095-2783.2015.23.017 The application of the analytic basis function expansion triangular-z nodal method for neutron diffusion calculation A. Hotta , H. Ninokata , A.J. Baratta , Development of parallel coupling system between three-dimensional nodal kinetic code ENTREE and two-fluid plant simulator TRAC / BF1 . J. Nucl. Sci. Technol , 37 ( 10 ), 840 - 854 ( 2000 ). https://doi.org/10.1080/18811248.2000.9714965 https://doi.org/10.1080/18811248.2000.9714965 Development of parallel coupling system between three-dimensional nodal kinetic code ENTREE and two-fluid plant simulator TRAC / BF1 10.1007/s41365-020-00781-y This work was supported by the National Key R&D Plan Project (No. 2019YFB1901200) and Fundamental Project of NPL, CAEP (No. 232). 2020-02-22 2020-05-11 2020-05-13 2020-06-21 2020-07-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00781-y http://dx.doi.org/10.1007/s41365-020-00781-y . 2020 Nuclear Science and Techniques 7 31