diff options
Diffstat (limited to 'Donjon/data/TC_EVO.x2m')
| -rw-r--r-- | Donjon/data/TC_EVO.x2m | 171 |
1 files changed, 171 insertions, 0 deletions
diff --git a/Donjon/data/TC_EVO.x2m b/Donjon/data/TC_EVO.x2m new file mode 100644 index 0000000..8e77940 --- /dev/null +++ b/Donjon/data/TC_EVO.x2m @@ -0,0 +1,171 @@ +*********************************************************** +* * +* Input file : TC_EVO.x2m * +* Purpose : Test-case for micro-depletion * +* calculation using DONJON4 * +* Author(s) : M. Guyot (2010) * +* * +* To run : rdonjon4 TC_EVO.x2m TC_EVO_CPOs * +* assertS.c2m * +* * +*********************************************************** +*********************************************************** +* A mixed method is used in this TC : the micro-depletion * +* is used to deplete the isotopes, except the residual * +* An interpolation method is used to compute the residual* +* properties at each step of time * +*********************************************************** + +LINKED_LIST GEOM GEOEVO GEOM2 GEOEVO2 MATEX MATEX2 TRACK FMAP MICROLIB + MICRO MACROLIB SYSTEM FLUX BURNUP CPO CPO2 MACRO MACRO3 + TRACKEVO POWER NEWFLUX ; +MODULE GEO: USPLIT: TRIVAT: RESINI: NCR: TRIVAA: FLUD: EVO: DELETE: + TINST: MACINI: END: GREP: FLPOW: UTL: ; +PROCEDURE assertS ; +REAL + Power Delt Timec Timei Timef := + 31.9713 1.0 1.0 0.0 0.0 ; +REAL Keff ; +REAL PTOT := 1.9039 ; ! power used for the interpolation +*--- +* Set the 3D Geometry +*--- +GEOM := GEO: :: CAR3D 2 2 1 + X- REFL X+ REFL MESHX 0.0 28.575 57.15 + SPLITX 1 1 + Y- REFL Y+ REFL MESHY 0.0 28.575 57.15 + SPLITY 1 1 + Z- REFL Z+ REFL MESHZ 0.0 49.53 SPLITZ 1 + MIX + PLANE 1 + 1 2 + 1 1 ; +GEOM2 MATEX := USPLIT: GEOM :: NGRP 2 MAXR 2000 + NREFL 1 RMIX 2 + NFUEL 1 FMIX 1 ; +TRACK := TRIVAT: GEOM2 :: MAXR 2000 MCFD 2 ; + +*--- +* Set the fuel-map and the burnup of the assembly +*--- +FMAP MATEX := RESINI: MATEX :: + ::: GEO: CAR3D 2 2 1 + X- REFL X+ REFL MESHX 0.0 28.575 57.15 + SPLITX 1 1 + Y- REFL Y+ REFL MESHY 0.0 28.575 57.15 + SPLITY 1 1 + Z- REFL Z+ REFL MESHZ 0.0 49.53 SPLITZ 1 + MIX + PLANE 1 + 1 0 + 1 1 ; + NXNAME '01' '02' + NYNAME 'A' 'B' + NCOMB ALL + FUEL WEIGHT 19.85 ; +FMAP := RESINI: FMAP :: EDIT 1 + BTYPE INST-BURN + INST-BVAL BUND + 2000.0 2000.0 2000.0 ; + +*--- +* Recover the fuel-map geometry to prepare the depletion calculation +*--- +GEOEVO := FMAP :: STEP UP 'GEOMAP' ; +GEOEVO2 MATEX2 := USPLIT: GEOEVO :: NGRP 2 MAXR 2000 + NFUEL 1 FMIX 1 ; +TRACKEVO := TRIVAT: GEOEVO2 :: MAXR 2000 MCFD 1 ; + +*--- +* Recover the reflector and fuel properties +*--- +SEQ_ASCII CPOFUEL CPOREFL ; +CPO := CPOFUEL ; +CPO2 := CPOREFL ; + +MICRO := NCR: CPO2 :: EDIT 0 + NMIX 2 MICRO CUBIC + COMPO CPO2 'moderator' + MIX 2 FROM 1 + SET 'BURN' 0. + MICRO ALL + ENDMIX ; +MICROLIB := NCR: CPO FMAP :: EDIT 0 + MICRO CUBIC ALLX 5 + TABLE CPO 'fuel' BURN + MIX 1 FROM 1 + INST-BURN + MICRO ALL + ENDMIX ; +MACRO := MICRO :: STEP UP 'MACROLIB' ; +MACROLIB := MICROLIB :: STEP UP 'MACROLIB' ; + +*---- +* Flux calculation and initialization of the BURNUP +*---- +MACRO3 MATEX := MACINI: MATEX MACRO MACROLIB :: EDIT 1 ; +SYSTEM := TRIVAA: MACRO3 TRACK :: EDIT 0 ; +FLUX := FLUD: SYSTEM TRACK :: EDIT 0 EXTE 200 1.E-6 ADI 5 ; +POWER NEWFLUX := FLPOW: FMAP FLUX TRACK MATEX :: + EDIT 0 PTOT <<PTOT>> BUND ; + +UTL: MICROLIB :: DUMP ; + +BURNUP MICROLIB := EVO: MICROLIB POWER :: + EDIT 2 FLUX_POW PIFI SAVE 0.0 DAY POWR <<Power>> ; + +GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 >>Keff<< ; +assertS FLUX :: "K-EFFECTIVE" 1 1.047273 ; + +WHILE Timei Timec < DO + EVALUATE Timef := Timei Delt + ; +* MICRO-DEPLETION + BURNUP MICROLIB := EVO: BURNUP MICROLIB POWER :: + EDIT 2 FLUX_POW PIFI DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; +* INTERPOLATION FOR THE MICROSCOPIC CROSS SECTIONS +* The RES keyword limits multiparameter interpolation to microscopic XS + FMAP := TINST: FMAP POWER :: EDIT 9 TIME <<Delt>> DAY ; + MICROLIB := NCR: MICROLIB CPO FMAP :: EDIT 0 + MICRO CUBIC ALLX 5 RES + TABLE CPO 'fuel' BURN + MIX 1 + INST-BURN + MICRO ALL + ENDMIX ; + + MACRO3 SYSTEM FLUX POWER NEWFLUX MACROLIB := DELETE: + MACRO3 SYSTEM FLUX POWER NEWFLUX MACROLIB ; +* FLUX CALCULATION + MACROLIB := MICROLIB :: STEP UP 'MACROLIB' ; + MACRO3 MATEX := MACINI: MATEX MACRO MACROLIB :: EDIT 1 ; + SYSTEM := TRIVAA: MACRO3 TRACK :: EDIT 0 ; + FLUX := FLUD: SYSTEM TRACK :: EDIT 0 EXTE 200 1.E-6 ADI 5 ; +* POWER CALCULATION + POWER NEWFLUX := FLPOW: FMAP FLUX TRACK MATEX :: + EDIT 0 PTOT <<PTOT>> BUND ; + + BURNUP MICROLIB := EVO: BURNUP MICROLIB POWER :: + FLUX_POW PIFI SAVE <<Timef>> DAY POWR <<Power>> ; +*---- +* change delta t for burnup and final time if required +*---- + IF Timef Timec = THEN + IF Timec 10.0 = THEN + EVALUATE Delt Timec := 10.0 50.0 ; + ENDIF ; + IF Timec 5.0 = THEN + EVALUATE Delt Timec := 5.0 10.0 ; + ENDIF ; + IF Timec 1.0 = THEN + EVALUATE Delt Timec := 4.0 5.0 ; + ENDIF ; + ENDIF ; + EVALUATE Timei := Timef ; +ENDWHILE ; + +GREP: FLUX :: GETVAL 'K-EFFECTIVE' 1 >>Keff<< ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.033521 ; + +ECHO "test TC_EVO completed" ; +END: ; +QUIT . |