def test_stagnation_pressure():
"""
This function tests the stagnation_pressure() function
"""
static_pressure = 6000000.0
stag_pres = Stagnation.stagnation_pressure(static_pressure, mach_number, gamma)
assert isclose(stag_pres, 6095032.0, rel_tol=1.0e-3)
def test_reactor_loop():
elec_power = 50000.0
species = 'nitrogen'
gas = Chemical(species)
heat_eff = 0.90
turb_eff = 0.90
comp_rat_guess = 1.18
comp_eff = 0.90
conv_eff = 0.25
inlet_stat_temp = 800.0 # Kelvins
inlet_stat_pres = 6000000.0 # Pascals
gas.calculate(T=inlet_stat_temp, P=inlet_stat_pres)
mw = gas.MW
gamma1 = gas.Cp / gas.Cvg
cp = gas.Cp
inlet_mach_number = 0.1
inlet_stag_temp = Stagnation.stagnation_temperature(inlet_stat_temp,
inlet_mach_number, gamma1)
inlet_stag_pres = Stagnation.stagnation_pressure(inlet_stat_pres,
inlet_mach_number, gamma1)
deltat = 50.0 # Kelvins
mdot = (elec_power / conv_eff) / (cp * deltat)
reac = ReactorLoop(heat_eff, turb_eff, species)
dicts = reac.performance(inlet_stag_temp, inlet_stag_pres, inlet_mach_number, inlet_stat_temp,
inlet_stat_pres, elec_power, comp_rat_guess, comp_eff,
conv_eff, mdot)
print('Compressor')
print(dicts[0])
print('Heat')
print(dicts[1])
print('Turbine')
print(dicts[2])
print('Rejection')
print(dicts[3])