1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
from rocketcea.cea_obj import CEA_Obj
from pylab import *
Pc = 500.0
ispIRFNA = CEA_Obj(propName='', oxName='IRFNA', fuelName="MHF3")
for e in [50.0,20.0,10.0]:
ispArr = []
MR = 1.1
mrArr = []
while MR < 3.5:
ispArr.append( ispIRFNA(Pc, MR, e ))
mrArr.append(MR)
MR += 0.05
plot(mrArr, ispArr, label='AreaRatio %g'%e)
legend(loc='best')
grid(True)
title( ispIRFNA.desc )
xlabel( 'Mixture Ratio' )
ylabel( 'Isp ODE (sec)' )
savefig('cea_plot.png', dpi=120)
show()
|
