Newer
Older
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
93
94
95
96
97
98
99
100
101
102
103
104
105
106
=%sampler
read parameters
n_samples=10
library="v7-56"
perturb_xs=yes
perturb_decay=yes
perturb_yields=yes
run_cases=yes
end parameters
read case [c1]
sequence=polaris
title "W17x17 pincell 80 GwD/MTU 8 wt%"
lib "v7-56"
system PWR
% increase pitch to account for GT and gap
opt ESSM NumAzim=1 NumPolar=1
opt KEFF PnOrder=0 TransportCorrection=no
EigSolver='power' PartitionUpscatter=no
UpscatterSolver='gauss_seidel' UpscatterInnerSolver='source'
geometry W17 : ASSM 1 1.32077
% compositions and materials
comp c_f50 : UOX enr=8.0
mat FUEL.1 : c_f50 dens=10.26
pin F : 0.4096 0.418 0.475
: FUEL.1 GAP.1 CLAD.1
mesh FUEL : nr=1
% state variables
state ALL : temp=590
COOL : dens=0.68 boron=800
FUEL : temp=900
CLAD : temp=700 TUBE : temp=580
pow 40
bu 1 5 10 20 30 40 50 60 70 80
end sequence
sequence=shell
mv ft71f001 discharge_iso
end sequence
sequence=origen
case
{
lib{ file=end7dec }
mat{ load{ file=discharge_iso pos=150 } }
time=[10L 1e-5 5000]
save=yes
}
end sequence
sequence=opus
data="ft71f001"
library="end7dec"
nrank=40
units=watts
end sequence
sequence=shell
rm ft71f001
end sequence
sequence=origen
case
{
lib{ file=end7dec }
mat{ load{ file=discharge_iso pos=154 } }
time=[10L 1e-5 5000]
save=yes
}
end sequence
sequence=opus
data="ft71f001"
library="end7dec"
nrank=40
units=watts
end sequence
end case
'response extracted from *.samplerfiles/c1_pert_*0.plt
read response[dh]
type=opus_plt
ndataset=0
nuclides=total end
end response
'response extracted from *.samplerfiles/c1_pert_*1.plt
read response[hbdh]
type=opus_plt
ndataset=1
nuclides=total end
end response
read analysis[pearson_steps2]
type = pearson_corr
targets= c1:dh.total(last) end
sources = c1:hbdh.total(last) end
end analysis
read analysis[pearson_set]
type = correlation_matrix
' targets= c1:dh.total(all) end
sources = c1:hbdh.total(all) end
end analysis
end