rp = RhoIspPlot(bipropL=[('LOX', 'LH2'), ('LOX', 'RP1')], Pc=1000., eps=16)
stg_obj = ReferenceStage(WtPayload=1600000.0,
volCuInRef=75500.0 * 1728,
WtInertRef=300000.0,
Name='Saturn V 1st stg')
rp.add_rho_isp_contours(label_frac_pos=0.2)
rp.add_stage_param_contours(
stg_obj,
set_param='GLOW',
param_value=5000000.0 + 1600000.0,
label_frac_posD={
'GLOW': 0.9,
'CubicFt': .4,
'MassFrac': .8
},
plot_param_valD={
'DeltaV': [12000, 13000, 14000, 15000],
'MassFrac': [0.86, 0.88, 0.9, 0.92, 0.94, 0.95],
'CubicFt': [50000, 75000, 100000, 125000, 150000, 175000, 200000]
},
plot_paramL=['DeltaV', 'CubicFt', 'MassFrac'],
num_ticks=6)
png_name = __file__[:-2] + 'png'
rp.savefig(png_name, dpi=120)
rp.show()
from rocketcea.biprop_utils.rho_isp_plot_obj import RhoIspPlot
from rocketcea.biprop_utils.veh_stage_obj import ConstMassFracStage
rp = RhoIspPlot(bipropL=[('LOX', 'LH2'), ('N2O4', 'MMH')],
nsteps_sg=90,
nsteps_isp=90)
stg_obj = ConstMassFracStage(mass_frac=0.8, WtPayload=1000.0)
rp.add_rho_isp_contours()
rp.add_stage_param_contours(stg_obj,
set_param='DeltaV',
param_value=7200.0,
plot_paramL=['GLOW'],
num_ticks=16)
rp.savefig('const_mass_frac1.png', dpi=120)
rp.show()
from rocketcea.biprop_utils.rho_isp_plot_obj import RhoIspPlot
from rocketcea.biprop_utils.veh_stage_obj import ReferenceStage
#rp = RhoIspPlot(bipropL=[('LOX','LH2'),('N2O4','MMH')], nsteps_sg=119, nsteps_isp=119)
rp = RhoIspPlot()
stg_obj = ReferenceStage(WtPayload=10000.0)
rp.add_rho_isp_contours(label_frac_pos=0.2)
rp.add_stage_param_contours(stg_obj,
set_param='DeltaV',
param_value=5000.0,
plot_param_valD={
'GLOW': [20000, 19000, 18000, 17000, 16000],
'MassFrac': [.65, .7, .75, .8, .85],
'CubicFt': [75, 100, 200, 300, 400]
},
label_frac_posD={
'GLOW': 0.9,
'CubicFt': .4
},
plot_paramL=['GLOW', 'CubicFt', 'MassFrac'],
num_ticks=6)
rp.savefig('rho_veh_1.png', dpi=120)
rp.show()
from rocketcea.biprop_utils.veh_stage_obj import ReferenceStage
stg_obj = ReferenceStage(volCuInRef=441.0,
WtInertRef=20.0,
WtPayload=30.0,
Name='N2H4 DACS')
rp = RhoIspPlot(bipropL=[('N2O4', 'MMH'), ('CLF5', 'N2H4')],
monopropL=['HYD40', 'HAN315', 'HAN269'])
rp.add_rho_isp_contours(label_frac_pos=0.2)
rp.add_stage_param_contours(
stg_obj,
set_param='GLOW',
param_value=64.0,
label_frac_posD={
'GLOW': 0.9,
'VolPropellant': .4
},
plot_param_valD={
'DeltaV': [1900, 2100, 2300, 2500, 2700, 2900, 3100, 3300],
'MassFrac': [.44, .46, .48, .5, .52],
'VolPropellant': [350, 375, 400, 425, 450]
},
plot_paramL=['DeltaV', 'VolPropellant'],
num_ticks=6)
rp.savefig('divert_veh.png', dpi=120)
rp.show()