blob: 134cc206f1709d63fd642087b7f2f66ee053c306 (
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
|
*----
* TEST CASE C2D20
* ACR ANNULAR CELL WITH NXT
*----
MODULE SHI: ASM: FLU: DELETE: END: ;
PROCEDURE GeoACRC2DN ;
PROCEDURE MicACRE620 MicACRE6UH MicACRE6UW ;
LINKED_LIST CFLXS.tds CFLXF.tds Flux Compo1 ;
SEQ_BINARY CFLXS.trk CFLXF.trk ;
XSM_FILE MicLib ConcIso EditDS Pij ;
STRING Modele Option Etat BCtrk DiscrOpt TypeMRG ;
REAL FuelSplit CoolSplit PTSplit CTSplit ModeSplit ;
PROCEDURE assertS ;
*----
* Get cross sections
*----
MicLib := MicACRE620 ;
MicLib := MicACRE6UH MicLib ;
MicLib := MicACRE6UW MicLib ;
*----
* Initialyze geometry options
*----
INTEGER NbAngles := 10 ;
REAL TrkDens := 20.0 ;
EVALUATE Etat BCtrk DiscrOpt := "Present" "TISO" "Thickness" ;
*----------
* Get geometry for resonance self shielding
*----------
EVALUATE Modele := "C2D6" ;
EVALUATE TypeMRG := "Supercell" ; ! "Supercell" OR "Reactor"
EVALUATE FuelSplit CoolSplit PTSplit CTSplit ModeSplit :=
0.0 0.0 0.0 0.0 0.0 ;
EVALUATE Option := "Shield" ;
CFLXS.tds CFLXS.trk := GeoACRC2DN ::
<<Modele>> <<Option>> <<Etat>> <<BCtrk>>
<<NbAngles>> <<TrkDens>>
;
*----------
* Get coarse geometry for flux calculations
* Uncomment EVALUATE command for fine mesh.
*----------
EVALUATE FuelSplit CoolSplit PTSplit CTSplit ModeSplit :=
-0.25 0.1 0.0 0.0 0.3 ;
EVALUATE Option := "Flux" ;
CFLXF.tds CFLXF.trk := GeoACRC2DN ::
<<Modele>> <<Option>> <<Etat>> <<BCtrk>> <<DiscrOpt>>
<<FuelSplit>> <<CoolSplit>> <<PTSplit>> <<CTSplit>> <<ModeSplit>>
<<NbAngles>> <<TrkDens>>
;
*----
* Self-shielding Calculation
*----
MicLib := SHI: MicLib CFLXS.tds CFLXS.trk ::
EDIT 0 LJ ;
*----
* Flux Calculation with SKIP
*----
Pij := ASM: MicLib CFLXF.tds CFLXF.trk ::
EDIT 0 PIJK SKIP PNOR HELI ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 SIGS ;
assertS Flux :: 'K-INFINITY' 1 1.165185 ;
Flux := DELETE: Flux ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 PNL ;
*assertS Flux :: 'K-INFINITY' 1 1.165915 1.165582 ;
Flux := DELETE: Flux ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 PNLR ;
*assertS Flux :: 'K-INFINITY' 1 1.165098 1.165501 ;
Flux := DELETE: Flux ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 TIBERE ;
assertS Flux :: 'K-INFINITY' 1 1.164401 ;
Pij Flux := DELETE: Pij Flux ;
Pij := ASM: MicLib CFLXF.tds CFLXF.trk ::
EDIT 0 PIJ SKIP PNOR HELI ECCO ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 ECCO ;
assertS Flux :: 'K-INFINITY' 1 1.164370 ;
Pij Flux := DELETE: Pij Flux ;
*----
* Flux Calculation without SKIP
*----
Pij := ASM: MicLib CFLXF.tds CFLXF.trk ::
EDIT 0 PIJK PNOR HELI ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 SIGS ;
assertS Flux :: 'K-INFINITY' 1 1.165184 ;
Flux := DELETE: Flux ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 PNL ;
assertS Flux :: 'K-INFINITY' 1 1.165599 ;
Flux := DELETE: Flux ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 PNLR ;
assertS Flux :: 'K-INFINITY' 1 1.165524 ;
Flux := DELETE: Flux ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 TIBERE ;
assertS Flux :: 'K-INFINITY' 1 1.164400 ;
Pij Flux := DELETE: Pij Flux ;
Pij := ASM: MicLib CFLXF.tds CFLXF.trk ::
EDIT 0 PIJ PNOR HELI ECCO ;
Flux := FLU: Pij MicLib CFLXF.tds ::
EDIT 1 TYPE B B1 ECCO ;
assertS Flux :: 'K-INFINITY' 1 1.164370 ;
*----
* Cleanup
*----
Pij Flux := DELETE: Pij Flux ;
CFLXS.tds CFLXF.tds CFLXS.trk CFLXF.trk := DELETE:
CFLXS.tds CFLXF.tds CFLXS.trk CFLXF.trk ;
MicLib := DELETE: MicLib ;
ECHO "test C2D20 completed" ;
END: ;
QUIT "LIST" .
|