critical-sphere.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "from IPython.display import Image\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import openmc"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "                               %%%%%%%%%%%%%%%\n",
      "                          %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                       %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                   %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                                    %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "                ###############      %%%%%%%%%%%%%%%%%%%%%%%%\n",
      "               ##################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "               ###################     %%%%%%%%%%%%%%%%%%%%%%%\n",
      "               ####################     %%%%%%%%%%%%%%%%%%%%%%\n",
      "               #####################     %%%%%%%%%%%%%%%%%%%%%\n",
      "               ######################     %%%%%%%%%%%%%%%%%%%%\n",
      "               #######################     %%%%%%%%%%%%%%%%%%\n",
      "                #######################     %%%%%%%%%%%%%%%%%\n",
      "                ######################     %%%%%%%%%%%%%%%%%\n",
      "                 ####################     %%%%%%%%%%%%%%%%%\n",
      "                   #################     %%%%%%%%%%%%%%%%%\n",
      "                    ###############     %%%%%%%%%%%%%%%%\n",
      "                      ############     %%%%%%%%%%%%%%%\n",
      "                         ########     %%%%%%%%%%%%%%\n",
      "                                     %%%%%%%%%%%\n",
      "\n",
      "                   | The OpenMC Monte Carlo Code\n",
      "         Copyright | 2011-2017 Massachusetts Institute of Technology\n",
      "           License | http://openmc.readthedocs.io/en/latest/license.html\n",
      "           Version | 0.9.0\n",
      "          Git SHA1 | 3b0c044974c59773dac4e3ce261a87e18fc53de5\n",
      "         Date/Time | 2017-08-13 17:56:38\n",
      "    OpenMP Threads | 4\n",
      "\n",
      " Reading settings XML file...\n",
      " Reading geometry XML file...\n",
      " Reading materials XML file...\n",
      " Reading cross sections XML file...\n",
      " Reading U235 from /home/blafleur/Sources/openmc/scripts/nndc_hdf5/U235.h5\n",
      " Reading He3 from /home/blafleur/Sources/openmc/scripts/nndc_hdf5/He3.h5\n",
      " Maximum neutron transport energy: 2.00000E+07 eV for U235\n",
      " Reading tallies XML file...\n",
      " Building neighboring cells lists for each surface...\n",
      " Initializing source particles...\n",
      "\n",
      " ====================>     K EIGENVALUE SIMULATION     <====================\n",
      "\n",
      "  Bat./Gen.      k            Average k         \n",
      "  =========   ========   ====================   \n",
      "        1/1    0.03689                       \n",
      "        2/1    0.03720                       \n",
      "        3/1    0.03818                       \n",
      "        4/1    0.03803                       \n",
      "        5/1    0.03781                       \n",
      "        6/1    0.03776                       \n",
      "        7/1    0.03772                       \n",
      "        8/1    0.03809                       \n",
      "        9/1    0.03855                       \n",
      "       10/1    0.03804                       \n",
      "       11/1    0.03806                       \n",
      "       12/1    0.03774    0.03790 +/- 0.00016\n",
      "       13/1    0.03779    0.03786 +/- 0.00010\n",
      "       14/1    0.03783    0.03786 +/- 0.00007\n",
      "       15/1    0.03778    0.03784 +/- 0.00006\n",
      "       16/1    0.03800    0.03787 +/- 0.00005\n",
      "       17/1    0.03726    0.03778 +/- 0.00010\n",
      "       18/1    0.03769    0.03777 +/- 0.00009\n",
      "       19/1    0.03820    0.03782 +/- 0.00009\n",
      "       20/1    0.03793    0.03783 +/- 0.00008\n",
      "       21/1    0.03818    0.03786 +/- 0.00008\n",
      "       22/1    0.03789    0.03786 +/- 0.00007\n",
      "       23/1    0.03814    0.03788 +/- 0.00007\n",
      "       24/1    0.03810    0.03790 +/- 0.00007\n",
      "       25/1    0.03811    0.03791 +/- 0.00006\n",
      "       26/1    0.03819    0.03793 +/- 0.00006\n",
      "       27/1    0.03765    0.03791 +/- 0.00006\n",
      "       28/1    0.03792    0.03791 +/- 0.00006\n",
      "       29/1    0.03802    0.03792 +/- 0.00005\n",
      "       30/1    0.03789    0.03792 +/- 0.00005\n",
      "       31/1    0.03813    0.03793 +/- 0.00005\n",
      "       32/1    0.03783    0.03792 +/- 0.00005\n",
      "       33/1    0.03774    0.03792 +/- 0.00005\n",
      "       34/1    0.03773    0.03791 +/- 0.00005\n",
      "       35/1    0.03811    0.03792 +/- 0.00004\n",
      "       36/1    0.03750    0.03790 +/- 0.00005\n",
      "       37/1    0.03802    0.03790 +/- 0.00004\n",
      "       38/1    0.03808    0.03791 +/- 0.00004\n",
      "       39/1    0.03807    0.03792 +/- 0.00004\n",
      "       40/1    0.03823    0.03793 +/- 0.00004\n",
      "       41/1    0.03755    0.03791 +/- 0.00004\n",
      "       42/1    0.03826    0.03793 +/- 0.00004\n",
      "       43/1    0.03825    0.03794 +/- 0.00004\n",
      "       44/1    0.03785    0.03793 +/- 0.00004\n",
      "       45/1    0.03805    0.03794 +/- 0.00004\n",
      "       46/1    0.03808    0.03794 +/- 0.00004\n",
      "       47/1    0.03818    0.03795 +/- 0.00004\n",
      "       48/1    0.03733    0.03793 +/- 0.00004\n",
      "       49/1    0.03768    0.03792 +/- 0.00004\n",
      "       50/1    0.03828    0.03793 +/- 0.00004\n",
      "       51/1    0.03744    0.03792 +/- 0.00004\n",
      "       52/1    0.03798    0.03792 +/- 0.00004\n",
      "       53/1    0.03823    0.03793 +/- 0.00004\n",
      "       54/1    0.03826    0.03794 +/- 0.00004\n",
      "       55/1    0.03797    0.03794 +/- 0.00004\n",
      "       56/1    0.03800    0.03794 +/- 0.00004\n",
      "       57/1    0.03806    0.03794 +/- 0.00004\n",
      "       58/1    0.03787    0.03794 +/- 0.00004\n",
      "       59/1    0.03816    0.03794 +/- 0.00004\n",
      "       60/1    0.03794    0.03794 +/- 0.00004\n",
      "       61/1    0.03809    0.03795 +/- 0.00003\n",
      "       62/1    0.03833    0.03796 +/- 0.00003\n",
      "       63/1    0.03775    0.03795 +/- 0.00003\n",
      "       64/1    0.03799    0.03795 +/- 0.00003\n",
      "       65/1    0.03789    0.03795 +/- 0.00003\n",
      "       66/1    0.03812    0.03795 +/- 0.00003\n",
      "       67/1    0.03792    0.03795 +/- 0.00003\n",
      "       68/1    0.03762    0.03795 +/- 0.00003\n",
      "       69/1    0.03762    0.03794 +/- 0.00003\n",
      "       70/1    0.03772    0.03794 +/- 0.00003\n",
      "       71/1    0.03799    0.03794 +/- 0.00003\n",
      "       72/1    0.03812    0.03794 +/- 0.00003\n",
      "       73/1    0.03790    0.03794 +/- 0.00003\n",
      "       74/1    0.03735    0.03793 +/- 0.00003\n",
      "       75/1    0.03789    0.03793 +/- 0.00003\n",
      "       76/1    0.03792    0.03793 +/- 0.00003\n",
      "       77/1    0.03755    0.03793 +/- 0.00003\n",
      "       78/1    0.03776    0.03792 +/- 0.00003\n",
      "       79/1    0.03784    0.03792 +/- 0.00003\n",
      "       80/1    0.03801    0.03792 +/- 0.00003\n",
      "       81/1    0.03828    0.03793 +/- 0.00003\n",
      "       82/1    0.03841    0.03793 +/- 0.00003\n",
      "       83/1    0.03763    0.03793 +/- 0.00003\n",
      "       84/1    0.03797    0.03793 +/- 0.00003\n",
      "       85/1    0.03811    0.03793 +/- 0.00003\n",
      "       86/1    0.03787    0.03793 +/- 0.00003\n",
      "       87/1    0.03784    0.03793 +/- 0.00003\n",
      "       88/1    0.03774    0.03793 +/- 0.00003\n",
      "       89/1    0.03766    0.03793 +/- 0.00003\n",
      "       90/1    0.03811    0.03793 +/- 0.00003\n",
      "       91/1    0.03783    0.03793 +/- 0.00003\n",
      "       92/1    0.03794    0.03793 +/- 0.00003\n",
      "       93/1    0.03776    0.03792 +/- 0.00003\n",
      "       94/1    0.03785    0.03792 +/- 0.00003\n",
      "       95/1    0.03799    0.03792 +/- 0.00003\n",
      "       96/1    0.03804    0.03793 +/- 0.00003\n",
      "       97/1    0.03797    0.03793 +/- 0.00003\n",
      "       98/1    0.03823    0.03793 +/- 0.00003\n",
      "       99/1    0.03780    0.03793 +/- 0.00002\n",
      "      100/1    0.03793    0.03793 +/- 0.00002\n",
      " Creating state point statepoint.100.h5...\n",
      "\n",
      " =======================>     TIMING STATISTICS     <=======================\n",
      "\n",
      " Total time for initialization     =  6.2905E-02 seconds\n",
      "   Reading cross sections          =  3.7924E-02 seconds\n",
      " Total time in simulation          =  9.6865E+00 seconds\n",
      "   Time in transport only          =  9.3313E+00 seconds\n",
      "   Time in inactive batches        =  2.7010E-01 seconds\n",
      "   Time in active batches          =  9.4164E+00 seconds\n",
      "   Time synchronizing fission bank =  4.7027E-02 seconds\n",
      "     Sampling source sites         =  3.0404E-02 seconds\n",
      "     SEND/RECV source sites        =  1.6397E-02 seconds\n",
      "   Time accumulating tallies       =  6.1912E-03 seconds\n",
      " Total time for finalization       =  1.4237E-01 seconds\n",
      " Total time elapsed                =  9.9469E+00 seconds\n",
      " Calculation Rate (inactive)       =  7.40471E+05 neutrons/second\n",
      " Calculation Rate (active)         =  1.91156E+05 neutrons/second\n",
      "\n",
      " ============================>     RESULTS     <============================\n",
      "\n",
      " k-effective (Collision)     =  0.03769 +/-  0.00009\n",
      " k-effective (Track-length)  =  0.03793 +/-  0.00002\n",
      " k-effective (Absorption)    =  0.03807 +/-  0.00025\n",
      " Combined k-effective        =  0.03794 +/-  0.00003\n",
      " Leakage Fraction            =  0.72989 +/-  0.00033\n",
      "\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Instantiate Nuclides\n",
    "u235 = openmc.Nuclide('U235')\n",
    "he3 = openmc.Nuclide('He3')\n",
    "\n",
    "# Create materials\n",
    "fuel = openmc.Material(name='fuel')\n",
    "fuel.set_density('g/cm3', 10.31341)\n",
    "fuel.add_nuclide(u235, 1.0)\n",
    "\n",
    "air = openmc.Material(name='fuel')\n",
    "air.set_density('g/cm3', 10.31341)\n",
    "air.add_nuclide(he3, 1.0)\n",
    "\n",
    "materials_file = openmc.Materials((fuel, air))\n",
    "materials_file.default_xs = '71c'\n",
    "materials_file.export_to_xml()\n",
    "\n",
    "radius = 0.5\n",
    "settings_file = openmc.Settings()\n",
    "settings_file.batches = 100\n",
    "settings_file.inactive = 10\n",
    "settings_file.particles = 20000\n",
    "# Create geometry\n",
    "sphere = openmc.Sphere(R=radius)\n",
    "min_x = openmc.XPlane(x0=-radius, boundary_type='vacuum') \n",
    "max_x = openmc.XPlane(x0=+radius, boundary_type='vacuum') \n",
    "min_y = openmc.YPlane(y0=-radius, boundary_type='vacuum') \n",
    "max_y = openmc.YPlane(y0=+radius, boundary_type='vacuum') \n",
    "min_z = openmc.ZPlane(z0=-radius, boundary_type='vacuum')\n",
    "max_z = openmc.ZPlane(z0=+radius, boundary_type='vacuum') \n",
    "\n",
    "# Create Universe\n",
    "universe = openmc.Universe(name='Universe')\n",
    "\n",
    "fuel_cell = openmc.Cell(name='fuel')\n",
    "fuel_cell.fill = fuel\n",
    "fuel_cell.region = -sphere\n",
    "universe.add_cell(fuel_cell)\n",
    "\n",
    "air_cell = openmc.Cell(name='air')\n",
    "air_cell.fill = air\n",
    "air_cell.region = +sphere\n",
    "universe.add_cell(air_cell)\n",
    "\n",
    "# Create root cell\n",
    "root_cell = openmc.Cell(name='root_cell')\n",
    "root_cell.fill = universe\n",
    "root_cell.region = +min_x & -max_x & +min_y & -max_y & +min_z & -max_z\n",
    "root_universe = openmc.Universe(universe_id=0, name='root universe')\n",
    "root_universe.add_cell(root_cell)\n",
    "\n",
    "geometry = openmc.Geometry()\n",
    "geometry.root_universe = root_universe\n",
    "geometry.export_to_xml()\n",
    "\n",
    "# Create source\n",
    "bounds = [-radius, -radius, -radius, radius, radius, radius]\n",
    "uniform_dist = openmc.stats.Box(bounds[:3], bounds[3:], only_fissionable=True)\n",
    "settings_file.source = openmc.source.Source(space=uniform_dist)\n",
    "settings_file.export_to_xml()\n",
    "\n",
    "# Run the calculation\n",
    "openmc.run()\n",
    "\n",
    "# Load the statepoint file\n",
    "# sp = openmc.StatePoint('statepoint.100.h5')\n",
    "# print(openmc.StatePoint.k_combined)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}