{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# TRIGA Mark II Reactor" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import os\n", "\n", "import numpy as np\n", "import openmc\n", "\n", "# remove previous model\n", "if os.path.exists('settings.xml'):\n", " os.system(\"rm *.xml\")\n", " os.system(\"rm summary.h5\")\n", " os.system(\"rm statepoint*\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Materials" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Nomenclature for thermal libraries: https://docs.openmc.org/en/stable/_modules/openmc/data/thermal.html\n", "\n", "T_ref_w = 294\n", "T_ref_f = 294\n", "T_ref_struct = 294\n", "\n", "# Water\n", "Water = openmc.Material(name='Water')\n", "Water.set_density('g/cm3',0.9985)\n", "Water.add_nuclide('H1', 2, 'ao')\n", "Water.add_nuclide('O16', 1, 'ao')\n", "Water.add_s_alpha_beta('c_H_in_H2O')\n", "Water.temperature = T_ref_w # (K)\n", "\n", "# Fuel 101 @ Pavia \n", "Fuel101 = openmc.Material(name='Fuel 101')\n", "Fuel101.set_density('g/cm3', 6.3)\n", "Fuel101.add_nuclide('U235', 1.58, 'wo')\n", "Fuel101.add_nuclide('U238', 6.42, 'wo')\n", "Fuel101.add_nuclide('H1', 1, 'wo')\n", "Fuel101.add_nuclide('Zr90', 46.1, 'wo')\n", "Fuel101.add_nuclide('Zr91', 10.2, 'wo')\n", "Fuel101.add_nuclide('Zr92', 15.7, 'wo')\n", "Fuel101.add_nuclide('Zr94', 16.3, 'wo')\n", "Fuel101.add_nuclide('Zr96', 2.68, 'wo')\n", "Fuel101.add_s_alpha_beta('c_H_in_ZrH')\n", "#Fuel101.add_s_alpha_beta('c_Zr_in_ZrH')\n", "Fuel101.temperature = T_ref_f\n", "\n", "# Fuel 103 @ Pavia \n", "Fuel103 = openmc.Material(name='Fuel 103')\n", "Fuel103.set_density('g/cm3', 5.95)\n", "Fuel103.add_nuclide('U235', 1.58, 'wo')\n", "Fuel103.add_nuclide('U238', 6.42, 'wo')\n", "Fuel103.add_nuclide('H1', 1.59, 'wo')\n", "Fuel103.add_nuclide('Zr90', 45.8, 'wo')\n", "Fuel103.add_nuclide('Zr91', 10.1, 'wo')\n", "Fuel103.add_nuclide('Zr92', 15.6, 'wo')\n", "Fuel103.add_nuclide('Zr94', 16.2, 'wo')\n", "Fuel103.add_nuclide('Zr96', 2.66, 'wo')\n", "Fuel103.add_s_alpha_beta('c_H_in_ZrH')\n", "#Fuel103.add_s_alpha_beta('c_Zr_in_ZrH')\n", "Fuel103.temperature= T_ref_f\n", "\n", "# Graphite\n", "Graphite = openmc.Material(name='Graphite')\n", "Graphite.set_density('g/cm3', 1.675)\n", "Graphite.add_element('C', 0.99, 'wo')\n", "Graphite.add_nuclide('Fe54', 0.00045164, 'wo')\n", "Graphite.add_nuclide('Fe56', 0.00735212, 'wo')\n", "Graphite.add_nuclide('Fe57', 0.00017283, 'wo')\n", "Graphite.add_nuclide('Fe58', 0.00002340, 'wo')\n", "Graphite.add_nuclide('O16', 0.0013, 'wo')\n", "Graphite.add_nuclide('Ti46', 7.92E-05, 'wo')\n", "Graphite.add_nuclide('Ti47', 7.30E-05, 'wo')\n", "Graphite.add_nuclide('Ti48', 7.38E-04, 'wo')\n", "Graphite.add_nuclide('Ti49', 5.53E-05, 'wo')\n", "Graphite.add_nuclide('Ti50', 5.40E-05, 'wo')\n", "Graphite.add_nuclide('Na23', 0.0001, 'wo')\n", "Graphite.add_nuclide('Mg24', 7.795E-05, 'wo')\n", "Graphite.add_nuclide('Mg25', 1.028E-05, 'wo')\n", "Graphite.add_nuclide('Mg26', 1.177E-05, 'wo')\n", "Graphite.add_nuclide('Cu63', 6.85E-04, 'wo')\n", "Graphite.add_nuclide('Cu65', 3.15E-04, 'wo')\n", "Graphite.add_nuclide('Co59', 0.00005, 'wo')\n", "Graphite.add_nuclide('H1', 0.000037, 'wo')\n", "Graphite.add_nuclide('F19', 0.000025, 'wo')\n", "Graphite.add_element('Cl', 0.000025, 'wo')\n", "Graphite.add_nuclide('B11', 0.0000005, 'wo')\n", "Graphite.add_s_alpha_beta('c_Graphite')\n", "Graphite.temperature = T_ref_struct\n", "\n", "# Cladding Al (from Chiesa PhD Thesis) \n", "CladdingAl = openmc.Material(name='Cladding Al') \n", "CladdingAl.set_density('g/cm3' , 2.713) \n", "CladdingAl.add_nuclide('Al27', 0.9957, 'wo')\n", "CladdingAl.add_nuclide('Cu63', 0.0006917, 'wo')\n", "CladdingAl.add_nuclide('Cu65', 0.0003083, 'wo')\n", "CladdingAl.temperature = T_ref_struct\n", "\n", "# Cladding SS\n", "CladdingSS = openmc.Material(name='Cladding SS') \n", "CladdingSS.set_density('g/cm3' , 8.03) \n", "CladdingSS.add_nuclide('Fe54' , 0.039980, 'wo')\n", "CladdingSS.add_nuclide('Fe56' , 0.650820, 'wo')\n", "CladdingSS.add_nuclide('Fe57' , 0.015299, 'wo')\n", "CladdingSS.add_nuclide('Fe58' , 0.002072, 'wo')\n", "CladdingSS.add_nuclide('Cr50' , 7.72E-03, 'wo')\n", "CladdingSS.add_nuclide('Cr52' , 1.49E-01, 'wo')\n", "CladdingSS.add_nuclide('Cr53' , 1.69E-02, 'wo')\n", "CladdingSS.add_nuclide('Cr54' , 4.19E-03, 'wo')\n", "CladdingSS.add_nuclide('Ni58' , 7.42E-02, 'wo')\n", "CladdingSS.add_nuclide('Ni60' , 2.86E-02, 'wo')\n", "CladdingSS.add_nuclide('Ni61' , 1.24E-03, 'wo')\n", "CladdingSS.add_nuclide('Ni62' , 3.95E-03, 'wo')\n", "CladdingSS.add_nuclide('Ni64' , 1.01E-03, 'wo')\n", "CladdingSS.add_nuclide('Si28' , 5.02E-03, 'wo')\n", "CladdingSS.add_nuclide('Si29' , 2.55E-04, 'wo')\n", "CladdingSS.add_nuclide('Si30' , 1.68E-04, 'wo')\n", "CladdingSS.temperature = T_ref_struct\n", "\n", "# Cladding for Control Rods (Aluminum)\n", "CladdingCR = openmc.Material(name='Cladding CR') \n", "CladdingCR.set_density('g/cm3' , 2.71)\n", "CladdingCR.add_nuclide('Al27', 0.9795, 'wo')\n", "CladdingCR.add_nuclide('Cu63', 0.00172925, 'wo')\n", "CladdingCR.add_nuclide('Cu65', 0.00077075, 'wo')\n", "CladdingCR.add_nuclide('Cr52', 0.0020, 'wo')\n", "CladdingCR.add_nuclide('Mg24' , 7.80E-03 , 'wo')\n", "CladdingCR.add_nuclide('Mg25' , 1.03E-03 , 'wo')\n", "CladdingCR.add_nuclide('Mg26' , 1.18E-03 , 'wo')\n", "CladdingCR.add_nuclide('Si28', 0.006, 'wo')\n", "CladdingCR.temperature = T_ref_struct \n", "\n", "# Zirconium thin rod\n", "Zirconium = openmc.Material(name='Zirconium')\n", "Zirconium.set_density('g/cm3', 6.52)\n", "Zirconium.add_nuclide('Zr90', 0.5145, 'ao')\n", "Zirconium.add_nuclide('Zr91', 0.1122, 'ao')\n", "Zirconium.add_nuclide('Zr92', 0.1715, 'ao')\n", "Zirconium.add_nuclide('Zr94', 0.1738, 'ao')\n", "Zirconium.add_nuclide('Zr96', 0.028, 'ao')\n", "Zirconium.temperature = T_ref_f \n", "\n", "# Poison disk\n", "Samarium = openmc.Material(name='Sm poison disk')\n", "Samarium.set_density('g/cm3', 2.4203) \n", "Samarium.add_nuclide('Sm147', 0.001035, 'wo')\n", "Samarium.add_nuclide('Sm148', 0.000782, 'wo')\n", "Samarium.add_nuclide('Sm149', 0.000968, 'wo')\n", "Samarium.add_nuclide('Sm150', 0.000520, 'wo')\n", "Samarium.add_nuclide('Sm152', 0.001911, 'wo')\n", "Samarium.add_nuclide('Sm154', 0.001647, 'wo')\n", "Samarium.add_nuclide('O16', 0.468019, 'wo')\n", "Samarium.add_nuclide('Al27', 0.524910, 'wo')\n", "Samarium.temperature = T_ref_f\n", "\n", "# Vacuum\n", "Vacuum = openmc.Material(name='Vacuum')\n", "Vacuum.set_density('g/cm3', 0.0012041)\n", "Vacuum.add_element('C', 0.000125, 'wo')\n", "Vacuum.add_nuclide('N14', 0.6869, 'wo')\n", "Vacuum.add_nuclide('O16', 0.301248, 'wo')\n", "Vacuum.temperature = 294\n", "\n", "# SHIM - REG\n", "ControlRodB4C = openmc.Material(name='Control rods (SHIM, REG)')\n", "ControlRodB4C.set_density('g/cm3' , 2.52)\n", "ControlRodB4C.add_nuclide('B11' , 64.2, 'ao')\n", "ControlRodB4C.add_nuclide('B10' , 15.8, 'ao')\n", "ControlRodB4C.add_element('C' , 20, 'ao')\n", "ControlRodB4C.add_s_alpha_beta('c_Graphite')\n", "ControlRodB4C.temperature = T_ref_struct\n", "\n", "# TRANS\n", "ControlRodBorated = openmc.Material(name='Control rods (TRANS)')\n", "ControlRodBorated.set_density('g/cm3' , 2.23)\n", "ControlRodBorated.add_nuclide('B11' , 21.7, 'ao')\n", "ControlRodBorated.add_nuclide('B10' , 5.35, 'ao')\n", "ControlRodBorated.add_element('C' , 72.3, 'ao')\n", "ControlRodBorated.add_s_alpha_beta('c_Graphite')\n", "ControlRodBorated.temperature = T_ref_struct\n", "\n", "\n", "# Instantiate a Materials collection and export to xml\n", "materials_file = openmc.Materials([CladdingAl, CladdingSS, CladdingCR, Graphite, Fuel101, Fuel103, Water, Vacuum, ControlRodB4C, ControlRodBorated, Samarium, Zirconium ])\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Geometry: planes and surfaces" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Control rod position: 0 = all in, 47 = all out\n", "z0_reg_value = 45 # cm\n", "z0_shim_value = 45 # cm\n", "z0_trans_value = 45 # cm" ] }, { "cell_type": "code", "execution_count": 27, "metadata": {}, "outputs": [], "source": [ "# Planes, cylinders and stuff\n", "\n", "r_fuel_101 = openmc.ZCylinder(r=1.791) #fuel radius\n", "r_clad_101 = openmc.ZCylinder(r=1.88) #cladding radius\n", "bot_active_101 = openmc.ZPlane(z0=10.20) #start active region 101\n", "top_active_101 = openmc.ZPlane(z0=45.76) #end active region 101\n", "Sm_bot_101 = openmc.ZPlane(z0=10.07) # samarium disk bot end \n", "Sm_top_101 = openmc.ZPlane(z0=45.89) # samarium disk top end \n", "reflector_bot_101 = openmc.ZPlane(z0=0) #bottom 101 reflector\n", "reflector_top_101 = openmc.ZPlane(z0=55.83) #top 101 reflector\n", "clad_bot_101 = openmc.ZPlane(z0=-1.27) #end bottom cladding 101\n", "clad_top_101 = openmc.ZPlane(z0=57.1) #end top cladding 101\n", "plug_top_101 = openmc.ZPlane(z0=64.26) #end plug 101\n", "plug_radius_101 = openmc.ZCylinder(r=0.395)\n", "\n", "\n", "# Geometry definitions for the Fuel Rod 103\n", "r_fuel_inner_103 = openmc.ZCylinder(r=0.285) #fuel radius\n", "r_fuel_outer_103 = openmc.ZCylinder(r=1.82) #fuel radius\n", "r_clad_103 = openmc.ZCylinder(r=1.88) #cladding radius\n", "bot_active_103 = openmc.ZPlane(z0=8.81) #start active region 103\n", "top_active_103 = openmc.ZPlane(z0=46.91) #end active region 103\n", "reflector_bot_103 = openmc.ZPlane(z0=0) #bottom 103 reflector\n", "reflector_top_103 = openmc.ZPlane(z0=55.65) #top 103 reflector\n", "clad_bot_103 = openmc.ZPlane(z0=-1.39) #end bottom cladding 103\n", "clad_top_103 = openmc.ZPlane(z0=57.31) #end top cladding 103\n", "plug_top_103 = openmc.ZPlane(z0=64.47) #end plug 103\n", "plug_radius_103 = openmc.ZCylinder(r=0.395)\n", "\n", "# Geometry for dummy element (graphite)\n", "r_dummy_graphite = openmc.ZCylinder(r=1.791)\n", "r_dummy_clad = openmc.ZCylinder(r=1.88)\n", "dummy_bot = openmc.ZPlane(z0=0) #bottom 101 reflector\n", "dummy_top = openmc.ZPlane(z0=55.83) #top 101 reflector\n", "dummy_bot_clad = openmc.ZPlane(z0=-1.27) #end bottom cladding 101\n", "dummy_top_clad = openmc.ZPlane(z0=57.1) #end top cladding 101\n", "\n", "# Geometry for central channel (with Al case for void feedback measurements)\n", "r_central_channel = openmc.ZCylinder(r=1.69)\n", "r_central_channel_clad = openmc.ZCylinder(r=1.905)\n", "\n", "r_Al_out = 3/2 # (cm)\n", "r_Al_in = 2.4/2 # (cm)\n", "l_Al_out = 28 # (cm)\n", "l_Al_in = 26.5258238 # (cm)\n", "midplane = 27.86\n", "plane_mat_separation = l_Al_in - 0.01 \n", "\n", "s_void_r_out = openmc.ZCylinder(r = r_Al_out)\n", "s_void_r_in = openmc.ZCylinder(r = r_Al_in)\n", "s_void_bot_in = openmc.ZPlane(z0 = midplane - l_Al_in/2)\n", "s_void_bot_out = openmc.ZPlane(z0 = midplane - l_Al_out/2)\n", "s_void_sep = openmc.ZPlane(z0 = midplane - l_Al_in/2 + plane_mat_separation)\n", "s_void_top_in = openmc.ZPlane(z0 = midplane + l_Al_in/2)\n", "s_void_top_out = openmc.ZPlane(z0 = midplane + l_Al_out/2)\n", "\n", "# Geometry for rabbit channel \n", "r_rabbit_inner = openmc.ZCylinder(r=1.804)\n", "r_rabbit_outer = openmc.ZCylinder(r=1.88)\n", "\n", "# Geometry for source\n", "r_source = openmc.ZCylinder(r=1.88)\n", "\n", "# Geometry for SHIM\n", "r_shim_inner = openmc.ZCylinder(r=1.425)\n", "r_shim_outer = openmc.ZCylinder(r=1.59)\n", "r_shim_bot_plug = openmc.ZCylinder(r=0.795)\n", "\n", "# Geometry for REG\n", "r_reg_inner = openmc.ZCylinder(r=0.965)\n", "r_reg_outer = openmc.ZCylinder(r=1.11)\n", "r_reg_bot_plug = openmc.ZCylinder(r=0.555)\n", "\n", "# Geometry for TRANS\n", "r_trans_inner = openmc.ZCylinder(r=1.105)\n", "r_trans_outer = openmc.ZCylinder(r=1.27)\n", "r_trans_bot_plug = openmc.ZCylinder(r=0.635)\n", "\n", "# CR heigths: length=45.47 cm, the starting position is zero.\n", "# I define three quantities, z_reg, z_shim and z_trans\n", "# such that they can be moved\n", "\n", "# CR position from their 0\n", "z0_reg= z0_reg_value\n", "z0_shim = z0_shim_value\n", "z0_trans = z0_trans_value\n", "\n", "# CALIBRAZIONE: cambiare le altezze delle barre di controllo per calcolare il k\n", "# (0 = rod tutta inserita, )\n", "# valori reg: 0.27 , 4.57 , 9.85 , 13.77, 17.69, 21.55, 24.77, 27.21, 38.1\n", "# valori shim: 19.92, 19.92, 19.38, 18.45, 16.93, 15.51, 13.93, 12.74, 11.81 \n", "# valori trans: sempre estratta \n", "\n", "# Model offset (z=+10.4 means that the rod is 'all in')\n", "offset= 10.40\n", "\n", "z0_reg = z0_reg + offset\n", "z0_shim = z0_shim + offset\n", "z0_trans = z0_trans + offset\n", "\n", "bot_reg = openmc.ZPlane(z0=z0_reg) #plane for the bottom reg\n", "bot_shim = openmc.ZPlane(z0=z0_shim) #plane for the bottom shim\n", "bot_trans = openmc.ZPlane(z0=z0_trans) #plane for the bottom trans\n", "\n", "bot_reg_clad = openmc.ZPlane(z0=z0_reg-0.8) #plane for the bottom clad\n", "bot_shim_clad = openmc.ZPlane(z0=z0_shim-0.8)\n", "bot_trans_clad = openmc.ZPlane(z0=z0_trans-0.8)\n", "\n", "bot_reg_plug = openmc.ZPlane(z0= z0_reg -1.6) #plane for the bottom of plug\n", "bot_shim_plug = openmc.ZPlane(z0= z0_shim -1.6) \n", "bot_trans_plug = openmc.ZPlane(z0= z0_trans -1.6)\n", "\n", "top_reg = openmc.ZPlane(z0=z0_reg+45.47) #plane for the bottom reg\n", "top_shim = openmc.ZPlane(z0=z0_shim+45.47) #plane for the bottom shim\n", "top_trans = openmc.ZPlane(z0=z0_trans+45.47) #plane for the bottom trans\n", "\n", "#top plug has width of 2.54 cm\n", "top_clad_reg = openmc.ZPlane(z0=z0_reg+45.47+2.54)\n", "top_clad_shim = openmc.ZPlane(z0=z0_shim+45.47+2.54)\n", "top_clad_trans = openmc.ZPlane(z0=z0_trans+45.47+2.54)\n", "\n", "top_CR = openmc.ZPlane(z0=45.47) #control rod length from 0\n", "start_position_CR = openmc.ZPlane(z0=0) #zero position for CR\n", "\n", "# Reflector dimensions planes\n", "bot_reflector = openmc.ZPlane(z0=1.27) #start of side reflector\n", "top_reflector = openmc.ZPlane(z0=57.17) #end of side reflector\n", "\n", "bot_clad_reflector = openmc.ZPlane(z0=0.0) #start of side reflector's clad\n", "top_clad_reflector = openmc.ZPlane(z0=59.04) #end of side reflector's clad\n", "\n", "#Grid at the top of the core\n", "top_grid_bot = openmc.ZPlane(z0=59.04)\n", "top_grid_top = openmc.ZPlane(z0=60.94)\n", "\n", "# Grid at the bottom of core\n", "bottom_grid_bot = openmc.ZPlane(z0=-6.995)\n", "bottom_grid_top = openmc.ZPlane(z0=-5.05) # true\n", "\n", "\n", "# Void in reflector (lazy susan)\n", "refl_void_inner = openmc.ZCylinder(r=30.5)\n", "refl_void_outer = openmc.ZCylinder(r=37)\n", "refl_void_b = openmc.ZPlane(z0=31.35)\n", "\n", "#IRRADIATION CHANNELS\n", "#--> Since OpenMC treat badly rotations, the channels' eqns have been defined as quadric (matlab script)\n", "# Also, a discriminator plane need to be defined since the cylinder is bi-directional. \n", "# Finally, a \"filling\" graphite region need to be defined to fill the water gap.\n", "\n", "channel_1_discriminator = openmc.YPlane(y0=0.0)\n", "channel_1_inner = openmc.Quadric(a=0.9963, b=0.0037, c=1, d=-0.1213, e=0, f=0, g=0, h=0, j=-42.7, k=396.5325) \n", "channel_1_outer = openmc.Quadric(a=0.9963, b=0.0037, c=1, d=-0.1213, e=0, f=0, g=0, h=0, j=-42.7, k=385.2625) \n", "\n", "channel_2_discriminator = openmc.YPlane(y0=0.0)\n", "channel_2_inner = openmc.Quadric(a=361, b=739.5735, c=1100.6, d=-1033.4, e=0, f=0, g=0, h=0, j=-46994, k=436413.1471) \n", "channel_2_outer = openmc.Quadric(a=361, b=739.5735, c=1100.6, d=-1033.4, e=0, f=0, g=0, h=0, j=-46994, k=424009.68) \n", "\n", "channel_3_discriminator = openmc.Plane(a=-19, b=+27.195, c=0, d=0) #3D plane which cut in half channel 2\n", "channel_3_inner = openmc.Quadric(a=67.1597, b=4.7337, c=71.8933, d=35.6602, e=0, f=0, g=4330.4, h=1149.7, j=-3069.8, k=98312, boundary_type = \"transmission\")\n", "channel_3_outer = openmc.Quadric(a=67.1597, b=4.7337, c=71.8933, d=35.6602, e=0, f=0, g=4330.4, h=1149.7, j=-3069.8, k=97502, boundary_type = \"transmission\")\n", "\n", "\n", "# Core delimitations\n", "core_radius = openmc.ZCylinder(r=22.3, boundary_type = 'transmission') #true\n", "inner_radius_reflector = openmc.ZCylinder(r=23.15)\n", "outer_radius_reflector = openmc.ZCylinder(r=53.35)\n", "s1 = openmc.ZCylinder(r=54.62, boundary_type='vacuum') #cylindrical surface delimiting the core\n", "end_domain_top = openmc.ZPlane(z0=125.7, boundary_type='vacuum')\n", "end_domain_bot = openmc.ZPlane(z0=-50, boundary_type='vacuum') \n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Geometry: cells and universes" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [], "source": [ "# Inline plot\n", "\n", "inline_plot = False\n", "cellid = 1" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Fuel element 101\n", "fuel_101_universe = openmc.Universe(name='UZrH Fuel 101 Universe')\n", "\n", "fuel_101_cell = openmc.Cell(name='Fuel 101', cell_id = cellid)\n", "fuel_101_cell.fill = Fuel101\n", "fuel_101_cell.region = -r_fuel_101 & +bot_active_101 & -top_active_101\n", "\n", "cellid += 1\n", "sm_top_101_cell = openmc.Cell(name='Sm disk top 101', cell_id = cellid)\n", "sm_top_101_cell.fill = Samarium\n", "sm_top_101_cell.region = -r_fuel_101 & +top_active_101 & -Sm_top_101\n", "\n", "cellid += 1\n", "graphite_top_101_cell = openmc.Cell(name='graphite top 101', cell_id = cellid)\n", "graphite_top_101_cell.fill = Graphite\n", "graphite_top_101_cell.region = -r_fuel_101 & +Sm_top_101 & -reflector_top_101\n", "\n", "cellid += 1\n", "sm_bot_101_cell = openmc.Cell(name='Sm disk bot 101', cell_id = cellid)\n", "sm_bot_101_cell.fill = Samarium\n", "sm_bot_101_cell.region = -r_fuel_101 & +Sm_bot_101 & -bot_active_101\n", "\n", "cellid += 1\n", "graphite_bot_101_cell = openmc.Cell(name='graphite bot 101', cell_id = cellid)\n", "graphite_bot_101_cell.fill = Graphite\n", "graphite_bot_101_cell.region = -r_fuel_101 & +reflector_bot_101 & -Sm_bot_101\n", "\n", "cellid += 1\n", "clad_101_cell = openmc.Cell(name='Al Cladding', cell_id = cellid)\n", "clad_101_cell.fill = CladdingAl\n", "clad_101_cell.region = -r_clad_101 & +r_fuel_101 & +reflector_bot_101 & -reflector_top_101\n", "\n", "cellid += 1\n", "clad_top_101_cell = openmc.Cell(name='Al Cladding top', cell_id = cellid)\n", "clad_top_101_cell.fill = CladdingAl\n", "clad_top_101_cell.region = -r_clad_101 & +reflector_top_101 & -clad_top_101\n", "\n", "cellid += 1\n", "clad_bot_101_cell = openmc.Cell(name='Al Cladding bot', cell_id = cellid)\n", "clad_bot_101_cell.fill = CladdingAl\n", "clad_bot_101_cell.region = -r_clad_101 & +clad_bot_101 & -reflector_bot_101\n", "\n", "cellid += 1\n", "plug_top_101_cell = openmc.Cell(name='top plug', cell_id = cellid)\n", "plug_top_101_cell.fill = CladdingAl\n", "plug_top_101_cell.region = -plug_radius_101 & +clad_top_101 & -plug_top_101\n", "\n", "cellid += 1\n", "plug_bot_101_cell = openmc.Cell(name='bot plug', cell_id = cellid)\n", "plug_bot_101_cell.fill = CladdingAl\n", "plug_bot_101_cell.region = -plug_radius_101 & +bottom_grid_top & -clad_bot_101\n", "\n", "cellid += 1\n", "plug_w_101_cell = openmc.Cell(name='top plug water', cell_id = cellid)\n", "plug_w_101_cell.fill = Water\n", "plug_w_101_cell.region = +plug_radius_101 & -r_clad_101 &+clad_top_101 & -plug_top_101\n", "\n", "cellid += 1\n", "water_top_fuel_101_cell=openmc.Cell(cell_id = cellid)\n", "water_top_fuel_101_cell.fill = Water\n", "water_top_fuel_101_cell.region = -r_clad_101 & +plug_top_101 &-end_domain_top\n", "\n", "cellid += 1\n", "water_bot_fuel_101_cell=openmc.Cell(cell_id = cellid)\n", "water_bot_fuel_101_cell.fill = Water\n", "water_bot_fuel_101_cell.region = -r_clad_101 & +end_domain_bot &-clad_bot_101\n", "\n", "\n", "fuel_101_universe.add_cell(fuel_101_cell)\n", "fuel_101_universe.add_cell(sm_top_101_cell)\n", "fuel_101_universe.add_cell(graphite_top_101_cell)\n", "fuel_101_universe.add_cell(sm_bot_101_cell)\n", "fuel_101_universe.add_cell(graphite_bot_101_cell)\n", "fuel_101_universe.add_cell(clad_101_cell)\n", "fuel_101_universe.add_cell(clad_top_101_cell)\n", "fuel_101_universe.add_cell(clad_bot_101_cell)\n", "fuel_101_universe.add_cell(plug_top_101_cell)\n", "fuel_101_universe.add_cell(plug_bot_101_cell)\n", "fuel_101_universe.add_cell(plug_w_101_cell)\n", "fuel_101_universe.add_cell(water_top_fuel_101_cell)\n", "fuel_101_universe.add_cell(water_bot_fuel_101_cell)\n", "\n", "if inline_plot:\n", " fuel_101_universe.plot(origin=(0,0,27), width=(6.0, 70.0), pixels=[400, 400], basis='xz', color_by='material')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Fuel element 103\n", "fuel_103_universe = openmc.Universe(name='UZrH Fuel 103 Universe')\n", "\n", "cellid += 1\n", "Zr_rod_103_cell = openmc.Cell(name='Zr thin rod', cell_id = cellid)\n", "Zr_rod_103_cell.fill = Zirconium\n", "Zr_rod_103_cell.region = -r_fuel_inner_103 & +bot_active_103 & -top_active_103\n", "\n", "cellid += 1\n", "fuel_103_cell = openmc.Cell(name='Fuel 103', cell_id = cellid)\n", "fuel_103_cell.fill = Fuel103\n", "fuel_103_cell.region = -r_fuel_outer_103 & +r_fuel_inner_103 & +bot_active_103 & -top_active_103\n", "\n", "cellid += 1\n", "graphite_top_103_cell = openmc.Cell(name='graphite top 103', cell_id = cellid)\n", "graphite_top_103_cell.fill = Graphite\n", "graphite_top_103_cell.region = -r_fuel_outer_103 & +top_active_103 & -reflector_top_103\n", "\n", "cellid += 1\n", "graphite_bot_103_cell = openmc.Cell(name='graphite bot 103', cell_id = cellid)\n", "graphite_bot_103_cell.fill = Graphite\n", "graphite_bot_103_cell.region = -r_fuel_outer_103 & +reflector_bot_103 & -bot_active_103\n", "\n", "cellid += 1\n", "clad_103_cell = openmc.Cell(name='SS Cladding 103', cell_id = cellid)\n", "clad_103_cell.fill = CladdingSS\n", "clad_103_cell.region = -r_clad_103 & +r_fuel_outer_103 & +reflector_bot_103 & -reflector_top_103\n", "\n", "cellid += 1\n", "clad_top_103_cell = openmc.Cell(name='SS Cladding top 103', cell_id = cellid)\n", "clad_top_103_cell.fill = CladdingSS\n", "clad_top_103_cell.region = -r_clad_103 & +reflector_top_103 & -clad_top_103\n", "\n", "cellid += 1\n", "clad_bot_103_cell = openmc.Cell(name='SS Cladding bot 103', cell_id = cellid)\n", "clad_bot_103_cell.fill = CladdingSS\n", "clad_bot_103_cell.region = -r_clad_103 & +clad_bot_103 & -reflector_bot_103\n", "\n", "cellid += 1\n", "plug_top_103_cell = openmc.Cell(name='top plug 103', cell_id = cellid)\n", "plug_top_103_cell.fill = CladdingSS\n", "plug_top_103_cell.region = -plug_radius_103 & +clad_top_103 & -plug_top_103\n", "\n", "cellid += 1\n", "plug_bot_103_cell = openmc.Cell(name='bot plug 103', cell_id = cellid)\n", "plug_bot_103_cell.fill = CladdingSS\n", "plug_bot_103_cell.region = -plug_radius_103 & +bottom_grid_top & -clad_bot_103\n", "\n", "cellid += 1\n", "plug_w_103_cell = openmc.Cell(name='top plug water 103', cell_id = cellid)\n", "plug_w_103_cell.fill = Water\n", "plug_w_103_cell.region = +plug_radius_103 & -r_clad_103 &+clad_top_103 & -plug_top_103\n", "\n", "cellid += 1\n", "water_top_fuel_103_cell=openmc.Cell(cell_id = cellid)\n", "water_top_fuel_103_cell.fill = Water\n", "water_top_fuel_103_cell.region = -r_clad_103 & +plug_top_103 &-end_domain_top\n", "\n", "cellid += 1\n", "water_bot_fuel_103_cell=openmc.Cell(cell_id = cellid)\n", "water_bot_fuel_103_cell.fill = Water\n", "water_bot_fuel_103_cell.region = -r_clad_103 & ~plug_bot_103_cell.region & +end_domain_bot &-clad_bot_103\n", "\n", "\n", "fuel_103_universe.add_cell(Zr_rod_103_cell)\n", "fuel_103_universe.add_cell(fuel_103_cell)\n", "fuel_103_universe.add_cell(graphite_top_103_cell)\n", "fuel_103_universe.add_cell(graphite_bot_103_cell)\n", "fuel_103_universe.add_cell(clad_103_cell)\n", "fuel_103_universe.add_cell(clad_top_103_cell)\n", "fuel_103_universe.add_cell(clad_bot_103_cell)\n", "fuel_103_universe.add_cell(plug_bot_103_cell)\n", "fuel_103_universe.add_cell(plug_top_103_cell)\n", "fuel_103_universe.add_cell(plug_w_103_cell)\n", "fuel_103_universe.add_cell(water_top_fuel_103_cell)\n", "fuel_103_universe.add_cell(water_bot_fuel_103_cell)\n", "\n", "if inline_plot:\n", " fuel_103_universe.plot(origin=(0,0,27), width=(6.0, 70.0), pixels=[400, 400], basis='xz', color_by='material')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Central Channel\n", "\n", "# Condition for filling the Al case for void feedback measurements. If void_water = True, fill the case with water. Else, air\n", "void_water = True \n", "\n", "if void_water: \n", " FillingMaterial = Water \n", "else:\n", " FillingMaterial = Vacuum\n", "\n", "central_channel_universe = openmc.Universe(name='CC Universe')\n", "\n", "cellid += 1\n", "central_channel_water = openmc.Cell(name='CC water below', cell_id = cellid)\n", "central_channel_water.fill = Water\n", "central_channel_water.region = -r_central_channel_clad & +end_domain_bot & -bottom_grid_bot\n", "\n", "cellid += 1\n", "central_channel_clad = openmc.Cell(name='CC clad', cell_id = cellid)\n", "central_channel_clad.fill = CladdingAl\n", "central_channel_clad.region = +r_central_channel & -r_central_channel_clad & +bottom_grid_bot & -end_domain_top\n", "\n", "cellid += 1\n", "central_channel_bot_cell = openmc.Cell(name='CC bot', cell_id = cellid)\n", "central_channel_bot_cell.fill = Vacuum\n", "central_channel_bot_cell.region = -r_central_channel & +bottom_grid_bot & -s_void_bot_out\n", "\n", "cellid += 1\n", "central_channel_top_cell = openmc.Cell(name='CC top', cell_id = cellid)\n", "central_channel_top_cell.fill = Vacuum\n", "central_channel_top_cell.region = -r_central_channel & +s_void_top_out & -end_domain_top\n", "\n", "\n", "# 6 celle per il barilotto di void:\n", "#0. and 1. = inner part, top and bottom, 2. = lateral clad, 3. and 4. = clad, top and bottom, 5. = lateral void\n", "cellid += 1\n", "void_cc_0 = openmc.Cell(name = 'void cc inner bot', cell_id = cellid)\n", "void_cc_0.fill = FillingMaterial\n", "void_cc_0.region = -s_void_r_in & +s_void_bot_in & -s_void_sep\n", "\n", "cellid += 1\n", "void_cc_1 = openmc.Cell(name = 'void cc inner top', cell_id = cellid)\n", "void_cc_1.fill = Vacuum \n", "void_cc_1.region = -s_void_r_in & +s_void_sep & -s_void_top_in\n", "\n", "cellid += 1\n", "void_cc_2 = openmc.Cell(name = 'void lateral clad', cell_id = cellid)\n", "void_cc_2.fill = CladdingAl\n", "void_cc_2.region = +s_void_r_in & -s_void_r_out & +s_void_bot_in & -s_void_top_in\n", "\n", "cellid += 1\n", "void_cc_3 = openmc.Cell(name = 'void bot clad', cell_id = cellid)\n", "void_cc_3.fill = CladdingAl\n", "void_cc_3.region = +s_void_bot_out & -s_void_bot_in & -s_void_r_out\n", "\n", "cellid += 1\n", "void_cc_4 = openmc.Cell(name = 'void top vlad', cell_id = cellid)\n", "void_cc_4.fill = CladdingAl\n", "void_cc_4.region = +s_void_top_in & -s_void_top_out & -s_void_r_out\n", "\n", "cellid += 1\n", "void_cc_5 = openmc.Cell(name = 'void lateral vacuum', cell_id = cellid)\n", "void_cc_5.fill = Vacuum\n", "void_cc_5.region = +s_void_bot_out & -s_void_top_out & +s_void_r_out & -r_central_channel\n", "\n", "\n", "central_channel_universe.add_cell(central_channel_bot_cell)\n", "central_channel_universe.add_cell(central_channel_top_cell)\n", "central_channel_universe.add_cell(central_channel_clad)\n", "central_channel_universe.add_cell(central_channel_water)\n", "central_channel_universe.add_cell(void_cc_0)\n", "central_channel_universe.add_cell(void_cc_1)\n", "central_channel_universe.add_cell(void_cc_2)\n", "central_channel_universe.add_cell(void_cc_3)\n", "central_channel_universe.add_cell(void_cc_4)\n", "central_channel_universe.add_cell(void_cc_5)\n", "\n", "\n", "if inline_plot:\n", " central_channel_universe.plot(origin=(0,0,27), width=(6.0, 90.0), pixels=[400, 400], basis='xz', color_by='material')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# TRANS rod\n", "transient_universe = openmc.Universe(name='TRANS Universe')\n", "\n", "cellid += 1\n", "TRANS_cell = openmc.Cell(name='TRANS Cell', cell_id = cellid)\n", "TRANS_cell.fill = ControlRodBorated #true\n", "TRANS_cell.region = -r_trans_inner & +bot_trans & -top_trans \n", "\n", "cellid += 1\n", "TRANS_clad = openmc.Cell(name='TRANS Clad Cell', cell_id = cellid)\n", "TRANS_clad.fill = CladdingCR\n", "TRANS_clad.region = -r_trans_outer & +r_trans_inner & +bot_trans & -top_trans\n", "\n", "cellid += 1\n", "TRANS_bot_clad_cell = openmc.Cell(name='TRANS clad bot cell', cell_id = cellid)\n", "TRANS_bot_clad_cell.fill = CladdingCR\n", "TRANS_bot_clad_cell.region = -r_trans_outer & +bot_trans_clad & -bot_trans\n", "\n", "cellid += 1\n", "TRANS_bot_plug_cell = openmc.Cell(name='TRANS clad bot plug', cell_id = cellid)\n", "TRANS_bot_plug_cell.fill = CladdingCR\n", "TRANS_bot_plug_cell.region = -r_trans_bot_plug & +bot_trans_plug & -bot_trans_clad\n", "\n", "cellid += 1\n", "TRANS_top_clad_cell = openmc.Cell(name='TRANS clad top cell', cell_id = cellid)\n", "TRANS_top_clad_cell.fill = CladdingCR\n", "TRANS_top_clad_cell.region = -r_trans_outer & +top_trans & -top_clad_trans\n", "\n", "cellid += 1\n", "TRANS_water_bot = openmc.Cell(name='TRANS water Cell', cell_id = cellid)\n", "TRANS_water_bot.fill = Water\n", "TRANS_water_bot.region = -r_clad_103 & +bottom_grid_top & -bot_trans_plug\n", "\n", "cellid += 1\n", "TRANS_water_bot_plug = openmc.Cell(name='TRANS water Cell bot plug', cell_id = cellid)\n", "TRANS_water_bot_plug.fill = Water\n", "TRANS_water_bot_plug.region = -r_clad_103 & +r_trans_bot_plug & -bot_trans_clad & +bot_trans_plug\n", "\n", "cellid += 1\n", "TRANS_water_top = openmc.Cell(name='TRANS water Cell', cell_id = cellid)\n", "TRANS_water_top.fill = Water\n", "TRANS_water_top.region = -r_clad_103 & +top_clad_trans & -end_domain_top\n", "\n", "cellid += 1\n", "TRANS_water_fill = openmc.Cell(name='TRANS water Cell for filling lattice', cell_id = cellid)\n", "TRANS_water_fill.fill = Water\n", "TRANS_water_fill.region = +r_trans_outer & -r_clad_103 & +bot_trans_clad & -top_clad_trans\n", "\n", "\n", "transient_universe.add_cell(TRANS_cell)\n", "transient_universe.add_cell(TRANS_clad)\n", "transient_universe.add_cell(TRANS_water_top)\n", "transient_universe.add_cell(TRANS_water_bot)\n", "transient_universe.add_cell(TRANS_water_bot_plug)\n", "transient_universe.add_cell(TRANS_water_fill)\n", "transient_universe.add_cell(TRANS_top_clad_cell)\n", "transient_universe.add_cell(TRANS_bot_clad_cell)\n", "transient_universe.add_cell(TRANS_bot_plug_cell)\n", "\n", "if inline_plot:\n", " transient_universe.plot(origin=(0,0,z0_trans), width=(6.0, 1.5*z0_trans), pixels=[400, 400], basis='xz', color_by='material')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# REG rod\n", "reg_universe = openmc.Universe(name='REG Universe')\n", "\n", "cellid += 1 \n", "REG_cell = openmc.Cell(name='REG Cell', cell_id = cellid)\n", "REG_cell.fill = ControlRodB4C\n", "REG_cell.region = -r_reg_inner & +bot_reg & -top_reg \n", "\n", "cellid += 1\n", "REG_clad = openmc.Cell(name='REG Clad Cell', cell_id = cellid)\n", "REG_clad.fill = CladdingCR\n", "REG_clad.region = -r_reg_outer & +r_reg_inner & +bot_reg & -top_reg\n", "\n", "cellid += 1\n", "REG_bot_clad_cell = openmc.Cell(name='REG clad bot cell', cell_id = cellid)\n", "REG_bot_clad_cell.fill = CladdingCR\n", "REG_bot_clad_cell.region = -r_reg_outer & +bot_reg_clad & -bot_reg\n", "\n", "cellid += 1\n", "REG_bot_plug_cell = openmc.Cell(name='REG clad bot plug', cell_id = cellid)\n", "REG_bot_plug_cell.fill = CladdingCR\n", "REG_bot_plug_cell.region = -r_reg_bot_plug & +bot_reg_plug & -bot_reg_clad\n", "\n", "cellid += 1\n", "REG_top_clad_cell = openmc.Cell(name='REG clad top cell', cell_id = cellid)\n", "REG_top_clad_cell.fill = CladdingCR\n", "REG_top_clad_cell.region = -r_reg_outer & +top_reg & -top_clad_reg\n", "\n", "cellid += 1\n", "REG_water_top = openmc.Cell(name='REG water Cell', cell_id = cellid)\n", "REG_water_top.fill = Water\n", "REG_water_top.region = -r_clad_103 & +top_clad_reg & -end_domain_top\n", "\n", "cellid += 1\n", "REG_water_bot = openmc.Cell(name='REG water Cell', cell_id = cellid)\n", "REG_water_bot.fill = Water\n", "REG_water_bot.region = -r_clad_103 & +bottom_grid_top & -bot_reg_plug\n", "\n", "cellid += 1\n", "REG_water_bot_plug = openmc.Cell(name='REG water Cell plug', cell_id = cellid)\n", "REG_water_bot_plug.fill = Water\n", "REG_water_bot_plug.region = -r_clad_103 & +r_reg_bot_plug & -bot_reg_clad & +bot_reg_plug\n", "\n", "cellid += 1\n", "REG_water_fill = openmc.Cell(name='REG water Cell for filling lattice', cell_id = cellid)\n", "REG_water_fill.fill = Water\n", "REG_water_fill.region = +r_reg_outer & -r_clad_103 & +bot_reg_clad & -top_clad_reg\n", "\n", "\n", "reg_universe.add_cell(REG_cell)\n", "reg_universe.add_cell(REG_clad)\n", "reg_universe.add_cell(REG_water_top)\n", "reg_universe.add_cell(REG_water_bot)\n", "reg_universe.add_cell(REG_water_bot_plug)\n", "reg_universe.add_cell(REG_water_fill)\n", "reg_universe.add_cell(REG_top_clad_cell)\n", "reg_universe.add_cell(REG_bot_clad_cell)\n", "reg_universe.add_cell(REG_bot_plug_cell)\n", "\n", "\n", "if inline_plot:\n", " reg_universe.plot(origin=(0,0,z0_reg), width=(6.0, 1.5*z0_reg), pixels=[400, 400], basis='xz', color_by='material')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# SHIM rod\n", "shim_universe = openmc.Universe(name='SHIM Universe')\n", "\n", "cellid += 1\n", "SHIM_cell = openmc.Cell(name='SHIM Cell', cell_id=cellid)\n", "SHIM_cell.fill = ControlRodB4C\n", "SHIM_cell.region = -r_shim_inner & +bot_shim & -top_shim \n", "\n", "cellid += 1\n", "SHIM_clad = openmc.Cell(name='SHIM Clad Cell', cell_id=cellid)\n", "SHIM_clad.fill = CladdingCR\n", "SHIM_clad.region = -r_shim_outer & +r_shim_inner & +bot_shim & -top_shim\n", "\n", "cellid += 1\n", "SHIM_bot_clad_cell = openmc.Cell(name='SHIM clad bot cell', cell_id=cellid)\n", "SHIM_bot_clad_cell.fill = CladdingCR\n", "SHIM_bot_clad_cell.region = -r_shim_outer & +bot_shim_clad & -bot_shim\n", "\n", "cellid += 1\n", "SHIM_bot_plug_cell = openmc.Cell(name='SHIM clad bot plug', cell_id=cellid)\n", "SHIM_bot_plug_cell.fill = CladdingCR\n", "SHIM_bot_plug_cell.region = -r_shim_bot_plug & +bot_shim_plug & -bot_shim_clad\n", "\n", "cellid += 1\n", "SHIM_top_clad_cell = openmc.Cell(name='SHIM clad top cell', cell_id=cellid)\n", "SHIM_top_clad_cell.fill = CladdingCR\n", "SHIM_top_clad_cell.region = -r_shim_outer & +top_shim & -top_clad_shim\n", "\n", "\n", "cellid += 1\n", "SHIM_water_top = openmc.Cell(name='SHIM water Cell top', cell_id=cellid)\n", "SHIM_water_top.fill = Water\n", "SHIM_water_top.region = -r_clad_103 & +top_clad_shim & -end_domain_top\n", "\n", "cellid += 1\n", "SHIM_water_bot = openmc.Cell(name='SHIM water Cell bot', cell_id=cellid)\n", "SHIM_water_bot.fill = Water\n", "SHIM_water_bot.region = -r_clad_103 & +end_domain_bot & -bot_shim_plug\n", "\n", "cellid += 1\n", "SHIM_water_bot_plug = openmc.Cell(name='SHIM water Cell bot plug', cell_id=cellid)\n", "SHIM_water_bot_plug.fill = Water\n", "SHIM_water_bot_plug.region = -r_clad_103 & +r_shim_bot_plug & +bot_shim_plug & -bot_shim_clad\n", "\n", "cellid += 1\n", "SHIM_water_fill = openmc.Cell(name='SHIM water Cell for filling lattice', cell_id=cellid)\n", "SHIM_water_fill.fill = Water\n", "SHIM_water_fill.region = +r_shim_outer & -r_clad_103 & +bot_shim_clad & -top_clad_shim\n", "\n", "shim_universe.add_cell(SHIM_cell)\n", "shim_universe.add_cell(SHIM_clad)\n", "shim_universe.add_cell(SHIM_water_top)\n", "shim_universe.add_cell(SHIM_water_bot)\n", "shim_universe.add_cell(SHIM_water_bot_plug)\n", "shim_universe.add_cell(SHIM_water_fill)\n", "shim_universe.add_cell(SHIM_top_clad_cell)\n", "shim_universe.add_cell(SHIM_bot_clad_cell)\n", "shim_universe.add_cell(SHIM_bot_plug_cell)\n", "\n", "\n", "if inline_plot:\n", " shim_universe.plot(origin=(0,0,z0_shim), width=(6.0, 3*z0_shim), pixels=[400, 400], basis='xz', color_by='material')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Dummy channel (graphite) \n", "graphite_channel_universe = openmc.Universe(name='Graphite Channel Universe')\n", "\n", "cellid += 1\n", "graphite_channel_cell = openmc.Cell(name='Graphite channel Cell', cell_id=cellid)\n", "graphite_channel_cell.fill = Graphite\n", "graphite_channel_cell.region = -r_dummy_graphite & +reflector_bot_101 & -reflector_top_101\n", "\n", "cellid += 1\n", "graphite_clad = openmc.Cell(name='Graphite Clad Cell', cell_id=cellid)\n", "graphite_clad.fill = CladdingCR\n", "graphite_clad.region = -r_dummy_clad & +r_dummy_graphite & +reflector_bot_101 & -reflector_top_101\n", "\n", "cellid += 1\n", "dummy_top_plug_cell = openmc.Cell(name='Graphite channel Cell plug top', cell_id=cellid)\n", "dummy_top_plug_cell.fill = CladdingCR\n", "dummy_top_plug_cell.region = -r_dummy_clad & +reflector_top_101 & -clad_top_101\n", "\n", "cellid += 1\n", "dummy_bot_plug_cell = openmc.Cell(name='Graphite channel Cell plug bot', cell_id=cellid)\n", "dummy_bot_plug_cell.fill = CladdingCR\n", "dummy_bot_plug_cell.region = -r_dummy_clad & +clad_bot_101 & -reflector_bot_101\n", "\n", "cellid += 1\n", "dummy_bot_water_cell = openmc.Cell(cell_id=cellid)\n", "dummy_bot_water_cell.fill = Water\n", "dummy_bot_water_cell.region = -r_dummy_clad & -clad_bot_101 & +end_domain_bot\n", "\n", "cellid += 1\n", "dummy_top_water_cell = openmc.Cell(cell_id=cellid)\n", "dummy_top_water_cell.fill = Water\n", "dummy_top_water_cell.region = -r_dummy_clad & +clad_top_101 & -end_domain_top\n", "\n", "\n", "graphite_channel_universe.add_cell(graphite_channel_cell)\n", "graphite_channel_universe.add_cell(graphite_clad)\n", "graphite_channel_universe.add_cell(dummy_top_plug_cell)\n", "graphite_channel_universe.add_cell(dummy_bot_plug_cell)\n", "graphite_channel_universe.add_cell(dummy_top_water_cell)\n", "graphite_channel_universe.add_cell(dummy_bot_water_cell)\n", "\n", "\n", "if inline_plot:\n", " graphite_channel_universe.plot(origin=(0,0,30), width=(7.0, 70), pixels=[400, 400], basis='xz', color_by='material')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Rabbit & Source empty channels\n", "rabbit_universe = openmc.Universe(name='Rabbit Universe')\n", "\n", "cellid += 1\n", "rabbit_cell = openmc.Cell(name='Rabbit', cell_id=cellid)\n", "rabbit_cell.fill = Vacuum\n", "rabbit_cell.region = -r_rabbit_inner & +bottom_grid_top & -end_domain_top\n", "\n", "cellid += 1\n", "rabbit_clad = openmc.Cell(name='Rabbit clad', cell_id=cellid)\n", "rabbit_clad.fill = CladdingCR\n", "rabbit_clad.region = +r_rabbit_inner & -r_rabbit_outer & +bottom_grid_top & -end_domain_top\n", "\n", "rabbit_universe.add_cell(rabbit_cell)\n", "rabbit_universe.add_cell(rabbit_clad)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# Reflector \n", "reflector_universe = openmc.Universe(name='Reflector')\n", "\n", "# Create reflector cell\n", "# --> 4 void zones inside the reflector: 1 cylindrical annulus and 3 channels\n", "\n", "# Cylindrical ring between 31.35 and 57\n", "cellid += 1\n", "reflector_cell_down = openmc.Cell(name='Reflector cell down', cell_id=cellid)\n", "reflector_cell_down.fill = Graphite\n", "reflector_cell_down.region = -outer_radius_reflector & +inner_radius_reflector & +bot_reflector & -refl_void_b & +channel_1_outer & +channel_2_outer & +channel_3_outer #aggiungere canali\n", "\n", "cellid += 1\n", "reflector_cell_fill_ch1 = openmc.Cell(name='Reflector cell fill ch1', cell_id=cellid)\n", "reflector_cell_fill_ch1.fill = Graphite\n", "reflector_cell_fill_ch1.region = -outer_radius_reflector & +inner_radius_reflector & +bot_reflector & -refl_void_b & -channel_1_outer & -channel_1_discriminator\n", "\n", "cellid += 1\n", "reflector_cell_fill_ch2 = openmc.Cell(name='Reflector cell fill ch2', cell_id=cellid)\n", "reflector_cell_fill_ch2.fill = Graphite\n", "reflector_cell_fill_ch2.region = -outer_radius_reflector & +inner_radius_reflector & +bot_reflector & -refl_void_b & -channel_2_outer & +channel_2_discriminator\n", "\n", "cellid += 1\n", "reflector_cell_fill_ch3 = openmc.Cell(name='Reflector cell fill ch3', cell_id=cellid)\n", "reflector_cell_fill_ch3.fill = Graphite\n", "reflector_cell_fill_ch3.region = -outer_radius_reflector & +inner_radius_reflector & +bot_reflector & -refl_void_b & -channel_3_outer & +channel_3_discriminator & +channel_2_outer\n", "\n", "cellid += 1\n", "reflector_cell_right_vacuum = openmc.Cell(name='Reflector right side from vacuum', cell_id=cellid)\n", "reflector_cell_right_vacuum.fill = Graphite\n", "reflector_cell_right_vacuum.region = -refl_void_inner & +inner_radius_reflector & +refl_void_b & -top_reflector\n", "\n", "cellid += 1\n", "reflector_cell_left_vacuum = openmc.Cell(name='Reflector left side from vacuum', cell_id=cellid)\n", "reflector_cell_left_vacuum.fill = Graphite\n", "reflector_cell_left_vacuum.region = +refl_void_outer & -outer_radius_reflector & +refl_void_b & -top_reflector \n", "\n", "# lazy susan\n", "cellid += 1\n", "reflector_cell_vacuum = openmc.Cell(name='Reflector vacuum', cell_id=cellid)\n", "reflector_cell_vacuum.fill = Vacuum\n", "reflector_cell_vacuum.region = -refl_void_outer & +refl_void_inner & +refl_void_b & -top_reflector #aggiungere canali\n", "\n", "cellid += 1\n", "reflector_clad_bot = openmc.Cell(name='Reflector clad bot', cell_id=cellid)\n", "reflector_clad_bot.fill = CladdingCR\n", "reflector_clad_bot.region = +core_radius & -s1 & +bot_clad_reflector & -bot_reflector\n", "\n", "cellid += 1\n", "reflector_clad_top = openmc.Cell(name='Reflector clad top', cell_id=cellid)\n", "reflector_clad_top.fill = CladdingCR\n", "reflector_clad_top.region = +core_radius & -s1 & +top_reflector & -top_clad_reflector\n", "\n", "cellid += 1\n", "reflector_clad_side_inner = openmc.Cell(name='Reflector clad side inner', cell_id=cellid)\n", "reflector_clad_side_inner.fill = CladdingCR\n", "reflector_clad_side_inner.region = +core_radius & -inner_radius_reflector & +bot_reflector & -top_reflector\n", "\n", "cellid += 1\n", "reflector_clad_side_outer = openmc.Cell(name='Reflector clad side outer', cell_id=cellid)\n", "reflector_clad_side_outer.fill = CladdingCR\n", "reflector_clad_side_outer.region = +outer_radius_reflector & -s1 & +bot_reflector & -top_reflector\n", "\n", "\n", "#CHANNELS\n", "cellid += 1\n", "reflector_ch1 = openmc.Cell(name='irradiation channel 1', cell_id=cellid)\n", "reflector_ch1.fill = Vacuum\n", "reflector_ch1.region = -channel_1_inner & +inner_radius_reflector & -outer_radius_reflector & +channel_1_discriminator\n", "\n", "cellid += 1\n", "reflector_ch1_clad = openmc.Cell(name='clad irradiation channel 1', cell_id=cellid)\n", "reflector_ch1_clad.fill = CladdingCR\n", "reflector_ch1_clad.region = +channel_1_inner & -channel_1_outer & +inner_radius_reflector & -outer_radius_reflector & +channel_1_discriminator\n", "\n", "cellid += 1\n", "reflector_ch2 = openmc.Cell(name='irradiation channel 2', cell_id=cellid)\n", "reflector_ch2.fill = Vacuum\n", "reflector_ch2.region = -channel_2_inner & +inner_radius_reflector & -outer_radius_reflector & -channel_2_discriminator\n", "\n", "cellid += 1\n", "reflector_ch2_clad = openmc.Cell(name='clad irradiation channel 2', cell_id=cellid)\n", "reflector_ch2_clad.fill = CladdingCR\n", "reflector_ch2_clad.region = +channel_2_inner & -channel_2_outer & +inner_radius_reflector & -outer_radius_reflector & -channel_2_discriminator\n", "\n", "cellid += 1\n", "reflector_ch3 = openmc.Cell(name='irradiation channel 3', cell_id=cellid)\n", "reflector_ch3.fill = Vacuum\n", "reflector_ch3.region = -channel_3_inner & +channel_2_outer & -outer_radius_reflector & -channel_3_discriminator\n", "\n", "cellid += 1\n", "reflector_ch3_clad = openmc.Cell(name='clad irradiation channel 3', cell_id=cellid)\n", "reflector_ch3_clad.fill = CladdingCR\n", "reflector_ch3_clad.region = +channel_3_inner & -channel_3_outer & +channel_2_outer & -outer_radius_reflector & -channel_3_discriminator\n", "\n", "cellid += 1\n", "water_reflector_top= openmc.Cell(cell_id=cellid)\n", "water_reflector_top.region= +top_clad_reflector & -end_domain_top & -s1 & +core_radius\n", "water_reflector_top.fill = Water\n", "\n", "cellid += 1\n", "water_reflector_bottom= openmc.Cell(cell_id=cellid)\n", "water_reflector_bottom.region= +end_domain_bot & -bottom_grid_bot & -s1\n", "water_reflector_bottom.fill = Water\n", "\n", "cellid += 1\n", "water_reflector_bottom2= openmc.Cell(cell_id=cellid)\n", "water_reflector_bottom2.region= -bot_clad_reflector & +bottom_grid_bot & -s1 & +core_radius\n", "water_reflector_bottom2.fill = Water\n", "\n", "\n", "reflector_universe.add_cell(reflector_cell_down)\n", "reflector_universe.add_cell(reflector_cell_vacuum)\n", "reflector_universe.add_cell(reflector_cell_right_vacuum)\n", "reflector_universe.add_cell(reflector_cell_left_vacuum)\n", "reflector_universe.add_cell(reflector_clad_bot)\n", "reflector_universe.add_cell(reflector_clad_top)\n", "reflector_universe.add_cell(reflector_clad_side_inner)\n", "reflector_universe.add_cell(reflector_clad_side_outer)\n", "reflector_universe.add_cell(reflector_ch1)\n", "reflector_universe.add_cell(reflector_cell_fill_ch1)\n", "reflector_universe.add_cell(reflector_ch1_clad)\n", "reflector_universe.add_cell(reflector_ch2)\n", "reflector_universe.add_cell(reflector_cell_fill_ch2)\n", "reflector_universe.add_cell(reflector_ch2_clad)\n", "reflector_universe.add_cell(reflector_ch3)\n", "reflector_universe.add_cell(reflector_cell_fill_ch3)\n", "reflector_universe.add_cell(reflector_ch3_clad)\n", "reflector_universe.add_cell(water_reflector_top)\n", "reflector_universe.add_cell(water_reflector_bottom)\n", "reflector_universe.add_cell(water_reflector_bottom2)\n", "\n", "\n", "if inline_plot:\n", " reflector_universe.plot(origin=(0,0,25), width=(130.0, 130), pixels=[500, 500], basis='xz', color_by='material')\n", " reflector_universe.plot(origin=(0,0,25), width=(130.0, 130), pixels=[500, 500], basis='xy', color_by='material')" ] }, { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [], "source": [ "# Concrete --> to do" ] }, { "cell_type": "code", "execution_count": 39, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Adding in a central channel rod...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a SHIM rod...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a transient rod...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a REG rod...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a graphite rod...\n", "Adding in a graphite rod...\n", "Adding in a water rod...\n", "Adding in a graphite rod...\n", "Adding in a graphite rod...\n", "Adding in a graphite rod...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a water rod...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n", "Adding in a 103 Fuel type...\n", "Adding in a 101 Fuel type...\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=91.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=92.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=93.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=94.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=95.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=96.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=97.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=98.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=99.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=100.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=101.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=102.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=103.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=104.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=105.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=106.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=107.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=108.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=109.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=300.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=110.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=400.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=111.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=401.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=112.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=402.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=113.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=403.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=114.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=404.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=115.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=405.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=116.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=500.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=117.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=501.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=118.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=502.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=119.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=503.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=120.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=504.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=121.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=505.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=122.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=506.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=123.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=507.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=124.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=508.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=125.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=509.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=126.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=510.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=127.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=511.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=128.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=600.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=129.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=601.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=130.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=602.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=131.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=603.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=132.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=604.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=133.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=605.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=134.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=606.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=135.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=607.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=136.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=608.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=137.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=609.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=138.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=610.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=139.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=611.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=140.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=612.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=141.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=613.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=142.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=614.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=143.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=615.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=144.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=616.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=145.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=617.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=146.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=700.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=147.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=701.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=148.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=702.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=149.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=703.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=150.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=704.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=151.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=705.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=152.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=706.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=153.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=707.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=154.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=708.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=155.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=709.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=156.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=710.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=157.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=711.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=158.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=712.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=159.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=713.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=160.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=714.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=161.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=715.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=162.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=716.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=163.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=717.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=164.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=718.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=165.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=719.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=166.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=720.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=167.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=721.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=168.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=722.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=169.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=723.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=170.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=800.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=171.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=801.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=172.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=802.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=173.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=803.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=174.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=804.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=175.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=805.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=176.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=806.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=177.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=807.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=178.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=808.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=179.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=809.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=180.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=810.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=181.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=811.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=182.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=812.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=183.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=813.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=184.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=814.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=185.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=815.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=186.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=816.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=187.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=817.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=188.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=818.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=189.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=819.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=190.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=820.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=191.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=821.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=192.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=822.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=193.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=823.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=194.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=824.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=195.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=825.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=196.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=826.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=197.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=827.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=198.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=828.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=199.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=829.\n", " warn(msg, IDWarning)\n" ] } ], "source": [ "# Core universe\n", "\n", "# Create surfaces that will divide rings in the circular lattice\n", "ring_radii = np.array([0.0, 4.0975, 8.2, 12.2, 16.3, 20.3])# From Riccardo Boccelli MSc thesis\n", "radial_surf = [openmc.ZCylinder(r=r) for r in\n", " (ring_radii[:-1] + ring_radii[1:])/2]\n", "\n", "water_cells_down = []\n", "water_cells_up = []\n", "top_grid_cells = []\n", "for i in range(ring_radii.size):\n", " # Create annular region\n", " if i == 0:\n", " water_region_down = -radial_surf[i] & -top_grid_bot \n", " water_region_up = -radial_surf[i] & +top_grid_top\n", " grid_region = -radial_surf[i] & +top_grid_bot & -top_grid_top\n", " elif i == ring_radii.size - 1:\n", " water_region_down = -top_grid_bot & +radial_surf[i-1] & -core_radius \n", " water_region_up = +top_grid_top & +radial_surf[i-1] & -core_radius\n", " grid_region = +radial_surf[i-1] & +top_grid_bot & -top_grid_top\n", " else:\n", " water_region_down = -top_grid_bot & +radial_surf[i-1] & -radial_surf[i] \n", " water_region_up = +top_grid_top & +radial_surf[i-1] & -radial_surf[i] \n", " grid_region = +radial_surf[i-1] & -radial_surf[i] & +top_grid_bot & -top_grid_top\n", " \n", " cellid +=1 \n", " water_cell_down = openmc.Cell(fill=Water, region=water_region_down, cell_id=cellid)\n", " \n", " cellid +=1 \n", " water_cell_up = openmc.Cell(fill=Water, region=water_region_up, cell_id=cellid)\n", " \n", " cellid +=1 \n", " top_grid_cell = openmc.Cell(fill=CladdingAl, region=grid_region, cell_id=cellid)\n", " top_grid_cells.append(top_grid_cell) \n", " \n", " water_cells_down.append(water_cell_down)\n", " water_cells_up.append(water_cell_up)\n", " \n", "\n", "core_universe = openmc.Universe()\n", "core_universe.add_cells(water_cells_down)\n", "core_universe.add_cells(water_cells_up)\n", "core_universe.add_cells(top_grid_cells)\n", "\n", "\n", "# Arrange the pins in the circular lattice \n", "num_pins = [1, 6, 12, 18, 24, 30]\n", "angles = [0, 0, 0, 0, 0, 0]\n", "#angles = [0, 30, 0, 30, 0, 30]\n", "\n", "# In order to insert pin in the lattice, it is needed to define first the parameters for the spatial location.\n", "# In particular: 'numPins' identifies the ring, and 'howleft' the angular degree\n", "\n", "centralChannel = {'numPins':[num_pins[0]], # Here the 'howLeft' parameter it is not necessary due to only one pin\n", " 'howLeftFrom3oclock':[0]}\n", "\n", "REG = {'numPins' :[num_pins[4]],\n", " 'howLeftFrom3oclock':[7]}\n", "\n", "TRANS = {'numPins' : [num_pins[3]],\n", " 'howLeftFrom3oclock':[11]}\n", "\n", "SHIM = {'numPins' : [num_pins[2]],\n", " 'howLeftFrom3oclock':[-1]}\n", "\n", "waterRods = {'numPins' :[num_pins[5], num_pins[5]],\n", " 'howLeftFrom3oclock': [1,11]}\n", "\n", "graphiteRods = {'numPins' : np.full(5, num_pins[5]),\n", " 'howLeftFrom3oclock':[-1,0,2,3,4]}\n", "\n", "# 6 pin in ring 1, 11 in ring 2 and so on\n", "Fuel_1033 = {'numPins': [num_pins[1],num_pins[1],num_pins[1],num_pins[1],num_pins[1],num_pins[1],\n", " \n", " num_pins[2], num_pins[2], num_pins[2],num_pins[2], num_pins[2], num_pins[2], num_pins[2],num_pins[2],\n", " num_pins[2],num_pins[2],num_pins[2],\n", "\n", " num_pins[3], \n", "\n", " num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5],\n", " num_pins[5], num_pins[5], num_pins[5], num_pins[5] \n", " ],\n", " 'howLeftFrom3oclock': [-1,0,1,2,3,4, 0,1,2,3,4,5,6,7,8,9,10, 14, 5,7,9,13,15,17,19,21,23,24,25,27 ] }\n", "\n", "\n", "#corona 3, ho (18-2) pins\n", "Fuel_1011 = {'numPins': [num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], num_pins[3], \n", "\n", " num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], num_pins[4], \n", "\n", " num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5], num_pins[5]\n", "\n", " ],\n", " 'howLeftFrom3oclock': [-1,0,1,2,3,4,5,6,7,8,9,10,12,13,15,16, -1,0,1,2,3,4,5,6,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22, 6,8,10,12,14,16,18,20,22,26,28] }\n", "\n", "\n", "\n", "# Here I define some variables in order to locate the pins in the next for iteration\n", "\n", "def CentralChannel(centralChannel,n,j):\n", " for irod in range(len(centralChannel['numPins'])): \n", " if n == centralChannel['numPins'][irod]:\n", " return True\n", " return False\n", "\n", "def REGChannel(REG,n,j):\n", " for irod in range(len(REG['numPins'])): \n", " if n == REG['numPins'][irod] and \\\n", " j-1 == REG['howLeftFrom3oclock'][irod]:\n", " return True\n", " return False\n", "\n", "def SHIMChannel(SHIM,n,j):\n", " for irod in range(len(SHIM['numPins'])): \n", " if n == SHIM['numPins'][irod] and \\\n", " j-1 == SHIM['howLeftFrom3oclock'][irod]:\n", " return True\n", " return False\n", "\n", "def TransientRod(TRANS,n,j):\n", " for irod in range(len(TRANS['numPins'])): \n", " if n == TRANS['numPins'][irod] and \\\n", " j-1 == TRANS['howLeftFrom3oclock'][irod]:\n", " return True\n", " return False\n", "\n", "def WaterRod(waterRods,n,j):\n", " for irod in range(len(waterRods['numPins'])): \n", " if n == waterRods['numPins'][irod] and \\\n", " j-1 == waterRods['howLeftFrom3oclock'][irod]:\n", " return True\n", " return False\n", "\n", "def GraphiteRod(graphiteRods,n,j):\n", " for irod in range(len(graphiteRods['numPins'])): \n", " if n == graphiteRods['numPins'][irod] and \\\n", " j-1 == graphiteRods['howLeftFrom3oclock'][irod]:\n", " return True\n", " return False\n", "\n", "def Fuel_103(Fuel_1033,n,j):\n", " for irod in range(len(Fuel_1033['numPins'])): \n", " if n == Fuel_1033['numPins'][irod] and \\\n", " j-1 == Fuel_1033['howLeftFrom3oclock'][irod]:\n", " return True\n", " return False\n", "\n", "\n", "def Fuel_101(Fuel_1011,n,j):\n", " for irod in range(len(Fuel_1011['numPins'])): \n", " if n == Fuel_1011['numPins'][irod] and \\\n", " j-1 == Fuel_1011['howLeftFrom3oclock'][irod]:\n", " return True\n", " return False\n", "\n", "\n", "\n", "# Filling the core\n", "for i, (r, n, a) in enumerate(zip(ring_radii, num_pins, angles)):\n", "\n", " for j in range(n):\n", " \n", " # Determine location of center of pin\n", " theta = (a + j/n*360.) * np.pi/180.\n", " x = r*np.cos(theta)\n", " y = r*np.sin(theta)\n", " cellid += 1\n", " \n", " pin_boundary = openmc.ZCylinder(x0=x, y0=y, r=r_clad_101.r)\n", " pin_boundary_CC = openmc.ZCylinder(x0=x, y0=y, r=r_central_channel_clad.r)\n", " water_cells_down[i].region &= +pin_boundary & -core_radius\n", " water_cells_up[i].region &= +pin_boundary & -core_radius\n", " top_grid_cells[i].region &= +pin_boundary & -core_radius\n", "\n", " # Create each fuel pin -- note that we explicitly assign an ID so \n", " # that we can identify the pin later when looking at tallies\n", " if CentralChannel(centralChannel,n,j):\n", " print('Adding in a central channel rod...')\n", " pin = openmc.Cell(fill=central_channel_universe, region=-pin_boundary_CC, cell_id=cellid)\n", " elif TransientRod(TRANS,n,j):\n", " print('Adding in a transient rod...')\n", " pin = openmc.Cell(fill=transient_universe, region=-pin_boundary, cell_id=cellid)\n", " elif WaterRod(waterRods,n,j):\n", " print('Adding in a water rod...')\n", " pin = openmc.Cell(fill=rabbit_universe, region=-pin_boundary, cell_id=cellid)\n", " elif SHIMChannel(SHIM,n,j):\n", " print('Adding in a SHIM rod...')\n", " pin = openmc.Cell(fill=shim_universe, region=-pin_boundary, cell_id=cellid)\n", " elif REGChannel(REG,n,j):\n", " print('Adding in a REG rod...')\n", " pin = openmc.Cell(fill=reg_universe, region=-pin_boundary, cell_id=cellid)\n", " elif GraphiteRod(graphiteRods,n,j):\n", " print('Adding in a graphite rod...')\n", " pin = openmc.Cell(fill=graphite_channel_universe, region=-pin_boundary, cell_id=cellid)\n", " elif Fuel_103(Fuel_1033,n,j):\n", " print('Adding in a 103 Fuel type...')\n", " pin = openmc.Cell(fill=fuel_103_universe, region=-pin_boundary, cell_id=cellid)\n", "\n", " elif Fuel_101(Fuel_1011,n,j):\n", " print('Adding in a 101 Fuel type...')\n", " pin = openmc.Cell(fill=fuel_101_universe, region=-pin_boundary, cell_id=cellid)\n", " \n", " pin.translation = (x, y, 0)\n", " pin.id = (i + 3)*100 + j \n", " core_universe.add_cell(pin)\n", "\n", "\n", "if inline_plot:\n", " core_universe.plot(origin=(0,0,25), width=(60.0, 100), pixels=[400, 400], basis='xz', color_by='material', seed = 3)\n", "\n", "if inline_plot:\n", " core_universe.plot(origin=(0,0,25), width=(60.0, 60), pixels=[600, 600], basis='xy', color_by='material', seed = 2)" ] }, { "cell_type": "code", "execution_count": 40, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=200.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=201.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=202.\n", " warn(msg, IDWarning)\n", "/home/renato/miniconda3/envs/mc/lib/python3.11/site-packages/openmc/mixin.py:70: IDWarning: Another Cell instance already exists with id=203.\n", " warn(msg, IDWarning)\n" ] } ], "source": [ "# Geometry definitions for the reactor\n", "\n", "cellid += 1\n", "reactor = openmc.Cell(cell_id = cellid)\n", "reactor.region = -core_radius & -end_domain_top & +bottom_grid_top\n", "reactor.fill = core_universe\n", "\n", "cellid += 1\n", "reflector = openmc.Cell(cell_id = cellid)\n", "reflector.region = +core_radius & -s1 & +end_domain_bot & -end_domain_top\n", "reflector.fill = reflector_universe\n", "\n", "cellid += 1\n", "bottom_grid_cell = openmc.Cell(name='bottom grid cell', cell_id = cellid)\n", "bottom_grid_cell.fill = CladdingCR\n", "bottom_grid_cell.region = -core_radius & +bottom_grid_bot & -bottom_grid_top\n", "\n", "cellid += 1\n", "bottom_grid_cell_water = openmc.Cell(name='bottom grid cell', cell_id = cellid)\n", "bottom_grid_cell_water.fill = Water\n", "bottom_grid_cell_water.region = -core_radius & -bottom_grid_bot & +end_domain_bot\n", "\n", "\n", "### --- ROOT UNIVERSE \n", "triga_universe = openmc.Universe(name='TRIGA universe')\n", "\n", "triga_universe.add_cell(reactor)\n", "triga_universe.add_cell(reflector)\n", "triga_universe.add_cell(bottom_grid_cell)\n", "triga_universe.add_cell(bottom_grid_cell_water)\n", "\n", "\n", "geometry = openmc.Geometry((triga_universe))\n", "\n", "if inline_plot:\n", " triga_universe.plot(origin=(0,0,25), width=(130.0, 130), pixels=[700, 700], basis='xz', color_by='material', seed = 3)\n", "\n", "if inline_plot:\n", " triga_universe.plot(origin=(0,0,25), width=(130.0, 130), pixels=[700, 700], basis='xy', color_by='material', seed = 3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plots" ] }, { "cell_type": "code", "execution_count": 41, "metadata": {}, "outputs": [], "source": [ "# Plot definition\n", "\n", "folder_path = './pictures/geometry/'\n", "if not os.path.exists(folder_path):\n", " os.makedirs(folder_path)\n", "\n", "colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76), }\n", "\n", "plot_xy1 = openmc.Plot()\n", "plot_xy1.basis = 'xy'\n", "plot_xy1.origin = (0.0, 0.0, 23.15-7) #23.15 , 28.15, 10.4 (start CR)\n", "plot_xy1.width = (110, 110)\n", "plot_xy1.pixels = (2000, 2000)\n", "plot_xy1.color_by = 'material'\n", "plot_xy1.filename = folder_path+'view_xy1'\n", "plot_xy1.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "plot_xy2 = openmc.Plot()\n", "plot_xy2.basis = 'xy'\n", "plot_xy2.origin = (0.0, 0.0, 23.15-4) #23.15 , 28.15, 10.4 (start CR)\n", "plot_xy2.width = (110, 110)\n", "plot_xy2.pixels = (2000, 2000)\n", "plot_xy2.color_by = 'material'\n", "plot_xy2.filename = folder_path+'view_xy2'\n", "plot_xy2.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "plot_xy3 = openmc.Plot()\n", "plot_xy3.basis = 'xy'\n", "plot_xy3.origin = (0.0, 0.0, 23.15) #23.15 , 28.15, 10.4 (start CR)\n", "plot_xy3.width = (110, 110)\n", "plot_xy3.pixels = (2000, 2000)\n", "plot_xy3.color_by = 'material'\n", "plot_xy3.filename = folder_path+'view_xy3'\n", "plot_xy3.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "plot_xy4 = openmc.Plot()\n", "plot_xy4.basis = 'xy'\n", "plot_xy4.origin = (0.0, 0.0, 23.15 + 2) #23.15 , 28.15, 10.4 (start CR)\n", "plot_xy4.width = (110, 110)\n", "plot_xy4.pixels = (2000, 2000)\n", "plot_xy4.color_by = 'material'\n", "plot_xy4.filename = folder_path+'view_xy4'\n", "plot_xy4.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "plot_xy5 = openmc.Plot()\n", "plot_xy5.basis = 'xy'\n", "plot_xy5.origin = (0.0, 0.0, 23.15+4) #23.15 , 28.15, 10.4 (start CR)\n", "plot_xy5.width = (110, 110)\n", "plot_xy5.pixels = (2000, 2000)\n", "plot_xy5.color_by = 'material'\n", "plot_xy5.filename = folder_path+'view_xy5'\n", "plot_xy5.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "plot_xy6 = openmc.Plot()\n", "plot_xy6.basis = 'xy'\n", "plot_xy6.origin = (0.0, 0.0, 23.15 + 7) #23.15 , 28.15, 10.4 (start CR)\n", "plot_xy6.width = (110, 110)\n", "plot_xy6.pixels = (2000, 2000)\n", "plot_xy6.color_by = 'material'\n", "plot_xy6.filename = folder_path+'view_xy6'\n", "plot_xy6.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "plot_vista_shim = openmc.Plot()\n", "plot_vista_shim.basis = 'xz'\n", "plot_vista_shim.origin = (0.0, 0, 0)\n", "plot_vista_shim.width = (110, 230)\n", "plot_vista_shim.pixels = (2000, 2000)\n", "plot_vista_shim.color_by = 'material'\n", "plot_vista_shim.filename = folder_path+'view_shim'\n", "plot_vista_shim.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "\n", "plot_vista_reg = openmc.Plot()\n", "plot_vista_reg.basis = 'xz'\n", "plot_vista_reg.origin = (0, 14.9, 0)\n", "plot_vista_reg.width = (110, 150)\n", "plot_vista_reg.pixels = (3000, 3000)\n", "plot_vista_reg.color_by = 'material'\n", "plot_vista_reg.filename = folder_path+'view_reg'\n", "plot_vista_reg.colors = {\n", " Water: (0,102,204),\n", " CladdingSS: (190,190,190), \n", " CladdingCR: (196,196,196),\n", " CladdingAl: (200,200,200),\n", " Fuel101: (255,30,30),\n", " Fuel103: (155,10,10),\n", " Graphite: (110,110,110),\n", " Vacuum: (249,255,227),\n", " ControlRodB4C: (152,251,152),\n", " ControlRodBorated: (152,251,152),\n", " Samarium: (255,255,51),\n", " Zirconium: (0,153,76),\n", "\n", "}\n", "\n", "\n", "plot_vista_trans = openmc.Plot()\n", "plot_vista_trans.basis = 'xz'\n", "plot_vista_trans.origin = (0, -11.5, 0)\n", "plot_vista_trans.width = (110, 150)\n", "plot_vista_trans.pixels = (2000, 2000)\n", "plot_vista_trans.color_by = 'material'\n", "plot_vista_trans.filename = folder_path+'view_trans'\n", "plot_vista_trans.colors = colors\n", "\n", "\n", "\n", "plots = openmc.Plots()\n", "plots.append(plot_xy1)\n", "plots.append(plot_xy2)\n", "plots.append(plot_xy3)\n", "plots.append(plot_xy4)\n", "plots.append(plot_xy5)\n", "plots.append(plot_xy6)\n", "plots.append(plot_vista_reg)\n", "plots.append(plot_vista_shim)\n", "plots.append(plot_vista_trans)\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Settings" ] }, { "cell_type": "code", "execution_count": 42, "metadata": {}, "outputs": [], "source": [ "# OpenMC simulation parameters \n", "\n", "entropy_calc = False\n", "resonance_calc = False\n", "uniform_source = True\n", "\n", "batches = 300\n", "inactive = 100\n", "particles = 120000\n", "\n", "settings_file = openmc.Settings()\n", "settings_file.batches = batches\n", "settings_file.inactive = inactive\n", "settings_file.particles = particles\n", "\n", "# uniform initial source\n", "if uniform_source:\n", " bounds = [-25, -25, 5, 25, 25, 50]\n", " uniform_dist = openmc.stats.Box(bounds[:3], bounds[3:], only_fissionable=True)\n", " settings_file.source = openmc.Source(space=uniform_dist)\n", "\n", "if resonance_calc:\n", " settings_file.resonance_scattering['enable'] = True\n", "\n", "if entropy_calc:\n", " entropy_mesh = openmc.Mesh()\n", " entropy_mesh.lower_left = (-60, -60, -5)\n", " entropy_mesh.upper_right = (60, 60, 60)\n", "\n", " entropy_mesh.dimension = (10, 10, 10)\n", " settings_file.entropy_mesh = entropy_mesh\n", "\n", "settings_file.temperature = {'method':'interpolation'}\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Tallies" ] }, { "cell_type": "code", "execution_count": 43, "metadata": {}, "outputs": [], "source": [ "# Measuring quantities during simulation\n", "tallies_file = openmc.Tallies()\n", "\n", "### Neutron energy flux \n", "energies = np.logspace(np.log10(1e-5), np.log10(20.0e6), 101)\n", "e_filter = openmc.EnergyFilter(energies)\n", "\n", "EnergySpectrum = openmc.Tally(name = \"Phi(E)\")\n", "EnergySpectrum.filters = [e_filter]\n", "EnergySpectrum.scores = ['flux']\n", "\n", "\n", "### Spatial flux\n", "mesh3D = openmc.RegularMesh()\n", "mesh3D.dimension = (60, 60, 100)\n", "mesh3D.lower_left = (-55, -55, -10)\n", "mesh3D.upper_right = (55, 55, 100)\n", "\n", "mesh3D_filter = openmc.MeshFilter(mesh3D)\n", "\n", "SpatialFlux = openmc.Tally(name = \"Phi(r)\")\n", "SpatialFlux.filters = [mesh3D_filter]\n", "SpatialFlux.scores = ['flux']\n", "\n", "\n", "### Deposited energy (for normalization)\n", "heating = openmc.Tally(name='factor-for-normalization')\n", "heating.scores = ['heating-local']\n", "\n", "tallies_file.append(EnergySpectrum)\n", "tallies_file.append(SpatialFlux)\n", "tallies_file.append(heating)\n" ] }, { "cell_type": "code", "execution_count": 44, "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.3\n", " Git SHA1 | 50e39a4e20dc9e0f3d7ccf07333f6a5e6c797c8c\n", " Date/Time | 2023-12-11 09:42:19\n", " OpenMP Threads | 16\n", "\n", " Reading settings XML file...\n", " Reading cross sections XML file...\n", " Reading materials XML file...\n", " Reading geometry XML file...\n", " Reading tallies XML file...\n", " Preparing distributed cell instances...\n", " Reading plot XML file...\n", "\n", " =======================> PLOTTING SUMMARY <========================\n", "\n", "Plot ID: 10\n", "Plot file: ./pictures/geometry/view_xy1.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 0 16.15\n", "Width: 110 110\n", "Coloring: Materials\n", "Basis: XY\n", "Pixels: 2000 2000\n", "\n", "Plot ID: 11\n", "Plot file: ./pictures/geometry/view_xy2.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 0 19.15\n", "Width: 110 110\n", "Coloring: Materials\n", "Basis: XY\n", "Pixels: 2000 2000\n", "\n", "Plot ID: 12\n", "Plot file: ./pictures/geometry/view_xy3.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 0 23.15\n", "Width: 110 110\n", "Coloring: Materials\n", "Basis: XY\n", "Pixels: 2000 2000\n", "\n", "Plot ID: 13\n", "Plot file: ./pictures/geometry/view_xy4.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 0 25.15\n", "Width: 110 110\n", "Coloring: Materials\n", "Basis: XY\n", "Pixels: 2000 2000\n", "\n", "Plot ID: 14\n", "Plot file: ./pictures/geometry/view_xy5.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 0 27.15\n", "Width: 110 110\n", "Coloring: Materials\n", "Basis: XY\n", "Pixels: 2000 2000\n", "\n", "Plot ID: 15\n", "Plot file: ./pictures/geometry/view_xy6.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 0 30.15\n", "Width: 110 110\n", "Coloring: Materials\n", "Basis: XY\n", "Pixels: 2000 2000\n", "\n", "Plot ID: 17\n", "Plot file: ./pictures/geometry/view_reg.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 14.9 0\n", "Width: 110 150\n", "Coloring: Materials\n", "Basis: XZ\n", "Pixels: 3000 3000\n", "\n", "Plot ID: 16\n", "Plot file: ./pictures/geometry/view_shim.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 0 0\n", "Width: 110 230\n", "Coloring: Materials\n", "Basis: XZ\n", "Pixels: 2000 2000\n", "\n", "Plot ID: 18\n", "Plot file: ./pictures/geometry/view_trans.png\n", "Universe depth: -1\n", "Plot Type: Slice\n", "Origin: 0 -11.5 0\n", "Width: 110 150\n", "Coloring: Materials\n", "Basis: XZ\n", "Pixels: 2000 2000\n", "\n", " Processing plot 10: ./pictures/geometry/view_xy1.png...\n", " Processing plot 11: ./pictures/geometry/view_xy2.png...\n", " Processing plot 12: ./pictures/geometry/view_xy3.png...\n", " Processing plot 13: ./pictures/geometry/view_xy4.png...\n", " Processing plot 14: ./pictures/geometry/view_xy5.png...\n", " Processing plot 15: ./pictures/geometry/view_xy6.png...\n", " Processing plot 17: ./pictures/geometry/view_reg.png...\n", " Processing plot 16: ./pictures/geometry/view_shim.png...\n", " Processing plot 18: ./pictures/geometry/view_trans.png...\n" ] } ], "source": [ "# Export to XML\n", "materials_file.export_to_xml()\n", "geometry.export_to_xml()\n", "plots.export_to_xml()\n", "settings_file.export_to_xml()\n", "tallies_file.export_to_xml()\n", "\n", "openmc.plot_geometry()" ] }, { "cell_type": "code", "execution_count": 45, "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.3\n", " Git SHA1 | 50e39a4e20dc9e0f3d7ccf07333f6a5e6c797c8c\n", " Date/Time | 2023-12-11 09:42:23\n", " OpenMP Threads | 16\n", "\n", " Reading settings XML file...\n", " Reading cross sections XML file...\n", " Reading materials XML file...\n", " Reading geometry XML file...\n", " Reading H1 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/H1.h5\n", " Reading O16 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/O16.h5\n", " Reading U235 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/U235.h5\n", " Reading U238 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/U238.h5\n", " Reading Zr90 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Zr90.h5\n", " Reading Zr91 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Zr91.h5\n", " Reading Zr92 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Zr92.h5\n", " Reading Zr94 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Zr94.h5\n", " Reading Zr96 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Zr96.h5\n", " Reading C0 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/C0.h5\n", " Reading Fe54 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Fe54.h5\n", " Reading Fe56 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Fe56.h5\n", " Reading Fe57 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Fe57.h5\n", " Reading Fe58 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Fe58.h5\n", " Reading Ti46 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ti46.h5\n", " Reading Ti47 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ti47.h5\n", " Reading Ti48 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ti48.h5\n", " Reading Ti49 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ti49.h5\n", " Reading Ti50 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ti50.h5\n", " Reading Na23 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Na23.h5\n", " Reading Mg24 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Mg24.h5\n", " Reading Mg25 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Mg25.h5\n", " Reading Mg26 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Mg26.h5\n", " Reading Cu63 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cu63.h5\n", " Reading Cu65 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cu65.h5\n", " Reading Co59 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Co59.h5\n", " Reading F19 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/F19.h5\n", " Reading Cl35 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cl35.h5\n", " Reading Cl37 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cl37.h5\n", " Reading B11 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/B11.h5\n", " Reading Al27 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Al27.h5\n", " Reading Cr50 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cr50.h5\n", " Reading Cr52 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cr52.h5\n", " Reading Cr53 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cr53.h5\n", " Reading Cr54 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Cr54.h5\n", " Reading Ni58 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ni58.h5\n", " Reading Ni60 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ni60.h5\n", " Reading Ni61 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ni61.h5\n", " Reading Ni62 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ni62.h5\n", " Reading Ni64 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Ni64.h5\n", " Reading Si28 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Si28.h5\n", " Reading Si29 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Si29.h5\n", " Reading Si30 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Si30.h5\n", " Reading Sm147 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Sm147.h5\n", " Reading Sm148 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Sm148.h5\n", " Reading Sm149 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Sm149.h5\n", " Reading Sm150 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Sm150.h5\n", " Reading Sm152 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Sm152.h5\n", " Reading Sm154 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/Sm154.h5\n", " Reading N14 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/N14.h5\n", " Reading B10 from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/B10.h5\n", " Reading c_H_in_H2O from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/c_H_in_H2O.h5\n", " Reading c_H_in_ZrH from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/c_H_in_ZrH.h5\n", " Reading c_Graphite from /home/renato/Desktop/xsdata/jeff-3.3-hdf5/c_Graphite.h5\n", " Minimum neutron data temperature: 294 K\n", " Maximum neutron data temperature: 600 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 H1\n", " Initializing source particles...\n", "\n", " ====================> K EIGENVALUE SIMULATION <====================\n", "\n", " Bat./Gen. k Average k\n", " ========= ======== ====================\n", " 1/1 0.96777\n", " 2/1 0.99735\n", " 3/1 1.01779\n", " 4/1 1.02098\n", " 5/1 1.02791\n", " 6/1 1.02535\n", " 7/1 1.02188\n", " 8/1 1.02741\n", " 9/1 1.02825\n", " 10/1 1.02687\n", " 11/1 1.02747\n", " 12/1 1.02886\n", " 13/1 1.02378\n", " 14/1 1.02980\n", " 15/1 1.02129\n", " 16/1 1.02872\n", " 17/1 1.01724\n", " 18/1 1.02503\n", " 19/1 1.02808\n", " 20/1 1.02649\n", " 21/1 1.02790\n", " 22/1 1.02462\n", " 23/1 1.03007\n", " 24/1 1.02630\n", " 25/1 1.03602\n", " 26/1 1.03232\n", " 27/1 1.03419\n", " 28/1 1.02863\n", " 29/1 1.03608\n", " 30/1 1.03671\n", " 31/1 1.02895\n", " 32/1 1.02679\n", " 33/1 1.02520\n", " 34/1 1.02888\n", " 35/1 1.02115\n", " 36/1 1.02557\n", " 37/1 1.02739\n", " 38/1 1.03153\n", " 39/1 1.03169\n", " 40/1 1.02154\n", " 41/1 1.02864\n", " 42/1 1.01558\n", " 43/1 1.02338\n", " 44/1 1.03512\n", " 45/1 1.02219\n", " 46/1 1.02973\n", " 47/1 1.02542\n", " 48/1 1.02512\n", " 49/1 1.02554\n", " 50/1 1.02434\n", " 51/1 1.02598\n", " 52/1 1.02526\n", " 53/1 1.02726\n", " 54/1 1.03077\n", " 55/1 1.02740\n", " 56/1 1.02722\n", " 57/1 1.02044\n", " 58/1 1.02819\n", " 59/1 1.02586\n", " 60/1 1.02692\n", " 61/1 1.02717\n", " 62/1 1.02398\n", " 63/1 1.03785\n", " 64/1 1.02292\n", " 65/1 1.03237\n", " 66/1 1.02918\n", " 67/1 1.02670\n", " 68/1 1.02403\n", " 69/1 1.02853\n", " 70/1 1.03072\n", " 71/1 1.02784\n", " 72/1 1.02805\n", " 73/1 1.02934\n", " 74/1 1.02622\n", " 75/1 1.02831\n", " 76/1 1.02457\n", " 77/1 1.02691\n", " 78/1 1.02340\n", " 79/1 1.02778\n", " 80/1 1.03008\n", " 81/1 1.02228\n", " 82/1 1.02696\n", " 83/1 1.03082\n", " 84/1 1.02785\n", " 85/1 1.02401\n", " 86/1 1.02613\n", " 87/1 1.02884\n", " 88/1 1.02725\n", " 89/1 1.02918\n", " 90/1 1.03261\n", " 91/1 1.02530\n", " 92/1 1.02692\n", " 93/1 1.02298\n", " 94/1 1.03199\n", " 95/1 1.03596\n", " 96/1 1.02777\n", " 97/1 1.02984\n", " 98/1 1.02997\n", " 99/1 1.03250\n", " 100/1 1.02953\n", " 101/1 1.02610\n", " 102/1 1.03681 1.03145 +/- 0.00535\n", " 103/1 1.03117 1.03136 +/- 0.00309\n", " 104/1 1.02337 1.02936 +/- 0.00296\n", " 105/1 1.02411 1.02831 +/- 0.00252\n", " 106/1 1.02606 1.02793 +/- 0.00209\n", " 107/1 1.03160 1.02846 +/- 0.00185\n", " 108/1 1.02655 1.02822 +/- 0.00162\n", " 109/1 1.02276 1.02761 +/- 0.00155\n", " 110/1 1.03534 1.02839 +/- 0.00159\n", " 111/1 1.02676 1.02824 +/- 0.00144\n", " 112/1 1.02659 1.02810 +/- 0.00132\n", " 113/1 1.02883 1.02816 +/- 0.00122\n", " 114/1 1.02528 1.02795 +/- 0.00115\n", " 115/1 1.02763 1.02793 +/- 0.00107\n", " 116/1 1.03310 1.02825 +/- 0.00105\n", " 117/1 1.03003 1.02836 +/- 0.00099\n", " 118/1 1.02453 1.02814 +/- 0.00096\n", " 119/1 1.03310 1.02841 +/- 0.00094\n", " 120/1 1.02153 1.02806 +/- 0.00096\n", " 121/1 1.02131 1.02774 +/- 0.00097\n", " 122/1 1.02138 1.02745 +/- 0.00097\n", " 123/1 1.02062 1.02715 +/- 0.00097\n", " 124/1 1.03122 1.02732 +/- 0.00094\n", " 125/1 1.02475 1.02722 +/- 0.00091\n", " 126/1 1.02481 1.02713 +/- 0.00088\n", " 127/1 1.03239 1.02732 +/- 0.00087\n", " 128/1 1.02615 1.02728 +/- 0.00084\n", " 129/1 1.02679 1.02726 +/- 0.00081\n", " 130/1 1.02996 1.02735 +/- 0.00079\n", " 131/1 1.02592 1.02731 +/- 0.00076\n", " 132/1 1.02891 1.02736 +/- 0.00074\n", " 133/1 1.02751 1.02736 +/- 0.00072\n", " 134/1 1.02770 1.02737 +/- 0.00070\n", " 135/1 1.02672 1.02735 +/- 0.00068\n", " 136/1 1.02679 1.02734 +/- 0.00066\n", " 137/1 1.03125 1.02744 +/- 0.00065\n", " 138/1 1.02852 1.02747 +/- 0.00063\n", " 139/1 1.01927 1.02726 +/- 0.00065\n", " 140/1 1.03625 1.02749 +/- 0.00067\n", " 141/1 1.02551 1.02744 +/- 0.00066\n", " 142/1 1.02082 1.02728 +/- 0.00066\n", " 143/1 1.02432 1.02721 +/- 0.00065\n", " 144/1 1.02108 1.02707 +/- 0.00065\n", " 145/1 1.03263 1.02720 +/- 0.00065\n", " 146/1 1.02775 1.02721 +/- 0.00063\n", " 147/1 1.03407 1.02735 +/- 0.00064\n", " 148/1 1.02499 1.02730 +/- 0.00062\n", " 149/1 1.02544 1.02727 +/- 0.00061\n", " 150/1 1.02688 1.02726 +/- 0.00060\n", " 151/1 1.02076 1.02713 +/- 0.00060\n", " 152/1 1.02576 1.02711 +/- 0.00059\n", " 153/1 1.01730 1.02692 +/- 0.00061\n", " 154/1 1.03201 1.02701 +/- 0.00060\n", " 155/1 1.02645 1.02700 +/- 0.00059\n", " 156/1 1.03299 1.02711 +/- 0.00059\n", " 157/1 1.02905 1.02714 +/- 0.00058\n", " 158/1 1.02989 1.02719 +/- 0.00057\n", " 159/1 1.03029 1.02724 +/- 0.00057\n", " 160/1 1.02166 1.02715 +/- 0.00057\n", " 161/1 1.02767 1.02716 +/- 0.00056\n", " 162/1 1.02701 1.02716 +/- 0.00055\n", " 163/1 1.03071 1.02721 +/- 0.00054\n", " 164/1 1.03338 1.02731 +/- 0.00054\n", " 165/1 1.03491 1.02743 +/- 0.00055\n", " 166/1 1.02669 1.02742 +/- 0.00054\n", " 167/1 1.02803 1.02743 +/- 0.00053\n", " 168/1 1.02544 1.02740 +/- 0.00052\n", " 169/1 1.02533 1.02737 +/- 0.00052\n", " 170/1 1.02373 1.02731 +/- 0.00051\n", " 171/1 1.02804 1.02732 +/- 0.00050\n", " 172/1 1.01677 1.02718 +/- 0.00052\n", " 173/1 1.03065 1.02723 +/- 0.00051\n", " 174/1 1.02981 1.02726 +/- 0.00051\n", " 175/1 1.03072 1.02731 +/- 0.00050\n", " 176/1 1.02618 1.02729 +/- 0.00050\n", " 177/1 1.03149 1.02735 +/- 0.00049\n", " 178/1 1.02972 1.02738 +/- 0.00049\n", " 179/1 1.02105 1.02730 +/- 0.00049\n", " 180/1 1.02514 1.02727 +/- 0.00048\n", " 181/1 1.03802 1.02740 +/- 0.00049\n", " 182/1 1.02964 1.02743 +/- 0.00049\n", " 183/1 1.02638 1.02742 +/- 0.00048\n", " 184/1 1.02662 1.02741 +/- 0.00048\n", " 185/1 1.02524 1.02738 +/- 0.00047\n", " 186/1 1.02059 1.02730 +/- 0.00047\n", " 187/1 1.02898 1.02732 +/- 0.00047\n", " 188/1 1.02537 1.02730 +/- 0.00046\n", " 189/1 1.02540 1.02728 +/- 0.00046\n", " 190/1 1.02328 1.02723 +/- 0.00046\n", " 191/1 1.02232 1.02718 +/- 0.00045\n", " 192/1 1.03147 1.02723 +/- 0.00045\n", " 193/1 1.02436 1.02720 +/- 0.00045\n", " 194/1 1.02453 1.02717 +/- 0.00044\n", " 195/1 1.02902 1.02719 +/- 0.00044\n", " 196/1 1.02302 1.02714 +/- 0.00044\n", " 197/1 1.02436 1.02712 +/- 0.00043\n", " 198/1 1.02165 1.02706 +/- 0.00043\n", " 199/1 1.02619 1.02705 +/- 0.00043\n", " 200/1 1.02406 1.02702 +/- 0.00043\n", " 201/1 1.02100 1.02696 +/- 0.00043\n", " 202/1 1.02359 1.02693 +/- 0.00042\n", " 203/1 1.02760 1.02693 +/- 0.00042\n", " 204/1 1.02722 1.02694 +/- 0.00041\n", " 205/1 1.03250 1.02699 +/- 0.00041\n", " 206/1 1.02241 1.02695 +/- 0.00041\n", " 207/1 1.02682 1.02695 +/- 0.00041\n", " 208/1 1.03065 1.02698 +/- 0.00041\n", " 209/1 1.02462 1.02696 +/- 0.00040\n", " 210/1 1.02561 1.02695 +/- 0.00040\n", " 211/1 1.02837 1.02696 +/- 0.00040\n", " 212/1 1.02689 1.02696 +/- 0.00039\n", " 213/1 1.03487 1.02703 +/- 0.00040\n", " 214/1 1.02529 1.02701 +/- 0.00039\n", " 215/1 1.02410 1.02699 +/- 0.00039\n", " 216/1 1.02708 1.02699 +/- 0.00039\n", " 217/1 1.02235 1.02695 +/- 0.00038\n", " 218/1 1.02515 1.02693 +/- 0.00038\n", " 219/1 1.02664 1.02693 +/- 0.00038\n", " 220/1 1.02326 1.02690 +/- 0.00038\n", " 221/1 1.02513 1.02689 +/- 0.00037\n", " 222/1 1.02468 1.02687 +/- 0.00037\n", " 223/1 1.02577 1.02686 +/- 0.00037\n", " 224/1 1.03120 1.02689 +/- 0.00037\n", " 225/1 1.01948 1.02683 +/- 0.00037\n", " 226/1 1.03144 1.02687 +/- 0.00037\n", " 227/1 1.02476 1.02685 +/- 0.00037\n", " 228/1 1.02805 1.02686 +/- 0.00036\n", " 229/1 1.02893 1.02688 +/- 0.00036\n", " 230/1 1.02462 1.02686 +/- 0.00036\n", " 231/1 1.02394 1.02684 +/- 0.00036\n", " 232/1 1.02065 1.02679 +/- 0.00036\n", " 233/1 1.02684 1.02679 +/- 0.00035\n", " 234/1 1.02661 1.02679 +/- 0.00035\n", " 235/1 1.03140 1.02683 +/- 0.00035\n", " 236/1 1.02658 1.02682 +/- 0.00035\n", " 237/1 1.02044 1.02678 +/- 0.00035\n", " 238/1 1.02493 1.02676 +/- 0.00035\n", " 239/1 1.02685 1.02677 +/- 0.00034\n", " 240/1 1.02868 1.02678 +/- 0.00034\n", " 241/1 1.02988 1.02680 +/- 0.00034\n", " 242/1 1.01644 1.02673 +/- 0.00034\n", " 243/1 1.02573 1.02672 +/- 0.00034\n", " 244/1 1.02592 1.02672 +/- 0.00034\n", " 245/1 1.03339 1.02676 +/- 0.00034\n", " 246/1 1.02469 1.02675 +/- 0.00034\n", " 247/1 1.02459 1.02673 +/- 0.00034\n", " 248/1 1.02925 1.02675 +/- 0.00033\n", " 249/1 1.02403 1.02673 +/- 0.00033\n", " 250/1 1.03013 1.02675 +/- 0.00033\n", " 251/1 1.02203 1.02672 +/- 0.00033\n", " 252/1 1.02794 1.02673 +/- 0.00033\n", " 253/1 1.02799 1.02674 +/- 0.00033\n", " 254/1 1.02791 1.02675 +/- 0.00032\n", " 255/1 1.03370 1.02679 +/- 0.00033\n", " 256/1 1.01831 1.02674 +/- 0.00033\n", " 257/1 1.03075 1.02676 +/- 0.00033\n", " 258/1 1.02733 1.02677 +/- 0.00032\n", " 259/1 1.02347 1.02675 +/- 0.00032\n", " 260/1 1.02169 1.02671 +/- 0.00032\n", " 261/1 1.02858 1.02673 +/- 0.00032\n", " 262/1 1.02725 1.02673 +/- 0.00032\n", " 263/1 1.02762 1.02673 +/- 0.00032\n", " 264/1 1.02390 1.02672 +/- 0.00032\n", " 265/1 1.02794 1.02672 +/- 0.00031\n", " 266/1 1.01782 1.02667 +/- 0.00032\n", " 267/1 1.02576 1.02667 +/- 0.00031\n", " 268/1 1.02614 1.02666 +/- 0.00031\n", " 269/1 1.02912 1.02668 +/- 0.00031\n", " 270/1 1.02348 1.02666 +/- 0.00031\n", " 271/1 1.02954 1.02667 +/- 0.00031\n", " 272/1 1.02743 1.02668 +/- 0.00031\n", " 273/1 1.02227 1.02665 +/- 0.00031\n", " 274/1 1.02673 1.02665 +/- 0.00030\n", " 275/1 1.02707 1.02666 +/- 0.00030\n", " 276/1 1.02359 1.02664 +/- 0.00030\n", " 277/1 1.02918 1.02665 +/- 0.00030\n", " 278/1 1.03028 1.02667 +/- 0.00030\n", " 279/1 1.03577 1.02672 +/- 0.00030\n", " 280/1 1.03646 1.02678 +/- 0.00031\n", " 281/1 1.03296 1.02681 +/- 0.00031\n", " 282/1 1.02302 1.02679 +/- 0.00030\n", " 283/1 1.02979 1.02681 +/- 0.00030\n", " 284/1 1.03217 1.02684 +/- 0.00030\n", " 285/1 1.02901 1.02685 +/- 0.00030\n", " 286/1 1.02483 1.02684 +/- 0.00030\n", " 287/1 1.02961 1.02685 +/- 0.00030\n", " 288/1 1.02835 1.02686 +/- 0.00030\n", " 289/1 1.02680 1.02686 +/- 0.00030\n", " 290/1 1.02673 1.02686 +/- 0.00029\n", " 291/1 1.01969 1.02682 +/- 0.00029\n", " 292/1 1.02122 1.02679 +/- 0.00029\n", " 293/1 1.02975 1.02681 +/- 0.00029\n", " 294/1 1.02754 1.02681 +/- 0.00029\n", " 295/1 1.02851 1.02682 +/- 0.00029\n", " 296/1 1.02779 1.02683 +/- 0.00029\n", " 297/1 1.02658 1.02682 +/- 0.00029\n", " 298/1 1.02947 1.02684 +/- 0.00029\n", " 299/1 1.02843 1.02685 +/- 0.00029\n", " 300/1 1.02923 1.02686 +/- 0.00028\n", " Creating state point statepoint.300.h5...\n", "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", " Total time for initialization = 5.4904e+00 seconds\n", " Reading cross sections = 5.3396e+00 seconds\n", " Total time in simulation = 5.7751e+02 seconds\n", " Time in transport only = 5.6785e+02 seconds\n", " Time in inactive batches = 1.3790e+02 seconds\n", " Time in active batches = 4.3961e+02 seconds\n", " Time synchronizing fission bank = 7.4705e+00 seconds\n", " Sampling source sites = 7.0682e+00 seconds\n", " SEND/RECV source sites = 4.0202e-01 seconds\n", " Time accumulating tallies = 4.4668e-01 seconds\n", " Time writing statepoints = 2.9013e-02 seconds\n", " Total time for finalization = 4.7723e-01 seconds\n", " Total time elapsed = 5.8362e+02 seconds\n", " Calculation Rate (inactive) = 87020.5 particles/second\n", " Calculation Rate (active) = 54593.9 particles/second\n", "\n", " ============================> RESULTS <============================\n", "\n", " k-effective (Collision) = 1.02699 +/- 0.00028\n", " k-effective (Track-length) = 1.02686 +/- 0.00028\n", " k-effective (Absorption) = 1.02758 +/- 0.00024\n", " Combined k-effective = 1.02733 +/- 0.00021\n", " Leakage Fraction = 0.10106 +/- 0.00006\n", "\n" ] } ], "source": [ "openmc.run()" ] }, { "cell_type": "code", "execution_count": 46, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0" ] }, "execution_count": 46, "metadata": {}, "output_type": "execute_result" } ], "source": [ "os.system(\"mv statepoint.300.h5 ./h5_files/\" + str(z0_reg-offset)+\"_\"+ str(z0_shim-offset) + \"_\" + str(z0_trans-offset)+\".h5\")" ] } ], "metadata": { "kernelspec": { "display_name": "mc", "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.11.9" } }, "nbformat": 4, "nbformat_minor": 2 }