summaryrefslogtreecommitdiff
path: root/Trivac/data/NodalTests_proc/prob5p3_nem.c2m
blob: f895037c47cc34f6cba8c9a9815fadb6e1204c26 (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
*----
*  TEST CASE prob5p3_nem
*
*  REF: A. Hebert, "Applied Reactor Physics", Presses Internationales
*       Polytechnique, Problem 5.3 (2009).
*
*----
*  Define STRUCTURES and MODULES used
*----
LINKED_LIST GEOM MACRO TRACK FLUX ;
MODULE GEO: MAC: NSST: NSSF: DELETE: END: ;
PROCEDURE assertS ;
*----
*  Macroscopic cross sections
*----
MACRO := MAC: ::
 EDIT 2 NGRO 2 NMIX 2 NIFI 1
 READ INPUT
 MIX     1
      DIFF  1.264E+00  0.9328E+00
     TOTAL  8.154E-03  4.1000E-03
    NUSIGF  0.000E+00  4.5620E-03
       CHI  1.000E+00  0.000E+00
  H-FACTOR  0.000E+00  4.5620E-03
      SCAT  1 1 0.0 2 2 0.0   7.368E-03
 MIX     2
      DIFF  1.310E+00  0.8695E+00
     TOTAL  1.018E-02  2.1170E-04
      SCAT  1 1 0.0 2 2 0.0   1.0180E-02
 ;
*----
*  Steady-state calculations with the Nodal Expansion Method
*----
GEOM := GEO: :: CAR1D 4
           X- ZERO X+ ZERO
           MIX  2 1 1 2
           MESHX 0.0 40.0 350.0 660.0 700.0
           SPLITX 1 4 4 1
           ;
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9977308 ;
TRACK FLUX := DELETE: TRACK FLUX ;
*
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 CMFD ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9976518 ;
GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
*
GEOM := GEO: :: CAR1D 4
           X- VOID X+ VOID
           MIX  2 1 1 2
           MESHX 0.0 40.0 350.0 660.0 700.0
           SPLITX 1 4 4 1
           ;
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9977973 ;
TRACK FLUX := DELETE: TRACK FLUX ;
*
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 CMFD ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9977070 ;
GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
*
GEOM := GEO: :: CAR1D 4
           X- ALBE 0.85 X+ ALBE 0.95
           MIX  2 1 1 2
           MESHX 0.0 40.0 350.0 660.0 700.0
           SPLITX 1 4 4 1
           ;
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9988235 ;
TRACK FLUX := DELETE: TRACK FLUX ;
*
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 CMFD ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9986206 ;
GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
*
GEOM := GEO: :: CAR1D 2
           X- VOID X+ REFL
           MIX  2 1
           MESHX 0.0 40.0 350.0
           SPLITX 1 4
           ;
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9977973 ;
TRACK FLUX := DELETE: TRACK FLUX ;
*
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 CMFD ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9977079 ;
GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
*
MACRO := MAC: MACRO ::
  ALBP 2 0.8 0.9 0.7 0.75
 ;
GEOM := GEO: :: CAR1D 4
           X- ALBE 1 X+ ALBE 2
           MIX  2 1 1 2
           MESHX 0.0 40.0 350.0 660.0 700.0
           SPLITX 1 4 4 1
           ;
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9983991 ;
TRACK FLUX := DELETE: TRACK FLUX ;
*
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 HYPE 2 CMFD ;
FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
assertS FLUX :: K-EFFECTIVE 1 0.9982761 ;
GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
*----
*  Steady-state calculations with the Analytic Nodal Method
*----
GEOM := GEO: :: CAR1D 3
           X- ZERO X+ ZERO
           MIX  2 1 2
           MESHX 0.0 40.0 660.0 700.0
           ;
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 99 MAXR 20 ANM ;
FLUX := NSSF: TRACK MACRO :: 
    EDIT 1 NUPD 100 1.0E-7 EXTE 100 1.0E-5 LEAK quadratic ;
assertS FLUX :: K-EFFECTIVE 1 0.99773398 ;
GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
*
GEOM := GEO: :: CAR1D 3
           X- ALBE 1 X+ ALBE 2
           MIX  2 1 2
           MESHX 0.0 40.0 660.0 700.0
           ;
TRACK := NSST: GEOM :: 
    TITLE 'BENCHMARK ENE6103'
    EDIT 2 MAXR 20 ANM ;
FLUX := NSSF: TRACK MACRO ::
    EDIT 1 NUPD 100 1.0E-7 EXTE 100 1.0E-5 LEAK quadratic ;
assertS FLUX :: K-EFFECTIVE 1 0.9984057 ;
ECHO "test prob5p3_nem completed" ;
END: ;