DV vs trip time demo for specific missions

See also: DVConstantThrust

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

Ps = [1 10]*1e6;
payload  = 1000;    % kg
uExhaust = 100;    % km/s
sigma    = 1000;   % W/kg
eta      = 0.4;    % fusion power to thrust
sR       = 0.05;

durations = linspace(1,10,20);
year = 365.25*86400;    % s

% loop over durations
DV = [];
for j = 1:2
  P = Ps(j);
  thrust = 2*P*eta/(uExhaust*1e3);
  mEngine = P/sigma;
  for k = 1:length(durations)
    duration = durations(k)*year;
    mFuel = thrust/(uExhaust*1e3)*duration;
    DV(j,k) = DVConstantThrust( thrust, payload+mEngine+sR*mFuel, mFuel, uExhaust, duration );
  end
end

Plot2D(durations,DV*1e-3,'Trip Time (years)','Delta-V (km/s)','DFD Capability')
legend('1 MW','10 MW')

yy = axis;
hold on
plot(yy(1:2),100*[1 1],'r');
plot(yy(1:2),150*[1 1],'r');
plot(yy(1:2),250*[1 1],'r');


%--------------------------------------
% $Date$
% $Id: 5e624be579b95c3e37bc56adcd2eb63c3a9322f4 $

Published with MATLAB® R2024a