def test_Anderson_example(self):
# Example 10.5 from Anderson2016 (recalculated manually due to different formula, calculation approach)
expected_result = 586.1 # [kg/s]
res = IRT.mass_flow(p_chamber=3.03e6,
A_throat=0.4,
R=520,
T_chamber=3500,
gamma=1.22)
self.assertEqual(round(res, 1), expected_result)
def test_zero_two(self):
expected_result = 0.769800358919501
res = IRT.mass_flow(p_chamber=1, A_throat=1, R=1, T_chamber=1, gamma=2)
self.assertEqual(res, expected_result)