def make_flfacts_alt_sweep(self,
model: BDF,
mach,
alts,
eas_limit: float = 1000.,
alt_units: str = 'm',
velocity_units: str = 'm/s',
density_units: str = 'kg/m^3',
eas_units: str = 'm/s') -> Tuple[Any, Any, Any]:
"""makes an altitude sweep"""
alts.sort()
alts = alts[::-1]
rho, mach, velocity = make_flfacts_alt_sweep(
mach,
alts,
eas_limit=eas_limit,
alt_units=alt_units,
velocity_units=velocity_units,
density_units=density_units,
eas_units=eas_units)
flfact_rho = self.sid + 1
flfact_mach = self.sid + 2
flfact_velocity = self.sid + 3
flfact_eas = self.sid + 4
flfact_alt = self.sid + 5
alts2 = alts[:len(rho)]
assert len(rho) == len(alts2)
comment = ' density: min=%.3e max=%.3e %s; alt min=%.0f max=%.0f %s' % (
rho.min(),
rho.max(),
density_units,
alts2.min(),
alts2.max(),
alt_units,
)
model.add_flfact(flfact_rho, rho, comment=comment)
model.add_flfact(flfact_mach, mach, comment=' Mach: %s' % mach.min())
comment = ' velocity: min=%.3f max=%.3f %s' % (
velocity.min(), velocity.max(), velocity_units)
model.add_flfact(flfact_velocity, velocity, comment=comment)
# eas in velocity units
rho0 = atm_density(0.,
alt_units=alt_units,
density_units=density_units)
eas = velocity * np.sqrt(rho / rho0)
kvel = _velocity_factor(velocity_units, eas_units)
eas_in_eas_units = eas * kvel
comment = ' EAS: min=%.3f max=%.3f %s' % (
eas_in_eas_units.min(), eas_in_eas_units.max(), eas_units)
model.add_flfact(flfact_eas, eas_in_eas_units, comment=comment)
comment = ' Alt: min=%.3f max=%.3f %s' % (alts2.min(), alts2.max(),
alt_units)
model.add_flfact(flfact_alt, alts2, comment=comment)
def test_sweep(self):
"""tests FLFACT sweeps"""
mach = 0.8
alts = np.linspace(-10000, 80000.)
rho, mach, vel = make_flfacts_alt_sweep(mach,
alts,
eas_limit=300.,
alt_units='m',
velocity_units='m/s',
density_units='kg/m^3',
eas_units='m/s')
del rho, mach, vel, alts
alt = 10000.
machs = np.linspace(0., 0.8)
rho, mach, vel = make_flfacts_mach_sweep(alt,
machs,
eas_limit=300.,
alt_units='m',
velocity_units='m/s',
density_units='kg/m^3',
eas_units='m/s')
del rho, mach, vel, alt
alt = 0.
machs = np.linspace(0.6, 0.8)
with self.assertRaises(RuntimeError):
rho, mach, vel = make_flfacts_mach_sweep(alt,
machs,
eas_limit=100.,
alt_units='m',
velocity_units='m/s',
density_units='kg/m^3',
eas_units='m/s')
del alt
alt = 10000.
eass = np.linspace(0., 300.)
rho, mach, vel = make_flfacts_eas_sweep(alt,
eass,
alt_units='m',
velocity_units='m/s',
density_units='kg/m^3',
eas_units='m/s')
del rho, mach, vel, alt
def test_sweep(self):
"""tests FLFACT sweeps"""
mach = 0.8
alts = np.linspace(-10000, 80000.)
rho, mach, vel = make_flfacts_alt_sweep(
mach, alts, eas_limit=300., alt_units='m',
velocity_units='m/s',
density_units='kg/m^3',
eas_units='m/s')
del rho, mach, vel, alts
alt = 10000.
machs = np.linspace(0., 0.8)
rho, mach, vel = make_flfacts_mach_sweep(
alt, machs, eas_limit=300.,
alt_units='m',
velocity_units='m/s',
density_units='kg/m^3',
eas_units='m/s')
del rho, mach, vel, alt