## Fick’s Law

The derivation of the diffusion equation depends on **Fick’s law**, which states that solute diffuses (**neutron current**) from high concentration (high flux) to low concentration. As can be seen we have to investigate the relationship between the **flux (φ)** and the **current (J)**. This relationship between the flux (φ) and the current (J) is identical in form to a law (**the Fick’s law**) used in the study of physical diffusion in liquids and gases.

**In chemistry, Fick’s law states that**:

*If the concentration of a solute in one region is greater than in another of a solution, the solute diffuses from the region of higher concentration to the region of lower concentration, with a magnitude that is proportional to the concentration gradient.*

In one (spatial) dimension, the law is:

where:

*J*is the diffusion flux,*D*is the**diffusion coefficient,***φ*(for ideal mixtures) is the concentration.

The use of this law in **nuclear reactor theory** leads to the **diffusion approximation**.

**The Fick’s law in reactor theory stated that**:

*The current density vector J is proportional to the negative of the gradient of the neutron flux. The proportionality constant is called the diffusion coefficient and is denoted by the symbol D.*

In one (spatial) dimension, the law is:

where:

(neutrons.cm*J*is the neutron current density^{-2}.s^{-1}) along x-direction, the net flow of neutrons that pass per unit of time through a unit area perpendicular to the x-direction.*D*is the**diffusion coefficient,**it has unit of cm and it is given by:*φ*is the neutron flux, which is the number of neutrons crossing through some arbitrary cross-sectional unit area in**all directions**per unit time.

The generalized Fick’s law (in three dimension) is:

where **J** denotes the **diffusion flux vector**. Note that the gradient operator turns the neutron flux, which is a **scalar quantity** into the neutron current, which is a **vector quantity**.

### Physical Interpretation

The physical interpretation is similar to fluxes of gases. The neutrons exhibit a net flow in the direction of least density. This is a natural consequence of **greater collision densities** at positions of **greater neutron densities**.

Consider neutrons passing through the plane at x=0 from left to right as the result of collisions to the left of the plane. Since the concentration of neutrons and the flux is larger for negative values of x, there are **more collisions per cubic centimeter on the left**. Therefore more neutrons are scattered from left to right, then the other way around. Thus the neutrons naturally diffuse toward the right.

**neutron flux**is often used, because it expresses the

**total path length**covered by

**all neutrons**. The total distance these neutrons can travel each second is determined by their velocity and therefore the

**neutron flux density**value is calculated as the

**neutron density (n)**multiplied by

**neutron velocity (v)**.

**Ф = n.v**

where:

**Ф – neutron flux (neutrons.cm ^{-2}.s^{-1})**

**n – neutron density (neutrons.cm**

^{-3})**v – neutron velocity (cm.s**

^{-1})**The neutron flux**, which is the number of neutrons crossing through some arbitrary cross-sectional unit area in **all directions** per unit time, is a **scalar quantity**. Therefore it is also known as the **scalar flux**. The expression **Ф(E).dE** is the total distance traveled during one second by all neutrons with energies between E and dE located in 1 cm^{3}.

The connection to the **reaction rate**, respectively the reactor power, is obvious. Knowledge of the **neutron flux** (the **total path length** of all the neutrons in a cubic centimeter in a second) and the** macroscopic cross sections** (the probability of having an interaction **per centimeter path length**) allows us to compute the rate of interactions (e.g. rate of fission reactions). **The reaction rate** (the number of interactions taking place in that cubic centimeter in one second) is then given by multiplying them together:

where:

**Ф – neutron flux (neutrons.cm ^{-2}.s^{-1})**

**σ – microscopic cross section (cm**

^{2})**N – atomic number density (atoms.cm**

^{-3})This flux of neutron flux is called the neutron current density. We have to distinguish between the** neutron flux** and the **neutron current density**. Although both physical quantities have the** same units**, namely, neutrons.cm^{-2}.s^{-1}, their physical interpretations are different. In contrast to the neutron flux, the neutron intensity is the number of neutrons crossing through some arbitrary cross-sectional unit area** in a single direction** per unit time (a surface is perpendicular to the direction of the beam). The **neutron current density** is a **vector quantity**.

The vector **J is defined as the following integral:**

A typical **thermal reactor** contains about **100 tons** of uranium with an average enrichment of **2%** (do not confuse it with the enrichment of the **fresh fuel**). If the reactor power is 3000MW_{th}, determine the **reaction rate** and the **average core thermal flux**.

**Solution:**

Multiplying the reaction rate per unit volume (RR = Ф . Σ) by the total volume of the core (V) gives us the** total number of reactions** occurring in the reactor core per unit time. But we also know the amount of energy released per one fission reaction to be about **200 MeV/fission**. Now, it is possible to determine the **rate of energy release** (power) due to the fission reaction. It is given by following equation:

**P = Ф . Σ _{f} . E_{r} . V = Ф . N_{U235} . σ_{f}^{235} . E_{r} . V**

where:

**P – reactor power (MeV.s ^{-1})**

**Ф – neutron flux (neutrons.cm**

^{-2}.s^{-1})**σ – microscopic cross section (cm**

^{2})**N – atomic number density (atoms.cm**

^{-3})**Er – the average recoverable energy per fission (MeV / fission)**

**V – total volume of the core (m**

^{3})The amount of fissile ^{235}U per the volume of the reactor core.

m_{235} [g/core] = 100 [metric tons] x 0.02 [g of ^{235}U / g of U] . 10^{6} [g/metric ton]
= **2 x 10 ^{6} grams of ^{235}U** per the volume of the reactor core

The atomic number density of ^{235}U in the volume of the reactor core:

N_{235} . V = m_{235} . N_{A} / M_{235}

= 2 x 10^{6} [g 235 / core] x 6.022 x 10^{23} [atoms/mol] / 235 [g/mol]
= **5.13 x 10 ^{27} atoms / core**

The microscopic fission cross-section of

^{235}U (for thermal neutrons):

**σ _{f}^{235} = 585 barns**

The average recoverable energy per ^{235}U fission:

**E _{r} = 200.7 MeV/fission**

## Validity of Fick’s Law

**It must be emphasized that Fick’s law is an approximation and was derived under the following conditions:**

**Infinite medium.**This assumption is necessary to allow integration over all space but flux contributions are negligible beyond a few mean free paths (about three mean free paths) from boundaries of the diffusive medium.**Sources or sinks.**Derivation of Fick’s law assumes that the contribution to the flux is mostly from elastic scattering reactions. Source neutrons contribute to the flux if they are more than a few mean free paths from a source.**Uniform medium.**Derivation of Fick’s law assumes that a uniform medium was used. There are different scattering properties at the boundary (interface) between two media.**Isotropic scattering**. Isotropic scattering occurs at low energies, but is not true in general. Presence of anisotropic scattering can be corrected by modification of the diffusion coefficient (based on transport theory).**Low absorbing medium**. Derivation of Fick’s law assumes (an expansion in a Taylor’s series) that the neutron flux,*φ,*is slowly varying. Large variations in φ occur when Σ_{a}(neutron absorption) is large (compared to Σ_{s}).**Σ**_{a}**<< Σ**_{s}**Time – independent flux.**Derivation of Fick’s law assumes that the neutron flux is independent of time.

To some extent, these limitations are valid in every practical reactor. Nevertheless Fick’s law gives a reasonable approximation. For more detailed calculations, higher order methods are available.

**Nuclear and Reactor Physics:**

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

**Advanced Reactor Physics:**

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