blob: 6103f55c625a0bf659fd793a9e3bd56dae482790 (
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
|
*----
* SmallCore_rt_mpo
*----
LINKED_LIST GEOM_OSC TRACK REFL MACRO SYS FLUX ;
MODULE GEO: TRIVAT: TRIVAA: FLUD: OUT: MCR: NCR: GREP: END: ;
HDF5_FILE CASEA :: FILE './Ass_3GLE_210-mpo.hdf' ;
HDF5_FILE CASEB :: FILE './Ass_3GLE_210_8Gd-mpo.hdf' ;
HDF5_FILE CASEC :: FILE './Ass_3GLE_320_16Gd-mpo.hdf' ;
HDF5_FILE CASED :: FILE './Ass_3GLE_420_16Gd-mpo.hdf' ;
SEQ_ASCII _REFL :: FILE './_OSC_reflector.txt' ;
PROCEDURE assertS ;
*----
* Local parameters
*----
REAL keff1 ;
REAL densB1 := 700.0 ; ! reference Boron concentration in ppm
REAL dens_mod := 0.7 ; ! reference coolant density (g/cc)
*
************************************************************************
* GEOMETRY AND TRACKING *
************************************************************************
REAL cote := 21.504 ;
REAL mesh1 := cote ;
REAL mesh2 := mesh1 cote + ;
REAL mesh3 := mesh2 cote + ;
REAL mesh4 := mesh3 cote + ;
REAL mesh5 := mesh4 cote + ;
REAL mesh6 := mesh5 cote + ;
REAL mesh7 := mesh6 cote + ;
REAL mesh8 := mesh7 cote + ;
REAL mesh9 := mesh8 cote + ;
REAL mesh10 := mesh9 cote + ;
REAL mesh11 := mesh10 cote + ;
*
GEOM_OSC := GEO: :: CAR2D 11 11
EDIT 0
X- VOID X+ VOID
Y- VOID Y+ VOID
MIX
0 0 5 5 5 5 5 5 5 0 0
0 5 5 4 4 4 4 4 5 5 0
5 5 4 4 2 2 2 4 4 5 5
5 4 4 3 2 3 2 3 4 4 5
5 4 2 2 3 2 3 2 2 4 5
5 4 2 3 2 1 2 3 2 4 5
5 4 2 2 3 2 3 2 2 4 5
5 4 4 3 2 3 2 3 4 4 5
5 5 4 4 2 2 2 4 4 5 5
0 5 5 4 4 4 4 4 5 5 0
0 0 5 5 5 5 5 5 5 0 0
MESHX 0.0 <<mesh1>> <<mesh2>> <<mesh3>> <<mesh4>> <<mesh5>> <<mesh6>>
<<mesh7>> <<mesh8>> <<mesh9>> <<mesh10>> <<mesh11>>
MESHY 0.0 <<mesh1>> <<mesh2>> <<mesh3>> <<mesh4>> <<mesh5>> <<mesh6>>
<<mesh7>> <<mesh8>> <<mesh9>> <<mesh10>> <<mesh11>>
SPLITX 2 2 2 2 2 2 2 2 2 2 2
SPLITY 2 2 2 2 2 2 2 2 2 2 2
;
TRACK := TRIVAT: GEOM_OSC ::
TITLE 'SMALL CORE 9x9 BENCHMARK'
EDIT 0 MAXR 3000 DUAL 3 1 ;
************************************************************************
* MULTICOMPO IMPORT OF THE REFLECTOR *
************************************************************************
REFL := _REFL ;
ECHO "Multicompo exportation completed" ;
************************************************************************
* CALCULATION OF MODERATOR DATA *
************************************************************************
REAL A := 6.022094E-1 ; ! Avogadro number
REAL P_B10 P_B11 := 0.202 0.798 ;
REAL M_H2O := 15.9994 2.0 1.00794 * + ;
REAL M_B := P_B10 10.012938 * P_B11 11.009276 * + ;
REAL f1 := 1.0 1.0E-06 densB1 * - ;
REAL f2 := 1.0E-06 densB1 * ;
REAL N_O16 := dens_mod f1 A * * M_H2O / ;
REAL N_H1 := 2.0 N_O16 * ;
REAL N_B10 := P_B10 dens_mod f2 A * * * M_B / ;
REAL N_B11 := P_B11 dens_mod f2 A * * * M_B / ;
ECHO "Mix_UOX N_O16=" N_O16 ;
ECHO "Mix_UOX N_H1=" N_H1 ;
ECHO "Mix_UOX B10=" N_B10 ;
ECHO "Mix_UOX B11=" N_B11 ;
************************************************************************
* MPO INTERPOLATION AT 700 ppm BORON *
************************************************************************
ECHO densB1 "ppm boron is" N_B10 "b10 nuclii concentration" ;
MACRO := MCR: CASEA CASEB CASEC CASED :: EDIT 3 MACRO NMIX 5
MPO CASEA output_0 MIX 1 SET 'BURN' 0.0 SET 'C_B10_WATER' <<N_B10>>
ENDMIX
MPO CASEB output_0 MIX 2 SET 'BURN' 0.0 SET 'C_B10_WATER' <<N_B10>>
ENDMIX
MPO CASEC output_0 MIX 3 SET 'BURN' 0.0 SET 'C_B10_WATER' <<N_B10>>
ENDMIX
MPO CASED output_0 MIX 4 SET 'BURN' 0.0 SET 'C_B10_WATER' <<N_B10>>
ENDMIX
;
************************************************************************
* MULTICOMPO INTERPOLATION AT 700 ppm BORON *
************************************************************************
MACRO := NCR: MACRO REFL :: EDIT 2 MACRO NMIX 5
COMPO REFL Radial MIX 5 FROM 2 SET 'C-BORE' <<densB1>> SET 'Method'
'ERM-NEM' SET 'Palier' 'FA3' ENDMIX
;
************************************************************************
* FLUX SOLUTION (RAVIART-THOMAS METHOD) *
************************************************************************
SYS := TRIVAA: MACRO TRACK ;
FLUX := FLUD: SYS TRACK :: EDIT 1 ADI 4 ACCE 5 3 ;
GREP: FLUX :: GETVAL 'K-EFFECTIVE ' 1 >>keff1<< ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.034929 ;
ECHO "test SmallCore_rt_mpo.x2m completed" ;
END: ;
|
