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
26
27
|
from rocketcea.cea_obj import CEA_Obj
from pylab import *
pcL = [ 2000., 500., 70.]
ispObj = CEA_Obj(propName='', oxName='LOX', fuelName="LH2")
for Pc in pcL:
cstarArr = []
MR = 2.0
mrArr = []
while MR < 8.0:
cstarArr.append( ispObj.get_Cstar( Pc=Pc, MR=MR) )
mrArr.append(MR)
MR += 0.05
plot(mrArr, cstarArr, label='Pc=%g psia'%Pc)
legend(loc='best')
grid(True)
title( ispObj.desc )
xlabel( 'Mixture Ratio' )
ylabel( 'Cstar (ft/sec)' )
savefig('cea_cstar_plot.png', dpi=120)
show()
|
