amm_dissociation.py

 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
28
from rocketcea.cea_obj import CEA_Obj
from pylab import *

Pc = 200.0
eps = 20.0

xL = []
ispL = []

for x in range(10, 100, 5):
    propName = 'HYD%g'%x

    ispObj = CEA_Obj(propName=propName)
    
    xL.append( x )
    
    IspFroz, Cstar, Tcomb = ispObj.getFrozen_IvacCstrTc( Pc=Pc, eps=eps, frozenAtThroat=1 )
    ispL.append( IspFroz )

plot(xL, ispL, linewidth=4)

grid(True)
title( 'Hydrazine Ideal Performance vs. Ammonia Dissociation\nPc=%g psia, Area Ratio=%g'%(Pc, eps) )
xlabel( '% Ammonia Dissociation' )
ylabel( 'IspVac (sec)' )
savefig('amm_dissociation.png', dpi=120)

show()