# Number of decay heat groups
n_dg = 0
# Fissioning Isotope
fission_iso = "sfr"
# Spectrum
spectrum = "fast"
# False to turn reactivity feedback off.
feedback = False
# External Reactivity
from reactivity_insertion import ReactivityInsertion
rho_ext = ReactivityInsertion(timer=ti)
# rho_ext = StepReactivityInsertion(timer=ti, t_step=1.0*units.seconds,
# rho_init=0.0*units.delta_k,
# rho_final=0.005*units.delta_k)
fuel = th.THComponent(name="fuel",
mat=SFRMetal(name="sfrfuel"),
vol=vol_fuel,
T0=t_fuel,
alpha_temp=alpha_f,
timer=ti,
heatgen=True,
power_tot=power_tot)
cool = th.THComponent(name="cool",
def init_rho_ext(self, rho_ext):
if rho_ext is None:
rho_ext = ReactivityInsertion(self._timer)
return rho_ext