Facebook Instagram Youtube Twitter

Nuclear Fission Chain Reaction

A nuclear fission chain reaction is a self-propagating sequence of fission reactions in which neutrons released in fission produce additional fission in at least one other nucleus. The chain reaction can take place only in the proper multiplication environment and only under proper conditions.

The fission process may produce 2, 3, or more free neutrons that are capable of inducing further fissions and so on. This sequence of fission events is known as the fission chain reaction, and it is important in nuclear reactor physics.

The chain reaction can take place only in the proper multiplication environment and only under proper conditions. Suppose one neutron causes two further fissions. In that case, the number of neutrons in the multiplication system will increase in time, and the reactor power (reaction rate) will also increase in time. To stabilize such a multiplication environment, it is necessary to increase the non-fission neutron absorption in the system (e.g.,, to insert control rods). Moreover, this multiplication environment (the nuclear reactor) behaves like an exponential system, which means the power increase is not linear, but it is exponential.

Nuclear chain reaction
The nuclear chain reaction occurs when one single nuclear reaction causes an average of one or more subsequent nuclear reactions.

On the other hand, if one neutron causes less than one further fission, the number of neutrons in the multiplication system will decrease in time, and the reactor power (reaction rate) will also decrease in time. It is necessary to decrease the non-fission neutron absorption in the system (e.g.,, to withdraw control rods) to sustain the chain reaction.

There is always a competition for the fission neutrons in the multiplication environment. Some neutrons will cause further fission reaction, some will be captured by fuel or non-fuel materials, and some will leak out of the system.

It is necessary to define the infinite and finite multiplication factors of a reactor to describe the multiplication system. The method of calculations of multiplication factors was developed in the early years of nuclear energy. It is only applicable to thermal reactors, where the bulk of fission reactions occurs at thermal energies. This method puts into context all the processes associated with the thermal reactors (e.g.,, neutron thermalization, neutron diffusion, or fast fission) because the most important neutron-physical processes occur in energy regions that can be clearly separated from each other. In short, the calculation of the multiplication factor gives a good insight into the processes that occur in each thermal multiplying system.

Fast vs. Thermal Flux Spectrum
thermal vs. fast reactor neutron spectrum
The spectrum of neutron energies produced by fission varies significantly with certain reactor designs. thermal vs. fast reactor neutron spectrum
Six Factor Formula - Fast Reactors
For fast reactors, in which neutrons cause the fission with a very broad energy distribution, such an analysis is inappropriate. The neutron flux in fast reactors has to be divided into many energy groups. Moreover, in fast reactors, neutron thermalization is an undesirable process, and therefore the four-factor formula does not make any sense. The resonance escape probability is insignificant because very few neutrons exist at energies where resonance absorption is significant. The thermal non-leakage probability does not exist because the reactor is designed to avoid the thermalization of neutrons.

Infinite Multiplication Factor – Four Factor Formula

In this section, the infinite multiplication factor, which describes all the possible events in the life of a neutron and effectively describes the state of an infinite multiplying system, will be defined.

The necessary condition for a stable, self-sustained fission chain reaction in a multiplying system (in a nuclear reactor) is that exactly every fission initiates another fission. The minimum condition is for each nucleus undergoing fission to produce, on average, at least one neutron that causes fission of another nucleus. Also, the number of fissions occurring per unit time (the reaction rate) within the system must be constant.

This condition can be expressed conveniently in terms of the multiplication factor. The infinite multiplication factor is the ratio of the neutrons produced by fission in one neutron generation to the number of neutrons lost through absorption in the preceding neutron generation. This can be expressed mathematically, as shown below.

Multiplication Factor

The infinite multiplication factor in a multiplying system measures the change in the fission neutron population from one neutron generation to the subsequent generation.

  • k < 1. If the multiplication factor for a multiplying system is less than 1.0. In this case, the number of neutrons decreases in time (with the mean generation time), and the chain reaction will never be self-sustaining. This condition is known as the subcritical state.
  • k = 1. If the multiplication factor for a multiplying system is equal to 1.0, then there is no change in neutron population in time, and the chain reaction will be self-sustaining. This condition is known as the critical state.
  • k > 1. If the multiplication factor for a multiplying system is greater than 1.0, then the multiplying system produces more neutrons than are needed to be self-sustaining. The number of neutrons is exponentially increasing in time (with the mean generation time). This condition is known as the supercritical state.
Example - Neutron Generation and Neutron Population
The number of neutrons (the neutron population) in the core at time zero is 1000 and k = 1.001 (~100 pcm).

Calculate the number of neutrons after 100 generations. Let say the mean generation time is ~0.1s.

Solution:
To calculate the neutron population after 100 neutron generations, we use following equation:

Nn=N0. (k)n
N1=N0.1.001 = 1001 neutrons after one generation
N2=N0.1.001.1.001 = 1002 neutrons after two generations
N3=N0.1.001.1.001.1.001 = 1003 neutrons after three generations
.
.

N50=N0. (k)50 = 1051 neutrons after fifty generations.

.

.

N100=N0. (k)100 = 1105 neutrons after hundred generations.

If we consider the mean generation time to be ~0.1s, the increase from 1000 neutrons to 1105 neutrons occurs within 10 seconds.

See also: Neutron Generation – Neutron Population

Effective Multiplication Factor in Reactor Kinetics
The simplest equation governing the neutron kinetics of the system with delayed neutrons is the point kinetics equation. This equation states that the time change of the neutron population is equal to the excess of neutron production (by fission) minus neutron loss by absorption in one mean generation time with delayed neutrons (ld). The role of ld is evident. Longer lifetimes give simply slower responses to multiplying systems.

If there are neutrons in the system at t=0, that is, if n(0) > 0, the solution of this equation gives the simplest point kinetics equation with delayed neutrons (similarly to the case without delayed neutrons):point kinetics equation with delayed neutronsLet us consider that the mean generation time with delayed neutrons is ~0.085 and k (k – neutron multiplication factor) will be step increased by only 0.01% (i.e., 10pcm or ~1.5 cents), that is k=1.0000 will increase to k=1.0001.

It must be noted such reactivity insertion (10pcm) is very small in the case of LWRs. The reactivity insertions of the order of one pcm are for LWRs practically unrealizable. In this case, the reactor period will be:

T = ld / (k-1) = 0.085 / (1.0001-1) = 850s

This is a very long period. In ~14 minutes, the reactor’s neutron flux (and power) would increase by a factor of e = 2.718. This is a completely different dimension of the response on reactivity insertion compared to the case without delayed neutrons, where the reactor period was 1 second.

Reactors with such kinetics would be quite easy to control. From this point of view, it may seem that reactor control will be quite a boring affair. It will not! The presence of delayed neutrons entails many specific phenomena that will be described in later chapters.

Interactive Chart - Reactor Kinetics
Press the “clear and run” button and try to increase the power of the reactor.

Compare the response of the reactor with the case of Infinite Multiplying System Without Source and without Delayed Neutrons (or set the β = 0).

But the infinite multiplication factor can also be defined in terms of the most important neutron-physical processes in the nuclear reactor. For simplicity, we will first consider a multiplying system that is infinitely large and therefore has no neutron leakage. In the infinite system. Four factors are completely independent of the size and shape of the reactor that gives the inherent multiplication ability of the fuel and moderator materials without regard to leakage:

Fast Fission Factor
The fast fission process is in the multiplication factor characterized by the fast fission factor, ε, which increases the fast neutron population in one neutron generation. The fast fission factor is the ratio of the fast neutrons produced by fissions at all energies to the number of fast neutrons produced in thermal fission.

See also: Fast Fission Factor

Resonance Escape Probability
The resonance escape probability, symbolized by p, is the probability that a neutron will be slowed to thermal energy and will escape resonance capture. This probability is defined as the ratio of the number of neutrons that reach thermal energies to the number of fast neutrons that slow down.

See also: Resonance Escape Probability

Thermal Utilisation Factor
The thermal utilization factor, f, is the fraction of the thermal neutrons that are absorbed in the nuclear fuel, in all isotopes of the nuclear fuel. It
describes how effectively (how well they are utilized) are thermal neutrons absorbed in the fuel.

The value of the thermal utilization factor is given by the ratio of the number of thermal neutrons absorbed in the fuel (all nuclides) to the number of thermal neutrons absorbed in all the material that makes up the core.

See also: Thermal Utilisation Factor

Reproduction Factor
The number of neutrons created in the new generation is determined by the neutron reproduction factor. The reproduction factor, η, is defined as the ratio of the number of fast neutrons produced by thermal fission to the number
of thermal neutrons absorbed in the fuel.

See also: Reproduction Factor

The infinite multiplication factor (k) may be expressed mathematically in terms of these factors by the following equation, usually known as the four-factor formula:

k = η.ε.p.f

In reactor physics, k or its finite form keff is the most significant parameter about reactor control. At any specific power level or condition of the reactor, keff is kept as near
to the value of 1.0 as possible. At this point in the operation, the neutron balance is kept to exactly one neutron completing the life cycle for each original neutron absorbed in the fuel.

From infinite to finite multiplication factor

The infinite multiplication factor is derived based on the assumption that no neutrons leak out of the reactor (i.e., a reactor is infinitely large). But in reality, each nuclear reactor is finite, and neutrons can leak out of the reactor core. The multiplication factor that takes neutron leakage into account is the effective multiplication factorkeff, which is the ratio of the neutrons produced by fission in one neutron generation to the number of neutrons lost through absorption and leakage in the preceding neutron generation.

The effective multiplication factor (keff) may be expressed mathematically in terms of the infinite multiplication factor (k) and two additional factors which account for neutron leakage during neutron thermalization (fast non-leakage probability) and neutron leakage during neutron diffusion (thermal non-leakage probability) by following equation, usually known as the six-factor formula:

keff = k . Pf . Pt

Nuclear Fission Chain Reaction
Neutron Life Cycle with keff = 1

Operational factors that affect the fission chain reaction in PWRs.

Detailed knowledge of all possible operational factors that may affect the multiplication factor of the system is of importance in reactor control. It was stated the keff is during reactor operation kept as near to the value of 1.0 as possible. Many factors influence the criticality of the reactor. For illustration, in an extreme case, the presence of humans (due to the water, carbon, which are good neutron moderators) near fresh uranium fuel assembly influences the multiplication properties of the assembly.
If any operational factor changes one of the contributing factors to keff (keff = η.ε.p.f.Pf.Pt), the ratio of 1.0 is not maintained, and this change in keff makes the reactor either subcritical or supercritical. Some examples of these operational changes that may take place in PWRs, are below and are described below:

Change in the control rods position
↓control rods ⇒ ↓keff = η.ε.p.  ↓f  .Pf.Pt

Control rods (insertion/withdrawal) influence the thermal utilization factor. For example, control rods insertion causes the addition of new absorbing material into the core, and this causes a decrease in thermal utilization.

The thermal utilization factor for heterogeneous reactor cores must be calculated in terms of reaction rates and volumes, for example, by the following equation:

thermal utilisation factor - equation2

where Σa is the macroscopic absorption cross-section, which is the sum of the capture cross-section and the fission cross-section, Σa = Σc + Σf. The superscripts U, M, P, CR, B, BA, and O, refer to uranium fuel, moderator, poisons, control rods, boric acid, burnable absorbers, etc. The presence of control rods, boric acid, or poisons causes a decrease in neutron utilization, which, in turn, causes a decrease in the multiplication factor.

Compared with the chemical shim, which offset positive reactivity excess in the entire core, with control rods, the unevenness of neutron-flux density in the reactor core may arise because they act locally.

Change in the boron concentration

↑boron ⇒ ↓keff = η.ε.p.  ↓f  .Pf.Pt

The concentration of boric acid diluted in the primary coolant influences the thermal utilization factor. For example, an increase in the concentration of boric acid (chemical shim) causes the addition of new absorbing material into the core, and this causes a decrease in the thermal utilization factor. 

The thermal utilization factor for heterogeneous reactor cores must be calculated in terms of reaction rates and volumes, for example, by the following equation:

thermal utilisation factor - equation2

where Σa is the macroscopic absorption cross-section, which is the sum of the capture cross-section and the fission cross-section, Σa = Σc + Σf. The superscripts U, M, P, CR, B, BA, and O, refer to uranium fuel, moderator, poisons, control rods, boric acid, burnable absorbers, etc. It is obvious that the presence of control rods, boric acid, or poisons causes a decrease in neutron utilization, which, in turn, causes a decrease in multiplication factor.

Compared with burnable absorbers (long-term reactivity control) or with control rods (rapid reactivity control), the boric acid avoids the unevenness of neutron-flux density in the reactor core because it is dissolved homogeneously in the coolant in the entire reactor core. On the other hand, high concentrations of boric acid may lead to a positive moderator temperature coefficient, which is undesirable. In this case, more burnable absorbers must be used.

Moreover, this method is slow in controlling reactivity. Normally, it takes several minutes to change the boric acid concentration (dilute or borate) in the primary loop. For rapid changes of reactivity, control rods must be used.

Change in the moderator temperature
↑TM ⇒ ↓keff = η.ε.  ↓p  . ↑f .  ↓Pf  .  ↓P (BOC)

↑TM ⇒ ↓keff = η.ε.  ↓p  .f.  ↓Pf  .  ↓P (EOC)

This operational change is very difficult to describe because changes in moderator temperature lead to almost all the coefficients. Major impacts on the multiplication of the system arise from the change of the resonance escape probability and the change of total neutron leakage (see thermal non-leakage probability and fast non-leakage probability).

  • moderator-to-fuel ratioChange of the resonance escape probability. It is known, the resonance escape probability is also dependent on the moderator-to-fuel ratio. All PWRs are designed as under moderated reactors. As the moderator temperature increases, the ratio of the moderating atoms (molecules of water) decreases due to the thermal expansion of water. Its density simply decreases. This, in turn, causes hardening of neutron spectrum in the reactor core resulting in higher resonance absorption (lower p). The decreasing density of the moderator causes that neutrons stay at a higher energy for a longer period, which increases the probability of non-fission capture of these neutrons. This process is one of two processes, which determine the moderator temperature coefficient (MTC). The second process is connected with the leakage probability of the neutrons. The moderator temperature coefficient must be for most PWRs negative, which improves the reactor stability because reactor core heating causes a negative reactivity insertion.
  • Boron letdown curve (chemical shim) and boron 10 depletion
    Boron letdown curve (chemical shim) and boron 10 depletion during a 12-month fuel cycle.

    Change of the thermal utilization factor. The impact on the thermal utilization factor depends strongly on the amount of boron that is diluted in the primary coolant (chemical shim). As the moderator temperature increases, the density of water decreases due to the thermal expansion of water. But along with the moderator also boric acid is expanded out of the core. Since boric acid is a neutron poison, expanding out of the core, positive reactivity is added. The positive reactivity addition due to the expansion of boron out of the core offsets the negative reactivity addition due to the expansion of the moderator out of the core. This effect is significant at the beginning of the cycle (BOC) and gradually loses significance as the boron concentration decreases.

  • Change of the neutron leakage. Since both (Pf and Pt) are affected by a change in moderator temperature in a heterogeneous water-moderated reactor and the directions of the feedbacks are the same, the resulting total non-leakage probability is also sensitive to the change in the moderator temperature. As a result, an increase in the moderator temperature causes that the probability of leakage to increase. In the case of the fast neutron leakage, the moderator temperature influences macroscopic cross-sections for elastic scattering reactionss.NH2O) due to the thermal expansion of water, which increases the moderation length. This, in turn, causes an increase in the leakage of fast neutrons.
    • For the thermal neutron leakage, there are two effects. Both processes have the same direction and together cause the increase in the thermal neutron leakage. This physical process is a part of the moderator temperature coefficient (MTC).

Neutron Moderators - Parameters

Change in the fuel temperature
↑Tf ⇒ ↓keff = η.ε.  ↓p  .f.Pf  .Pt

Change in the fuel temperature affects primarily the resonance escape probability, which is connected with the phenomenon usually known as the Doppler broadening. The Doppler effect is generally considered to be the most important effect, which improves the reactor stability. Especially in the case of reactivity-initiated accidents (RIA), the Doppler coefficient of reactivity would be the first in the compensation of the inserted positive reactivity. The time for heat to be transferred to the moderator is usually measured in seconds, while the Doppler coefficient is effective almost instantaneously. The Doppler broadening with the self-shielding process causes the Doppler coefficient (or the fuel temperature coefficient) for all power reactors to always be negative. Therefore an increase in the fuel temperature promptly causes an increase in the resonance integral (Ieff), which, in turn, causes a negative reactivity insertion. It is of the highest importance in reactor safety.

See also: Doppler Broadening.
See also: Self-shielding

Change in the pressure
↓pressure ⇒ ↓keff = η.ε.  ↓p  .f.  ↓Pf  .  ↓Pt

Although water is considered incompressible, in reality, it is slightly compressible (especially at 325°C (617°F)). The effect of pressure in the primary circuit has similar consequences as the moderator temperature. Compared with the effects of moderator temperature changes, changes in pressure have a lower-order impact on reactivity. The causes are only in the density of the moderator, not in the change of microscopic cross-sections.

The pressure coefficient of reactivity has a slightly positive effect on reactivity as the system’s pressure is increased if. At high boron concentrations, the pressure coefficient may reach negative values, but for many PWRs, it is prohibited to operate under such conditions. Therefore burnable absorbers are usually added to the fuel. They lower the initial boron concentration.

Note: Effects of the nuclate boiling of the primary coolant are not discussed here.

Change in the coolant flow rate
↓flow rate ⇒ ↑TM(average) ⇒ ↓keff = η.ε.  ↓p  .f.  ↓Pf  .  ↓Pt

The effect of change in the flow rate through the primary circuit has identical consequences as the effects of the moderator temperature. In reality, when there is an abrupt change (e.g.,, due to a disconnection of the reactor coolant pump) in the flow rate and the reactor power remains the same, the difference between inlet and outlet temperatures must increase. It follows from the basic energy equation of reactor coolant, which is below:

P=↓ṁ.c.↑∆t

flow rate decreaseThe pressure determines the inlet temperature in the steam generators. Therefore the inlet temperature changes minimally during the transient. It follows the outlet temperature must change significantly as the flow rate changes. When the inlet temperature remains almost the same and the outlet changes significantly, it stands to reason, the average temperature of coolant (moderator) will also change significantly. Therefore the effect of change in the flow rate through the primary circuit has identical consequences as the effects of the moderator temperature.

The decrease in flow rate is associated with negative reactivity insertion. Special attention is needed in case of an abrupt increase in the flow rate. At normal operation, such an increase in the flow rate can not occur, except for the controlled reactor coolant pump connection, which can be connected only under specific conditions.

Presence of boiling of the coolant
boiling ⇒ ↓keff = η.ε.  ↓p  .f.Pf  .Pt

In pressurized water reactors, nucleate boiling may occur even during operational conditions. Nucleate boiling occurs when any surface of fuel cladding reaches the saturation temperature (e.g.,, 350°C), which is determined by the pressure in the pressurizer (e.g.,, 16MPa). Such local nucleate boiling does not pose any problem for the reactor operation.

On the other hand, during abnormal conditions, boiling in the reactor core is one of the most important phenomena that may take place in the core. From the reactivity point of view, nucleate boiling has very important consequences on the reactivity of the reactor core. Boiling affects reactivity in the same manner as voids, and therefore it is characterized by the void coefficient.

The formation of voids in the core has the same effect as the change in the moderator temperature (change in the density of the moderator). Compared with the change in the moderator temperature, boiling minimally affects the neutron leakage. It is unlikely that local boiling occurs at the periphery of the reactor core, where the local power drops significantly.

Presence of burnable absorbers

↑burnable absorbers ⇒ ↓keff = η.ε.p.  ↓f  .Pf.Pt

The number of burnable absorbers in the nuclear fuel influences the thermal utilization factor. In some cases (especially in the case of gadolinium absorbers), the presence of burnable absorbers influences all the factors in the four-factor formula due to very high self-shielding effects. But at this place, we consider only the change in the thermal utilization factor. An increase in the number of burnable absorbers causes the addition of new absorbing material into the core, and this causes a decrease in the thermal utilization factor. 

The thermal utilization factor for heterogeneous reactor cores must be calculated in terms of reaction rates and volumes, for example, by the following equation:

thermal utilisation factor - equation2

effect of gadolinium absorbers
The effect of gadolinium burnable absorbers (BA) can be demonstrated on boron letdown curves. At the beginning of a specific fuel cycle, the critical concentration of boric acid in the reactor core without burnable absorbers (blue curve) significantly differs from the critical concentration of boric acid in the reactor core with burnable absorbers (red curve). The difference is dependent on the amount of BA used.

where Σa is the macroscopic absorption cross-section, which is the sum of the capture cross-section and the fission cross-section, Σa = Σc + Σf. The superscripts U, M, P, CR, B, BA, and O, refer to uranium fuel, moderator, poisons, control rods, boric acid, burnable absorbers, etc. The presence of control rods, boric acid, or burnable absorbers causes a decrease in neutron utilization, which, in turn, causes a decrease of multiplication factors.

Fuel burnup
↑burnup ⇒ ↓keff = ↓η  .ε.p.  ↓f  .  ↑↓Pf.  ↑↓Pt

It is hard to describe the effects of fuel burnup on the six-factor formula. It must be noted the criticality must be maintained for a long period. Therefore all the negative effects must be compensated by increasing the thermal utilization factor (boron dilution or compensating rods withdrawal).

Thermal Utilization Factor

The thermal utilization factor slightly changes with the fuel burnup. The fresh fuel at the beginning of the cycle comprises only the absorption by the 235U. As the amount of 239Pu and other higher transuranic elements increases because of the radiative capture of a neutron by the 238U in the core, it is necessary to consider the change of fuel composition in determining the value of f at different times of the fuel cycle.

In general, the thermal utilization factor decreases in time as the total content of fissile isotopes decreases and the total content of neutron poisons (fission products with high absorption cross-sections) increases. But in the power reactors, in which the criticality must be maintained for a long period (e.g.,, 12-month or up to 24-month) without refueling, the thermal utilization factor may not decrease. Such a decrease would imply an inevitable reactor shutdown. The continuous decrease in ΣaU must be offset by the continuous decrease in ΣaB, which means the concentration of boric acid (in the case of PWRs) must be continuously decreased as the fuel loses its reactivity (kinf). For reactors in which the chemical shim can not be used, the excess of reactivity is compensated by compensating rods.

Reproduction Factor

There is essentially a small change in η over the lifetime of the reactor core (decreases). This is because there is a continuous decrease in ΣfU, but on the other hand, this decrease is partially offset by the increase in ΣfPu. As the fuel burnup increases, the 239Pu begins to contribute to the neutron economy of the core.

See also: Nuclear Breeding

Neutron Leakage

In power reactors, the total non-leakage probability also significantly changes with fuel burnup. This dependency is not associated with any of the parameters like the diffusion coefficient or the geometrical buckling. In power reactors, the total non-leakage probability depends on the certain fuel loading pattern, and the reload strategy. Neutron leakage is one of the key parameters in the neutron and fuel economy.

To enhance the neutron and fuel economy, core designers design the low leakage loading patterns. Fresh fuel assemblies are not situated in the peripheral positions of the reactor core. The peripheral positions are loaded with the fuel with the highest fuel burnup. Compared to the average assemblies, these “high” burnup assemblies have inherently lower relative power (due to the lower kinf and the fact they feel the presence of a non-multiplying environment). In short, this parameter is significantly dependent on a certain loading pattern. During fuel burnup, the neutron leakage usually increases, especially in low leakage loading patterns. This process is caused by reducing the differences in kinf between fresh fuel assemblies and peripheral high-burnup assemblies.

 
References:
Nuclear and Reactor Physics:
  1. J. R. Lamarsh, Introduction to Nuclear Reactor Theory, 2nd ed., Addison-Wesley, Reading, MA (1983).
  2. J. R. Lamarsh, A. J. Baratta, Introduction to Nuclear Engineering, 3d ed., Prentice-Hall, 2001, ISBN: 0-201-82498-1.
  3. W. M. Stacey, Nuclear Reactor Physics, John Wiley & Sons, 2001, ISBN: 0- 471-39127-1.
  4. Glasstone, Sesonske. Nuclear Reactor Engineering: Reactor Systems Engineering, Springer; 4th edition, 1994, ISBN: 978-0412985317
  5. W.S.C. Williams. Nuclear and Particle Physics. Clarendon Press; 1 edition, 1991, ISBN: 978-0198520467
  6. G.R.Keepin. Physics of Nuclear Kinetics. Addison-Wesley Pub. Co; 1st edition, 1965
  7. Robert Reed Burn, Introduction to Nuclear Reactor Operation, 1988.
  8. U.S. Department of Energy, Nuclear Physics and Reactor Theory. DOE Fundamentals Handbook, Volume 1 and 2. January 1993.

Advanced Reactor Physics:

  1. K. O. Ott, W. A. Bezella, Introductory Nuclear Reactor Statics, American Nuclear Society, Revised edition (1989), 1989, ISBN: 0-894-48033-2.
  2. K. O. Ott, R. J. Neuhold, Introductory Nuclear Reactor Dynamics, American Nuclear Society, 1985, ISBN: 0-894-48029-4.
  3. D. L. Hetrick, Dynamics of Nuclear Reactors, American Nuclear Society, 1993, ISBN: 0-894-48453-2. 
  4. E. E. Lewis, W. F. Miller, Computational Methods of Neutron Transport, American Nuclear Society, 1993, ISBN: 0-894-48452-4.

See previous:

See above:

Neutron Reactions

See next: