blob: 4fca2b2f321bf67ba030991eeb8ba255f877d457 (
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
|
* Regression test for SENS: module based on a pin cell calculation.
* jef2.2 library used
* Running time = ~ 30sec
* C. Laville
*
LINKED_LIST
GEOMETRY ASSEMBLY TRACKING FLUX MACLIB MICLIB EDIT
ADJ_ASSEMBLY ADJ_FLUX ADJ_MACLIB ADJ_EDIT ;
MODULE
LIB: GEO: EDI: SYBILT: USS: ASM: FLU: END: T: DELETE: GREP:
SENS: ;
PROCEDURE
assertS ;
REAL kiso b2 ;
SEQ_ASCII SENS.sdf :: FILE './SENS.sdf' ;
*
*--------------------------------------------------------
* Calcul flux direct
*--------------------------------------------------------
*
MICLIB := LIB: ::
NMIX 8
CTRA APOL
ANIS 6
ADED 4 NELAS N4N N2N N3N
SUBG
MIXS LIB: DRAGON FIL: DLIB_J2
MIX 1 293.0
Pu238 7.38570E-06 1
Pu239 4.20190E-03 1
Pu240 5.60600E-04 1
Pu241 8.75230E-05 1
Pu242 1.69720E-05 1
U235 1.22230E-04 1
U238 1.68760E-02 1
Am241 3.42580E-05 1
O16 4.37130E-02
MIX 2 COMB 1 1.0
MIX 3 COMB 1 1.0
MIX 4 COMB 1 1.0
MIX 5 293.0
O16 1.00000E-06
MIX 6 293.0
Fe54 3.20914E-03
Fe56 5.07487E-02
Fe57 1.21726E-03
Fe58 1.54924E-04
Cr50 6.20381E-04
Cr52 1.20564E-02
Cr53 1.37463E-03
Cr54 3.42577E-04
Ni58 6.89729E-03
Ni60 2.66523E-03
Ni61 1.20916E-04
Ni62 3.71893E-04
Ni64 1.09739E-04
Mo0 2.11680E-03
Mn55 1.69350E-03
C0 6.77390E-05
Si0 8.46730E-04
MIX 7 293.0
O16 3.33680E-02
H1_H2O 6.67330E-02
MIX 8 COMB 7 1.0
;
*
GEOMETRY := GEO: :: TUBE 8
R+ REFL
RADIUS
0.00000
0.174666
0.220937
0.24076
0.247015
0.254
0.2921
0.418214
0.5374
MIX 1 2 3 4 5 6 7 8 ;
*
TRACKING := SYBILT: GEOMETRY :: ;
MACLIB := USS: MICLIB TRACKING :: TRAN PASS 2 ;
*
ASSEMBLY := ASM: MACLIB TRACKING :: PIJ ;
*
FLUX := FLU: ASSEMBLY MACLIB TRACKING :: TYPE K ;
GREP: FLUX :: GETVAL 'K-INFINITY' 1 1 1 >>kiso<< ;
IF kiso 1.0 > THEN
FLUX := DELETE: FLUX ;
FLUX := FLU: ASSEMBLY MACLIB TRACKING :: TYPE B B0TR ;
ENDIF ;
GREP: FLUX :: GETVAL 'B2 B1HOM' 1 1 1 >>b2<< ;
*
*--------------------------------------------------------
* Calcul flux adjoint
*--------------------------------------------------------
*
ADJ_MACLIB := T: MACLIB ;
ADJ_ASSEMBLY := ASM: ADJ_MACLIB TRACKING :: PIJ ;
*
ADJ_FLUX := FLU: ADJ_ASSEMBLY ADJ_MACLIB TRACKING ::
TYPE K B0TR PNL BUCK <<b2>> ;
*
*--------------------------------------------------------
* Calcul de sensibilité
*--------------------------------------------------------
*
SENS.sdf := SENS: MACLIB TRACKING FLUX ADJ_FLUX ::
EDIT 10 ANIS 1 ;
*
*
assertS FLUX :: 'K-EFFECTIVE' 1 1.0 ;
END: ;
|
