Diffusion Equation – Infinite Slab Reactor

Solutions of the Diffusion Equation – Multiplying Systems

In previous section it has been considered that the environment is non-multiplying. In non-multiplying environment neutrons are emitted by a neutron source situated in the center of coordinate system and then they freely diffuse through media. We are now prepared to consider neutron diffusion in multiplying system, which contains fissionable nuclei (i.e. Σf ≠ 0). In this multiplying system we will also study spatial distribution of neutrons, but in contrast to non-multiplying environment these neutrons can trigger nuclear fission reaction.
Multiplying systems from criticality point of view
The required condition for a stable, self-sustained fission chain reaction in a multiplying system (in a nuclear reactor) is that exactly every fission initiate another fission. The minimum condition is for each nucleus undergoing fission to produce, on the 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

It is obvious the infinite multiplication factor in a multiplying system is a measure of 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, then the number of neutrons is decreasing 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.
Reactor criticality
Reactor criticality. A – a supercritical state; B – a critical state; C – a subcritical state
Diffusion Theory - Multiplying Systems
In this section, we will solve the following diffusion equation

diffusion equation - multiplying system

in various geometries that satisfy the boundary conditions. In this equation ν is number of neutrons emitted in fission and  Σf is macroscopic cross-section of fission reaction. Ф denotes a reaction rate. For example a fission of 235U by thermal neutron yields 2.43 neutrons.

It must be noted that we will solve the diffusion equation without any external source. This is very important, because such equation is a linear homogeneous equation in the flux. Therefore if we find one solution of the equation, then any multiple is also a solution. This means that the absolute value of the neutron flux cannot possibly be deduced from the diffusion equation. This is totally different from problems with external sources, which determine the absolute value of the neutron flux.

We will start with simple systems (planar) and increase complexity gradually. The most important assumption is that all neutrons are lumped into a single energy group, they are emitted and diffuse at thermal energy (0.025 eV). Solutions of diffusion equations in this case provides an illustrative insights, how can be the neutron flux distributed in a reactor core.

Solution for the Infinite Slab Reactor

Diffusion Theory - Slab Reactor-minLet assume a uniform reactor (multiplying system) in the shape of a slab of physical width a in the x-direction and infinite in the y- and z-directions. This reactor is situated in the center at x=0. In this geometry the flux does not vary in y and z allowing us to eliminate the y and z derivatives from ∇2. The flux is then a function of x only, and therefore the Laplacian and diffusion equation can be written as:

solution of diffusion equation - infinite slab

The quantity Bg2 is called the geometrical buckling of the reactor and depends only on the geometry. This term is derived from the notion that the neutron flux distribution is somehow ‘‘buckled’’ in a homogeneous finite reactor. It can be derived the geometrical buckling is the negative relative curvature of the neutron flux (Bg2 = ∇2Ф(x) / Ф(x)). In a small reactor the neutron flux have more concave downward or ‘‘buckled’’ curvature (higher Bg2) than in a large one. This is a very important parameter and it will be discussed in following sections.

For x > 0, this diffusion equation has two possible solutions sin(Bgx) and cos(Bgx), which give a general solution:

Φ(x) = A.sin(Bgx) + C.cos(Bg x)

From finite flux condition (0≤ Φ(x) < ∞), that required only reasonable values for the flux, it can be derived, that A must be equal to zero. The term sin(Bgx) goes to negative values as x goes to negative values and therefore it cannot be part of a physically acceptable solution. The physically acceptable solution must then be:

Φ(x) =  C.cos(Bg x)

where Bg can be determined from vacuum boundary condition.

extrapolated length - boundary condition
Neutron flux as a function of position near a free surface according to diffusion theory and transport theory.

The vacuum boundary condition requires the relative neutron flux near the boundary to have a slope of -1/d, i.e., the flux would extrapolate linearly to 0 at a distance d beyond the boundary. This zero flux boundary condition is more straightforward and is can be written mathematically as:

infinite slab reactor - boundary condition

If d is not negligible, physical dimensions of the reactor are increased by d and extrapolated boundary is formulated with dimension ae/2 = a/2 + d and this condition can be written as Φ(a/2 + d) =  Φ(ae/2) = 0.

Therefore, the solution must be  Φ(ae/2) =  C.cos(Bg .ae/2) = 0 and the values of geometrical buckling, Bg, are limited to Bg = /a_e, where n is any odd integer. The only one physically acceptable odd integer is n=1, because higher values of n would give cosine functions which would become negative for some values of x. The solution of the diffusion equation is:

infinite slab reactor - solution

It must be added the constant C cannot be obtained from this diffusion equation, because this constant gives the absolute value of neutron flux. In fact the neutron flux can have any value and the critical reactor can operate at any power level. It should be noted the cosine flux shape is only in theoretical case in a uniform homogeneous reactor at low power levels (at “zero power criticality”).

In power reactor core (at full power operation), the neutron flux can reach, for example, about 3.11 x 1013 neutrons.cm-2.s-1, but this values depends significantly on many parameters (type of fuel, fuel burnup, fuel enrichment, position in fuel pattern, etc.).

The power level does not influence the criticality (keff) of a power reactor unless thermal reactivity feedbacks act (operation of a power reactor without reactivity feedbacks is between 10E-8% – 1% of rated power).

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 above:

Diffusion Theory