openmc.deplete.abc.SIIntegrator¶
-
class
openmc.deplete.abc.
SIIntegrator
(operator, timesteps, power=None, power_density=None, n_steps=10)[source]¶ Abstract class for the Stochastic Implicit Euler integrators
Does not provide a
__call__
method, but scales and resets the number of particles used in initial transport calculationParameters: - operator (openmc.deplete.TransportOperator) – The operator object to simulate on.
- timesteps (iterable of float) – Array of timesteps in units of [s]. Note that values are not cumulative.
- power (float or iterable of float, optional) – Power of the reactor in [W]. A single value indicates that
the power is constant over all timesteps. An iterable
indicates potentially different power levels for each timestep.
For a 2D problem, the power can be given in [W/cm] as long
as the “volume” assigned to a depletion material is actually
an area in [cm^2]. Either
power
orpower_density
must be specified. - power_density (float or iterable of float, optional) – Power density of the reactor in [W/gHM]. It is multiplied by
initial heavy metal inventory to get total power if
power
is not speficied. - n_steps (int, optional) – Number of stochastic iterations per depletion interval. Must be greater than zero. Default : 10
Variables: - operator (openmc.deplete.TransportOperator) – Operator to perform transport simulations
- chain (openmc.deplete.Chain) – Depletion chain
- timesteps (iterable of float) – Size of each depletion interval in [s]
- power (iterable of float) – Power of the reactor in [W] for each interval in
timesteps
- n_steps (int) – Number of stochastic iterations per depletion interval
-
__call__
(conc, rates, dt, power, i)¶ Perform the integration across one time step
Parameters: - conc (numpy.ndarray) – Initial concentrations for all nuclides in [atom]
- rates (openmc.deplete.ReactionRates) – Reaction rates from operator
- dt (float) – Time in [s] for the entire depletion interval
- power (float) – Power of the system in [W]
- i (int) – Current depletion step index
Returns: - proc_time (float) – Time spent in CRAM routines for all materials in [s]
- conc_list (list of numpy.ndarray) – Concentrations at each of the intermediate points with the final concentration as the last element
- op_results (list of openmc.deplete.OperatorResult) – Eigenvalue and reaction rates from intermediate transport simulations
-
__iter__
()¶ Return pairs of time steps in [s] and powers in [W]
-
__len__
()¶ Return integer number of depletion intervals