Chemical mission to Pluto

%--------------------------------------------------------------------------
% See also: VOrbit, RPRA2AE, Plot2D
%--------------------------------------------------------------------------

%--------------------------------------------------------------------------
%   Copyright 2017 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------

d             = 6.09e9; % km to Pluto (at 40.7 AU)

uE            = 324*9.806/1000;
fS            = 0.02;
muPluto       = Constant('mu pluto');
rPluto        = Constant('equatorial radius pluto');
rA            = rPluto + 1000;
rP            = rPluto + 200;
a             = RPRA2AE(rP,rA);
vOrbit        = VOrbit(rA,a,muPluto);

dVMax         = uE*log((1+fS)/fS);
vMax          = dVMax - vOrbit;
vTransfer     = linspace(5,0.99*vMax);
tTransfer     = d./vTransfer/86400/365.25;
gamma         = exp((vTransfer-vOrbit)/uE);
mR            = (gamma-1)./(1+fS*(1-gamma));

Plot2D(tTransfer,mR,'Years','Mass Ratio' ,'Chemical Mission')


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

% $Id: 1042112fa2e7170d4b8c329e7b214cbf14ec4634 $

Published with MATLAB® R2024a