## 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()