class Transition_L2G():
t_max = 10.0
initial = dict(T=350, p=5.0e5, phase="Liquid")
params = dict(k_p=1.0e-2, k_T=1.0e7, Dt=t_max / 100.0)
@staticmethod
def schedule(sim, t):
sim.set_params(T_inf=350, p_inf=5.0e3)
answer = dict(T=350,
p=5.0e3,
rho=iapws97.density_region2(350.0, 5.0e3),
h=iapws97.enthalpy_region2(350.0, 5.0e3))
class Hot_Gas_Fill():
t_max = 1000.0
initial = dict(T=450, p=5.0e5, phase="Gas")
params = dict(mass_source=0.1, k_p=0.0, k_T=0.0, Dt=t_max / 100.0)
@staticmethod
def schedule(sim, t):
pass
answer = dict(T=450,
p=5.0e5,
rho=iapws97.density_region2(450.0, 5.0e5),
h=iapws97.enthalpy_region2(450.0, 5.0e5))
class Hot_Gas():
t_max = 1000.0
initial = dict(T=450, p=5.0e5, phase="Gas")
params = dict(k_p=1.0e-4, k_T=1.0e4, Dt=100.0)
@staticmethod
def schedule(sim, t):
sim.set_params(T_inf=550, p_inf=5.0e5)
answer = dict(T=550,
p=5.0e5,
rho=iapws97.density_region2(550.0, 5.0e5),
h=iapws97.enthalpy_region2(550.0, 5.0e5))
class Liquid_Drain():
t_max = 2000.0
initial = dict(T=373.15, p=3.0e5, phase="Liquid")
params = dict(mass_source=-0.1, k_p=0.0, k_T=0.0, Dt=t_max / 100.0)
@staticmethod
def schedule(sim, t):
pass
answer = dict(T=350,
p=5.0e3,
rho=iapws97.density_region2(350.0, 5.0e3),
h=iapws97.enthalpy_region2(350.0, 5.0e3))
class Cycle_sgclg():
t_max = 100.0
initial = dict(T=250, p=5.0e3, phase="Solid")
params = dict(k_p=1.0e-4, k_T=1.0e3)
@staticmethod
def schedule(sim, t):
if t < 1000.0:
sim.set_params(T_inf=800, p_inf=5.0e4)
elif t < 2000.0:
sim.set_params(T_inf=800, p_inf=3.0e7)
elif t < 3000.0:
sim.set_params(T_inf=400, p_inf=3.0e7)
elif t < 4000.0:
sim.set_params(T_inf=400, p_inf=5.0e3)
answer = dict(T=400,
p=5.0e3,
rho=iapws97.density_region2(400.0, 5.0e3),
h=iapws97.enthalpy_region2(400.0, 5.0e3))