Neutron density, flux and reaction rates

In reactor physics it is customary [48] to define the number of neutrons per unit volume, per velocity bin and solid angle unit n(r, v, fi, t) such that the number of neutrons in volume d3r, at a position r, with a velocity between v and v + dv pointing in direction fi (unit vector) within solid angle d2^ is n(r, v, fi, t) d3r dv d2^. The flux of neutrons is defined as

ф(г, v, fi, t)=vn(r, v, fi, t). (39)

The number of neutrons per time unit with velocity v and direction fi which cross a planar unit surface at position r and unit normal vector u is:

image030

ф(г, v, fi, t)fi • u. The case where ф is isotropic is particularly interesting, since it is nearly satisfied in reactors. Then, if measuring angles with respect to the normal vector u, we obtain the total number of neutrons crossing the surface, regardless of their direction, by

Подпись: 4:ф(г, v, fi, t)■»::/2

cos(9) sin(0) d9 = 2:ф(г, v, fi, t).

0

image032 image033

We now consider a thin slab with unit surface and an atomic thickness of ns identical nuclei per unit surface. The nuclei have reaction cross-section a. The number of reactions in the slab per unit time reads

The total flux is the directional flux integrated over angle, thus

‘(r, v, t) = 4:ф(г, v, fi, t).

This quantity is usually called the neutron flux. In term of this quantity, the number of reactions per time unit is thus

image034

Figure 3.4. Energy dependence of q for the principal fissile nuclei.

 

while the number of neutrons crossing a unit surface plane per time unit is

1 Mr, V; t)).[12]

Instead of computing the number of reactions per unit time, it is instructive to compute the total length travelled by all neutrons traversing a thin unit surface slab, which we assume to have thickness l, that is very small compared with the transverse dimensions of the slab. Then the total length travelled by the neutrons is

Подпись: L = 4гкф(г, v, fi, t)l^n/2

sin(h) dh = ‘(r, v, t)l

0

Подпись: '(r, v, t) Подпись: L V Подпись: (3.11)

where the volume V of the slab is simply l, since it has unit surface. Thus we obtain an expression for the flux:

The formula was demonstrated for an infinitely thin slab. However, for iso­tropic neutron fluxes, it can be generalized to any arbitrary volume. Indeed, any volume can be subdivided, at will, into n small, thin, parallelograms. For each of the elementary volumes we have the flux value ‘ = Lj/V. The average flux over the volume is the average of the ‘ weighted by the volume,

Подпись: (3.12)и = V" = Li = L

(‘) E V’ v v

since the total length travelled by the neutrons is clearly the sum of the elementary lengths. This is an important formula since it is the one used in all Monte Carlo simulations. It can be shown that this formula also holds for anisotropic neutron fluxes.

Подпись: J(r, V; t) Подпись: Нф(г, V, H, t) d2H. (4w) Подпись: (3.13)

Note that the definition of the neutron flux is different from the usual definition of flux in other fields of physics. For example, in thermodynamics, a finite heat flux through a surface requires a temperature gradient across this surface. The analogy of the heat flux in neutronics is the neutron current defined as

image043 image044 image045

This current is only different from 0 for anisotropic neutron fluxes. One-sided currents are also frequently used and defined by

where N is the direction in which the flux is measured.

The number of interactions of type (a) per cm3 is S(a)p where ‘ is the neutron flux expressed in n/cm2/s. The most important macroscopic cross­sections are the scattering cross-section XS, the absorption cross-section Xa and the fission cross-section Xf = XaPf where Pf is the fission probability.