Research Article | Open Access
Wenzhen Chen, Jianli Hao, Ling Chen, Haofeng Li, "Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method", Science and Technology of Nuclear Installations, vol. 2013, Article ID 261327, 6 pages, 2013. https://doi.org/10.1155/2013/261327
Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method
The singularly perturbed method (SPM) is proposed to obtain the analytical solution for the delayed supercritical process of nuclear reactor with temperature feedback and small step reactivity inserted. The relation between the reactivity and time is derived. Also, the neutron density (or power) and the average density of delayed neutron precursors as the function of reactivity are presented. The variations of neutron density (or power) and temperature with time are calculated and plotted and compared with those by accurate solution and other analytical methods. It is shown that the results by the SPM are valid and accurate in the large range and the SPM is simpler than those in the previous literature.
The analysis of variation of neutron density (or power) and reactivity with time under the different conditions is an important content of nuclear reactor physics or neutron kinetics [1–7]. Some important achievements on the supercritical transient with temperature feedback with big () or small () reactivity inserted have been approached through the effort of many scholars [7–12]. The studies on the delayed supercritical transient with small reactivity inserted and temperature feedback are introduced in the related literature [13–15], in which the explicit function of density (or power) and reactivity with respect to time is derived mainly with decoupling method, power prompt jump approximation, precursor prompt jump approximation, temperature prompt jump approximation [10, 16], and so forth. From the detailed analysis and comparison of the results in the early and recent literature [7, 12, 14], it is found that some results have certain limit and rather big error under the particular conditions. In present work, the variation law of power, reactivity, and precursor density with respect to time at any level of initial power is obtained by the singularly perturbed method (SPM). All the results are compared with those obtained by the numerical solution which tend to the accurate solution under very small time step size . It is proved that the SPM is correct and reliable and is simpler than the analytical methods by the related literature.
2. Theoretical Derivation
The point reactor neutron kinetics equations with one group of delayed neutrons are [3, 4] where is the average neutron density, is the time, is the reactivity, is the total fraction of the delayed neutron, is the prompt neutron lifetime, is the radioactive decay constant of delayed neutron precursor, and is the average density of delayed neutron precursor. When multiplied with a certain coefficient, represents the power. It is assumed that the reactor has a negative temperature coefficient of reactivity () when a small step reactivity (β) is inserted. Consider the temperature feedback, and the real reactor reactivity is where is the temperature increment of the reactor, namely, , where and are the instantaneous temperature and initial temperature, respectively. After the reactivity is inserted into the reactor the adiabatic model is still employed [3, 15]; then we have where is the reciprocal of thermal capacity of reactor.
Substituting obtained from (1) into (6) and simplifying it yields The transient process is supposed to begin at and [15, 17], so the initial conditions can be given as , , , , where is the initial neutron density (or power).
For , , (1) and (2) are stiff equations. According to the singularly perturbed method , solution of includes the inner solution in inner part and the outer solution in outer part. In this paper both inner and outer solutions are approximated to be zero order: The initial conditions are where is the initial value of inner solution and is the initial value of outer solution.
In the outer part the inner solution attenuates to zero, and and can be neglected, namely, and ; therefore, in the outer part (12) can be simplified as follows:
Because varies slowly and s, compared to other terms, the term on the left side of (13) can be neglected, and then we have
Integrating (16) subjective to the initial conditions and yields
With the initial conditions , the solution of (18) is where .
In the inner part is assumed to be constant , and , , so (12) can be simplified as follows, respectively:
In addition, the temperature feedback is not fast enough to affect the reactivity, and in the inner part it is assumed that , . From (20) we can get
The solution of (22) is
The fast varying part is assumed to attenuate to zero in the inner part, so and
Then we can get the fast varying part in the inner part as follows:
From (25) we have
Integrating (29) with the initial conditions , and yields
3. Calculation and Analysis
The PWR with fuel 235U is taken as an example with parameters , s, 1/s, K/MW·s, and 1/K [8, 13]. For the reactor with the initial power 1 MW, while reactivity and is inserted, respectively, the variations of reactivity, temperature, power with time, and power with reactivity are presented in Figures 1, 2, 3, and 4. The curves with smaller change are for and the curves with larger change are for . The solid line notes the accurate solution by the best basic function method with very small step size . The results of this paper and the accurate solution as well as the temperature prompt jump (TPJ) method in the literature  are almost the same and hard to distinguish in Figures 1–4. The short dashed-dot line and long dashed-dot line represent the results of precursor prompt jump (PrPJ) method in the literature  and the small parameter (SmP) method in the literature , respectively. The dashed line notes the results of power prompt jump (PPJ) method in the literature . The difference is caused by the approximate treatment to obey the methods of PrPJ, SmP, PPJ, and SmP. Furthermore the variation in the vicinity of prompt supercritical process is also calculated and is shown in Figures 5 and 6. The correct results cannot be obtained by the small parameter method (SmP) in the vicinity of prompt supercritical process and are not shown in Figures 5 and 6.
From Figures 1–6 it can be concluded that very good results cannot be obtained by the precursor prompt jump (PrPJ) method to calculate the delayed supercritical progress with small step reactivity and temperature feedback.
For small step reactivity, the results by the small parameter (SmP) method are close to those by the power prompt jump (PPJ) method and are better than those by the precursor prompt jump (PrPJ) method, but the accuracy of results by the small parameter method decreases with the increase of the reactivity inserted. The power is negative when the small parameter method is used to calculate the transient process in the vicinity of prompt supercritical state. From Figures 1–4, it can be seen that the small parameter method is more suitable for the calculation of reactivity and temperature increase than for that of power.
The results are quite precise using the power prompt jump (PPJ) method for the delayed supercritical process, but the main problem compared to the accurate solution is that some displacement exists along time axis. Furthermore it should be pointed out that each power peak value obtained by the precursor prompt jump (PrPJ) method, power prompt jump (PPJ) method, or small parameter (SmP) method is lower than that obtained by the accurate solution or singularly perturbed method (SPM) see Figure 1.
From Figures 1–4 it can be also found that the temperature prompt jump method (TPJ) and the singularly perturbed method (SPM) in this paper are the two most precise methods for the delayed supercritical process with small step reactivity and temperature feedback. However from Figures 5 and 6 it can be seen that as the reactivity inserted increases to the vicinity of prompt supercritical process, the total discrepancy of power by the TPJ method is larger than that by the SPM or PPJ method, and the irrelevant phenomena that the power jumps at first and then decreases monotonously from the peak will appear in the TPJ method as shown in Figure 6.
The analytical expressions of power (or neutron density), reactivity, the precursor power (or density), and temperature increase with respect to time are derived for the delayed supercritical process with small reactivity () and temperature feedback by the singularly perturbed method. Compared with the results by the accurate solution and other methods in the literature, it is shown that the singularly perturbed method (SPM) in this paper is valid and accurate in the large range and is simpler than those in the previous literature. The method in this paper can provide a new theoretical foundation for the analysis of reactor neutron dynamics.
This research is supported by the National Natural Science Foundation of China (Project no. 11301540) and the Natural Science Foundation of Naval University of Engineering.
- H. P. Gupta and M. S. Trasi, “Asymptotically stable solutions of point-reactor kinetics equations in the presence of Newtonian temperature feedback,” Annals of Nuclear Energy, vol. 13, no. 4, pp. 203–207, 1986.
- H. Van Dam, “Dynamics of passive reactor shutdown,” Progress in Nuclear Energy, vol. 30, no. 3, pp. 255–264, 1996.
- D. L. Hetrick, Dynamics of Nuclear Reactors, American Nuclear Society, La Grange Park, Ill, USA, 1993.
- W. M. Stacey, Nuclear Reactors Physics, Wiley-Interscience, New York, NY, USA, 2001.
- A. A. Nahla and E. M. E. Zayed, “Solution of the nonlinear point nuclear reactor kinetics equations,” Progress in Nuclear Energy, vol. 52, no. 8, pp. 743–746, 2010.
- G. Espinosa-Paredes, M.-A. Polo-Labarrios, E.-G. Espinosa-Martínez, and E. D. Valle-Gallegos, “Fractional neutron point kinetics equations for nuclear reactor dynamics,” Annals of Nuclear Energy, vol. 38, no. 2-3, pp. 307–330, 2011.
- A. E. Aboanber and A. A. Nahla, “Solution of the point kinetics equations in the presence of Newtonian temperature feedback by Padé approximations via the analytical inversion method,” Journal of Physics A, vol. 35, no. 45, pp. 9609–9627, 2002.
- A. E. Aboanber and Y. M. Hamada, “Power series solution (PWS) of nuclear reactor dynamics with newtonian temperature feedback,” Annals of Nuclear Energy, vol. 30, no. 10, pp. 1111–1122, 2003.
- W. Z. Chen, B. Zhu, and H. F. Li, “The analytical solution of point-reactor neutron-kinetics equation with small step reactivity,” Acta Physica Sinica, vol. 50, no. 8, pp. 2486–2489, 2004 (Chinese).
- H. Li, W. Chen, F. Zhang, and L. Luo, “Approximate solutions of point kinetics equations with one delayed neutron group and temperature feedback during delayed supercritical process,” Annals of Nuclear Energy, vol. 34, no. 6, pp. 521–526, 2007.
- T. Sathiyasheela, “Power series solution method for solving point kinetics equations with lumped model temperature and feedback,” Annals of Nuclear Energy, vol. 36, no. 2, pp. 246–250, 2009.
- Y. M. Hamada, “Confirmation of accuracy of generalized power series method for the solution of point kinetics equations with feedback,” Annals of Nuclear Energy, vol. 55, pp. 184–193, 2013.
- Q. Zhu, X.-L. Shang, and W.-Z. Chen, “Homotopy analysis solution of point reactor kinetics equations with six-group delayed neutrons,” Acta Physica Sinica, vol. 61, no. 7, Article ID 070201, 2012.
- S. D. Hamieh and M. Saidinezhad, “Analytical solution of the point reactor kinetics equations with temperature feedback,” Annals of Nuclear Energy, vol. 42, pp. 148–152, 2012.
- A. A. Nahla, “An analytical solution for the point reactor kinetics equations with one group of delayed neutrons and the adiabatic feedback model,” Progress in Nuclear Energy, vol. 51, no. 1, pp. 124–128, 2009.
- W. Chen, L. Guo, B. Zhu, and H. Li, “Accuracy of analytical methods for obtaining supercritical transients with temperature feedback,” Progress in Nuclear Energy, vol. 49, no. 4, pp. 290–302, 2007.
- H. Li, W. Chen, L. Luo, and Q. Zhu, “A new integral method for solving the point reactor neutron kinetics equations,” Annals of Nuclear Energy, vol. 36, no. 4, pp. 427–432, 2009.
- Z. Q. Huang, Kinetics Base of Nuclear Reactor, Peking University Press, Beijing, China, 2007 (Chinese).
Copyright © 2013 Wenzhen Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.