Abstract
The environment effect arises when pin-cell homogenized parameters are generated with reflective boundary conditions. To treat it in whore-core pin-by-pin calculation, two works are summarized in this article. Firstly, by analyzing the relative errors of pin-cell homogenized group constants and the relative importance of pin-cell discontinuity factors (PDF) in each group, the importance of correcting the PDF of the thermal group is recognized. Secondly, the least-squares method for a multivariate polynomial is utilized to functionalize the relation of the thermal group PDF and the core parameters, including diffusion coefficient, removal cross-section, neutron source, and normalized surface flux. The C5G7 and KAIST benchmarks are employed to evaluate the performance of the PDF predication. Numerical results indicate its effectiveness in reducing the errors of eigenvalue and pin power, especially for the cases with the fuel pins located near the interface between different assemblies.
Introduction
For Pressurized Water Reactor (PWR) simulation, the computational cost of one-step calculation with fully detailed description is still too expensive for routine calculation using either stochastic or deterministic methods, even with the currently most advanced computing powers. Therefore, approximations in spatial, energy, and angular spaces have been developed to provide efficient solutions with acceptable accuracy. With the advantages of small storage and high computing speed, the two-step calculation schemes have successfully been applied to PWR analysis for decades.
In the last few decades, the nodal diffusion calculation accompanied by assembly homogenization has been widely employed (). It contains 2D lattice neutron-transport calculation and 3D whole-core neutron diffusion calculation. The 2D lattice neutron-transport calculations are carried out for each type of assemblies with reflective boundary conditions to provide assembly-homogenized parameters (including cross-sections, diffusion coefficients, and discontinuity factors) and pin-power form functions. The whole-core diffusion solver determines the nodal power shape and then reconstructs the 3D pin-power distribution (). One of the most important error sources in the scheme is applying the reflective boundary condition in the lattice calculation. For most large-scale PWRs, the neutron leakage model applied in the assembly homogenization can lead to an accurate result without any other corrections for the reflective boundary condition. There are several reasons. First, the environmental effect is not so strong for large-scale commercial PWRs. Second, fuel assemblies are optically thick which contain more than five neutron-free paths.
With the development of computing technology, an improved scheme named pin-by-pin calculation has become popular in recent years (; ). Different from the traditional two-step calculation, only the heterogeneous structure within each pin is homogenized, leaving the assembly heterogeneous during the three-dimensional whole-core calculation. It eliminates the assembly homogenization and the pin-power reconstruction. Alone with the development of the pin-cell homogenization technology (; ), the reflective boundary condition used in the process of the pin-cell homogenized parameters becomes the main error source. In principle, the homogenized cross-sections, diffusion coefficients, and the pin-cell discontinuity factors can preserve the neutron leakage and reaction rates only for the boundary conditions at which they were created, i.e., for the reflective boundary condition. However, the exact boundary condition of an assembly in the active core cannot be known a priori for the two-step scheme. What is more, the size of the homogenization region in the pin-by-pin calculation is almost equal to the averaged neutron-free path, making the pin-cell homogenized parameters more dependent on the assembly environment compared with the assembly-homogenized parameters.
In order to achieve further accuracy of the pin-by-pin calculation, several embedding homogenization methods have been estimated (; ). These methods usually utilize the embedded heterogeneous assembly calculations to update the group constants on-line. The artificial neural network (ANN) method has also been developed to fit the group constants (; ). However, the methods are unstable because of the unpredictable ANN training and uncertainty in the quality of the fitted function. Focused on the environment effect, first in Analysis of Pin-Cell Homogenized Parameters, this article analyzed the errors of the pin-cell homogenized parameters and the relative importance of pin-cell discontinuity factors in each group (). Secondly, the least-squares method (LSM) utilized for the functionalization of PDF () is introduced in the Method for the PDF Predication section. In order to numerically analyze the performance of the method, C5G7 () and KAIST () benchmarks are evaluated in this summary.
Analysis of Pin-Cell Homogenized Parameters
A typical checkerboard assembly problem based on the KAIST benchmark (), as shown in Figure 1, is utilized to analyze the pin-cell homogenized data. In the pin-cell homogenization, the flux-volume-weight method is used to generate the homogenized cross-sections:where h, g stands for the energy group and stands for the homogenized cross-section.
FIGURE 1
Because of the reflective boundary condition performed in the assembly calculations, neither the axial nor radial leakage rate is considered within the heterogeneous assembly calculation. The heterogeneous leakage model () is applied in this work; the heterogeneous leakage equations are shown as follows:
The pin-cell diffusion coefficient for core pin-by-pin calculation would be determined by the space-dependent leakage coefficient from the heterogeneous leakage model.
The GET method was applied in this work. For single-assembly calculation with a specific boundary condition (zero current), the heterogeneous information can be provided by lattice code, while the homogeneous information is supposed to be obtained from the homogeneous neutron diffusion calculation. To get rid of the nonlinear relationship between the homogeneous surface flux and the PDF, a procedure is used to estimate the homogeneous surface flux accurately for the SP3 method using the heterogeneous information (). The discontinuity factor is defined as the ratio of the heterogeneous over the homogeneous surface fluxes:where and stand for the pin-cell heterogeneous surface flux and the homogeneous one.
Two sets of 7-group pin-cell homogenized parameters are generated for the colorset problem:
• The first set is obtained from the whole-core checkerboard calculation in order to eliminate the single-assembly zero net current boundary conditions assumption and incorporate the correct interface spectrum effect between UO2 and MOX assemblies into the homogenized parameters.
• The pin-cell homogenized parameters in the second case are determined by single-assembly calculations with reflective boundary conditions.
Pin-Power Distribution
Figure 2 shows the pin-power relative error distribution of the pin-by-pin SP3 calculation utilized in the second library set. The reference solution is provided by the 2D whole-core one-step transport calculation. It can be found that large pin-power errors always occur in the pins close to the interface between two assemblies due to the large spectral mismatch between the single-assembly and full-core conditions in those locations. The maximum of the relative errors is more than 4% at the interface, while that of other locations would be about 1%.
FIGURE 2
Pin-Cell Homogenized Parameters
Figure 3 shows the neutron flux profile between the UO2 and MOX assemblies. Dramatic changes can be found close to the interface of the thermal energy group. The spectrum interference effect between UO2 and MOX is very strong for the fast groups (5–7 groups) while being weak for the thermal group (1–4 groups).
FIGURE 3
Part of the relative error distributions of the cross-sections is shown in Figure 4. It can be found that the errors of the pin-cell homogenized cross-sections are small, and the spectrum interference effect of the thermal group is so strong. The maximum error is less than 1.5%.
FIGURE 4
The relative error distributions of the PDF ratio of the first and seventh energy groups are shown in Figure 5. It can be found that the errors of the PDF ratios of fast groups are negligible. It is because that the spectrum interference effect of the fast groups is very weak. In contrast, the errors of the PDF ratios of the thermal groups are very large because of the strong spectrum interference effect between different assemblies. The maximum PDF error of the thermal groups would be more than 10%.
FIGURE 5
Relative Importance of PDF in Each Group
Single-assembly calculations of four different assemblies mentioned in the checkerboard problem were carried out to determine the relative importance of PDF in each group and analyze the interaction effect of the PDF between different energy groups.
Several calculations are carried out as follows:
• The reference solutions are obtained from the transport calculation by the house-developed Bamboo lattice code, which uses the MOC method as the transport solver (; ; ; ).
• The second one is the pin-by-pin SP3 calculation without PDF (; ; ; ).
• The third results come from the calculation with exact PDF in fast groups or thermal groups and then turned off (set to 1.0) in all other groups.
From the results presented in Table 1, several important conclusions can be made. First, the results with the correct PDF of the thermal group are more accurate than the calculation with the correct PDF of the fast group. It is indicated that the PDF of the thermal group are more important and contribute more to the error than the fast group ones. Secondly, the sum of the eigenvalue differences of each exact PDF equals the calculation results without PDF. It is indicated that the PDF satisfy the neutron balance equation in a group-by-group way. Thirdly, correct PDF in one group with those of other groups set as 1.0 might not improve the solution accuracy; however, they do not make the results worse.
TABLE 1
| Case | kinf | Eigenvalue difference (pcm) pin-power % RMS error | ||
|---|---|---|---|---|
| Without PDF | PDF effect | |||
| Fast group = 1.0 Thermal group = exact | Fast group = exact Thermal group = 1.0 | |||
| UOX-1 | 1.49161 | 141 | 28 | 113 |
| 0.57 | 0.09 | 0.59 | ||
| UOX-2 | 1.44302 | 453 | −9 | 461 |
| 1.09 | 0.19 | 1.04 | ||
| MOX-1 | 1.40060 | 847 | −31 | 881 |
| 1.86 | 0.56 | 1.47 | ||
| MOX-1 CR | 1.35576 | 1,429 | 18 | 1,410 |
| 1.63 | 0.34 | 1.45 | ||
Single-assembly calculation results.
Method for the PDF Predication
The thermal group PDF have the largest relative errors and the thermal group is more important. Thus, attention should be paid to the correction of the thermal group PDF. In this work, the least-squares method (LSM) for multivariate polynomial is utilized to functionalize the relation of the thermal group PDF and the core parameters. The PDF predication function would be created with heterogeneous solutions obtained from several single-assembly calculations and several spectral geometry calculations. It is expected to predicate PDF with a whole-core homogeneous solution.
The least-squares method (LSM) for multivariate polynomial is written in the way by assuming that x is a vector variable with two or more components:where stands for the function variable, which is the pin-cell homogenized parameters in our study, such as the cross-section and the pin-cell discontinuity factors; stands for the core-parameter vector x, is the undetermined coefficient, which can be determined using the discrete state points (xn, yn).
The selection of the core parameters for the function of PDF predication plays an important role in the accuracy of the solutions. In this article, several elements are taken into consideration to pick up the core parameters. First, neutron behaviors are different in different energy groups. Therefore, the core parameters should be selected in a group-by-group manner. Second, the parameters should be associated with the method used in the core nodal solver. Different core methods would lead to different PDF. Third, the parameters should be suitable for any case that appeared in the problem.
In this article, the choice of core parameters is based on the physics underlying core neutronics behavior: neutron leakage, removal, and generation. Thus, there are seven core parameters as follows:where subscripts i and g stand for the node and the energy group; represents the diffusion coefficient (cm); represents the removal cross-section (cm−1); S is the sum of the fission source and the scatter source (cm−3); is the average neutron flux (s−1 cm−2); is the nodal surface flux (s−1 cm−2); superscript s represents pin-cell surface in the different directions, .
Based on the selected seven inputs, the coefficient of the PDF predication function can be determined using the discrete state points (xn, yn) via the least-squares method. For the whole-core pin-by-pin calculation, the state points (xn, yn) would be provided by the single-assembly transport calculations and colorset problem calculation.
In order to determine the order of the PDF predication function, several calculations of different orders are carried out. The relative fitting errors of the PDF predication function are presented in Table 2. It can be found that when the polynomial order of the function is three, the maximum root-mean-square (RMS) error is less than 1%.
TABLE 2
| Polynomial order | Number of coefficients | RMS error/% | ||||||
|---|---|---|---|---|---|---|---|---|
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | ||
| 0 | 1 | 2.09 | 1.21 | 1.82 | 5.69 | 5.77 | 8.33 | 8.17 |
| 1 | 8 | 1.42 | 0.29 | 0.41 | 1.25 | 2.14 | 3.11 | 3.27 |
| 2 | 36 | 0.50 | 0.15 | 0.22 | 0.69 | 0.95 | 1.66 | 0.95 |
| 3 | 120 | 0.30 | 0.09 | 0.05 | 0.10 | 0.30 | 0.26 | 0.23 |
Relative fitting errors of the different energy-group PDF by different predication function.
In this work, the function of PDF predication is selected as follows:
Taking into consideration the different neutron behaviors between the internal and peripheral locations of the core, two different functions of PDF predication were established for the internal fuel pins and the peripheral fuel, respectively. The PDF predication of reflector pins is considered separately.
Numerical Results
The 2D C5G7 benchmark () shown in Figure 6 is carried out and the results are summarized in Table 3. The pin-power error distributions are shown in Figure 7. The “MC” in Table 3 represents the reference solutions of Monte Carlo calculation obtained from the benchmark. The “SA” represents that the whole-core SP3 calculation is carried out with the cross-sections and the PDF generated by the single-assembly calculation. The “LSM” represents the PDF of thermal groups obtained by the predication function and the other parameters are just the same as the “SA.” The reference keff is 1.18655 ± 0.00006 within the confidence interval of 98%.
FIGURE 6
TABLE 3
| Solutions | keff | Eigenvalue difference/pcm | Pin-power Max/RMS error (%) |
|---|---|---|---|
| MC | 1.18655 | — | — |
| SA | 1.18493 | −162 | 9.62/2.18 |
| LSM | 1.18772 | 117 | 3.09/0.57 |
Summary of the results.
FIGURE 7
Takingg into consideration the different neutron behaviors between the internal and peripheral locations of the core, two functions of PDF predication were established for the internal fuel pins and the peripheral fuel or reflector pins, respectively. Numerical results indicate that the application of LSM can accurately predicate the PDF of the thermal group and reduce the error of eigenvalue and pin power effectively. Compared with the “SA” calculation, which causes about −162 pcm error in keff and 9.62% maximum pin-power error, the “LSM” calculation reduces these errors, respectively, to 117 pcm and 3.09%. The root-mean-square (RMS) pin-power error is also reduced from 2.18 to 0.57%.
The 2D whole-core problem of the KAIST benchmark () shown in Figure 8A is also carried out to analyze the performance of the method. Figure 8B is the pin-power distribution of the reference solution, which comes from the one-step transport calculation. The “SA” and “LSM” have the same meaning as those in the C5G7 benchmark problem. Results are summarized in Table 4. The pin-power error distributions are shown in Figure 9. It leads to the same conclusion. The function of thermal PDF predication works effectively. It can reduce the error and improve the accuracy of the pin-power distribution, especially for the interface locations between different assemblies.
FIGURE 8
TABLE 4
| Solutions | keff | Eigenvalue difference/pcm | Pin-power max/RMS error (%) |
|---|---|---|---|
| MOC | 0.97948 | — | — |
| SA | 0.97914 | −34 | 3.29/0.70 |
| LSM | 0.97896 | −52 | −2.13/0.70 |
Summary of the results.
FIGURE 9
Conclusion
All the pin-cell homogenized parameters are analyzed based on a checkerboard problem. It can be found that the errors of the PDF ratios of the thermal groups are very large, while others are negligible. What is more, the relative importance of the PDF of thermal groups is larger than that of other parameters. Attention would be paid to the correction of the thermal group PDF.
The PDF of the thermal group is functionalized using the least-squares method. The PDF predication function is created with given solutions. It has been proved that the predicated PDF of thermal groups can lead to a more accurate core solution.
Results of the C5G7 and KAIST benchmarks indicate that the PDF predication of the thermal group works effectively in reducing the error of eigenvalue and pin power, especially for the pins located near the interface between different assemblies.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
YL contributed significantly to analysis and manuscript preparation; HW helped perform the analysis with constructive discussions.
Funding
This research was carried out under financial support by the National Natural Science Foundation of China (No.11775171).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
SmithM. A.LewisL. L.NaB-C. (2004). Benchmark on Deterministic 2-D MOX Fuel Assembly Transport Calculations Without Spatial Homogenization. Progress Nuclear Energy45, 107–118. 10.1016/j.pnueene.2004.09.003
2
ChenQ.WuH.CaoL. (2008). Auto MOC-A 2D Neutron Transport Code for Arbitrary Geometry Based on the Method of Characteristics and Customization of AutoCAD. Nucl. Eng. Des.238, 2828–2833. 10.1016/j.nucengdes.2008.04.014
3
ChoN. Z. (2000). Benchmark Problem in Reactor and Particle Transport Physics. Daejeon, Korea. Korea Advanced Institute of Science and Technology.
4
HeL.WuH.CaoL. (2014). Improvements of the Subgroup Resonance Calculation Code SUGAR. Ann. Nucl. Energ.66, 5–12. 10.1016/j.anucene.2013.11.029
5
HuangKai.WuH. C.CaoL. Z. (2016). Improvements to the Transmutation Trajectory Analysis of Depletion Evaluation. Ann. Nucl. Energ.87, 637–647. 10.1016/j.anucene.2015.10.013
6
KozlowskiT.LeeD. J.DownerT. (2004). The Use of an Artificial Neural Network for On-Line Predication of Pin-Cell Discontinuity Factors in PARCS. Chicago, Illinois, USA: PHYSOR2004.
7
KozlowskiT. (2005). Spatial Homogenization Methods for Pin-By-Pin Neutron Transport Calculations, PhD Thesis. Westville, USA: Purdue University.
8
KozlowskiT. (2004). The Use of an Artificial Neural Network for On-Line Prediction of Pin-Cell Discontinuity Factors in PARCS. Chicago, Illinois, USA: PHYSOR2004.
9
KozlowskiT.XuY.DownarT. J.LeeD. (2011). Cell Homogenization Method for Pin-By-Pin Neutron Transport Calculations. Nucl. Sci. Eng.169, 1–18. 10.13182/nse08-85
10
LiY. Z.CaoL. Z.YuanX. B. (2014). High Order Source Approximation for the EFEN Method. Japan: PHYSOR.
11
LiY.ZhangB.WuH.ShenW. (2017). Heterogeneous Neutron-Leakage Model for PWR Pin-By-Pin Calculation. Ann. Nucl. Energ.110, 443–452. 10.1016/j.anucene.2017.07.002
12
LiY. Z.WuH. C.CaoL. Z. (2013). Unstructured Triangular Nodal-SP3 Method Based on an Exponential Function Expansion. Nucl. Sci. Eng.174 (2), 163–171. 10.13182/nse11-111
13
LiY. Z.ZhangB.TianC. (2015). NECP-CACTI: Pressurized Water Reactor Lattice Code Development. Tran. Nucl. Soc.112, 814–816.
14
ParkK.ParkC. (2001). Pin-cell Homogenization via Generalized Equivalence Theory and Embedding Assembly Calculation. Trans. Am. Nucl. Soc.85, 334.
15
SmithK. S. (1986). Assembly Homogenization Techniques for Light Water Reactor Analysis. Prog. Nucl. Energ.17, 303–335. 10.1016/0149-1970(86)90035-1
16
SmithK. (1980). Spatial Homogenization Methods for Light Water Reactors, PhD Thesis. MIT.
17
SugimuraN.UshioT. (2006). Calculation Models of AEGIS/SCOPE2, A Core Calculation System of Next Generation. Vancouver, BC, Canada: PHYSOR2006.
18
YamajiK.MatsumotoH.NakanoM. (2006). Development of the New Pin-By-Pin Core Calculation Method with Embedded Heterogeneous Assembly Calculation. Vancouver, BC, Canada: PHYSOR2006.
19
YangW.ZhengY. Q.WuH. C.CaoL. Z.LiY. Z. (2014). High-Performance Whole Core Pin-by-Pin Calculation Based on EFEN-SP3 Method. Nucl. Power Eng.35, 164–167. 10.13832/j.jnpe.2014.05.0164
20
YangW.LiY. Z.WuH. C. (2016). Multi-group GMRES Algorithm for the Exponential Function Expansion Nodal SP3 Method. Idaho, USA: PHYSOR 2016, Sun Valley.
21
YuH.RehmanH.KimY. H. (2017). Artificial Neural Network Modeling for 2-group Pin-wise Group Constants. San Francisco, CA, USA: ANS Annual Meeting.
22
ZhangB.WuH.LiY.CaoL.ShenW. (2017). Evaluation of Pin-Cell Homogenization Techniques for PWR Pin-By-Pin Calculation. Nucl. Sci. Eng.186 (2), 134–146. 10.1080/00295639.2016.1273018
23
ZhangS. H.TangC. T.HuangH. (2008). Feasibility of Embedding Nodal Homogenization in Next Generation Methods for 3D Pin-By-Pin Core Simulation. Interlaken, Switzerland: PHYSOR2008.
Summary
Keywords
PWR, pin-by-pin, environment effect, least-squares method, functionalization
Citation
Zhang B, Li Y and Wu H (2021) Environment Effect Treatments in PWR Whole-Core Pin-by-Pin Calculation. Front. Energy Res. 9:760339. doi: 10.3389/fenrg.2021.760339
Received
18 August 2021
Accepted
07 September 2021
Published
01 October 2021
Volume
9 - 2021
Edited by
Tengfei Zhang, Shanghai Jiao Tong University, China
Reviewed by
Kun Zhuang, Nanjing University of Aeronautics and Astronautics, China
Xubo Ma, North China Electric Power University, China
Zhuo Li, Sun Yat-sen University, China
Updates
Copyright
© 2021 Zhang, Li and Wu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bin Zhang, jsszwtzb@163.com
This article was submitted to Nuclear Energy, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.