blob: a6df25ed5ddf0d9ca92c14f7728c63b64550ad0c (
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
|
*---
* Raviart-Thomas Reflector DF-RT equivalence prodedure
*---
LINKED_LIST Geom FVol Param MACRO ;
*** STRING nameFeed := "assbly_caseA_mpo_boron.hdf" ;
STRING nameFeed := "assbly_caseA_mpo_draglib_feed.hdf" ;
STRING nameRefl := "Reflector_beavrs_DF-RT_Jeff3p1p1Apolib99_mpo.hdf" ;
STRING Palier := "BEAVRS" ;
STRING htype := "DF-RT" ;
STRING LibType := "JEFF311" ;
STRING NuclData := "DLIB_99" ;
STRING OptSph := "SPH" ;
STRING OptNoal := "ALBE" ;
STRING OptSS := "TONE" ;
MODULE GEO: UTL: MCR: END: ;
HDF5_FILE MpoRefl :: FILE <<nameRefl>> ;
PROCEDURE REFL_MPO ;
!
Geom := GEO: ::
CAR1D 12
X- REFL X+ VOID
MESHX -43.0 -21.5 -5.0 0.0 0.005 2.2225 3.9825 9.6975 10.1975 15.9125 21.5 21.505
35.9125
MIX 1 1 1 2 3 4 3 5 6 5
! fuel gap ss ss+h2o ss h2o ss_shield h2o
7 8
! gap h2o
SPLITX 20 15 20 1 25 24 10 5 6 5
1 5
;
* Define volume fractions in 8 mixtures
* mix=1 : feeding assembly (first fraction is always set to 0)
* mix=2 : fuel-reflector gap
* mix=7 : right-most gap to compute albedo
* mix=8 : residual right reflector
INTEGER nbZones := 8 ;
REAL fvacier := 0.05 ;
REAL fveau := 1.0 fvacier - ;
FVol := UTL: ::
! MIX 1 2 3 4 5 6 7 8
CREA 'H2O' <<nbZones>> = 0. 0. 0. <<fveau>> 1. 0. 1. 1.
CREA 'SS ' <<nbZones>> = 0. 1. 1. <<fvacier>> 0. 1. 0. 0.
;
UTL: FVol :: DUMP ;
* Define feeding assembly parameters (0 to 3 parameters allowed)
Param := UTL: ;
INTEGER nlf := 1 ;
EVALUATE nbZones := nbZones 1 - ; ! decreases nbZones by one
MpoRefl := REFL_MPO Geom FVol Param :: <<nameFeed>> <<nameRefl>> <<Palier>> <<htype>>
<<nlf>> <<LibType>> <<NuclData>> <<OptSph>> <<OptNoal>>
<<OptSS>> <<nbZones>> ;
************************************************************************
* MPO extrapolation of the reflector at 350 ppm Boron
************************************************************************
REAL densB1 := 350.0 ;
MACRO := MCR: MpoRefl :: EDIT 3 MACRO NMIX 1
MPO MpoRefl output_0
MIX 1 FROM 2 SET 'C-BORE' <<densB1>> SET 'Method' <<htype>> SET 'Palier' <<Palier>> ENDMIX
;
UTL: MACRO :: DIR DUMP ;
ECHO "test Reflector_beavrs_DF-RT_Jeff3p1p1Apolib99.x2m completed" ;
END: ;
|
