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 %--------------------------------------