Geometrical and Material Buckling

In reactor physics, the concept of buckling is used to describe the relationship between requirements on fissile material inside a reactor core and dimensions and shape of that core. In general, criticality is achieved when the rate of neutron production is equal to the rate of neutron losses, including both neutron absorption and neutron leakage.

Material Buckling

In earlier text, the term Bm first appeared in the following equation:

material-buckling-equation

This parameter is known as the material buckling and it describes the characteristics of the fuel material in an infinite medium. For example let 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

Geometrical Buckling

geometric buckling of various shapes of reactors

Geometrical buckling and solutions of various shapes of reactors

Geometrical buckling is a measure of neutron leakage, while material buckling is a measure of neutron production minus absorption. With this terminology the criticality condition may also be stated as the material and geometric buckling being equal:

Bm = Bg

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.

The value of geometrical buckling for infinite slab reactor can be derived, when the vacuum boundary condition is applied on the solution of diffusion equation. The physically acceptable solution for infinite slab reactor is:

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

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

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

Criticality Condition

The basic classification of states of a reactor is according to the multiplication factor as eigenvalue which is a measure of the change in the fission neutron population from one neutron generation to the subsequent generation.

But these three basic states may be defined also according to the material and geometrical bucklings:

  • Bm < Bg. When a reactor is smaller (i.e. higher Bg and higher relative curveture) than the critical size for a given material, Bm < Bg, then the reactor is subcritical.
  • Bm = Bg. When a reactor size matches the critical size for a given material, Bm = Bg, then the reactor is critical.
  • Bm > Bg. When a reactor is larger than the critical size for a given material, Bm > Bg, then the reactor is supercritical.

It must be added, for any positive value of materials buckling, there is a unique critical size for each reactor geometry. For reactors of shapes other than spheres the geometrical buckling takes the form Bg = C/R, where the coefficient C is determined by the solution of diffusion equation with vacuum boundary condition. R is a characteristic dimension. Generally, the multiplication of a uniform reactor of any shape and size is given by keff = k.PNL, with the non/leakage probability written as (for large reactors):

fast non-leakage probability_4

where M2 is the migration area and the subscript is dropped from B, the geometric buckling. As can be seen the total non-leakage probability of large reactors is primarily a function of migration area and the relative curvature of the neutron flux given by the geometrical buckling.

Example: Calculate the geometrical buckling

Let assume a bare square cylinder (i.e., height = diameter; H = 2R). Assume that the material buckling of this reactor, which is given using one-group cross sections, is:

material-buckling-equation2

Calculate the critical radius (Bm = Bg) using one-group diffusion theory.

The geometrical buckling with extrapolated distance is:

geometric-buckling-equation

From this equation we can get Re= 30.5 cm.

Geometric Buckling of Reflected Reactor

Adding a reflector to a reactor allows either the volume of the reactor or the requirements on fuel to be reduced, or some combination of the two. If the reflector savings is known, the calculation of the critical dimensions of a reflected reactor needs only the solution of the bare reactor, which is simpler problem. For example, for the cylindrical reactor, it is only necessary to determine the bare critical radius R0 and the reflected radius is simply R = R0 – δ, where R0  is critical diameter of a bare reactor. The geometrical buckling of infinite cylindrical reactor will then be:

geometrical-buckling-savings

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: