Cryocooler model comparison

Compare the two models at 60 K, by computing input power.

%--------------------------------------------------------------------------
% See also: Cryocooler, CryocoolerBrake
%--------------------------------------------------------------------------

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


Tc = 60; % cold tip temperature
qC = logspace(-1,3); % lift power

[mass,eta,cOP] = Cryocooler( qC, Tc ); % stirling coolers, >50k
Pin = qC./cOP;
[massB,etaB,cOPB] = CryocoolerBrake( qC, Tc ); % more types of coolers
PinB = qC./cOPB;

Plot2D(qC,[mass;massB],'Qc (W)','Mass (kg)','Mass Comparison (40 K)','log')
legend('Ladner','ter Brake')

% Note: this seems much too low
Plot2D(Pin,mass,'Input Power','Mass (kg)','Ladner Model','log')

% note: this matches the paper
Plot2D(PinB,massB,'Input Power','Mass (kg)','ter Brake Model','log')
hold on
massB2 = 0.0711*PinB.^0.905;
plot(PinB,massB2,'r')
legend('Temp. Dependent','Power Only')

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

Published with MATLAB® R2020a