Repository navigation
Add reactivity estimation for continuous load following #81
Description
Current behavior:
pyrk only accepts a reactivity insertion and returns a power response.
Desired behavior:
pyrk should accept
- initial or current power level
- target power level
- estimation method (e.g. prompt jump approximation)
and initialize the reactivity insertion based on the estimation.
How to close:
This issue can be closed when neutronics.py has the new parameters:
-
p0 -
pt -
rho_extaccepts a string"pja"for the prompt jump approximation
class Neutronics(object):
"""This class handles calculations and data related to the
neutronics subblock
"""
def __init__(self, iso="u235", e="thermal", n_precursors=6, n_decay=11,
n_fic=0,
timer=Timer(),
rho_ext=None,
feedback=False,
p0=None,
pt =None):
"""
:param p0: The starting power of the reactor.
:param pt: The target power of the reactor.
"""- and there is a function that estimates the reactivity via the requested method.
Guidance
- Short post about the prompt jump approximation.
- Recall that reactor power is proportional neutron population (P=fission_xs * neutron_density * flux * volume * chi)
Activity