Calculation code for spherical nuclear reactor

The governing equation being used is the steady state neutron diffusion equation:

eq1

eq2

Numerical schemes being used are:

  • Central finite difference for flux calculation
  • Gauss-Siedel and S.O.R for flux calculation
  • Power method for criticality calculation

Code package:

Flowchart of the code:

powermethod

Some previews:

snr01snr02snr03

Leave a Reply