**World Journal of Nuclear Science and Technology**

Vol.05 No.01(2015), Article ID:53644,14 pages

10.4236/wjnst.2015.51006

A New Formulation to the Point Kinetics Equations Considering the Time Variation of the Neutron Currents

Anderson Lupo Nunes^{1}, Aquilino Senra Martinez^{1}, Fernando Carvalho da Silva^{1}, Daniel Artur Pinheiro Palma^{2}

^{1}Department of Nuclear Engineering, COPPE/UFRJ, Rio de Janeiro, Brazil

^{2}Brazilian Commission for Nuclear Energy, Rio de Janeiro, Brazil

Email: anunes@con.ufrj.br, aquilino@lmp.ufrj.br, fernando@con.ufrj.br, dpalmaster@gmail.com

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 6 January 2015; accepted 27 January 2015; published 29 January 2015

ABSTRACT

The system of point kinetics equations describes the time behaviour of a nuclear reactor, assuming that, during the transient, the spatial form of the flux of neutrons varies very little. This system has been largely used in the analysis of transients, where the numerical solutions of the equations are limited by the stiffness problem that results from the different time scales of the instantaneous and delayed neutrons. Its derivation can be done directly from the neutron transport equation, from the neutron diffusion equation or through a heuristics procedure. All of them lead to the same functional form of the system of differential equations for point kinetics, but with different coefficients. However, the solution of the neutron transport equation is of little practical use as it requires the change of the existent core design systems, as used to calculate the design of the cores of nuclear reactors for different operating cycles. Several approximations can be made for the said derivation. One of them consists of disregarding the time derivative for neutron density in comparison with the remaining terms of the equation resulting from the P1 approximation of the transport equation. In this paper, we consider that the time derivative for neutron current density is not negligible in the P1 equation. Thus being, we obtained a new system of equations of point kinetics that we named as modified. The innovation of the method presented in the manuscript consists in adopting arising from the P1 equations, without neglecting the derivative of the current neutrons, to derive the modified point kinetics equations instead of adopting the Fick’s law which results in the classic point kinetics equations. The results of the comparison between the point kinetics equations, modified and classical, indicate that the time derivative for the neutron current density should not be disregarded in several of transient analysis situations.

**Keywords:**

Reactor Point-Kinetics, Neutron Current Density, Nuclear Power Density

1. Introduction

In order to determine the nuclear power distribution in a reactor core, one should investigate the neutron transport in a heterogeneous medium and with strong neutron absorption, where these neutrons can also be scattered or escape from the active part of the reactor. Notwithstanding that, the advances in computer processing and of the countless methods to solve the neutron transport equation, in practice the approximation of the neutron diffusion is largely used in stationary calculations to predict the distribution of neutrons and of the critical concentration of boron. To deal with the movement of the neutrons in a way similar to that of heat diffusion, one needs to make several approximations in the transport equation, which include a weak angular dependency of the angular distribution of the neutrons, isotropic sources of neutrons, and the disregarding of the derivative for neutron current density, in comparison with other terms that appear in the neutron transport equation [1] .

Once the spatial distribution of the neutrons in the nuclear reactor is known, it is also important to predict the time behaviour of this distribution, induced as it is by the variation in nuclear reactivity due to the variation of fuel temperature, the variation of the material composition of the reactor core, the variation in moderator density, amongst others. The simplest way to determine the time behaviour for the nuclear power is through the solution of point kinetics equations. These equations include approximations that are added to those made to obtain the equations for neutron diffusion in the structure of multi-groups of neutron energy [2] . Point kinetics equations consist on a system for the calculation of the nuclear power and the concentrations of the delayed neutrons precursors. They are first-order differential equations, coupled and non-linear in their more general form.

Though quite questionable, the approximations made in the development of the classical point kinetics equations have already been widely analysed and discussed in the literature. However, the influence on these equations from not considering the time derivative in neutron current density did not deserve, until now, a systematic and specific evaluation. Due to this, in this paper we developed, from the neutron transport equation, the modified point kinetics equations, that are different from the classical ones, as they include the time derivative of neutron current density.

In Section 2, it is presented the development of the modified point kinetics equations. Section 3 presents the calculation to obtain their analytical solutions. Section 4 presents the results of the analytical solutions of the classical and modified point kinetics equations. And Section 5 discusses the results obtained and provides the conclusions of this paper.

2. The Modified Point Kinetics Equations

The theory of neutron transport is the wide model to describe neutron distribution in a nuclear reactor. It is described in [1] and [3] in terms of the angular flux of neutrons,:

(1)

and

(2)

where and the operators, , and are defined as follows:

(3)

(4)

(5)

and

(6)

The scattering cross section can be expanded in terms of the polynomials of Legendre up to the second term, that is, the expansion is done for and. It consists of the approximation

(7)

Equation (3) can be re-written thus:

(8)

We apply the operator to Equations (1) and (2) and considering the following definitions:

and

it results that:

(9)

and

(10)

where the operators, and are defined:

(11)

(12)

and

(13)

In replacing Equation (8) in Equation (1), multiplying the resulting equation by and after that integrating in the solid angle, it results that:

(14)

Considering the approximation described in [1] :

and defining the transport cross section:

,

from Equation (14), it is possible to write:

(15)

There is a similarity of Equation (15) with the Telegrapher Equation, as can be seen in [4] and [5] . However the methodology used to arrive at Equation (15) is totally different, as well as the results that follow until we get to the Equation (28).

In dividing Equation (15) by the transport cross section, using the definition of diffusion coefficient, applying

the diverging operator and disregarding the term, it results that:

(16)

where the diffusion coefficient is defined as:

After that, in deriving Equation (9) in relation to time, and multiplying by, and adding to the

very Equation (9) and replacing Equation (16), one obtains:

(17)

In the Equations (10) and (17), when in a stationary regime, all of its time derivatives are disregarded. Then, it results that:

(18)

Considering the adjoint flux of neutrons from Equation (18), this adjoint equation is:

(19)

where,

(20)

(21)

(22)

and

(23)

Equation (18) is integrated in the volume and in the energy, and re-written:

(24)

In multiplying Equation (17) by the adjoint flux of neutrons and integrating in the volume and in the energy, and subtracting of the result by Equation (24), one obtains:

(25)

where, by simplicity, the functional dependence both of the diffusion coefficient as of the speed was omitted. The neutron flux is written as the product of an amplitude factor, which is dependent on time only, and a shape factor (or shape function); thus,

. (26)

In writing the neutron flux as the product of the two factors in Equation (26), the intent is that the amplitude factor, , should describe most of the time dependence whereas the shape factor, , will change very little with time.

In defining the integral:

(27)

one obtains, after replacing Equation (26) in Equation (25) and the division of the result by, that

(28)

It was considered that the time variation of the diffusion coefficient and of the neutron cross sections are negligible, which immediately implies that their derivatives are null. The several kinetic parameters that appear in Equation (28) are defined as:

(29)

(30)

(31)

(32)

(33)

(34)

and

(35)

In multiplying Equation (10) by the adjoint flux of neutrons and after that integrating in the volume and in the energy, after using Equation (26), one obtains that:

(36)

Dividing Equation (36) by and using Equations (32) and (35), one obtains:

(37)

where.

Equations (28) and (37) form the new model for the point kinetics, called modified point kinetics.

3. Equivalency of the Modified Point Kinetics Equations with Respect the Classic Point Kinetics Equations

The classical point kinetics Equations can be obtained in many ways. Following the development which resulted in the modified point kinetic Equations, the model for the classical point kinetics is obtained by making the approximation known as Fick’s Law in Equation (15),

(38)

From there on, some terms of the first and second derivative of and the first-order derivative of appear in Equation (17) and do not appear in the corresponding equation for the classical kinetics. When tends to zero, Equation (28) results to the classical point kinetics equation.

(39)

The new parameters and can be called as neutron transport frequency and neutron absorption frequency, respectively. When tends to zero the Equation (28) lies in Equation (39), in other words, the point classical kinetic equation. Equivalent to state that the neutron transport frequency is much larger than the other parameters.

Considering a homogeneous medium Equations (30) and (34) are simplified. First, the diffusion coefficient constant and the speed are considered as the medium is homogeneous, which gives:

(40)

Note that appears in the definition of Equation (29) into Equation (40).

(41)

Therefore,

(42)

In Equation (34) replace the operators, and.

(43)

Considering the homogeneous medium rewrite the Equation (43):

(44)

Replacing the Equations (27), (29) and (42) in Equation (44):

(45)

once:

(46)

follows:

(47)

4. Analytical Solution of the Modified Point Kinetics Equations for One Group of Precursors with Constant Reactivity

The solution for point kinetics equations can be obtained in several ways, as in [6] -[8] . Point kinetics equations without the approximation related to the time derivative of neutron current density, for one group of precursors and constant reactivity, according to Equations (28) and (37), are:

(48)

and

(49)

The initial conditions are:

(50)

(51)

and

(52)

Applying the conditions (50), (51) and (52) in Equation (48) it is possible to write:

(53)

One therefore chooses

(54)

It is possible to verify that the initial conditions (50), (51) and (52) as applied to Equation (38), that is, to the classical point kinetics, result exactly in Equation (54). So, in replacing Equation (54) in Equation (53) it results:

(55)

It adds to the Equation (48) with Equation (49):

(56)

You can also replace defined by Equation (49) in Equation (48):

(57)

The Equation (57) is derived in the time. Then, it results that:

(58)

Multiplying the Equation (56) by:

(59)

It adds to the Equation (58) with the Equation (59):

(60)

The Equation (60) is a third-order homogeneous linear differential equation and will be solved by the Laplace Transform technique. One can simplify Equation (60) as:

(61)

In applying the Laplace transform [9] in Equation (61), it results that:

(62)

being,

(63)

and

(64)

In Equation (62) we isolate the Laplace transform of the neutron density:

(65)

Applying the reverse Laplace transform on both sides of Equation (65) one have that:

(66)

To solve the reverse Laplace transform one should factor the polynomials so to obtain:

(67)

Note that the Equation (67) corresponds exactly to Equation (66), in other words, can rewrite (67):

Comparing the above equations it is possible to deduce the number system:

(68)

(69)

(70)

(71)

(72)

and

(73)

The Equation (73) can be solved directly to obtain the value of. Note that it is an equation of the third de-

gree. Initially it replaces the variable by:

(74)

Note that in Equation (74) the term corresponding to is reset. Thus, it can be rewritten:

(75)

The Equation (75) was resolved according [10] and the solution is known as Cardano-Tartaglia’s formula:

(76)

Note that Equation (76) corresponds only one of the three solutions of the Equation (85). The other two solutions can be obtained by dividing Equation (75) by which allow us to obtain:

therefore,

Follows:

(77)

The roots of Equation (77) are the other two solutions of Equation (76). Thus three possible solutions which satisfy the Equation (73) are obtained. The solutions are substituted in Equation (72). Equation (72) becomes an Equation of the second degree variable and can be trivially solved. With solutions and can be obtained from Equation (71). Equations (68), (69) and (70) are all expressed only in terms of, and and so the values of, and are determined directly. From the resolution of the inverse Laplace transform of Equation (67), results that:

(78)

5. Results

Modified point kinetics equations for one group of precursors with constant reactivity are solved with the method described in Section 3. Numerical values for the nuclear parameters are considered from references found in the literature, i.e., [1] [11] -[15] as listed in Table 1.

The values considered in Table 1 for the absorption cross section and the diffusion coefficient correspond to the average cross sections of a typical core in a PWR nuclear reactor, according to [1] and [12] .

The initial conditions (50), (51), (52) are used in the calculations of classical and modified equations and the initial conditions (54) and (55) are used only in the modified point kinetics equations. We consider and the values of Table 1. The results are shown in the Figures 1-6 and in the Table 2 to Table 3.

Figures 1-6 show the results of the calculations made with the solution of classical point kinetics equation, modified point kinetics equation with transport frequency equal to 104 s^{−1} and modified point kinetics equation with transport frequency equal to 103 s^{−1}. In the first three figures the time interval is from 0 to 100 s and in the last three it goes from 0 to 10 s. Note that in the graph contained in Figure 3 the order of magnitude for neutron density of the classical kinetics and of the modified kinetics are quite different for a reactivity of 0.007.

The variation in the neutron density as a function of the reactivity is seen through a comparison between the graphs. It is possible to see that, for a reactivity equal to the fraction of neutrons delayed by the total of neutrons, the neutron density obtained by the classical point kinetics equations for a time corresponding to 100 s is of the

Figure 1. Neutron density as a function of the time of 0 s at 100 s for a reactivity of 0.001.

Table 1. Parameters used in the tests.

Table 2. Calculation of n(t) (cm^{−3}) with point kinetics for a group of precursors with a neutron transport frequency of 10^{4} s^{−1}.

Table 3. Calculation of n(t) (cm^{−3}) with point kinetics for a group of precursors with a neutron transport frequency of 10^{3} s^{−1}.

Figure 2. Neutron density as a function of the time of 0 s at 100 s for a reactivity of 0.003.

Figure 3. Neutron density as a function of the time of 0 s at 40 s for a reactivity of 0.006.

Figure 4. Neutron density as a function of the time of 0 s at 10 s for a reactivity of 0.001.

Figure 5. Neutron density as a function of the time of 0 s at 10 s for a reactivity of 0.003.

Figure 6. Neutron density as a function of the time of 0 s at 10 s for a reactivity of 0.001.

order of 1022, whilst the use of the modified point kinetics equations is of the order of 109, for a neutron transport frequency of 103 s^{−1}. Thus, the difference between the results of the models is more significant for high- reactivity situations.

6. Conclusions

The objective of this paper is to obtain a new system of equations called equations of point kinetics modified in which is considered the effect of the time derivative for neutron current density in the Equation (15). In general, the time derivative of the density of neutrons is neglected for the obtainment of the classical model.

The results presented in this article show that the difference between the neutron density obtained from classical point kinetics equations and that obtained from modified point kinetics equations is relevant. With the neutron transport frequency equal to 104 s^{−1} the difference between the neutron density obtained from classical point kinetics equations and that obtained from point kinetics equations without the approximation for the time derivative of neutron current density is relevant. With a neutron transport frequency equal to 103 s^{−1} the difference between them it is quite significative.

Modified point kinetics equations imply a significant difference in results, in relation to those obtained with the classical point kinetics. Table 1 and Table 2 show that the results from the classical kinetics have an important difference in relation to the model of the modified point kinetics that increases when the frequency is smaller.

References

- Duderstadt, J.J. and Hamilton, L.J. (1976) Nuclear Reactor Analysis. John Wiley & Sons Ltd., New York.
- Bell, G.I. and Glasstone (1970) Nuclear Reactor Theory. Van Nostrand Reinhold Ltd., New York.
- Henry, A.F. (1975) Nuclear Reactor Analysis. The MIT Press, Cambridge and London.
- Heizler, S.I. (2010) Asymptotic Telegrapher’s Equation (P1) Approximation for the Transport Equation. Nuclear Science and Engineering, 166, 17-35. http://physics.biu.ac.il/files/physics/shared/staff/u47/nse_166_17.pdf
- Espinosa-Paredes, G., Polo-Labarrios, M.A., Espinosa-Martinez, E.G. and del Valle-Gallegos, E. (2011) Fractional Neutron Point Kinetics Equations for Nuclear Reactor Dynamics. Annals of Nuclear Energy, 38, 307-330. http://www.sciencedirect.com/science/article/pii/S0306454910003816 http://dx.doi.org/10.1016/j.anucene.2010.10.012
- Akcasu, Z., Lellouche, G. and Shotkin, L.M. (1971) Mathematical Methods in Nuclear Reactor Dynamics. Academic Press, New York and London.
- Chao, Y.A. and Attard, A. (1985) A Resolution of the Stiffness Problem of Reactor Kinetics. Nuclear Science and Engineering, 90, 40-46. http://www.ans.org/pubs/journals/nse/a_17429 http://dx.doi.org/10.13182/NSE85-7
- Zhang, F., Chen, W.Z. and Gui, X.W. (2008) Analytic Method Study of Point-Reactor Kinetic Equation When Cold Start-Up. Annals of Nuclear Energy, 35, 746-749. http://www.sciencedirect.com/science/article/pii/S0306454907002368 http://dx.doi.org/10.1016/j.anucene.2007.08.015
- Hoogenboom, J.E. (1985) The Laplace Transformation of Adjoint Transport Equations. Annals of Nuclear Energy, 12, 151-152. http://www.sciencedirect.com/science/article/pii/030645498590091X http://dx.doi.org/10.1016/0306-4549(85)90091-X
- Fuchs, D. and Tabachnikov, S. (2000) Mathematical Omnibus: Thirty Lectures on Classic Mathematics. American Mathematical Society, Rhode Island.
- Hetrick, D.L. (1971) Dynamics of Nuclear Reactor. The University of Chicago Press Ltd., Chicago and London.
- Stacey, W.M. (2007) Nuclear Reactor Analysis. 2nd Edition, Wiley-VCH GmbH & CO KGaA, Weinheim.
- Kinard, M. and Allen, E.J. (2003) Efficient Numerical Solution of the Point Kinetics Equations in Nuclear Reactor Dynamics. Annals of Nuclear Energy, 31, 1039-1051. http://www.sciencedirect.com/science/article/pii/S0306454904000027 http://dx.doi.org/10.1016/j.anucene.2003.12.008
- Palma, D.A.P., Martinez, A.S. and Gonçalves, A.C. (2009) Analytical Solution of Point Kinetics Equations for Linear Reactivity Variation. Annals of Nuclear Energy, 36, 1469-1471. http://www.sciencedirect.com/science/article/pii/S0306454909001947 http://dx.doi.org/10.1016/j.anucene.2009.06.016
- Jahanbin, A. and Malmir, H. (2012) Kinetic Parameters Evaluation of PWRs Using Static Cell and Core Calculation Codes. Annals of Nuclear Energy, 41, 110-114. http://www.sciencedirect.com/science/article/pii/S0306454911004464