Home

A USER GUIDE FOR DONJON VERSION5 A. Hébert, D. Sekki, and

1. mSORot 100000 quadr 1800 1 E 5 1 E 5 multi compo file k SEQ_ASCII SCPO MCFD 1 05 gt Nowe N 1 0 DANN FILE MultiCompo geometry construction x GEOM Pgeom x reactor material index x GEOM MATEX x numerical discretization USPLIT GEOM NREFL lt lt nbRef1 gt gt RMIX lt lt mRef11 gt gt lt lt mRef12 gt gt NFUEL lt lt nbFuel gt gt FMIX lt lt mFueli gt gt lt lt mFuel2 gt gt EDIT 1 NGRP 2 MAXR lt lt MaxReg gt gt 163 IGE 344 kh IF Method TRACK MCFD THEN TRIVAT GEOM EDIT 1 MAXR lt lt MaxReg gt gt MCFD lt lt degree gt gt ELSEIF Method PRIM THEN TRACK TRIVAT GEOM EDIT 1 MAXR lt lt MaxReg gt gt PRIM lt lt degree gt gt ELSEIF Method DUAL THEN TRACK ENDIF kh TRIVAT GEOM EDIT 1 MAXR lt lt MaxReg gt gt DUAL lt lt degree gt gt lt lt quadr gt gt macrolib for reflector and devices k LCPO SCPO MACRO1 x NCR LCPO EDIT 1 MACRO NMIX lt lt nbMix gt gt COMPO LCPO MODE MIX lt lt mRef11 gt gt SET MTYPE CELL20 ENDMIX MIX lt lt mRef12 gt gt SET MTYPE CELL18 ENDMIX COMPO LCPO FUEL MIX lt lt mZCRin gt gt SET FTYPE ZCU SET RDTPOS 1 SET RDDPOS ENDMIX MIX lt lt mZCRot gt gt SET FTYPE ZCU SET RDTPOS 1 SET RDD
2. 2nd initial size of the hypercube for MAP method default 0 1 3rd Eext limit for external convergence default 1074 4th Eint limit for internal convergence default 1074 5th Equad limit for convergence of the quadratic constraint default 1074 The other value of the record are not used and set to 0 0 The optional array DLEAK STATE contains integer values related to the definition of mixture and group indices in module DLEAK e The number of energy groups in macrolib S The number of mixtures in macrolib S The type of leakage parameters S3 where SI 1 use diffusion coefficients 2 use P weighted macroscopic total cross sections The type of control variables S where SI 1 use leakage parameter itself 2 usea correction factor The minimum group index S with 1 lt S lt Sf e The maximum group index Sf with S lt SZ lt SJ The minimum mixture index S with 1 lt S lt SZ e The maximum mixture index Sg with S lt Sf lt S 8 8 1 The sub directory OLD VALUE in optimize Table 104 Main records and sub directories in OLD VALUE Condition UnitsComment VAR VALUE uu Noar The values of the decision variables of the last valid iteration FOBJ CST VAL Nest 1 The values of the objective function and the contraints of the last valid iteration GRADIENT uuu Novar Nest 1 The gradients of the objective function and the constraints of the last valid
3. AXIS youu uuu The number used to identify the rod mouve ment direction 1 along x axis 2 along y axis 3 along z axis FROM ou The number used to identify the side of ge ometry from which the controller rod is in serted into the reactor core along its direction of mouvement l if a rod is inserted from the highest position e g from the top 1 if a rod is inserted from the lowest position e g from the bottom LENGTH yuu The initial and final position coordinates of the full inserted rod along its direction of movement The rod length is the distance be tween these two coordinates CORE LIMITS The initial and final position coordinates of the full core along each Cartesian direction MAX POSuuuuu The limiting 3D Cartesian coordinates of the full inserted rod This data is given for each part of the rod continued on next page IGE 344 135 Records in DEV ROD sub directories continued from last page Condition Units Comment ROD MIX youu The rod type mixture indices The first num ber corresponds to the inserted rod position and the second to the withdrawn rod position This data is given for each part of the rod LEVEL uuuuuu The actual insertion level of the controller rod This value must be between 0 0 for the full withdrawn rod and 1 0 for the full inserted rod SPEED uuu uu The speed of rod movement insertion or ex traction TIME uuuuuuu 3 The time for the full rod insertion or extrac tion
4. means change line Generated by GenPMAXS V6 1 1 Block BRANCHES information optional This blocks identifies the state variables used for the branching and the correponding states for all branches STA_VAR 4 CR DC PC TF BRANCHES 1 2 6 18 54 RE 1 0 00000 0 71100 1000 00000 900 00000 CR 1 1 00000 0 71100 1000 00000 900 00000 CR 2 2 00000 0 71100 1000 00000 900 00000 DC 1 0 00000 0 66100 1000 00000 900 00000 DC 2 0 00000 0 75200 1000 00000 900 00000 DC 3 1 00000 0 66100 1000 00000 900 00000 DC 4 1 00000 0 75200 1000 00000 900 00000 DC 5 2 00000 0 66100 1000 00000 900 00000 DC 6 2 00000 0 75200 1000 00000 900 00000 PC 1 0 00000 0 66100 0 00000 900 00000 In this example the PMAXS cross sections depend on 4 state variables CR DC PC and TF There are 2 branches for control rods 6 for coolant density 18 for boron concentration and 54 for fuel temperature Block BURNUP information optional It contains the number of burnup sets and burnup points Each cross sections will be repeated for each burnup points IGE 344 99 BURNUPS 1 1 35 0 00000 0 00900 0 01900 0 07500 0 15000 0 50000 1 00000 2 00000 3 00000 4 00000 6 00000 8 00000 10 0000 12 0000 14 0000 16 0000 18 0000 20 0000 24 0000 28 0000 32 0000 36 0000 40 0000 44 0000 48 0000 52 0000 56 0000 60 0000 64 0000 68 0000 72 0000 76 0000 80 0000 84 0000 86 0000 In the above example one can observe a set of 35 burnup points from 0 to 86 GWd t Block XS SET identification mand
5. FULL HALF QUARTER CHECKER CHECKER 1 2 CHECKER 1 4 54 integer index used to control the printing on screen 0 for no print 1 for minimum printing 2 modified fuel indices and coolant densities are printed per bundle over each channel 3 modified fuel indices are printed per each radial plane for larger values of iprint everything will be printed keyword used to specify mixF integer fuel type mixture number of the non perturbed fuel cell This number must be specified for each fuel type as been recorded in the MATEX object see Section 3 2 1 keyword used to specify mixV integer new mixture number assigned to the voided fuel cell Note that this number must be specified for each fuel type and it must be different from any other reactor material mixtures keyword used to specify PNAME This information is required only for the interpola tion of fuel properties using the NCR module character 12 user defined identification name of local parameter associated with the relative coolant density The recommended name is D COOL This parameter name and the unperturbed densities values should be previously recorded in the FMAP object see Section 3 1 2 The same PNAME will be set for the coolant density in the perturbed FMAPV but the actual values of coolant densities throughout the core will be reordered by the CVR module according to the specified voiding pattern keyword used to specify the value dcoolV
6. NAMDET NHEX nhex HEX ihex POSITION pos RESP Tep ENDN 24 keyword to specify the detector name character 12 name of the detector The different names in alphabetical order must fit their usual numbering in the core Ex PLATNO1 CHIONO1C keyword to set the number of hexagons where the detector is placed number of hexagons keyword to set the hexagon numbers corresponding to the detector position array containing the hexagon numbers where the detector is present as ordered in the geometry definition keyword to specify the detector coordinates array containing the positions of the specified detector The positions must be read as X X Y Y Z Z For 2 D geometry Z coordinates must be 0 0 and a value greater than 1 0 For hexagonal geometry only Z coordinates are used in 3 D representation keyword to specify the detector initial responses array containing the initial responses of the detector To use the current detector models in DONJON responses are given as e For vanadium detectors current response last response e For platinum detectors current response reference flux last detector slow re sponses e For ion chamber detectors current logarithmic response current log rate re sponse reference flux keyword to specify the end of the detector informations IGE 344 25 3 6 The LZC module The LZC module is used for the modeling of liquid zone controllers which are no
7. ROD P0Sovuuu The actual 3D Cartesian coordinates of the rod This data is given for each part of the rod 8 3 4 The ROD GROUP sub directories Inside each ROD GROUP sub directory the following records will be found Table 90 Records in ROD GROUP sub directories Type Condition Units Comment GROUP ID uuu The identification number of the rod group NUM ROD uuu The total number Nya of rod devices in the group ROD ID uuuuu An array of identification numbers of rods which belong to the same group 8 3 5 The DEV LZC sub directories Inside each DEV LZC sub directory the following records will be found IGE 344 136 Table 91 Records in DEV LZC sub directories Type Condition Units Comment LZC IDuuuuuu The identification number of the liquid zone controller MAX P0Suuuuu The limiting 3D Cartesian coordinates of the whole liquid controller including its empty and full parts AXIS youu uuu The number used to identify the water filling direction 1 along x axis 2 along y axis 3 along z axis HEIGHT yoy uuu The water height of the full filled controller along its direction of filling LEVEL uuuuuu The actual water level of the liquid controller device This value must be between 0 0 for the empty state and 1 0 for the full filled state EMPTY POSLuu The actual 3D Cartesian coordinates of the empty part of liquid contoller FULL POS juju The actual 3D Cartesian coordinates of the full part of liquid
8. descthm structure describing the input data to the THM module 5 1 1 Input data to the THM module Table 71 Structure descthm EDIT iprint RELAX relax TIME caltype timestep timeiter time FPUISS fract CRITFL cflux continued on next page IGE 344 108 Structure descthm continued from last page CWSECT sect flow SPEED velocity ASSMB sass nbf nbg INLET poutlet tinlet RADIUS rl 12 r3 r4 POROS poros PUFR pufr CONDF ncond kcond k k 0 ncond INV inv ref unit CONDC ncond kcond k k 0 ncond unit HGAP hgap HCONV hconv TEFF wtef CONV maxitl maxit2 maxit3 ermaxt ermaxc RODMESH nb1 nb2 FORCEAVE BOWR SAHA SET PARAM PNAME pvalue 3 where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing larger values produce increasing amounts of output RELAX keyword used to set the relaxation parameter relax relax relaxation parameter selected in the interval 0 lt relax lt 1 and used to update the fuel average and surface temperature coolant temperature and coolant density The updated value is taken equal to 1 relax times the previous iteration value plus relax times the actual iteration value The default value is relax 1 TIME keyword used to specify the type of calculation steady state or transient performed by the THM module and the temporal pa
9. integer total number of combustion zones This value must be greater than or equal to 1 and less than or equal to the total number of reactor channels keyword used to specify icz integer array of combustion zone indices specified for every channel A reactor channel can belong to only one combustion zone however a combustion zone can be specified for several channels keyword used to indicate that the total number of combustion zones equals to the number of reactor channels In this particular case each channel will have a unique combustion zone number Hence an explicit specification of the combustion zone in dices can be omitted Non number of fuel channels in the radial plane Np number of fuel bundles or assembly subdivisions in the axial plane after NCOMB keyword to specify that one combustion zone per assembly is to be defined keyword used to specify a basic assembly layout for the SIM PWR refuelling module see Section 3 14 number of assemblies along the X axis Typical values are 15 or 17 number of assemblies along the Y axis character 3 identification name corresponding to the basic naval coordinate position of an assembly naval i is the concatenation of a letter generally chosen between A and T and of an integer generally chosen between 01 and 17 An assembly may occupies four positions in the fuel map in order to be represented by four radial burnups In this case the same naval coordinate value
10. ivarty index of the combustion zone for differentiation of cross section information IGE 344 SET DELTA ADD LINEAR CUBIC PARKEY vall val2 MAP REF valref SAMEASREF ENDREF MICRO ALL ONLY HISO conc 70 keyword used to indicate a simple interpolation at vall or an averaging between vall and val2 The result Gef is also used as the reference value when the ADD is used Note see at the ending note of this section for a detailed description and examples keyword used to indicate a delta sigma calculation between val2 and vall i e Aorep Oval2 Oval IS computed Note see at the ending note of this section for a detailed description and examples keyword used to indicate a delta sigma calculation between val2 and vall is added to the reference value i e Ao Oval2 Oval is used as contribution Gref Ao or Aref Ao is returned Note see at the ending note of this section for a detailed description and examples keyword indicating that interpolation of the MULTICOMPO for parameter PARKEY uses linear Lagrange polynomials It is possible to set different interpolation modes to different parameters By default the interpolation mode is set in Sect 4 2 1 keyword indicating that interpolation of the MULTICOMPO for parameter PARKEY uses the Ceschino method with cubic Hermite polynomials as presented in Ref 17 By default the interpolation mode is set in Sect 4 2 1 character
11. 1 for minimum editing default value 2 only channel powers in radial plane are printed 3 only bundle powers per each radial plane are printed 10 only bundle powers per each channel are printed Any combination of the values 2 3 and 10 is possible for example 5 243 Note that any other value of iprint behaves as the first lower possible value for example 7 gives the same output as 5 Moreover channel and bundle powers can be printed only if the FMAP object was provided in the calling specification PTOT keyword used to specify the input of power By default a power is recovered from the KINET object IGE 344 power FSTH fsth INIT P NEW PRINT MAP DISTR FLUX RATIO POWER ALL NORM BUND 38 real total reactor power given in MW This value must correspond to the reactor nominal conditions keyword to specify the thermal to fission power ratio thermal to fission power ratio By default this value is not used and the total power is the one given after the PTOT keyword keyword used to save the actual power distribution in the BUND PW INI record of the fuel map object FMAP It is used by the AFM module to apply power feedback during a fast transient using the initial power distribution instead of the actual power keyword used to indicate that a new total reactor power is to be recomputed based on the previously calculated flux normalization factor The flux and power distributions over
12. ENDROD ROD 23 ROD NAME ZCR12A LEVEL 0 0 AXIS Y FROM H MAXPOS 205 0 227 0 0 0 260 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 24 ROD NAME ZCR12B LEVEL 0 0 AXIS Y FROM H MAXPOS 205 0 227 0 260 0 520 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD IGE 344 174 ROD 25 ROD NAME ZCR13A LEVEL 0 0 AXIS Y FROM H MAXPOS 249 0 271 0 0 0 260 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 26 ROD NAME ZCR13B LEVEL 0 0 AXIS Y FROM H MAXPOS 249 0 271 0 260 0 520 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 27 ROD NAME ZCR14A LEVEL 0 0 AXIS Y FROM H MAXPOS 293 0 315 0 0 0 260 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 28 ROD NAME ZCR14B LEVEL 0 0 AXIS Y FROM H MAXPOS 293 0 315 0 260 0 520 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 29 ROD NAME ZCR15A LEVEL 0 0 AXIS Y FROM H MAXPOS 337 0 359 0 0 0 260 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 30 ROD NAME ZCR15B LEVEL 0 0 AXIS Y FROM H MAXPOS 337 0 359 0 260 0 520 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD oe SOR S ROD 31 ROD NAME SORO1 LEVEL 0 0 AXIS Y FROM H MAXPOS 117 0 139 0 53 75 466 25 49 53 99 06 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD IGE 344
13. FLMAP MATEX RESINI MATEX COMPO descresinil FLMAP RESINI FLMAP FLMAP2 descresini2 where IGE 344 9 FLMAP character 12 name of the RESINI object that will contain the fuel lattice information If FLMAP appears on both LHS and RHS it will be updated otherwise it is created MATEX character 12 name of the MATEX object specified in the modification mode MATEX is required only when FLMAP is created COMPO character 12 name of the MULTICOMPO data structure L_COMPO signature where the detailed subregion geometry at assembly level is stored FLMAP2 character 12 name of the RESINI object that contains the fuel lattice information to recover from descresinil structure describing the main input data to the RESINI module Note that this input data is mandatory and must be specified only when FLMAP is created descresini2 structure describing the input data for global and local parameters This data is permitted to be modified in the subsequent calls to the RESINI module 3 1 1 Main input data to the RESINI module Note that the input order must be respected Table 3 Structure descresinil EDIT iprint SPLIT NAP GEO descgeo NAP descnap ASSEMBLY na nax nay A ZONE iza i i 1 nch ASBLY AXNAME XNAMEA i i 1 nax AYNAME YNAMEA i i 1 nay NXNAME XNAME i i 1 nx NYNAME YNAME i i 1 ny NCOMB ncomb B ZONE icz i
14. HEXZ 22 HGAP 108 110 hgap 108 110 HISO 65 66 69 71 80 82 History 60 history 140 141 HISTORY 5 56 61 HISTORY CASE identification mandatory 99 HISTORY DATA 147 HistXSM 61 HM_Dens 153 hm_dens 104 HOUR 43 46 47 HST 4 5 56 59 141 hstbrn 56 58 hstdim 56 57 hstpar 56 59 hx 45 46 hy 45 46 I BURNUP 65 66 1135 101 102 ialch 12 14 ialch i 14 ibm1 112 113 ibm2 112 113 ibun 56 57 ibun 61 icha 56 57 icha 61 IGE 344 icz 9 11 id 20 21 26 28 33 34 idfuel 56 57 idfuel 61 ifuel 46 IFX 121 ifs 121 122 igc 94 igrp 21 28 33 34 ihex 23 24 ilg 95 ilgrp 29 30 ilzc 29 30 IMET 89 91 imet 89 91 imix 43 44 69 80 88 89 101 102 123 imixold 69 80 in 18 inam 101 102 indcycle 45 46 indcyclel 46 indcycle2 46 INF 100 101 NFO 22 23 NFOR 88 89 NIT 37 38 50 88 89 NLET 108 109 NN STEP EPS 115 116 NSR 33 34 NST BURN 12 65 69 80 81 NST BVAL 12 NV 108 110 inv 108 110 INVCONST 22 23 iph 101 iphase 100 101 iphase 1 101 iprint 9 12 15 19 26 29 31 33 35 37 39 40 43 45 46 49 50 53 54 57 64 67 68 78 79 88 89 94 100 101 107 108 112 115 118 121 123 iprt 22 51 irgrp 29 30 irod 29 30 ishift 12 13 ISOTOPES 101 102 ispx 123 ispy 123 IUPS 101 106 jups 101 106 ivarty 65 66 69 80 81 iza 9 10
15. MICRO HISO conc ALL ENDMIX Table 51 Structure descdata2 MIX mix NAMDIR DERIV UPS TIMAV BURN INST BURN AVG EX BURN ivarty MICRO HISO conc y ALL ENDMIX where MIX keyword used to set the material mixture mix mix integer identifier for the material mixture that will be included in the MACROLIB The maximum number of identifiers permitted is nmix and the maximum value that mix may have is nmix Note that if TABLE is the selected option then mix identifies the fuel type as defined in the reactor geometry NAMDIR character 12 directory name in the CPO object from which the nuclear properties for material mixture mix are to be recovered DERIV keyword used to compute the derivative of the MACROLIB information with respect to burn or burn value By default the MACROLIB information is not differentiated UPS keyword used to compute properties with no up scattering contribution TIMAV BURN keyword used to compute time averaged cross section information This option is available only if TABLE is the selected option By default the type of calculation TIMAV BURN or INST BURN is recovered from the FMAP object INST BURN keyword used to compute cross section information at specific bundle burnups This option is available only if TABLE is the selected option By default the type of calcu lation TIMAV BURN or INST BURN is recovered from the FMAP object AVG EX
16. They are many options on how to use the module AFM for its different purposes A compact summary is presented on Tab 65 The Rozon correlation for fuel temperature as a function of bundle power is True La 0 476 P 2 267 P x 1074 where P is in kW and temperatures are in Kelvin IGE 344 MAP with local parameters MAP local ters without parame Table 65 AFM options summary REFT REFT TFUEL TCOOL sat R REFT TFUEL TCOOL di REFT REFT TFUEL TCOOL i REFT TFUEL TCOOL Nominal values Nominal values except for specified parameters Nominal values except for TFUEL parameter which is computed according to the Rozon cor relation using nominal power Same as above except for specified parameters which will have a constant value Nominal values except for local parameters in cluded in MAP Same as above except for specified parameters which will have a constant value Nominal values except for TFUEL parameter which is computed according to the Rozon cor relation if power distribution is available Same as above except for specified parameters which will have a constant value 92 IGE 344 93 4 5 The T16CPO module The WIMS AECL Tape16 file is a FORTRAN seguential binary file which is used to transfer the results of a WIMS AECL calculation to other applications The explicit contents of this file may vary from application to application since the output of most records to thi
17. fuel map macrolib x MACFL x CRE LFUEL1 LFUEL2 FMAP EDIT O READ TABLE LFUEL1 MIX lt lt mFuel1 gt gt lt lt NamFuel1 gt gt ENDMIX TABLE LFUEL2 MIX lt lt mFuel2 gt gt lt lt NamFuel2 gt gt ENDMIX extended macrolib k MACRO2 MACFL kh MATEX MACINI MATEX MACRO1 MACFL DELETE MACFL complete macrolib x MACRO MATEX MACR02 dini NEWMAC MATEX MACRO2 DEVICE DELETE MACRO2 numerical solution Sini SYSTEM TRIVAA MACRO TRACK EDIT O EDIT 1 ENDMIX ENDMIX ENDMIX ENDMIX EDITO EDIT O 161 IGE 344 162 MACRO DELETE MACRO IF iter 1 THEN FLUX FLUD SYSTEM TRACK EDIT O ACCE 3 3 ADI 4 EXTE 1000 lt lt Precf gt gt THER 1000 ELSE FLUX FLUD FLUX SYSTEM TRACK EDIT O ACCE 3 3 ADI 4 EXTE 1000 lt lt Precf gt gt THER 1000 ENDIF SYSTEM DELETE SYSTEM flux and power POWER FLPOW FMAP FLUX TRACK MATEX EDIT 0 PTOT lt lt Power gt gt burnups integration limits B FMAP TAVG FMAP POWER EDIT 0 AX SHAPE RELAX 0 5 B EXIT POWER DELETE POWER current parameters GREP FLUX GETVAL K EFFECTIVE 1 gt gt Keff lt lt GREP FMAP GETVAL EPS AX 1 gt gt Eps lt lt ECHO Iteration No iter ECHO AXIAL SHAPE ERROR Eps ECHO RESULTING K EFF Keff ENDWHILE edit resulting fluxes and powers po
18. 2563 4136 3173 3801 3614 4046 3705 3990 3666 3910 3826 3889 LINKED LIST FMAP END RESINI gt 2 EDIT 2 REF SHIFT 8 61 79 27 72 73 47 69 41 12 48 31 96 02 17 07 1095 41 1656 2556 2898 3111 3953 3764 4119 3960 4037 4115 3915 4232 34 62 28 80 67 97 31 65 40 83 26 46 3890 81 1137 1756 2822 3261 3212 3813 4030 4241 4380 92 62 28 49 21 93 43 15 19 182 IGE 344 END QUIT 183 IGE 344 184 10 11 12 13 14 15 16 17 18 19 References A H bert Applied Reactor Physics Presses Internationales Polytechnique ISBN 978 2 553 01436 9 424 p Montr al 2009 A H bert Development Procedures for Version4 of Reactor Physics Codes Report IGE 287 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2006 A H bert and R Roy The GAN Generalized Driver Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2000 A H bert and R Roy The Ganlib5 kernel guide 64 bit clean version Report IGE 332 Ecole Polytechnique de Montr al January 2013 G Marleau A H bert and R Roy A User Guide for DRAGON Version4 Report IGE 294 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2007 A H bert A User Gu
19. 26 27 MAX MIX GEO 123 124 MAXIMIZE 115 116 maxit1 108 111 maxit2 108 111 maxit3 108 111 MAXPOS 20 21 26 27 MAXR 15 16 maxreg 15 16 MB 95 MCC 32 MCC 4 31 32 MCO 100 101 MCO 102 103 153 MCR 88 90 MD 95 96 MEMORY 78 79 METH 122 METHOD 115 MIC 100 101 MICLIB 42 MICLIB2 42 MICLIB3 42 MICRO 43 65 70 78 81 MICROLIB 42 MICROLIB 49 MICROLIB 4 MICROLIB_XS 156 157 MINIMIZE 115 116 MINUTE 43 46 47 MIX 10 65 69 80 94 95 101 102 192 mix 65 66 MIX FUEL 53 54 MIX VOID 53 54 mix1 20 21 mix2 20 21 MIXASS 123 mixdir 101 102 mixE 26 27 mixF 26 27 53 54 mixf 15 16 MIXMAX 112 113 MIXMIN 112 113 MIXNAM 94 96 mixr 15 16 mixV 53 54 MLIB 67 68 78 79 MLIB2 67 78 mmix 88 89 MODNAME 6 MODULE 6 module 6 7 MOVDEV 4 21 33 133 MOVE 18 19 MP 95 MT 95 96 MTMD 94 96 MTS 96 MULTICOMPO 5 9 120 121 MW day per tonne 13 N2N 82 N3N 82 NAN 83 NA 83 na 9 10 12 namburn 67 68 78 79 NAMCHA 43 44 NAMDB 88 89 NAMDET 23 24 NAMDIR 65 67 68 NAMDIR 101 102 NAMDPL 82 NAME 20 NAME 23 24 name_detect 138 name_detect 138 name_type 137 138 name_type 137 namedir 121 123 NAMPAR 59 82 83 NAMPER 94 96 NAMTYP 22 23 52 88 89 NAP 120 IGE 344 NAP 4 9 10 16 120 124 naval 9 11 nax 9 10 nay 9
20. B MIX 1 SET A lt lt va0 gt gt ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF SCR SAP FMAP ENDMIX Table 62 SCR inputs for TA cases IGE 344 88 4 4 The AFM module The AFM module is used to create an extended MACROLIB containing set of interpolated nuclear properties from a feedback model database The DATABASE information are obtained by previous DRAGON calculations using module CFC 1 There are two possible utilizations e Construction of an extended MACROLIB for fuel properties directly from DATABASE information with respect to local parameters contained in the fuel map object or directly input e Construction of an extended MACROLIB containing only one set of cross sections derivated from the DATABASE information Properties can be obtained for fuel or reflector The calling specifications are Table 63 Structure AFM MACRO AFM MACRO DBASE MAPFL descafm where MACRO character 12 name of the extended MACROLIB The MACROLIB can be in modification mode DBASE character 12 name of the DATABASE object containing fuel properties with respect to local parameters MAPFL character 12 name of the MAP object containing fuel regions description and burnup informations This file is only required when a MACRO is created for fuel area descafm structure containing the data to module AFM 4 4 1 Input data to the AFM module Table 64
21. FMAP character 12 name of a FMAP object that will be updated by the TINST module The FMAP object must contain the instantaneous burnups for each fuel bundle and the weight of each fuel mixture POWER character 12 name of a POWER object containing the channel and bundle powers previously computed by the FLPOW module The channel and bundle powers are used by the TINST module to compute the new burn up of each bundle If bundle powers are previously specified with the module RESINI you can refuel your core without a POWER object MICLIB3 character 12 name of a LIBRARY object that will be created by the TINST module This MICROLIB contains the fuel properties after refueling when keyword MICRO is used in desctinst MICLIB2 character 12 name of a LIBRARY object that will be read by the TINST module This must be a fuel map LIBRARY created either created by the NCR or the EVO module MICLIB character 12 name of a LIBRARY object that will be read by the TINST module This MICROLIB contains the new fuel properties that should be used for the refueling desctinst structure describing the input data to the TINST module 3 13 1 Input data to the TINST module Note that the input order must be respected IGE 344 43 Table 31 Structure desctinst EDIT iprint BURN STEP rburn TIME rtime DAY HOUR MINUTE SECOND REFUEL MICRO CHAN NAMCHA nsh NEWFUEL CHAN NAMCHA nsh SOME imix i i
22. In addition the POWER object will store several parameters that can be used as power and criticity constraints for the optimization and fuel management purposes namely the maximum channel and bundle powers the channel and bundle power form factors the ef fective multiplication factor recovered from the FLUX or KINET data structure The FLPOW module specifications are Table 26 Structure FLPOW POWER NRMFLUX FMAP FLPOW POWOLD FMAP FLUX KINET TRACK MATEX descflpow POWER FLPOW POWOLD FLUX KINET TRACK MACRO descflpow where POWER character 12 name of the POWER object that will be created by the module It will contain the information related to the reactor fluxes and powers NRMFLUX character 12 name of the FLUX object in creation mode According to the chosen option this object contains either the fluxes normalized to the given total reactor power or the fluxes per bundle Is it useful if you want to compute the detectors readings with the DETECT module POWOLD character 12 name of the read only POWER object It must contain the previously computed flux normalization factor which corresponds to the reactor nominal or equi librium conditions FMAP character 12 name of the FMAP object containing the fuel lattice specification When FMAP is specified on the RHS the fluxes and powers calculations are performed over IGE 344 37 the fuel lattice as well as over the whole r
23. POWER FLPOW FMAP FLUX TRACK MATEX EDIT lt lt iEdit gt gt PTOT lt lt Power gt gt Haz last parameters Hu GREP FLUX GETVAL K EFFECTIVE 1 gt gt Keff lt lt GREP FMAP GETVAL EPS AX 1 gt gt Eps lt lt GREP FMAP GETVAL B EXIT 1 gt gt Bexit lt lt ECHO Number of Iterations iter ECHO AXIAL SHAPE ERROR Eps p ECHO CORE AVERAGE EXIT BURNUP Bexit ECHO RESULTING K EFFECTIVE Keff assertS FLUX K EFFECTIVE 1 1 050128 END QUIT IGE 344 9 2 Example2 Multicompo related example Input data for test case Example2 x2m DR KK ARE FKK ARE K AE AE AA A FK FK ARE AE AE FK K AAA AA FK o FK K FK K FK FK AA KK K K K K K K OK Purpose XK XA XA A dd OK K K K K ok PROCEDURE MODULE LINKED_LIST x Input file Author s Example2 x2m Test for non regression using DONJON 4 D Sekki 2007 11 kk X assertS Pgeom Pfmap Pburn Pdevc DELETE GREP END NCR MACINI FLUD USPLIT DSET TAVG FLPOW TRIVAT TRIVAA NEWMAC GEOM TRACK MATEX FMAP FLUX POWER MACFL DEVICE MACRO1 MACRO2 MACRO SYSTEM LCPO variables INTEGER INTEGER INTEGER INTEGER INTEGER INTEGER INTEGER STRING INTEGER INTEGER REAL REAL REAL REAL fos nbMix nbRefl mFuel1 mRef11 mZCRin mSORin MaxReg Method degree iter iEdit Power epsil Precf Eps Keff Bexit 8 nbFuel mFuel2 mRef12 mZCRot
24. SET AVGB avburn SET FUEL ifuel FROM hcold2 AT asmb2 BURN indcycle burncycle B DIST AX asmb1 SET axn i i 1 nb FROM hcold2 AT asmb2 BURN indcycle burncycle I BURN STEP rburn TIME rtime DAY HOUR MINUTE SECOND ENDCYCLE COMPARE hcl BURN indcyclel burncyclel hc2 BURN indcycle2 burncycle2 DIST BURN gt gt epsburn lt lt DIST POWR gt gt epspowr lt lt SET PARAM PNAME pvalue where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value for larger values of iprint everything will be printed CYCLE keyword defining operations based on the actual fuel cycle hcnew character 12 identification name of the specific fuel cycle FROM keyword defining the previous fuel cycle in case that some information needs to be transmitted to the actual fuel cycle heold character 12 identification name of the previous fuel cycle BURN keyword defining the burnup at which the assembly is recycled in the previous fuel cycle By default the last burnup step is used indcycle integer index of the burnup step in the previous fuel cycle burncycle real value of the burnup in the previous fuel cycle MAP keyword defining the assembly layout in naval coordinate positions in the actual fuel cycle Here lx and ly values are those defined in the fue
25. XK XA XA XX XX X A dd K OK K K KK ok PARAMETER GEOM LINKED_LIST GEOM MODULE END GEO GEOM GEO CAR3D 24 24 12 167 IGE 344 Figure 13 Geometry definition plane 1 EDIT 1 X VOID X VOID Y VOID Y VOID Z VOID Z VOID MIX PLANE 1 k 123456789012345678 000000003333 333300000000 000000333333 333333000000 000003333333 333333300000 000333311111 111113333000 111111133000 1 111 111 11330 0 111111113330 111111111330 111111111333 1411144141113323 111111111333 MOA 331111111 31111 111 1 311111111 0 0 03 33 0 0 0 033111111111 333111111111 333111111111 333111111111 168 IGE 344 111111111333 333111111111 K L M N 0 P Q R S 111111111333 LAZ 1 11 1 398 393 111111111333 11 11 11 111333 111111111330 111111113330 111111113300 111111133000 111113333000 333111111111 333111111111 333111111111 33 911 111111 1 033111111111 033311111111 003311111111 000331111111 000333311111 000003333333 333333300000 000000333333 333333000000 000000003333 333300000000 PLANE 2 SAME 1 PLANE 3 SAME 1 PLANE 4 SAME 1 PLANE 5 123456789 012345678 000000004444 444400000000 000000444444 444444000000 000004444444 444444400000 000444422222 222224444000 000442222222 222222244000 004422222222 222222224400 044422222222 222222224440 044222222222 222222222440 444222222222 222222222444 444222222222 222222222444 444
26. XX XX XX A dd 2K K K KK ok PARAMETER DEVICE MATEX LINKED_LIST DEVICE MATEX MODULE END DEVINI INTEGER mZCRin mZCRout mSORin mSORout Read arguments 1 gt gt mZCRin lt lt gt gt mZCRout lt lt gt gt mSORin lt lt gt gt mSORout lt lt 169 IGE 344 170 visi kai 15 eneen MZA L EET m E T E bi CER Li Lo i cc 4 13 PA O IL pizzo sissi ie 12 Cyne Stree IA Meee S N 0 a 0 Q N SOR Shutoff Units 26 LISU Liquid Injection Shutdown Nozzles 6 ZCU Zone Control Units 15 VFD Vertical Flux Detector Units 28 Possible CAU Control Absorber Units 8 Sites 6 FCBR Fission Chamber Units 3 VP View Port 1 aa Fuel Bundle Positions Figure 14 Top View of ACR Benchmark Core Model DEVICE MATEX DEVINI MATEX EDIT 1 NUM ROD 56 FADE x ZCR x ROD 1 ROD NAME ZCRO1A LEVEL 0 0 AXIS Y FROM H MAXPOS 161 0 183 0 0 0 260 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 2 IGE 344 171 ROD NAME ZCRO1B LEVEL 0 0 AXIS Y FROM H MAXPOS 161 0 183 0 260 0 520 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 3 ROD NAME ZCRO2A LEVEL 0 0 AXIS Y FROM H MAXPOS 205 0 227 0 0 0 260 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 4 ROD NAME ZCRO2B LEVEL 0 0 AXIS Y FROM H MAXPOS 205 0 227 0 260 0 5
27. for more details on TRIVAC modules specification refer to its user guide Because the code execution is controlled by the GAN generalized driver it is also possible to use its utility modules A brief description of each module that can be executed using DONJON is given below A short description of each data structure that can be used in DONJON is given in Section 2 2 e The following DRAGON modules can be executed using DONJON GEO module used to create or modify a reactor geometry The spatial locations of the reactor material mixtures must also be defined using the GEO module Only 3 D Cartesian reactor geometries are allowed with DONJON MAC module used to create or modify a MACROLIB containing the material properties by directly specifying the group ordered macroscopic cross sections for each selected ma terial mixture e The following TRIVAC modules can be executed using DONJON TRIVAT module used to perform a 3 D numerical discretization or tracking of the reactor geometry TRIVAA module used to compute the set of system matrices with respect to the previously obtained tracking information FLUD module used to compute the numerical solution to an eigenvalue problem correspond ing to a previously obtained set of system matrices e The following are short descriptions of utility modules that can be executed using DONJON UTL module used to perform several utility actions on a data structure DEL
28. izae 82 HWSHHHHKHHHHEHH 191 JOB_OPT 101 105 JOB_TIT 101 105 jobtit 101 105 ke 51 52 kcond 108 110 KELVIN 110 kg 14 KINET 36 37 KINET 36 37 KINSOL 36 37 L_POWER 127 ladf 103 105 ladf Ixes Ided Ij1f Ichi Ichd linv Idet lyld Icdf Igff lbet lamb Idec 101 lamb 105 lbet 105 ledf 105 Ichi 105 Ichid 105 Idec 105 Ided 105 Idet 105 LEAK 67 68 78 79 LEMKE 115 length 10 LEVEL 20 26 27 29 30 33 34 LFLAG 153 Igff 105 LIBRARY 42 LIBRARY 42 49 LINEAR 67 70 78 81 LINKED LIST 6 7 linv 105 LIST 94 95 1j1 105 LOCAL 12 13 57 58 Ix 9 11 45 46 Ixes 105 ly 9 11 45 46 lyld 105 LZC 26 27 29 30 LZC GR 26 LZC GROUP 29 30 136 Izc group 26 28 LZC ID 28 LZC 4 18 25 30 133 MAC 3 5 17 MACD 22 MACEL 17 63 MACGED 123 MACINI 4 5 17 35 40 IGE 344 MACRES 120 MACRO 88 MACRO 17 36 37 63 64 88 112 113 MACRO 67 68 78 79 MACRO2 17 35 MACRO3 35 MACROLIB 71 83 macrolib 132 MACROLIB 1 3 5 11 17 35 37 40 63 66 88 89 120 MACROLIB_XS 156 MAP 31 MAP 120 122 MAP 37 38 45 46 69 70 80 81 88 89 115 MAP 56 61 88 91 107 120 122 124 MAPO 60 MAPFL 67 70 73 74 78 81 85 86 88 89 107 109 MATEX 8 9 15 17 18 25 35 37 120 MATEX 15 matex 131 132 MATEX 5 9 14 15 17 18 25 35 37 54 120 MAX FULL
29. nunber of evolution 3 nglo nunber of global parameters 1 nloc nunber of local parameters 2 bunl bundle length in cm 49 53 cm goose x x Xx Xx INTEGER ncha nbun nevo nglo nloc 93312 REAL bunl 49 53 Esla Initialize History using MAPO si Reseau MAPO History HST Reseau DIMENSIONS GLOBAL lt lt nglo gt gt LOCAL lt lt nloc gt gt BUNDLES lt lt nbun gt gt lt lt bunl gt gt CHANNELS lt lt ncha gt gt Here the HISTORY data structure will be stored in the XSM file History One global and two local parameters are considered No information about the name or the value of the global and local parameters will be available This initialization procedure stores information only on the main level of the HISTORY data structure if the MAP data structure is not available In this case the HISTORY is updated using a MAP data structure in sequential ASCII file MAPO The number of channels and bundles per channel are stored and compared with the same information in the MAP structure For each bundle in the MAP cell type and fuel type directories are constructed The bundle powers and burnups available in MAP are used to generate the power rates in kW kg and the depletion time in days required to reach the specified burnups These values are stored in the HISTORY in the PARAMBURNTAR record The fuel mass is mandatory for such calculation thus the fuel weight is recovered from the MAP If the HISTO
30. see Section 4 2 2 MACRO keyword indicating that MLIB is a MACROLIB default option MICRO keyword indicating that MLIB is a MICROLIB Object MLIB contains an embedded MACROLIB but the CPU time required to obtain it is longer LINEAR keyword indicating that interpolation of the MULTICOMPO uses linear Lagrange poly nomials CUBIC keyword indicating that interpolation of the MULTICOMPO uses the Ceschino method with cubic Hermite polynomials as presented in Ref 17 default option LEAK keyword used to introduce leakage in the embedded MACROLIB This option should only be used for non regression tests b2 the imposed buckling corresponding to the leakage NMIX keyword used to define the maximum number of material mixtures This information is required only if MLIB is created nmixt the maximum number of mixtures a mixture is characterized by a distinct set of macroscopic cross sections the MACROLIB may contain The default value is nmixt 0 or the value recovered from MLIB if it appears on the RHS COMPO keyword used to set CPONAM and to define each global and local parameter TABLE keyword used to set CPONAM and to recover some global and local parameter from a MAP object named MAPFL CPONAM character 12 name of the LCM object containing the MULTICOMPO data structure where the interpolation is performed This name must be set in the RHS of the NCR data structure NAMDIR access the MULTICOMPO structure of CPONAM from the
31. seses K Ey OO k we 65 IGE 344 vi Structure desedata2 o o sa cee k HR ee HE EAM ea 65 Atticus INCHES e ee hb ee 67 Structure neridabay i ss e e a SS 67 Strchurs descinth cansada Ree ee ee a ee e E 69 NCR inputs for instantaneous cases 6 a e 74 NCR Inputs for TA gases oscila a Ae AR a We AAA A ee ooo 75 Structime SCR ail ee Ea TE 78 Bimuichune Sr dataj oa RARA 78 Structure deseints serora i AR 80 Structure descdepl ks 82 SCR inputs for instantaneous cases e e 86 SCR inputs lor TA Cases LL e ee A ee hee a ea a 87 DUCA i au Se S Rae ee Rd es a ee 88 Structure deseaiii o s i ire reggea adala eee ee eed bb ele 88 APM options summary peress asera Aone eS So ee ee ni 92 purucpore TRSCPOR Pas eh Lee eRe E eae oes 93 Structure desetl6cpo ic ee rd ee ee ee we 94 Structure DIPI oe ee a ae A ee E 100 Siruchire descphasel bic sesame uw a en Deni aa na AAA Be Mes 101 SUTILES TEME so ks eh ee e a ee ee ee A e 107 Biructine desethbi isa ci 0g oe Aw RS ee ee ea es Le 107 Structure DLBAK oe eas we dace dE DR Be a eS eG ee ees 112 Structure dlealedata o sacca k bee AA 112 Druce GRADE LR AMA REALI ELARIO 114 ee Bram Mate iui sc una e oo a aria i a ad ge ed 114 Sirene PLUS gt ariano ha he 4 ERASER ER ER oe Oe ere 118 Structure pla dat Lk e A a ee ee oe E 118 Sucre NAP iin ole KOR AD Ree a oe hee Be Oe BS 120 Structure desenapl e s esre te See eee A ew Ee be ee a Oe eS 121 Structure desena
32. the grid is represented by the thick gray lines on the axis on Fig 6 and 7 If we use the notations of Fig 6 and 7 the best estimate interpolated values f we can get are given by f fV fo f Va f Vo f V8 F Vo F Ve F Vo f V8 f Va f Vo F Ve f Vo for instataneous f fF V V f g Va f Va f W f Vo F Vo for TA Note that the reference point Vo in the example does not have to be the same for all parameters Database structures such as represented on Fig 8 can also been used In this case we even have two choices for the Af computation on axis A The last case is in fact a mix of cases 2 and 3 The gray rectangle and the gray line on Fig 9 and 10 reprensent where all the lattice computations have been performed With the notations used on those figures one can write that the best estimate interpolated values f we can get are given by f fV fVe fVac f V_ f Va f Vo HV Bo F Va J Vo for instataneous F VV fh Va F Vic Vac Va Vs H Va Fo Vos Vac F Va f Vo for TA Note once again that the reference point Vo in the example does not have to be the same for all parameters Database structures such as represented on Fig 11 can also been used IGE 344 85 The input files will actually reflect the previous eguations However they are different if the param eters are stored in a MAP object MAPFL or provi
33. 12 user defined keyword associated to a global or local parameter to be set value of a global or local parameter used to interpolate vall is the initial value of this parameter in case an average is required vall can be an integer real or string value value of the final global or local parameter By default a simple interpolation is performed so that val2 vall val2 is always a real value with val2 gt vall keyword used to indicate that the value of parameter vall or the second value for the Ao calculation is recovered from MAPFL i e the MAP object containing fuel regions description keyword only available together with the ADD option It is used to set all the other variable values when a A contribution is performed for one variable value of the reference parameter when it is directly given by the user Note that there is no default value keyword used to specify that the reference value will be the same as in the refence case i e for the Cret computation keyword only available together with the ADD option It is used to specify that all the other variable values which are required are given keyword used to set the number densities of some isotopes present in the MULTICOMPO object The data statement MICRO ALL is used by default keyword to indicate that all the isotopes present in the MULTICOMPO object will be used in the MICROLIB and MACROLIB objects Concentrations of these isotopes will be recovered from the
34. 39 40 59 AXIS 20 26 27 axn 46 47 AXNAME 9 10 AYNAME 9 10 B EXIT 39 40 B ZONE 9 11 b2 67 68 78 79 BARR 101 103 153 BFAC 101 106 bfac 101 106 BLIN 89 91 BORON 89 90 BOWR 108 111 BRANCH_INFO 155 BRANCHES information optional 98 BREFL 56 57 62 brnpar 57 BTYPE 12 BUNBLES 58 BUND 12 14 37 38 BUNDLE POW 12 13 BUNDLES 57 bunl 57 58 BURN 45 46 58 89 90 102 103 153 burn 65 66 BURN STEP 43 46 47 burn0 65 66 burn1 65 66 burncycle 45 46 burncyclel 46 IGE 344 burncycle2 46 BURNUP 56 57 59 61 BURNUP information optional 98 BURNUP BEG 59 bval 89 90 bvalue 12 13 ByPass 153 bypass 104 CALCUL DX 118 caltype 107 108 CBOR 102 103 153 CD 95 CELL 12 14 CELLAV 94 96 CELLDIR 141 celldir 61 CELLID 59 CELLID 56 57 61 CELSIUS 110 CFC 88 cflux 107 109 CHAIN 82 CHAN 12 14 43 44 CHAN VOID 53 55 CHANNEL 147 CHANNELS 57 58 CHECKER 53 54 CHECKER 1 2 53 54 CHECKER 1 4 53 54 Chi inV YLD Bet Lam 99 COMB 12 13 COMMENT 101 105 comment 101 105 COMPARE 46 47 COMPO 8 10 120 121 124 COMPO 64 67 68 COMPO 5 22 63 64 66 compo 102 COMPO 5 68 conc 65 66 69 70 80 82 CONDC 108 110 CONDF 108 110 CONDG 94 control 115 116 CONV 108 111 CONV TEST 118 119 COST EXTRAP 118 CP 95 CPO 63 65 CPO 93 96 CPO 5 CPONAM 67 68 188 CPONAM1
35. 90 reaction 82 83 READ 31 32 64 REC 31 32 recl 31 32 rec2 31 32 REF 51 52 69 70 80 81 ref 108 110 REF SHIFT 12 13 REFLECTOR 101 104 refnam 102 refnam i 101 REFT 88 89 REFUEL 43 REGION 94 96 RELAX 39 40 107 108 relax 107 108 relval 39 40 rep 23 24 RES 67 68 78 79 RESINI 9 RESINI 4 5 8 9 39 40 42 45 53 125 127 RESP 23 24 REWIND 155 RMIX 15 16 ROD 20 29 33 34 ROD GR 18 19 ROD GROUP 29 30 135 rod group 18 19 21 ROD ID 21 ROD NAME 20 RODMESH 108 111 RP 95 RT 95 rtime 43 46 47 SAHA 108 111 SAM 89 91 IGE 344 SAME 12 14 31 32 SAMEASREF 69 70 80 81 SAP 100 102 104 SAP 101 103 153 SAP 5 SAPHYB 78 79 SAPHYB 5 SAPHYB_INFO 152 153 SAPNAM 78 79 sapnam 102 sapnam i 101 SAPNAMI 78 SAPNAM2 78 sass 108 109 SCR 78 79 SCR 78 100 151 scr_data 78 80 scr_data 78 SECOND 43 46 47 sect 108 109 SEL 101 104 153 SEQUASCII 6 7 SEQ_BINARY 6 7 SET 46 47 53 54 69 71 80 81 83 121 SET PARAM 12 13 46 48 108 111 SFAC 101 106 sfac 101 106 SHUFF 43 44 side 10 SIM 9 11 SIM 4 45 SIMEX 52 SIMPLEX 115 SLPITX ASS 123 SLPITY ASS 123 SM149 101 102 smnam 101 102 SMOOTH 12 SOME 43 SPC 47 SPEC 46 47 spec 22 23 SPECTRAL 22 23 SPEED 20 29 30 33 34 108 109 speed 20 29 30 33 34 SPLINE 52 SPLIT NAP 9 sr 115 S
36. BURN keyword used to compute the derivatives of cross section information relative to the exit burnup of a single combustion zone The derivatives are computed using Eq 3 3 of Ref 16 written as 05 1 a ee S Sa pee ral dB NS B e99 Esl SR Bre oc DE B o OB Bi Bie Bee il Boe j k where Be B and By are the beginning of cycle burnup of bundle j k end of cycle burnup of bundle i j k and exit burnup of channel j This option is available only if TABLE is the selected option IGE 344 ivarty I BURNUP burn T BURNUP burn0 burnl MICRO HISO conc ALL ENDMIX 66 index of the combustion zone for differentiation of cross section information keyword used to perform a single interpolation and to set the burnup interpolation value burn real interpolation value of the burnup given in MW day per tonne of initial heavy elements keyword used to perform a time average MACROLIB evaluation between the burnup values burn0 and burn1 real initial value of the burnup given in MW day per tonne of initial heavy elements real final value of the burnup given in MW day per tonne of initial heavy elements keyword used to set the number densities of the extracted isotopes present in the COMPO linked list or xSM file By default the extracted isotopes are not added to the resulting MACROLIB character 12 name of an extracted isotope user defined real number density of the extracted isoto
37. Center Benchmark Problem Book ANL 7416 Supp 2 ID11 A2 Argonne National Laboratory 1977 S Loubi re R Sanchez M Coste A H bert Z Stankovski C Van Der Gucht and I Zmi jarevic APOLLO2 Twelve Years Later paper presented at the Int Conf on Mathematics and Computation Reactor Physics and Environmental Analysis in Nuclear Applications Madrid Spain September 27 30 1999 J Griffiths WIMS AECL Users Manual Report RC 1176 Atomic Energy of Canada Limited Chalk River Ontario 1994 J V Donnelly Post Processing WIMS AECL Results with Proc16 Report FCC RCP 001 Atomic Energy of Canada Limited Chalk River Ontario 1997 J V Donnelly Wrfsp Post Processing Data from WIMS AECL to RESP Report FCC RCP 006 Atomic Energy of Canada Limited Chalk River Ontario 1997 See the home page at http freesteam sourceforge net D Rozon A H bert and D McNabb The application of generalized perturbation theory and mathematical programming to equilibrium refueling studies of a CANDU reactor Nucl Sci Eng 78 211 1981 R Chambon E Varin and D Rozon CANDU fuel management optimization using alternative gradient methods Ann Nucl Energy 34 1002 2007 J A Ferland A linear programming problem with an additional quadratic constraint solved by parametric linear complementarity Publication number 497 D partement d informatique et de recherche op rationnelle Universit de
38. DETINI 4 5 22 137 DEV LZC 135 136 dev lzc 26 DEV ROD 134 135 dev rod 18 20 DEVICE 18 25 29 33 35 device 133 134 DEVICE 5 18 25 29 33 35 DEVINI 4 5 18 25 30 34 133 DFLUX 114 DIFF 112 113 DIMENSIONS 57 DIRGEO 123 DIST AX 46 47 DIST BURN 46 48 DIST POWR 46 48 DISTR 37 38 DIVERS 157 DIVERS 155 157 DLEAK 112 DLEAK 112 dleak_data 112 DMACRO 112 DMIX 20 21 DMOD 102 103 153 dmod 89 90 DNAME 88 89 189 DONCPO 93 94 DONJON 6 DRA 101 103 153 154 DSET 4 18 21 25 29 133 dt 51 52 ecost 118 119 EDI 68 EDIT 9 12 15 19 22 26 29 31 33 35 37 39 40 43 45 46 49 51 53 57 64 67 78 79 88 89 94 100 101 107 108 112 115 118 121 123 emin 104 EMPTY MIX 26 27 END 29 30 END 3 6 7 ENDADD 101 ENDBARR 101 ENDCHAIN 82 83 ENDCYCLE 46 47 ENDGLOBAL 101 ENDISOTOPES 101 ENDMIX 65 66 69 71 80 82 ENDN 23 24 ENDPKEY 101 ENDREF 69 70 80 81 ENDROD 20 21 energy 82 83 ENRICH 12 14 epsburn 46 48 epsilon4 116 epspowr 46 48 EQUI 78 79 ermaxc 108 111 ermaxt 108 111 EVO 42 58 EXTR 33 34 F 105 FACTOR 112 113 FADE 18 19 FILE_CONT_1 101 104 153 FILE_CONT_2 101 104 153 FILE_CONT_3 101 104 153 154 FILE_CONT_4 101 104 153 FILE_NAME 101 105 FIXP 88 89 FLMAP 8 9 13 FLMAPI 31 32 FLMAP2 8 9 12 14 31 32 flow
39. E 5 REAL Eps Keff Bexit x compo files k SEQ_ASCII SFUEL1 FILE CpoFuell SEQ_ASCII SFUEL2 FILE CpoFuel2 SEQ_ASCII SREFL1 FILE CpoModel SEQ_ASCII SREFL2 FILE CpoMode2 SEQ_ASCII SZCRin FILE CpoZCRin SEQ_ASCII SZCRot FILE CpoZCRot SEQ_ASCII SSORin FILE CpoSORin SEQ_ASCII SSORot FILE CpoSORot x compo directories STRING NamFuel1 FUEL1 Lit STRING NamFuel2 FUEL2 1 3 STRING NamRef11 MODE1 gws STRING NamRef12 MODE2 prz 159 IGE 344 160 STRING NamZCRin ZCRIN da STRING NamZCRot ZCROT de STRING NamSORin SORIN os STRING NamSORot SOROT 13 A ii he ee sete FULL CORE CALCULATION FAIDA geometry construction GEOM Pgeom COLE reactor material index si GEOM MATEX USPLIT GEOM EDIT O NGRP 2 MAXR lt lt MaxReg gt gt NREFL lt lt nbRef1 gt gt RMIX lt lt mRef11 gt gt lt lt mRef12 gt gt NFUEL lt lt nbFuel gt gt FMIX lt lt mFuel1 gt gt lt lt mFuel2 gt gt ds numerical discretization IF Method MCFD THEN TRACK TRIVAT GEOM EDIT 1 MAXR lt lt MaxReg gt gt MCFD lt lt degree gt gt ELSEIF Method PRIM THE TRACK TRIVAT GEOM EDIT 1 MAXR lt lt MaxReg gt gt PRIM lt lt degree gt gt ELSEIF Method DUAL THE TRACK TRIVAT GEOM EDIT 1 MAXR lt lt MaxReg gt gt DUAL lt
40. HELIOS dra keyword used to set der character T F type of data in non reference branches of output PMAXS file If der T data are partial derivatives otherwise it is raw cross sections Default der T keyword used to set vers the version of PARCS which will be used If vers gt 2 705 generate PMAXS for PARCS 2 71 or later versions otherwise it is for PARCS 2 7 or earlier versions Default vers 3 0 keyword used to set a comment line character 40 Comment line for the user Default comment PWR CASE UOX MOX CORE FUEL keyword use to set logical flags it indicates write or not write correponding data into PMAXS file If the flag is F default values given in Ref 40 will be used in PARCS For reflector case all flags will be forced to F except for ladf and linv character T F assembly discontinuity factor Default ladf F character T F microscopic cross section of Xe and Sm Default lxes F character T F direct energy deposition fraction Default Ided F character T F J1 factor for minimal critical power ratio Default j1f F character T F fission spectrum Default Ichi F character T F delay neutron fission spectrum Default Ichid F character T F inverse neutron velocity Must be T for transient Default linv F character T F Detector response Default Idet F character T F yield values of I Xe and P
41. IGE 344 141 Main records and sub directories in history continued from last page Type Condition Units Comment FUELDIR Dir list of sub directories FUEL that contain the prop sd erties associated with the fuel type Fi j CELLDIR Dir list of sub directories CELL that contain the prop J erties associated with the cell Ci The signature for this data structure is SIGNA L_HISTORY ouu The array Sh contains the following information e S N contains the number of global parameters e Si N contains the number of local parameters Sy N contains the number of bundles per channel St N contains the number of channels in the core e S N contains the number of bundle shift o Se T contains the type of depletion solution used e S T contains the type of burnup considered Sh Nr contains the number of isotopes o Se G contains the number of transport groups e Si N contains the number of regions e S Nr contains the number of fuel types The fuel directory name FUEL associated with fuel type F j is composed using the following FOR TRAN instruction WRITE FUEL A4 18 8 FUEL Fi This directory will contain the initial isotopic content of this fuel type The cell directory name CELL associated with C j is composed using the following FORTRAN instruction WRITE CELL A4 18 8 CELL Ciy This directory will contain the value of the local parameters
42. MACRO TABLE CPO default B MIX 1 SET A lt lt va gt gt SET C lt lt vc gt gt ENDMIX MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET A lt lt va gt gt SET C lt lt vc0 gt gt ADD C lt lt vc0 gt gt lt lt vc gt gt REF PA lt lt va0 gt gt B SAMEASREF ENDREF ENDMIX MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 ENDMIX MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET C lt lt vc0 gt gt ADD C lt lt vc0 gt gt MAP REF A lt lt va0 gt gt B SAMEASREF ENDREF SET A lt lt va0 gt gt ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B SAMEASREF ENDREF ENDMIX continued on next page IGE 344 MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET A lt lt va0 gt gt SET C lt lt vc0 gt gt ADD A lt lt va0 gt gt lt lt va gt gt REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF ADD C lt lt vc0 gt gt lt lt vc gt gt REF PA lt lt va0 gt gt B lt lt vb0 gt gt ENDREF ENDMIX PLANE AXE Fig 10 MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET A lt lt va0 gt gt SET C lt lt vc gt gt ADD A lt lt va0 gt gt lt lt va gt gt REF C lt lt vc0 gt gt B lt lt vb0
43. PARKEY vall val2 MAP REF valref SAMEASREF ENDREF MICRO ALL ONLY 81 only if a MAPFL object is set By default the type of calculation TIMAV BURN or INST BURN is recovered from the MAPFL object index of the combustion zone for differentiation of cross section information keyword used to indicate a simple interpolation at vall or an averaging between vall and val2 The result 0 ef is also used as the reference value when the ADD is used Note see at the ending note of this section for a detailed description and examples keyword used to indicate a delta sigma calculation between val2 and vall i e Aoret Oval2 Oval is computed Note see at the ending note of this section for a detailed description and examples keyword used to indicate a delta sigma calculation between val2 and vall is added to the reference value i e Ao Oyaj2 Oval is used as contribution Gref Ao or Adyer Ao is returned Note see at the ending note of this section for a detailed description and examples keyword indicating that interpolation of the SAPHYB for parameter PARKEY uses linear Lagrange polynomials It is possible to set different interpolation modes to different parameters By default the interpolation mode is set in Sect 4 3 1 keyword indicating that interpolation of the SAPHYB for parameter PARKEY uses the Ceschino method with cubic Hermite polynomials as presented in Ref 17 By default the interp
44. RECNAME WRITE RECNAME 5HFINF_ 13 3 IFX Thus for example for ifx 2 RECNAME would be FINF_002 SET keyword to specify assembly calculations at which pho have been calculated Re peated as many times as there are parameters in the MULTICOMPO pname name of the parameter in the MULTICOMPO pvalue value of the parameter in the MULTICOMPO IGE 344 122 Note that in the case of heterogeneously homogenized assembly the pin wise projected diffusion flux is stored in mixture 1 7 1 2 Pin power reconstruction EDIT iprint PPR NZASS nzass METH GPPR ifx POWER pow ti where EDIT iprint PPR NZASS nzass METH GPPR ifx POWER pow Table 80 Structure descnap2 keyword used to modify the print level iprint index used to control the printing of this module The iprint parameter is important for adjusting the amount of data that is printed by this calculation step e iprint 0 results in no output e iprint 1 keyword to perform pin power reconstruction and stored results in the MAP data structure MAP keyword to specify nzass number of mesh in Z direction along assemblies keyword to select the type of methodology used for pin power reconstruction keyword to select the generalized pin power reconstruction from an heterogeneous assembly definition several mixtures per assembly Note that if there is only one mixture homogeneous assembly this method is actually the regular
45. Rha REA ws 166 9 3 1 Input file for geometry Li o 166 9 3 2 Inpub file for devices cocer a a eae E a 169 9 3 3 Input file for fuel map s ssu i LE a A A ee 178 9 3 4 Input file for exit burnups o 181 a s easa ra aa es A ESB RE MR A Se ee S eo BAS SSO 184 ie a we oe a A E ese ees PRL ES es a ee ES 187 IGE 344 v oOo N D Ubs WIM oP SBP i je fea a 49 49 12 13 0 WAN D 00 N N NINNINNNNNDN_r_rrrrrrrro CSCO ON SD TU 0 N HO O AN OOBWNF Oo 0 XS IA WIN OS 6021 00 0 EO0O0ONR O List of Tables Structite DONJON k e I n BOS 6 MORIRE 8 Structure deseresinill 0 4 06 68 bee ee eee n a 9 Structure descresini2 ee eee e 12 Structure USPLITS S eo Re A we ae BR ee a 15 Structure deselink i s gt su ba cde ee hn S a 15 Structure MACINES oa Gad So e ae e a A a A 17 Sirueturs DEVINIS se picio aa OL A A de re Ak a EEN E a 18 Structure desedey o e 244s 4 860 PREG ee SENS a ee ER 18 Structure der eea a Li A a a ee ee wo Owais 20 Structure rod sroup ii 44 504 ee nnd bane a aa ee we 21 Stricture DETIENE 2 44 64 amp fe AAA A a a oe Aka we ede 22 Structure descinid t lt icare k e a G 22 Structure desedet i i O eee RR 23 Structure MACY Lili e pie e e n ne e Pad n a 25 Structure deselac gt i ia 26 Structure dA6YIZ L ili bee E e td i ai A EN 26 Structure E BEOUP i e RI a E 28 Sirene DSET o 4 5 gos eee oR K IO a Sa 29 Structure desedset k k k k Ee eae SEER aa 29 Strictu
46. Structure descafm MAP MCR mmix INFOR NAMDB DNAME ntyp NAMTYP i i 1 ntyp REFT imix i NAMTYP i 1 ntyp EDIT iprint FIXP INIT pow y PWF NPWF TFUEL tfuel TCOOL tcool TMOD tmod continued on next page IGE 344 89 Structure descafm continued from last page BORON nB RDCL dcool RDMD dmod PUR purity BURN bval XENON nXe XEREF NEP nNp NREF y SAM nSm IMET imet BLIN where MAP MCR mmix INFOR NAMDB DNAME ntyp NAM TYP i REFT imix i EDIT iprint FIXP INIT pow PWF NPWF keyword to specify that a MACROLIB for fuel properties will be computed keyword to specify that a MACROLIB containing only one non zero mixture will be created maximum number of mixtures in the MACROLIB keyword to specify the data base name character 72 title of the database as it has been created keyword to specify the number of fuel types and their names as stored in the data base number of fuel types For MCR option ntyp must be 1 character 12 name of the directory where each fuel type information has been stored keyword to specify a number associated with a fuel type name fuel type index as specified for the fuel map or a non zero mixture number for the single property sc macrolib keyword used to set iprint index used to control the printing in module AFM 0 for no print default value 1 for minimum printing l
47. VAR WEIGHT 115 116 varmax 115 116 varmin 115 116 VChan 153 154 vchan 104 VClad 153 154 vclad 104 VCnRd 153 154 venrd 104 VCool 153 154 vcool 104 vcool vwatr vmodr venrd vfuel vclad vchan 101 vecmax 115 116 vecmin 115 116 velocity 108 109 vers 101 105 VERSION 101 105 VFCM 154 vicm 105 VFuel 153 154 vfuel 104 VModr 153 154 vmodr 104 vnorm 51 52 VOID PATTERN 53 54 VWatR 153 154 vwatr 104 WARNING ONLY 118 WEIGHT 12 14 weight 115 116 width 10 196 WIMS16 93 94 wt 14 wteff 108 110 X 20 26 27 xbe 104 XE135 101 102 xenam 101 102 XENON 89 90 XENON 49 50 XEREF 89 91 XFAC 35 xfac 35 XNAME 9 10 53 55 XNAMEA 9 10 XS CONTROL 99 XS CONTROL information mandatory 98 XS data mandatory 99 XS SET identification mandatory 99 XS_CONT 101 104 154 XSM 66 XSM_FILE 6 7 Y 20 26 27 ybe 104 yield 82 83 YNAME 9 11 53 55 YNAMEA 9 10 2 390 30 97
48. array containing the pin power for each pin on each mesh plane of the assembly ASS POWER uuu total assembly power SIGF PHl uu array containing the integral of Xp times the flux for each pin on each mesh plane of the assembly continued on next page IGE 344 131 Assembly type sub directory continued from last page Condition UnitsComment PU sueina NN Ng array containing the flux of each group for pin each pin on each mesh plane of the assembly 8 2 Contents of matex data structure A matex data structure is used to store several information related to the reactor extended material index and geometry This object has a signature L_MATEX it is created using the USPLIT module The information contained in this data structure can be used and updated in other DONJON modules 8 2 1 The state vector content The dimensioning parameters S which are stored in the state vector for this data structure represent e The number of energy groups Ng Si e The maximum number of material mixtures Nm S2 Nm equals to the total number of material regions plus the number of device mixtures e The number of reflector types N 3 e The number of fuel types Nr Sa e The total number of mixtures indices Nio Ss Nioz equals to the total number of mesh splitted volumes plus the number of device mixtures e The type of reactor geometry I S only I 7 for 3D Cartesian geometry or I 9 for 3D Hexagonal geomet
49. contoller EMPTY MIX uuu The empty part mixture number and the ref erence mixture number FULL MIX uuu The full part mixture number and the refer ence mixture number RATE uuuuuuu The water filling rate TIME ou The water filling time 8 3 6 The LZC GROUP sub directories Inside each LZC GROUP sub directory the following records will be found Table 92 Records in LZC GROUP sub directories Type Condition Units Comment GROUP ID uuu The identification number of the lzc group NUM LZC uuu The total number Na of lzc devices in the group LZC ID unio An array of identification numbers of liquid controllers which belong to the same group IGE 344 137 8 4 Contents of a detect data structure The detect data structure is used to store detector positions and responses This object has a signature L_DETECT it is created using the DETINI module The information contained in this data structure can be used and updated in other DONJON modules which are related to the detectors namely DETECT and DETINI modules 8 4 1 The state vector content The dimensioning parameters S which are stored in the state vector for this data structure represent e The number of energy groups Ng Si e The total number of detectors X Z1 S2 e Flag for hexagonal detector definition S3 1 for hexagonal detector definition 0 otherwise The dimensioning parameters for a specific detector type which are stored in the vecto
50. cross sections respectively These indices will be used to compute the incremental cross sections in the NEWMAC module IGE 344 28 3 6 3 Description of lzc group input structure The partition of lzc devices into groups is similar to that of rod devices Table 18 Structure lzc group GROUP ID grp LZC ID id ALL where GROUP ID keyword used to set igrp number igrp integer identification number of a group to be created Each controllers group must be assigned a unique identification number given in ascending order ranging from 1 to ngrp LZC ID keyword used to set the controllers id numbers id integer identification numbers of the liquid controllers which belong to the same group igrp A particular controller or several devices may belong to different groups but it could not be repeated inside the same group The total number of liquid controllers in any group must be between 1 and nlzc ALL keyword used to specify that all liquid controllers will belong to the same group igrp IGE 344 29 3 7 The DSET module The DSET module is used to set or to update some of the devices parameters The new parameters can be applied for the rod type devices and or for the liquid zone controllers such as the new inser tion level for the rods or water filling level for the lzc type devices etc It is possible to apply the new parameters to the individual user selected devices as well as to the user selected groups of d
51. desccre1 READ COMPO CPO descdatal EDIT iprint Table 49 Structure desccre2 READ TABLE CPO descdata2 3 where EDIT iprint NMIX nmix READ COMPO TABLE CPO descdatal descdata2 keyword used to set iprint integer index used to control the printing of information on screen 0 for no print 1 for minimum printing larger values will produce increasing amounts of output keyword used to define the number of material mixtures nmix This data must be given only if MACRO is created and the FMAP object is not specified integer maximum number of reactor material mixtures as defined in the reactor ge ometry keyword used to read the MACROLIB specification from the input data file keyword used to indicate a simple MACROLIB creation i e according to the first calling specification when FMAP object is not specified keyword used to indicate a fuel map MACROLIB creation i e according to the second calling specification with FMAP object specified character 12 name of the selected COMPO object This name must appear in the calling specification to the CRE module structure containing the interpolation specification if COMPO is the selected option structure containing the interpolation specification if TABLE is the selected option IGE 344 65 Table 50 Structure descdata1 MIX mix NAMDIR DERIV UPS I BURNUP burn T BURNUP burn0 burn1
52. eee 53 31S The BST mde II dee es Khaw SSA A a eS eR a 56 IGE 344 111 S 18 1 Example aaa ea ee ane ee Rnb ee heb baa ee 59 4 CROSS SECTION INTERPOLATION MODULES 63 4 1 The GRE modale ooo ek E we e ee e RE ES 63 4 1 1 Input data for the CRES module i Re ee 64 4 2 The NGR nodale i e o a a de ed 67 4 2 1 Interpolation data input for module NCR Ls 67 4 2 2 Defining local and global parameters o o o 69 4 2 3 Interpolation in the parameter grid 71 4 3 The sori Module oo a hee BS ba ee B ha ae eee 78 4 3 1 Interpolation data input for module SCR 78 4 3 2 Defining global parameters o o ee eevee 80 4 34 Depletion data Structure k RT 82 4 3 4 Interpolation in the parameter grid saaua coeno anaso 83 4 4 The APM modale cs o a semea a Anda eee pa edb bois 88 4 4 1 Input data to the AFM module o e 88 4 5 The TIGEPO modde gt coso a al a oes 93 4 5 1 Input data for the TI6CPO module 94 4 6 The D2P Module e na bee a ee ee lado aa ee ed 97 4 6 1 The PMAXS format oscars a RRR EE k 97 4 6 2 General format of the module o o 100 5 THERMAL HYDRAULICS MODULES o e 107 5 1 The THM adult k a S we Oe a E 107 5 1 1 Input data to the THM module aoaaa a a 107 6 OPTIMIZATION MODULES socorrer LEED ERS 112 6 1 The DLEAK module s sa soa k 84
53. emin 6 2506E 01 1E 04 keyword used to set the FILE_CONT_3 block volume of regions See Ref 40 volume of coolant Default vcool 2 4921E 02 volume of water By default vwatr 0 0000E 00 volume of moderator Default vmodr 2 4921E 02 volume of control rods By default venrd 2 3020E 01 volume of fuel By default vfuel 1 4407E 02 volume of cladding By default vclad 4 5099E 01 volume of channel By default vchan 4 5099E 01 keyword used to set the FILE_CONT_4 block See Ref 40 rod lattice pitch cm Default pitch 1 44270E 00 position of first column rods cm By default xbe 7 21350E 01 position of first row rods cm By default ybe 7 21350E 01 keyword used to set the XS_CONT block See Ref 40 number of sides in assembly Default nside 4 IGE 344 ncorner v cm GENPMAXS JOB_TIT jobtit FILE_NAME fname DERIVATIVE der VERSION vers COMMENT comment JOB_OPT ladf Ixes lded Ij1f Ichi Ichid linv Idet lyld ledf Igff lbet lamb Idec 105 number of corners in assembly By default ncorner 4 value of vfem By default vfem 5 32151E 01 keyword used to indicate that the input data for the GEN file will be set by the user keyword used to set jobtit character 16 name of the PMAXS file created by the D2P module Default jobtit D2P PMAX keyword used to set fname character 12 name of the HELIOS like file HEL created by the D2P module De fault fname
54. gt gt ENDREF ENDMIX 75 only the burnup in MAP MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET A lt lt va0 gt gt SET C lt lt vc0 gt gt ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF ADD C lt lt vc0 gt gt MAP REF A lt lt va0 gt gt B lt lt vb0 gt gt ENDREF ENDMIX MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET A lt lt va0 gt gt ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF ENDMIX Table 56 NCR inputs for TA cases IGE 344 76 The following pictures correspond to the previous different examples Figure 2 Complete grid one point case Figure 5 Partial grid complete planes TA case AZ Figure 3 Complete grid TA case Figure 6 Partial grid complete axis one point case Figure 4 Partial grid complete planes one point case Figure 7 Partial grid complete axis TA case IGE 344 77 Figure 8 Partial grid complete axis with another Figure 10 Partial grid one complete plane and one configuration one point case complete axis TA case VV Figure 11 Partial grid one complete plane and one Figure 9 Partial grid one complete plane and one complete axis with another configuration one point complete axis one point case case IGE 344 78 4 3 The SCR module This compo
55. if present lzc type device The computing algorithm is based on the determination of the volumic fraction occupied by each device the incremental cross sections are then adjusted accordingly Note that the NEWMAC module must be executed each time the devices positions are modified from the previously computed ones The NEWMAC module specification is Table 25 Structure NEWMAC MACRO3 MATEX NEWMAC MATEX MACRO2 DEVICE EDIT iprint XFAC xfac where MACRO3 character 12 name of the MACROLIB to be created by the module It will contain the updated properties of each material region with respect to the current position of each device MATEX character 12 name of the MATEX object containing the complete reactor material index including devices MATEX must be specified in the modification mode it will store the updated h factors computed per each fuel region with respect to the devices positions MACRO2 character 12 name of the read only extended MACROLIB previously created by the MACINI module DEVICE character 12 name of the read only DEVICE object containing the devices information and parameters EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing larger values produce increasing amounts of output The default value is iprint 1 XFAC keyword used to specify the number of cells on which incremental cross sections were computed in
56. in the modification mode it will store the recovered h factors per each fuel region MACRO character 12 name of a MACROLIB created using either MAC CRE NCR or AFM module for the evolution independent material properties see structure desccrel or refer to the user guidel MACFL character 12 name of a fuel map MACROLIB created using either CRE NCR or AFM module for the interpolated fuel properties see structure desccre2 or refer to the user guidel EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing larger values produce increasing amounts of output The default value is iprint 1 IGE 344 18 3 4 The DEVINI module The DEVINI module is used for the modeling of reactivity mecanisms based on the devices specifi cations which are read from the input data file The module will create a new DEVICE object that will store the devices specifications and parameters see Section 8 3 Note that only the rod type i e solid devices are considered using the DEVINI module the liquid zone controllers can be added subsequently using the LZC module see Section 3 6 A rod type device is a reactivity controller rod or plate such as a zone control rod ZCR a shutoff rod SOR etc Several devices parameters can be modified using the DSET module see Section 3 7 A device specification includes several controller rod parameters
57. iteration VAR VALUE2 jj Naar The values of the decision variables of the second last valid iteration BEST VARuuuu The values of the decision variables corre sponding to the best valid solution ever found BEST FCT uuu The value of the objective function corre sponding to the best valid solution ever found IGE 344 151 8 9 Contents of d2p_info data structure A d2p_info data structure is used to store cross sections and Saphyb information such as keff kinf assembly discontinuity factors This object has a signature L INFO it is created using the D2P module 8 9 1 The state vector content The dimensioning parameters S which are stored in the state vector for this data structure represent The number of energy groups Ngr S The number of pkey for the generation of the PMAXS file Npk S2 The number of cross sections recovered using the SCR module Nx S3 The number of points with respect to burnup Npu Sa The type of branching calculation requested for the generation of PMAXS tree Isria S5 default meshing according to the saphyb content user defined with global otpion user defined with add otpion user defined with add new otpion Isrid BWIN The number of control rod composition Neomp S6 The number of delay neutron groups Nagel S7 The number of columns in assembly Ncols Sg The number of rows in assembly Nrows Sy The computed part of asembly assembly part S10 0 whole Es
58. lt degree gt gt lt lt quadr gt gt ENDIF k macrolib for reflector kh LREFL1 SREFL1 LREFL2 SREFL2 MACRO1 CRE LREFL1 LREFL2 EDIT 1 NMIX lt lt nbMix gt gt READ COMPO LREFL1 MIX lt lt mRef11 gt gt lt lt NamRef11 gt gt ENDMIX COMPO LREFL2 MIX lt lt mRef12 gt gt lt lt NamRef12 gt gt ENDMIX k device specification S DEVICE MATEX Pdevc MATEX lt lt mZCRin gt gt lt lt mZCRot gt gt lt lt mSORin gt gt lt lt mSORot gt gt es full insertion of ZCR devices fos DEVICE DSET DEVICE EDIT 1 ROD GROUP 1 LEVEL 1 0 END ROD GROUP 2 LEVEL 1 0 END ROD GROUP 3 LEVEL 1 0 END En update macrolib for devices IGE 344 yes LZCRin LZCRot LSORin LSORot MACRO1 x SZCRin SZCRot SSORin SSORot CRE MACRO1 LZCRin LZCRot LSORin LSORot READ COMPO LZCRin COMPO LZCRot COMPO LSORin COMPO LSORot MIX lt lt mZCRin gt gt lt lt NamZCRin gt gt MIX lt lt mZCRot gt gt lt lt NamZCRot gt gt MIX lt lt mSORin gt gt lt lt NamSORin gt gt MIX lt lt mSORot gt gt lt lt NamSORot gt gt fuel map specification x FMAP MATEX x average exit burnups x Pfmap MATEX FMAP Pburn FMAP k initialization kh LFUEL1 LFUEL2 SFUEL1 SFUEL2 EVALUATE Eps epsil 1 WHILE Eps epsil gt iter 20 lt DO EVALUATE iter iter 1 x
59. modify the print level iprint iprint index used to control the printing in this module It must be set to 0 if no printing on the output file is required while values lt 10 will print general information about each record requested on Tape16 as well as other generic information pertinent to the T16CPO module Finally for values of iprint gt 10 additional information required for debugging will be printed The default value is iprint 1 NMIX optional keyword used to define the number of mixtures created on the CPO data structure nmixt the maximum number of mixtures created The default value is nmixt 1 CONDG optional keyword used to define the group structure for condensation In the case where the CPO is to be updated the information following CONDG must yield an energy group structure compatible with that already available on this data structure If it is absent the code will first try to use the CPO group structure if available Then it will try to use the editing group structure corresponding to NGREAC on the following Tape16 record REACTION Gy FLUX yu uuu NEL Finally if everything else fails it will select the main transport group structure corre sponding to NGMTR on the following Tape16 record IGE 344 ngcond LIST MIX MIXNAM CELLAV REGION noreg RC nburn frstrec NAMPER valref 95 WIMS uuuuu CONSTANT jj NEL the number of condensed groups required the last group number associated
60. name of isotopes in the saphyb ob ject for depletion chain recovery of fission yield and number densities 1 Xe135 2 1135 3 Sm 149 suname Values for each state variable specified by PKEY Ny is the number of value taken for the k PKEY continued on next page IGE 344 Records and sub directories in SAPHYB_INFO Condition Units PKIDX uuu I Npx ADFuvuuvuuuy 0x3 153 continued from last page Comment Correspondance of indices between SAP MCO and GenPMAXS state vari ables Type of ADF recovered GET or SEL Each component of the list PKEY_INFO is a directory containing for all possible state variables 1 BARR 2 DMOD 3 CBOR 4 TCOM 5 TMOD 6 BURN Inside each groupwise directory the following records will be found Table 107 Records in PKEY_INFO Type Condition Units LFLAG uuuuuu I 1 NAME Sou C 12 Iflag Comment Indicates if the corresponding state variable can be found in the SAP or MCO object Name of the state parameter in the SAP or MCO object For BARR and BURN state vari ables this record exists even if LFLAG is false Table 108 Records and sub directories in HELIOS_HEAD Type Condition Units FILE_CONT_1 FILE_CONT_2 FILE_CONT_3 FILE_CONT_4 Comment Set of data for FILE_CONT_1 block in DRA file Heavy metal Density HM_Dens Bypass Den sity ByPass Set of data for FILE_CONT_2 block in DRA file Lower Energy Limits of Neutron
61. names only a fixed normalisation can be performed INFO keyword to specify the information associated with the detector type ndetect number of detectors of the specified type nrep number of detector response components for the specified type It must be greater or equal to 2 corresponding to a response in fraction and the reference flux value SPECTRAL keyword to specify the energy spectral of a detector type spec array containing the energy spectral of a detector type DEFAULT keyword to specify the energy spectral will be initialized as 1 0 for the highest energy group and 0 0 for other groups INVCONST keyword to specify the inverse time constants of the detector type model This option is only valid for platinum VAMTYP 1 5 PLATN detector type tinv array containing the inverse time constants of the detector model FRACTION keyword to specify the fractions corresponding to each delayed or prompt reponse of the detector type model This option is only valid for platinum NAMTYP 1 5 PLATN detector type frac array containing the detector type model fractions descdet structure describing the format used to read detector information 3 5 2 Description of the detector data Note that the information input order must be respected Table 14 Structure descdet NAME NAMDET NHEX nhex HEX ihex i i 1 nhex POSITION pos i i 1 6 RESP rep i i 1 nrep ENDN IGE 344 where NAME
62. perturbed cross sections second of two integer rod mixture indices Index mix2 corresponds to the reference cross sections Indices mix and mix2 will be used to compute the incremental cross sections in the NEWMAC module keyword used to end the rod description 3 4 3 Description of rod group input structure The partition of devices into groups is very useful when the same action is to be applied to several rods e g setting of new parameters using the DSET module or rods moving using the MOVDEV module Table 11 Structure rod group GROUP ID igrp ROD ID id ALL where GROUP ID igrp ROD ID id ALL keyword used to set igrp number integer identification number of a group to be created Each rods group must be assigned a unique identification number given in ascending order ranging from 1 to ngrp keyword used to set the rod id numbers integer identification numbers of rods which belong to the same group igrp A partic ular rod or several rods may belong to different groups but it could not be repeated inside the same group The total number of rods in any group must be between 1 and nrod keyword used to specify that all rods will belong to the same group igrp IGE 344 22 3 5 The DETINI module The DETINI module is used to read and store detector information A detector is represented by a 2 D or 3 D Cartesian Hexagonal geometry The DETINI module specification is Table 1
63. r acteurs CANDU refroidis leau l g re Ph D Thesis Ecole Polytechnique de Montr al 2006 A H bert Revisiting the Ceschino Interpolation Method in MATLAB A Ubiquitous Tool for the Practical Engineer Clara M Ionescu Ed InTech Open Access Publisher ISBN 978 953 307 907 3 Croatia 2011 A H bert G Marleau and R Roy A description of the DRAGON Data Structures Report IGE 295 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2007 E Varin and G Marleau CANDU reactor core simulations using fully coupled DRAGON and DONJON calculations Ann Nucl Energy 33 682 2006 IGE 344 185 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 J Tajmouati Optimisation de la gestion du combustible enrichi d un r acteur CANDU avec prise en compte des param tres locaux Ph D thesis Ecole Polytechnique de Montr al 1993 M T Sissaoui G Marleau and D Rozon CANDU Reactor Simulations Using the Feedback Model with Actinide Burnup History Nucl Technology 125 197 1999 R Roy The CLE 2000 Tool box Report IGE 163 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 1999 A H bert and R Roy A Programmer s Guide for the GAN Generalized Driver FORTRAN 77 version Report IGE 158 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 1994 Argonne Code
64. set of global and local parameters For example the fuel bundles have a different evolution of the fuel properties according to the given burnup distribution which is a global parameter Consequently the homogenized cell properties will differ from one fuel region to another i e they are not uniform over the fuel lattice Thus the realistic modeling of a reactor core requires the fuel properties to be interpolated with respect to global and local parameters which must be specified in the fuel map Note that the above considerations correspond to the typical core modeling of CANDU or PWR reactors The RESINI module will create a new FMAP object that will store the information related to the fuel lattice specification and properties see Section 8 1 In PWR cases each channel correspond to an assembly Using heterogeneous mixtures in one assembly increases the complexity of the geometry However two levels geometries embedded geometry are not possible in the DONJON code The general idea is then to define one channel per mixture for all assemblies All these channels have then to be regrouped by assembly to impose the same burnup This process could be done manually but if the heterogeneity of the cross section is large ex one mixture per pin within a complete core the geometry definition may be too complex This task can be performed automatically by the module NAP The RESINI module specifications are Table 2 Structure RESINI
65. stored on the main CPO directory will corre spond to a catenation of MIXNAM and a perturbation name and an index i describing the perturbation order It is created using the following FORTRAN instructions for the reference mixture WRITE MIXDIR A6 A6 MIXNAM RC oui while for the i perturbed state associated with NAMPER J we will use WRITE MIXDIR gt A6 A2 A2 12 MIXNAM NAMPER J uu i Finally for the i perturbed state associated with the MTMD perturbation we will use WRITE MIXDIR A6 A4 12 gt MIXNAM MTMD i Typically if the desct16cpo structure takes the form EDIT 0 NMIX 2 MIX Candu RC 15 1 FT 900 0 2 1100 0 16 1300 0 46 Maple RC 70 RP 1 0 1 0 5 71 Then the first 15 cases stored on the Tape16 file will correspond to a reference CANDU fuel with burnup The reference fuel temperature is 900 0 K The next 15 cases are for a fuel temperature of 1100 0 K Finally cases 46 to 60 are for a fuel temperature of 1300 0 K The Maple mixture will have no burnup The reference Maple cross sections correspond to case 70 while case 71 contains the effect on the Maple fuel mixture cross sections of a 50 reduction in reflector purity As a result we will end up with a CPO data structure which contains 5 mixtures called respectively Candu RCuuuu Candu FTuuu1t Candu FTuuu2 MapleuRCouuu Maple RPuuut The beginning of a new case on Tape16 will be identified by the presence of the record CELLAV uuu MODERATO
66. structure L_TRACK signature containing the tracking character 12 name of the read only FLUXUNK data structure L_FLUX signature con taining a transport solution character 12 name of the MAP data structure L_MAP signature containing fuel re gions description global and local parameter information burnup fuel coolant tem peratures coolant density etc Keyword PPR is expected in descnap2 character 12 name of the read only MATEX data structure L_MATEX signature The object corresponds to the heterogeneously splited geometry Keyword PPR is expected in descnap2 character 12 name of the read only MACROLIB data structure L_MACROLIB signature containing a cross section for the fuel The MACROLIB data structure must have been created with a MULTICOMPO data structure with pin level properties transport flux H factor infinite domain diffusion flux Keyword PPR is expected in descnap2 character 12 name of the created GEOMETRY data structure L_GEOM signature con taining the detailed core geometry definition at heterogeneous assembly level character 12 name of the read only GEOMETRY data structure L_GEOM signature containing the core geometry definition with homogeneous assembly only 1 mesh per assembly mandatory structure containing the input data to this module to compute additional properties for subregions see Section 7 1 1 structure containing the input data to this module to perform pin power reconstru
67. sub directory named NAMDIR This value must be set equal to default if not previously defined by a STEP UP key word in module COMPO namburn name of the parameter for burnup or irradiation in the sub directory named NAMDIR This value is defined if option TABLE is set and if burnup or irradiation is to be con sidered as parameter descintf input data structure containing interpolation information relative to the MULTICOMPO data structure named CPONAM see Section 4 2 2 IGE 344 69 4 2 2 Defining local and global parameters If a MAP object is defined on the RHS of structure ner_data and if the TABLE keyword is set some information required to set the interpolation points is found in this object In this case the NCR operator search the MULTICOMPO object for global or local parameters having an arbitrary name specified in the MAP object or set directly in this module Note that any parameter s value set directly in this module prevails on a value stored in the MAP object Each instance of descintf is a data structure specified as Table 54 Structure descintf MIX imix FROM imixold USE TIMAV BURN INST BURN AVG EX BURN ivarty SET DELTA ADD LINEAR CUBIC PARKEY vall MAP 4 val2 MAP REF PARKEY valref SAMEASREF ENDREF MICRO ALL ONLY HISO conc ENDMIX where MIX keyword used to set imix Discontinuity factor infor
68. the branching generation of the PMAXS file keyword used to indicate that the meshing is the one used in the SAP or MCO object Default option IGE 344 USER GLOBAL ADD NEW pkey nval val DEF ADF DRA GET 103 keyword used to indicate that the meshing is defined by the user keyword used to set a global meshing by defining for each desired state variables a number of points for the branching calculation keyword used to add a set of new points for the branching calculation The new points are added to the meshing contained in the SAP or MCO object keyword used to indicate that the points contained in the SAP or MCO object are ignored conseguently only the set of points indicated using ADD will be considered for the branching calculation name of the state variable If pkey does not correpond to any name in the SAP or MCO object it will be ignored It is not necessary to set all state variable contained in SAP or MCO if a state variable is missing the SAP or MCO meshing for this state variable will be considered NB the BARR parameter cannot be modified by the user number of points for the state variable pkey In the case GLOBAL the nval points are obtained by splitting the pkey range from the first to the last values contained in the SAP or MCO object otherwise it corresponds to the number of new points to be introduced in the meshing value to be added in the branching calculati
69. the reactor core are updated accordingly Note that this option is valid only if a read only POWOLD object is provided keyword used to indicate the printing on files Note that all produced files will have the same extension res keyword used to specify the printing of the average fluxes and flux ratios per fuel bundle The normalized bundle fluxes are computed and printed for each reactor channel and per each energy group The flux ratios are computed with respect to the thermal energy group fluxes they are printed on the same file keyword used to indicate the printing of data computed over the whole reactor geom etry keyword used to specify the printing of flux distribution The normalized fluxes are printed in separated files one file per energy group the number of produced files will then equal to the total number of energy groups The flux values are printed for each mesh splitted volume in X Y and Z planes the virtual regions will have the fluxes values set to 0 keyword used to specify the printing of flux ratio distribution The flux ratios are com puted with respect to the thermal energy group fluxes per each mesh splitted volume They are printed in separated files the number of produced files will equal to the total number of energy groups less one keyword used to specify the printing of power distribution The power values are printed for each mesh splitted volume in X Y and Z planes the non fuel regions w
70. to different types of interpolation module used for the core voiding simulations module used to manage a full reactor execution in DONJON using explicit DRAGON calculations for each cell see Section 3 18 pin power reconstruction module l 6 7 IGE 344 2 2 Data structures The transfer of information between the modules is performed by means of well defined data struc tures also called objects The objects can be defined in either create read only or modification mode Each object has its own specific signature that can be easily recognized by a module A detailed descrip tion of DONJON data structures is given in Section 8 For more details on DRAGON and TRIVAC data structures refer to their guidel A brief description of all data structures that can be used in DONJON is given below GEOMETRY MACROLIB COMPO MULTICOMPO SAPHYB FMAP MATEX DEVICE DETECT TRACK SYSTEM FLUX POWER HISTORY data structure containing the geometry information This object has a signature L_GEOM it is created using DRAGON module GEO data structure containing the multigroup macroscopic properties it has a signature L_MACROLIB This object can be created in several modules namely using DRAGON modules MAC and NCR or using DONJON modules CRE MACINI and NEWMAC data structure containing the mono parameter database generated by the lattice code This object has a signature L_COMPO it is created usi
71. vc0 gt gt SET C lt lt vc0 gt gt ADD C lt lt vc0 gt gt MAP ADD C lt lt vc0 gt gt lt lt vc gt gt REF A lt lt va0 gt gt REF PA lt lt va0 gt gt 7B SAMEASREF ENDREF B SAMEASREF ENDREF SET A lt lt va0 gt gt ENDMIX ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B SAMEASREF ENDREF ENDMIX continued on next page IGE 344 MACROLIB NMIX 1 TABLE SAP MIX 1 SET A SET G ADD A REF ADD C REF ENDMIX PLANE AXE Fig MACROLIB 10 NMIX 1 TABLE SAP MIX 1 SET A SET G ADD A REF ENDMIX SCR SAP FMAP B lt lt va0 gt gt lt lt vc0 gt gt lt lt va0 gt gt lt lt va gt gt 2C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF lt lt vc0 gt gt lt lt vc gt gt PA lt lt va0 gt gt B lt lt vb0 gt gt ENDREF SCR SAP FMAP B lt lt va0 gt gt lt lt vc gt gt lt lt va0 gt gt lt lt va gt gt C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF 87 only the burnup in MAP MACROLIB NMIX 1 TABLE SAP B MIX 1 SET A SCR SAP FMAP lt lt va0 gt gt SET C lt lt vc0 gt gt ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF ADD C lt lt vc0 gt gt MAP REF A lt lt va0 gt gt B lt lt vb0 gt gt ENDREF ENDMIX MACROLIB NMIX 1 TABLE SAP
72. were provided at the database computation time keyword used to set tmod moderator temperature in K If this data is omitted the reference value in the data base is used The reference value is 345 66 K if none were provided at the database computation time keyword used to set nB Boron concentration in ppm If this data is omitted the reference value in the data base is used The reference value is 0 0 ppm See note below for inside equations keyword used to set dcool coolant density in g cm If this data is omitted the reference value in the data base is used The reference value is 0 81212 g cm if none were provided at the database computation time keyword used to set dmod moderator density in g cm3 If this data is omitted the reference value in the data base is used The reference value is 1 082885 g cm if none were provided at the database computation time keyword used to set purity moderator purity in atm If this data is omitted the reference value in the data base is used The reference value is 99 911 atm if none were provided at the database computation time keyword used to set bval This option is valid only when MCR is used and can not be omitted fuel burnup in MWd t This value must be positive keyword used to set nXe Xenon concentration in 10 at cm This concentration will be applied to every bun dle keyword used to specify that the Xenon concentrations as computed with DRAGON wil
73. will appear at four different i indices 3 1 2 Input of global and local parameters The information with respect to the fuel burnup is required for the fuel map MACROLIB construction using either the CRE NCR or AFM module The fuel region properties related to other local or global parameters can be interpolated only using the NCR module IGE 344 12 Table 4 Structure descresini2 EDIT iprint BTYPE TIMAV BURN INST BURN TIMAV BVAL bvalue i i 1 ncomb INST BVAL SAME bvalue CHAN bvalue i i 1 nch BUND bvalue i i 1 nch nb SMOOTH ASBLY bvalue i i 1 na OLDMAP BUNDLE POW SAME bvalue CHAN bvalue i i 1 nch BUND bvalue i i 1 nch nb REF SHIFT ishift COMB ishift i i 1 ncomb ADD PARAM PNAME PNAME PARKEY PARKEY GLOBAL LOCAL SET PARAM PNAME pvalue OLDMAP TIMES PNAMEREF SAME pvalue CHAN pvalue i i 1 nch BUND pvalue i i 1 nch nb y FUEL WEIGHT ENRICH POISON fvalue i i 1 nfuel CELL ialch i i 1 nch where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value 2 to print the channels refuelling schemes if they are new or modified 3 initial burnup limits per each channel are also printed if the axial power shape has been reinitialized BTY
74. 0 NZASS 122 nzass 122 OLDMAP 12 14 ONLY 69 70 80 81 OPTIM 112 114 118 optimize 148 149 OPTIMIZE 114 118 or 30 33 53 out 18 OUT STEP EPS 115 116 OUT STEP LIM 114 115 P NAME 32 P NEW 37 38 IGE 344 PARAB 52 PARABOLIC 115 116 PARAM 129 130 PARAMBURNTAR 60 61 PARKEY 12 13 69 70 80 81 PARKEY 12 13 part 104 pcore 31 32 PENAL METH 115 period 58 PHASE 100 101 pitch 104 pitch xbe ybe 101 PKEY 101 102 152 pkey 103 PKEY_INFO 153 PKEY 152 pkey i 103 PKEY_INFO 153 PLQ 118 plq_data 118 plq_data 118 PMAXS 102 PNAME 12 14 46 48 53 54 108 111 PNAME 12 13 48 111 pname 121 PNAMEREF 12 14 POISON 12 14 POOL 43 44 POROS 108 109 poros 108 109 pos 20 21 23 24 26 27 POSITION 23 24 poutlet 108 109 pow 88 89 122 POWER 36 37 39 42 43 45 47 49 POWER 36 38 122 power 138 140 POWER 5 36 39 42 45 49 power 37 38 58 POWOLD 36 38 PPR 120 122 PRINT 37 38 PROJECTION 121 PTOT 37 38 PUFR 108 109 pufr 108 109 PUR 89 90 purity 89 90 PUT 56 57 59 pvalue 12 14 46 48 108 111 121 PWF 88 89 PWR CASE UOX MOX CORE FUEL 105 194 QLP 149 QMAP 45 47 QUARTER 53 54 r1 108 109 r2 108 109 r3 108 109 r4 108 109 RADIUS 108 109 RATE 26 27 rate 26 27 RATIO 37 38 rburn 43 46 47 RC 94 95 RD 95 RDCL 89 90 RDMD 89
75. 05 106 107 108 109 110 111 112 113 114 115 Main records and sub directories in OLD VALUE LL 150 Records and sub directories in d2p_info data structure o oo o 152 Records and sub directories in SAPHYBINFO LL 152 Recordsin SPRENINEDF y era sole aa aaa a oe na 153 Records and sub directories in HELIOS_HEAD LL 153 Records and sub directories in GENPMAXS_INP o o oo 154 Records in sub directory TH DATA o oo ooo o 154 Records and sub directories in sub directory BRANCH_INFO 155 Sub directories in CROSSSECT osca be o a ee 155 Records in the sub directory MACROLIB_XS LL 156 Records in the sub directory MICROLIB_XS LL 156 Records in the sub directory DIVERS LL 157 IGE 344 vill List of Figures 1 Presentation of fully and partially inserted 3 part control rods a oaoa aaa 19 2 Complete grid one point case LL 76 3 Complete grid TA Case o ooo a a ee ee a ao 76 4 Partial grid complete planes one point case ooo a ee 76 5 Partial grid complete planes TA case oo oo ocea torota ana a ee 76 6 Partial grid complete axis one point case a ee 76 T Partial grid complete axis TA case kas 76 8 Partial grid complete axis with another configuration one point case 77 9 Partial grid one complete plane and one complete axis one point case 77 10 Partial grid one complete plane and one complet
76. 05 470 535 DMIX lt lt mSORin gt gt lt lt mSDRout gt gt ENDROD ROD 51 ROD NAME SOR21 LEVEL 0 0 AXIS Y FROM H MAXPOS 117 0 139 0 53 75 466 25 495 3 544 83 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 52 ROD NAME SOR22 LEVEL 0 0 AXIS Y FROM H MAXPOS 183 0 205 0 24 5 495 5 495 3 544 83 DMIX lt lt mSORin gt gt lt lt mSDRout gt gt ENDROD ROD 53 ROD NAME SOR23 LEVEL 0 0 AXIS Y FROM H MAXPOS 227 0 249 0 24 5 495 5 495 3 544 83 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 54 ROD NAME SOR24 LEVEL 0 0 AXIS Y FROM H MAXPOS 271 0 293 0 24 5 495 5 495 3 544 83 DMIX lt lt mSORin gt gt lt lt mSORout gt gt IGE 344 x ENDROD ROD 55 ROD NAME SOR25 LEVEL 0 0 AXIS Y FROM H MAXPOS 315 0 337 0 24 5 495 5 495 3 544 83 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 56 ROD NAME SOR26 LEVEL 0 0 AXIS Y FROM H MAXPOS 381 0 403 0 53 75 466 25 495 3 544 83 DMIX lt lt mSORin gt gt lt lt mSDRout gt gt ENDROD create rod devices groups x END CREATE ROD GR 5 GROUP ID 1 ROD ID 1 2 15 16 19 20 23 24 27 28 GROUP ID 2 ROD ID 3 4 7 8 9 10 13 14 GROUP ID 3 ROD ID 5 6 11 12 17 18 21 22 25 26 29 30 GROUP ID 4 ROD ID 31 36 37 38 39 42 43 44 45 48 49 50 51 52 53 54 55 GROUP ID 5 ROD ID 32 33 34 35 40 41 46 47 53 56 QUIT 9 3 3 Input file for fuel map Input data for test case Pfmap c2m DK KKK KK K KK ARE ARE K K K AE
77. 1 ABS nsh ALL imix SHUFF CHAN NMCHA1 TO NMCHA2 POOL J where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value 2 only the burnup limits over each channel are printed 3 only the axial power shape values over each channel are printed 4 only the channel refueling rates are printed for larger values of iprint everything will be printed BURN STEP keyword used to indicate an increase of core average burn up rburn keyword used to indicate in MWd t the average increase of burn up in the core TIME keyword used to indicate the time of combustion at the power specified in POWER structure rtime keyword used to set the time combustion value in DAY or HOUR or MINUTE or SECOND DAY keyword used to specify that rtime is a number of days HOUR keyword used to specify that rtime is a number of hours MINUTE keyword used to specify that rtime is a number of minutes SECOND keyword used to specify that rtime is a number of seconds REFUEL key word to specify a channel refueling MICRO keyword used to perform a microscopic refueling In this case three libraries have to be provided when TINST is called CHAN key word to specify the refueled channel information NAMCHA channel name as defined by NXNAME and NYNAME NAMCHA is a character 4 variable constructed as WRITE NAMCHA A1 A3 NYNAME 1 1 NXNAME 1 2 nsh refuel
78. 10 nB 89 90 nb 11 12 46 nb1 108 111 nb2 108 111 nbf 108 109 nbfuel 67 68 nbg 108 109 nbun 57 58 nburn 94 95 Nc 59 nch 9 11 12 ncha 57 58 ncols 104 ncols nrows part hm_dens bypass 101 NCOMB 9 11 ncomb 9 11 12 ncond 108 110 ncorner 105 NCR 67 68 NCR 4 5 11 12 17 40 42 43 53 54 63 67 ncr_data 67 69 ncr_data 67 nest 114 117 ndetect 22 23 NEP 89 90 NEW 46 101 103 NEWFUEL 43 NEWMAC 4 5 17 21 27 35 40 133 NFTOT 83 NFUEL 15 16 nfuel 12 14 16 53 NG 82 ngcond 94 95 nglo 57 58 nglob 58 NGMTR 94 ngr1 112 113 ngr2 112 113 NGREAC 94 NGRP 15 16 22 ngrp 15 16 18 19 21 22 26 28 30 34 NHEX 23 24 nhex 23 24 nk 43 nloc 57 58 nlzc 26 28 30 NMCHA1 43 44 NMCHA2 43 44 NMIX 15 16 64 67 68 78 79 94 193 nmix 64 65 123 nmixt 15 16 67 68 78 79 94 nmxgeo 123 124 nNp 89 90 NO STORE OLD 118 noreg 94 95 NORM 37 38 51 52 NP 83 npert 94 96 NPWF 88 89 NREF 89 NREFL 15 16 nrefl 15 16 nrep 22 23 NRMFLUX 36 nrod 18 21 30 34 nrows 104 nsh 43 44 nside 104 nside ncorner vfcm 101 nSm 89 91 NTOT1 112 113 ntyp 88 89 NUM LZC 26 NUM ROD 18 19 nval 103 nvar 114 116 nvoid 53 55 nx 9 10 nxa 10 nxass 123 nXe 89 90 NXNAME 43 44 NXNAME 9 10 ny 9 11 nya 10 nyass 123 NYNAME 43 44 NYNAME 9 1
79. 108 109 FLPOW 4 5 36 37 39 40 42 45 49 51 89 138 FLUD 3 5 36 37 40 51 122 IGE 344 FLUNAM 120 FLUX 36 37 51 FLUX 37 38 flux 140 FLUX 5 36 37 51 114 FLUXUNK 120 FMAP 36 39 42 45 53 54 63 65 FMAP 8 13 fmap 125 128 FMAP 5 36 39 42 45 53 63 64 FMAPYV 53 54 fmax 26 27 FMIX 15 16 fname 101 105 FOBJ CST VAL 115 116 FORCEAVE 108 111 FPUISS 107 108 frac 23 fract 22 107 108 FRACTION 22 23 FROM 20 45 47 69 80 82 83 frstrec 94 95 FSTH 37 38 fsth 37 38 FT 95 FUEL 12 14 46 47 101 102 128 129 fuel 47 FUELDIR 141 FULL 53 54 FULL MIX 26 27 funct 115 116 fvalue 12 14 GEN 100 105 GENPMAXS 101 105 GENPMAXS_INP 154 GEO 3 5 9 10 15 51 GEOM 15 16 51 GEOMAP 128 geometry 125 127 128 GEOMETRY 5 15 16 51 120 124 GEOMOLD 15 16 GEONEW 120 GEOOLD 120 GET 56 57 59 101 103 153 GFF 123 GLOBAL 12 13 57 58 101 103 GPPR 122 GPT 114 GPT 114 GRAD 114 grad_data 114 grad_data 114 115 190 GREP 3 GRID 101 102 GROUP 33 34 GROUP ID 21 28 GRPMAX 112 113 GRPMIN 112 113 H 20 H 20 HALF 53 54 115 116 hel 46 47 he2 46 47 hcase 45 46 hcnew 45 46 hcold 45 47 hcold2 46 47 HCONV 108 110 hconv 108 110 height 10 HEL 100 104 105 HELIOS 101 104 HELTOS dra 105 HELIOS_HEAD 153 154 HEX 23 24
80. 14 24 36 44 uonmoniwns Nn Y OTVOZ ZTOY 247 29 218 23 22 14 13 23 35 45 65 210 21 22 34 46 16 22 21 33 46 180 IGE 344 END QUIT 9 3 4 Input file for exit burnups AE popapepe pap pe a Me B alel efio njijijnjiolo s 7zfs lilejslals slalsjela Figure 15 Combustion zones definition 65 68 67 ss 57 s8 47 48 49 38 39 40 28 30 91 ale 43 14 ls ae 181 IGE 344 Input data for test case Pburn c2m Lu KKK K K oo K RE FK K OK K dd FK K kk ok CALL XK XA XA XX XX X Procedure Purpose Author s Pburn c2m Provide average exit burnups D Sekki 2007 11 FMAP Pburn FMAP gt XK Xx XA XX XX XX X FARR k kak ak k kak 3K k kak I k 3K 3K AO AA AAA K K FKK K FK FK K 3K FKK K FK PARAMETER FMAP MODULE FMAP BTYPE TIMAV BURN TIMAV BVAL A 1008 B 1365 2014 C 1806 3250 D 2279 3409 E 2258 3456 F 2678 3857 G 2792 4025 H 1693 4134 J 2268 4121 98 60 25 12 93 49 15 14 56 92 62 21 20 96 99 81 05 1057 1586 2364 2019 3391 2267 3715 2345 48 3705 3263 3925 3087 3978 2811 4025 3065 3995 69 43 85 33 85 03 98 30 24 74 13 10 67 20 97 54 RESINI FMAP 1071 1643 2425 47 3761
81. 2 Structure DETINI DETECT DETINI DETECT descdet where DETECT character 12 name of the DETECT object that will be created by the module it will contain the detector informations If DETECT appear on RHS it is updated otherwise it is created descdev structure describing the input data to the DETINI module 3 5 1 Input data to the DETINI module Note that the input order must be respected Table 13 Structure descinidet EDIT iprt HEXZ NGRP ngrp TYPE NAMTYP INFO ndetect nrep SPECTRAL spec i i 1 ngrp DEFAULT INVCONST tinv i i 1 nrep 2 FRACTION fract i i 1 nrep 1 descdet i 1 ndetect where EDIT keyword used to set iprt iprt index used to control the printing in module INIDET 1 2 for no print default value 3 for printing the contents of the output DETECT HEXZ keyword to specify that only hexagonal detectors will be defined If this keyword is absent Cartesian detectors will be defined NGRP keyword used to set ngrp ngrp number of energy groups in the calculation It must be equal to the number set in the MACD module or by the COMPO files IGE 344 23 TYPE keyword to specify the detector type NAMTYP character 12 name of the detector type To correspond to the actual detector re sponse model encoded the type of detector must be in this list e PLATN_REGUL PLATN_SAU e VANAD_REGUL CHION_SAU CHION_REGUL For other type
82. 20 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 5 ROD NAME ZCRO3A LEVEL 0 0 AXIS Y FROM H MAXPOS 249 0 271 0 0 0 260 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 6 ROD NAME ZCRO3B LEVEL 0 0 AXIS Y FROM H MAXPOS 249 0 271 0 260 0 520 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 7 ROD NAME ZCRO4A LEVEL 0 0 AXIS Y FROM H MAXPOS 293 0 315 0 0 0 260 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 8 ROD NAME ZCRO4B LEVEL 0 0 AXIS Y FROM H MAXPOS 293 0 315 0 260 0 520 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 9 ROD NAME ZCRO5A LEVEL 0 0 AXIS Y FROM H MAXPOS 337 0 359 0 0 0 260 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt IGE 344 172 ENDROD ROD 10 ROD NAME ZCRO5B LEVEL 0 0 AXIS Y FROM H MAXPOS 337 0 359 0 260 0 520 0 123 825 173 355 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 11 ROD NAME ZCRO6A LEVEL 0 0 AXIS Y FROM H MAXPOS 161 0 183 0 0 0 260 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 12 ROD NAME ZCRO6B LEVEL 0 0 AXIS Y FROM H MAXPOS 161 0 183 0 260 0 520 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 13 ROD NAME ZCRO7A LEVEL 0 0 AXIS Y FROM H MAXPOS 205 0 227 0 0 0 260 0 272 415 321 945 DMIX lt lt mZCRi
83. 2222 2222222000 0022222222 2222222200 A B Cc IGE 344 NX 1 10 17 23 29 32 o ooo NNNNNO NNNNNDN NNNNNDN NNNNINN NNNNNN NNNNNN NNNNNN NNNNNN NNNNNN 6 O O o 5 5 0 5 oo NINNNN o Oo NNNNNNN O ONNNNNNNNIN DONNNNNNNN ONNNNNNNNN ONNNNNNNNN ONNNNNNNNN ONNNNNNNNN ONNNNNNNNN PLANE PLANE PLANE PLANE PLANE PLANE PLANE MESHX MESHY MESHZ 0 0 194 348 0 0 297 18 346 NAME NYNAME 2 11 18 24 30 31 3 3 2147 2 2 K 3 12 19 25 31 30 2 3 Hom M PPR RB fH ta a al 62 16 70 84 238 392 84 238 392 NNNNNN NNNNNN NNNNNN NNNNNN NNNNNN NNNNNN NNNNNN NNNNNN NNNNNO o oo ONNNNNNNNN ONNNNNNNNN ONNNNNNNNDN ONNNNNNNNN ONNNNNNNNN DONNNNNNNN O ONNNNNNNN DOONNNNNNN DOOODOONNNNN Oo0DO0OO0O0O0OOoOoOoOoOo 106 260 414 106 260 414 49 53 99 06 71 396 24 445 77 495 30 544 83 594 1 21 2p L 4 12 20 26 32 29 22 2 B 3 73 213 2 20 M N NCOMB 73 B ZONE 5 11 20 27 33 38 5 10 19 28 34 39 24 214 D 20 18 28 35 40 128 282 436 128 282 436 148 252 215 E P w 17 27 36 41 59 26 2 F Q N 16 26 37 42 150 304 458 150 304 458 172 326 520 12 326 520 198 12 2 x7 16 al 7G 2 R 15 25 37 43 28 Lif 7H 7S 13
84. 222222222 222222222444 444222222222 222222222444 A B Cc D E F G H K L M N 0 P Q R S 444222222222 222222222444 444222222222 222222222444 444222222222 222222222444 444222222222 222222222444 044222222222 222222222440 044422222222 222222224440 004422222222 222222224400 000442222222 222222244000 000444422222 222224444000 000004444444 444444400000 000000444444 444444000000 000000004444 444400000000 PLANE 6 SAME 5 PLANE 7 SAME 5 PLANE 8 SAME 5 IGE 344 PLANE 9 SAME 1 PLANE 10 SAME 1 PLANE 11 SAME 1 PLANE 12 SAME 1 MESHX 0 0 20 40 62 84 106 128 150 172 194 216 238 260 282 304 326 348 370 392 414 436 458 480 500 520 MESHY 0 0 20 40 62 84 106 128 150 172 194 216 238 260 282 304 326 348 370 392 414 436 458 480 500 520 MESHZ 0 0 49 53 99 06 148 59 198 12 247 65 297 18 346 71 396 24 445 77 495 30 544 83 594 36 SPLITX NIN NN NN NN NN N N NN NN NIN NN SPLITY NN NN NN NN NN NN NN NN NN NN SPLITZ 222222222222 END QUIT 9 3 2 Input file for devices Input data for test case Pdevc c2m DR ARE KKK KKK KK K FK K K K AE FKK FK EE dd K FK K 2K FK AA AA K K K K K K K K OK Procedure Pdevc c2m Purpose Reactor rod devices specification Author s D Sekki 2007 06 CALL DEVICE MATEX Pdevc MATEX lt lt mZCRin gt gt lt lt mZCRout gt gt lt lt mSORin gt gt lt lt mSORout gt gt XK XA Xx XX XK XA Xx
85. 4445 08 wee eee E 112 61 1 Data input for module DLEAK ok ee k k os 112 6 2 The GRAD module a ca eb dad A Oo hee tl ee oe 114 621 Data input for module GRAD o e e cocs c s sorea a tee RE EE wes 114 6 3 The PLS mode ug Deo Seay wees A a he e pe a ee we 118 Gi Data input for module PLQ a a k uoi a fl a 118 7 PIN POWER RECONSTRUCTION MODULES 120 Tal The NAP Module escran segad a ee AAA A 120 TL Additional properties calculations LL 121 TaL2 Pin power reconstruction i 122 LS Heterogeneous assembly geometry definition 123 8 DONJON DATA STRUCTURES s ee e pi a k ee ae 125 8 1 Contents of fmap data structure lt s ce ee 125 8 1 1 The state vector content i Lei pi 125 812 Tih inain fmap directora a a Be e S 126 8 1 3 The FUEL sub directories o o o 128 8 1 4 The hcycle sub directories o o oo e eee eee 129 8 1 5 The PARAM sub directories 2 o ee 129 8 1 6 The ASSEMBLY sub directory o o e e 130 8 2 Contents of matex data structure 2 00 ee 131 8 2 1 The state vector content 131 8 2 2 The matex directory 4 cio b ra ee DL 131 8 3 Contents of device data structure ks 133 8 3 1 The state vector content 2 ke A ee 133 8 3 2 The main device directory lt c coos pacte oo toek ki i 133 8 5 3 The DEV ROD sub directories o e e 134 8 3 4 The ROD GROUP sub direct
86. 67 CPONAM2 67 CRE 4 5 11 12 17 40 43 63 64 CREATE 18 19 26 CRITFL 107 109 CROSS_SECT 155 CROSS_SECT 155 156 CST OBJ 115 116 CST QUAD EPS 115 116 CST TYPE 115 116 CST WEIGHT 114 115 117 cstval 115 117 cstw 115 117 CT 95 CUBIC 67 70 78 81 CVR 4 53 22 CWSECT 108 109 CYCLE 45 46 D2P PMAX 105 D2P 97 100 101 105 151 d2p info 151 152 DATABASE 88 DAY 43 46 47 DBASE 88 dcool 89 90 dcoolV 53 54 DECAY 82 83 DEF 101 103 DEFAULT 22 23 DELETE 3 DELH 33 34 delh 33 34 DELT 33 34 delt 33 34 DELTA 69 71 80 81 83 112 113 DENS COOL 53 54 der 101 105 DERIV 65 DERIVATIVE 101 105 descafm 88 89 desccre1 17 63 64 desccre2 17 63 64 descdatal 64 65 descdata2 64 65 descdepl 78 80 82 descdet 22 23 descdetect 51 descdev 18 22 descdset 29 descflpow 36 37 descgeo 9 10 descinidet 22 IGE 344 descintf 67 69 descintf 68 69 descints 78 80 descints 79 80 desclink 15 desclzc 25 26 descmcc1 31 descmove 33 descnap 9 descnap1 120 121 descnap2 120 122 descnap3 120 121 123 descrcvr 53 descresinil 8 9 descresini2 8 9 12 descsim 45 46 desct16cpo 93 94 96 desctavg 39 descthm 107 108 desctinst 42 43 descxenon 19 50 DETEC 51 DETECT 22 detect 137 DETECT 5 22 51 52 DETECT 4 5 36 38 51 52 137
87. 6cpo where DONCPO name of data structure where the output CPO is stored This can be a new data structure or an old data structure which will be updated desct16cpo input specifications for the execution of the T16CP0 module i end of record keyword This keyword is used to delimit the part of the input data stream associated the current module In the following dataset MODULE T16CPO SEQ_BINARY WIMS16 LINKED_LIST DONCPO DONCPO T16CPO WIMS16 means that that the module will read the sequential binary file WIMS16 file in readonly mode and create the CPO data structure DONCPO while the dataset MODULE T16CPO SEQ_BINARY WIMS16 LINKED_LIST DONCPO IGE 344 94 DONCPO T16CPO DONCPO WIMS16 means that the data structure DONCPO will be updated The input instructions replaced by here should indicate what part of the information located on WIMS16 should be transferred to DONCPO and in what order 4 5 1 Input data for the T16CP0 module The input data structure desct16cpo will take the form Table 67 Structure desct16cpo EDIT iprint NMIX nmixt CONDG ngcond igc i i 1 ngcond LIST MIX MIXNAM CELLAV REGION noreg RC nburn frstrec MTMD valreft valrefd npert valpert i valperd i frstrec i i 1 npert NAMPER valref npert valper i frstrec i i 1 npert 11 where EDIT optional keyword used to
88. By default the rod insertion level is left undefined value real positive value of the rod insertion level This value is used to compute the actual rod position in the reactor core The rod insertion level is minimal value 0 0 when the rod is completely withdrawn and it is maximal value 1 0 when the rod is fully inserted For the partially inserted rod the insertion level must be 0 0 lt value lt 1 0 SPEED keyword used to specify speed By default the speed is left undefined speed real positive value of the rod movement speed given in cm s This value is needed only for the reactor regulating purpose IGE 344 TIME time MAXPOS pos DMIX mix mix2 ENDROD 21 keyword used to specify time By default the insertion time is left undefined real value of time for the rod insertion or extraction given in sec This value is needed only for the reactor regulating purpose keyword used to specify the full inserted coordinates of a rod part The sequence of MAXPOS and DMIX data structures is repeated for each part making the rod real array containing 3 D Cartesian coordinates of the full inserted rod This is the limiting rod position in the reactor core which may or may not be the same as the actual rod position These coordinates must be given in the order X X Y Y Z and Z keyword used to specify mix1 and mix2 first of two integer rod mixture indices Index mix1 corresponds to the
89. CTER 2 4 General input structure DONJON is built around the GAN generalized driver 1 Accordingly all the modules that will be used during the current execution must be first identified It is also necessary to define the format of each object data structure that will be processed by these modules Then the modules required for the specific DONJON calculation are called successively information being transferred from one module to the next via the objects Finally the execution of DONJON is terminated when it encounters the END module even if it is followed by additional data records in the input data stream The general input data structure therefore follows the calling specifications given below Table 1 Structure DONJON MODULE MODNAME LINKED_LIST STRNAME XSM_FILE STRNAME SEQ_BINARY STRNAME SEQ_ASCII STRNAME module where MODULE keyword used to specify the names of all modules that will be used in the current DONJON execution MODNAME character 12 name of a DONJON or DRAGON or TRIVAC or utility module The list of modules that can be executed using DONJON code is provided in Section 2 1 IGE 344 7 LINKED LIST keyword used to specify the names of data structure that will be stored as linked lists XSM_FILE keyword used to specify the names of all data structure that will be stored on XSM format files SEQ_BINARY keyword used to specify the names of al
90. E control i i 1 nvar VAR WEIGHT weight i i 1 nvar VAR VAL MIN vecmin i i 1 nvar ALL varmin VAR VAL MAX vecmax i i 1 nvar ALL varmax FOBJ CST VAL funct i i 1 ncst 1 CST TYPE type i i 1 ncst CST OBJ cstval i i 1 ncst CST WEIGHT cstw i i 1 ncst keyword used to set iprint index used to control the printing in module keyword used to define the guasi linear programming method Note If the general Lemke method is used the guadratic constraint must be active The strategy consists to proceed in two steps e At first step the linear programming problem i e without the quadratic con traint is solved and the control variable displacement is computed If this dis placement is less than the radius of the guadratic constraint the step one solution is accepted and step two is not performed If this displacement is greater than the radius of the guadratic constraint the step one solution is rejected and step two is performed Step one can be solved with the SIMPLEX method or with the linear LEMKE method e At step two the general LEMKE method is used to find the correct solution The general Lemke method is based on a parametric linear complementarity principle as explained in Ref 32 keyword used to specify that the SIMPLEX method will be used at step one and the general LEMKE method at step two keyword used to specify that the linear LEMKE method will be
91. ETE module used to delete one or many data structures GREP module used to extract a single value from a data structure END module used to delete all the local linked lists to close all the remaining local files and to return from a procedure or to terminate the overall DONJON execution controlled by the GAN generalized driver e The following are short descriptions of DONJON modules IGE 344 CRE NCR AFM USPLIT RESINI MACINI DEVINI DETINI LZC DSET MCC MOVDEV NEWMAC FLPOW TAVG TINST SIM DETECT CVR HST NAP 4 module used to create a MACROLIB containing the material properties by interpolating the nuclear properties from a mono parameter database previously generated in the lattice code module used to create a MICROLIB or a MACROLIB containing the material properties by interpolating the nuclear properties from a multi parameter database previously generated in the lattice code module used to create a MACROLIB containing the material properties by interpolating the nuclear properties from a multi parameter feedback model database previously generated in the lattice code module used to create an extended reactor material index over the whole mesh splitted reactor geometry module used to define the fuel lattice to create the fuel map geometry and to specify the global and local parameters module used to create an extended MACROLIB in whi
92. FKK FK K FK FK FK FKK FKK FK K FK FK FK FK K FK K FK K FK FK K FK K FK FK K K K FK K K K K 2K K K OK XK XA Xx XX XX X Procedure Pfmap c2m Purpose Reactor fuel map specification Author s D Sekki 2007 11 CALL FMAP MATEX Pfmap MATEX x Xx Xx XX XX X 178 179 IGE 344 Lu KKK KK K FKK K FKK OK K K FKK FK K FK FK FK FK FK FK FK FK K FK FK FK FK K FK K FK K K FK FK FK K FK FK K K K FK K FK K K 2K K K ok PARAMETER FMAP MATEX gt LINKED_LIST FMAP MATEX END RESINI MODULE EDIT 0 RESINI MATEX FMAP MATEX GEO CAR3D 20 20 12 EDIT 0 X VOID X VOID Y VOID Y VOID Z VOID Z VOID MIX PLANE 1 123456789 012345678 0000000000 0000000000 0000011111 A B C D E F G H J 1111100000 1111111000 111111100 1111111100 1111111110 1111111110 1111111110 e ie es i LELLO 1111111110 0001111111 0 0 Led 1 14111 1 0011111111 001111141111 0111111111 0111111111 0 171 EL LL dd dl 0111111111 K L M N 0 P Q R S 11111111 11 10 VELIA LA L O 1111111110 1111111110 1111111110 1111111100 1111111100 1111111000 1111100000 0111111111 O At 11d LL 0111111111 01111141111 034 1 1 1 1 11 1 1 0011111111 0011111111 0001111111 0000011111 0000000000 0000000000 PLANE 2 SAME 1 PLANE 3 SAME 1 PLANE 4 SAME 1 PLANE 5 123456789 012345678 0000000000 0000000000 0000022222 2222200000 000222
93. Groups Set of data for FILE_CONT_3 block in DRA file Volume of coolant VCoo1 moderator VModr control rods VCnRd fuel VFuel cladding VClad channels VChan water VWatR Set of data for FILE_CONT_4 block in DRA file Cell Pitch and X Y Position of First Cell continued on next page IGE 344 Records and sub directories in HELIOS_HEAD XS_ CONTunuvua FLAG uuuuuua JOE TTT II FILE NAME jj DERIVATIVE VERSION ui COMMENT uuu SOB OPT iis DAT SRC jain Type R 1 Condition Units 154 continued from last page Comment Set of data for XS_CONT block in DRA file VFCM Table 109 Records and sub directories in GENPMAXS_INP Type Condition Units Comment Indicates the end of a branch calculation Name of final PMAXS file Name of output DRA file Set of data for FILE_CONT_3 block in DRA file Volume of coolant VCool moderator VModr control rods VCnRd fuel VFuel cladding VClad channels VChan water VWatR Version of PARCS used Free comment Options for GenPMAXS running see Data source information see 10 Each component of the list TH_DATA is a directory containing TH data for a single burnup point Inside each groupwise directory the following records associated will be found CHI ou uu OVER V uuuuuu LAMBDA juuuuu BETA y uvuuuuua YLDPm uuuuuva YLDXeLuuuuuu YLDI uuuuuuva Type Table 110 Records in sub direct
94. ISOTOPES XE135 xenam SM149 smnam 1135 inam ENDISOTOPES REFLECTOR HELIOS FILE_CONT_1 ncols nrows part hm_dens bypass FILE_CONT_2 emin g g 1 ngroup FILE_CONT_3 vcool vwatr vmodr venrd vfuel vclad vchan FILE_CONT_4 pitch xbe ybe XS_CONT nside ncorner vfcm GENPMAXS JOB_TIT jobtit FILE NAME fname DERIVATIVE der VERSION vers COMMENT comment JOB_OPT ladf Ixes Ided Ij1f Ichi Ichd linv Idet lyld Icdf Igff lbet lamb Idec IUPS iups SFAC sfac BFAC bfac IGE 344 where NAMDIR mixdir MIX imix PKEY refnam sapnam FUEL BARR unrodded aicg aicn compo ISOTOPES XE135 xenam SM149 smnam 1135 inam GRID SAP 102 keyword used to set mixdir name of sub directory in Multicompo containing information to be recovered Default mixdir default keyword used to set imix index of the mixture in the SAP object which will be considered by the module Default imix 1 keyword used to associate a name of PKEY in the SAP object to a type of state variable type of state variable The possible refnam are BARR for the control rod DMOD for the density of coolant g cc CBOR for the boron concentration ppm TCOM for the fuel temperature C TMOD for the moderator temperature C BURN for the burnup exposure MWd T It is not necessary to specify all state variable names only state variables with a different nam
95. IT O ACCE 3 3 ADI 4 EXTE 1000 lt lt Precf gt gt THER 1000 ENDIF SYSTEM DELETE SYSTEM flux and power POWER FLPOW FMAP FLUX TRACK MATEX EDIT 0 PTOT lt lt Power gt gt burnups integration limits FMAP TAVG FMAP POWER EDIT O AX SHAPE RELAX 0 5 B EXIT POWER DELETE POWER dini current parameters IGE 344 166 GREP FLUX GETVAL K EFFECTIVE 1 gt gt Keff lt lt GREP FMAP GETVAL EPS AX 1 gt gt Eps lt lt ECHO Iteration No iter ECHO AXIAL SHAPE ERROR Eps ECHO RESULTING K EFF Keff ENDWHILE edit resulting fluxes and powers pe POWER FLPOW FMAP FLUX TRACK MATEX EDIT lt lt iEdit gt gt PTOT lt lt Power gt gt os last parameters Ho GREP FLUX GETVAL K EFFECTIVE 1 gt gt Keff lt lt GREP FMAP GETVAL EPS AX 1 gt gt Eps lt lt GREP FMAP GETVAL B EXIT 1 gt gt Bexit lt lt ECHO Number of Iterations iter ECHO AXIAL SHAPE ERROR Eps p ECHO CORE AVERAGE EXIT BURNUP Bexit ECHO RESULTING K EFFECTIVE gt M Keff assertS FLUX K EFFECTIVE 1 1 050102 END QUIT 9 3 Procedures 9 3 1 Input file for geometry Input data for test case Pgeom c2m FE GAA kak kak k I IG kkk kkk LI ICA ACI A A I k kkk k ak Procedure Pgeom c2m Purpose Reactor geometry definition Author s D Sekki 2007 06 CALL GEOM Pgeom XK XA Xx XX X
96. IX keyword used to set imix Discontinuity factor information present in the Saphyb is interpolated as mixture 1 values imix index of the mixture that is to be created in the MICROLIB and MACROLIB FROM keyword used to set the index of the mixture in the SAPHYB object imixold index of the mixture that is recovered in the SAPHYB object By default imixold 1 USE keyword used to set the index of the mixture in the SAPHYB object equal to imix TIMAV BURN keyword used to compute time averaged cross section information This option is avail able only if a MAPFL object is set By default the type of calculation TIMAV BURN or INST BURN is recovered from the MAPFL object INST BURN keyword used to compute cross section information at specific bundle burnups This option is available only if a MAPFL object is set By default the type of calculation TIMAV BURN or INST BURN is recovered from the MAPFL object AVG EX BURN keyword used to compute the derivatives of cross section information relative to the exit burnup of a single combustion zone The derivatives are computed using Eq 3 3 of Ref 16 written as 0 1 Bie e pe eoc oc o dB De B ik Er Ik _ B el re OBS Be Beye Bes f j j i j b Bie where BR 5 and B are the beginning of cycle burnup of bundle j k end of cycle burnup of bundle j k and exit burnup of channel j This option is available IGE 344 ivarty SET DELTA ADD LINEAR CUBIC
97. K DMACRO OPTIM DLEAK MACRO dleak_data where DMACRO character 12 name of a LCM object type L_MACROLIB containing the delta MACROLIB information DMACRO is created by the module A STEP heteroneneous list is present in DMACRO OPTIM character 12 name of a second LCM object type L_OPTIMIZE created by the module Leakage related parameters are saved in the the control variable record VAR VALUE of OPTIM object Input data defined in Sect 6 1 1 is also saved in OPTIM object MACRO character 12 name of the LCM object type L_MACROLIB containing the input MACROLIB dleak_data structure containing the data to module DLEAK see Sect 6 1 1 6 1 1 Data input for module DLEAK Table 73 Structure dleak_data EDIT iprint TYPE DIFF NTOT1 DELTA VALUE FACTOR MIXMIN ibm MIXMAX ibm2 GRPMIN ngr1 GRPMAX ngr2 3 IGE 344 where EDIT iprint TYPE DIFF NTOT1 DELTA VALUE FACTOR MIXMIN ibm1 MIXMAX ibm2 GRPMIN ngrl GRPMAX ngr2 113 keyword used to set iprint index used to control the printing in module DLEAK keyword used to set the leakage parameter that is differentiated differentiation with respect to diffusion coefficients differentiation with respect to P weighted macroscopic total cross sections keyword used to set the type of differentiation differentiation with respect to the leakage parameter itself differentiation with respect to th
98. MAP 20 lt lt vc0 gt gt B SAMEASREF ENDREF SCR SAP FMAP B lt lt va0 gt gt lt lt vc0 gt gt lt lt va0 gt gt MAP 20 lt lt vc0 gt gt 7B lt lt vb0 gt gt ENDREF lt lt vc0 gt gt MAP gt A lt lt va0 gt gt 7B lt lt vb0 gt gt ENDREF continued on next page IGE 344 86 all parameters explicitly set all parameters in MAP PLANE o Fig MACROLIB SCR SAP MACROLIB SCR SAP FMAP NMIX NMIX 1 SAPHYB SAP TABLE SAP B MIX 1 MIX 1 SET A lt lt va0 gt gt SET A lt lt va0 gt gt SET B lt lt vb gt gt ADD A lt lt va0 gt gt MAP SET C lt lt vc gt gt REF C lt lt vc0 gt gt ADD A lt lt va0 gt gt lt lt va gt gt B lt lt vb0 gt gt ENDREF REF C lt lt vc0 gt gt ENDMIX B lt lt vb0 gt gt ENDREF ENDMIX Table 61 SCR inputs for instantaneous cases For the TA the burnup variable has no other choice than to be stored in the MAP object MAPFL Then the input files will be only the burnup in MAP all parameters in MAP MACROLIB SCR SAP FMAP MACROLIB SCR SAP FMAP NMIX 1 NMIX 1 TABLE SAP B TABLE SAP B MIX 1 MIX 1 SET A lt lt va gt gt ENDMIX SET C lt lt vc gt gt ENDMIX MACROLIB SCR SAP FMAP MACROLIB SCR SAP FMAP NMIX 1 NMIX 1 TABLE SAP B TABLE SAP B MIX 1 MIX 1 SET A lt lt va gt gt SET C lt lt
99. MAXPOS 271 0 293 0 24 5 495 5 198 12 247 65 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 42 ROD NAME SOR12 LEVEL 0 0 AXIS Y FROM H MAXPOS 381 0 403 0 53 75 466 25 222 885 272 415 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 43 ROD NAME SOR13 LEVEL 0 0 AXIS Y FROM H MAXPOS 117 0 139 0 53 75 466 25 272 415 321 945 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 44 ROD NAME SOR14 LEVEL 0 0 AXIS Y FROM H MAXPOS 381 0 403 0 53 75 466 25 272 415 321 945 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 45 ROD NAME SOR15 LEVEL 0 0 AXIS Y FROM H MAXPOS 117 0 139 0 53 75 466 25 321 945 371 475 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 46 ROD NAME SOR16 LEVEL 0 0 AXIS Y FROM H MAXPOS 227 0 249 0 24 5 495 5 346 71 396 24 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 47 IGE 344 177 ROD NAME SOR17 LEVEL 0 0 AXIS Y FROM H MAXPOS 271 0 293 0 24 5 495 5 346 71 396 24 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 48 ROD NAME SOR18 LEVEL 0 0 AXIS Y FROM H MAXPOS 381 0 403 0 53 75 466 25 321 945 371 475 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 49 ROD NAME SOR19 LEVEL 0 0 AXIS Y FROM H MAXPOS 117 0 139 0 53 75 466 25 421 005 470 535 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 50 ROD NAME SOR20 LEVEL 0 0 AXIS Y FROM H MAXPOS 381 0 403 0 53 75 466 25 421 0
100. MULTICOMPO object or set using the HISO conc data statement keyword to indicate that only the isotopes set using the HISO conc data statement will be used in the MICROLIB and MACROLIB objects character 8 name of an isotope 024 user defined value of the number density in 1024 particles per cm of the isotope IGE 344 71 the value of the number density for isotope HISO is recovered from the MULTICOMPO object ENDMIX end of specification keyword for the material mixture 4 2 3 Interpolation in the parameter grid The following example corresponds to a delta sigma computation in mixture 1 corresponding to a perturbation Note that in this case the MACROLIB object may content negative cross section MACROLIB NCR CPO EDIT 40 NMIX 1 MACRO COMPO CPO default MIX 1 delta sigma contribution SET CELL 3D DELTA PITCH 0 0 1 0 ENDMIX When the number of parameters used for the interpolation is increased all the lattice computations corresponding to all the combinations of parameters may not be done for computation time reasons In this case some approximations may be required The choice for the SET DELTA and ADD is then dependent of the structure of the database i e how the database grid of possibilities is filled When a MAP object containing fuel regions description is used the problem become even more complex because values have to be automatically changed for all bundles In order to clar
101. Montr al January 1984 W H Jens P A Lottes Analysis of heat transfer burnout pressure drop and density data for high pressure water ANL 4627 Argonne National Laboratory 1951 R W Bowring Physical model based on bubble detachment and calculation of steam voidage in the subcooled region of a heated channel Norway Institutt for Atomenergi OECD Halden Reaktor Prosjekt 1962 R T Lahey and F J Moody The Thermal Hydraulics of a Boiling Water Nuclear Reactor American Nuclear Society Publications U S A 1977 R Chambon Specifications and User Guide for NAP module in DRAGON DONJON VERSION5 Pin Power Reconstruction module Report IGE 345 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2014 M Fliscounakis E Girardi and T Courau A Generalized Pin Power Reconstruction Method for Arbitrary Heterogeneous Geometries M amp C 2011 Rio de Janero Brasil May 8 12 2011 IGE 344 186 38 G Rowlands Resonance Absorption and Non Uniform Temperature Distributions Journal of Nuclear Energy Parts A B 16 pp 235 236 1962 39 T Downar Y Xu and V Seker PARCS v3 0 U S NRC Core Neutronics Simulator User Manual Department of Nuclear Engineering and Radiological Science University of Michigan Ann Arbor MI 2012 40 A Ward Y Xu and T Downar GenPMAXS v6 1 1 Code for Generating the PARCS Cross Section Interface File PMAXS University of Michi
102. OLIB in which the properties are stored per each material region over the whole mesh splitted reactor geometry This MACROLIB is obtained by combining the material properties which are contained in the two distinct MACROLIB objects e The first MACROLIB contains the material properties which are evolution independent such as reflector and device properties It is created using either MAC CRE NCR or AFM module e The second is a fuel map MACROLIB created using either CRE NCR or AFM module It must contain the interpolated fuel properties per each fuel bundle The resulting MACROLIB will contain the properties that are stored for each reactor material and per each mesh splitted volume When the devices are not present in the reactor core then the resulting MACROLIB can be considered as a complete reactor MACROLIB and it can be directly used for the numer ical solving However when the devices are inserted into the reactor core the resulting MACROLIB is not yet complete it must be subsequently updated with respect to the device properties using the NEWMAC module see Section 3 10 The MACINI module specification is Table 7 Structure MACINI MACRO2 MATEX MACINI MATEX MACRO MACEL EDIT iprint where MACRO2 character 12 name of the extended MACROLIB to be created by the module MATEX character 12 name of the MATEX object containing an extended material index over the reactor geometry MATEX must be specified
103. OVE 3 level lt 1 2 L FADE x el n o ww Y b i Na gt level 1 level lt 1 part 2 FADE level lt 1 MOVE Figure 1 Presentation of fully and partially inserted 3 part control rods where EDIT iprint NUM ROD nrod FADE MOVE CREATE ROD GR ngrp dev rod rod group keyword used to set iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value larger values produce increasing amounts of output keyword used to specify nrod integer total number of the reactor rod type devices This number must be greater than 0 fading rod keyword A fraction of the fully inserted rod vanishes default option moving rod keyword The complete rod is moving DONJON3 type movement keyword used to create the rod groups of devices The creation of groups is optional keyword used to set ngrp integer total number of the rod groups to be created This number must be greater than 0 structure describing the input data for each individual rod structure describing the input data for each group of rods 3 4 2 Description of dev rod input structure A rod position is referred by its 3 D Cartesian coordinates only Note that the devices positions can not overlap The input order of data must be respected IGE 344 20 Table 10 Structure dev rod ROD id ROD NAME NAME AXIS X Y Z FROM H H LEVEL value
104. PE keyword used to specify the type of interpolation with respect to burnup data This information will be used during the execution of CRE NCR or AFM module TIMAV BURN keyword used to indicate the burnups interpolation according to the time average model This option is not available in 3 D Hexagonal geometry INST BURN keyword used to indicate the burnups interpolation according to the instantaneous model TIMAV BVAL keyword used to indicate the input of average exit burnup values per each combustion zone Note that the axial power shape and the first burnup limits will be reinitialized each time the average exit burnups are modified by the user These data are required for the time average calculation see Section 3 12 This option is not available with 3 D Hexagonal geometry INST BVAL keyword used to specify the instantaneous burnup values for each fuel bundle SMOOTH keyword used to level fuel mixtures burnup If the burnup is supposed to be the same at each occurence of every fuel mixture for symetry reasons SMOOTH will make sure they share the exact same value the first one in the burnup map Purpose is only to correct calculation noise in historic calculation ASBLY keyword to specify that one burnup value per assembly is to be defined OLDMAP keyword to specify that the burnup value is recovered from FLMAP2 The recovered burnup distribution is either from a previous calculation e with the same geometry but different initial
105. POS ENDMIX MIX lt lt mSORin gt gt SET FTYPE SOR SET RDTPOS 1 SET RDDPOS ENDMIX MIX lt lt mSORot gt gt SET FTYPE SOR SET RDTPOS 1 SET RDDPOS ENDMIX device specification x DEVICE MATEX Pdevc MATEX lt lt mZCRin gt gt lt lt mZCRot gt gt lt lt mSORin gt gt lt lt mSORot gt gt x full insertion of ZCR devices kh DEVICE x DSET DEVICE EDIT 1 ROD GROUP 1 LEVEL 1 0 END ROD GROUP 2 LEVEL 1 0 END ROD GROUP 3 LEVEL 1 0 END fuel map specification x FMAP MATEX k average dini Pfmap MATEX exit burnups FMAP Pburn FMAP 164 IGE 344 165 EVALUATE Eps epsil 1 WHILE Eps epsil gt iter 10 lt DO EVALUATE iter iter 1 SE fuel map macrolib H MACFL NCR LCPO FMAP EDIT O MACRO TABLE LCPO FUEL BURN MIX lt lt mFuel1 gt gt SET FTYPE CELL20 ENDMIX TABLE LCPO FUEL BURN MIX lt lt mFuel2 gt gt SET FTYPE CELL18 ENDMIX extended macrolib MACRO2 MATEX MACINI MATEX MACRO1 MACFL EDITO MACFL DELETE MACFL Pee complete macrolib MACRO MATEX NEWMAC MATEX MACRO2 DEVICE EDITO MACRO2 DELETE MACRO2 sa numerical solution Hz SYSTEM TRIVAA MACRO TRACK EDITO MACRO DELETE MACRO IF iter 1 THEN FLUX FLUD SYSTEM TRACK EDIT O ACCE 3 3 ADI 4 EXTE 1000 lt lt Precf gt gt THER 1000 ELSE FLUX FLUD FLUX SYSTEM TRACK ED
106. R in a Tape16 file Accordingly the keyword CELLAV should be used in the WIMS AECL run creating this file In addition if the REGION option is used in the T16CPO input data structure then it should also be used in the WIMS AECL run creating this file IGE 344 97 4 6 The D2P module The objective of the D2P module is to produce a file containing the macroscopic cross sections generated by the DRAGONS lattice code and readable by the GenPMAXS software This module makes possible the use of DRAGON integrated XS into the PARCS core code PARCS Purdue Advanced Reactor Core Simulator from the U S NRC 9 is a full 3D core code for the simulation of nuclear reactor steady state and transient behavior at a specific burnup state The main objective of the D2P module is to produce an output file with which can be accepted by GenPMAXS to produce the PMAXS file In order to minimize the development in the GenPMAXS code the choice has been made to reproduce an existing format already accepted by GenPMAXS The HELIOS output format has been selected A Microlib is extracted from a Saphyb or Multicompo obtained by DRAGON or APOLLO in a previous calculation The D2P module extracts cross sections contained in this microlib and creates two files e an input file needed by GenPMAXS to produce a PMAXS extention inp e a file containing data cross sections in HELIOS like format extention dra Note that the D2P module is compati
107. ROD 32 ROD NAME SORO2 LEVEL 0 0 AXIS Y FROM H MAXPOS 183 0 205 0 24 5 495 5 49 53 99 06 DMIX lt lt mSORin gt gt lt lt mSDRout gt gt ENDROD ROD 33 ROD NAME SORO3 LEVEL 0 0 AXIS Y FROM H MAXPOS 227 0 249 0 24 5 495 5 49 53 99 06 DMIX lt lt mSORin gt gt lt lt mSDRout gt gt ENDROD ROD 34 ROD NAME SORO4 LEVEL 0 0 AXIS Y FROM H MAXPOS 271 0 293 0 24 5 495 5 49 53 99 06 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 35 ROD NAME SORO5 LEVEL 0 0 AXIS Y FROM H MAXPOS 315 0 337 0 24 5 495 5 49 53 99 06 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 36 ROD NAME SORO6 LEVEL 0 0 AXIS Y FROM H MAXPOS 381 0 403 0 53 75 466 25 49 53 99 06 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 37 ROD NAME SORO7 LEVEL 0 0 AXIS Y FROM H MAXPOS 117 0 139 0 53 75 466 25 123 825 173 355 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 38 ROD NAME SOR08 LEVEL 0 0 AXIS Y FROM H MAXPOS 381 0 403 0 53 75 466 25 123 825 173 355 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 39 ROD NAME SORO9 LEVEL 0 0 AXIS Y FROM H MAXPOS 117 0 139 0 53 75 466 25 222 885 272 415 175 IGE 344 176 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 40 ROD NAME SOR10 LEVEL 0 0 AXIS Y FROM H MAXPOS 227 0 249 0 24 5 495 5 198 12 247 65 DMIX lt lt mSORin gt gt lt lt mSORout gt gt ENDROD ROD 41 ROD NAME SOR11 LEVEL 0 0 AXIS Y FROM H
108. RY is in modification mode the fuel weight is computed using the bundle lenght and the initial fuel density Now let us assume that a DRAGON calculation was performed for the cell located in bundle j 1 of channel i 1 We will also assume that these cells contain a single type of fuel Here the moderator temperature TMod is a global parameter while the fuel TComb and coolant TCalo temperatures are considered local parameters We assume that after the cell flux calculation a BURNUP data structure was generated using the following instructions k Procedures for cell calculation CellCalc k PROCEDURE CellCalc k Global parameter Tmod for moderator temperature Local parameters TComb for fuel temperature TCalo for coolant temperature k REAL TMod 345 66 REAL TComb TCalo 941 29 560 66 k Initial burnup options for cell calculation k REAL Power DeltaT 31 9713 5 0 k Local data structures k IGE 344 61 LINKED LIST Burnup Edition poe Execution control parameters icha channel number 1 ibun bundle number 1 us S INTEGER icha ibun 1 1 pann Perform cell calculation EL Burnup Edition CellCalc Burnup lt lt TComb gt gt lt lt TCalo gt gt lt lt TMod gt gt lt lt Power gt gt lt lt DeltaT gt gt Then assuming that the history structure HistXSM was created using the options above we can use go LOS Update history s
109. SPEED speed TIME time MAXPOS pos i i 1 6 DMIX mix1 mix2 ENDROD where ROD keyword used to specify the rod id number id integer identification number of the current rod Each rod type device must be assigned a unique id number given in an ascending order ranging from 1 to nrod ROD NAME keyword used to specify the rod NAME NAME character 12 name of the current rod In general this name is composed by the rod specific type e g SOR ZCR etc followed by its sequential number e g 01 02 etc AXIS keyword used to specify the rod movement axis A rod can be displaced along only one of the axis X keyword used to specify that a rod is displaced along X axis Y keyword used to specify that a rod is displaced along Y axis Z keyword used to specify that a rod is displaced along Z axis FROM keyword used to specify the insertion side of geometry The rod devices can be inserted into the reactor core from only one side of geometry For example some vertically moving devices can be inserted only from the top whereas other only from the bottom H keyword used to specify that a rod will be inserted into reactor core from the highest position e g from the top for vertically moving rod device H keyword used to specify that a rod will be inserted into reactor core from the lowest position e g from the bottom for vertically moving rod device LEVEL keyword used to specify the actual rod insertion level value
110. STORY data structure the values of these parameters The general form of this structure is IGE 344 59 Table 46 Structure hstpar NAMPAR valpar where NAMPAR name of a local or global parameter to process The parameters specified before the keyword CELLID is read will be considered global otherwise they will be considered local valpar real value for the local or global parameter to process In the case where the GET option is activated the history module will extract this parameter from the input data stream In the case where the PUT option is activated the history module will try to transfer this information into a real CLE 2000 variable 3 18 1 Example The history interface between the codes DRAGON and DONJON has been written as a new module in order to facilitate the access to the GANLIB utilities that manage the required hierarchical data structures The resulting HST module can be called both by DRAGON and DONJON The reactor model we will consider as an example is a 3 D model with an x 3 y 3 and z 3 mesh Here we will assume that the x y plane describes fuel channels The z plane will be associated with the so called fuel bundles This choice is somewhat arbitrary however it is useful if the refueling takes place in a specific direction as in a CANDU reactor Here a 2 bundle shift fueling strategy will be considered To each fresh fuel cell introduced in the core the HST module will associate a unique c
111. TABLE 82 83 STEP 121 STEP REDUCT 115 116 STOP 155 STRNAME 6 7 195 SYSTEM 5 40 T 105 T BURNUP 65 66 T H invariant variables mandatory 99 T16CPO 93 94 96 TABLE 64 65 67 69 78 80 Tape16 93 96 TAVG 4 39 40 TCalo 60 TCOM 102 103 153 TComb 60 TCOOL 88 90 tcool 88 90 TEFF 108 110 TEXT80 78 79 TFUEL 88 90 tfuel 88 90 TH_DATA 154 TH_DATA 154 THERMO 107 THERMO 107 thm 143 144 THM 107 108 144 147 TIMAV BURN 12 65 69 80 81 TIMAV BVAL 12 TIME 20 21 26 27 29 30 43 46 47 51 52 107 108 time 20 21 26 27 29 30 107 108 timeiter 107 108 TIMES 12 14 times 10 timestep 107 108 tinlet 108 109 TINST 4 42 43 tinv 22 23 TMOD 88 90 102 103 153 TMod 60 tmod 88 90 TO 43 44 TRACK 36 37 51 TRACK 5 37 51 TRACKING 120 TRIVAA 3 5 40 TRIVAT 3 5 37 TRKNAM 120 TTD 31 32 TYPE 22 23 112 113 type 115 116 UNI 31 32 unit 108 110 unrodded 101 102 UPS 65 IGE 344 USE 69 80 USER 101 103 USPLIT 4 5 14 15 51 131 UTL 3 val 103 vall 69 70 80 81 val2 69 70 80 81 valpar 59 valper 94 96 valperd 94 96 valpert 94 96 valref 69 70 80 81 94 95 valrefd 94 96 valreft 94 96 VALUE 112 113 value 20 26 27 29 30 33 34 valuel 31 32 value2 31 32 VAR VAL MAX 115 116 VAR VAL MIN 115 116 VAR VALUE 115 116
112. TECHNICAL REPORT IGE 344 A USER GUIDE FOR DONJON VERSION5 A HEBERT D SEKKI AND R CHAMBON in Institut de g nie nucl aire 64 D partement de g nie m canique cole Polytechnique de Montr al 64 bit September 24 2015 IGE 344 i Contents CON nas PL AAA ee DA es ES aaa i E ol oe a a as E TTT TTT v List Of FIQUES 6 4 5 4 26444 Db RAS a ee ee e GER RA viii 1 INTRODUCCI N lodo aid he LL LEE 1 2 GENERAL SPECIFICATION OF DONJON oo o e 3 2 1 A O E Rea ek See wed 3 2 2 Data structures corra e ionu e a Re he e ee 5 2 3 Syntactic rules for input specification 2 2 2 0 2 0 0000 eee 6 2 4 General input StructUre e co k A ee ee eae ee 6 3 GENERAL CORE DESCRIPTION MODULES o e e 8 3 1 The RESTN module k m e Pita tae ws 8 DLE Main input data to the RESINI module 9 3 1 2 Input of global and local parameters o 11 3 2 The USPLIT modules lt lt csi pe LA AG a a a eee 15 2 2 Input data to the USPLIT module 15 3 3 The MAGINI module so o uce wa Le E 17 3 4 The DEVINE modula sa eA RAR RR Re ee RD AR we 18 3 4 1 Input data to the DEVINI module 18 3 4 2 Description of dev rod input structure ks 19 34 3 Description of rod group input structure 21 3 5 The DETINI module i ra nce sel a eee thane dba ea a ds 22 dde l Input data to the DETINI module c ev eee ees 22 52 Description of
113. Transport cross section Absorption cross section Nu fission cross section Kappa fission cross section Microscopic capture cross section of Xenon Microscopic capture cross section of Samarium Fission cross section Detector response parameter cattering cross section Assembly discontinuity factor Direct energy deposition T Factors Corner discontinuity factor Group Wise form function Y 4 6 2 General format of the module The D2P module can perform a sequence of phases related to the generation of a cross section format readable by GenPMAXS e PHASE 1 recover input data from Saphyb and create output files 1 recover information from a Saphyb file 2 store general information in output file 3 generate the GenPMAXS input file e PHASE 2 recover crosss section from microlib thanks to the SCR module and store in memory e PHASE 3 store cross sections in output file The general format of the data for the D2P module is the following Table 68 Structure D2P HEL GEN INF D2P INF SAP MCO PHASE 1 EDIT iprint descphase1 GEN INF D2P MIC INF GEN SAP MCO PHASE 2 EDIT iprint HEL GEN INF D2P INF GEN HEL PHASE 3 EDIT iprint In the DRAGON formalism the Left Hand Side LHS is dedicated to the objects created or modified by the module the Right Hand Side RHS is used for input objects all parameters are passed to the module after the delimiter where HEL ascii fil
114. Units Comment Signature of the power data structure SIGNA L_POWER uu Vector describing the various parameters as sociated with this data structure S The total reactor power The total reactor volume The flux normalization factor Fuel mixture indices per fuel bundle cm s The normalized fluxes over the whole reactor geometry recorded per each mesh splitted vol ume and per each energy group The flux val ues over the virtual regions are set to 0 cm The volume of each fuel bundle em s7 The normalized average fluxes recorded per each fuel bundle and per each energy group cm s_ The normalized flux distribution over the whole reactor geometry recorded per each X Y Z planes and per each energy group The fluxes ratios with respect to the thermal energy group fluxes The bundle powers The channel powers The power distribution over the reactor core recorded per each X Y Z planes The power values over the non fuel regions are set to 0 The maximum channel power The maximum bundle power The radial power form factor defined as maximum to average channel power in core continued on next page IGE 344 140 Records in power data structure continued from last page Condition Units Comment FORM BUNDouu R 1 The overall power form factor defined as maximum to average bundle power in core K EFFECTIVE R 1 The effective multiplication factor recovered from the flux data structure All stored f
115. VE Ve for TA where f represents the average value between two points The third case corresponds to a minimal grid where the lattice computations have been perfomed only for one parameter variation at a time In this case the grid is represented by the thick gray lines on the axis on Fig 6 and 7 If we use the notations of Fig 6 and 7 the best estimate interpolated values f we can get are given by I FV gt FVo Va f a Ve FO HF Ve Fe fF Va Va F Vo f Vc f Vo for instataneous f FVV FOB Ve S Va f Vo S Vo f Vo for TA Note that the reference point Vo in the example does not have to be the same for all parameters Database structures such as represented on Fig 8 can also been used In this case we even have two choices for the Af computation on axis A The last case is in fact a mix of cases 2 and 3 The gray rectangle and the gray line on Fig 9 and 10 reprensent where all the lattice computations have been performed With the notations used on those IGE 344 72 figures one can write that the best estimate interpolated values f we can get are given by f fV f Ve f Vec f Ve f Va f W f Vac f Va f Vo for instataneous KVV f Vg Va Vbo Vac FV Va FVA V0 Vbo Vee f Va f Vo for TA Note once again that the reference point Vo in the example does not have to be the same for all parameters D
116. ace for its production yield branching ratio or production yield expressed in fraction NAMPAR character 12 name of the a parent isotope or isomer that appears as a particularized isotope of the Saphyb ENDCHAIN keyword to specify the end of the depletion chain 4 3 4 Interpolation in the parameter grid The following example corresponds to a delta sigma computation in mixture 1 corresponding to a perturbation Note that in this case the MACROLIB object may content negative cross section MACROLIB SCR SAP EDIT 40 NMIX 1 SAPHYB SAP MIX 1 delta sigma contribution SET CELL 3D DELTA PITCH 0 0 1 0 ENDMIX When the number of parameters used for the interpolation is increased all the lattice computations corresponding to all the combinations of parameters may not be done for computation time reasons In this case some approximations may be required The choice for the SET DELTA and ADD is then dependent of the structure of the database i e how the database grid of possibilities is filled When a MAP object containing fuel regions description is used the problem become even more complex because values have to be automatically changed for all bundles In order to clarify all the different possibilities and limitations dependently of the database structure we will use a 3 parameter case The paramaters are referenced by A B and C But before we explain the different cases we want to remind tha
117. ard full voiding pattern checkerboard half voiding pattern checkerboard quarter voiding pattern user defined voiding pattern Ivoida NOOR WNrF CO IGE 344 e The type of the geometry F Sia f 7 Cartesian 3 D geometry 9 Hexagonal 3 D geometry 126 The naval coordinate layout used by the SIM module Lim S13 The number of assemblies along X and Y axis are given using Isim Ly 100 and Ly mod Jsim 100 The total number of assemblies Nass S14 The number of assemblies along X direction Nya S15 e The number of assemblies along X direction Nya S16 8 1 2 The main fmap directory The number of plane of the mesh along Z direction where assemblies are situated Nz ass S17 The following records and sub directories will be found on the first level of fmap directory SIGNATURE uu STATE VECTOR FLMIX uuuuvy FLMIX INI uuu S ZONE LULU BMIX uuuuuuva XNAME y uuuuua YNAME uuuuua Table 82 Records and sub directories in fmap data structure Condition Units Comment Signature of the fmap data structure SIGNA L_MAP uuu Vector describing the various parameters as sociated with this data structure S Fuel mixture indices per bundle or assembly subdivisions for each reactor channel Fuel mixture indices per bundle or assembly subdivisions for each reactor channel as de fined by user in RESINI module identification name corresponding to the basic naval coordinate p
118. arger values produce increasing amounts of output keyword used to set the power used for cross section interpolation a distributed beginning of transient bundle power in kW is used This power distribu tion has to be pre calculated in the FLPOW module using the INIT keyword uniform bundle power in kW If this data is omitted the reference value in the data base is used or the bundle powers present in a MAP The reference value is 615 kW if none were provided at the database computation time keyword used to activate power bundle feedback on fuel properties using powers recov ered from BUND PW record in MAPFL This is the default option if MAP is selected keyword used to desactivate PWF feedback This is the only possible option if MCR is selected IGE 344 TFUEL tfuel TCOOL tcool TMOD tmod BORON nB RDCL dcool RDMD dmod PUR purity BURN bval XENON nXe XEREF NEP nNp 90 keyword used to set tfuel fuel temperature in K If this data is omitted and the bundle powers present in a MAP fuel temperatures are computed with respect to powers If this data is omitted and there is no bundle power the reference value in the data base is used where it is 941 29 K if none were provided at the database computation time keyword used to set tcool coolant temperature in K If this data is omitted the reference value in the data base is used The reference value is 560 66 K if none
119. associated with cell C as well as the current isotopic content of this cell The identification number C associated with channel j and bundle i can be seen as the serial number of the bundle located at a position in space identified by i j It is automatically managed by the HST module l For a fresh core Ci j n where n represents the cell order definition in the input file Upon refueling some bundles in channel k of the core are displaced from region I k to m k new bundles are introduced at location I k and old bundles removed from location m k If one assumes that CNEW and COLD represents the value of C after and before refueling then we will have NEW _ OLD La Ci k NEW _ FRESH Chic T Cri k IGE 344 where CELERE 142 represent a fresh fuel cell The local parameters and burnup power density of the fuel cell previously located at m k are preserved and the fresh fuel isotopic densities is that provided in Fm k Hu FRESH the fuel type associated with CY 8 6 2 The fuel type sub directory Each fuel sub directory FUEL contains the following information Table 98 Fuel type sub directory Type Condition Units FUELDEN INIT ISOTOPESUSED ISOTOPESMIX ISOTOPESDENS 8 6 3 The cell type sub directory Comment array containing the initial density of heavy element in the fuel pp in g cm and the initial linear density of heavy element in the fuel my in g cm array containi
120. atabase structures such as represented on Fig 11 can also been used IGE 344 73 The input files will actually reflect the previous equations However they are different if the param eters are stored in a MAP object MAPFL or provided directly by the user For the case of one point interpolation i e instantaneous the input files will be all parameters explicitly set all parameters in MAP MACROLIB NCR CPO NMIX 1 MACRO COMPO CPO default MIX 1 SET A lt lt va gt gt SET B lt lt vb gt gt SET C lt lt vc gt gt ENDMIX MACROLIB NCR CPO NMIX 1 MACRO COMPO CPO default MIX 1 SET A lt lt va gt gt SET B lt lt vb gt gt SET C lt lt vc0 gt gt ADD C lt lt vc0 gt gt lt lt vc gt gt REF PA lt lt va0 gt gt B lt lt vb gt gt ENDREF lor B SAMEASREF ENDREF ENDMIX MACROLIB NCR CPO NMIX 1 MACRO COMPO CPO default MIX 1 SET A lt lt va0 gt gt SET SET ADD 7B IC 2A REF ADD 22 REF ENDMIX lt lt vb gt gt lt lt vc0 gt gt lt lt va0 gt gt lt lt va gt gt 2C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF lt lt vc0 gt gt lt lt vc gt gt 7A lt lt va0 gt gt B lt lt vb0 gt gt ENDREF MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 ENDMIX MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET C l
121. atory In this block the geometrical configuration of core reactor is specified and also the number of ADF in each group the number of rod rows and columns in whole assembly the rod lattice pitch etc Some of these parameters have default values XS_SET 00000001 1 1 1 17 17 3 1 44270 0 72135 0 72135 2 78613 0 73659 0 00016 0 00000 0 00000 Block HISTORY CASE identification mandatory This block describes the state variables values for all histories contained in the PMAXS file The first parameter refers to the burnup set HISTORYC 1 0 00000 0 71100 1000 00000 900 00000 Block T H invariant variables mandatory It contains invariant variables if the corresponding logical flag in block XS CONTROL is set to T The list of invariant variables repeated for each burnup points is H NO e Chi spectra Yield of I Xe and Sm Beta of delayed neutron Lambda of delayed neutron Decay heat data 00000E 00 0 00000E 00 5 13949E 08 2 17697E 06 87406E 14 1 13375E 13 4 25558E 15 68628E 04 1 43419E 03 1 39641E 03 3 23740E 03 1 43931E 03 5 99082E 04 33535E 02 3 26045E 02 1 21056E 01 3 05531E 01 8 60559E 01 2 89034E 00 This block contains the necessary information for 1 burnup point and for each energy group Chi inV YLD Bet Lam Block XS data mandatory Cross sections in PMAXS file are listed for each burnup point and for each neutron energy group Some cross sections are optional see table below IGE 344 100 XS data block
122. ble with the last version of GenPMAXS v6 1 3 and PARCS v32m13 etc 4 6 1 The PMAXS format The depletion capabilty of the PARCS code is reachable thanks to a specific format of cross section file named PMAXS Purdue Macrosocopic Cross Section File This format is generated using the GenPMAXS code based on output files of several lattice codes such as HELIOS CASMO TRITON WIMS CONDOR and SERPENT This module intend to add DRAGON in this list The macroscopic cross sections are stored in the PMAXS file using partial derivatives as a function of state variables Consequently the PMAXS format is a multi dimentional table including burnup dependance This format is a flexible way to obtain a more or less accurate meshing of cross sections In addition of burnup the list of state variables around which PMAXS is built is the following 1 CR control rod fraction 2 DC density of coolant PC soluble poison concentration in coolant TF temperature of fuel TC temperature of coolant IC impurity of coolant DM density of moderator PM soluble poison concentration in moderator A oS eb A TM temperature of moderator 10 IM impurity of moderator 11 DN density difference between neighbor and current assembly 12 BN burnup difference between neighbor and current assembly IGE 344 98 These variables should be specified in this order With the exception of burnup each variable is optional The following equation is
123. bly subdivisions for each channel MW d t7 Low burnup integration limits according to the time average model MW d t Upper burnup integration limits according to kW kW 2 s71 the time average model Bundle powers set in RESINI module or re covered from L_POWER object Beginning of transient bundle powers recov ered from L_POWER object The normalized average fluxes recorded per each fuel bundle and for each energy group recovered from L_POWER object MW d t Core average discharge burnup kg dot d d kW Bundle shifts per combustion zone A bundle shift corresponds to the number of displaced fuel bundles during the refueling operation Refueling pattern vector per combustion zone Refueling scheme of each channel it corre sponds to the positive or negative bundle shift number according to the flow direction Channel refueling rates Time values at which channels are refueled in side a refueling time period Refueling time period in days The power of the bundles shifted the i th time MWdaT The burnup of the bundles shifted the i th time The number of shifts per bundle during refu eling Normalized axial power shape values over the fuel bundles Equal to fuel bundle powers di vided by channel powers Convergence factor for the axial power shape calculation it is defined as a relative error be tween the two successives calculations Sub directory containing the embedded 3D Cart
124. burnup integration limits are calculated approximately using the RESINI module 2 A time average integration is performed and a new fuel map MACROLIB is created using either NCR CRE or AFM module 3 An extended MACROLIB over the whole reactor geometry is created using the MACINI module 4 If the devices are inserted into the reactor core then the previously created MACROLIB is to be updated for the devices properties using the NEWMAC module 5 The complete MACROLIB is subsequently used by the TRIVAA module in order to create a matrix SYSTEM 6 The full core numerical solution i e fluxes and effective multiplication factor is computed using the FLUD module 7 The channel and bundle powers are next calculated using the FLPOW module 8 Finally the new axial power shape and burnup limits are computed using the TAVG module IGE 344 41 Note that the steps from 2 to 8 are to be repeated until the reguired precision for the axial power shape convergence is satisfied IGE 344 42 3 13 The TINST module The TINST module is used to compute the instantaneous burnup for each fuel bundle You can also use TINST to refuel your reactor according to a refueling scheme The scheme can be either specified with RESINI or directly in TINST The TINST module specification is Table 30 Structure TINST FMAP TINST FMAP POWER MICLIB3 FMAP TINST FMAP MICLIB2 MICLIB y desctinst where
125. burnup options for cell calculation A is after refueling B is before refueling k IGE 344 REAL TCombA TCaloA TCombB TCaloB REAL PowerA DeltaTA PowerB DeltaTB Burnup HST History PUT TMod gt gt TMod lt lt CELLID lt lt icha gt gt lt lt ibun gt gt PUT BREFL BURN gt gt DeltaTB lt lt gt gt PowerB lt lt TComb gt gt TCombB lt lt TCalo gt gt TCaloB lt lt AREFL BURN gt gt DeltaTA lt lt gt gt PowerA lt lt TComb gt gt TCombA lt lt TCalo gt gt TCaloA lt lt IF DeltaTB 0 0 gt THEN goes Burn before refueling eae Burnup Edition CellCalc Burnup lt lt TCombB gt gt lt lt TCaloB gt gt lt lt TMod gt gt lt lt PowerB gt gt lt lt DeltaTB gt gt Edition DELETE Edition ENDIF pets Burn after refueling pecca Burnup Edition CellCalc Burnup lt lt TCombA gt gt lt lt TCaloA gt gt lt lt TMod gt gt lt lt PowerA gt gt lt lt DeltaTA gt gt a Update History goose History HST History Burnup CELLID lt lt icha gt gt lt lt ibun gt gt 62 Note that here there are two sets of local parameters that can be extracted from the history data structure namely the before BREFL and the after AREFL refueling information In the case of fresh fuel single fuel description or a refueled bundle extracting the before information is not required However if one uses the general procedure described above to extrac
126. ch calcualtion Each component of the list CROSS_SECT is a directory containing cross sections for all burnup point of a specific branch Cross sections for a single burnup points are filled segentially during the branching calculation Inside each groupwise directory the following sub directories will be found Table 112 Sub directories in CROSS_SECT Type Condition Units Comment MACROLIBXS Dir Directory for macroscopic cross sections for a given calculated branch cross sections are overwrited after each branch calculation continued on next page IGE 344 Sub directories in CROSS_SECT MICROLIB_XS XTRo oo ABSORPTION KAPPA FI uu X_NU_FI LULU SCATLULLLULU DET uuuuuuuva SFI ooo uu ADF uvuuuuua Name XENGuuuuuuua SMNG Guo Type Dir Condition lxes Units 156 continued from last page Comment Directory for microscopic cross sections for a given calculated branch cross sections are overwrited after each branch calculation Table 113 Records in the sub directory MACROLIB_XS Condition Units Comment The transport cross section 2 Lf 1 3 Dy The absorption cross section VI DY Ef Ny BP The product of DA the fission cross section with the energy emitted by this reaction K KOS The product of DI the fission cross section with the averaged number of fissionemitted delayed neutron v VMS Scattering cross section elements fro
127. ch channel is represented using 1D convection equations along the channel and 1D cylindrical equations for a single pin cell A two fluid homogeneous model is used The THM module is built around freesteam an open source implementation of IAPWS IF97 steam tables for light water The THM module works both in steady state and in transient conditions and includes a subcooled flow boiling model based on the Jens amp Lottes correlation 33 and on Bowring s model for two phase homogeneous flows P4 The 1D thermal hydraulics equations are solved in each channel as a fonction of two fixed inlet conditions for the coolant velocity and temperature and one fixed outlet condition for the pressure The THM module specification is Table 70 Structure THM THERMO MAPFL THM THERMO MAPFL descthm where THERMO character 12 name of the THERMO object that will be created or updated by the THM module Object THERMO contains thermal hydraulics information set or computed by THM in transient or in permanent conditions such as the distribution of the enthalpy the pressure the velocity the density and the temperatures of the coolant for all the channels in the geometry It also contains all the values of the fuel temperatures in transient or in permanent conditions according to the discretisation chosen for the fuel rods MAPFL character 12 name of the MAP object containing fuel regions description and local parameter informations
128. ch the properties are stored per each material region over the whole mesh splitted reactor geometry module used for 3 D modeling of rod type devices in the reactor core module used to read and store detector information module used for 3 D modeling of liquid zone controllers in the reactor core module used to set the new devices parameters that can be used for the reactivity worth studies module used to modify the fuel map in order to compute a Doppler or general reactivity coefficient module used to compute the time dependent positions of the moving rod type devices module used to create an extended MACROLIB that will contain the updated material properties computed with respect to the actual devices positions module used to compute and print powers and normalized fluxes over the reactor core module used to perform burnups calculation according to the time average model compute burnups integration limits core average exit burnup axial power shapes and channel refuelling rates module used to perform burnups calculation according to the time linear model and compute instantaneous burnups values This module is specific to Candu reactor refuelling module used to perform burnups calculation according to the time linear model and compute instantaneous burnups values This module is specific to PWR reactor refu elling module used to compute the mean flux at each detector site and the response of each detector according
129. clad thermal conductivity as a function of local clad temper ature Teraa Clad conductivity is computed with the following polynomial ncond Aaa Y kcond k Tetaa k 0 with Acad in W m K and Teraa in the selected unit Kelvin or Celsius By default a built in model is used keyword used to set the heat exchange coefficient of the gap as a constant By default a built in model is used real value set to the constant heat exchange coefficient of the gap in W m K keyword used to set the heat transfer coefficient between clad and fluid as a constant By default this coefficient is computed using a built in correlation real value set to the constant heat transfer coefficient between clad and fluid in W m K keyword used to set the weighting factor in the effective fuel temperature approxima tion The effective fuel temperature is used for the cross sections interpolations on fuel temperature real value Wieg set to the weighting factor in the effective fuel temperature The effective fuel temperature is computed as fuel fuel fuel Tot Weert tada 1 Wier Toen center where 0 lt Wieg lt 1 T is the temperature at the surface of the fuel pellet K surface and T uel is the temperature at the center of the fuel pellet K 38 By default the Rowlands weighting factor Wieg 2 is used IGE 344 111 CONV keyword used to set the convergence criteria for solving the conduction and the con serva
130. containing fuel regions description global and local parameter information burnup fuel coolant temperatures coolant density etc Keyword TABLE is expected in ncr_data ncr data input data structure containing interpolation information see Section 4 2 1 4 2 1 Interpolation data input for module NCR Table 53 Structure ncr_data EDIT iprint ALLX nbfuel RES MACRO MICRO LINEAR CUBIC LEAK b2 NMIX nmixt COMPO CPONAM NAMDIR descintf TABLE CPONAM NAMDIR namburn descintf gt where EDIT keyword used to set iprint IGE 344 68 iprint index used to control the printing in module NCR 0 for no print 1 for minimum printing default value ALLX keyword used to register the region number of each isotope before merging This option is useful if the same keyword has been specified in EDI and COMPO before nbfuel number of fuel rings used for micro depletion calculations RES keyword indicating that the interpolation is done only for the microscopic cross sections and not for the isotopic densities In this case a RHS MICROLIB must be defined and the number densities are recovered from it This option is useful for micro depletion applications Important note It is possible to force interpolation of some isotopic densities with RES option if these isotopes are explicitely specified with a flag after MICRO keyword in descintf input data structure
131. cretization of the reactor geometry and to provide the numerical solution according to the user selected numerical procedure 7 The UTILIB library provides the utility and linear algebra libraries Finally the GANLIB computer code provides CLE 2000 capabilities to control data flows and to implement computational schemes GANLIB also provide LCM data structures to exchange information between modules The DONJON code is divided into several modules each module is designed to perform some partic ular tasks The transfer of information between the modules is achieved by means of well defined data structure Several design features data structure and computing algorithms were recovered revised and adapted from the previous DONJON versionl One of the main concerns of the DONJON developers is to ensure the code reliability and extensibility The DONJON modules are first designed for the reactor full core modeling in 3 D Cartesian geometry These modules are built around the reactor fuel lattice specification corresponding to the common design features of CANDU reactors The modules related to the modeling of reactivity mechanisms which are normally presented in the reactor core also constitute an important part of code The DONJON code can perform several full core calculations and can be used to determine some important core characteristics such as the power and normalized flux distributions over the reactor core All full core calculations us
132. ction see Section 7 1 2 IGE 344 121 structure containing the input data to this module to automatically define the core geometry with heterogeneous assembly see Section 7 1 3 descnap3 7 1 1 Additional properties calculations Table 79 Structure descnap1 EDIT iprint PROJECTION STEP namedir IFX ifx SET pname pvalue 2 where EDIT keyword used to modify the print level iprint iprint integer index used to control the printing in module NAP 0 for no print 1 for min imum printing default value larger values of iprint will produce increasing amounts of output PROJECTION keyword to specify that additional properties for subregions will be computed and stored in the MULTICOMPO data structure COMPO STEP keyword to specify namedir namedir name of the directory containing the homogenized cross sections homogeneous or heterogeneous IFX keyword to specify ifx ifx number used to create the name of the flux record representing we This flux rep resents the results of calculation in an infinite domain computed in diffusion with cross sections homogenized either homogeneously or heterogeneously One record is associated with each type of homogenization when the MULTICOMPO is created Thus MULTICOMPO can be enriched several times each time using a different homogeniza tion of the assembly at the end of the transport calculations The following format is used for the flux record name
133. ctors can be recall in case of error of the mathematical programming CALCUL DX keyword used to specify that the new step will be calculated NO STORE OLD keyword used to specify that the old value of decision variables and gradients will not be stored in the L OPTIMIZE 0LD VALUE directory COST EXTRAP keyword used to calculate the extrapolated objective constant ecost IGE 344 119 ecost extrapolated objective constant CONV TEST keyword used to calculate if the external convergence has been reached u 1 means that external convergence has been reached 0 otherwise IGE 344 120 7 PIN POWER RECONSTRUCTION MODULES This section is related to pin power reconstruction capabilities available in Donjon The corresponding theory is explained in 36 37 7 1 The NAP module The NAP module supplies the main transport diffusion equivalence options to DRAGON and DON JON It can be used to perform the pin power reconstruction 97 The calling specifications are Table 78 Structure NAP COMPO NAP COMPO TRKNAM FLUNAM descnap1 MAP NAP MAP TRKNAM FLUNAM MATEX MACRES descnap2 GEONEW NAP GEOOLD COMPO descnap3 where COMPO TRKNAM FLUNAM MAP MATEX MACRES GEONEW GEOOLD descnap1 descnap2 character 12 name of the MULTICOMPO data structure L_COMPO signature where the detailed subregion properties will be stored character 12 name of the read only TRACKING data
134. d hexagons in the first X Y plane POSITION uu R 6 The coordinates of the detector RESPON uu R Z2 The responses of the detector 8 5 Contents of power data structure A power data structure is used to store the information related to the powers and fluxes over the reactor core This object has a signature L_POWER it is created using the FLPOW module The reactor fluxes and powers are recorded using several data formats 8 5 1 The state vector content The dimensioning parameters S which are stored in the state vector for this data structure represent e The number of energy groups Ng S IGE 344 139 e The total number of mesh splitted volumes Nr S2 e The number of mesh splitted volumes along x axis Ly S3 e The number of mesh splitted volumes along y axis Ly S4 The number of reactor channels Nen S e The number of bundles per channel NV S7 8 5 2 The power directory The number of mesh splitted volumes along z axis L S5 The following records will be found on the power directory SIGNATURE uu STATE VECTOR PTOT Guu UU VTOT Goo UK NORM ivuuuuua FLMIX uuuuva FLUX uuuuuuva VOLU BUND uu FLUX BUND uu FLUX DISTRuu FLUX RATIO u POWER BUND u POWER CHAN u POWER DISTR PMAX CHAN uu PMAX BUND uu FORM CHAN uu Table 96 Records in power data structure Nen No Nen No Nor Talla Ng Ly Ly La Nor ES 1 Nen No Nen La Lilia Condition
135. d in the HISTORY data structure nbun the number of bundles per channels for the reactor model Note that if nbun is different from the value already defined on the HISTORY data structure or the MAP data structure the execution will be aborted bunl bundle length in cm This information is required to compute inital fuel weight CHANNELS keyword used to specify the number of fuel channels for the reactor model considered in the HISTORY data structure ncha the number of fuel channels for the reactor model Note that if ncha is different from the value already defined on the HISTORY data structure or the MAP data structure the execution will be aborted The hstbrn serves a unique purpose mainly to extract from the HISTORY file the information required to process a burnup evaluation in DRAGON using the EVO module The information must be stored inside CLE 2000 variables The general form of this output structure is Table 45 Structure hstbrn BURN period power where BURN keyword to indicate that burnup information follows period the burnup period in days that will be transferred to a real CLE 2000 variable power the power density in kW kg that will be transferred to a real CLE 2000 variable The hstpar serves two purposes First it is used to define the names of the local and global parameters that may be used in our calculations as well as the values of these local parameters In can also be used to extract from a HI
136. d used to indicate the initialization of the library for a recursive calculation using the XENON module The Xenon concentration is set to zero for all the bundles IGE 344 51 3 16 The DETECT module The DETECT module is used to compute the mean flux at each detector site and the response of each detector The DETECT module specifications are Table 36 Structure DETECT DETEC DETECT DETEC FLUX TRACK GEOM descdetect where DETEC character 12 name of the DETECT containing the detector positions and responses FLUX character 12 name of the FLUX containing the flux solution computed by the FLUD or FLPOW modules To obtain a correct result the best is to use a normalized flux coming from the FLPOW module In this case the fluxes are normalized to the reactor power TRACK character 12 name of the TRACK containing the TRIVAC tracking GEOM character 12 name of the GEOMETRY containing the mesh splitting geometry created by the USPLIT or GEO modules descdetect structure containing the data to module DETECT 3 16 1 Input data to the DETECT module Note that the fuel lattice power distribution can be printed only on the screen Table 37 Structure descdetect EDIT iprt TIME dt REF kc NORM vnorm SIMEX SPLINE PARAB where EDIT key word used to set iprt iprt index used to control the printing in module DETECT 0 for no print 1 for minimum printing default value 4 for pri
137. ded directly by the user For the case of one point interpolation i e instantaneous the input files will be all parameters explicitly set all parameters in MAP MACROLIB NMIX 1 SAPHYB MIX 1 SET SET SET ENDMIX MACROLIB NMIX 1 SAPHYB MIX 1 SET SET SET ADD REF lor ENDMIX MACROLTB NMIX 1 SAPHYB MIX 1 SET SET SET ADD REF ADD REF ENDMIX SCR SAP SAP PA lt lt va gt gt 7B lt lt vb gt gt 2C lt lt vc gt gt SCR SAP SAP 2A B 262 63 lt lt va gt gt lt lt vb gt gt lt lt vc0 gt gt lt lt vc0 gt gt lt lt vc gt gt 7A lt lt va0 gt gt B lt lt vb gt gt ENDREF B SAMEASREF ENDREF SCR SAP SAP 2A 7B 20 A lt lt va0 gt gt lt lt vb gt gt lt lt vc0 gt gt lt lt va0 gt gt lt lt va gt gt C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF lt lt vc0 gt gt lt lt vc gt gt 2A lt lt va0 gt gt B lt lt vb0 gt gt ENDREF 20 MACROLIB NMIX 1 TABLE SAP MIX 1 ENDMIX MACROLIB NMIX 1 TABLE SAP MIX 1 SET ADD 262 262 REF SET ADD A 2A REF ENDMIX MACROLIB NMIX 1 TABLE SAP MIX 1 SET SET ADD DA 262 A REF ADD C REF ENDMIX SCR SAP FMAP B SCR SAP FMAP B lt lt vc0 gt gt lt lt vc0 gt gt MAP gt A lt lt va0 gt gt 7B SAMEASREF ENDREF lt lt va0 gt gt lt lt va0 gt gt
138. dex numt and sorted channel by channel numt can take values between 1 and 9999 In each of the Nen CHANNEL numc sub directories the following records will be found Table 102 Records in each CHANNEL directory Condition Units Comment VINLET uuu TINLET uuu PINLET uuu VELOCITIES u inlet velocity inlet temperature inlet pressure velocity in each of the N bundles of the chan nel numbered numc pressure in each bundle of the channel enthalpy in each bundle of the channel density in each bundle of the channel density of liquid phase in each bundle of the channel distribution of the temperature in the fuel pin for each bundle of the channel center fuel pellet temperature in each bundle of the channel PRESSURE uuu ENTHALPY uuu DENSITY Luv LIQUID DENS TEMPERATURES R R R R R R R R R R CENTER TEMPS IGE 344 148 8 8 Contents of a optimize data structure The optimize specification is used to store the optimization variables and functions values and definitions limits and options In any case the signature variable for this data structure must be SIGNA L OPTIMIZE jj The dimensioning parameters for this data structure which are stored in the state vector S represents e The number of decision variables Nyar S e The number of constraints Nest e The type of optimization S where di 1 minimization 3 1 maximization e The result of a test for external convergenc
139. ds and sub directories in device data structure Condition Units Comment SIGNATURE uu Signature of the device data structure SIGNA L_DEVICE uuu STATE VECTOR Vector describing the various parameters as sociated with this data structure S DEV_ROD uuuu i Sub directories for each controller rod A sub directory is created for each controller rod ac cording to the rod identification number ROD_GROUP uu i Sub directories for each group of rod type de vices A sub directory is created for each group of rod type devices according to the rod group identification number continued on next page IGE 344 134 Records and sub directories in device data structure continued from last page Condition Units Comment DEV_LZC uuu Dir Nzc Sub directories for each liquid zone controller A sub directory is created for each liquid con troller according to the liquid controller iden tification number LZC_ GROUP uu Dir Nigrp Sub directories for each group of lzc type de vices A sub directory is created for each group of lzc type devices according to the lzc group identification number 8 3 3 The DEV ROD sub directories Inside each DEV ROD sub directory the following records will be found Table 89 Records in DEV ROD sub directories Condition Units Comment ROD IDuuuuuu The identification number of the rod ROD NAME uuu The identification name of the rod ROD PARTS vu The number of parts in the rod Npart gt 1
140. e Output file with HELIOS like format compulsory if iphase 1 in creation mode or if iphase 3 in modification mode GEN ascii file Input file for running GenPMAXS compulsory if iphase 1 in creation mode or if iphase 2 or 3 in modification mode IGE 344 101 INF LCM object Block of data for the dialogue between different seguence of operations SAP Saphyb object with cross sections to be extracted compulsory if iphase 1 or 2 MCO Multicompo object with cross sections to be extracted compulsory if iphase 1 or 2 MIC microlib object with cross sections for one burnup point compulsory if iphase 2 PHASE keyword used to set iphase iph integer index used to control the current phase of D2P module EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print default value 1 for minimum printing for larger values of iprint everything will be printed descphase1 input data structure for PHASE 1 of this module In the case where the current PHASE of the D2P module is iphase 1 the descphase1 takes the form Table 69 Structure descphase1 NAMDIR mixdir MIX imix PKEY refnam i sapnam i i 1 npkey ENDPKEY FUEL BARR DEF unrodded aicg aicn USER unrodded compo i i 1 ncompo ENDBARR GRID SAP DEF USER GLOBAL pkey i nval i i 1 npkey ENDGLOBAL NEW ADD pkey i nval i val j j 1 nval i i 1 npkey ENDADD ADF DRA GET SEL
141. e axis TA case 77 11 Partial grid one complete plane and one complete axis with another configuration one point case 77 12 Face View of ACR Benchmark Core Model 292 Channels 158 13 Geometry definition plane 1 ooo e e 167 14 Top View of ACR Benchmark Core Model lt 170 15 Combustion zones definition io semet 2 bb bee Oe ee ee ea A we 181 IGE 344 1 1 INTRODUCTION DONJON is a full core modelization code designed around solution techniques of the neutron diffusion or simplified P equation The current DONJON package is an evolution version released as an attempt to introduce the innovative capabilities for the full core modeling and simulations of different types of nuclear reactors sush as Pressurized Water Reactors PWRs legacy CANDU reactors and Advanced CANDU Reactors ACRs The computer code DONJON Release 4 0 is part of Version4 distribution built around the GAN generalized driver The current DONJON package DONJON Version5 is a rewrite of the code around the GANLIB5 kernel intended to be 64 bit clean DONJON execution depends on other computer codes components of Version4 namely GANLIB UTILIB DRAGON and TRIVAC codes The DRAGON modules are used with DONJON code to define the reactor geometry to provide the macroscopic cross section libraries and to perform micro depletion calculations The TRIVAC solver modules are used to perform a spatial dis
142. e compared to refnam are expected name of state variable in the SAP object associated to refnam Default values are sapnam refnam NB If a state variable name is not correctly associated an error will occur during processing keyword used to indicate that the input SAP object contains cross sections for fuel assembly keyword used to associate an index of control rod in the SAP object to an index composition in PMAXS index of control rod in the SAP object for the unrodded cross section No default index of control rod in the SAP object for the aicg cross section No default index of control rod in the SAP object for the aicn cross section No default index of control rod in the SAP object for the composition i No default keyword used to associate a name of isotope in the SAP object to a specific isotope keyword used to indicate that the following record correponds to the name of Xe 135 in the SAP object name of Xe 135 isotope in the SAP or MCO object Default xenam XE135PF keyword used to indicate that the following record correponds to the name of Sm 149 in the SAP or MCO object name of Xe 135 isotope in the SAP or MCO object Default smnam SM149PF keyword used to indicate that the following record correponds to the name of I 135 in the SAP or MCO object name of Xe 135 isotope in the SAP or MCO object Default inam T1135PF keyword used to select the grid of state variables used for
143. e correction factor ju keyword used to set the first mixture where leakage parameters are differentiated By default the first mixture index is used minimum mixture index where leakage parameters are differentiated keyword used to set the last mixture where leakage parameters are differentiated By default the total number of mixtures in MACRO is used maximum mixture index where leakage parameters are differentiated keyword used to set the first energy group where leakage parameters are differentiated By default the first energy group index is used minimum energy group index where leakage parameters are differentiated keyword used to set the last energy group where leakage parameters are differentiated By default the total number of energy groups in MACRO is used maximum energy group index where leakage parameters are differentiated IGE 344 114 6 2 The GRAD module The GRAD module is designed to perform the following tasks compute the gradients of the system characteristics using solutions of direct or adjoint fixed source eigenvalue problems Here we assume an optimization problem with nvar control variables and with ncst constraints The total number of system characteristics is therefore equal to ncst 1 define options and parameters for the different method to solve the optimization problem The non linear optimization problem can be solved as a converging sequence of linear optimization problems with a
144. e of the quadratic constraint 57 where so 0 not converged 4 1 converged e The number of iterations relative to the quadratic constraint 2 e The type of reduction for the radius if the quadratic constraint Sg where o 1 half fa T parabolic e The number of inner iterations 57 e The number of outer iterations Sg e The resolution s method for the linear problem with quadratic constraint 58 where SIMPLEX LEMKE LEMKE LEMKE MAP Augmented Lagragian Penalty Method Sy UNIT e The number of outer iterations without step back Sfo e S not used e 52 not used A flag for unsuccessful resolution in module PLQ 53 where So 0 successful at last iteration 15 gt 1 number of iteration with unsuccessful resolution IGE 344 SIGNATURE uu STATE VECTOR DLEAK STATE VAR VALVE uu VAR MAX VAL VAR MIN VAL VAR WEIGHT uu CST 0BJ uuuuu CST TYPE uuu CST WEIGHTLU FOBJ CST VAL OPT PARAM R GRADIENT uuu OLD VALUE uu 149 Table 103 Main records and sub directories in optimize Nest I Nest R Nest R Nest 1 R 40 R Noar Nest ze 1 UnitsComment Signature of the data structure SIGNA Vector describing the various parameters as sociated with data structure S Vector describing the various parameters as sociated with data structure S7 This array is available if the OPTIMIZE object has been created using module DLEAK The val
145. e that the values of the D COOL record are to be updated They will be computed from the moderator temperature stored in the T C00L folder and the core pressure using the water tables Note that the position of the TTD keyword matters the density being calculated when TTD is called and not at the end of MCC execution core pressure in Pa IGE 344 33 3 9 The MOVDEV module The MOVDEV module can be used for the transient simulations and reactor control studies which are related to the time dependent rod devices displacement in the reactor core The rods can be inserted into or extracted from the reactor core at constant or at variable speed of movement The rod positions are recomputed at every given time step of movement The new rod positions can be computed in several ways based on either current time increment and movement speed relative change in rod positions or current rod insertion level The MOVDEV module allows the rod devices to be displaced individually or simultaneously in groups The MOVDEV module specification is Table 23 Structure MOVDEV DEVICE MOVDEV DEVICE descmove where DEVICE character 12 name of the DEVICE object that will be modified by the module The rods positions are updated according to the current time step of movement descmove structure describing the input data to the MOVDEV module 3 9 1 Input data to the MOVDEV module It is possible to move several individual rod
146. eactor channel BURN INSTouu R Nen No Instantaneous burnups per assembly subdivi sions for each channel POWER BUND y R Nen Np Powers per assembly subdivisions for each channel 8 1 5 The PARAM sub directories Each PARAM sub directory contains the information corresponding to a single local or global parameter excluding burnups Inside a such sub directory the following records will be found IGE 344 130 Table 85 Records in PARAM sub directories Condition Units Comment P NAME uuuuu Unique identification name of this parame ter This name is user defined however it is recommended to use the following pre defined values Boron concentration Averaged fuel temperature Surfacic fuel temperature Averaged coolant temperature Averaged coolant density CANDU only parameters T MODE Averaged moderator temperature PARKEY yuu Corresponding name of this parameter as recorded in a multi parameter Compo file P TYPE uuuuu Number associated to the type of recorded parameter ptype 1 for global parameter ptype 2 for local parameter P VALUE uu R 1 ptype 1 Recorded single value for global parameter R Ven No ptype 2 Recorded values for local parameter per each fuel bundle for every channel 8 1 6 The ASSEMBLY sub directory Table 86 Assembly type sub directory Condition UnitsComment LABEL uuuuuu Label of the assembly It consists a 8 char acter name composed of the AYNAME and AXNAME PIN POWER uuu
147. eactor geometry If FMAP is specified on the LHS its records BUND PW and FLUX AV will be set according to the information present in POWER FLUX character 12 name of the FLUX object previously created by the FLUD module The numerical flux solution contained in FLUX is recovered and all flux are normalized to the given total reactor power KINET character 12 name of the KINET object previously created by the KINSOL module The numerical flux solution contained in KINET is recovered TRACK character 12 name of the TRACK object created by the TRIVAT module The infor mation stored in TRACK is recovered and used for the average flux calculation MATEX character 12 name of the MATEX object containing the reactor material index and the h factors that will be recovered and used for the power calculation MACRO character 12 name of the MACROLIB object containing the h factors that will be recovered and used for the power calculation descflpow structure describing the input data to the FLPOW module 3 11 1 Input data to the FLPOW module Note that the fuel lattice power distribution can be printed only on the screen Table 27 Structure descflpow EDIT iprint PTOT power P NEW FSTH fsth INIT NORM BUND PRINT MAP DISTR FLUX RATIO POWER ALL where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print
148. ell calculations in DRAGON The HISTORY data structure can also be used to update the MAP data structure Finally the module HST can be used to create an initial BURNUP data structure that can be used to evolve the cell another time step in DRAGON The HST module can be used to create or update an HISTORY data structure The possible options are Table 40 Updating an HISTORY structure using a MAP structure HISTORY HST HISTORY MAP hstdim GET hstpar Table 41 Updating an HISTORY structure using a BURNUP structure HISTORY HST HISTORY BURNUP hstdim GET hstpar CELLID icha ibun idfuel GET hstpar It can also be used to create a BURNUP data structure from the information available on an HISTORY data structure Table 42 Updating a BURNUP structure using an HISTORY structure BURNUP HST HISTORY hstdim PUT hstpar CELLID icha ibun PUT BREFL hstbrn hstpar AREFL hstbrn hstpar AREFL hstbrn hstpar al It can also be used to update a MAP data structure from the information available on an HISTORY data structure Table 43 Updating an HISTORY structure using a MAP structure MAP HST MAP HISTORY where IGE 344 57 HISTORY character 12 name of an HISTORY data structure BURNUP character 12 name of a BURNUP data structure MAP character 12 name of a MAP data structure hstdim structure containing the dimen
149. ell number between 1 and Nc the maximum number of cells in the reactor Most of the information associated with the fresh fuel cells will be extracted from a DRAGON BURNUP file or defined using variable local parameters Each fresh fuel cell inserted in the core will also be associated with a specific fuel type Each fuel type is defined as a unique initial fuel composition The fuel management for the reactor including burnup and refueling will be performed by the DONJON code Here the HST will interact with this code via the MAP data structure Typically each cell in the reactor will be burned inside DRAGON using the power provided in the AX SHAPE record and the depletion time provided in the BURNUP BEG record stored in the MAP structure When refueling takes place some of the fuel cells will be extracted other will be displaced from one position to another and finally new fresh fuel cells inserted The fresh fuel cells properties will be extracted from the fuel types properties available on the HISTORY data structure In a coupled DRAGON DONJON execution the HST module will be called at various points and for various reasons The first call to HST can be performed using MODULE HST Enola Map data structure for initialization MAPO History data structure History E SEG ASCII MAPO XSM_FILE History XSM_FILE Reseau gui Reactor parameters IGE 344 60 ncha nunber of channels 9 nbun nunber of bundles 3 nevo
150. er 2 array of horizontal channel names A horizontal channel name is iden tified by the channel column using numerical characters 1 2 3 and so on Note that the total number of X names must equal to the total number of subdivisions along the X direction in the fuel map geometry All non fuel regions are to be assigned a single character This option is not available for 3 D Hexagonal geometry When assembly are defined and split several names can be the same integer total number of subdivisions along the X direction in the fuel map geometry Not used for 3 D Hexagonal geometry keyword used to specify YNAME Not used for 3 D Hexagonal geometry IGE 344 YNAME ny NCOMB ncomb B ZONE icz ALL nch nb ASBLY SIM Ix ly naval 11 character 2 array of vertical channel names A vertical channel name is identified by the channel row using alphabetical letters A from the top B C and so on The total number of Y names must equal to the total number of subdivisions along the Y direction in the fuel map geometry All non fuel regions are to be assigned a single character This option is not available for 3 D Hexagonal geometry When assembly are defined and split several names can be the same integer total number of subdivisions along the Y direction in the fuel map geometry Not used for 3 D Hexagonal geometry keyword used to specify the number of combustion zones
151. er of global and local parameters used in a HISTORY data structure can be increased at all time The number of channels bundles and the refueling scheme must be defined at the creation of the HISTORY data structure This information can be provided manually or extracted from a MAP data structure The general form of the hstdim input structure follows Table 44 Structure hstdim EDIT iprint DIMENSIONS GLOBAL nglo LOCAL nloc BUNDLES nbun bunl CHANNELS ncha where EDIT keyword used to modify the print level iprint iprint index used to control the printing in this module It must be set to 0 if no printing on the output file is required DIMENSIONS keyword used to indicate that the general dimensioning of the HISTORY data struc ture will be modified IGE 344 58 GLOBAL keyword used to modify the number of global parameters on the HISTORY data structure nglo the number of global parameters Note that the history module will use the maxi mum value between the current nglob and the value if any defined on the HISTORY data structure LOCAL keyword used to modify the number of local parameters on the HISTORY data struc ture nloc the number of local parameters Note that the history module will use the maximum value between the current nloc and the value if any defined on the HISTORY data structure BUNBLES keyword used to specify the number of bundles per channels for the reactor model considere
152. erformed This name must be set in the RHS of the SCR data structure name of the parameter for burnup or irradiation This value is defined if option TABLE is set and if burnup or irradiation is to be considered as parameter input data structure containing interpolation information relative to the SAPHYB data structure named SAPNAM see Section 4 3 2 IGE 344 80 descdepl input structure describing the depletion chain see Section 4 3 3 This input structure reguires option MICRO By default the depletion chain data is not written in the output MICROLIB 4 3 2 Defining global parameters If a MAP object is defined on the RHS of structure scr_data and if the TABLE keyword is set some information required to set the interpolation points is found in this object In this case the SCR operator search the SAPHYB object for global parameters having an arbitrary name specified in the MAP object or set directly in this module Note that any parameter s value set directly in this module prevails on a value stored in the MAP object Each instance of descints is a data structure specified as Table 59 Structure descints MIX imix FROM imixold USE TIMAV BURN INST BURN AVG EX BURN ivarty SET DELTA ADD LINEAR CUBIC PARKEY vall MAP 4 val2 MAP REF PARKEY valref SAMEASREF ENDREF MICRO ALL ONLY HISO 4 conc ENDMIX where M
153. eriod of time delt keyword used to specify that a particular rod or a group of rods will be extracted from the reactor core during the period of time delt keyword used to specify the new level value real positive value of the rod insertion level at current time step This value will be used to compute the new rod position in the reactor core The insertion level is minimal value 0 0 when the rod is completely withdrawn and it is maximal value 1 0 when the rod is fully inserted For the partially inserted rod the insertion level must be 0 0 lt value lt 1 0 keyword used to specify the value delh real positive absolute value of the relative change in the rod position during the period of time delt This is a time dependent rod displacement along the rod movement axis which must be given in cm keyword used to set the current value of speed real positive absolute value of the rod movement speed given in cm s The rod speed can be kept constant or it can be modified at any time step delt The devices could also have the different speeds of movement IGE 344 35 3 10 The NEWMAC module The NEWMAC module is used to create a complete MACROLIB with respect to the devices parameters The resulting MACROLIB will contain the exact properties for every material region over the whole mesh splitted reactor geometry The material properties of each region are recomputed with respect to the actual position of each rod type and
154. erturbed FMAP object in which the fuel type mixtures indices are modified according to the specified core voiding pattern The information with respect to the relative coolant den sities is required only for the subsequent interpolation of fuel properties using the NCR module These data will also be reordered by the CVR module according to the specified voiding pattern and recorded as local parameter in the perturbed fuel map object see Section 3 1 2 The CVR module specification is Table 38 Structure CVR CVR FMAP descrevr where FMAP character 12 name of a read only FMAP object created in the RESINI module This object must contain the non perturbed fuel cell properties FMAPV character 12 name of a new FMAP object that will contain the modified fuel type indices and reordered coolant densities according to the specified core voiding pattern descrcvr structure describing the input data to the CVR module 3 17 1 Input data to the CVR module Note that the input order must be respected Table 39 Structure descrevr EDIT iprint MIX FUEL mixF i MIX VOID mixV i i 1 nfuel DENS COOL PNAME SET dcoolV VOID PATTERN FULL HALF QUARTER CHECKER CHECKER 1 2 CHECKER 1 4 CHAN VOID nvoid YNAME i XNAME i i 1 nvoid 3 where EDIT keyword used to set iprint IGE 344 iprint MIX FUEL mixF MIX VOID mixV DENS COOL PNAME SET dcoolV VOID PATTERN
155. esian geometry of the fuel lattice continued on next page IGE 344 Records and sub directories in fmap data structure Condition Units FUEL uuuuuuu Dir Neue hcycle Dir Nburn Tsim 0 PARAM junuuua Dir Nparm Nparm gt 0 128 continued from last page Comment List of fuel type sub directories Each com ponent of the list is a directory containing the information relative to a single fuel type Sub directory containing information related to a fuel cycle in a PWR Nhurn is the number of burnup steps used during the simulation of the cycle These burnup steps may not be of increasing values List of parameter type sub directories Each component of the list is a directory containing the information relative to a single parame ter The total number of sub directories cor responds to the total number of recorded pa rameters Nparm excluding burnups The contents of the GEOMAP sub directory correspond to the typical contents of the geometry data structure The dimensioning parameters Nx Ny and N represent the number of volumes along the corresponding axis in the fuel map geometry The shifting information records pshift bshift and ishift will be composed using the following FORTRAN instructions respectively as WRITE pshift A6 12 PSHIFT ell WRITE bshift A6 12 BSHIFT ell WRITE ishift A6 12 ISHIFT ell for 1 lt ell lt Nene Each ti
156. et in a previous call to module GRAD IGE 344 117 cstval array containing ncst real values CST WEIGHT keyword to specify the weights or penalties of the constraints These weights are not used with Lemke or MAP methods These values can also be set in a previous call to module GRAD cstw array containing ncst real values IGE 344 118 6 3 The PLQ module The PLQ module is used to solve the linear programming problem with a quadratic constraint The gradients of the system characteristics are calculated with module GRAD The options and parameters for the different method to solve the optimization problem are also defined in module GRAD The calling specifications are Table 76 Structure PLQ OPTIM PLQ OPTIM plq_data where OPTIM character 12 name of the OPTIMIZE object L_OPTIMIZE signature containing the optimization informations Object OPTIM must appear on the RHS to be able to updated the previous values plq_data structure containing the data to the module PLQ see Sect 6 3 1 6 3 1 Data input for module PLQ Table 77 Structure plq_data EDIT iprint WARNING ONLY CALCUL DX NO STORE OLD COST EXTRAP gt gt ecost lt lt CONV TEST gt gt leony lt lt where EDIT keyword used to set iprint iprint index used to control the printing in module WARNING ONLY keyword used to specify that only a warning will be used when no valid previous decision ve
157. evices If the device rod insertion or lzc filling level is selected for the modification then a new device position is recomputed accordingly The DSET module can be used to perform the device reactivity studies and also to predict the reactivity worth of the rod devices The DSET module specification is Table 19 Structure DSET DEVICE DSET DEVICE descdset where DEVICE character 12 name of the DEVICE object that will be updated by the module descdset structure describing the input data to the DSET module 3 7 1 Input data to the DSET module It is possible to set or to modify the parameters for several individual devices and or for several groups of devices simultaneously Table 20 Structure descdset EDIT iprint ROD irod ROD GROUP irgrp LZC ilzc LZC GROUP ilgrp LEVEL value SPEED speed TIME time END 2 where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing larger values produce increasing amounts of output ROD keyword used to specify the rod irod number irod integer identification number of a rod to be modified Each rod type device has a IGE 344 ROD GROUP irgrp LZC ilzc LZC GROUP ilgrp LEVEL value SPEED speed TIME time END 30 unique irod number ranging from 1 to nrod as been defined in the DEVINI module see Section 3 4 2 ke
158. eyword used to indicate the input of poison load in a fuel real value of the fuel type parameter specified for each fuel type in the same order as the fuel mixture indices have been recorded in the MATEX object see Section 3 2 1 integer total number of the fuel types as been defined in the USPLIT module keyword used to specify that a patterned age distribution will be input and used to compute instantaneous bundle burnup real array containing the refueling sequence numbers This channel is refueled the ialch i th one The channels are ordering from the top left to the bottom right of the core The expression of the resulting bundle burnups are given in Ref 20 IGE 344 15 3 2 The USPLIT module The USPLIT module is used to create a MATEX object that will provide a link between the reactor geometry and material index The 3 D Cartesian or 3 D Hexagonal reactor geometry which is previously produced in the GEO module is analyzed and the material mixture indices are recomputed in order to provide a unique mixture number for each material sub volume Such renumbering permits a complex reactor core modeling A MATEX object is also used to store some additional information that will be required and updated by other DONJON modules see Section 8 2 The USPLIT module specification is Table 5 Structure USPLIT GEOM MATEX USPLIT GEOM GEOMOLD desclink where GEOM character 12 name of a GEOMETRY object This objec
159. for the whole geometry MAX T FUEL ji MAX D C00Luu MAX T C00Luu HISTORY DATA J v yy Y R The signature for this data structure is SIGNA L_THM The array S contains the following informa tion e S contains the number of active fuel rods e St contains the number of guide tubes e S contains the maximum number of iterations for computing the conduction integral e S contains the maximum number of iterations for computing the center pellet temperature e S contains the maximum number of iterations for computing the coolant parameters velocity pressure enthapy density in case of a transient calculation e St contains the number of discretisation points in fuel e S contains the number of total discretisation points in the whole fuel rod fuel cladding S contains the integer setting the type of calculation steady state or transient performed by the THM module e Si contains the current time index e S flag to set the gap correlation Sth 0 built in correlation is used 10 1 set the heat exchange coefficient of the gap as a user defined constant e S flag to set the heat transfer coefficient between the clad and fluid sth 0 built in correlation is used 11 1 set the heat exchange coefficient between the clad and fluid as a user defined constant e St flag to set the fuel thermal conductivity Sth 0 built in correlation is used 12 1 set the fuel thermal conducti
160. fy that the step will be reduced by a factor of 2 keyword used to specify that the step will be reduced with the parabolic method keyword to specify the values of the control variables These values can also be set in a previous call to module GRAD or set in another module array containing nvar real values keyword to specify the values of the control variable weights in the quadratic constraint All weights are set to 1 0 by default array containing nvar real values keyword to specify the minimum values of the control variables These values can also be set in a previous call to module GRAD array containing nvar real values single real value used for all control variables keyword to specify the maximum values of the control variables These values can also be set in a previous call to module GRAD array containing nvar real values single real value used for all control variables keyword to specify the value of the objective function followed by the actual values of the constraints These values can also be set in a previous call to module GRAD or set in another module array containing ncst 1 real values keyword to specify the relation types of the constraints These values can also be set in a previous call to module GRAD array containing ncst integer values These values are 1 for gt 0 for equalily and 1 for lt keyword to specify the RHS values of the constraints These values can also be s
161. gan 2013 IGE 344 Index ext 115 116 Einn 115 116 Equad 115 116 leonv 118 119 BUND PW 37 gt FLUX AV 37 1135PF 102 SM149PF 102 XE135PF 102 default 102 159 163 descphase1 100 101 compo i i 1 ncompo 101 emin g g 1 ngroup 101 pkey i nval i i 1 npkey 101 pkey i nval i val j j 1 nval i i 1 npkey 101 65 66 68 69 71 79 80 82 46 DRA 152 dra 97 inp 97 res 38 OLD VALUE 150 22 8 15 17 18 22 25 28 31 33 25 36 38 42 45 49 51 53 56 63 67 78 88 93 100 107 112 114 118 120 O 8 15 17 18 22 25 29 31 33 85 36 38 42 45 49 51 53 56 63 67 78 88 93 107 112 114 118 120 D2P 100 6 7 51 93 BURNUP 56 57 HISTORY 56 57 MAP 56 57 D2P 100 descphase1 101 3 D 1 3 4 8 10 13 15 16 18 19 25 126 3D 127 131 133 136 1 100 0001 147 2 100 3 100 DETECT 51 INIDET 22 A ZONE 9 10 ADD 31 32 69 71 80 81 83 101 103 ADD PARAM 12 13 187 ADF 101 103 AFM 4 11 12 17 38 40 43 88 aicg 101 102 aicn 101 102 ALL 9 11 21 28 37 38 43 44 65 66 69 70 80 81 115 ALLX 67 68 and or 36 AREFL 56 57 62 ASBLY 9 12 15 16 asmbl 46 47 asmb2 46 47 ASSEMBLY 9 10 ASSMB 108 109 AT 46 47 AUG LAGRANG 115 avburn 46 47 AVG EX BURN 65 69 80 AVGB 46 47 AX SHAPE
162. given in m s The value of speed is required only for the reactor regulating purpose keyword used to specify a new value for time real value of time either for the rod insertion or extraction or for the liquid controller filling given in sec The value of time is required only for the reactor regulating purpose keyword used to indicate the end of input of the new parameters for the current device or group of devices IGE 344 31 3 8 The MCC module The MCC module supplies tools to edit selectively one or several parameters of a fuel map object These parameters are the ones defined when the fuel map is created according to the calculation grid chosen to compute the macroscopic cross sections libraries Their selective edition is useful to study their impact on the core reactivity For instance this module enables to increase uniformly the fuel tempera ture in each cell without modifying any other parameter such as the moderator temperature It grants access to the Doppler coefficient even at hot full power when the fuel temperature is different in every cell of the core This module enables the computation of non exhausive list the Doppler coefficient the power Doppler coefficient Doppler during a power transient the moderator temperature coefficient the boron concentration coefficient the reactivity with the fuel temperature set to the moderator one etc Note that this module does not perform any reactivity calculatio
163. hat MLIB is a MACROLIB default option keyword indicating that MLIB is a MICROLIB Object MLIB contains an embedded MACROLIB but the CPU time required to obtain it is longer keyword indicating that interpolation of the SAPHYB uses linear Lagrange polynomials keyword indicating that interpolation of the SAPHYB uses the Ceschino method with cubic Hermite polynomials as presented in Ref 17 default option keyword used to introduce leakage in the embedded MACROLIB This option should only be used for non regression tests the imposed buckling corresponding to the leakage keyword used to select a SPH factor set in the Saphyb By default the cross sections and diffusion coefficients are not SPH corrected character 80 name of the SPH factor set selected in the Saphyb These sets are stored as local parameter information in the Saphyb keyword used to define the maximum number of material mixtures This information is required only if MLIB is created the maximum number of mixtures a mixture is characterized by a distinct set of macroscopic cross sections the MACROLIB may contain The default value is nmixt 0 or the value recovered from MLIB if it appears on the RHS keyword used to set SAPNAM and to define each global parameter keyword used to set SAPNAM and to recover some global parameter from a MAP object named MAPFL character 12 name of the LCM object containing the SAPHYB data structure where the interpolation is p
164. he nucleus 0 for a nucleus in its ground state 1 for an isomer in its first exited state etc For example U in its ground state will be represented by izae 922380 reaction character 6 identification of a neutron induced reaction that takes place either for pro duction of this isotope its depletion or for producing energy Example of reactions are following NG indicates that a radiative capture reaction takes place either for produc tion of this isotope its depletion or for producing energy N2N indicates that the following reaction is taking place n 4 Xz gt 2n Xz N3N indicates that the following reaction is taking place n 4 Xz gt 3n 4 Xz IGE 344 83 N4N indicates that the following reaction is taking place n 4 Xz gt 4n 4 3 Xz NP indicates that the following reaction is taking place pits phe Veg NA indicates that the following reaction is taking place n LS Xe 4 Hey L NFTOT indicates that a fission is taking place energy energy in MeV recoverable per neutron induced reaction of type reaction If the energy associated to radiative capture is not explicitely given it should be added to the energy released per fission By default energy 0 0 MeV STABLE non depleting isotope Such an isotope may produces energy by neutron induced reactions such as radiative capture FROM indicates that this isotope is produced from decay or neutron induced reactions DECAY indicates that a decay reaction takes pl
165. his number must be greater than 0 dev lzc structure describing the input data for each individual liquid controller lzc group structure describing the input data for each group of liquid controllers 3 6 2 Description of dev lzc input structure Note that the devices positions can not overlap in the reactor core The input order of data must be respected Table 17 Structure dev lzc LZC id MAXPOS pos i i 1 6 MAX FULL fmax ants X Y 2 LEVEL value RATE rate TIME time EMPTY MIX mixE n n 1 2 FULL MIX mixF n n 1 2 IGE 344 where LZC id MAXPOS pos MAX FULL fmax AXIS X Y Z LEVEL value RATE rate TIME time EMPTY MIX mixE FULL MIX mixF 27 keyword used to specify the liquid controller id number integer identification number of the current liquid controller Each controller must be assigned a unique id number given in an ascending order ranging from 1 to nlzc keyword used to specify the entire position of a liquid zone controller including its empty and full parts real array containing 3 D Cartesian coordinates of the liquid zone controller position in the reactor core These coordinates must be given in the order X X Y Y Z Z keyword used to specify fmax real value of the limiting coordinate along the controller filling axis which corresponds to the maximum full filling level for the current liquid controlle
166. i 1 nch ALL ASBLY SIM Ix ly naval i i 1 nch descresini2 where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value larger values produce increasing amounts of output keyword used to indicate the call to an embedded module SPLIT NAP keyword to specify that the embedded geometry will be split by the embedded NAP module IGE 344 GEO descgeo NAP ASSEMBLY na nax nay A ZONE iza ASBLY AXNAME XNAMEA AYNAME YNAMEA NXNAME XNAME nx NYNAME 10 keyword used to call the GEO module The fuel map geometry differs from the com plete reactor geometry in the sense that it must be defined as a coarse geometry i e without mesh splitting over the fuel bundles Consequently the mesh spacings over the fuel regions must correspond to the bundle dimensions e g hz width hy height h length or in 3 D Hexagonal geometry h side h height Note that the total number of non virtual regions in the embedded geometry must equal to the number of fuel channels times the number of fuel bundles per channel This means that only the fuel type mixture indices are to be provided in the data input to the GEO module for MIX record Other material regions e g reflector must be declared as virtual i e with the mixtures indices set to 0 structure describing the input data to the GEO
167. ich can be used with DONJON code is given in Section 2 2 In this section a detailed description of the DONJON data structures is presented 8 1 Contents of fmap data structure A fmap data structure is used to store fuel assembly or bundle map and fuel information such as powers average fluxes control zones burnup or refueling scheme The fuel bundle location are given in an embedded sub directory which contains the records as a geometry data structure This object has a signature L_MAP it is created using the RESINI module 8 1 1 The state vector content The dimensioning parameters S which are stored in the state vector for this data structure represent e The number of fuel bundles per channel N Si e The number of fuel channels Nen Sa e The number of combustion zones Neomp S3 e The number of energy groups Ner S4 e The type of interpolation with respect to burnup Iptyp Ss 0 interpolation type is not provided Ibtyp 1 according to the time average model 2 according to the instantaneous model e The number of bundle shift Nani S e The number of fuel types Niue S7 e The number of recorded parameters Nparm Sg e The total number of fuel bundles Miot S9 e The number of voided reactor channels Nyoia S10 e The option with respect to the core voiding pattern Lua S11 voiding pattern not provided full core voiding pattern half core voiding pattern quarter core voiding pattern checkerbo
168. ide for TRIVAC Version4 Report IGE 293 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2007 A H bert TRIVAC A Modular Diffusion Code for Fuel Management and Design Applications Nucl J of Canada Vol 1 No 4 325 1987 A H bert Application of the Hermite Method for Finite Element Reactor Calculations Nucl Sci Eng 91 34 1985 A H bert Variational Principles and Convergence Acceleration Strategies for the Neutron Diffusion Equation Nucl Sci Eng 91 414 1985 A H bert Preconditioning the Power Method for Reactor Calculations Nucl Sci Eng 94 1 1986 A H bert Development of the Nodal Collocation Method for Solving the Neutron Diffusion Equa tion Ann Nucl Energy 14 527 1987 A H bert Application of a Dual Variational Formulation to Finite Element Reactor Calculations Ann nucl Energy 20 823 1993 E Varin A H bert J Koclas and R Roy A User Guide for DONJON Report IGE 208 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2005 E Varin A H bert Data Structures for DONJON Report IGE 226 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2003 D Rozon Gestion du combustible nucl aire Notes de cours Ene6109 Report IGE 298 Ecole Polytechnique de Montr al Institut de G nie Nucl aire 2007 R Chambon Optimisation de la gestion du combustible dans les
169. ify all the different possibilities and limitations dependently of the database structure we will use a 3 parameter case The paramaters are referenced by A B and C But before we explain the different cases we want to remind that the interpolation factors are computed on each axis seperatly The first case corresponds to a complete grid represented by a gray paralepiped on Fig 2 and 3 The figure 2 shows that the interpolated value in point V can be obtained directly without MAP object For time average TA computation lets assume that the parameter B represents the burnup and keep this convention for other database structure also In this case the figure 3 shows also that the direct interpolation can be done to compute an average value between the points V and V Note that the TA burnups are stored in the MAP object and are then recovered automatically The second case corresponds to a partial grid where all the lattice computations have been perfomed for several pairs of parameters which are represented as the gray rectangles on Fig 4 and 5 If we use the notations of Fig 4 and 5 the best estimate interpolated values f we can get are given by f F V f Va VBa f Ve WVac f VB f Vsc f Vea f VB Vea f Vec f Va for instataneous f fV V F Vg Va HVga Vea f Vp Va HVgo Vac MVE Ve Vac Vac f VBa VBa S Vg Va Vba Vea F Vbo Vec
170. ill have the power values set to 0 keyword used to indicate the printing of all available information i e without partic ular selection of data keyword to specify that the output flux object will contain a value per mesh splitted element normalized to the given power as required by the DETECT module This is the default option keyword to specify that the output flux object will contain a value per bundle nor malized to the given power IGE 344 39 3 12 The TAVG module The TAVG module is used to compute the burnup integration limits for each fuel bundle the axial power shape over the fuel lattice the channel refuelling rates and the reactor core average exit burnup All calculations using the TAVG module are performed according to the time average model for the equilibrium core conditions The computing algorithm is based on bidirectional refuelling schemes of channels and average exit burnups specified over the fuel lattice which should be recorded in the fuel map using the RESINI module Note that the complete time average calculation is a complex and iterative procedure requiring of several full core calculations external iterations to be performed The main steps of the time average calculation using DONJON are briefly described at the end of this section The TAVG module can also be used to compute the instantaneous fuel burnups according to the channel patterned age model for the fuel management and optimization purp
171. inate position of an assembly identified with a SPC keyword Up to 30 coordinates can be set aside if many assemblies have the same specification keyword indicating that a user defined value will be assigned to the assembly keyword indicating that an averaged burnup will be assigned to the assembly real value of the average burnup in MWd t keyword indicating that a new fuel assembly will be used integer index of the fuel type corresponding to the new fuel assembly Fuel type indices are those used in the RESINI PLANE descriptions of Sect 3 1 keyword indicating that a value recovered from another assembly will be assigned to the current assembly character 12 identification name of a previous fuel cycle keyword indicating that the naval coordinate position of the other assembly will be given character 3 naval coordinate position of the other assembly in cycle hcold2 keyword used to impose an axial burnup distribution to the assembly The burnup distribution is recovered from an existing assembly or is set to user suppled values real values of the axial burnup distribution keyword used to indicate an increase of core average burn up keyword used to indicate in MWd t the average increase of burn up in the core keyword used to indicate the time of combustion at the power specified in POWER structure keyword used to set the time combustion value in DAY or HOUR or MINUTE or SECOND keyword used to specify that rtime is a nu
172. ing current version of DONJON correspond to the reactor static conditions The modeling of the reactor fuel lattice using DONJON is made in considering that the fuel lattice is composed of a well defined number of fuel channels and bundles All reactor channels contain the same number of fuel bundles and are identified by their specific names The fuel bundles have a distinct set of properties that are recovered and interpolated according to the specified global and local parameters The interpolation of fuel properties with respect to burnup distribution can be performed according to the time average or instantaneous models The time average calculation is performed in considering the bidirectional refuelling scheme of reactor channels and assuming that all channels have the same bundle shift The modeling of the reactivity mechanisms is based on their specified parameters which include the devices position rods insertion level water filling level direction of movement etc The rod devices insertion level can be set according to their nominal positions or they can be displaced in and out of core The devices can also be divided into several groups so that they can be manipulated displaced or moved simultaneously The time dependent behaviour of the moving devices can be modeled and used for the transient simulations or reactor control studies The reactivity worth of devices can also be studied and predicted using DONJON The reactor material pr
173. ing records and sub directories will be found in the first level of a thm directory SIGNATURE uu STATE VECTOR REAL PARAM ji KCONDF Guu UCONDF uu UCONDC LULU ERROR T FUEL ERROR D COOL Table 100 Main records and sub directories in thm Condition Units Comment K parameter SIGNA containing the signature of the data structure array S containing various integer parameters that are required to describe this data structure array R containing various floating point parame ters that are required to describe this data structure coefficients of the user defined correlation for the fuel thermal conductivity string variable set to KELVIN or to CELSIUS coefficients of the user defined correlation for the clad thermal conductivity string variable set to KELVIN or to CELSIUS absolute error in fuel temperature g cc absolute error in coolant density continued on next page IGE 344 144 Main records and sub directories in thm continued from last page Condition Units Comment ERROR T COOL MIN T FUEL MIN D CO0L jj MIN T CO0L absolute error in coolant temperature minimum fuel temperature minimum coolant density minimum coolant temperature maximum fuel temperature maximum coolant density maximum coolant temperature sub directory containing the historical values taken by the thermal hydraulics parameters mass flux density pressure enthalpy temperature in the coolant and in the fuel rod
174. ing scheme The absolute value of nsh is the number of fuel bundles inserted in the channel NAMCHA The sign of nsh define the refueling direction positive direction is from the first to the nk th bundle and negative is from the nk th to the first bundle NEWFUEL key word to specify that a channel will be refueled with a different type of fuel SOME key word to specify that the nsh values of fuel types can be different imix i index number of a fuel type with respect to the values defined in module NCR CRE or AFM IGE 344 ALL SHUFF CHAN NMCHA1 TO NMCHA2 POOL 44 key word to specify that the nsh values of fuel types will be identical to imix key word to specify that a specified channel will move into an other one or discharge into the pool key word to specify the moved channel name channel name as defined by NXNAME and NYNAME It is constructed as NAMCHA key word to specify the bundle destination channel name as defined by NXNAME and NYNAME It is constructed as NAMCHA key word to specify that the channel referenced by NMCHA1 is discharged into the pool IGE 344 45 3 14 The SIM module The SIM module can perform a sequence of operations related to fuel management in PWRs e simulate a refuelling and shuffling scheme and update the burnup distribution accordingly The refuelling scheme is specified directly in SIM e increase the burnup using the power available in the POWER object and compute the fi
175. ization values Example homoge neous calculation followed by a pin power reconstruction where assemblies were not defined in the first place IGE 344 BUNDLE POW bvalue REF SHIFT COMB ishift ADD PARAM PNAME PNAME PARKEY PARKEY GLOBAL LOCAL SET PARAM SAME 13 e with a different geometry In this case the assembly geometry of the new FLMAP and the geometry of the FLMAP2 must match Example homogeneous calcula tion followed by a heterogeneous calculation or pin power reconstruction keyword used to specify the power values for each fuel bundle This option is not available in 3 D Hexagonal geometry real array containing the burnups powers values given in MW day per tonne MW of initial heavy elements The fuel burnup is considered as a global parameter keyword used to specify ishift Note that the axial power shape and the first burnup limits will be reinitialized each time the channel refuelling schemes are modified by the user This option is not avaialble in 3 D Hexagonal geometry keyword used to indicate the input of bundle shift numbers per combustion zone integer array or single value of the bundle shift numbers A single ishift value means that the same bundle shift will be applied for all combustion zones Note that the bundle shift value must be positive it corresponds to the number of displaced fuel bundles during each channel refuelling keyword used to indicate the input of infor
176. l be taken If this option is omitted and MAP contains bundle fluxes new Xenon concentrations will be computed and used keyword used to set nNp Neptunium concentration in 10 4at em IGE 344 XEREF SAM nSm IMET imet BLIN 91 keyword used to specify that the Neptunium concentrations as computed with DRAGON will be taken If this option is omitted and MAP contains bundle fluxes new Neptunium concentrations will be computed and used keyword used to set nSm Samarium concentration in 10 at cm If this data is omitted bundle concentrations as computed by DRAGON is used keyword used to set imet interpolation type for time average calculations imet 1 using Lagrange approxima tions imet 2 using spline approximations imet 3 using Hermite approximations default value keyword used to linear interpolation for burnup instead of the Lagrangian interpolation method Note The concentration of boron is provided in terms of 10 4at cm in the database However the usual units are ppm wt of Boron Thus the input asks for ppm of Boron ng and automatically transform the units into 10 at cm using the following equations pe g em NB Pwater g cm and 3 3 3 N 3 Pwaterlat cm 3Pwater molecule cm gg Prater 9 em water 3 3 N 3 ps at fem pp molecule em ma Pig em B thus Mwa er pg 0 at cm np a Pwater 102 at cm3 3 Mp where M molar mass and N the Avogadro number
177. l data structure that will be stored on sequential binary files SEQ_ASCII keyword used to specify the names of all data structures that will be stored on sequen tial ASCII files STRNAME character 12 name of a data structure The list of data structure that can be used in DONJON is presented in Section 2 2 module input specification for a module that will be executed For DONJON specific modules these input structures are described in Section 3 to Section 6 END keyword to call the normal end of execution utility module i keyword to specify the end of record This keyword is used to delimit the part of the input data stream associated with each module Generally the user has the choice to declare the most of data structure in the format of a linked list to reduce CPU times or as a XSM file to reduce memory resources In general the data structure are stored on the sequential ASCII files only for the backup purposes The input data normally ends with a call to the END module However the GAN driver will insert automatically the END module even if it was not provided upon reaching an end of file in the input stream Each module calling specification contains a module execution description and its associated input structure All these modules except the END module may be called more than once IGE 344 8 3 GENERAL CORE DESCRIPTION MODULES 3 1 The RESINI module The RESINI module is used for modeling of the reactor fue
178. l lattice in 3 D Cartesian geometry or 3 D Hexagonal geometry This modeling is based on the following considerations e For 3 D Cartesian geometry the reactor fuel lattice is composed of a well defined number of fuel channels Each channel is composed of a well defined number of fuel bundles or assembly subdivi sions All channels contain the same number of fuel bundles or assembly subdivisions Each reactor channel is identified by its specific name which corresponds to its position in the fuel lattice In a Candu reactor the channels are refuelled according to the bidirectional refuelling scheme The refuelling scheme of a channel corresponds to the number of displaced fuel bundles bundle shift during each channel refuelling The direction of refuelling corresponds to the direction of coolant flow along the channel In a PWR a basic assembly layout can be projected over the fuel map using a naval coordinate position system Assembly refuelling and shuffling will be possible using the ad hoc module SIM see Section 3 14 e For 3 D Hexagonal geometry the reactor fuel lattice is composed of a well defined number of fuel channels and each channel is composed of a well defined number of fuel bundle All fuel bundles have the same volume All channels contain the same number of fuel bundles Refuelling is not available during the calculation The lattice indexation is kept to identify the hexagons e The fuel regions generally have a different
179. l map see Section 3 1 2 QMAP keyword defining the assembly layout in naval coordinate positions using quarter core symmetry conditions Here the lower right quarter is defined The full map is recon structed through rotations around the center hx ordered list of available character 1 prefixes for the X oriented naval coordinate po sitions Values are generally chosen between A and T hy ordered list of available character 2 suffixes for the Y oriented naval coordinate po sitions Values are generally chosen between 01 and 17 hase character 4 identification value for the i j position Accepted values are e or for a position outside the core e NEW for a new assembly at zero burnup selected according to the fuel map specified in Sect 3 1 IGE 344 SPEC asmbl SET AVGB avburn FUEL fuel FROM hcold2 AT asmb2 DIST AX axn BURN STEP rburn TIME rtime DAY HOUR MINUTE SECOND ENDCYCLE COMPARE hel he2 47 e SPC for an assembly described later in the dataset using a SPEC specification e or a naval coordinate position referring to the position of an assembly in cycle hcold that is recycled in the current cycle keyword defining specifications related to all assemblies previously identified with the SPC keyword If QMAP keyword has been used with SPC values the 4 equivalent assem blies must be specified i e not only the lower right quarter assembly character 3 naval coord
180. long x axis of the reactor geometry The mesh splitted coordinates along y axis of the reactor geometry The mesh splitted coordinates along z axis of the reactor geometry The h factors per each mixture and per each energy group as recovered from the extended macrolib data structure IGE 344 133 8 3 Contents of device data structure A device data structure is used to store several information related to the reactor devices This object has a signature L_DEVICE it is created using the DEVINI module The information contained in this data structure can be used and updated in other DONJON modules which are related to the devices namely LZC DSET MOVDEV and NEWMAC modules 8 3 1 The state vector content The dimensioning parameters S which are stored in the state vector for this data structure represent e The type of reactor geometry I S only I 7 for 3D Cartesian geometry allowed The total number of rod type devices N ogq Sa The total number of the rod type groups Ny grp 93 e The total number of lzc type devices Nize S4 The total number of the lzc type groups Nigrp 5 Type of control rod movement Imoy Sg where ae 1 Fading A fraction of the fully inserted rod vanishes mov 2 Moving The complete rod is moving DONJON3 type movement 8 3 2 The main device directory The following records and sub directories will be found on the first level of device directory Table 88 Recor
181. lt value is 0 05 keyword used to set the plutonium mass enrichment of fuel Plutonium enrichment affects some built in correlations used to represent the heat conduction phenomenon in fuel real value set to the plutonium mass enrichment The default value is 0 0 IGE 344 CONDF ncond kcond unit INV inv ref CONDC HGAP hgap HCONV hconv TEFF wtelf 110 keyword used to set the fuel thermal conductivity as a function of local fuel temperature True Fuel conductivity is computed as ncond inv A fuel kcond k T puei fuel gt cond k T fuel Fy ret with Afue in W m K and True in the selected unit Kelvin or Celsius By default built in models are used taking into account oxyde porosity and pluto nium mass enrichment Note that oxyde porosity and plutonium mass enrichment are ignored if this keyword is used integer value set to the degree of the conductivity polynomial real value set to the coefficient of the conductivity polynomial ncond 1 coefficients are expected string value set to the unit of temperature T in the conductivity function Can be either CELSIUS or KELVIN keyword used to add an inverse term in the fuel conductivity function real value set to the coefficient in the inverse term of fuel conductivity The default value is 0 0 i e no inverse term real value set to the reference in the inverse term of fuel conductivity keyword used to set the
182. luxes are normalized either to the given total reactor power or using the previously recorded normalization factor The recorded values of the maximum channel and bundle powers the channel and bundle power form factors and the effective multiplication factor can be used as power and criticity constraints for the optimization and fuel management purposes 8 6 Contents of history data structure This data structure contains the information required to ensure a smooth coupling of DRAGON with DONJON when a history based full reactor calculation is to be performed 8 6 1 The main directory The following records and sub directories will be found in the first level of a history directory Table 97 Main records and sub directories in history Condition Units Comment SIGNATURE uu parameter SIGNA containing the signature of the data structure STATE VECTOR array S containing various parameters that are re quired to describe this data structure BUNDLELENGTH parameter L containing the fuel bundle length NAMEGLOBAL jj array G containing the names of the global parame ters PARAMGLOBAL array G containing the value of the global parame ters NAMELOCAL uu array containing the names of the local parame ters CEELI Deiana array Ci j containing an identification number asso ciated with bundle and channel j FUELIDuuuuuu array F containing the fuel type associated with bundle and channel j continued on next page
183. m Default lyld F character T F corner discontinuity factor Default Icdf F character T F Beta of delayed neutron Default Ibet F T F T F T F character T F group wise power form function Default lgff F T F character T F Lambda of delayed neutron Default lamb F T F character T F Decay heat beta and lambda Default Idec F IGE 344 IUPS iups SFAC sfac BFAC bfac 106 keyword used to set the treatment for up scattering 0 1 2 0 keep up scatter XS 1 remove up scatter XS modify down scatter XS with DRAGON spectrum 2 remove up scatter XS modify down scatter with infinite medium spectrum Default iups 0 y gt Dg Ys gg Ys gg T Ys geg X dg for g lt g where g g are the spectra flux either provided by DRAGON or infinite spectra computed in GenPMAXS keyword used to set sfac See Ref 40 the scattering cross section factor If sfac is different from 1 then the scattering cross section will be multiplied by sfac Default sfac 1 0 keyword used to set bfac See Ref 40 the multiplier for betas If bfac is different from 1 then the betas will be multiplied by bfac Default bfac 1 0 IGE 344 107 5 THERMAL HYDRAULICS MODULES 5 1 The THM module The THM module is a simplified thermal hydraulics module where the reactor is represented as a collection of independent channels with no cross flow between them Ea
184. m group gtog be The detector response parameter it is product of cross section and local flux ratio The fission cross section 3 Assembly Discontinuity Factors Table 114 Records in the sub directory MICROLIB_XS Type R Ngr R Ner Condition Units Comment The microscopic absorption cross section of ad al gg Xenon Ca Xe Otot Xe _ Bri Di Xe E The microscopic absorption cross section of A 3 gg Samarium o Olot Sm a Sm s Sm continued on next page IGE 344 Records in the sub directory MICROLIB_XS Type Condition Units XEND uuuuuuva R Ner SMND uuuuuuu R Ner 157 continued from last page Comment The Xenon number density nxe The Samarium number density ngm Each component of the list DIVERS is a directory containing informations for all burnup points of a specific branch recovered from the DIVERS block of the saphyb object Inside each groupwise directory the following records will be found Table 115 Records in the sub directory DIVERS Name Type Condition Units B2unuuuuuuua KEFF uuuuuuva XINF uuuuuuva Comment Critical buckling Effective multiplication factor Infinite multiplication factor IGE 344 158 9 EXAMPLES The following examples of input files represent a typical core modeling using DONJON The main characteristics of a simplified design for the ACR 700 benchmark core model are given below upper ZCU abs
185. mation for a new global or local parameter For more information about the parameter data organization on FMAP data structure see Section 8 1 5 keyword used to specify PNAME character 12 identification name of a given parameter This name is user defined so that it is arbitrary however it must be unique so that it can be used for the search of parameter information and interpolation purpose Moreover it is recommended to use the following pre defined values Boron concentration Averaged fuel temperature Surfacic fuel temperature Averaged coolant temperature Averaged coolant density CANDU only parameters T MODE Averaged moderator temperature D MODE Averaged moderator density keyword used to specify PARKEY character 12 corresponding name of a given parameter as it is recorded in the par ticular multi parameter compo file The PARKEY name of a parameter may not be same as its PNAME and can also differ from one multi compo file to another keyword used to indicate that a given parameter is global which will have a single and constant parameter s value keyword used to indicate that a given parameter is local In this case the total number of recorded parameter s values will be set to Nen x Np keyword used to indicate the input or modification of the actual values for a param eter specified using its PNAME keyword used to indicate that a core average value of a local parameter will be pro vided If the ke
186. mation present in the Multicompo is interpolated as mixture 1 values imix index of the mixture that is to be created in the MICROLIB and MACROLIB FROM keyword used to set the index of the mixture in the MULTICOMPO object imixold index of the mixture that is recovered in the MULTICOMPO object By default imixold 1 USE keyword used to set the index of the mixture in the MULTICOMPO object equal to imix TIMAV BURN keyword used to compute time averaged cross section information This option is avail able only if a MAPFL object is set By default the type of calculation TIMAV BURN or INST BURN is recovered from the MAPFL object INST BURN keyword used to compute cross section information at specific bundle burnups This option is available only if a MAPFL object is set By default the type of calculation TIMAV BURN or INST BURN is recovered from the MAPFL object AVG EX BURN keyword used to compute the derivatives of cross section information relative to the exit burnup of a single combustion zone The derivatives are computed using Eq 3 3 of Ref 16 written as Da _ 1 o dB Y B E SY coe BEL Biz OB Be Bese Bboc i Bee j k where Be B5 and By are the beginning of cycle burnup of bundle j k end of cycle burnup of bundle j j k and exit burnup of channel j This option is available only if a MAPFL object is set By default the type of calculation TIMAV BURN or INST BURN is recovered from the MAPFL object
187. mber of days keyword used to specify that rtime is a number of hours keyword used to specify that rtime is a number of minutes keyword used to specify that rtime is a number of seconds keyword indicating the end of data specific to the actual fuel cycle keyword for obtaining a CLE 2000 variable that is a measure of the discrepancy be tween two cycles character 12 identification name of the first fuel cycle to compare character 12 identification name of the second fuel cycle to compare IGE 344 DIST BURN epsburn DIST POWR epspowr SET PARAM PNAME PNAME pvalue 48 keyword used to recover the discrepancy on burnup distribution in a CLE 2000 vari able character 12 CLE 2000 variable name in which the extracted burnup discrepancy expressed in MW day tonne will be placed keyword used to recover the relative error on power distribution in a CLE 2000 variable character 12 CLE 2000 variable name in which the extracted power relative error will be placed keyword used to indicate the input or modification of the actual values for a param eter specified using its PNAME keyword used to specify PNAME character 12 name of a parameter single real value containing the actual parameter s values Note that this value will not be checked for consistency by the module It is the user responsibility to provide the valid parameter s value which should be consistent with those recorded in the multicompo
188. me a bundle is shifted and stay in the reactor its burnup and power will be saved in the records bshift and pshift For example bshift i and pshift i will contain all the burnups and powers of bundles that have been shifted th time 8 1 3 The FUEL sub directories Each FUEL sub directory contains the information corresponding to a single fuel type Inside each sub directory the following records will be found IGE 344 129 Table 83 Records in FUEL sub directories Type Condition Units Comment MIX uuuuuuua Fuel type mixture number TOT Goo Total number of fuel bundles for this fuel type MIX VOID uuu Voided cell mixture number for this fuel type WEIGHT jou Fuel weight in a bundle for this fuel type ENRICH uuu Fuel enrichment for this fuel type POISON LULU Poison load for this fuel type 8 1 4 The hcycle sub directories Each hcycle sub directory contains the information corresponding to a single PWR fuel cycle Inside each sub directory the following records will be found Table 84 Records in hcycle sub directories Condition Units Comment TIME ou uuu Depletion time corresponding to instanta neous burnup values BURNAVG uuuu Average burnup of the assembly NAME Sou Names of each assembly or of each quart of assembly during a refuelling cycle All quart of assembly belonging to the same assembly have the same name FLMIX uuu I Nen No Fuel mixture indices per assembly subdivi sions for each r
189. module see the user guide Only 3 D Cartesian or 3 D Hexagonal fuel map geometry is allowed keyword used to call the NAP module The heterogeneous assembly geometry defini tion will be called using the geometry defined previously with the embedded module GEO and the COMPO data structure See section 7 1 3 for important note on the coarse geometry requirement keyword to specify that assembly related information are provided number of assemblies number of assemblies along x direction number of assemblies along y direction keyword to specify the assembly number iza of each channels assembly belonging number after A ZONE keyword to automatically compute the assembly number of each chan nel A call to the embedded module NAP is required previously keyword to specify the assembly position names along x direction XNAMEA character 2 array of horizontal channel names A horizontal channel name is iden tified by the channel column using numerical characters 1 2 3 and so on Note that the total number of X names must equal to nxa keyword to specify the assembly position names along y direction YNAMEA character 2 array of horizontal channel names A horizontal channel name is identi fied by the channel column using numerical characters A B C and so on Note that the total number of Y names must equal to nya keyword used to specify XNAME Not used for 3 D Hexagonal geometry charact
190. n only fuel map edition The MCC module specifications are Table 21 Structure MCC FLMAP1 MCC FLMAP1 FLMAP2 descmcec1 where FLMAPI character 12 name of the MAP object that will contain the updated fuel lattice information If FLMAPI appears on both LHS and RHS it will be updated if it only appears on RHS it will only be read to display its contents FLMAP2 character 12 name of the MAP object that contains information to be recovered to update FLMAP1 If FLMAP2 exists data to update FLMAPI will be taken in it If not data to update FLMAPI will be taken in FLMAPI descmcc1 structure describing the main input data to the MCC module Note that this input data is mandatory and must be specified either if FLMAPI is updated or only read 3 8 1 Main input data to the MCC module Table 22 Structure descmcc1 EDIT iprint REC rec1 UNI valuel ADD value2 SAME READ rec2 TTD pcore where IGE 344 EDIT iprint REC recl UNI valuel ADD value2 SAME READ rec2 TTD pcore 32 keyword used to set iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value larger values produce increasing amounts of output A value of 5 or higher will display the contents of FLMAPI object after edition and a value of 6 or higher will display the contents before edition keyword used to indicate tha
191. n gt gt lt lt mZCRout gt gt ENDROD ROD 14 ROD NAME ZCRO7B LEVEL 0 0 AXIS Y FROM H MAXPOS 205 0 227 0 260 0 520 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 15 ROD NAME ZCRO8A LEVEL 0 0 AXIS Y FROM H MAXPOS 249 0 271 0 0 0 260 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 16 ROD NAME ZCRO8B LEVEL 0 0 AXIS Y FROM H MAXPOS 249 0 271 0 260 0 520 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 17 ROD NAME ZCRO9A IGE 344 173 LEVEL 0 0 AXIS Y FROM H MAXPOS 293 0 315 0 0 0 260 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 18 ROD NAME ZCRO9B LEVEL 0 0 AXIS Y FROM H MAXPOS 293 0 315 0 260 0 520 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 19 ROD NAME ZCR10A LEVEL 0 0 AXIS Y FROM H MAXPOS 337 0 359 0 0 0 260 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 20 ROD NAME ZCR10B LEVEL 0 0 AXIS Y FROM H MAXPOS 337 0 359 0 260 0 520 0 272 415 321 945 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 21 ROD NAME ZCR11A LEVEL 0 0 AXIS Y FROM H MAXPOS 161 0 183 0 0 0 260 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt ENDROD ROD 22 ROD NAME ZCR11B LEVEL 0 0 AXIS Y FROM H MAXPOS 161 0 183 0 260 0 520 0 421 005 470 535 DMIX lt lt mZCRin gt gt lt lt mZCRout gt gt
192. nal instan taneous burnup of each assembly subdivision e modify a local parameter such as the Boron concentration in the coolant The SIM module specification is Table 32 Structure SIM FMAP SIM FMAP POWER descsim where FMAP character 12 name of a FMAP object that will be updated by the SIM module The FMAP object must contain the instantaneous burnups for each assembly subdivision a basic naval coordinate assembly layout and the weight of each assembly subdivision POWER character 12 name of a POWER object containing the channel and powers of the assembly subdivisions previously computed by the FLPOW module The channel and powers of the assembly subdivisions are used by the SIM module to compute the new burn up of each assembly subdivision If the powers of the assembly subdivisions are previously specified with the module RESINI you can burn your core without a POWER object descsim structure describing the input data to the SIM module 3 14 1 Input data to the SIM module Note that the input order must be respected Table 33 Structure descsim EDIT iprint CYCLE hcnew FROM hcold BURN indcycle burncycle MAP hx i i 1 Ix hease i j i 1 Ix j 1 ly hy j QMAP hx i i lx 2 1 Ix hy j hease i j i 1x 2 1 Ix j 1y 2 1 ly continued on next page IGE 344 46 Structure descsim continued from last page SPEC asmb1
193. nent of DONJON is dedicated to the interpolation of MACROLIB data from a SAPHYB ob ject the reactor database produced by module SAP in DRAGON or by module SAPHYB in APOLLO2 A set of global parameters are defined for each material mixture and used as multi dimensional interpo lation variables The calling specifications are Table 57 Structure SCR MLIB SCR MLIB MLIB2 SAPNAMI SAPNAM2 MAPFL scr data where MLIB character 12 name of a MICROLIB type L_LIBRARY or MACROLIB type L_MACROLIB containing the interpolated data If this object also appears on the RHS it is open in modification mode and updated A MACROLIB object cannot be specified on the RHS MLIB2 character 12 name of an optional MICROLIB object whose content is copied on MLIB SAPNAMI character 12 name of the LCM object containing the SAPHYB data structure L_SAPHYB signature SAPNAM2 character 12 name of an additional LCM object containing an auxiliary SAPHYB data structure L_SAPHYB signature This object is optional MAPFL character 12 name of the MAP object containing fuel regions description global pa rameter information burnup fuel coolant temperatures coolant density etc Keyword TABLE is expected in scr_data scr_data input data structure containing interpolation information see Section 4 3 1 Note SAPHYB files generated by APOLLO2 don t have a signature If such a SAPHYB is given as input to module SCR a signatu
194. ng DRAGON module CPO data structure containing the multi parameter database generated by the lattice code This object has a signature L_MULTICOMPO it is created using DRAGON module COMPO data structure containing the multi parameter database generated by the lattice code This object has a signature L_SAPHYB it is created using the APOLLO2 lattice code or the DRAGON module SAP data structure containing the fuel lattice specification This object has a signature L_MAP it is created using DONJON module RESINI data structure containing the extended reactor material index This object has a signature L_MATEX it is created using DONJON module USPLIT data structure containing the devices specification This object has a signature LDEVICE it is created using DONJON module DEVINI data structure containing detector positions and responses This object has a signature L_DETECT it is created using DONJON module DETINI and can be modified by the modules DETINI and DETECT data structure containing a tracking information of the reactor geometry This object has a signature L_TRACK it is created using TRIVAC module TRIVAT data structure containing a set of system matrices This object has a signature L_SYSTEM it is created using TRIVAC module TRIVAA data structure containing the numerical solution to an eigenvalue problem This object has a signature L_FLUX it is created using TRIVAC module FLUD data structu
195. ng the name of isotopes used in this fuel type array containing the mixture associated with each isotopes in this fuel type array p containing the density of each iso topes Each cell isotopic sub directory CELL contains the following information Table 99 Cell sub directory Type Condition Units FUELDEN INIT PARAMLOCALBR PARAMLOCALAR Comment array containing the initial density of heavy element in the fuel pp in g cm and the initial linear density of heavy element in the fuel my in g cm array y containing the value of the local pa rameters before refueling array y containing the value of the local pa rameters after refueling continued on next page IGE 344 Cell sub directory PARAMBURNTBR PARAMBURNTAR DEPL PARAM ji ISOTOPESDENS Type Condition Units 8 7 Contents of thm data structure 143 continued from last page Comment array containing the depletion time TP in days and the burnup power rate P in kW kg be fore refueling array containing the depletion time T4 in days and the burnup power rate PA in kW kg after refueling array containing the time step T in days the burnup B in kWd kg and the irradiation w in n kb currently reached by the fuel in this cell array p containing the density of each iso topes This data structure contains the thermal hydraulics information required in a multi physics calculation 8 7 1 The main thm directory The follow
196. nting each detector name 5 for finite element numbers and total number of finite elements for each detector IGE 344 TIME dt REF kc NORM vnorm SIMEX SPLINE PARAB 52 key word used to set dt time step between two calls to the DETECT module key word used to set kc index used to control the type of calculation 0 for reference calculation 1 normal calculation The reference responses are used to obtain detector current responses in full power fractions key word used to set vnorm value used to normalized responses of all the detectors present in DETECT key word used to specify that a polynomial interpolation of detector fluxes according to HQSIMEX method This interpolation will be applied only for vanadium detectors under NAMTYP of value VANAD_REGUL key word to specify that the flux at detector site will be computed with a spline method key word to specify that the flux at detector site will be computed with a parabolic method IGE 344 53 3 17 The CVR module The CVR module is used to update the fuel type index and the coolant densities throughout the reac tor core as required for the voiding simulations A particular core voiding pattern is either selected from the several pre defined patterns or directly defined by the user in an arbitrary fashion In the last case the user may specify the individual voided channels by indicating their identification names The CVR module will create a new p
197. o control the axial shape conver gence over the external time average iterations The optimal value which corresponds to the minimal total number of such iterations can be found by performing several runs at different relval The default value of the relaxation parameter is set to 0 5 B EXIT keyword used to indicate the calculation of the core average exit burnup and the chan nel refuelling rates 3 12 2 Time average calculation using DONJON When the average exit burnups are provided for each channel the exact burnup integration limits for each fuel bundle are unknown and need to be determined The burnups integration limits are function of the normalized axial power shape which in turn depends on the flux solution over the fuel lattice Moreover the flux solution depends on the fuel map macrolib i e fuel properties which in turn depends on the burnups integration limits for each fuel bundle Consequently the time average calculation is an iterative procedure that consists to repeat all the steps required for the axial power shape computation This repetition is to be made until the relative error between the two successives axial power shape calculations becomes as small as required for the precision The axial power shape computing scheme is composed of several steps each step is performed using an appropriate DONJON or TRIVAC module 1 An initial axial power shape is set as a flat distribution over the fuel lattice and the first
198. obers parked above core ZCU 1 2 amp 3 ZCU Upper 7 889 upper lower lower ZCU absobers parked below core Figure 12 Face View of ACR Benchmark Core Model 292 Channels e Number of reactor channels 292 e Number of fuel bundles per channel 12 e Core radius 260 cm e Core length 594 36 cm e Lattice pitch 22 cm e Reactor thermal power 1800 MW th IGE 344 9 1 Examplel Compo related example Input data for test case Examplel x2m DR KK K K FK FKK K K AE AE AAA A FK FK ARE ARE AE AE AAA AA AAA K FK K FK FK K FK FK FK K K K K K K K K OK Purpose Author s x XA Input file Examplei x2m Test for non regression using DONJON 4 D Sekki 2007 11 x XA XX RE RE RE KK K KK K KK K FK K OK K dd K K ok OK PROCEDURE MODULE LINKED_LIST LINKED_LIST kh variables kh assertS Pgeom Pfmap Pburn Pdevc DELETE GREP END CRE MACINI FLUD USPLIT DSET TAVG FLPOW TRIVAT TRIVAA NEWMAC GEOM TRACK MATEX FMAP FLUX POWER MACFL DEVICE MACRO1 MACRO2 MACRO SYSTEM LREFL1 LREFL2 LFUEL1 LFUEL2 LZCRin LZCRot LSORin LSORot INTEGER nbMix 8 INTEGER nbRefl nbFuel INTEGER mFueli mFuel2 INTEGER mRefli mRef12 INTEGER mZCRin mZCRot INTEGER mSORin mSORot 200 N On ANN 100000 INTEGER MaxReg STRING Method MCFD INTEGER degree quadr 11 INTEGER iter iEdit 0 5 REAL Power 1800 REAL epsil 1 E 5 REAL Precf 1
199. olation mode is set in Sect 4 3 1 character 12 user defined keyword associated to a global parameter to be set value of a global parameter used to interpolate vall is the initial value of this param eter in case an average is required vall can be an integer real or string value value of the final global parameter By default a simple interpolation is performed so that val2 vall val2 is always a real value with val2 gt vall keyword used to indicate that the value of parameter vall or the second value for the Ao calculation is recovered from MAPFL i e the MAP object containing fuel regions description keyword only available together with the ADD option It is used to set all the other variable values when a A contribution is performed for one variable value of the reference parameter when it is directly given by the user Note that there is no default value keyword used to specify that the reference value will be the same as in the refence case i e for the Cret computation keyword only available together with the ADD option It is used to specify that all the other variable values which are required are given keyword used to set the number densities of some isotopes present in the SAPHYB object The data statement MICRO ALL is used by default keyword to indicate that all the isotopes present in the SAPHYB object will be used in the MICROLIB and MACROLIB objects Concentrations of these isotopes will be recove
200. on corresponding to the pkey i In the case where pkey i is TCOM or TMOD the temperature must be in Celsius keyword used to call a default meshing the values for BARR and BURN are extracted from Saphyb four default values are considered for DMOD CBOR TCOM and TMOD if exists These values correspond to the first mid and last values of the initial SAP or MCO meshing This otpion is used if the number of banches in the Saphyb or defined by the user exceeds 1000 keyword used to set the type of Assembly Discontonuity Factor to be recovered from the SAP or MCO object NB the ladf flag must be set to true Assembly discontinuity factors are generated using the DRAGON V5 procedure De fault option NB This option is available with Saphyb produced by the DRAGON V5 lattice code using a 2 level flux calculation with the Method Of Characteristics a ADE Hom g where g is the energy group f the assembly surface and on surfacic homogene and heterogene fluxes in asssembly et item are the average Assembly discontinuity factors are generated using the Generalized Equivalence The ory NB This option is available with Saphyb produced by the APOLLO2 lattice code using a 2 level flux calculation with the Method of Characteristics Het g 17Ne txh ADF Hom dg where g is the energy group f the assembly surface di ia ani are the homogene and heterogene fluxes in asssembly D the diffusion coefficien
201. onds to the pin by pin geometry Note that another multicompo with all pin wise properties is needed to be able to use the automatically generated geometry keyword to specify mixtures corresponding to assembly in the coarse geometry number of type of assembly mixture number of the assemblies keyword to specify the mesh splitting at assembly level split along x direction for each mesh of the heterogeneous assembly number of mesh of the heterogeneous assembly along x direction keyword to specify the mesh splitting at assembly level split along y direction for each mesh of the heterogeneous assembly number of mesh of the heterogeneous assembly along y direction IGE 344 124 MAX MIX GEO keyword to specify the number of mixtures in the original core geometry i e before the core geometry is splited by the NAP module This keyword is mandatory if there is a reflector in the geometry otherwise the numbers for fuel mixtures will not match between the split core geometry GEOMETRY and the split fuel geometry GEOMETRY embedded in MAP nmxgeo number of mixtures in the original core geometry Note The included geometry in the COMPO has to be unfolded even if the transport calculations are done on a 1 8th assembly Moreover no split can be defined in the geometry one mesh ONLY per heterogeneous mixture is mandatory IGE 344 125 8 DONJON DATA STRUCTURES A brief description of each DRAGON DONJON and TRIVAC data structures wh
202. only with the mono parameter COMPO objects and the nuclear properties can be interpolated only with respect to the burnup data In case of the MACROLIB con struction from a multi parameter database the NCR module should be used instead In this case the interpolation of nuclear properties can be made with respect to global and local parameters if they were previously specified in the fuel map see Section 3 1 2 The CRE module specifications are Table 47 Structure CRE MACRO CRE MACRO CPO desccre1 MACEL CRE CPO FMAP desccre2 where MACRO character 12 name of the MACROLIB object to be created or updated for the few reactor material properties Note that if MACRO appears on the RHS the information previously stored in MACRO is kept CPO character 12 name of the COMPO object containing the mono parameter database from transport calculations MACFL character 12 name of the fuel map MACROLIB that will be created only for the fuel properties over the fuel lattice FMAP character 12 name of the FMAP object containing the fuel map specification and burnup informations desccrel structure describing the input data to the CRE module when the FMAP object is not specified IGE 344 desccre2 64 structure describing the input data to the CRE module for the fuel map MACROLIB construction 4 1 1 Input data for the CRE module EDIT iprint NMIX nmix Table 48 Structure
203. operties are essentially recovered from the reactor database obtained from the lattice calculations using DRAGON code The two distinct macroscopic cross section libraries can be constructed using DONJON The first MACROLIB is constructed only for the material properties which are evolution independent such as reflector and devices properties The second MACROLIB is constructed only for the fuel properties defined per each fuel bundle over the fuel lattice The two libraries are next combined and updated according to the devices insertion level The produced extended MACROLIB is subsequently used to obtain the numerical solution using TRIVAC modules Finally it should be noted that the DONJON code development is permanently in progress The IGE 344 2 future updates will provide several extended capabilities for the reactor design and calculations they will be gradually added to the subseguent DONJON versions IGE 344 3 2 GENERAL SPECIFICATION OF DONJON 2 1 Modules Reactor calculations using DONJON are performed by means of seguential execution of several user selected modules according to the user defined computing scheme Each module is designed to perform some particular tasks The detailed description of DONJON modules is given in Section 3 to Section 6 In order to perform the reactor calculations it is also reguired to use some DRAGON and TRIVAC modules For more details on DRAGON modules specification refer to its user guidel
204. or Saphyb database IGE 344 49 3 15 The XENON module The XENON module is used to correct the Xenon distribution coming from an interpolation calculation This module computes the new densities according to the bundle flux and the equation providing the balance concentration of Xenon 135 Yr Yx 210 Nx _ 3 1 ds Ax 0x9 ey where e Yr is the fission yield of 1135 e Yx is the fission yield of Xe135 e ox is the capture cross section of Xel35 e Ax is the decay constant of Xe135 e V is the total fission cross section e is the bundle flux The XENON module specification is Table 34 Structure XENON MICROLIB XENON MICROLIB POWER descxenon where MICROLIB character 12 name of a LIBRARY object that will be updated by the XENON module The Xenon should be extracted in this library for the use of this module POWER character 12 name of a POWER object containing the bundle fluxes previously com puted by the FLPOW module The fluxes should be normalized to the reactor power descxenon structure describing the input data to the XENON module 3 15 1 Input data to the XENON module Note that the input order must be respected Table 35 Structure descxenon EDIT iprint continued on next page IGE 344 50 Structure descxenon continued from last page INIT where EDIT keyword used to set iprint iprint integer index used to control the printing on screen INIT keywor
205. ories lt o ee 135 8 3 5 The DEV LZC sub directories LL 135 8 3 6 The LZC GROUP sub directories 136 IGE 344 iv 8 4 Contents of a detect data structure a 137 8 4 1 The state vector content lt o oo o es 137 8 4 2 The main detect directory o oooooooo oo m o 137 8 4 3 The name_type sub directories lt a 137 8 4 4 The name detect sub directories 20 0 20 0 00000 138 8 5 Contents of power data structure o oo ee 138 8 5 1 The state vector content 138 8 5 2 The pow t directory oo ia ea cesera a 139 8 6 Contents of history data structure aooaa a a 140 8 6 1 The main directory cor a die he A We e 140 8 6 2 The fuel type sub directory lt 142 8 6 3 The cell type sub directory o eee eevee 142 87 Contents of thm data structure il lle a eh wa we 143 8 7 1 The main thm directory lt o ec 6 eee de es 143 2 The HISTORY DATA sub directory lt 147 8 8 Contents of a optimize data structure LL 148 8 8 1 The sub directory OLD VALUE in optimize 150 8 9 Contents of d2p_info data structure lt 151 5 9 1 Whe st te vector CONES o us aa e AR E DE 151 9 EXAMPLES lt lt sac ie eke DA wae e g a AS 158 9 1 Examplel Compo related example 0 0000 eee 159 9 2 Example2 Multicompo related example o o o o 163 9 3 o AI Ra eRe eee dae PERE
206. ory TH DATA Condition Units Comment The steady state energy spectrum of the neu tron emitted by fission x The average of the inverse neutron velocity Decay constant of delayed neutron Effective delayed neutron fraction Effective yield value of Pm 149 Effective yield value of Xe 135 Effective yield value of 1 135 IGE 344 155 Table 111 Records and sub directories in sub directory BRANCH_INFO Condition Units Comment BRANCH_IT uu Index of the current branch type calculation BRANCH uuu Current BRANCH type calculation BRANCH NB uu Number of branches to be created in PMAXS files REF_STATE uu Values of state variables for reference state burnup excepted HST_STATE uu Values of state variables for history state bur nup excepted STATE_INDEX Current index value for each state variables BRANCH_INDEX Index of current branch calculation STATE uuuuuu State variable values for the current branch calculation REWIND uuuuu Indicates the end of a branch calculation REWIND 1 STOP iii Indicates the end of calculation STOP 1 CROSS SECT Dir Nbu Each component of the list is a directory con taining cross sections for all burnup point of a specific branch cross sections are overwrited after each branch calcualtion DIVERS uuu Dir Vpu Igria lt 1 Each component of the list is a directory con taining the DIVERS data block recovered from the saphyb object inforamtions are overwrited after each bran
207. oses The TAVG module specification is Table 28 Structure TAVG TAVG FMAP POWER desctavg where FMAP character 12 name of a FMAP object that will be updated by the TAVG module The FMAP object must contain the average exit burnups and refuelling schemes of channels POWER character 12 name of a POWER object containing the channel and bundle powers previously computed by the FLPOW module The channel and bundle powers are used by the TAVG module to compute the normalized axial power shape over each channel desctavg structure describing the input data to the TAVG module 3 12 1 Input data to the TAVG module Note that the input order must be respected Table 29 Structure desctavg EDIT iprint AX SHAPE RELAX relval B EXIT IGE 344 40 where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value 2 only the burnup limits over each channel are printed 3 only the axial power shape values over each channel are printed 4 only the channel refuelling rates are printed for larger values of iprint everything will be printed AX SHAPE keyword used to indicate the calculation of the new axial power shape and correspond ing burnups limits over each reactor channel RELAX keyword used to set the relaxation parameter relval relval real value of the relaxation parameter generally used t
208. osition of an assembly as defined by user in RESINI module Renumbered mixture indices per each fuel re gion over the fuel map geometry for the non fuel regions these indices are set to 0 Channel identification names with respect to their horizontal position Channel identification names with respect to their vertical position Name of the assembly on X direction 4 char acter name per assembly continued on next page IGE 344 Records and sub directories in fmap data structure AYNAME uuu ASSEMBLY uuu B ZONE uuuuu BURN AVG uuu BURN INST uu BURN BEG uuu BURN END uuu BUND PWuuuuu BUND PW INI FLUX AV uuu B EXITuuuuvu REF SHIFT uuu REF VECTOR yu REF SCHEME REF RATE is REF CHAN uu DEPL TIME uu pshift bshift ishift AX SHAPE uuu EPS AXuuuuuu GEOMAP uuu z SFE B a 2 ui E B D da Dea Z Z n No Nor R 1 I Neomb I Neomb No I Nen R Nen R 1 R Na Np R Nen Np I Ncn No R Nen Np R 1 Dir Condition Nsht gt I Nsht 2 1 Nsht gt 1 Ibtyp 1 Typ 1 Units 127 continued from last page Comment Name of the assembly on Y direction 4 char acter name per assembly List of assembly directories Each component of this list follows the specification presented in Section 8 1 6 Combustion zone indices per channel MW dt Average exit burnups per combustion zone MW d t Instantaneous burnups per bundle or assem
209. p2 i ss 44 Gases BAe EI ae ee ee he es PURE 122 Structure desenap3 2 4 66 44 eR e 123 Records and sub directories in fmap data structure 2 2 2 o o o 126 Records in FUEL sub directories s ew 129 Records in hcycle sub directories 2 ee 129 Records in PARAM sub directories k k ws 130 Assembly type sub directory 2 soa e e ee 130 Records in matex data structure lt ks 132 Records and sub directories in device data structure 22 2 133 Records in DEV ROD sub directories i ee aa 134 Records in ROD GROUP sub directories 2 o e k aeea 135 Records in DEV LZC sub directories o e ee 136 Records in LZC GROUP sub directories lt a be ss 136 Records and sub directories in detect data structure lt o e mo 137 Records in name_type sub directories ks 138 Records in name_detect sub directories ks 138 Records in power data structure o 139 Main records and sub directories in history o ooo 140 Pusl EDSON cosa ee RA EB nw eee a AA AR 142 Cell BuB dIEEEEGES o eo poco le he AA a ee ee a Pek bk eh eee 142 Main records and sub directories in thm ooo o e 143 Sub directories in HISTORY DATA directory e eee 147 Records in each CHANNEL directory ee E a 147 Main records and sub directories in optimize 2 0 2 ooo 149 IGE 344 vii 104 1
210. pe given in 1074 particles per cm keyword used to indicate that the number density for the isotope HISO will be recov ered from the COMPO object keyword used to indicate that all the number densities are to be recovered from the COMPO object keyword used to indicate the end of data specification for the material mixture mix IGE 344 67 4 2 The NCR module This component of DONJON is dedicated to the interpolation of MICROLIB and MACROLIB data from a MULTICOMPO object the reactor database produced by COMPO A set of global and or local parameters are defined for each material mixture and used as multi dimensional interpolation variables The calling specifications are Table 52 Structure NCR MLIB NCR MLIB MLIB2 CPONAMI CPONAM2 MAPFL ncr_data where MLIB character 12 name of a MICROLIB type L_LIBRARY or MACROLIB type L_MACROLIB containing the interpolated data If this object also appears on the RHS it is open in modification mode and updated A MACROLIB object cannot be specified on the RHS MLIB2 character 12 name of an optional MICROLIB object whose content is copied on MLIB CPONAMI character 12 name of the LCM object containing the MULTICOMPO data structure L_MULTICOMPO signature CPONAM2 character 12 name of an additional LCM object containing an auxiliary MULTICOMPO data structure L_MULTICOMPO signature This object is optional MAPFL character 12 name of the MAP object
211. pe 1 half part 2 quarter 3 eighth The number of sides in assembly Nsurf S11 The number of corners in assembly Neorn 612 The number of assembly discontinuity factors per energy groups Naaf S13 IGE 344 152 The following records and sub directories will be found on the first level of d2p_info directory Table 105 Records and sub directories in d2p info data struc ture Condition Units Comment SIGNATURE uu Cx12 Signature of the d2p info data structure SIGNA L_INFO uuuuu STATE VECTOR 1 40 Vector describing the various parameters as sociated with this data structure BARR_ INFO uu I Ncomp Meaning of control rod in the saphyb object Neomp is the number of composition for con trol rods including the unrodded case SAPHYB_INFO Dir Information related to the saphyb content HELIOS_HEAD Information related to the header of output DRA file GENPMAXS_INP i Information related to the GenPMAXS input file BRANCH_INFO i Information related to the current branch cal culation TH_DATALLULU Information related to the invariant T H data block Each componant of the list is a direc tory containing TH data for a single burnup point Table 106 Records and sub directories in SAPHYB_INFO Condition Units Comment PKEY_INFOLuL Each component of the list is a directory con taining information for all possible state vari ables STATE VAR uuu Name of state variables PKEY ISOTOPES uuu Reference
212. pin power recon struction number used to create the named of the flux record representing y for more details d m ae See Sect 7 1 1 keyword used to normalize the flux If this keyword is not used the flux directly computed by the FLUD is used to perform the pin power reconstruction then normal ization has to be performed by the user indepently power used to normalize the flux MW IGE 344 123 7 1 3 Heterogeneous assembly geometry definition EDIT iprint Table 81 Structure descnap3 DIRGEO namedir MACGEO MIXASS nmix imix i i 1 nmix SLPITX ASS ispx i i 1 nxass SLPITY ASS ispy i i 1 nyass MAX MIX GEO nmxgeo 3 where EDIT iprint DIRGEO namedir MACGEO MIXASS nmix imix SLPITX ASS ispx nxass SLPITY ASS ispy nyass keyword used to modify the print level iprint index used to control the printing of this module The iprint parameter is important for adjusting the amount of data that is printed by this calculation step e iprint 0 results in no output e iprint 1 keyword to specify namedir name of the directory containing the heterogeneously or homogeneously homogenized cross sections keyword to specify that the macro geometry stored in the GFF record of the macrolib will be used instead of the macro geometry of the heterogeneously or homogeneously homogenized cross sections Usually this other macro geometry of the assembly cor resp
213. quadratic constraint of the form nvar 2 2 does aaf lt 5 i 1 where w is a weight defined after keyword CST WEIGHT and Ag is a displacement for i th control variable at iteration n The initial value of radius SU is defined after keyword OUT STEP LIM reduces the radius SV of the quadratic constraint The calling specifications are Table 74 Structure GRAD OPTIM GRAD OPTIM DFLUX GPT grad_data where OPTIM character 12 name of the OPTIMIZE object L_OPTIMIZE signature containing the optimization informations Object OPTIM must appear on the RHS to be able to updated the previous values DFLUX character 12 name of the FLUX object L_FLUX signature containing a set of solutions of fixed source eigenvalue problems GPT character 12 name of the GPT object L_GPT signature containing a set of direct or adjoint sources grad_data structure containing the data to the module GRAD see Sect 6 2 1 6 2 1 Data input for module GRAD Table 75 Structure grad_data EDIT iprint continued on next page IGE 344 Structure grad data where EDIT iprint METHOD SIMPLEX LEMKE MAP AUG LAGRANG PENAL METH OUT STEP LIM ST 115 continued from last page METHOD SIMPLEX LEMKE MAP AUG LAGRANG PENAL METH OUT STEP LIM sr OUT STEP EPS 64 INN STEP EPS inn CST QUAD EPS quad MAXIMIZE MINIMIZE STEP REDUCT HALF PARABOLIC VAR VALU
214. r keyword used to specify the controller filling axis A liquid controller can be filled along only one vertical axis keyword used to specify that a liquid controller is filled along X axis keyword used to specify that a liquid controller is filled along Y axis keyword used to specify that a liquid controller is filled along Z axis keyword used to specify the actual filling level real positive value of the water level This value is minimal value 0 0 when the controller is empty and it is maximal value 1 0 when the controller is full filled For the partially filled controller the water level must be 0 0 lt value lt 1 0 keyword used to specify rate real positive value of the water filling rate given in m s This value is needed only for the reactor regulating purpose keyword used to specify time real value of the filling time given in sec This value is needed only for the reactor regulating purpose keyword used to specify mixE two integer mixture indices specified for the empty part of liquid controller The first and the second mixture indices correspond to the perturbed and the reference cross sections respectively These indices will be used to compute the incremental cross sections in the NEWMAC module keyword used to specify mixF two integer mixture indices specified for the full part of liquid controller The first and the second mixture indices correspond to the perturbed and the reference
215. r Ti represents e The number of detectors of type name_type Zi e The number of delayed responses 2 Za 8 4 2 The main detect directory The following records and sub directories will be found on the first level of detect directory Table 93 Records and sub directories in detect data structure Type Condition Units Comment SIGNATURE uu Signature of the detect data structure SIGNA L_DETECT 111 STATE VECTOR Vector describing the various parameters as sociated with this data structure S name_type i Detector type sub directory contains informa tions for each detector of this type 8 4 3 The name_type sub directories Inside each name_type sub directory the following records will be found IGE 344 138 Table 94 Records in name_type sub directories Type Condition Units Comment INFORMATIONS Record containing describing the various pa rameters associated with a detector type Z INV CONST uuu The inverse time constant of the delayed de tector responses FRACTION uuu The delayed and prompt fractions of the de tector responses SPECTRAL uuu The energy spectrum of the detector name_detect i Detector information sub directory 8 4 4 The name_detect sub directories Inside each name_detect sub directory the following records will be found Table 95 Records in name_detect sub directories Type Condition Units Comment NHEX uuu H nhex her 1 The numbers of affecte
216. r reflector temperature gt RD for reflector density 11 RP for reflector purity Note that these keywords are identical to those used in the Proc16 program Here the moderator coolant and reflector can be D20 H20 or any other mixture since DONJON is not aware of the compositions of these mixtures In the case where many different Tape16 files contains the reference and the individual perturbation effects one must first define the reference case before updating the CPO using the Tape16 files containing the perturbations reference value of the associated local parameter IGE 344 96 npert number of local parameter perturbations valper perturbed values of the local parameter MTMD character 4 name of perturbation associated with combine temperature and density changes effects Note that this keyword is equivalent to the MTS keyword used in the Proc16 program In principle any combined perturbations effects could be built from the catenation of two individual perturbations given in NAMPER valreft reference temperature This is required if either the MT or the MD perturbation is not defined valrefd reference density This is required if either the MT or the MD perturbation is not defined npert number of simultaneous perturbations in moderator temperature and density valpert perturbed values of the moderator temperature valperd perturbed value of the moderator density The explicit name of the mixtures MIXDIR that will be
217. rameters in case of a transient calculation By default a steady state calculation is performed caltype integer value set to control the type of calculation that will be performed by the THM module 0 for steady state 1 for transient The default value is 0 timestep real value set to the time step in case of a transient calculation The default value is 0 0 timeiter integer value of the current time step index used for transient calculations The default value is 0 time real value of time in second used for transient calculations The default value is 0 0 FPUISS keyword used to specify the fraction of the power released in fuel The remaining fraction is assumed to be released in coolant The default value is 0 974 fract real value set to the fraction f Power densities released in coolant and fuel are computed as Veool Viuel Paresh cool 1 a i 1 Veool Vinesh Veool Viuel Pesh Veuel Vinesh Qiuel f IGE 344 CRITFL cflux CWSECT sect flow SPEED velocity ASSMB Sass nbf nbg INLET poutlet tinlet RADIUS rl r2 r3 r4 POROS poros PUFR pufr 109 where Veoo and Vruel are coolant and fuel area computed from sass nbf nbg r3 and r4 The mesh power Pmesh and volume Vmesh are recovered from MAPFL object keyword used to specify the critical heat flux real value set to the critical heat flux in W m The default value is 2 0 x 10 W m keyword used to specify
218. re MOC bce a ea da Sao GS bee ee ea wee de 31 Structure desemicell o cui a a e E eG E 31 Structure MOVDEVS o cade eho rasa AA A 33 Structure desemove k k k a ee 33 Structure NEWMACS ciar oe Ba A A Be A ADR 35 AICA L a G RA R i ela 36 Structure descABOW i pe ee eR ee 37 DEF TAVOS Las is AM LTL ek ee RR a ee 39 Structure desetavg i E ee eee ee 39 atico TINGIS ss a bia Lilia A Ke ee ee Rebate LESAGE 42 Structure desetinst 4 4 4004 8 008 BN eee dee ee ee e 43 SEUCIUTE SIM osorno ER 45 Structure desesimi 0 2 ee GERD ESR e ee ee e 45 Structure AE ONS i a a moa es OG AA a eee De BSR ie 49 Structure desexenom e ss coca k oe a E DA 49 Birches DETECTS o REA ai GE Gwe as Be e EO 51 Structure desedete6t pk kn ha eRe Eau a 51 e CIRE i hk ie se E A A a 53 Structure deserevr kia k heh SD m REL EE LESSEE 53 Updating an HISTORY structure using a MAP structure Ls 56 Updating an HISTORY structure using a BURNUP structure 004 56 Updating a BURNUP structure using an HISTORY structure 56 Updating an HISTORY structure using a MAP structure ks 56 Structures ASA se sa bu wewa da A RI 57 Sfrtichurs OSTBER o o o soma ka a O k R ee ee he a es de ai 58 Structure hstpar lt u ea baad a k bb ALG eee Bai eee ad eb ee 59 Structure CRET c esei RAGE OA a a ee eR Raa ae 63 Structure AEseCREL ini sk See EERE Mee ee ALLE 64 Structure desecr 2 os a eaae e Re ee eS 64 Structure desedatal
219. re containing the powers and normalized fluxes over the reactor core This object has a signature L_POWER it is created using DONJON module FLPOW This data structure contains the information required to ensure a smooth coupling of DRAGON with DONJON when an history based full reactor calculation is to be performed It is used only by the HST module IGE 344 6 2 3 Syntactic rules for input specification The input data to any module is read in free format using the subroutine REDGET CLE 2000 variables 2 are also allowed The user guide for DONJON is written using the following convention e the parameters surrounded by single square brackets denote an optional input e the parameters surrounded by double square brackets denote an input which may be repeated as many times as needed e the parameters in braces separated by vertical bars V denote a choice where one and only one input is mandatory e the parameters in bold face and in brackets denote an input structure e the parameters in italics and in brackets with an index data i i 1 n denote a set of n inputs e the words using the typewriter font KEYWORD are character constants used as keywords e the words in italics denote the user defined variables they are lower case and of integer type when starting from to n or of real type when starting from a to h or from o to z or they are upper case and of character type CHARA
220. re must be included before using it The following instruction can do the job Saphyb UTL Saphyb CREA SIGNATURE 3 L SA PHYB 4 3 1 Interpolation data input for module SCR Table 58 Structure scr_data EDIT iprint MEMORY RES MACRO MICRO LINEAR CUBIC LEAK b2 EQUI TEXTS0 NMIX nmixt SAPHYB SAPNAM descints TABLE SAPNAM namburn descints descdepl IGE 344 where EDIT iprint MEMORY RES MACRO MICRO LINEAR CUBIC LEAK b2 EQUI TEXT80 NMIX nmixt SAPHYB TABLE SAPNAM namburn descints 79 keyword used to set iprint index used to control the printing in module SCR 0 for no print 1 for minimum printing default value keyword activating a copy of the Saphyb into memory before performing interpolation In some cases this operation may reduce CPU resources in SCR keyword indicating that the interpolation is done only for the microscopic cross sections and not for the isotopic densities In this case a RHS MICROLIB must be defined and the number densities are recovered from it This option is useful for micro depletion applications Important note It is possible to force interpolation of some isotopic densities with RES option if these isotopes are explicitely specified with a flag after MICRO keyword in descints input data structure see Section 4 3 2 keyword indicating t
221. real value of the relative coolant density with respect to the nominal or unperturbed conditions associated with the voided reactor channels In general this value equals to 0 0 for the complete voiding of a channel and to 1 0 for an unperturbed channel Intermediate values of dcoolV will then correspond to the partially voided channels It is supposed that all voided channels will have the same dcoolV value keyword used to specify the core voiding pattern which will be used for a particular voiding simulation keyword used to specify the full core voiding pattern According to this pattern the fuel mixtures will be modified for all reactor channels keyword used to specify the half core voiding pattern According to this pattern the fuel mixtures will be modified only for the upper half of reactor channels keyword used to specify the quarter core voiding pattern According to this pattern the fuel mixtures will be modified only for the upper left quarter of reactor channels keyword used to specify the checkerboard full voiding pattern According to this pat tern the fuel mixtures will be modified for all reactor channels in which the direction of coolant flow is positive keyword used to specify the checkerboard half voiding pattern According to this pattern the fuel mixtures will be modified only for the upper half of reactor channels in which the direction of coolant flow is positive keyword used to specify the checkerboard q
222. red from the SAPHYB object or set using the HISO conc data statement keyword to indicate that only the isotopes set using the HISO conc data statement will be used in the MICROLIB and MACROLIB objects IGE 344 82 HISO character 8 name of an isotope conc user defined value of the number density in 1024 particles per cm of the isotope the value of the number density for isotope HISO is recovered from the SAPHYB object ENDMIX end of specification keyword for the material mixture 4 3 3 Depletion data structure Part of the depletion data used in the isotopic depletion calculation the fission yields and the ra dioactive decay constants is recovered from the Saphyb file Remaining depletion data is recovered from the input data structure descdepl This data describes the heredity of the radioactive decay and the neutron activation chain Table 60 Structure descdepl CHAIN NAMDPL izae reaction energy STABLE FROM DECAY reaction yield NAMPAR J ENDCHAIN with CHAIN keyword to specify the beginning of the depletion chain NAMDPL character 12 name of an isotope or isomer of the depletion chain that appears as a particularized isotope of the Saphyb izae optional six digit integer representing the isotope The first two digits represent the atomic number of the isotope the next three indicate its mass number and the last digit indicates the excitation level of t
223. rmally presented in the CANDU6 type reactor core The liquid zone controllers specifications are read from the input data file Note that this modeling can be made after the rod type devices have been previously defined using the DEVINI module see Section 3 4 In this case the previously created DEVICE object will be updated by the LZC module it will store the additional and separate information with respect to the liquid controllers see Section 8 3 The liquid zone controller specification includes several device parameters such as the whole de vice position water filling level direction of filling etc Note that a liquid zone controller is normally composed of two parts one part is empty and the second part is full filled The water level can be adjusted according to the control reactivity requirements The controllers positions are referred using 3 D Cartesian coordinates Several devices parameters can be modified using the DSET module see Section 3 7 The liquid controllers can also be divided into the several user defined groups so that they can be manipulated simultaneously The LZC module specification is Table 15 Structure LZC DEVICE MATEX LZC DEVICE MATEX desclzc where DEVICE character 12 name of the DEVICE object Note if the rod type devices are not present in the reactor core then DEVICE object must appear only on the LHS i e in create mode it will contain the information only with respect
224. ry are allowed e The total number of mesh splitted volumes Nr S7 e The number of mesh splitted volumes along x axis Ly Sg e The number of mesh splitted volumes along y axis Ly Sy e The number of mesh splitted volumes along z axis L S10 8 2 2 The matex directory The following records will be found on the matex directory IGE 344 SIGNATURE uu STATE VECTOR RMIX uuuuuuva RTOT Goo uu PMT X Guu FTOTLUULLLLL MAT uuuuuuuua INDEX yuuuuuu MESHX youuu MESHY uuuuua MESHZ y uuuuua H FACTOR uuu 132 Table 87 Records in matex data structure Condition Units Comment Signature of the matex data structure SIGNA L_MATEX uuu Vector describing the various parameters as sociated with this data structure S The reflector type mixture indices as defined in the reactor geometry The total number of reflector regions per each reflector type The fuel type mixture indices as defined in the reactor geometry The total number of fuel regions per each fuel type The material mixture indices per each region and including the device mixtures The fuel type indices are set negative the device in dices are appended at the end of vector the virtual region indices are set to 0 The renumbered mixture indices A unique number is associated with each mesh splitted volume The device indices are not included the virtual region indices are set to 0 The mesh splitted coordinates a
225. s and or several groups of rods simultaneously A user must be aware that a particular device will not be displaced more than once during the same time step Note that the input order of data to the module must be respected Table 24 Structure descmove EDIT iprint DELT delt ROD id GROUP grp INSR EXTR LEVEL value DELH delh SPEED speed 2 where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value larger values produce increasing amounts of output IGE 344 DELT delt ROD id GROUP igrp INSR EXTR LEVEL value DELH delh SPEED speed 34 keyword used to set delt real value of the time increment for the current time step given in sec keyword used to specify the rod id number integer identification number of a rod type device to be displaced Each rod has a unique id number ranging from 1 to nrod as been defined in the DEVINI module see Section 3 4 2 keyword used to specify a rod group igrp number integer number of a group of rods that will be displaced simultaneously with the same parameters of movement Each group of rod devices has a unique igrp number ranging from 1 to ngrp as been defined in the DEVINI module see Section 3 4 3 keyword used to specify that a particular rod or a group of rods will be inserted into the reactor core during the p
226. s file is controlled by the user who can activate specific keywords in the WIMS AECL input file The standard CPO data structure used by the code DONJON is generally generated by the cell code DRAGON This data structure can be stored on a FORTRAN direct access binary file in the form of a hierarchical data base There is also the possibility to keep the contents of this data structure in memory with the same hierarchical structure for faster access The structure of the data base is in the form of a list of material directories which contain burnup sub directories Inside each of these burnup sub directories the isotopic contents of a mixture is described and the multigroup cross sections associated with a specific isotope are stored in individual sub directories Note that in this database the macroscopic cross sections associated with a mixture are stored in a default isotopic sub directory The interface between the Tape16 file and the CPO data structure should be written as a new module of the code DONJON in order to facilitate the access to the GANLIB utilities which manage the hierarchical data structures This module will be called T16CPO The transfer of information from a Tape16 format file to a CPO data structure will reguire the following DONJON instructions The T16CPO module specifications for creating or updating a CPO data structure from a Tape16 file are Table 66 Structure T16CPO DONCPO T16CPO DONCPO WIMS16 desct1
227. sions for the HISTORY data structure CELLID keyword to identify the cell for which history information is to be processed icha channel number for which history information is to be processed ibun bundle number for which history information is to be processed idfuel fuel type number associated with this cell One can associate to each fuel cell a different fuel type By default a single fuel type is defined and it fills every fuel cell Only the initial properties of each fuel type are saved These properties are used for refueling GET keyword to specify that the values of the parameters selected in brnpar will be read from the input stream or CLE 2000 local variables and stored on the HISTORY data structure PUT keyword to specify that the values of the parameters selected in brnpar will be read from the HISTORY data structure and transferred to local CLE 2000 variables BREFL to specify that the information to extract from the HISTORY data structure is related to the properties of the cell before refueling takes place AREFL to specify that the information to extract from the HISTORY data base is related to the properties of the cell after refueling took place hstbrn structure containing the burnup options hstpar structure containing the local parameters options The hstdim input structure is required for general dimensioning purpose It is generally used only when creating the HISTORY data structure However the numb
228. such as a rod position rod inser tion level direction of movement etc The devices positions can not overlap in the reactor core they are referred using 3 D Cartesian coordinates The insertion level of rods can be set according to their nominal positions or they can be displaced in or out of core The rods can also be divided into the several user defined groups so that they can be manipulated displaced or moved simultaneously The DEVINI module specification is Table 8 Structure DEVINI DEVICE MATEX DEVINI MATEX descdev where DEVICE character 12 name of the DEVICE object that will be created by the module it will contain the devices information MATEX character 12 name of the MATEX object that will be updated by the module The rod devices material mixtures are appended to the previous material index and the rod devices indices are also modified accordingly descdev structure describing the input data to the DEVINI module 3 4 1 Input data to the DEVINI module The DEVINI module allows the definition of rod type devices made of one or many up to 10 parts as depicted in Fig 1 Table 9 Structure descdev EDIT iprint NUM ROD nrod FADE MOVE dev rod i 1 nrod CREATE ROD GR ngrp rod group i 1 ngrp E IGE 344 level 1 Upper insertion H 19 Lower insertion H EEES ee ew ee amp a Lsup ay level lt 1 3 M
229. t h the mesh dimension and a et the net average surfacic current IGE 344 SEL REFLECTOR HELIOS FILE_CONT_1 ncols nrows part hm_dens bypass FILE_CONT_2 emin FILE_CONT_3 vcool vwatr vmodr venrd vfuel vclad vchan FILE_CONT_4 pitch xbe ybe XS_CONT nside 104 Assembly discontinuity factors are generated using the Selengut normalization This option is available with Saphyb prduced by the APOLLOZ2 lattice code using a 2 level flux calculation with the Method of Characteristics IX o ADF EINEAN Hom dg where g is the energy group f the assembly surface pom asssembly Dg the diffusion coefficient h the mesh dimension and Jj J7 JN the incoming outgoing and net average surfacic currents is the homogene fluxe in keyword used to indicate that the input SAP object contains cross sections for reflector keyword used to indicate that the input data for the HEL file will be set by the user keyword used to set the FILE_CONT_1 block See Ref 40 number of rod columns Default ncols 15 number of rod rows Default nrows 15 index for computed part of assembly 0 1 2 3 whole half quarter eight By default part 3 initial heavy metal density g cm 3 By default hm_dens 2 78613 the saturated moderator density g cm 3 By default bypass 0 73659 keyword used to set the lower energy limits of neutron groups lower energy limits of neutron groups Default
230. t lt vc0 gt gt ADD C lt lt vc0 gt gt MAP REF A lt lt va0 gt gt B SAMEASREF ENDREF SET A lt lt va0 gt gt ADD 7A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B SAMEASREF ENDREF ENDMIX MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET A lt lt va0 gt gt SET C lt lt vc0 gt gt ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF ADD C lt lt vc0 gt gt MAP REF A lt lt va0 gt gt B lt lt vb0 gt gt ENDREF ENDMIX continued on next page IGE 344 74 all parameters explicitly set all parameters in MAP PLANE AXE Fig 9 MACROLIB NCR CPO NMIX 1 MACRO COMPO CPO default MIX 1 SET A lt lt va0 gt gt SET B lt lt vb gt gt SET C lt lt vc gt gt ADD A lt lt va0 gt gt lt lt va gt gt REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF ENDMIX MACROLIB NCR CPO FMAP NMIX 1 MACRO TABLE CPO default B MIX 1 SET A lt lt va0 gt gt ADD A lt lt va0 gt gt MAP REF C lt lt vc0 gt gt B lt lt vb0 gt gt ENDREF ENDMIX Table 55 NCR inputs for instantaneous cases For the TA the burnup variable has no other choice than to be stored in the MAP object MAPFL Then the input files will be only the burnup in MAP all parameters in MAP MACROLIB NCR CPO FMAP NMIX 1
231. t is defined in creation appears only on LHS or modification appears on both LHS and RHS mode An existing geometry previously created in the GEO module is modified Only 3 D Cartesian or 3 D Hexagonal reactor geometries are allowed MATEX character 12 name of a MATEX object to be created by the module GEOMOLD character 12 name of a GEOMETRY object previously created in the GEO module This object must be specified in read only mode appears only on RHS It is copied into GEOM at the beginning of USPLIT module Only 3 D Cartesian or 3 D Hexag onal reactor geometries are allowed desclink structure describing the input data to the USPLIT module 3 2 1 Input data to the USPLIT module Note that the fuel type and reflector type mixture indices are need to be specified explicitly and the input order must be respected Table 6 Structure desclink EDIT iprint NGRP ngrp MAXR maxreg NMIX nmixt NREFL nrefl RMIX mixr i i 1 nrefl NFUEL nfuel FMIX mixf i i 1 nfuel ASBLY IGE 344 where EDIT iprint NGRP ngrp MAXR maxreg NMIX nmixt NREFL nrefl RMIX mixr NFUEL nfuel FMIX mixf ASBLY 16 keyword used to set iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value larger values produce increasing amounts of output keyword used to specify ngrp integer total number of energy gro
232. t the interpolation factors are computed on each axis seperatly IGE 344 84 The first case corresponds to a complete grid represented by a gray paralepiped on Fig 2 and 3 The figure 2 shows that the interpolated value in point V can be obtained directly without MAP object For time average TA computation lets assume that the parameter B represents the burnup and keep this convention for other database structure also In this case the figure 3 shows also that the direct interpolation can be done to compute an average value between the points V and V Note that the TA burnups are stored in the MAP object and are then recovered automatically The second case corresponds to a partial grid where all the lattice computations have been perfomed for several pairs of parameters which are represented as the gray rectangles on Fig 4 and 5 If we use the notations of Fig 4 and 5 the best estimate interpolated values f we can get are given by f fV gt f Ve Vea f VB f Vec f Ve f Vac f Vea f Va Vea f Vec f Va for instataneous f KVV Vb Va HF a Vea F Vb Va f Vee Vac F Vk Va FV bo Vac HV 24 Visa Vh Va F Vba Vaa f Vic Veo VB Va for TA where f represents the average value between two points The third case corresponds to a minimal grid where the lattice computations have been perfomed only for one parameter variation at a time In this case
233. t the before and after information one will be able to identify the new fuel bundles as well as the bundle that have not been moved in the core by the fact that At 0 for burnup before refueling For bundles that have been displaced in the core during refueling then At gt 0 IGE 344 63 4 CROSS SECTION INTERPOLATION MODULES 4 1 The CRE module The CRE module is used for the recovering and interpolation of nuclear properties from one or many COMPO objects originated from the transport calculations using lattice code DRAGON A resulting MACROLIB will be created or updated by the CRE module it will contain the nuclear properties of some selected reactor materials Two types of MACROLIB can be constructed using the CRE module e A MACROLIB that will be constructed for the few reactor materials namely for the devices and or reflector properties It can also be created for the few fuel regions defined in the reactor core This MACROLIB is permitted to be updated for the new properties in the subsequent calls to the CRE module e A fuel map MACROLIB that will be constructed over the fuel lattice only This MACROLIB will contain a set of interpolated fuel properties with respect to the burnup distribution over the fuel lattice and according to the interpolation option defined in the FMAP object The total number of mixtures in the resulting MACROLIB will equal to the total number of fuel bundles Note that the CRE module can be used
234. t the name of the record to be updated will follow name of the record to be updated The authorised values are defined in the table 8 1 5 page 129 P NAME variable To be allowed these values have to be previously defined in the fuel map before calling MCC keyword to specify that an uniform and absolute value is to be set With this param eter the recl of every region will have the same value absolute value in Kelvin if it is a temperature that will be set for rec1 in every region keyword to specify that an uniform variation of rec1 value is to be applied With this parameter the value of rec1 will be altered in every region value of recl variation in Kelvin if it is a temperature that will be applied to every region keyword used to indicate that the values to edit rec1 are to be recovered from rec2 in the same fuel map FLMAP1 For instance it enables the user to set the fuel temperature with the moderator temperature keyword used to indicate that the values to edit recl are to be recovered from rec2 in an other fuel map FLMAP2 For instance it enables the user to set the fuel temperature in FLMAPI with the fuel temperature of FLMAP2 name of the record where the data to perform the update of FLMAP1 is to be recov ered The authorised values are defined in the table 8 1 5 page 129 P NAME variable To be allowed these values have to be previously defined in the fuel map before calling MCC keyword used to indicat
235. the core coolant section and the coolant inlet flow 2 real value set to the core coolant section in m real value set to the coolant flow in m3 hr This value doesn t include the by pass flow The inlet coolant velocity in m s is computed as flow ds 3600 cwsect keyword used to specify the inlet coolant velocity real value set to the inlet coolant velocity in m s keyword used to specify the assembly characteristics 2 real value set to the assembly surface in m This value is equal to the square of an assembly side including the water gap integer value set to the number of active fuel rods in a single assembly integer value set to the number of active guide tubes in a single assembly keyword used to specify the outlet pressure and inlet absolute temperature real value set to the outlet coolant pressure in Pa The pressure along each channel is assumed to be constant and equal to poutlet in permanent conditions real value set to the inlet coolant absolute temperature in K keyword used to set the pin cell radii real value set to the fuel pellet radius in m real value set to the internal clad rod radius in m real value set to the external clad rod radius in m real value set to the guide tube radius in m keyword used to set the oxyde porosity of fuel Porosity affects some built in correla tions used to represent the heat conduction phenomenon in fuel real value set to the oxyde porosity The defau
236. the detector data 23 3 6 The LAC module soso as i i A E A ote 25 310 1 Input data to the LZC modul o o e 25 3 6 2 Description of dev lzc input structure 26 3 6 3 Description of lzc group input structure 2 2 28 dl Whe BRET Modales ceros ti A eh ee a ae e n 29 Sl Input data to the DSET module 29 3 8 The MCC modules Lui RR ada a d a 31 3 91 Main input data to the MCC module 0000 31 3 9 The MOVDEVE modile a corro RE ls oe E 33 2 9 1 Input data to the MOVDEV module o 33 310 The NEWNAGE module gt cua ca aa a ROS ES 35 SLI The FEPGWs module ooo ee OG Pee CER ee EW bo ewe Ia 36 SLI Inputdata tothe PLPOW module ui iii 37 212 The TAVG module 6 eae LE ee ee es a a 39 3 12 1 Input data to the TAVG module aa 0 200005 ae 39 3 12 2 Time average calculation using DONJON o o 40 3 13 The TINST module ie A ee ee a p A a 42 3 13 1 Input data to the TINST module 0 200200 2 42 3 14 The SEM module sss Cee bad ee BA eae dd eed eee ee ww ee 45 3 14 1 Input data to the SIM module o aaa ee kanan 45 Silo The XENON module coco cosita aa a eee Hae es 49 3 15 1 Input data to the XENON module o o 49 3 16 The DETECT module io essa e ma e edha ed e A ane 51 3 16 1 Input data to the DETECT module 51 Sole The CYR module Lv a a at ADLER 53 3 17 1 Input data to the CVR module 000
237. the supercell code xfac corrective factor for delta sigmas real number For DRAGON code xfac is generally set to 2 0 and for MULTICELL code set to 1 0 The default value is 2 0 IGE 344 36 3 11 The FLPOW module The FLPOW module is used to compute and print the flux and power distributions over the reactor core It also computes and prints some additional information for example the fluxes ratios with respect to the thermal energy group fluxes the mean power density the power and flux form factors etc The computed fluxes and powers are printed either on files or on the screen Note that the calculation using the FLPOW module can be performed once the numerical solution has been previously established using the FLUD or KINSOL module According to the user selected module specification the average fluxes and powers can be computed per each fuel region over the fuel lattice and or per each material region over the whole reactor geometry In either case all fluxes are normalized to the given total reactor power corresponding to the reactor nominal conditions at core equilibrium If the reactor is perturbed from its initial state then a new total reactor power can be recomputed and accordingly the flux and power distributions will be updated using the previously computed normalization factor The FLPOW module will create a new POWER object that will store the information related to the reactor fluxes and powers see Section 8 5
238. thermal confuctivity contains the temperature maximum absolute error in K allowed in the solution of the con duction equations e R contains the maximum relative error allowed in the matrix resolution of the conservation equations of the coolant Re contains the relaxation parameter for the multiphysics parameters temperature of fuel and coolant and density of coolant e Rit contains the time in s IGE 344 146 e Ri contains the heat transfer coefficient of the gap if Sf 1 e R A contains the heat transfer coefficient between the clad and fluid if Sf 1 e Ri contains the surface temperature weighting factor of effective fuel temperature for the Rowlands approximation IGE 344 147 8 7 2 The HISTORY DATA sub directory In the HISTORY DATA directory the following sub directories will be found Table 101 Sub directories in HISTORY DATA directory Type Condition Units Comment STATIC PARAM Dir Ncn sub directory containing all the values of the thermal hydraulics parameters computed by the THM module in steady state conditions and sorted channel by channel Each chan nel is identified by an integer numc that can take values between 1 and 9999 For exam ple the first channel is identified by the string character CHANNEL 0001 TIMESTEPnumt Dir Nen sub directories containing all the values of the thermal hydraulics parameters computed by the THM module in transient conditions at a given time in
239. tion eguation maxitl integer value set to the maximum number of iterations for computing the conduction integral The default value is 50 maxit2 integer value set to the maximum number of iterations for computing the center pellet temperature The default value is 50 maxit3 integer value set to the maximum number of iterations for computing the coolant pa rameters mass flux pressure enthalpy and density in case of a transient calculation The default value is 50 ermaxt real value set to the maximum temperature error in K The default value is 1 K ermaxc real value set to the maximum relative error for parameters given by the resolution of flow conservation equations pressure velocity and enthalpy The default value is 1073 RODMESH keyword used to set the radial discretization of pin cells nbl integer value set to the number of discretisation points in fuel The default value is 5 nb2 integer value set to the number of discretisation points in the whole pin cell fuel cladding The default value is 8 FORCEAVE keyword used to force the use of the average approximation during the fuel conductivity evaluation By default a rectangle quadrature approximation is used BOWR keyword used to set a subcooling model based on the Jens amp Lottes correlation with the Bowring model default option SAHA keyword used to set a subcooling model based on the Saha Zuber correlation This option is recommended for BWR applica
240. tions SET PARAM keyword used to indicate the input or modification of the actual values for a param eter specified using its PNAME PNAME keyword used to specify PNAME PNAME character 12 name of a parameter pvalue single real value containing the actual parameter s values Note that this value will not be checked for consistency by the module It is the user responsibility to provide the valid parameter s value which should be consistent with those recorded in the multicompo or Saphyb database IGE 344 112 6 OPTIMIZATION MODULES This section is related to optimization capabilities available in Donjon and based on generalized perturbation theory 3 General information about the generalized perturbation theory can be found in Sect 5 3 of Ref 1 6 1 The DLEAK module The DLEAK module is used to create a delta MACROLIB type L_MACROLIB with respect to leakage information Derivatives of leakage related information recovered from the input MACROLIB are stored in the STEP heteroneneous list components present in the output MACROLIB Derivatives can be taken with respect to a leakage parameter itself Dg or 214 4 or relative to factor u in Dy or 21 9 Note that factor y is not a SPH factor because it multiplies only leakage related parameters One component of the STEP heteroneneous list is created for each value of energy group g and for each value of mixture The calling specifications are Table 72 Structure DLEA
241. to the liquid zone controllers However if the rod type devices are present in the reactor core then they must be specified first i e before the liquid controllers using the DEVINI module see Section 3 4 In the last case the DEVICE object must also appear on the RHS i e in modification mode it will contain the additional and separate information with respect to the liquid zone controllers MATEX character 12 name of the MATEX object that will be updated by the module The lzc devices material mixtures are appended to the previous material index and the Izc devices indices are also modified accordingly desclzc structure describing the input data to the LZC module 3 6 1 Input data to the LZC module Note that the input order must be respected IGE 344 26 Table 16 Structure desclzc EDIT iprint NUM LZC nlzc dev lzc i 1 nlzc CREATE LZC GR ngrp lzc group i 1 ngrp 3 where EDIT keyword used to set iprint iprint integer index used to control the printing on screen 0 for no print 1 for minimum printing default value larger values produce increasing amounts of output NUM LZC keyword used to specify nlzc nlzc integer total number of liguid zone controllers This number must be greater than 0 CREATE keyword used to create the lzc groups of devices The creation of groups is optional LZC GR keyword used to set ngrp ngrp integer total number of the lzc groups to be created T
242. tructure pica History HST History Burnup GET TMod lt lt TMod gt gt CELLID lt lt icha gt gt lt lt ibun gt gt GET TComb lt lt TComb gt gt TCalo lt lt TCalo gt gt where no idfuel is given see Table 41 thus we have used the default value for idfuel 1 to store in HISTORY the general information associated with fuel channel 1 and bundle 1 Here the initial properties associated with fuel type 1 will be generated from the initial isotope densities in the BURNUP For the CELLID here icha 1 and ibun 1 the burnup information isotope densities depletion parameters and initial fuel density are stored in a celldir directory Moreover the power rate 31 9713 kW kg and the depletion time 5 0 days are kept in the PARAMBURNTAR record A HISTORY data structure that contains the initial cell information can be updated using a MAP data structure duna Map data structure for refueling MAP1 pat SEQ_ASCII MAP1 S Refuel goes Reseau MAP1 History HST History Reseau Here new burnup power ratings will be stored in the HISTORY data structure reflecting the power distribution in the DONJON calculation The refueling information available in the MAP structure will also be used to redistribute the fuel in the HISTORY structure at various cell location Finally the last option is to recover this information in DRAGON to perform a new series of cell calculations k Local parameters Initial
243. uarter voiding pattern According to this pattern the fuel mixtures will be modified only for the upper left quarter of reactor channels in which the direction of coolant flow is positive IGE 344 CHAN VOID nvoid YNAME XNAME 55 keyword used to specify the user defined voiding pattern Each voided channel must be identified by its YNAME name followed by its XNAME integer total number of the voided channels This number must be greater than 0 and less than or equal to the total number of reactor channels character 2 vertical name of the voided channel A vertical channel name is identified by the channel row using an alphabetical letter A B C etc The total number of the specified Y names must equal to the total number of voided channels nvoid character 2 horizontal name of the voided channel A horizontal channel name is identified by the channel column using a numerical character 1 2 3 etc The to tal number of the specified X names must equal to the total number of voided channels nvoid IGE 344 56 3 18 The HST module The HST module has been designed to manage a full reactor execution in DONJON using explicit DRAGON calculations for each cell This module saves in an HISTORY data structure the informa tion available in BURNUP data structures generated by DRAGON It can also read MAP data structure generated by DONJON to prepare the HISTORY data structure for a new series of c
244. ues of the decision variables The maximum values of the decision variables can be The minimum values of the decision variables can be The weight of the decision variables w in the quadratic constraint The limit value of the contraints The units depends with the type of the constraint type The type of the contraints 1 for gt 0 for 1 for lt The weight of the constraint 7 and y for the duals and meta heuristic methods The actual values of the objective function first value and the contraints the following values The number of the constraints are assigned in the order they have been defined The different limits and values for the iterative calculations of the optimization problem The gradients of the objective function and the constraints The gradients of the objec tive for all the decision variables are in first position then follow the gradients of the con straints Directory containing differents informations of the previous iterations the values of the deci sion variables the objective function the con straints and the gradients of these functions for the previous external iteration This reper tory will be created by the module QLP unless it is specified to not do The array OPT PARAM R contains real values related with the different limits and values for the iterative calculations of the optimization problem IGE 344 150 lst S initial radius of the quadratic constraint default 1 0
245. ups This value must be greater than 0 keyword used to specify maxreg integer maximum number of mesh splitted regions in the reactor geometry In 3 D Hezagonal geometry it corresponds to the total number of prismatic blocks In l keyword used to extend number of material mixtures in case new fuels are going to be inserted in the fuel map in upcoming fuel cycles By default nmixt is set to the maximum mixture index in RHS geometry GEOM or GEOMOLD the maximum fuel mixture index in the complete life of the reactor This number must be greater than the maximum mixture index in RHS geometry GEOM or GEOMOLD keyword used to specify nrefl integer total number of reflector types A reactor should have at least one reflector material keyword used to specify mixr integer array of the reflector type mixture indices Each reflector type is assigned a distinct mixture number as previously defined in the GEOMETRY object keyword used to specify nfuel integer total number of fuel types A reactor should have at least one fuel type keyword used to specify mixf integer array of the fuel type mixture indices Each fuel type is assigned a distinct mixture number as previously defined in the GEOMETRY object keyword used to compute automatically nfuel and mixf i This option is only available when the geometry has been split by the NAP module IGE 344 17 3 3 The MACINI module The MACINI module is used to construct an extended MACR
246. used at step one and the general LEMKE method at step two keyword used to specify that the MAP method will be used The guadratic constraint is linearized and a converging seguence of SIMPLEX calculations is performed keyword used to specify that the augmented Lagrangian method will be used keyword used to specify that the penalty method will be used keyword used to set the initial radius of the guadratic constraint default value is sr 1 0 initial radius of the guadratic constraint real IGE 344 DUT STEP EPS Eext INN STEP EPS Einn CST QUAD EPS Equad MAXIMIZE MINIMIZE STEP REDUCT HALF PARABOLIC VAR VALUE control VAR WEIGHT weight VAR VAL MIN vecmin varmin VAR VAL MAX vecmax varmax FOBJ CST VAL funct CST TYPE type CST OBJ 116 keyword used to set the tolerance of outer iteration convergence inside module PLQ tolerance value real keyword used to set the tolerance used within the SIMPLEX or LEMKE method tolerance value real keyword to set the convergence parameter epsilon4 for the radius of the quadratic constraint inside module GRAD tolerance for convergence of the radius of the quadratic constraint real keyword used to specify that the optimization problem will be a maximization keyword used to specify that the optimization problem will be a minimization default keyword used to define the method of the reduction of the outer step keyword used to speci
247. used to compute a cross section in the PMAXS formalism including 4 state variables with r the reference state and m the mid point between the reference state and the current node state CR DC VTF TC S CR DC VTF TC S DC VTF TC ADC a EN 0DC CR DCm YTF TC aya E ATC DE OVTF CR DC TF TC OTC oR DC VTF TC The PMAXS formalism and the procedure branching generation are completely described in the GenPMAXS manual 9 The PMAXS file contains eight blocks few of them are optional and others mandatory A precise description of each block is proposed in the APPENDIX A of GenPMAXS manual In this section a short description of blocks is given Block XS CONTROL information mandatory The first block stores data reflecting the conditions in which cross sections are obtained and what kind of informations is contained in the PMAXS file It is composed of five integers for the number of energy groups of delay neuton groups etc Then fifteen logical flags indicates if the PMAXS contains the correponding data such as assembly discontinuity factor ADF Xe and Sm microscopic cross sections The block is ended by five lines of comments to be filled by the user GLOBAL_V 1266114517 FFFFFFFFFFFFFFT Contents of T H Invariant Variables TIV block and Cross Sections XS block TIV XS tr ab nf kf sct 2 Group value of each variable are put together in a line Some variables separated by share a line
248. vity as a function of a simple user defined correlation IGE 344 e S flag to set the clad thermal conductivity th 0 built in correlation is used 33 145 1 set the clad thermal conductivity as a function of a simple user defined correlation S type of approximation used during the fuel conductivity evaluation sth 0 use a rectangle quadrature approximation 14 1 use an average approximation Si type of subcooling model 0 use the Jens Lottes correlation and Bowring s model th _ 15 1 use the Saha Zuber subcooling model if Sh 1 if Sh 1 The array Ri contains the following information e Rih e R e RE e R P e RE e Re e RE e Re e Re e Ri e Ri Ris Ri e Rii contains the current time step in s contains the fraction of reactor power released in fuel contains the critical heat flux in W m contains the inlet coolant velocity in m s contains the outlet coolant pressure in Pa contains the inlet coolant temperature in K contains the Plutonium mass fraction in fuel contains the fuel porosity contains the fuel pellet radius contains the internal clad rod radius in m contains the external clad rod radius in m contains the guide tube radius in m contains the assembly surface in m S contains the number of terms in the user defined correlation for the fuel thermal confuctivity S contains the number of terms in the user defined correlation for the clad
249. with each condensed group keyword to specify that the complete contents of Tape16 must be listed on the output file keyword to specify that the remaining information will be associated with mixture properties definition character 6 name of the mixture to create or update on the CPO optional keyword to specify that cell averaged data will be taken from Tape16 This is the default option optional keyword to specify that regional data will be taken from Tape16 The default option is CELLAV region number associated with this material in Tape16 optional keyword to specify that the cross section taken from Tape16 are at reference value This information must be defined at least once for each mixture It must also precede the definition of perturbation parameters number of consecutive burnup steps associated with mixture The default value is nburn 1 We will assume that the same number of burnup steps is also available for the nuclear properties associated with the perturbed local parameters first Tape16 record number associated with this mixture character 2 name of the perturbation Each perturbation is associated with a single local parameter The values permitted for NAMPER are the following FT for fuel temperature MT for moderator temperature MD for moderator density MP for moderator purity MB for moderator boron CT for coolant temperature CD for coolant density CP for coolant purity 5 00 SS SO iS a RT fo
250. yword SAME is specified then this average value will be set for all fuel bundles for every reactor channel IGE 344 CHAN BUND TIMES PNAMEREF pvalue OLDMAP FUEL WEIGHT ENRICH POISON fvalue nfuel CELL ialch 14 keyword used to indicate that the values of a local parameter will be provided per each reactor channel If the keyword CHAN is specified then the channel averaged parameter s value will be set for all fuel bundles containing in the same reactor channel keyword used to indicate that the values of a local parameter will be specified per each fuel bundle for every channel keyword used to indicate that the values of the local parameter PNAME is a translation of the local parameter PNAMEREF via a multiplication of the constant indicated by SAME character 12 identification name of a given parameter real array or a single value containing the actual parameter s values Note that these values will not be checked for consistency by the module It is the user responsibility to provide the valid parameter s values which should be consistent with those recorded in the multicompo database keyword to specify that the pvalue value s is are recovered from FLMAP2 keyword used to indicate the input of data which will be specified for each fuel type keyword used to indicate the input of fuel weight in a bundle given in kg keyword used to indicate the input of fuel enrichment values given in wt k
251. yword used to specify the rod group irgrp number integer identification number of a rod group of devices that will be modifed with the same parameters Each rod group has a unique irgrp number ranging from 1 to ngrp as been defined in the DEVINI module see Section 3 4 3 keyword used to specify the liquid controller ilzc number integer identification number of a liquid controller to be modified Each lzc type device has a unique ilzc number ranging from 1 to nlzc as been defined in the LZC module see Section 3 6 2 keyword used to specify the lzc group ilgrp number integer identification number of a lzc group of devices that will be modifed with the same parameters Each lzc group has a unique ilgrp number ranging from 1 to ngrp as been defined in the LZC module see Section 3 6 3 keyword used to specify a new level value real positive value of the new device level For the rod type devices this value must correspond to the new rod insertion level see Section 3 4 3 For the lzc type devices this value must correspond to the new water filling level see Section 3 6 2 In any case the new level value must be 0 0 lt value lt 1 0 keyword used to specify a new value for speed real positive value of the device speed For the rod type devices this value must correspond to the speed of rod movement insertion or extraction given in cm s For the lzc type devices this value must correspond to the water filling rate

A USER GUIDE FOR DONJON VERSION5 A. Hébert, D. Sekki, and

Contents

Download Pdf Manuals

Related Search

Related Contents

A USER GUIDE FOR DONJON VERSION5 A. H&eacute;bert, D. Sekki, and

Contents

Download Pdf Manuals

Related Search

Related Contents

A USER GUIDE FOR DONJON VERSION5 A. Hébert, D. Sekki, and