spectrum = "thermal" # two-point model n_ref = 0 Lambda_ref = 0 ref_lambda = [] ref_rho = [] # Feedbacks, False to turn reactivity feedback off. True otherwise. feedback = True # External Reactivity from reactivity_insertion import RampReactivityInsertion rho_ext = RampReactivityInsertion(timer=ti, t_start=60.0 * units.seconds, t_end=70.0 * units.seconds, rho_init=0.0 * units.delta_k, rho_rise=600.0 * units.pcm, rho_final=600.0 * units.pcm) # maximum number of internal steps that the ode solver will take nsteps = 5000 k_mod = random.gauss(17, 17 * 0.05) * units.watt / (units.meter * units.kelvin) cp_mod = random.gauss(1650.0, 1650.0 * 0.05) * \ units.joule / (units.kg * units.kelvin) rho_mod = DensityModel(a=1740. * units.kg / (units.meter**3), model="constant") Moderator = Material('mod', k_mod, cp_mod, dm=rho_mod) k_fuel = random.uniform(15.0, 19.0) * units.watt / (units.meter * units.kelvin) cp_fuel = random.gauss( 1818.0, 1818 * 0.05) * units.joule / units.kg / units.kelvin # [J/kg/K]
ref_rho = [0.084349, 0.168983] # Feedbacks, False to turn reactivity feedback off. True otherwise. feedback = True # External Reactivity from reactivity_insertion import RampReactivityInsertion # from reactivity_insertion import StepReactivityInsertion # rho_ext = StepReactivityInsertion(timer=ti, # t_step=t_feedback + 10.0*units.seconds, # rho_init=0.0*units.delta_k, # rho_final=600.0*units.pcm) rho_ext = RampReactivityInsertion(timer=ti, t_start=t_feedback + 10.0 * units.seconds, t_end=t_feedback + 15.0 * units.seconds, rho_init=0.0 * units.delta_k, rho_rise=650.0 * units.pcm, rho_final=650.0 * units.pcm) # maximum number of internal steps that the ode solver will take nsteps = 5000 k_mod = 17 * units.watt / (units.meter * units.kelvin) cp_mod = 1650.0 * units.joule / (units.kg * units.kelvin) rho_mod = DensityModel(a=1740. * units.kg / (units.meter**3), model="constant") Moderator = Material('mod', k_mod, cp_mod, rho_mod) k_fuel = 15 * units.watt / (units.meter * units.kelvin) cp_fuel = 1818.0 * units.joule / units.kg / units.kelvin rho_fuel = DensityModel(a=2220.0 * units.kg / (units.meter**3), model="constant") Fuel = Material('fuel', k_fuel, cp_fuel, rho_fuel)