Updates to unstructured mesh example notebooks.

This commit is contained in:
Patrick Shriwise 2020-03-27 11:39:12 -05:00
parent 1716d35ce5
commit 25fe1ef4f3
2 changed files with 97 additions and 960 deletions

View file

@ -35,7 +35,7 @@
},
{
"cell_type": "code",
"execution_count": 3,
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
@ -70,7 +70,7 @@
},
{
"cell_type": "code",
"execution_count": 4,
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
@ -86,7 +86,7 @@
},
{
"cell_type": "code",
"execution_count": 5,
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
@ -109,7 +109,7 @@
},
{
"cell_type": "code",
"execution_count": 6,
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
@ -130,16 +130,16 @@
},
{
"cell_type": "code",
"execution_count": 7,
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7fd60775e080>"
"<matplotlib.image.AxesImage at 0x7fe9650cd320>"
]
},
"execution_count": 7,
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
@ -169,16 +169,16 @@
},
{
"cell_type": "code",
"execution_count": 8,
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7fd609966588>"
"<matplotlib.image.AxesImage at 0x7fe967252e48>"
]
},
"execution_count": 8,
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
},
@ -213,7 +213,7 @@
},
{
"cell_type": "code",
"execution_count": 9,
"execution_count": 8,
"metadata": {},
"outputs": [
{
@ -248,9 +248,9 @@
" Copyright | 2011-2020 MIT and OpenMC contributors\n",
" License | http://openmc.readthedocs.io/en/latest/license.html\n",
" Version | 0.12.0-dev\n",
" Git SHA1 | 34da5a1e2a3942428dfca6e8dfc2ce02d04333dc\n",
" Date/Time | 2020-03-19 09:53:33\n",
" OpenMP Threads | 8\n",
" Git SHA1 | e7eceff2aa3a9bbfd4dd553f585d1a31b051ddb3\n",
" Date/Time | 2020-03-27 11:18:32\n",
" OpenMP Threads | 2\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
@ -272,7 +272,6 @@
" Maximum neutron transport energy: 20000000.000000 eV for U235\n",
" Minimum neutron data temperature: 294.000000 K\n",
" Maximum neutron data temperature: 294.000000 K\n",
" Reading tallies XML file...\n",
" Preparing distributed cell instances...\n",
" Writing summary.h5 file...\n",
" Initializing source particles...\n",
@ -292,24 +291,23 @@
" 9/1 0.29029 0.23536 +/- 0.01935\n",
" 10/1 0.20094 0.22848 +/- 0.01649\n",
" Creating state point statepoint.10.h5...\n",
" Writing unstructured mesh tally_2.10.vtk...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 2.1114e+02 seconds\n",
" Reading cross sections = 9.1607e-01 seconds\n",
" Total time in simulation = 4.1149e+00 seconds\n",
" Time in transport only = 2.6868e-02 seconds\n",
" Time in inactive batches = 4.2519e-03 seconds\n",
" Time in active batches = 4.1106e+00 seconds\n",
" Time synchronizing fission bank = 4.4463e-05 seconds\n",
" Sampling source sites = 3.7916e-05 seconds\n",
" SEND/RECV source sites = 3.5140e-06 seconds\n",
" Time accumulating tallies = 9.2627e-03 seconds\n",
" Total time for finalization = 5.3702e-01 seconds\n",
" Total time elapsed = 2.1580e+02 seconds\n",
" Calculation Rate (inactive) = 117594.0 particles/second\n",
" Calculation Rate (active) = 121.637 particles/second\n",
" Total time for initialization = 6.8458e-01 seconds\n",
" Reading cross sections = 6.7039e-01 seconds\n",
" Total time in simulation = 1.9462e-02 seconds\n",
" Time in transport only = 1.7538e-02 seconds\n",
" Time in inactive batches = 9.2816e-03 seconds\n",
" Time in active batches = 1.0181e-02 seconds\n",
" Time synchronizing fission bank = 5.2594e-05 seconds\n",
" Sampling source sites = 4.6704e-05 seconds\n",
" SEND/RECV source sites = 3.1580e-06 seconds\n",
" Time accumulating tallies = 2.0290e-06 seconds\n",
" Total time for finalization = 1.0320e-06 seconds\n",
" Total time elapsed = 7.0439e-01 seconds\n",
" Calculation Rate (inactive) = 53869.9 particles/second\n",
" Calculation Rate (active) = 49113.3 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
@ -324,10 +322,10 @@
{
"data": {
"text/plain": [
"0.2070678896718187+/-0.019653037754288158"
"0.20706788967181938+/-0.019653037754288005"
]
},
"execution_count": 9,
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
@ -345,7 +343,7 @@
},
{
"cell_type": "code",
"execution_count": 10,
"execution_count": 9,
"metadata": {},
"outputs": [
{
@ -355,7 +353,7 @@
"<IPython.core.display.Image object>"
]
},
"execution_count": 10,
"execution_count": 9,
"metadata": {
"image/png": {
"width": 600
@ -379,7 +377,7 @@
},
{
"cell_type": "code",
"execution_count": 11,
"execution_count": 10,
"metadata": {},
"outputs": [
{
@ -434,7 +432,7 @@
},
{
"cell_type": "code",
"execution_count": 12,
"execution_count": 11,
"metadata": {},
"outputs": [],
"source": [
@ -452,7 +450,7 @@
},
{
"cell_type": "code",
"execution_count": 13,
"execution_count": 12,
"metadata": {},
"outputs": [],
"source": [
@ -474,218 +472,16 @@
},
{
"cell_type": "code",
"execution_count": 14,
"execution_count": 13,
"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-2020 MIT and OpenMC contributors\n",
" License | http://openmc.readthedocs.io/en/latest/license.html\n",
" Version | 0.12.0-dev\n",
" Git SHA1 | 34da5a1e2a3942428dfca6e8dfc2ce02d04333dc\n",
" Date/Time | 2020-03-19 09:57:45\n",
" OpenMP Threads | 8\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
" Reading materials XML file...\n",
" Reading geometry XML file...\n",
" Reading U234 from /home/shriwise/opt/openmc/xs/nndc_hdf5/U234.h5\n",
" Reading U235 from /home/shriwise/opt/openmc/xs/nndc_hdf5/U235.h5\n",
" Reading U238 from /home/shriwise/opt/openmc/xs/nndc_hdf5/U238.h5\n",
" Reading O16 from /home/shriwise/opt/openmc/xs/nndc_hdf5/O16.h5\n",
" Reading Zr90 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Zr90.h5\n",
" Reading Zr91 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Zr91.h5\n",
" Reading Zr92 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Zr92.h5\n",
" Reading Zr94 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Zr94.h5\n",
" Reading Zr96 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Zr96.h5\n",
" Reading H1 from /home/shriwise/opt/openmc/xs/nndc_hdf5/H1.h5\n",
" Reading B10 from /home/shriwise/opt/openmc/xs/nndc_hdf5/B10.h5\n",
" Reading B11 from /home/shriwise/opt/openmc/xs/nndc_hdf5/B11.h5\n",
" Reading c_H_in_H2O from /home/shriwise/opt/openmc/xs/nndc_hdf5/c_H_in_H2O.h5\n",
" Maximum neutron transport energy: 20000000.000000 eV for U235\n",
" Minimum neutron data temperature: 294.000000 K\n",
" Maximum neutron data temperature: 294.000000 K\n",
" Reading tallies XML file...\n",
" Preparing distributed cell instances...\n",
" Writing summary.h5 file...\n",
" Initializing source particles...\n",
"\n",
" ====================> K EIGENVALUE SIMULATION <====================\n",
"\n",
" Bat./Gen. k Average k\n",
" ========= ======== ====================\n",
" 1/1 0.22539\n",
" 2/1 0.23030\n",
" 3/1 0.23180\n",
" 4/1 0.23343\n",
" 5/1 0.22940\n",
" 6/1 0.22765\n",
" 7/1 0.23238\n",
" 8/1 0.23083\n",
" 9/1 0.23204\n",
" 10/1 0.23189\n",
" 11/1 0.23555\n",
" 12/1 0.23220\n",
" 13/1 0.22907\n",
" 14/1 0.23016\n",
" 15/1 0.23055\n",
" 16/1 0.23062\n",
" 17/1 0.22656\n",
" 18/1 0.23380\n",
" 19/1 0.23268\n",
" 20/1 0.23233\n",
" 21/1 0.23256\n",
" 22/1 0.22965 0.23111 +/- 0.00145\n",
" 23/1 0.22846 0.23022 +/- 0.00121\n",
" 24/1 0.22936 0.23001 +/- 0.00089\n",
" 25/1 0.23054 0.23012 +/- 0.00069\n",
" 26/1 0.22708 0.22961 +/- 0.00076\n",
" 27/1 0.23159 0.22989 +/- 0.00070\n",
" WARNING: No tet found for location between triangle hits\n",
" 28/1 0.23703 0.23078 +/- 0.00108\n",
" 29/1 0.23227 0.23095 +/- 0.00097\n",
" 30/1 0.23114 0.23097 +/- 0.00086\n",
" 31/1 0.23133 0.23100 +/- 0.00078\n",
" 32/1 0.23143 0.23104 +/- 0.00072\n",
" 33/1 0.23125 0.23105 +/- 0.00066\n",
" 34/1 0.23218 0.23113 +/- 0.00061\n",
" 35/1 0.23140 0.23115 +/- 0.00057\n",
" 36/1 0.22911 0.23102 +/- 0.00055\n",
" 37/1 0.23143 0.23105 +/- 0.00052\n",
" 38/1 0.23342 0.23118 +/- 0.00051\n",
" 39/1 0.23186 0.23121 +/- 0.00048\n",
" 40/1 0.23029 0.23117 +/- 0.00046\n",
" 41/1 0.23132 0.23118 +/- 0.00043\n",
" 42/1 0.23167 0.23120 +/- 0.00042\n",
" 43/1 0.23244 0.23125 +/- 0.00040\n",
" 44/1 0.23101 0.23124 +/- 0.00038\n",
" 45/1 0.23225 0.23128 +/- 0.00037\n",
" 46/1 0.22945 0.23121 +/- 0.00036\n",
" 47/1 0.22978 0.23116 +/- 0.00035\n",
" 48/1 0.23335 0.23124 +/- 0.00035\n",
" 49/1 0.23298 0.23130 +/- 0.00034\n",
" 50/1 0.23095 0.23129 +/- 0.00033\n",
" 51/1 0.23724 0.23148 +/- 0.00037\n",
" 52/1 0.22973 0.23142 +/- 0.00037\n",
" 53/1 0.23066 0.23140 +/- 0.00035\n",
" 54/1 0.22838 0.23131 +/- 0.00036\n",
" 55/1 0.23262 0.23135 +/- 0.00035\n",
" 56/1 0.23593 0.23148 +/- 0.00036\n",
" 57/1 0.23358 0.23153 +/- 0.00036\n",
" 58/1 0.23050 0.23151 +/- 0.00035\n",
" 59/1 0.23273 0.23154 +/- 0.00034\n",
" 60/1 0.22842 0.23146 +/- 0.00034\n",
" 61/1 0.23344 0.23151 +/- 0.00033\n",
" 62/1 0.23333 0.23155 +/- 0.00033\n",
" 63/1 0.22987 0.23151 +/- 0.00032\n",
" 64/1 0.23117 0.23150 +/- 0.00032\n",
" 65/1 0.23197 0.23151 +/- 0.00031\n",
" 66/1 0.23379 0.23156 +/- 0.00031\n",
" 67/1 0.23461 0.23163 +/- 0.00031\n",
" 68/1 0.23109 0.23162 +/- 0.00030\n",
" 69/1 0.22916 0.23157 +/- 0.00030\n",
" 70/1 0.23008 0.23154 +/- 0.00029\n",
" 71/1 0.23157 0.23154 +/- 0.00029\n",
" 72/1 0.23126 0.23153 +/- 0.00028\n",
" 73/1 0.23377 0.23157 +/- 0.00028\n",
" 74/1 0.23105 0.23157 +/- 0.00028\n",
" 75/1 0.23654 0.23166 +/- 0.00029\n",
" 76/1 0.23198 0.23166 +/- 0.00028\n",
" 77/1 0.23390 0.23170 +/- 0.00028\n",
" 78/1 0.23455 0.23175 +/- 0.00028\n",
" 79/1 0.23245 0.23176 +/- 0.00027\n",
" 80/1 0.23121 0.23175 +/- 0.00027\n",
" 81/1 0.23183 0.23175 +/- 0.00026\n",
" 82/1 0.23496 0.23181 +/- 0.00027\n",
" 83/1 0.22763 0.23174 +/- 0.00027\n",
" 84/1 0.23184 0.23174 +/- 0.00027\n",
" 85/1 0.23074 0.23173 +/- 0.00026\n",
" 86/1 0.23178 0.23173 +/- 0.00026\n",
" 87/1 0.23135 0.23172 +/- 0.00025\n",
" 88/1 0.23117 0.23171 +/- 0.00025\n",
" 89/1 0.22815 0.23166 +/- 0.00025\n",
" 90/1 0.22852 0.23162 +/- 0.00025\n",
" 91/1 0.22910 0.23158 +/- 0.00025\n",
" 92/1 0.23143 0.23158 +/- 0.00025\n",
" 93/1 0.23097 0.23157 +/- 0.00024\n",
" 94/1 0.23348 0.23160 +/- 0.00024\n",
" 95/1 0.23068 0.23158 +/- 0.00024\n",
" 96/1 0.23089 0.23157 +/- 0.00024\n",
" 97/1 0.23373 0.23160 +/- 0.00024\n",
" 98/1 0.23336 0.23163 +/- 0.00023\n",
" 99/1 0.23084 0.23162 +/- 0.00023\n",
" 100/1 0.23116 0.23161 +/- 0.00023\n",
" Creating state point statepoint.100.h5...\n",
" Writing unstructured mesh tally_1.100.vtk...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 2.0470e+02 seconds\n",
" Reading cross sections = 6.4733e-01 seconds\n",
" Total time in simulation = 3.5637e+02 seconds\n",
" Time in transport only = 3.5039e+02 seconds\n",
" Time in inactive batches = 1.8166e+01 seconds\n",
" Time in active batches = 3.3821e+02 seconds\n",
" Time synchronizing fission bank = 6.9102e-01 seconds\n",
" Sampling source sites = 6.0920e-01 seconds\n",
" SEND/RECV source sites = 8.1712e-02 seconds\n",
" Time accumulating tallies = 2.2190e-01 seconds\n",
" Total time for finalization = 8.3563e-01 seconds\n",
" Total time elapsed = 5.6214e+02 seconds\n",
" Calculation Rate (inactive) = 110098.0 particles/second\n",
" Calculation Rate (active) = 23654.1 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
" k-effective (Collision) = 0.23129 +/- 0.00020\n",
" k-effective (Track-length) = 0.23161 +/- 0.00023\n",
" k-effective (Absorption) = 0.23099 +/- 0.00019\n",
" Combined k-effective = 0.23119 +/- 0.00017\n",
" Leakage Fraction = 0.79477 +/- 0.00014\n",
"\n"
]
},
{
"data": {
"text/plain": [
"0.23118721242921098+/-0.0001701893231378034"
"0.2311872124292096+/-0.00017018932313754055"
]
},
"execution_count": 14,
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
@ -694,7 +490,7 @@
"model.settings.particles = 100_000\n",
"model.settings.inactive = 20\n",
"model.settings.batches = 100\n",
"model.run()"
"model.run(output=False)"
]
},
{
@ -718,14 +514,14 @@
},
{
"cell_type": "code",
"execution_count": 15,
"execution_count": 14,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"manifold_flux.vtk tally_1.100.vtk tally_1.200.vtk tally_2.10.vtk\r\n"
"tally_1.100.vtk\r\n"
]
}
],
@ -742,7 +538,7 @@
},
{
"cell_type": "code",
"execution_count": 16,
"execution_count": 15,
"metadata": {},
"outputs": [
{
@ -752,7 +548,7 @@
"<IPython.core.display.Image object>"
]
},
"execution_count": 16,
"execution_count": 15,
"metadata": {
"image/png": {
"width": 600
@ -775,7 +571,7 @@
},
{
"cell_type": "code",
"execution_count": 17,
"execution_count": 16,
"metadata": {},
"outputs": [
{
@ -785,7 +581,7 @@
"<IPython.core.display.Image object>"
]
},
"execution_count": 17,
"execution_count": 16,
"metadata": {
"image/png": {
"width": 600
@ -807,7 +603,7 @@
},
{
"cell_type": "code",
"execution_count": 18,
"execution_count": 17,
"metadata": {},
"outputs": [],
"source": [
@ -825,20 +621,15 @@
},
{
"cell_type": "code",
"execution_count": 19,
"execution_count": 18,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[1.43381086e-04]\n",
" [1.48043747e-04]\n",
" [1.60408339e-04]\n",
" ...\n",
" [7.04197023e-05]\n",
" [7.04197023e-05]\n",
" [7.04197023e-05]]\n",
"[1.43381086e-04 1.48043747e-04 1.60408339e-04 ... 7.04197023e-05\n",
" 7.04197023e-05 7.04197023e-05]\n",
"[[ 2.88485691 -2.55429784 9.97768184]\n",
" [ 2.87565092 -2.60469781 9.8884092 ]\n",
" [ 2.85832254 -2.65291228 9.97768184]\n",
@ -858,19 +649,14 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"The combination of these values can provide for an appoxmiate visualization of the unstructured mesh without its explicit representation or use of an additional mesh library."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The combination of these values can provide for an appoxmiate visualization of the unstructured mesh without its explicit representation or use of an additional mesh library.\n",
"\n",
"We hope you've found this example notebook useful!"
]
},
{
"cell_type": "code",
"execution_count": 20,
"execution_count": 19,
"metadata": {},
"outputs": [
{
@ -880,7 +666,7 @@
"<IPython.core.display.Image object>"
]
},
"execution_count": 20,
"execution_count": 19,
"metadata": {
"image/png": {
"width": 600

View file

@ -40,8 +40,8 @@
"source": [
"import urllib.request\n",
"\n",
"manifold_geom_url = 'https://tinyurl.com/rp7grox'\n",
"manifold_mesh_url = 'https://tinyurl.com/wojemuh'\n",
"manifold_geom_url = 'https://tinyurl.com/rp7grox' # 99 MB\n",
"manifold_mesh_url = 'https://tinyurl.com/wojemuh' # 5.4 MB\n",
"\n",
"def download(url, filename='dagmc.h5m'):\n",
" \"\"\"\n",
@ -108,53 +108,41 @@
"source": [
"air = openmc.Material(name='air')\n",
"air.set_density('g/cc', 0.001205)\n",
"air.add_nuclide('N14',0.781557629247)\n",
"air.add_nuclide('N15',0.002873370753)\n",
"air.add_nuclide('O16',0.210668126508)\n",
"air.add_nuclide('O17',7.9873492e-05)\n",
"air.add_nuclide('Ar36',1.53456e-05)\n",
"air.add_nuclide('Ar38',2.8934e-06)\n",
"air.add_nuclide('Ar40',0.004581761)\n",
"air.add_element('N', 0.784431)\n",
"air.add_element('O', 0.210748)\n",
"air.add_element('Ar',0.0046)\n",
"\n",
"steel = openmc.Material(name='steel')\n",
"steel.set_density('g/cc', 8.0)\n",
"steel.add_nuclide('Si28',0.0092672382464)\n",
"steel.add_nuclide('Si29',0.00047056391679999997)\n",
"steel.add_nuclide('Si30',0.00031019783679999996)\n",
"steel.add_element('Si', 0.010048)\n",
"steel.add_element('S', 0.00023)\n",
"steel.add_element('Fe', 0.669)\n",
"steel.add_element('Ni', 0.12)\n",
"steel.add_element('Mo', 0.025)\n",
"steel.add_nuclide('P31',0.00023)\n",
"steel.add_nuclide('S32',0.000218593702)\n",
"steel.add_nuclide('S33',1.721987e-06)\n",
"steel.add_nuclide('S34',9.650777000000001e-06)\n",
"steel.add_nuclide('S36',3.3534e-08)\n",
"steel.add_nuclide('Mn55',0.011014)\n",
"steel.add_nuclide('Fe54',0.03910305)\n",
"steel.add_nuclide('Fe56',0.6138342600000001)\n",
"steel.add_nuclide('Fe57',0.01417611)\n",
"steel.add_nuclide('Fe58',0.0018865800000000001)\n",
"steel.add_nuclide('Ni58',0.08169227999999999)\n",
"steel.add_nuclide('Ni60',0.03146772)\n",
"steel.add_nuclide('Ni61',0.00136788)\n",
"steel.add_nuclide('Ni62',0.0043614000000000005)\n",
"steel.add_nuclide('Ni64',0.00111072)\n",
"steel.add_nuclide('Mo100',0.0024360000000000002)\n",
"steel.add_nuclide('Mo92',0.0036622500000000006)\n",
"steel.add_nuclide('Mo94',0.0022967499999999997)\n",
"steel.add_nuclide('Mo95',0.00396825)\n",
"steel.add_nuclide('Mo96',0.00416825)\n",
"steel.add_nuclide('Mo97',0.0023955)\n",
"steel.add_nuclide('Mo98',0.006073)\n",
"\n",
"materials = openmc.Materials([air, steel])\n",
"materials.export_to_xml()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now let's download the geometry and mesh files.\n",
"(This may take some time.)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"download(manifold_geom_url)\n",
"# get the manifold DAGMC geometry file\n",
"download(manifold_geom_url) \n",
"# get the manifold tet mesh\n",
"download(manifold_mesh_url, 'manifold.h5m')"
]
},
@ -246,9 +234,9 @@
" Copyright | 2011-2020 MIT and OpenMC contributors\n",
" License | http://openmc.readthedocs.io/en/latest/license.html\n",
" Version | 0.12.0-dev\n",
" Git SHA1 | 34da5a1e2a3942428dfca6e8dfc2ce02d04333dc\n",
" Date/Time | 2020-03-19 10:10:21\n",
" OpenMP Threads | 8\n",
" Git SHA1 | e7eceff2aa3a9bbfd4dd553f585d1a31b051ddb3\n",
" Date/Time | 2020-03-27 10:25:39\n",
" OpenMP Threads | 2\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
@ -269,12 +257,10 @@
" Reading Si28 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si28.h5\n",
" Reading Si29 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si29.h5\n",
" Reading Si30 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si30.h5\n",
" Reading P31 from /home/shriwise/opt/openmc/xs/nndc_hdf5/P31.h5\n",
" Reading S32 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S32.h5\n",
" Reading S33 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S33.h5\n",
" Reading S34 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S34.h5\n",
" Reading S36 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S36.h5\n",
" Reading Mn55 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mn55.h5\n",
" Reading Fe54 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe54.h5\n",
" Reading Fe56 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe56.h5\n",
" Reading Fe57 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe57.h5\n",
@ -291,13 +277,13 @@
" Reading Mo96 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo96.h5\n",
" Reading Mo97 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo97.h5\n",
" Reading Mo98 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo98.h5\n",
" Reading P31 from /home/shriwise/opt/openmc/xs/nndc_hdf5/P31.h5\n",
" Reading Mn55 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mn55.h5\n",
" Maximum neutron transport energy: 20000000.000000 eV for N15\n",
" Minimum neutron data temperature: 294.000000 K\n",
" Maximum neutron data temperature: 294.000000 K\n",
" Reading tallies XML file...\n",
" Preparing distributed cell instances...\n",
" Writing summary.h5 file...\n",
" Initializing source particles...\n",
"\n",
" ===============> FIXED SOURCE TRANSPORT SIMULATION <===============\n",
"\n",
@ -312,24 +298,22 @@
" Simulating batch 9\n",
" Simulating batch 10\n",
" Creating state point statepoint.10.h5...\n",
" Writing unstructured mesh tally_1.10.vtk...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 2.2053e+02 seconds\n",
" Reading cross sections = 2.1899e+00 seconds\n",
" Total time in simulation = 1.4710e+01 seconds\n",
" Time in transport only = 6.4195e+00 seconds\n",
" Time in active batches = 1.4710e+01 seconds\n",
" Time sampling source = 3.6842e+00 seconds\n",
" Time accumulating tallies = 1.9801e-02 seconds\n",
" Total time for finalization = 6.9415e-01 seconds\n",
" Total time elapsed = 2.3623e+02 seconds\n",
" Calculation Rate (active) = 3399.10 particles/second\n",
" Total time for initialization = 2.0935e+01 seconds\n",
" Reading cross sections = 2.3075e+00 seconds\n",
" Total time in simulation = 7.0885e+00 seconds\n",
" Time in transport only = 7.0872e+00 seconds\n",
" Time in active batches = 7.0885e+00 seconds\n",
" Time accumulating tallies = 3.8640e-06 seconds\n",
" Total time for finalization = 1.2420e-06 seconds\n",
" Total time elapsed = 2.8024e+01 seconds\n",
" Calculation Rate (active) = 7053.66 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" Leakage Fraction = 0.97422 +/- 0.00091\n",
" Leakage Fraction = 0.97472 +/- 0.00083\n",
"\n"
]
}
@ -342,7 +326,7 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"Now let's setup the unstructured mesh tally. We'll do this the same way we did in the previous notebook."
"Now let's setup the unstructured mesh tally. We'll do this the same way we did in the [previous notebook](./unstructured-mesh-part-i.ipynb)."
]
},
{
@ -379,332 +363,16 @@
"cell_type": "code",
"execution_count": 11,
"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-2020 MIT and OpenMC contributors\n",
" License | http://openmc.readthedocs.io/en/latest/license.html\n",
" Version | 0.12.0-dev\n",
" Git SHA1 | 34da5a1e2a3942428dfca6e8dfc2ce02d04333dc\n",
" Date/Time | 2020-03-19 10:14:19\n",
" OpenMP Threads | 8\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
" Reading materials XML file...\n",
" Reading DAGMC geometry...\n",
"Loading file dagmc.h5m\n",
"Initializing the GeomQueryTool...\n",
"Using faceting tolerance: 0.001\n",
"Building OBB Tree...\n",
" Reading N14 from /home/shriwise/opt/openmc/xs/nndc_hdf5/N14.h5\n",
" Reading N15 from /home/shriwise/opt/openmc/xs/nndc_hdf5/N15.h5\n",
" Reading O16 from /home/shriwise/opt/openmc/xs/nndc_hdf5/O16.h5\n",
" Reading O17 from /home/shriwise/opt/openmc/xs/nndc_hdf5/O17.h5\n",
" Reading Ar36 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ar36.h5\n",
" WARNING: Negative value(s) found on probability table for nuclide Ar36 at 294K\n",
" Reading Ar38 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ar38.h5\n",
" Reading Ar40 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ar40.h5\n",
" Reading Si28 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si28.h5\n",
" Reading Si29 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si29.h5\n",
" Reading Si30 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si30.h5\n",
" Reading P31 from /home/shriwise/opt/openmc/xs/nndc_hdf5/P31.h5\n",
" Reading S32 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S32.h5\n",
" Reading S33 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S33.h5\n",
" Reading S34 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S34.h5\n",
" Reading S36 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S36.h5\n",
" Reading Mn55 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mn55.h5\n",
" Reading Fe54 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe54.h5\n",
" Reading Fe56 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe56.h5\n",
" Reading Fe57 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe57.h5\n",
" Reading Fe58 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe58.h5\n",
" Reading Ni58 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni58.h5\n",
" Reading Ni60 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni60.h5\n",
" Reading Ni61 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni61.h5\n",
" Reading Ni62 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni62.h5\n",
" Reading Ni64 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni64.h5\n",
" Reading Mo100 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo100.h5\n",
" Reading Mo92 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo92.h5\n",
" Reading Mo94 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo94.h5\n",
" Reading Mo95 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo95.h5\n",
" Reading Mo96 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo96.h5\n",
" Reading Mo97 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo97.h5\n",
" Reading Mo98 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo98.h5\n",
" Maximum neutron transport energy: 20000000.000000 eV for N15\n",
" Minimum neutron data temperature: 294.000000 K\n",
" Maximum neutron data temperature: 294.000000 K\n",
" Reading tallies XML file...\n",
" Preparing distributed cell instances...\n",
" Writing summary.h5 file...\n",
" Initializing source particles...\n",
"\n",
" ===============> FIXED SOURCE TRANSPORT SIMULATION <===============\n",
"\n",
" Simulating batch 1\n",
" Simulating batch 2\n",
" Simulating batch 3\n",
" Simulating batch 4\n",
" Simulating batch 5\n",
" Simulating batch 6\n",
" Simulating batch 7\n",
" Simulating batch 8\n",
" Simulating batch 9\n",
" Simulating batch 10\n",
" Simulating batch 11\n",
" Simulating batch 12\n",
" Simulating batch 13\n",
" Simulating batch 14\n",
" Simulating batch 15\n",
" Simulating batch 16\n",
" Simulating batch 17\n",
" Simulating batch 18\n",
" Simulating batch 19\n",
" Simulating batch 20\n",
" Simulating batch 21\n",
" Simulating batch 22\n",
" Simulating batch 23\n",
" Simulating batch 24\n",
" Simulating batch 25\n",
" Simulating batch 26\n",
" Simulating batch 27\n",
" Simulating batch 28\n",
" Simulating batch 29\n",
" Simulating batch 30\n",
" Simulating batch 31\n",
" Simulating batch 32\n",
" Simulating batch 33\n",
" Simulating batch 34\n",
" Simulating batch 35\n",
" Simulating batch 36\n",
" Simulating batch 37\n",
" Simulating batch 38\n",
" Simulating batch 39\n",
" Simulating batch 40\n",
" Simulating batch 41\n",
" Simulating batch 42\n",
" Simulating batch 43\n",
" Simulating batch 44\n",
" Simulating batch 45\n",
" Simulating batch 46\n",
" Simulating batch 47\n",
" Simulating batch 48\n",
" Simulating batch 49\n",
" Simulating batch 50\n",
" Simulating batch 51\n",
" Simulating batch 52\n",
" Simulating batch 53\n",
" Simulating batch 54\n",
" Simulating batch 55\n",
" Simulating batch 56\n",
" Simulating batch 57\n",
" Simulating batch 58\n",
" Simulating batch 59\n",
" Simulating batch 60\n",
" Simulating batch 61\n",
" Simulating batch 62\n",
" Simulating batch 63\n",
" Simulating batch 64\n",
" Simulating batch 65\n",
" Simulating batch 66\n",
" Simulating batch 67\n",
" Simulating batch 68\n",
" Simulating batch 69\n",
" Simulating batch 70\n",
" Simulating batch 71\n",
" Simulating batch 72\n",
" Simulating batch 73\n",
" Simulating batch 74\n",
" Simulating batch 75\n",
" Simulating batch 76\n",
" Simulating batch 77\n",
" Simulating batch 78\n",
" Simulating batch 79\n",
" Simulating batch 80\n",
" Simulating batch 81\n",
" Simulating batch 82\n",
" Simulating batch 83\n",
" Simulating batch 84\n",
" Simulating batch 85\n",
" Simulating batch 86\n",
" Simulating batch 87\n",
" Simulating batch 88\n",
" Simulating batch 89\n",
" Simulating batch 90\n",
" Simulating batch 91\n",
" Simulating batch 92\n",
" Simulating batch 93\n",
" Simulating batch 94\n",
" Simulating batch 95\n",
" Simulating batch 96\n",
" Simulating batch 97\n",
" Simulating batch 98\n",
" Simulating batch 99\n",
" Simulating batch 100\n",
" Simulating batch 101\n",
" Simulating batch 102\n",
" Simulating batch 103\n",
" Simulating batch 104\n",
" Simulating batch 105\n",
" Simulating batch 106\n",
" Simulating batch 107\n",
" Simulating batch 108\n",
" Simulating batch 109\n",
" Simulating batch 110\n",
" Simulating batch 111\n",
" Simulating batch 112\n",
" Simulating batch 113\n",
" Simulating batch 114\n",
" Simulating batch 115\n",
" Simulating batch 116\n",
" Simulating batch 117\n",
" Simulating batch 118\n",
" Simulating batch 119\n",
" Simulating batch 120\n",
" Simulating batch 121\n",
" Simulating batch 122\n",
" Simulating batch 123\n",
" Simulating batch 124\n",
" Simulating batch 125\n",
" Simulating batch 126\n",
" Simulating batch 127\n",
" Simulating batch 128\n",
" Simulating batch 129\n",
" Simulating batch 130\n",
" Simulating batch 131\n",
" Simulating batch 132\n",
" Simulating batch 133\n",
" Simulating batch 134\n",
" Simulating batch 135\n",
" Simulating batch 136\n",
" Simulating batch 137\n",
" Simulating batch 138\n",
" Simulating batch 139\n",
" Simulating batch 140\n",
" Simulating batch 141\n",
" Simulating batch 142\n",
" Simulating batch 143\n",
" Simulating batch 144\n",
" Simulating batch 145\n",
" Simulating batch 146\n",
" Simulating batch 147\n",
" Simulating batch 148\n",
" Simulating batch 149\n",
" Simulating batch 150\n",
" Simulating batch 151\n",
" Simulating batch 152\n",
" Simulating batch 153\n",
" Simulating batch 154\n",
" Simulating batch 155\n",
" Simulating batch 156\n",
" Simulating batch 157\n",
" Simulating batch 158\n",
" Simulating batch 159\n",
" Simulating batch 160\n",
" Simulating batch 161\n",
" Simulating batch 162\n",
" Simulating batch 163\n",
" Simulating batch 164\n",
" Simulating batch 165\n",
" Simulating batch 166\n",
" Simulating batch 167\n",
" Simulating batch 168\n",
" Simulating batch 169\n",
" Simulating batch 170\n",
" Simulating batch 171\n",
" Simulating batch 172\n",
" Simulating batch 173\n",
" Simulating batch 174\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
" Simulating batch 175\n",
" Simulating batch 176\n",
" Simulating batch 177\n",
" Simulating batch 178\n",
" Simulating batch 179\n",
" Simulating batch 180\n",
" Simulating batch 181\n",
" Simulating batch 182\n",
" Simulating batch 183\n",
" Simulating batch 184\n",
" Simulating batch 185\n",
" Simulating batch 186\n",
" Simulating batch 187\n",
" Simulating batch 188\n",
" Simulating batch 189\n",
" Simulating batch 190\n",
" Simulating batch 191\n",
" Simulating batch 192\n",
" Simulating batch 193\n",
" Simulating batch 194\n",
" Simulating batch 195\n",
" Simulating batch 196\n",
" Simulating batch 197\n",
" Simulating batch 198\n",
" Simulating batch 199\n",
" Simulating batch 200\n",
" Creating state point statepoint.200.h5...\n",
" Writing unstructured mesh tally_1.200.vtk...\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 4.1748e+01 seconds\n",
" Reading cross sections = 2.1717e+00 seconds\n",
" Total time in simulation = 1.8229e+02 seconds\n",
" Time in transport only = 1.0991e+02 seconds\n",
" Time in active batches = 1.8229e+02 seconds\n",
" Time sampling source = 7.1488e+01 seconds\n",
" Time accumulating tallies = 4.7612e-02 seconds\n",
" Total time for finalization = 1.4539e-01 seconds\n",
" Total time elapsed = 2.2446e+02 seconds\n",
" Calculation Rate (active) = 5485.91 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" Leakage Fraction = 0.97394 +/- 0.00017\n",
"\n"
]
}
],
"outputs": [],
"source": [
"openmc.run()"
"openmc.run(output=False)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Again we should see that `tally_1.100.vtk` file which we can use to visualize our results in VisIt or another tool of your choice that supports VTK files."
"Again we should see that `tally_1.200.vtk` file which we can use to visualize our results in VisIt or another tool of your choice that supports VTK files."
]
},
{
@ -716,7 +384,7 @@
"name": "stdout",
"output_type": "stream",
"text": [
"tally_1.10.vtk\ttally_1.200.vtk\r\n"
"tally_1.200.vtk\r\n"
]
}
],
@ -820,326 +488,9 @@
"cell_type": "code",
"execution_count": 15,
"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-2020 MIT and OpenMC contributors\n",
" License | http://openmc.readthedocs.io/en/latest/license.html\n",
" Version | 0.12.0-dev\n",
" Git SHA1 | 34da5a1e2a3942428dfca6e8dfc2ce02d04333dc\n",
" Date/Time | 2020-03-19 10:18:04\n",
" OpenMP Threads | 8\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
" Reading materials XML file...\n",
" Reading DAGMC geometry...\n",
"Loading file dagmc.h5m\n",
"Initializing the GeomQueryTool...\n",
"Using faceting tolerance: 0.001\n",
"Building OBB Tree...\n",
" Reading N14 from /home/shriwise/opt/openmc/xs/nndc_hdf5/N14.h5\n",
" Reading N15 from /home/shriwise/opt/openmc/xs/nndc_hdf5/N15.h5\n",
" Reading O16 from /home/shriwise/opt/openmc/xs/nndc_hdf5/O16.h5\n",
" Reading O17 from /home/shriwise/opt/openmc/xs/nndc_hdf5/O17.h5\n",
" Reading Ar36 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ar36.h5\n",
" WARNING: Negative value(s) found on probability table for nuclide Ar36 at 294K\n",
" Reading Ar38 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ar38.h5\n",
" Reading Ar40 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ar40.h5\n",
" Reading Si28 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si28.h5\n",
" Reading Si29 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si29.h5\n",
" Reading Si30 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Si30.h5\n",
" Reading P31 from /home/shriwise/opt/openmc/xs/nndc_hdf5/P31.h5\n",
" Reading S32 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S32.h5\n",
" Reading S33 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S33.h5\n",
" Reading S34 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S34.h5\n",
" Reading S36 from /home/shriwise/opt/openmc/xs/nndc_hdf5/S36.h5\n",
" Reading Mn55 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mn55.h5\n",
" Reading Fe54 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe54.h5\n",
" Reading Fe56 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe56.h5\n",
" Reading Fe57 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe57.h5\n",
" Reading Fe58 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Fe58.h5\n",
" Reading Ni58 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni58.h5\n",
" Reading Ni60 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni60.h5\n",
" Reading Ni61 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni61.h5\n",
" Reading Ni62 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni62.h5\n",
" Reading Ni64 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Ni64.h5\n",
" Reading Mo100 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo100.h5\n",
" Reading Mo92 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo92.h5\n",
" Reading Mo94 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo94.h5\n",
" Reading Mo95 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo95.h5\n",
" Reading Mo96 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo96.h5\n",
" Reading Mo97 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo97.h5\n",
" Reading Mo98 from /home/shriwise/opt/openmc/xs/nndc_hdf5/Mo98.h5\n",
" Maximum neutron transport energy: 20000000.000000 eV for N15\n",
" Minimum neutron data temperature: 294.000000 K\n",
" Maximum neutron data temperature: 294.000000 K\n",
" Reading tallies XML file...\n",
" Preparing distributed cell instances...\n",
" Writing summary.h5 file...\n",
" Initializing source particles...\n",
"\n",
" ===============> FIXED SOURCE TRANSPORT SIMULATION <===============\n",
"\n",
" Simulating batch 1\n",
" Simulating batch 2\n",
" Simulating batch 3\n",
" Simulating batch 4\n",
" Simulating batch 5\n",
" Simulating batch 6\n",
" Simulating batch 7\n",
" Simulating batch 8\n",
" Simulating batch 9\n",
" Simulating batch 10\n",
" Simulating batch 11\n",
" Simulating batch 12\n",
" Simulating batch 13\n",
" Simulating batch 14\n",
" Simulating batch 15\n",
" Simulating batch 16\n",
" Simulating batch 17\n",
" Simulating batch 18\n",
" Simulating batch 19\n",
" Simulating batch 20\n",
" Simulating batch 21\n",
" Simulating batch 22\n",
" Simulating batch 23\n",
" Simulating batch 24\n",
" Simulating batch 25\n",
" Simulating batch 26\n",
" Simulating batch 27\n",
" Simulating batch 28\n",
" Simulating batch 29\n",
" Simulating batch 30\n",
" Simulating batch 31\n",
" Simulating batch 32\n",
" Simulating batch 33\n",
" Simulating batch 34\n",
" Simulating batch 35\n",
" Simulating batch 36\n",
" Simulating batch 37\n",
" Simulating batch 38\n",
" Simulating batch 39\n",
" Simulating batch 40\n",
" Simulating batch 41\n",
" Simulating batch 42\n",
" Simulating batch 43\n",
" Simulating batch 44\n",
" Simulating batch 45\n",
" Simulating batch 46\n",
" Simulating batch 47\n",
" Simulating batch 48\n",
" Simulating batch 49\n",
" Simulating batch 50\n",
" Simulating batch 51\n",
" Simulating batch 52\n",
" Simulating batch 53\n",
" Simulating batch 54\n",
" Simulating batch 55\n",
" Simulating batch 56\n",
" Simulating batch 57\n",
" Simulating batch 58\n",
" Simulating batch 59\n",
" Simulating batch 60\n",
" Simulating batch 61\n",
" Simulating batch 62\n",
" Simulating batch 63\n",
" Simulating batch 64\n",
" Simulating batch 65\n",
" Simulating batch 66\n",
" Simulating batch 67\n",
" Simulating batch 68\n",
" Simulating batch 69\n",
" Simulating batch 70\n",
" Simulating batch 71\n",
" Simulating batch 72\n",
" Simulating batch 73\n",
" Simulating batch 74\n",
" Simulating batch 75\n",
" Simulating batch 76\n",
" Simulating batch 77\n",
" Simulating batch 78\n",
" Simulating batch 79\n",
" Simulating batch 80\n",
" Simulating batch 81\n",
" Simulating batch 82\n",
" Simulating batch 83\n",
" Simulating batch 84\n",
" Simulating batch 85\n",
" Simulating batch 86\n",
" Simulating batch 87\n",
" Simulating batch 88\n",
" Simulating batch 89\n",
" Simulating batch 90\n",
" Simulating batch 91\n",
" Simulating batch 92\n",
" Simulating batch 93\n",
" Simulating batch 94\n",
" Simulating batch 95\n",
" Simulating batch 96\n",
" Simulating batch 97\n",
" Simulating batch 98\n",
" Simulating batch 99\n",
" Simulating batch 100\n",
" Simulating batch 101\n",
" Simulating batch 102\n",
" Simulating batch 103\n",
" Simulating batch 104\n",
" Simulating batch 105\n",
" Simulating batch 106\n",
" Simulating batch 107\n",
" Simulating batch 108\n",
" Simulating batch 109\n",
" Simulating batch 110\n",
" Simulating batch 111\n",
" Simulating batch 112\n",
" Simulating batch 113\n",
" Simulating batch 114\n",
" Simulating batch 115\n",
" Simulating batch 116\n",
" Simulating batch 117\n",
" Simulating batch 118\n",
" Simulating batch 119\n",
" Simulating batch 120\n",
" Simulating batch 121\n",
" Simulating batch 122\n",
" Simulating batch 123\n",
" Simulating batch 124\n",
" Simulating batch 125\n",
" Simulating batch 126\n",
" Simulating batch 127\n",
" Simulating batch 128\n",
" Simulating batch 129\n",
" Simulating batch 130\n",
" Simulating batch 131\n",
" Simulating batch 132\n",
" Simulating batch 133\n",
" Simulating batch 134\n",
" Simulating batch 135\n",
" Simulating batch 136\n",
" Simulating batch 137\n",
" Simulating batch 138\n",
" Simulating batch 139\n",
" Simulating batch 140\n",
" Simulating batch 141\n",
" Simulating batch 142\n",
" Simulating batch 143\n",
" Simulating batch 144\n",
" Simulating batch 145\n",
" Simulating batch 146\n",
" Simulating batch 147\n",
" Simulating batch 148\n",
" Simulating batch 149\n",
" Simulating batch 150\n",
" Simulating batch 151\n",
" Simulating batch 152\n",
" Simulating batch 153\n",
" Simulating batch 154\n",
" Simulating batch 155\n",
" Simulating batch 156\n",
" Simulating batch 157\n",
" Simulating batch 158\n",
" Simulating batch 159\n",
" Simulating batch 160\n",
" Simulating batch 161\n",
" Simulating batch 162\n",
" Simulating batch 163\n",
" Simulating batch 164\n",
" Simulating batch 165\n",
" Simulating batch 166\n",
" Simulating batch 167\n",
" Simulating batch 168\n",
" Simulating batch 169\n",
" Simulating batch 170\n",
" Simulating batch 171\n",
" Simulating batch 172\n",
" Simulating batch 173\n",
" Simulating batch 174\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
" Simulating batch 175\n",
" Simulating batch 176\n",
" Simulating batch 177\n",
" Simulating batch 178\n",
" Simulating batch 179\n",
" Simulating batch 180\n",
" Simulating batch 181\n",
" Simulating batch 182\n",
" Simulating batch 183\n",
" Simulating batch 184\n",
" Simulating batch 185\n",
" Simulating batch 186\n",
" Simulating batch 187\n",
" Simulating batch 188\n",
" Simulating batch 189\n",
" Simulating batch 190\n",
" Simulating batch 191\n",
" Simulating batch 192\n",
" Simulating batch 193\n",
" Simulating batch 194\n",
" Simulating batch 195\n",
" Simulating batch 196\n",
" Simulating batch 197\n",
" Simulating batch 198\n",
" Simulating batch 199\n",
" Simulating batch 200\n",
" Creating state point statepoint.200.h5...\n",
" WARNING: Skipping unstructured mesh writing for tally 1. More than one filter\n",
" is present on the tally.\n",
"\n",
" =======================> TIMING STATISTICS <=======================\n",
"\n",
" Total time for initialization = 4.0432e+01 seconds\n",
" Reading cross sections = 1.9152e+00 seconds\n",
" Total time in simulation = 1.7878e+02 seconds\n",
" Time in transport only = 1.0807e+02 seconds\n",
" Time in active batches = 1.7878e+02 seconds\n",
" Time sampling source = 7.0570e+01 seconds\n",
" Time accumulating tallies = 9.5415e-02 seconds\n",
" Total time for finalization = 2.9952e-01 seconds\n",
" Total time elapsed = 2.1979e+02 seconds\n",
" Calculation Rate (active) = 5593.58 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
" Leakage Fraction = 0.97394 +/- 0.00017\n",
"\n"
]
}
],
"outputs": [],
"source": [
"openmc.run()"
"openmc.run(output=False)"
]
},
{
@ -1273,7 +624,7 @@
"name": "stdout",
"output_type": "stream",
"text": [
"manifold_flux.vtk tally_1.10.vtk tally_1.200.vtk\r\n"
"manifold_flux.vtk tally_1.200.vtk\r\n"
]
}
],