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.
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 [
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 [
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.
The theory of neutron transport is the wide model to describe neutron distribution in a nuclear reactor. It is described in [
and
where
and
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
Equation (3) can be re-written thus:
We apply the operator
and
it results that:
and
where the operators
and
In replacing Equation (8) in Equation (1), multiplying the resulting equation by
Considering the approximation described in [
and defining the transport cross section:
from Equation (14), it is possible to write:
There is a similarity of Equation (15) with the Telegrapher Equation, as can be seen in [
In dividing Equation (15) by the transport cross section, using the definition of diffusion coefficient, applying
the diverging operator and disregarding the term
where the diffusion coefficient is defined as:
After that, in deriving Equation (9) in relation to time, and multiplying by
very Equation (9) and replacing Equation (16), one obtains:
In the Equations (10) and (17), when in a stationary regime, all of its time derivatives are disregarded. Then, it results that:
Considering the adjoint flux of neutrons from Equation (18), this adjoint equation is:
where,
and
Equation (18) is integrated in the volume and in the energy, and re-written:
In multiplying Equation (17) by the adjoint flux of neutrons and integrating in the volume and in the energy
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
In writing the neutron flux as the product of the two factors in Equation (26), the intent is that the amplitude factor,
In defining the integral:
one obtains, after replacing Equation (26) in Equation (25) and the division of the result by
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:
and
In multiplying Equation (10) by the adjoint flux of neutrons and after that integrating in the volume and in the energy
Dividing Equation (36) by
where
Equations (28) and (37) form the new model for the point kinetics, called modified point kinetics.
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),
From there on, some terms of the first and second derivative of
The new 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:
Note that appears in the definition of Equation (29) into Equation (40).
Therefore,
In Equation (34) replace the operators
Considering the homogeneous medium rewrite the Equation (43):
Replacing the Equations (27), (29) and (42) in Equation (44):
once:
follows:
The solution for point kinetics equations can be obtained in several ways, as in [
and
The initial conditions are:
and
Applying the conditions (50), (51) and (52) in Equation (48) it is possible to write:
One therefore chooses
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:
It adds to the Equation (48) with Equation (49):
You can also replace
The Equation (57) is derived in the time. Then, it results that:
Multiplying the Equation (56) by
It adds to the Equation (58) with the Equation (59):
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:
In applying the Laplace transform [
being,
and
In Equation (62) we isolate the Laplace transform of the neutron density:
Applying the reverse Laplace transform on both sides of Equation (65) one have that:
To solve the reverse Laplace transform one should factor the polynomials so to obtain:
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:
and
The Equation (73) can be solved directly to obtain the value of
gree. Initially it replaces the variable
Note that in Equation (74) the term corresponding to
The Equation (75) was resolved according [
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
therefore,
Follows:
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
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., [
The values considered in
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
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
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
Parameter | Symbol | Value |
---|---|---|
Decay constant | λ | 0.0810958 s−1 |
Mean generation time | Λ | 0.002 s |
Absorption cross section | Σa | 0.14 cm−1 |
Diffusion coefficient | D | 10 cm |
Neutron velocity | v | 3 × 106 cm/s |
Fraction of delayed neutrons | β | 0.007 |
Absorption frequency | fA | 4.167 s−1 |
Model and Reactivity | t = 0.4 s | t = 1 s | t = 10 s | t = 20 s | t = 40 s | t = 100 s |
---|---|---|---|---|---|---|
Classical kinetics ρ = 0.001 | 1.1173 | 1.1648 | 1.3186 | 1.5033 | 1.9537 | 4.2888 |
Modified Kinetics ρ = 0.001 | 1.0979 | 1.1528 | 1.3129 | 1.4946 | 1.9368 | 4.2151 |
Classical kinetics ρ = 0.003 | 1.4161 | 1.6770 | 2.9278 | 5.1718 | 16.138 | 490.25 |
Modified Kinetics ρ = 0.003 | 1.3341 | 1.5954 | 2.8417 | 4.9561 | 15.076 | 424.36 |
Classical kinetics ρ = 0.006 | 2.0958 | 3.4712 | 56.554 | 947.92 | 2.663 × 105 | 5.905 × 1012 |
Modified Kinetics ρ = 0.006 | 1.8212 | 2.9008 | 38.432 | 487.09 | 78236 | 3.242 × 1011 |
Model and Reactivity | t = 0.4 s | t = 1 s | t = 10 s | t = 20 s | t = 40 s | t = 100 s |
---|---|---|---|---|---|---|
Classical kinetics ρ = 0.001 | 1.1173 | 1.1648 | 1.3186 | 1.5033 | 1.9537 | 4.2888 |
Modified kinetics ρ = 0.001 | 1.0379 | 1.0798 | 1.2688 | 1.4279 | 1.8083 | 3.6727 |
Classical kinetics ρ = 0.003 | 1.4161 | 1.6770 | 2.9278 | 5.1718 | 16.138 | 490.25 |
Modified kinetics ρ = 0.003 | 1.1185 | 1.2640 | 2.3225 | 3.7127 | 9.4734 | 157.38 |
Classical kinetics ρ = 0.006 | 2.0958 | 3.4712 | 56.554 | 947.92 | 2.663 × 105 | 5.905 × 1012 |
Modified kinetics ρ = 0.006 | 1.2524 | 1.6157 | 9.1475 | 42.311 | 898.86 | 8.617 × 106 |
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.
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.