Skip to main content

Convert OpenMC tallies into user friendly units

Project description

This Python package aims to help convert OpenMC tallies to user specified units.

Installation

pip install openmc_tally_unit_converter

Usage

OpenMC tally results are save into a statepoint h5 file without units.

This package ascertains the base units of common tallies by inspecting their tally filters and scores.

The following worked example is for a heating tally. Flux, effective dose and DPA / damage-energy tallies are also supported.

import openmc_tally_unit_converter as otuc
import openmc

# loads up tally from an openmc statepoint output file
statepoint = openmc.StatePoint(filepath="statepoint.2.h5")
my_tally = statepoint.get_tally(name="my_cell_heating_tally")

# gets the base units of the tally
tally = otuc.process_tally(tally=my_tally)
print(tally)
>>> 218559.22320927 electron_volt / source_particle

The package then allows users to scale the base tally units to different units. For example the package can easily convert cm to meters or electron volts to joules.

converted_tally = otuc.process_tally(
    tally=my_tally,
    required_units = "joules / source_particle"
)

print(converted_tally)
>>> 3.50170481e-14 Joules / source_particle

Additional inputs can be applied to normalize the the tallies and convert the units further:

  • The source strength of the source in particles per second can be specified with the strength_strength argument. This allows the tally results to be converted from the units of score per simulated particle to score per unit time (e.g seconds, hours, days etc).
converted_tally = otuc.process_tally(
    tally=my_tally,
    source_strength=1e20,  # input units for this argument are particles per second
    required_units = "joules / minute"
)

print(converted_tally)
>>> 2.10102288e+08 joules / source_particle
  • The volume of the cell in cm3 can also be specified with the volume argument. This allows the tally result to be converted from the base units to base units per unit volume.
converted_tally = otuc.process_tally(
    tally=my_tally,
    source_strength=13458.3,  # input units for this argument are particles per second
    volume=12,  # input units are in cm3
    required_units = "joules / second / meter **3"
)

print(converted_tally)
>>> 3.92724948e-05 Joules / meter ** 3 / second

:point_right: Further examples

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for openmc-tally-unit-converter, version 0.3.0
Filename, size File type Python version Upload date Hashes
Filename, size openmc_tally_unit_converter-0.3.0.tar.gz (8.1 kB) File type Source Python version None Upload date Hashes View
Filename, size openmc_tally_unit_converter-0.3.0-py3-none-any.whl (8.0 kB) File type Wheel Python version py3 Upload date Hashes View

Supported by

AWS AWS Cloud computing Datadog Datadog Monitoring Facebook / Instagram Facebook / Instagram PSF Sponsor Fastly Fastly CDN Google Google Object Storage and Download Analytics Huawei Huawei PSF Sponsor Microsoft Microsoft PSF Sponsor NVIDIA NVIDIA PSF Sponsor Pingdom Pingdom Monitoring Salesforce Salesforce PSF Sponsor Sentry Sentry Error logging StatusPage StatusPage Status page