Syeilendra Pramuditya

Calculation code for spherical nuclear reactor

Posted by Syeilendra Pramuditya on June 13, 2009

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:

 

Leave a Reply

XHTML: You can use these tags: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <pre> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>