crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-26 21:25:36 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

Updating unstructured mesh examples to include the library argument in Python classes.

Browse source
This commit is contained in:
Patrick Shriwise 2021-01-06 10:25:56 -06:00
parent 6c7c0a6816
commit e9a13d67ab
2 changed files with 61 additions and 68 deletions
Download patch file Download diff file Expand all files Collapse all files

62
examples/jupyter/unstructured-mesh-part-i.ipynb
View file

File diff suppressed because one or more lines are too long

67
examples/jupyter/unstructured-mesh-part-ii.ipynb
View file

@ -244,20 +244,24 @@
"\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 | c9cbdb7c70b202e847c7169f9e5602f110a43853\n",
" Date/Time | 2020-04-24 09:25:02\n",
" OpenMP Threads | 8\n",
" License | https://docs.openmc.org/en/latest/license.html\n",
" Version | 0.12.1-dev\n",
" Git SHA1 | 8ae407c90e927af62bfdc8f150e96bfd5d7b2ec2\n",
" Date/Time | 2021-01-06 09:17:05\n",
" MPI Processes | 1\n",
" OpenMP Threads | 2\n",
"\n",
" Reading settings XML file...\n",
" Reading cross sections XML file...\n",
" Reading materials XML file...\n",
" Reading DAGMC geometry...\n",
"Set overlap thickness = 0\n",
"Set numerical precision = 0.001\n",
"Loading file dagmc.h5m\n",
"Initializing the GeomQueryTool...\n",
"Using faceting tolerance: 0.001\n",
"Building OBB Tree...\n",
"Building acceleration data structures...\n",
"Implicit Complement assumed to be Vacuum\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",
@ -291,12 +295,11 @@
" 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",
" Minimum neutron data temperature: 294.0 K\n",
" Maximum neutron data temperature: 294.0 K\n",
" Preparing distributed cell instances...\n",
" Writing summary.h5 file...\n",
" Maximum neutron transport energy: 20000000.0 eV for N15\n",
"\n",
" ===============> FIXED SOURCE TRANSPORT SIMULATION <===============\n",
"\n",
@ -311,20 +314,19 @@
" Simulating batch 9\n",
" Simulating batch 10\n",
" Creating state point statepoint.10.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.3651e+01 seconds\n",
" Reading cross sections = 2.1907e+00 seconds\n",
" Total time in simulation = 3.4743e-01 seconds\n",
" Time in transport only = 2.9555e-01 seconds\n",
" Time in active batches = 3.4743e-01 seconds\n",
" Time accumulating tallies = 7.6575e-03 seconds\n",
" Total time for finalization = 2.2273e-01 seconds\n",
" Total time elapsed = 4.4224e+01 seconds\n",
" Calculation Rate (active) = 2878.27 particles/second\n",
" Total time for initialization = 2.4196e+01 seconds\n",
" Reading cross sections = 2.1428e+00 seconds\n",
" Total time in simulation = 1.9512e-01 seconds\n",
" Time in transport only = 1.9269e-01 seconds\n",
" Time in active batches = 1.9512e-01 seconds\n",
" Time accumulating tallies = 2.0220e-06 seconds\n",
" Time writing statepoints = 2.2461e-03 seconds\n",
" Total time for finalization = 1.8940e-06 seconds\n",
" Total time elapsed = 2.4401e+01 seconds\n",
" Calculation Rate (active) = 5125.02 particles/second\n",
"\n",
" ============================> RESULTS <============================\n",
"\n",
@ -350,7 +352,7 @@
"metadata": {},
"outputs": [],
"source": [
"unstructured_mesh = openmc.UnstructuredMesh(\"manifold.h5m\")\n",
"unstructured_mesh = openmc.UnstructuredMesh(\"manifold.h5m\", library='moab')\n",
"\n",
"mesh_filter = openmc.MeshFilter(unstructured_mesh)\n",
"\n",
@ -399,8 +401,7 @@
"name": "stdout",
"output_type": "stream",
"text": [
"manifold_flux.vtk tally_1.100.vtk tally_1.200.vtk\r\n",
"manifold.vtk\t tally_1.10.vtk\r\n"
"tally_1.200.vtk\r\n"
]
}
],
@ -495,7 +496,7 @@
"text": [
"<?xml version='1.0' encoding='utf-8'?>\r\n",
"<tallies>\r\n",
" <mesh id=\"1\" type=\"unstructured\">\r\n",
" <mesh id=\"1\" library=\"moab\" type=\"unstructured\">\r\n",
" <filename>manifold.h5m</filename>\r\n",
" </mesh>\r\n",
" <filter id=\"1\" type=\"mesh\">\r\n",
@ -562,16 +563,7 @@
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/shriwise/.pyenv/versions/3.7.3/lib/python3.7/site-packages/vtk/util/numpy_support.py:137: FutureWarning: Conversion of the second argument of issubdtype from `complex` to `np.complexfloating` is deprecated. In future, it will be treated as `np.complex128 == np.dtype(complex).type`.\n",
" assert not numpy.issubdtype(z.dtype, complex), \\\n"
]
}
],
"outputs": [],
"source": [
"data_dict = {'Flux 0 - 1 MeV' : thermal_flux,\n",
" 'Flux 1 - 5 MeV' : fast_flux,\n",
@ -596,8 +588,7 @@
"name": "stdout",
"output_type": "stream",
"text": [
"manifold_flux.vtk tally_1.100.vtk tally_1.200.vtk\r\n",
"manifold.vtk\t tally_1.10.vtk\r\n"
"manifold.vtk tally_1.200.vtk\r\n"
]
}
],
@ -647,7 +638,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.3"
"version": "3.7.3"
}
},
"nbformat": 4,