EEOUT files (Unstructured Meshes)

The complete API can be found at f4enix.output.eeout.EEOUT

F4Enix can also parse outputs of unstrutured meshes MCNP simulation (i.e., eeout files) and export it to .vtu format.

The feature is limited for the moment to 1st and 2nd order tetraheadra with no mixed formulation. The removing of this limitations though would not take much effort and can be considered in the future if deemed necessary.

from f4enix.output.eeout import EEOUT
from pprint import pprint  # just to print dictionaries in a nicer way

# Parse the eeout file
eeout = EEOUT('cyl_tetra1.eeout')
pprint(eeout.info)

/home/docs/checkouts/readthedocs.org/user_builds/f4enix/envs/stable/lib/python3.10/site-packages/tqdm/auto.py:21: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html
  from .autonotebook import tqdm as notebook_tqdm

{'NUMBER OF 1st HEXS': 0,
 'NUMBER OF 1st PENTS': 0,
 'NUMBER OF 1st TETS': 1306,
 'NUMBER OF 2nd HEXS': 0,
 'NUMBER OF 2nd PENTS': 0,
 'NUMBER OF 2nd TETS': 0,
 'NUMBER OF EDITS': 2,
 'NUMBER OF HISTORIES': 1000000,
 'NUMBER OF INSTANCES': 2,
 'NUMBER OF MATERIALS': 2,
 'NUMBER OF NODES': 391,
 'NUMBER OF PARTICLES': 1}

# A PyVista object is created containing all available data in the eeout file
eeout.grid

UnstructuredGrid (0x70608579fdc0)
  N Cells:    1306
  N Points:   391
  X Bounds:   -1.000e+01, 1.000e+01
  Y Bounds:   -1.000e+01, 1.000e+01
  Z Bounds:   0.000e+00, 5.000e+01
  N Arrays:   7

UnstructuredGrid 391 points 1,306 cells 180 KiB

Bounds⧉X [-1.000e+01, 1.000e+01]Y [-1.000e+01, 1.000e+01]Z [0.000e+00, 5.000e+01]

• Cell Data: (7)
  density active
  • material⧉
    scalar
    <U9
  • density⧉
    scalar
    float64
    [9.450e-01, 7.300e+00]
    active
  • volume⧉
    scalar
    float64
    [9.480e-01, 3.142e+01]
  • Tally1_par1_FLUX_14 - value⧉
    scalar
    float64
    [0.000e+00, 1.658e-01]
  • Tally1_par1_FLUX_14 - error⧉
    scalar
    float64
    [1.492e-03, 1.000e+00]
  • Tally2_par1_ENERGY_16 - value⧉
    scalar
    float64
    [0.000e+00, 3.240e-02]
  • Tally2_par1_ENERGY_16 - error⧉
    scalar
    float64
    [2.007e-03, 1.000e+00]

# import pyvista as pv

# # As usual it is possible to directly work with the PyVista object
# plotter = pv.Plotter()
# plotter.add_mesh(eeout.grid.clip(), scalars='Tally2_par1_ENERGY_16 - value',
#                  log_scale=True, cmap='jet', show_edges=True,
#                  scalar_bar_args={'vertical': True})
# plotter.show(jupyter_backend='static')

# or export it to .vtu to use it for instance in ParaView
import tempfile  # only for tmp directory

outfolder = tempfile.gettempdir()
eeout.export(outfolder)