Compare engines

See also: Turbojet, StdAtm, RamjetKerrebrock, RDE

%--------------------------------------------------------------------------
%   Copyright (c) 2000 Princeton Satellite Systems, Inc. All rights reserved.
%--------------------------------------------------------------------------


leg = {'Turbojet' 'Ramjet' 'Turbojet/Afterburner' 'RDE'};

m   = linspace(0,8);
n 	= length(m);
h   = 20000*ones(1,n);
p 	= StdAtm(h);

iSp = zeros(4,n);
tS  = zeros(4,n);

d = Turbojet;
d.qR = 121e6;
[~,iSp(1,:),tS(1,:)]	= Turbojet( m, d, p );
[~,iSp(2,:),tS(2,:)]	= RamjetKerrebrock( m, d, p );
d.tT7 = 2000;
[~,iSp(3,:),tS(3,:)]	= Turbojet( m, d, p );
d = RDE;
[~,iSp(4,:),tS(4,:)]    = RDE(m,d,p);

%
% Plot2D(m,[tS;iSp], 'Mach', {'Specific Thrust','I_{sp} (s)'},'Engines with H2 Fuel','lin',{'1:3' '4:6'});
%  legend(leg,'location','best')
NewFig('EngineComparison');
subplot(2,1,1)
plot(m,tS)
xlabel('Mach')
ylabel('Specific Thrust')
legend(leg,'location','best')
subplot(2,1,2)
plot(m, iSp)
xlabel('Mach')
ylabel('I_{sp}')
legend('Engines with H2 Fuel','lin',{'1:3' '4:6'})
legend(leg,'location','best')

% PrintFig(1,4,1,'Engines'); % make a nice image


%--------------------------------------