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

    Pc = 500.0
    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()