Main_Sim.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "20b39670-179f-4da3-8dab-cd2935858f33",
   "metadata": {},
   "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",
      "                 | The OpenMC Monte Carlo Code\n",
      "       Copyright | 2011-2023 MIT, UChicago Argonne LLC, and contributors\n",
      "         License | https://docs.openmc.org/en/latest/license.html\n",
      "         Version | 0.13.4-dev\n",
      "        Git SHA1 | 6c7f857e1e98ce4fd5b39d8be034308596628d7e\n",
      "       Date/Time | 2023-08-10 15:11:42\n",
      "  OpenMP Threads | 8\n",
      "\n",
      " Reading settings XML file...\n",
      " Reading cross sections XML file...\n",
      " Reading materials XML file...\n",
      " Reading geometry XML file...\n",
      " Reading Cu63 from /nuclear_data/ENDFB-8.0-NNDC_Cu63.h5\n",
      " Reading Cu65 from /nuclear_data/ENDFB-8.0-NNDC_Cu65.h5\n",
      " Reading Sn112 from /nuclear_data/ENDFB-8.0-NNDC_Sn112.h5\n",
      " Reading Sn114 from /nuclear_data/ENDFB-8.0-NNDC_Sn114.h5\n",
      " Reading Sn115 from /nuclear_data/ENDFB-8.0-NNDC_Sn115.h5\n",
      " Reading Sn116 from /nuclear_data/ENDFB-8.0-NNDC_Sn116.h5\n",
      " Reading Sn117 from /nuclear_data/ENDFB-8.0-NNDC_Sn117.h5\n",
      " Reading Sn118 from /nuclear_data/ENDFB-8.0-NNDC_Sn118.h5\n",
      " Reading Sn119 from /nuclear_data/ENDFB-8.0-NNDC_Sn119.h5\n",
      " Reading Sn120 from /nuclear_data/ENDFB-8.0-NNDC_Sn120.h5\n",
      " Reading Sn122 from /nuclear_data/ENDFB-8.0-NNDC_Sn122.h5\n",
      " Reading Sn124 from /nuclear_data/ENDFB-8.0-NNDC_Sn124.h5\n",
      " Reading N14 from /nuclear_data/ENDFB-8.0-NNDC_N14.h5\n",
      " Reading N15 from /nuclear_data/ENDFB-8.0-NNDC_N15.h5\n",
      " Reading O16 from /nuclear_data/ENDFB-8.0-NNDC_O16.h5\n",
      " Reading O17 from /nuclear_data/ENDFB-8.0-NNDC_O17.h5\n",
      " Reading O18 from /nuclear_data/ENDFB-8.0-NNDC_O18.h5\n",
      " Reading Ar36 from /nuclear_data/ENDFB-8.0-NNDC_Ar36.h5\n",
      " WARNING: Negative value(s) found on probability table for nuclide Ar36 at 294K\n",
      " Reading Ar38 from /nuclear_data/ENDFB-8.0-NNDC_Ar38.h5\n",
      " Reading Ar40 from /nuclear_data/ENDFB-8.0-NNDC_Ar40.h5\n",
      " Reading C12 from /nuclear_data/ENDFB-8.0-NNDC_C12.h5\n",
      " Reading C13 from /nuclear_data/ENDFB-8.0-NNDC_C13.h5\n",
      " Minimum neutron data temperature: 294 K\n",
      " Maximum neutron data temperature: 294 K\n",
      " Reading tallies XML file...\n",
      " Preparing distributed cell instances...\n",
      " Reading plot XML file...\n",
      " Writing summary.h5 file...\n",
      " Maximum neutron transport energy: 20000000 eV for Sn112\n",
      "\n",
      " ===============>     FIXED SOURCE TRANSPORT SIMULATION     <===============\n",
      "\n",
      " Simulating batch 1\n",
      " Simulating batch 2\n",
      " Simulating batch 3\n",
      " Simulating batch 4\n",
      " Simulating batch 5\n",
      " Simulating batch 6\n",
      " Simulating batch 7\n",
      " Simulating batch 8\n",
      " Simulating batch 9\n",
      " Simulating batch 10\n",
      " Creating state point statepoint.10.h5...\n",
      "\n",
      " =======================>     TIMING STATISTICS     <=======================\n",
      "\n",
      " Total time for initialization     = 6.2979e-01 seconds\n",
      "   Reading cross sections          = 6.0642e-01 seconds\n",
      " Total time in simulation          = 1.0285e-01 seconds\n",
      "   Time in transport only          = 9.5595e-02 seconds\n",
      "   Time in active batches          = 1.0285e-01 seconds\n",
      "   Time accumulating tallies       = 4.2690e-03 seconds\n",
      "   Time writing statepoints        = 2.4154e-03 seconds\n",
      " Total time for finalization       = 8.7210e-05 seconds\n",
      " Total time elapsed                = 7.3753e-01 seconds\n",
      " Calculation Rate (active)         = 972306 particles/second\n",
      "\n",
      " ============================>     RESULTS     <============================\n",
      "\n",
      " Leakage Fraction            = 0.89522 +/- 0.00066\n",
      "\n",
      " ==============>     TALLY 1: FLUX OF TOP TEST CYLINDER     <===============\n",
      "\n",
      " Cell 3\n",
      "   Total Material\n",
      "     Flux                                 0.195518 +/- 0.000955274\n",
      " =============>     TALLY 2: FLUX OF BOTTOM TEST CYLINDER     <=============\n",
      "\n",
      " Cell 2\n",
      "   Total Material\n",
      "     Flux                                 0.430725 +/- 0.0013198\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQsAAAEHCAYAAAC5l2MpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAWxElEQVR4nO3dfVAU9/0H8Pdh4DjhuAzKg1QBjTGKD8ShmEgcPAyDUI0hMDZTTYP1sZ0DRdpMg4kjHaelNRIzYVBMoxAnQRySoCSO1IhwxBSaEWMqqfGpmiiP0cgdYIPU2/7hj/vl5Ol7cMcu8H7N3Iy3993dDzG83f3e7n5UkiRJICLqh4vcBRDR8MCwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhD8ldgLNZLBbU19dDq9VCpVLJXQ6RokiShNbWVgQEBMDFpe9jhxEfFvX19Zg0aZLcZRAp2vXr1zFx4sQ+x4z4sNBqtQCAtz/8DGM9PGWuhkhZ7rS3YW3CU9bfk76M+LDoOvUY6+GJsR79/wchGo1ETtE5wUlEQhgWRCSEYUFEQhgWRCSEYUFEQhgWRCSEYUFEQmQNi8zMTISHh0Or1cLX1xfx8fG4cOGC9fPvv/8eKSkpeOyxx6DRaBAYGIiNGzfCZDLJWDXR6CRrWBiNRhgMBlRXV+OTTz5BZ2cnYmJi0N7eDuD+pdr19fXYuXMnamtrkZ+fj9LSUqxZs0bOsolGJZWSnu793XffwdfXF0ajEZGRkT2OKSoqwgsvvID29nY89FD/F6CazWbodDoU/O1LXsFJ9IA77a1YsTgUJpMJXl5efY5V1OXeXacX3t7efY7x8vLqNSg6OjrQ0dFhfW82mx1bJNEopZiwsFgsSE1NxVNPPYVZs2b1OObmzZvYvn071q9f3+t2MjMz8Yc//MFZZQ5LzwaPkbsExZIkCSXfWOQuY1hQzLchBoMBtbW1KCws7PFzs9mMJUuWICQkBBkZGb1uJz09HSaTyfq6fv26kyomGl0UcWSRnJyMjz/+GJWVlT3eU9/a2orY2FhotVoUFxfD1dW1122p1Wqo1Wpnlks0Ksl6ZCFJEpKTk1FcXIyTJ09i8uTJ3caYzWbExMTAzc0NJSUlcHd3l6FSIpL1yMJgMKCgoABHjhyBVqtFY2MjAECn00Gj0ViD4s6dO3j33XdhNputE5Y+Pj4YM4bn4kRDRdaw2LNnDwBAr9fbLM/Ly8OqVatw5swZ/OMf/wAATJ061WbM1atXERwcPBRlEhFkDov+LvHQ6/X9jiGioaGYb0OISNkYFkQkhGFBREIYFkQkhGFBREIYFkQkhGFBREIYFkQkhGFBREIYFkQkhGFBREIYFkQkhGFBREIYFkQkhGFBREIYFkQkhGFBREIU3esUAN566y3o9Xp4eXlBpVKhpaVFnmKJRjlF9zoFgDt37iA2NhZbtmyRsVIikvUZnKWlpTbv8/Pz4evri5qaGmuv09TUVABARUXFEFdHRD+mqDkLkV6nRCQPRXQkA8R6nYpgY+Tu+IR0cgTFhEVXr9NTp04NajtsjPwAScLlSr3cVSiWBABBJ+UuY1hQxGlIV6/T8vLyHnud2oONkYmcQ/YmQykpKSguLkZFRUWPvU7txcbIRM6h6F6nANDY2IjGxkZcvnwZAHDu3DlotVoEBgZyIpRoCMl6GrJnzx6YTCbo9XpMmDDB+jp06JB1TG5uLubOnYt169YBACIjIzF37lyUlJTIVTbRqCT7aUh/MjIykJGR4fxiiKhPipjgJCLlY1gQkRCGBREJYVgQkRCGBREJYVgQkRCGBREJYVgQkRCGBREJYVgQkRCGBREJYVgQkRCGBREJYVgQkRCGBREJYVgQkRCGBREJYVgQkRDFN0b+4YcfYDAYMG7cOHh6eiIxMRFNTU0yVUw0eim+MfLmzZvx0UcfoaioCEajEfX19UhISJCxaqLRSdGNkU0mE/bt24eCggIsWrQIAJCXl4cZM2aguroaTz75pBxlE41KipqzeLAxck1NDTo7OxEdHW0dM336dAQGBqKqqqrHbXR0dMBsNtu8iGjwFNPrtKfGyI2NjXBzc8PDDz9sM9bPz8/akOhB7HX6AJUKjy40yl2FYkmShPPfWOQuY1hQzJFFV2PkwsLCQW2HvU6JnEMRRxZdjZErKyttGiP7+/vj7t27aGlpsTm6aGpqgr+/f4/bYq9TIueQ9chCkiQkJyejuLgYJ0+e7NYYOSwsDK6urigrK7Muu3DhAr799lvMnz9/qMslGtUU3RhZp9NhzZo1SEtLg7e3N7y8vJCSkoL58+fzmxCiISZrWOzZswcAoNfrbZbn5eVh1apVAIBdu3bBxcUFiYmJ6OjowOLFi7F79+4hrpSIFN8Y2d3dHTk5OcjJyRmCioioN4r5NoSIlI1hQURChE5DBnIvRm5uLnx9fe1ej4iUSSgsDh8+jJ///OfQaDRCGy0oKEBbWxvDgmgEEZ7gfPPNN4V/+d9///0BF0REyiQ0Z1FeXm69uUvEsWPH8JOf/GTARRGR8ggdWSxcuNCujS5YsGBAxRCRcg34Oovm5mY0NzfDYrG9Y2/OnDmDLoqIlMfusKipqUFSUhLOnz9vvahKpVJBkiSoVCrcu3fP4UUSkfzsDovVq1dj2rRp2LdvH/z8/KBSqZxRFxEpjN1h8e9//xsffPABpk6d6ox6iEih7L6C8+mnn8aXX37pjFqISMHsPrJ4++23kZSUhNraWsyaNQuurq42ny9btsxhxRGRctgdFlVVVfjss89w7Nixbp9xgpNo5LL7NCQlJQUvvPACGhoaYLFYbF4MCqKRy+6wuHXrFjZv3gw/Pz9n1ENECmV3WCQkJKC8vNwZtRCRgtk9ZzFt2jSkp6fj1KlTmD17drcJzo0bNzqsOCJSjgF9G+Lp6Qmj0Qij0bZ5jUqlsissKisr8dprr6GmpgYNDQ0oLi5GfHy89fOmpib8/ve/x/Hjx9HS0oLIyEhkZ2fj0UcftbdsIhoku8Pi6tWrDtt5e3s7QkNDsXr16m4P2JEkCfHx8XB1dcWRI0fg5eWF119/HdHR0fjXv/4FDw8Ph9VBRP2T9YG9cXFxiIuL6/GzS5cuobq6GrW1tZg5cyaA+08D9/f3x8GDB7F27dqhLJVo1LN7gjMxMRF/+ctfui3fsWMHli9f7pCigPsNjoH7T/fu4uLiArVajVOnTvW5Hhsj25Ikia9eXiTO7iOLyspKZGRkdFseFxeHrKwsR9QE4P+7paenp2Pv3r3w8PDArl27cOPGDTQ0NPS6Hhsjd1fCxr/kAHYfWbS1tcHNza3bcldXV4f+K+7q6ooPP/wQFy9ehLe3N8aOHYvy8nLExcXBxaX3stkYmcg57A6L2bNn49ChQ92WFxYWIiQkxCFFdQkLC8PZs2fR0tKChoYGlJaW4tatW5gyZUqv66jVanh5edm8iGjw7D4N2bp1KxISEnDlyhUsWrQIAFBWVoaDBw+iqKjI4QUC93ufAvcnPU+fPo3t27c7ZT9E1Du7w+KZZ57B4cOH8ac//Qnvv/8+NBoN5syZgxMnTtj9rM62tjZcvnzZ+v7q1as4e/YsvL29ERgYiKKiIvj4+CAwMBDnzp3Dpk2bEB8fj5iYGHvLJqJBGtBXp0uWLMGSJUsGvfPTp08jKirK+j4tLQ0AkJSUhPz8fDQ0NCAtLQ1NTU2YMGECXnzxRWzdunXQ+yUi+6mkEf79kdlshk6nQ8HfvsRYD63c5RApyp32VqxYHAqTydTv/J7QBKe3tzdu3rwpXEBgYCC++eYb4fFEpHxCpyEtLS04duyYdaKxP7du3eKzLYhGGOE5i6SkJGfWQUQKJxQWDzYSIqLRx+6LsohodGJYEJEQhgURCWFYEJEQhgURCbE7LBYtWtTj8yJu375tvbGMiEYeu+8NqaiowLlz5/DFF1/gvffesz4L8+7du90e4EtEI8eATkNOnDiBxsZGPPnkk7h27ZqDSyIiJRpQWEyYMAFGoxGzZ89GeHg4KioqHFwWESmN3WGhUqkA3H8iVUFBATZt2oTY2Fjs3r3b4cURkXLYPWfx4B3tr776KmbMmMF7R4hGuAE1GfLx8bFZlpiYiOnTp+P06dMOK4yIlMXusAgKCupx+cyZM63NgIho5JH1oqzKyko888wzCAgIgEqlwuHDh20+b2trQ3JyMiZOnAiNRoOQkBDk5ubKUyzRKCdrWHT1Os3Jyenx87S0NJSWluLdd9/F+fPnkZqaiuTkZJSUlAxxpUSk2F6nAPD3v/8dSUlJ0Ov1AID169dj7969+Pzzz7Fs2bIhqpKIAIXfGxIREYGSkhLU1dVBkiSUl5fj4sWLfbYCYK9TIudQdFhkZ2cjJCQEEydOhJubG2JjY5GTk4PIyMhe18nMzIROp7O+Jk2aNIQVE41cig+L6upqlJSUoKamBllZWTAYDDhx4kSv67DXKZFzyDpn0Zf//Oc/2LJlC4qLi60NjebMmYOzZ89i586diI6O7nE9tVoNtVo9lKUSjQqKPbLo7OxEZ2dnt47pY8aM4QOEiWQg65FFf71OFy5ciJdeegkajQZBQUEwGo04cOAAXn/9dRmrJhqdZA2L/nqdFhYWIj09HStXrsT333+PoKAg/PGPf8Svf/1ruUomGrVkDQu9Xt/txrQf8/f3R15e3hBWRES9UeycBREpC8OCiIQwLIhICMOCiIQwLIhICMOCiIQwLIhICMOCiIQwLIhICMOCiIQwLIhICMOCiIQwLIhICMOCiIQwLIhICMOCiIQwLIhICMOCiIQoujGySqXq8fXaa6/JUzDRKKboxsgNDQ02r/3790OlUiExMXGIKyUiRTdG9vf3t3l/5MgRREVFYcqUKc4ujYgeoNiOZA9qamrC0aNH8c477/Q5rqOjAx0dHdb3bIxM5BjDZoLznXfegVarRUJCQp/j2BiZyDmGTVjs378fK1euhLu7e5/j2BiZyDmGxWnIp59+igsXLuDQoUP9jmVjZCLnGBZHFvv27UNYWBhCQ0PlLoVo1FJ0Y2Tg/gRlUVERsrKy5CqTiKDwxsgAUFhYCEmS8Itf/EKOEono/6ikvjoTjwBmsxk6nQ4Ff/sSYz20cpdDpCh32luxYnEoTCYTvLy8+hw7LOYsiEh+DAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhiu51CgDnz5/HsmXLoNPp4OHhgfDwcHz77bdDXyzRKKfoXqdXrlzBggULMH36dFRUVOCf//wntm7d2m/vECJyPEX3On3llVfws5/9DDt27LAue+SRR4aiNCJ6gGLnLCwWC44ePYpp06Zh8eLF8PX1xRNPPNHjqcqPdXR0wGw227yIaPAUGxbNzc1oa2vDn//8Z8TGxuL48eN47rnnkJCQAKPR2Ot67HVK5ByKDQuLxQIAePbZZ7F582Y8/vjjePnll7F06VLk5ub2uh57nRI5h2J7nY4fPx4PPfQQQkJCbJbPmDEDp06d6nU99jolcg7FHlm4ubkhPDwcFy5csFl+8eJFBAUFyVQV0eil6F6nL730Ep5//nlERkYiKioKpaWl+Oijj1BRUSFf0USjlKJ7nT733HPIzc1FZmYmNm7ciMceewwffPABFixYIFfJRKMWe50SjWLsdUpEDsewICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEqLYe0McpesykjvtbTJXQqQ8Xb8XIpdbjfiLsm7cuMHb1In6cf36dUycOLHPMSM+LCwWC+rr66HVaqFSqeQup09msxmTJk3C9evX+72abrjiz6gskiShtbUVAQEBcHHpe1ZixJ+GuLi49JuYSuPl5aX4/8kGiz+jcuh0OqFxnOAkIiEMCyISwrBQELVajW3bto3oJ33xZxy+RvwEJxE5Bo8siEgIw4KIhDAsiEgIw4KIhDAsFOjatWtYs2YNJk+eDI1Gg0ceeQTbtm3D3bt35S5t0HJychAcHAx3d3c88cQT+Pzzz+UuyWEyMzMRHh4OrVYLX19fxMfHd2tlMZwxLBTo66+/hsViwd69e/HVV19h165dyM3NxZYtW+QubVAOHTqEtLQ0bNu2DWfOnEFoaCgWL16M5uZmuUtzCKPRCIPBgOrqanzyySfo7OxETEwM2tvb5S7NMSQaFnbs2CFNnjxZ7jIGZd68eZLBYLC+v3fvnhQQECBlZmbKWJXzNDc3SwAko9EodykOwSOLYcJkMsHb21vuMgbs7t27qKmpQXR0tHWZi4sLoqOjUVVVJWNlzmMymQBgWP+9/RjDYhi4fPkysrOzsWHDBrlLGbCbN2/i3r178PPzs1nu5+eHxsZGmapyHovFgtTUVDz11FOYNWuW3OU4BMNiCL388stQqVR9vr7++mubderq6hAbG4vly5dj3bp1MlVO9jIYDKitrUVhYaHcpTjMiL9FXUl++9vfYtWqVX2OmTJlivXP9fX1iIqKQkREBN566y0nV+dc48ePx5gxY9DU1GSzvKmpCf7+/jJV5RzJycn4+OOPUVlZOewej9AXhsUQ8vHxgY+Pj9DYuro6REVFISwsDHl5ef0+mETp3NzcEBYWhrKyMsTHxwO4f6heVlaG5ORkeYtzEEmSkJKSguLiYlRUVGDy5Mlyl+RQDAsFqqurg16vR1BQEHbu3InvvvvO+tlw/lc4LS0NSUlJ+OlPf4p58+bhjTfeQHt7O371q1/JXZpDGAwGFBQU4MiRI9Bqtda5GJ1OB41GI3N1DiD31zHUXV5engSgx9dwl52dLQUGBkpubm7SvHnzpOrqarlLcpje/s7y8vLkLs0heIs6EQkZ3ifCRDRkGBZEJIRhQURCGBZEJIRhQURCGBZEJIRhQURCGBbkUNeuXbPeFPf44487dV/5+fnWfaWmpjp1X8SwICc5ceIEysrKnLqP559/Hg0NDZg/f75T90P38d4Qcopx48Zh3LhxTt2HRqOBRqOBm5ubU/dD9/HIgnp14MABjBs3Dh0dHTbL4+Pj8ctf/tLu7e3fvx8zZ86EWq3GhAkTbO42ValU2Lt3L5YuXYqxY8dixowZqKqqwuXLl6HX6+Hh4YGIiAhcuXJl0D8XDQzDgnq1fPly3Lt3DyUlJdZlzc3NOHr0KFavXm3Xtvbs2QODwYD169fj3LlzKCkpwdSpU23GbN++HS+++CLOnj2L6dOnY8WKFdiwYQPS09Nx+vRpSJI0Ym5nH5ZkvpGNFO43v/mNFBcXZ32flZUlTZkyRbJYLD2Ov3r1qgRA+uKLL2yWBwQESK+88kqv+wEgvfrqq9b3VVVVEgBp37591mUHDx6U3N3du627cOFCadOmTYI/EQ0UjyyoT+vWrcPx48dRV1cH4P43EKtWrYJKpRLeRnNzM+rr6/H000/3OW7OnDnWP3c9q3P27Nk2y3744QeYzWZ7fgRyEE5wUp/mzp2L0NBQHDhwADExMfjqq69w9OhRu7Yh+uAXV1dX65+7wqinZRaLxa79k2PwyIL6tXbtWuTn5yMvLw/R0dGYNGmSXetrtVoEBwc7/atUci6GBfVrxYoVuHHjBv7617/aPbHZJSMjA1lZWXjzzTdx6dIlnDlzBtnZ2Q6ulJyJYUH90ul0SExMhKenp/Vhu/ZKSkrCG2+8gd27d2PmzJlYunQpLl265NhCyak4Z0FC6urqsHLlSqjV6gFvY8OGDb02SpIeeLpjcHBwt2V6vb7bMho6PLKgPt2+fdv6aHuDwSC8XkREBCIiIpxYGfDee+/B09MTn376qVP3Q/fxgb3Up+DgYNy+fRtbt27F7373u37H//e//8W1a9cAAGq12u7JUHu0trZamxY9/PDDGD9+vNP2RQwLIhLE0xAiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEsKwICIhDAsiEvI/NldlLKMguFUAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 258.065x259.74 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "\"\"\"\n",
    "CATHERINE DICKSTEIN\n",
    "JACK SZYPROWSKI\n",
    "8-10-23\n",
    "\"\"\"\n",
    "\n",
    "%matplotlib inline\n",
    "import openmc\n",
    "import os\n",
    "import numpy as np\n",
    "\n",
    "#MATERIALS\n",
    "#Creating bronze material for cylinder\n",
    "castBronze = openmc.Material(1, \"castBronze\")\n",
    "castBronze.add_element(\"Cu\", 0.90)\n",
    "castBronze.add_element(\"Sn\", 0.10)\n",
    "castBronze.temperature = 293 #K\n",
    "castBronze.set_density(\"g/cm3\", 8.77)\n",
    "\n",
    "#Creating nitrogen\n",
    "nitrogen = openmc.Material(2, \"nitrogen\")\n",
    "nitrogen.set_density(\"g/cm3\", 0.0012506)\n",
    "nitrogen.add_element('N', 2.0)\n",
    "nitrogen.temperature = 293 #K\n",
    "\n",
    "#Creating oxygen\n",
    "oxygen = openmc.Material(3, \"oxygen\")\n",
    "oxygen.set_density(\"g/cm3\", 0.001429)\n",
    "oxygen.add_element('O', 2.0)\n",
    "oxygen.temperature = 293 #K\n",
    "\n",
    "#Creatiing argon\n",
    "argon = openmc.Material(4, \"argon\")\n",
    "argon.set_density(\"g/cm3\", 0.0017837)\n",
    "argon.add_element('Ar', 1.0)\n",
    "argon.temperature = 293 #K\n",
    "\n",
    "#Creating Carbon Dioxide\n",
    "CO2 = openmc.Material(5, \"CO2\")\n",
    "CO2.set_density(\"g/cm3\", 0.001976)\n",
    "CO2.add_element('C', 1.0)\n",
    "CO2.add_element('O', 2)\n",
    "CO2.temperature = 293 #K\n",
    "\n",
    "#Adds all gasses into one material\n",
    "airmix = openmc.Material.mix_materials([nitrogen, oxygen, argon, CO2], [0.7808, 0.2095, 0.0093, 0.0004], 'ao')\n",
    "\n",
    "#Exporting materials\n",
    "mats = openmc.Materials([castBronze, nitrogen, oxygen, argon, CO2, airmix])\n",
    "mats.export_to_xml()\n",
    "\n",
    "# GEOMETRY\n",
    "cylWidth = 1.01092\n",
    "cylDiam = 2.51968\n",
    "\n",
    "#Making cylinder\n",
    "cylinderSurface = openmc.ZCylinder(r=(cylDiam/2), boundary_type='transmission')\n",
    "upperSurface = openmc.ZPlane(z0=(20+cylWidth), boundary_type='transmission')\n",
    "lowerSurface = openmc.ZPlane(z0=20, boundary_type='transmission')\n",
    "\n",
    "#Creates cylinder region\n",
    "cylinderInside = -cylinderSurface & -upperSurface & +lowerSurface\n",
    "\n",
    "# Create cells, mapping materials to cells\n",
    "cylinder = openmc.Cell(name='cylinder')\n",
    "cylinder.fill = castBronze\n",
    "cylinder.region = cylinderInside\n",
    "\n",
    "#Create cylinder behind to count flux\n",
    "BtestcylinderSurface = openmc.ZCylinder(r=(cylDiam/2), boundary_type='transmission')\n",
    "BtestupperSurface = openmc.ZPlane(z0=20, boundary_type='transmission')\n",
    "BtestlowerSurface = openmc.ZPlane(z0=19, boundary_type='transmission')\n",
    "\n",
    "BtestcylinderInside = -BtestcylinderSurface & -BtestupperSurface & +BtestlowerSurface\n",
    "\n",
    "# Create cells, mapping materials to cells\n",
    "Btestcylinder = openmc.Cell(name='Btestcylinder')\n",
    "Btestcylinder.fill = airmix\n",
    "Btestcylinder.region = BtestcylinderInside\n",
    "\n",
    "#Create cylinder in front to count flux\n",
    "TtestcylinderSurface = openmc.ZCylinder(r=(cylDiam/2), boundary_type='transmission')\n",
    "TtestupperSurface = openmc.ZPlane(z0=(21+cylWidth), boundary_type='transmission')\n",
    "TtestlowerSurface = openmc.ZPlane(z0=(20+cylWidth), boundary_type='transmission')\n",
    "\n",
    "TtestcylinderInside = -TtestcylinderSurface & -TtestupperSurface & +TtestlowerSurface\n",
    "\n",
    "# Create cells, mapping materials to cells\n",
    "Ttestcylinder = openmc.Cell(name='Ttestcylinder')\n",
    "Ttestcylinder.fill = airmix\n",
    "Ttestcylinder.region = TtestcylinderInside\n",
    "\n",
    "#Create a box of air for cylinders to be in\n",
    "xmin = openmc.XPlane(-60.0, boundary_type='vacuum')\n",
    "xmax = openmc.XPlane(60.0, boundary_type='vacuum')\n",
    "ymin = openmc.YPlane(-60.0, boundary_type='vacuum')\n",
    "ymax = openmc.YPlane(60.0, boundary_type='vacuum')\n",
    "zmin = openmc.ZPlane(-60.0, boundary_type='vacuum')\n",
    "zmax = openmc.ZPlane(60.0, boundary_type='vacuum')\n",
    "box = +xmin & -xmax & +ymin & -ymax & +zmin & -zmax\n",
    "\n",
    "#Making regions and assigning to cells for box\n",
    "air_region = ~cylinderInside & ~BtestcylinderInside & ~TtestcylinderInside & box\n",
    "air = openmc.Cell(name='air')\n",
    "air.fill = airmix\n",
    "air.region = air_region\n",
    "\n",
    "#Create universe, export geometry\n",
    "root_universe = openmc.Universe(cells=[cylinder, Ttestcylinder, Btestcylinder, air])\n",
    "geom = openmc.Geometry(root_universe)                    \n",
    "geom.export_to_xml()  \n",
    "\n",
    "#Graphs Geometry\n",
    "root_universe.plot(width=(7,7), basis='yz', origin=(0,0,19), color_by=\"cell\", colors={Ttestcylinder: 'mistyrose', Btestcylinder: 'mistyrose', cylinder: 'goldenrod', air: 'lightsteelblue'})\n",
    "\n",
    "#SETTINGS\n",
    "settings = openmc.Settings()\n",
    "settings.particles = 10000\n",
    "settings.batches = 10\n",
    "settings.inactive = 0\n",
    "\n",
    "#Creates source\n",
    "r = openmc.stats.PowerLaw(0.0, 2, 1.0)\n",
    "phi = openmc.stats.Uniform(0.0, 2 * np.pi)\n",
    "z = openmc.stats.Discrete([1.0], [0.0])\n",
    "spatial_dist = openmc.stats.CylindricalIndependent(r, phi, z)\n",
    "source = openmc.Source(space=spatial_dist)\n",
    "\n",
    "#Sets source type and energy\n",
    "source.angle = openmc.stats.Monodirectional([0, 0, 1])\n",
    "source.energy = openmc.stats.Discrete([0.025], [1.0])\n",
    "source.strength = 1\n",
    "settings.source = source\n",
    "\n",
    "#Exports settings\n",
    "settings.run_mode = 'fixed source'\n",
    "settings.export_to_xml()\n",
    "\n",
    "#TALLIES\n",
    "#Creates a filter on the top of the sample\n",
    "Ttest_filter = openmc.CellFilter(Ttestcylinder)\n",
    "\n",
    "#Counts the amount of initial flux\n",
    "Ttest_tally = openmc.Tally(name=\"flux of Top test cylinder\")\n",
    "Ttest_tally.filters = [Ttest_filter]\n",
    "Ttest_tally.scores = ['flux']\n",
    "\n",
    "#Creates a filter on the bottom of the sample\n",
    "Btest_filter = openmc.CellFilter(Btestcylinder)\n",
    "\n",
    "#Counts the amount of  flux transmitted\n",
    "Btest_tally = openmc.Tally(name=\"flux of Bottom test cylinder\")\n",
    "Btest_tally.filters = [Btest_filter]\n",
    "Btest_tally.scores = ['flux']\n",
    "\n",
    "#Export Tallies\n",
    "tallies = openmc.Tallies([Ttest_tally, Btest_tally])\n",
    "tallies.export_to_xml()\n",
    "\n",
    "#Run Code\n",
    "model=openmc.model.Model(geom, mats, settings, tallies)\n",
    "openmc.run()\n",
    "\n",
    "#Prints results\n",
    "!cat tallies.out"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}