In [1]:

    

import isa
from isentropic_gas import IsentropicGas


    

In [2]:

    

gas_isoentropico_real = IsentropicGas(selected_cp_air_model='naca', selected_gamma_air_model='naca')
gas_isoentropico_ideal = IsentropicGas(selected_cp_air_model='ideal', selected_gamma_air_model='ideal')


    

In [3]:

    

gas_isoentropico_ideal?


    

In [4]:

    

h0 = 3500 # m
temp0 = isa.temp_isa(h0) #K
pres0 = isa.pres_isa(h0)
v0 = 163.65 #m/s
gas_isoentropico_ideal.mach_number?
mach0 = gas_isoentropico_ideal.mach_number(v0, temp0)
print(mach0)


    

    




0.500121411084

In [5]:

    

gas_isoentropico_ideal.stagnation_temperature_from_mach?
T0t = gas_isoentropico_ideal.stagnation_temperature_from_mach(mach0, temp0)
print(T0t)


    

    




278.737312002

In [6]:

    

gas_isoentropico_real.stagnation_pressure_from_mach?
P0t = gas_isoentropico_real.stagnation_pressure_from_mach(mach0, pres0, temp0)
print(P0t)


    

    




78103.5225407

In [7]:

    

gas_isoentropico_ideal.velocity_from_stagnation_temperature?
mach0_ = gas_isoentropico_ideal.velocity_from_stagnation_temperature(T0t, temp0)
v0_ = gas_isoentropico_ideal.mach_from_stagnation_temperature(T0t, temp0)
print(mach0_, v0_)


    

    




163.65 0.500121411084

In [ ]:

    

gas_isoentropico_ideal.velocity_from_stagnation_temperature?
mach0_ = gas_isoentropico_real.velocity_from_stagnation_temperature(T0t, temp0)
v0_ = gas_isoentropico_real.mach_from_stagnation_temperature(T0t, temp0)
print(mach0_, v0_)


    

In [ ]: