Contents

Explore D-3He fusion power model

Explore combinations of D and 3He temperature and fuel ratios.

%--------------------------------------------------------------------------
% See also PowerDHe3
%--------------------------------------------------------------------------

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

nE = linspace(5e19,7e20,25);

d = PowerDHe3;

tE = 30;
ratios = [1 2 3];
temps = [140 120 100;70 90 100];
for i = 1:length(ratios)
  for j = 1:size(temps,2)
    d.tHe3 = temps(1,j);
    d.tD   = temps(2,j);
    for k = 1:length(nE)
      d.nD = nE(k)/(2*ratios(i)+1);
      d.nHe3 = nE(k)/(2+1/ratios(i));
      [pF(i,j,k), pN(i,j,k), pB(i,j,k), nNd(i,j,k), nNt(i,j,k)] = PowerDHe3( d );
    end
  end
end

Plot with subplots per

NewFig('Fusion Power with Ion Temp')
subplot(3,1,1)
loglog(nE,squeeze(pF(:,1,:)))
grid on
ylabel('Temp 140/70')
title('Fusion Power MW/m3')
subplot(3,1,2)
loglog(nE,squeeze(pF(:,2,:)))
grid on
ylabel('Temp 120/90')
subplot(3,1,3)
loglog(nE,squeeze(pF(:,3,:)))
grid on
xlabel('Electron Density #/m3')
legend('Ratio 1','Ratio 2','Ratio 3')
ylabel('Temp 100/100')

2nd plot

NewFig('Fusion Power with Ratios')
subplot(3,1,1)
loglog(nE,squeeze(pF(1,:,:)))
grid on
ylabel('Ratio 1')
title('Fusion Power MW/m3')
subplot(3,1,2)
loglog(nE,squeeze(pF(2,:,:)))
grid on
ylabel('Ratio 2')
subplot(3,1,3)
loglog(nE,squeeze(pF(3,:,:)))
grid on
ylabel('Ratio 3')
xlabel('Electron Density #/m3')
legend('T 140/70','T 120/90','T 100/100')

Combined plot

NewFig('Fusion Power - Combined Effects')
loglog(nE,squeeze(pF(1,1,:)),'linewidth',2)
hold on
loglog(nE,squeeze(pF(1,2,:)),'linewidth',1.5)
loglog(nE,squeeze(pF(1,3,:)),'linewidth',1)
colors = get(gca,'colororder');
loglog(nE,squeeze(pF(2,1,:)),'--','color',colors(1,:),'linewidth',2)
loglog(nE,squeeze(pF(2,2,:)),'--','color',colors(2,:),'linewidth',1.5)
loglog(nE,squeeze(pF(2,3,:)),'--','color',colors(3,:),'linewidth',1)
xlabel('Electron Density #/m3')
ylabel('Fusion Power MW/m3')
loglog(nE,squeeze(pF(3,1,:)),'-.','color',colors(1,:),'linewidth',2)
loglog(nE,squeeze(pF(3,2,:)),'-.','color',colors(2,:),'linewidth',1.5)
loglog(nE,squeeze(pF(3,3,:)),'-.','color',colors(3,:),'linewidth',1)
legend('rho=1, T=140/70',...
       'rho=1, T=120/90',...
       'rho=1, T=100/100',...
       'rho=2, T=140/70',...
       'rho=2, T=120/90',...
       'rho=2, T=100/100',...
       'rho=3, T=140/70',...
       'rho=3, T=120/90',...
       'rho=3, T=100/100')

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

Published with MATLAB® R2020a