import math
# given data
delG=-39.59 # kJ/mol
delH=-56.83 # change in enthalpy in kJ/mol
mdotmethanol=32.0 # in g/s
mdotoxygen=48.0 # in g/s
Wmax=166.3 # -delG in kJ
flowmethanol=mdotmethanol*100*3600/(Wmax*1000) # in kg/h
flowoxygen=mdotoxygen*100*3600/(Wmax*1000) # in kg/h
print "The required flow rate of methanol is ",round(flowmethanol,2),"kg/h"
print "The required flow rate of oxygen is ",round(flowoxygen,2),"kg/h"
delQ=delH-delG # using eq 12.7
fuelrate=-delQ*19.24/mdotmethanol # in kcal/s
print "The required heat removal rate is ",round(fuelrate,3),"kcal/s"
# The answer is wrong in textbook because of incorrect conversion from g/s to kg/h