Burnup(MWd/kg) vs Time graph

User Support
1

Hello everyone,
I’m new to openMC. So, I do not know most of its functions. Currently, I’m running a depletion analysis on my BWR model. The timestep_unit for my depletion is ‘d’. So, How can I get the burnup(MWd/kg) for the specified days in the depletion to draw a graph for the time in the X axis and burnup(MWd/kg) in the Y axis?
Thanks,
Niloy

1 Like
2

use time_step units MWd/kg

Example:
import numpy as np

Create depletion “operator”

chain_file = ‘chain_endfb71_pwr.xml’
op = openmc.deplete.Operator(geometry, settings_file, chain_file)
burnup_step= np.array([6.0, 12.0, 18.0, 24.0, 30.0, 36.0, 42.0, 48.0, 54.0, 60.0]) #MWd/kg
burnup = np.diff ( burnup_step, prepend =0.0 )
power = 174 # W/cm, for 2D simulations only
integrator = openmc.deplete.PredictorIntegrator(op, burnup, power, timestep_units=‘MWd/kg’)
integrator.integrate()