1
2
3
4
5
6
7
8
9
10
11
12
13
14
|
from rocketcea.cea_obj import CEA_Obj
ispObj = CEA_Obj( oxName='LOX', fuelName='LH2')
for frozen, frozenAtThroat in [(0,0), (1,0), (1,1)]:
molWtD, massFracD = ispObj.get_SpeciesMassFractions( Pc=1000.0, MR=6.0, eps=40.0,
frozen=frozen, frozenAtThroat=frozenAtThroat )
print(' ROCKETCEA MASS FRACTIONS (frozen=%i, frozenAtThroat=%i)'%(frozen, frozenAtThroat) )
for species, mfL in massFracD.items():
s = ' '.join( [ '%7.5f'%mf for mf in mfL ] )
print( ' %-15s '%species, s, ' MW=%g'%molWtD[species] )
print('='*55)
|
