blob: 5615cf9741a60e175d7c2a6398d28833e68c069e (
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
|
*----
* TEST CASE TCM51
* 3-REGION PIJ AND SN WITH DOUBLE HETEROGENEITY (BIHET) MIXTURES
* ONE-GROUP MACROLIB-DEFINED CROSS SECTIONS
*----
* Define STRUCTURES and MODULES used
*----
LINKED_LIST MACRO ASSMB TRACK SYS FLUX EDIT ;
MODULE GEO: SYBILT: SNT: MAC: ASM: FLU: UTL: DELETE: GREP: END: ;
REAL keff ;
SEQ_BINARY TRACK_FIL ;
PROCEDURE assertS ;
*
MACRO := MAC: ::
EDIT 2 NGRO 1 NMIX 4 NIFI 1
READ INPUT
MIX 1
TOTAL 0.41 SCAT 1 1 0.3
NUSIGF 0.1 CHI 1.0
MIX 2
TOTAL 0.45 SCAT 1 1 0.347
NUSIGF 0.17 CHI 1.0
MIX 3
TOTAL 0.53 SCAT 1 1 0.5
MIX 4
TOTAL 0.3 SCAT 1 1 0.2
NUSIGF 0.1 CHI 1.0
;
ASSMB := GEO: :: TUBE 3
R+ REFL
MIX 5 4 6
RADIUS 0.0 1.0 3.0 5.0
BIHET SPHE (*NG=*) 2 (* NMILG= *) 2 (* SPHERICAL MICRO-STRUCTURE *)
(* NS= *) 1 1
(* M-S-1 *) 0.0 0.3 (* M-S 2 *) 0.0 0.5
(* COMPOSITE MIXTURES *) 5 6
(* MIXTURES SURROUNDING M-S *) 3 1
(* COMPOSITE MIXTURE 5 FRACT *) 0.4 0.0
(* REAL MIXTURE CONTENT M-S-1 *) 3
(* COMPOSITE MIXTURE 6 FRACT *) 0.0 0.3
(* REAL MIXTURE CONTENT M-S-1 *) 2
;
*----
* Collision probability method
*----
TRACK := SYBILT: ASSMB ::
EDIT 1
MAXR 20
HEBE QUAB 7
TITLE 'DOUBLE-HETEROGENEITY TEST'
;
*
SYS := ASM: MACRO TRACK :: EDIT 99 PIJ ;
FLUX := FLU: MACRO TRACK SYS :: TYPE K ;
GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>keff<< ;
ECHO "HEBE K-EFFECTIVE=" keff ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.074210 ;
TRACK SYS FLUX := DELETE: TRACK SYS FLUX ;
*----
* Discrete ordinates method HEBE
*----
TRACK := SNT: ASSMB ::
EDIT 2
DIAM 0 SN 8 SCAT 1
HEBE QUAB 7
TITLE 'DOUBLE-HETEROGENEITY TEST'
;
*
SYS := ASM: MACRO TRACK :: EDIT 99 ARM ;
FLUX := FLU: MACRO TRACK SYS :: EDIT 99 TYPE K ;
GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>keff<< ;
ECHO "HEBE K-EFFECTIVE=" keff ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.073958 ;
TRACK SYS FLUX := DELETE: TRACK SYS FLUX ;
*
*----
* Collision probability method SAPO
*----
TRACK := SYBILT: ASSMB ::
EDIT 1
MAXR 20
SAPO QUAB 7
TITLE 'DOUBLE-HETEROGENEITY TEST'
;
*
SYS := ASM: MACRO TRACK :: EDIT 99 PIJ ;
FLUX := FLU: MACRO TRACK SYS :: TYPE K ;
GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>keff<< ;
ECHO "SAPO K-EFFECTIVE=" keff ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.074108 ;
TRACK SYS FLUX := DELETE: TRACK SYS FLUX ;
*
*----
* Discrete ordinates method SAPO
*----
TRACK := SNT: ASSMB ::
EDIT 2
DIAM 0 SN 8 SCAT 1
SAPO QUAB 7
TITLE 'DOUBLE-HETEROGENEITY TEST'
;
*
SYS := ASM: MACRO TRACK :: EDIT 99 ARM ;
FLUX := FLU: MACRO TRACK SYS :: EDIT 99 TYPE K ;
GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>keff<< ;
ECHO "SAPO K-EFFECTIVE=" keff ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.073781 ;
TRACK SYS FLUX := DELETE: TRACK SYS FLUX ;
*
*----
* Collision probability method SLSI
*----
TRACK := SYBILT: ASSMB ::
EDIT 1
MAXR 20
SLSI 0.65 QUAB -100
TITLE 'DOUBLE-HETEROGENEITY TEST'
;
*
SYS := ASM: MACRO TRACK :: EDIT 99 PIJ ;
FLUX := FLU: MACRO TRACK SYS :: TYPE K ;
GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>keff<< ;
ECHO "SLSI K-EFFECTIVE=" keff ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.068906 ;
TRACK SYS FLUX := DELETE: TRACK SYS FLUX ;
*
*----
* Discrete ordinates method SLSI
*----
TRACK := SNT: ASSMB ::
EDIT 2
DIAM 0 SN 8 SCAT 1
SLSI 0.65 QUAB -100
TITLE 'DOUBLE-HETEROGENEITY TEST'
;
*
SYS := ASM: MACRO TRACK :: EDIT 99 ARM ;
FLUX := FLU: MACRO TRACK SYS :: EDIT 99 TYPE K ;
GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 1 1 >>keff<< ;
ECHO "SLSI K-EFFECTIVE=" keff ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.068608 ;
TRACK SYS FLUX := DELETE: TRACK SYS FLUX ;
ECHO "test TCM51 completed" ;
END: ;
QUIT "LIST" .
|