-
Notifications
You must be signed in to change notification settings - Fork 3
/
all_ana.m
92 lines (82 loc) · 2.58 KB
/
all_ana.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
input_name = "infinite";
run(input_name + "_det0.m")
run(input_name + "_his0.m");
% ----------------------------------------------------------------
% Plotting Total and Material Specific Flux
% ----------------------------------------------------------------
path_images = "./Flux_Images/";
% Scale the energy integrated flux to a maximum of 1.0
DETEnergyDetector(:,11) = DETEnergyDetector(:,11)/max(DETEnergyDetector(:,11));
%DETspectrum(:,11) = DETspectrum(:,11)/max(DETspectrum(:,11));
figure('visible','off');
errorbar(DETEnergyDetectorE(:,3), DETEnergyDetector(:,11), 2*DETEnergyDetector(:,11).*DETEnergyDetector(:,12), 'k.');
%plot(DETEnergyDetectorE(:,3), DETEnergyDetector(:,11));
legend({'U-235 Infinite System'}, 'Location', 'southwest')
set(gca, 'XScale', 'log');
set(gca, 'YScale', 'log');
set(gca, 'XTick', [1e-12,1e-10,1e-8,1e-6,1e-4,1e-2,1e0,1e2]);
set(gca, 'FontSize', 16);
xlabel('Energy (MeV)')
ylabel('Energy integrated neutron flux (a.u.)')
grid on
box on
ylim([0,1.1]);
print(gcf, (path_images + "Flux_Log"), '-dpng', '-r300')
set(gca,'YScale','linear');
print(gcf, (path_images + "Flux_Linear"), '-dpng', '-r300')
% ----------------------------------------------------------------
% Plotting K-eff values per cycle and relative error
% ----------------------------------------------------------------
path_images = "./Convergance_Images/";
figure('visible','off');
hold on
errorbar(HIS_ANA_KEFF(:, 1), HIS_ANA_KEFF(:,2), ...
2*HIS_ANA_KEFF(:,2).*HIS_ANA_KEFF(:,4), 'K.');
set(gca, 'FontSize', 16);
title('Total K-eff (Analog)')
xlabel('Cycle Number')
ylabel('K-eff')
grid on
box on
hold off
print(gcf,(path_images + "total_keff_ana"),'-dpng','-r300')
clf
figure('visible','off');
hold on
errorbar(HIS_ANA_KEFF(:, 1), HIS_ANA_KEFF(:,5), ...
2*HIS_ANA_KEFF(:,5).*HIS_ANA_KEFF(:,7), 'K.');
set(gca,'FontSize',16);
title('Total Prompt')
xlabel('Cycle Number')
ylabel('K-eff')
grid on
box on
hold off
print(gcf,(path_images + "prompt_keff"),'-dpng','-r300')
clf
figure('visible','off');
hold on
errorbar(HIS_ANA_KEFF(:, 1), HIS_ANA_KEFF(:,8), ...
2*HIS_ANA_KEFF(:,8).*HIS_ANA_KEFF(:,10), 'K.');
set(gca,'FontSize',16);
title('Total Delayed')
xlabel('Cycle Number')
ylabel('K-eff')
grid on
box on
hold off
print(gcf,(path_images + "delayed_keff"),'-dpng','-r300')
clf
figure('visible','off');
hold on
errorbar(HIS_IMP_KEFF(:, 1), HIS_IMP_KEFF(:,2), ...
2*HIS_IMP_KEFF(:,2).*HIS_IMP_KEFF(:,4), 'K.');
set(gca,'FontSize',16);
title('Total K-eff (Implicit)')
xlabel('Cycle Number')
ylabel('K-eff')
grid on
box on
hold off
print(gcf,(path_images + "total_keff_imp"),'-dpng','-r300')
clf