In a nuclear reactor, criticality is achieved when the rate of neutron production is equal to the rate of neutron losses, including both neutron absorption and neutron leakage. Geometric buckling is a measure of neutron leakage, while material buckling is a measure of neutron production minus absorption. Thus, in the simplest case of a bare, homogeneous, steady state reactor, the geometric and material buckling must be equal.
Both buckling terms are derived from the diffusion equation:
S+DΔ2Φ-ΣaΦ = 0
where S is the source, and from diffusion theory D=1/3Σtr and L2=D/Σa. For thermal neutrons the source is:
S=k∞lfΣaΦ.
where k∞ is from the four factor formula and lf is the fast neutron non-leakage probability. Rearranging the diffusion equation becomes
-Δ2Φ/Φ=(k∞lf-1)/L2
The left side of the equation is the geometric buckling and the right side is the material buckling.
The geometric buckling is an eigenvalue problem that can be solved for different geometries. The table below lists the geometric buckling for some common geometries.
| Geometry | Geometric Buckling Bg2 |
|---|---|
| Sphere of radius R | (π/R)2 |
| Cylinder of height H and radius R | (π/H)2+(2.405/R)2 |
| Parallelepided with lengths 2a, 2b, and 2c | (π/2a)2+(π/2b)2+(π/2c)2 |
Since the diffusion theory calculations overpredict the critical dimensions, an extrapolation distance δ must be subtracted to obtain an estimate of actual values. The buckling could also be calculated using actual dimensions and extrapolated distances using the following table.
Expressions for Geometric Buckling in Terms of Actual Dimensions and Extrapolated Distances[1]
| Geometry | Geometric Buckling Bg2 |
|---|---|
| Sphere of radius R | (π/R+δ)2 |
| Cylinder of height H and radius R | (π/H+2δ)2+(2.405/R+δ)2 |
| Parallelepided with lengths a, b, and c | (π/a+2δ)2+(π/b+2δ)2+(π/c+2δ)2 |
Defining a fast diffusion area or neutron age τ then the thermal non-leakage probability and fast non-leakage probability are respectively
lth1/(1+L2B2)
lf1/(1+τB2)
The effective multiplication factor then becomes
keff=k∞lthlf=1/((1+L2B2)(1+τB2))
In the case of a large reactor, the B4 term can be neglected and we are left with
keff=1/(1+M2B2)
where M2=L2+τ. For a critical reactor keff=1, so solving for B2, the material buckling becomes
BM2=(k∞-1)/M2
By equating the geometric and material buckling, one can determine the critical dimensions of a nuclear reactor.
|
|