blob: 4eaaeb2ba8529bebb9309ed2d6216b71cf2316ea (
plain)
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
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
|
*----
* TEST CASE TCM61
* 4-group eighth UOX assembly with AIC in DRAGON (ALLG option)
* Author: A. Hebert
*
* Ref: A. Hebert, "High Order Linear Discontinuous and Diamond
* Differencing Schemes along Cyclic Characteristics,"
* Nucl. Sci. Eng. 184, 591-603 (2016).
*----
*----
* Define STRUCTURES and MODULES used
*----
LINKED_LIST GEON2 MACRO4G TRACKN2 SYS FLUX2 EDIT ;
MODULE GEO: SNT: ASM: FLU: MAC: EDI: ABORT: UTL: OUT: ERROR: END: ;
SEQ_ASCII RCOND2 :: FILE './AIC_SN_CCOND2.ref' ;
REAL erra maxa ;
PROCEDURE assertS ;
*----
* Input geometry data
*----
*
REAL mesh1 := 0.37 ;
REAL mesh2 := mesh1 0.52 + ;
REAL mesh3 := mesh2 0.74 + ;
REAL mesh4 := mesh3 0.52 + ;
REAL mesh5 := mesh4 0.74 + ;
REAL mesh6 := mesh5 0.52 + ;
REAL mesh7 := mesh6 0.74 + ;
REAL mesh8 := mesh7 0.52 + ;
REAL mesh9 := mesh8 0.74 + ;
REAL mesh10 := mesh9 0.52 + ;
REAL mesh11 := mesh10 0.74 + ;
REAL mesh12 := mesh11 0.52 + ;
REAL mesh13 := mesh12 0.74 + ;
REAL mesh14 := mesh13 0.52 + ;
REAL mesh15 := mesh14 0.74 + ;
REAL mesh16 := mesh15 0.52 + ;
REAL mesh17 := mesh16 0.74 + ;
REAL mesh18 := mesh17 0.26 + ;
*
GEON2 := GEO: :: CAR2D 18 18
EDIT 0
X- DIAG X+ REFL
Y- REFL Y+ DIAG
MESHX 0.0 <<mesh1>> <<mesh2>> <<mesh3>> <<mesh4>> <<mesh5>> <<mesh6>>
<<mesh7>> <<mesh8>> <<mesh9>> <<mesh10>> <<mesh11>> <<mesh12>>
<<mesh13>> <<mesh14>> <<mesh15>> <<mesh16>> <<mesh17>>
<<mesh18>>
MESHY 0.0 <<mesh1>> <<mesh2>> <<mesh3>> <<mesh4>> <<mesh5>> <<mesh6>>
<<mesh7>> <<mesh8>> <<mesh9>> <<mesh10>> <<mesh11>> <<mesh12>>
<<mesh13>> <<mesh14>> <<mesh15>> <<mesh16>> <<mesh17>>
<<mesh18>>
SPLITX 2 2 4 2 4 2 4 2 4 2 4 2 4 2 4 2 4 1
SPLITY 2 2 4 2 4 2 4 2 4 2 4 2 4 2 4 2 4 1
MIX
1 1 2 1 2 1 4 1 2 1 2 1 4 1 2 1 3 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 1 3 1 2 1 3 1 3 1 2 1 3 1 3 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
3 1 2 1 3 1 3 1 2 1 3 1 3 1
1 1 1 1 1 1 1 1 1 1 1 1 1
4 1 2 1 2 1 4 1 2 1 3 1
1 1 1 1 1 1 1 1 1 1 1
3 1 2 1 2 1 3 1 3 1
1 1 1 1 1 1 1 1 1
4 1 2 1 3 1 3 1
1 1 1 1 1 1 1
3 1 3 1 3 1
1 1 1 1 1
3 1 3 1
1 1 1
3 1
1
;
*----
* Input Macrolib
*----
MACRO4G := MAC: ::
EDIT 1 NGRO 4 NMIX 4 NIFI 1 ANIS 2
READ INPUT
MIX 1 (* WATER *)
TOTAL 5.316E-01 9.352E-01 9.931E-01 1.5931E+00
TRANC 2.880E-01 5.529E-01 4.936E-01 4.4340E-01
SCAT 1 1 (*1->1*) 4.973E-01
2 2 (*2->2*) 6.570E-01 (*1->2*) 3.134E-02
4 4 (*4->3*) 1.163E-02 (*3->3*) 7.176E-01 (*2->3*) 2.528E-01
(*1->3*) 2.304E-03
4 4 (*4->4*) 1.565E+00 (*3->4*) 2.707E-01 (*2->4*) 2.412E-02
(*1->4*) 2.228E-04
1 1 (*1->1*) 2.735E-01
2 2 (*2->2*) 4.253E-01 (*1->2*) 1.410E-02
4 4 (*4->3*) 5.064E-03 (*3->3*) 4.286E-01 (*2->3*) 1.203E-01
(*1->3*) 3.050E-04
4 4 (*4->4*) 4.384E-01 (*3->4*) 6.426E-02 (*2->4*) 3.373E-03
(*1->4*) 1.245E-05
MIX 2 (* FUEL 1 *)
TOTAL 4.097E-01 6.286E-01 4.594E-01 7.508E-01
TRANC 4.841E-02 7.939E-03 8.080E-03 8.038E-03
NUSIGF 1.632E-02 9.259E-02 1.129E-01 6.508E-01
CHI 1.000E+00 0.000E+00 0.000E+00 0.000E+00
SCAT 1 1 (*1->1*) 3.908E-01
2 2 (*2->2*) 4.407E-01 (*1->2*) 8.546E-04
3 4 (*4->3*) 5.318E-03 (*3->3*) 3.785E-01 (*2->3*) 8.057E-03
2 4 (*4->4*) 3.983E-01 (*3->4*) 1.203E-02
1 1 (*1->1*) 4.793E-02
2 2 (*2->2*) 5.419E-03 (*1->2*) -2.640E-04
3 4 (*4->3*) -7.338E-04 (*3->3*) 9.014E-03 (*2->3*) -2.453E-03
2 4 (*4->4*) 8.769E-03 (*3->4*) -2.572E-03
MIX 3 (* FUEL 2 *)
TOTAL 4.072E-01 6.223E-01 4.599E-01 7.504E-01
TRANC 4.927E-02 7.872E-03 8.166E-03 8.040E-03
NUSIGF 1.652E-02 9.136E-02 1.143E-01 6.500E-01
CHI 1.000E+00 0.000E+00 0.000E+00 0.000E+00
SCAT 1 1 (*1->1*) 3.885E-01
2 2 (*2->2*) 4.387E-01 (*1->2*) 8.562E-04
3 4 (*4->3*) 5.236E-03 (*3->3*) 3.780E-01 (*2->3*) 8.591E-03
2 4 (*4->4*) 3.984E-01 (*3->4*) 1.229E-02
1 1 (*1->1*) 4.981E-02
2 2 (*2->2*) 3.100E-03 (*1->2*) -2.624E-04
3 4 (*4->3*) -7.231E-04 (*3->3*) 1.142E-02 (*2->3*) -2.610E-03
2 4 (*4->4*) 8.764E-03 (*3->4*) -2.604E-03
MIX 4 (* AIC *)
TOTAL 4.240E-01 1.0476E+00 2.7503E+00 1.1194E+01
TRANC 0.000E+00 0.000E+00 0.000E+00 0.000E+00
SCAT 1 1 (*1->1*) 3.738E-01
2 2 (*2->2*) 4.102E-01 (*1->2*) 2.439E-04
3 4 (*4->3*) 4.985E-03 (*3->3*) 4.956E-01 (*2->3*) 7.370E-04
2 4 (*4->4*) 2.996E-01 (*3->4*) 3.248E-03
1 1 (*1->1*) 0.000E+00
1 2 (*2->2*) -6.161E-02
1 3 (*3->3*) -2.326E-01
1 4 (*4->4*) 0.000E+00
;
TRACKN2 := SNT: GEON2 ::
TITLE 'SECOND LEVEL TRACKING'
EDIT 2 MAXR 25000 DIAM 0 SN 16 SCAT 2
MAXI 50 EPSI 1.E-5 GMRES 10
QUAD 2 (*level symmetric*)
ALLG ;
SYS := ASM: MACRO4G TRACKN2 :: ARM EDIT 1 ;
FLUX2 := FLU: MACRO4G SYS TRACKN2 ::
EDIT 1 TYPE K EXTE 100 1E-5 ;
assertS FLUX2 :: 'K-EFFECTIVE' 1 0.9171573 ;
EDIT := OUT: FLUX2 TRACKN2 MACRO4G GEON2 ::
EDIT 2 INTG
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 107 108 109 110 111 112 113 114 115 116
117 118 119 120 121 122 123 124 125 126
127 128 129 130 131 132 133 134 135
136 137 138 139 140 141 142 143
144 145 146 147 148 149 150
151 152 153 154 155 156
157 158 159 160 161
162 163 164 165
166 167 168
169 170
171
;
UTL: EDIT :: DIR IMPR STATE-VECTOR * ;
ERROR: RCOND2 EDIT :: PICK ERASUM >>erra<< ERAMAX >>maxa<< ;
ECHO "erra=" erra "maxa=" maxa ;
*----
* Assert statistic on erra
*----
REAL REFVALUE := 0.2793 ;
REAL DELTA := erra REFVALUE - REFVALUE / ABS ;
IF DELTA 1.0E-2 < THEN
ECHO "TEST SUCCESSFUL; DELTA=" DELTA ;
ELSE
ECHO "------------" ;
ECHO "TEST FAILURE" ;
ECHO "------------" ;
ECHO "REFERENCE=" REFVALUE " CALCULATED=" erra ;
ABORT: ;
ENDIF ;
ECHO "test TCM61 completed" ;
END: ;
|
