def test_02(self):
# Truth values from NASA RP 1046
Value = SA.alt_temp2density_ratio(8000, -0.8496)
Truth = 0.06011 / 0.076474
self.failUnless(RE(Value, Truth) <= 1e-5)
def test_02(self):
# Truth values from NASA RP 1046
Value = SA.alt_temp2density_ratio(8000, -0.8496)
Truth = 0.06011 / 0.076474
self.assertTrue(RE(Value, Truth) <= 1e-5)
def test_03(self):
# truth value calculated from program at:
# http://www.sworld.com.au/steven/space/atmosphere/
Value = SA.alt_temp2density_ratio(1999, 275.1565, alt_units="m", temp_units="K")
Truth = 0.82168
self.failUnless(RE(Value, Truth) <= 5e-5)
def test_03(self):
# truth value calculated from program at:
# http://www.sworld.com.au/steven/space/atmosphere/
Value = SA.alt_temp2density_ratio(1999,
275.1565,
alt_units='m',
temp_units='K')
Truth = .82168
self.assertTrue(RE(Value, Truth) <= 5e-5)
def power_alt(ff, n, Hp, OAT, k=.00283, k2=.007, ff_units='USG/h', temp_units = 'F'):
"""
Alternate, simplified method to determine power. Assume that the fuel
flow for peak EGT varies linearly with density ratio and rpm. Then use
the main power method to calculate power for cases where there is not a
good indication of fuel flow at peak EGT. Only applicable to Lycoming
IO-360A series engines. Based on a very limited set of flight test data,
so the accuracy of the fuel flow at peak EGT has not been established
under all conditions.
"""
oat = U.temp_conv(OAT, from_units='F', to_units='R')
dr = SA.alt_temp2density_ratio(Hp, OAT, temp_units='F')
ff_peak_EGT = k * dr * n + k2 * oat
p = power(ff, ff_peak_EGT, n)
return p
def power_alt(ff,
n,
Hp,
OAT,
k=.00283,
k2=.007,
ff_units='USG/h',
temp_units='F'):
"""
Alternate, simplified method to determine power. Assume that the fuel
flow for peak EGT varies linearly with density ratio and rpm. Then use
the main power method to calculate power for cases where there is not a
good indication of fuel flow at peak EGT. Only applicable to Lycoming
IO-360A series engines. Based on a very limited set of flight test data,
so the accuracy of the fuel flow at peak EGT has not been established
under all conditions.
"""
oat = U.temp_conv(OAT, from_units='F', to_units='R')
dr = SA.alt_temp2density_ratio(Hp, OAT, temp_units='F')
ff_peak_EGT = k * dr * n + k2 * oat
p = power(ff, ff_peak_EGT, n)
return p
def test_01(self):
PR = SA.alt_temp2density_ratio(0, 59, temp_units='F')
self.assertEqual(PR, 1)
def test_01(self):
PR = SA.alt_temp2density_ratio(0, 59, temp_units="F")
self.assertEqual(PR, 1)
def ffp(Hp, OAT, n, k=.00283, k2=.007):
oat = U.temp_conv(OAT, from_units='F', to_units='R')
dr = SA.alt_temp2density_ratio(Hp, OAT, temp_units='F')
ff_peak_EGT = k * dr * n + k2 * oat
return ff_peak_EGT
def ffp(Hp, OAT, n, k=.00283, k2=.007):
oat = U.temp_conv(OAT, from_units='F', to_units='R')
dr = SA.alt_temp2density_ratio(Hp, OAT, temp_units='F')
ff_peak_EGT = k * dr * n + k2 * oat
return ff_peak_EGT
