# Initialization of Variable
from math import pi
from math import atan
from numpy import *
T1=80.0 #F
T2=460.0 #F
T3=1300.0 #F
T4=780.0 #F
#calculations
#"from mollier charts"
h1=129.1 #Btu/lb
h2 = 221.2 #Btu/lb
h3= 438.8 #Btu/lb
h4 = 301.5 #Btu/lb
wcom=h2-h1
wcob=h3-h2
wtur=h3-h4
eta=(wtur-wcom)/wcob *100
#results
print " work done by compressor in btu input as work per lb of air compressed",round(wcom,3)
print " Heat supplied in the combustor in Btu supplied per lb of air ",round(wcob,3)
print " work done in the turbine in Btu output as work per lb of air",round(wtur,3)
print " Cycle efficiency in percent",round(eta,3)
# Initialization of Variable
from math import pi
from math import atan
from numpy import *
T1=80 #F
T2=460 #F
T=700 #F
T3=1300 #F
T4=780 #F
#"from mollier charts"
h1=129.1 #Btu/lb
h2 = 221.2 #Btu/lb
h3= 438.8 #Btu/lb
h4 = 301.5 #Btu/lb
h=281.1 #Btu/lb
#calculations
wcom=h2-h1
wcob=h3-h2
wtur=h3-h4
output=-wcom+wtur
Q=h3-h
eff=output/Q *100
#results
print " Heat supplied in the combustor in Btu supplied per lb of air ",round(Q,3)
print " Cycle efficiency in percent",round(eff,3)