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])