diff options
Diffstat (limited to 'Dragon/data/lumpSS.x2m')
| -rw-r--r-- | Dragon/data/lumpSS.x2m | 268 |
1 files changed, 268 insertions, 0 deletions
diff --git a/Dragon/data/lumpSS.x2m b/Dragon/data/lumpSS.x2m new file mode 100644 index 0000000..b070c3a --- /dev/null +++ b/Dragon/data/lumpSS.x2m @@ -0,0 +1,268 @@ +*---- +* +* TEST CASE lumpSS +* DISTRIBUTED SELF-SHIELDING ON LUMPED ASSEMBLIES +* 295-GROUP JEFF3.1.1 DRAGLIB +* +* Author: A. Hebert +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST GEOM1 GEOM2 GEOM12 TRACK_SS1 TRACK_SS2 TRACK_2L LIBRARY LIBRARY2A LIBRARY2B + LIBRARY2 ASB FLUX EDIBU BURNUP MCOMPO StepList ; +MODULE LIB: GEO: SYBILT: ASM: FLU: COMPO: EDI: EVO: UTL: GREP: DELETE: END: ; +PROCEDURE assertS ; +PROCEDURE mixA1_lumpSS ; +REAL P_Zr90 P_Zr91 P_Zr92 P_Zr94 P_Zr96 := + 51.45E-02 11.22E-02 17.15E-02 17.38E-02 2.80E-02 ; +REAL N_Zr0 := 4.3241E-2 ; +REAL N_Zr90 N_Zr91 N_Zr92 N_Zr94 N_Zr96 := + N_Zr0 P_Zr90 * N_Zr0 P_Zr91 * N_Zr0 P_Zr92 * N_Zr0 P_Zr94 * + N_Zr0 P_Zr96 * ; + +************************************************************************ +* BEGIN DECLARATION * +************************************************************************ +INTEGER maxstep := 3 ; +StepList := UTL: :: + CREA + ListBU <<maxstep>> = + 0.0 10.0 30.0 + ; +INTEGER maxautop := 1 ; +StepList := UTL: StepList :: + CREA + ListAutop <<maxautop>> = + 30.0 + ; +************************************************************************ +* END DECLARATION * +************************************************************************ + +************************************************************************ +* COMPO INITIALIZATION * +************************************************************************ +MCOMPO := COMPO: :: + EDIT 10 + STEP UP 'EDIBU_HOM' + MAXCAL 5 + COMM 'Evolution burnup - Condensed - Assembly scale homogenized' ENDC + PARA 'Burnup' IRRA + ISOT 11 U235 U236 U238 Pu239 Pu240 Pu241 Pu242 Gd155 Gd157 Xe135 Sm149 + INIT +; + +************************************************************************ +* DRAGLIB PROCESSING * +************************************************************************ +LIBRARY := LIB: :: + EDIT 1 + NMIX 11 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 (*PROCESS SCATTERING DATA TO P1 ORDER*) + PT (*CALENDF TYPE PROBABILITY TABLES*) +* + DEPL LIB: DRAGON FIL: DLIB_295 + MIXS LIB: DRAGON FIL: DLIB_295 + MIX 1 293.0 + O16 = O16 4.6624E-2 + U235 = U235 7.0803E-4 1 + U238 = U238 2.2604E-2 1 + Pu239 = Pu239 1.0000E-9 1 + MIX 2 COMB 1 1.0 + MIX 3 COMB 1 1.0 + MIX 4 COMB 1 1.0 + MIX 5 293.0 + Zr90 = Zr90 <<N_Zr90>> 2 IRSET PT NONE + Zr91 = Zr91 <<N_Zr91>> 2 IRSET PT NONE + Zr92 = Zr92 <<N_Zr92>> 2 IRSET PT NONE + Zr94 = Zr94 <<N_Zr94>> + Zr96 = Zr96 <<N_Zr96>> + MIX 6 293.0 NOEV + H1 = H1_H2O 4.6892E-2 + O16 = O16 2.3446E-2 + MIX 7 293.0 + O16 = O16 4.6624E-2 + U235 = U235 7.0803E-4 1 + U238 = U238 2.2604E-2 1 + Pu239 = Pu239 1.0000E-9 1 + MIX 8 COMB 7 1.0 + MIX 9 COMB 7 1.0 + MIX 10 COMB 7 1.0 + MIX 11 COMB 5 1.0 + ; +************************************************************************ +* GEOMETRY DEFINITION AND TRACKING * +************************************************************************ +GEOM1 := GEO: :: CARCEL 6 + X- REFL X+ REFL Y- REFL Y+ REFL + MESHX 0.0 1.2 MESHY 0.0 1.2 + RADIUS 0.0 0.288712 0.365195 0.397962 0.4083 0.45 0.5748331 + MIX 1 2 3 4 5 6 6 ; +TRACK_SS1 := SYBILT: GEOM1 :: + MAXR 20 MAXZ 15000 DP00 QUA2 20 3 ; + +GEOM2 := GEO: :: CARCEL 6 + X- REFL X+ REFL Y- REFL Y+ REFL + MESHX 0.0 1.2 MESHY 0.0 1.2 + RADIUS 0.0 0.288712 0.365195 0.397962 0.4083 0.45 0.5748331 + MIX 7 8 9 10 11 6 6 ; +TRACK_SS2 := SYBILT: GEOM2 :: + MAXR 20 MAXZ 15000 DP00 QUA2 20 3 ; + +GEOM12 := GEO: :: CAR2D 2 1 + X- REFL X+ REFL Y- REFL Y+ REFL + CELL GEOM1 GEOM2 + ::: GEOM1 := GEO: CARCEL 6 + MESHX 0.0 1.2 MESHY 0.0 1.2 + RADIUS 0.0 0.288712 0.365195 0.397962 0.4083 0.45 0.5748331 + MIX 1 2 3 4 5 6 6 ; + ::: GEOM2 := GEO: CARCEL 6 + MESHX 0.0 1.2 MESHY 0.0 1.2 + RADIUS 0.0 0.288712 0.365195 0.397962 0.4083 0.45 0.5748331 + MIX 7 8 9 10 11 6 6 ; +; +TRACK_2L := SYBILT: GEOM12 :: + MAXR 20 MAXZ 30000 DP00 QUA2 20 3 ; + +************************************************************************ +* NOMINAL REFERENCE CALCULATION * +************************************************************************ +INTEGER istepNext ; +REAL BUbeg Tbeg ; +REAL BU BUautop BUcompo := 0.0 0.0 0.0 ; +REAL Norm_f2 := 39.0 ; ! SPECIFIC POWER MW/t +REAL BUend Tend := 0. 0. ; +INTEGER istep := 0 ; +INTEGER iauto := 1 ; + +********************************************** +* BEGIN DEPLETION * +********************************************** +ECHO "start burnup loop" ; +WHILE istep maxstep < DO + + EVALUATE istep := istep 1 + ; + ECHO "burnup step: " istep "/" maxstep ; + + EVALUATE BUbeg Tbeg := BUend Tend ; + GREP: StepList :: GETVAL 'ListBU' <<istep>> >>BUend<< ; + EVALUATE Tend := BUend Norm_f2 / ; + ECHO "Nominal depletion between" Tbeg " and" Tend "days" ; + +**************************** +* BEGIN SELF-SHIELDING * +**************************** + IF istep 1 = THEN + ECHO "perform initial self-shielding at" Tbeg "days" ; + + ECHO "$$$ ------------- begin : mixAUTOP1 " ; + LIBRARY := UTL: LIBRARY :: DIR ; + LIBRARY2 LIBRARY2A LIBRARY2B := mixA1_lumpSS LIBRARY TRACK_SS1 TRACK_SS2 :: + <<istep>> ; + ECHO "$$$ ------------- end : mixAUTOP1 " ; + + ELSE + GREP: StepList :: GETVAL "ListAutop" <<iauto>> >>BUautop<< ; + GREP: StepList :: GETVAL "ListBU" <<istep>> >>BU<< ; + IF BUautop BU = THEN + + ECHO "perform self-shielding at" Tbeg "days" ; + ECHO "$$$ ----------- begin : mixAUTOP2 " ; + LIBRARY2 LIBRARY2A LIBRARY2B := + mixA1_lumpSS LIBRARY2 LIBRARY2A LIBRARY2B LIBRARY TRACK_SS1 TRACK_SS2 :: + <<istep>> ; + ECHO "$$$ ----------- end : mixAUTOP2 " ; + + EVALUATE iauto := iauto 1 + ; + + ENDIF ; + ENDIF ; + +**************************** +* END OF SELF-SHIELDING * +**************************** + +************************************************************ +* BEGIN FLUX CALCULATION * +************************************************************ + ECHO "flux calculation at" Tbeg "days" ; + ASB := ASM: LIBRARY2 TRACK_2L :: PIJ ; + IF istep 1 = THEN + FLUX := FLU: LIBRARY2 TRACK_2L ASB :: TYPE K ; + ELSE + FLUX := FLU: FLUX LIBRARY2 TRACK_2L ASB :: TYPE K ; + ENDIF ; + ASB := DELETE: ASB ; + +************************************************************ +* NON-REGRESSION * +************************************************************ + IF istep 1 = THEN + assertS FLUX :: K-EFFECTIVE 1 1.329496 ; + ELSEIF istep 2 = THEN + assertS FLUX :: K-EFFECTIVE 1 1.314796 ; + ELSEIF istep 3 = THEN + assertS FLUX :: K-EFFECTIVE 1 1.292124 ; + ENDIF ; + +************************************************************ +* BEGIN FLUX NORMALIZATION * +************************************************************ + IF istep 1 = THEN + BURNUP LIBRARY2 := EVO: LIBRARY2 FLUX TRACK_2L :: + EDIT 3 + SAVE <<Tbeg>> DAY POWR <<Norm_f2>> + ; + ELSE + BURNUP LIBRARY2 := EVO: BURNUP LIBRARY2 FLUX TRACK_2L :: + EDIT 3 + SAVE <<Tbeg>> DAY POWR <<Norm_f2>> + ; + ENDIF ; + +************************************************************ +* EDITION * +************************************************************ + EDIBU := EDI: FLUX LIBRARY2 TRACK_2L :: + EDIT 0 + MICR 11 U235 U236 U238 Pu239 Pu240 Pu241 Pu242 Gd155 Gd157 Xe135 Sm149 + MERG MIX + 1 1 1 1 0 0 1 1 1 1 0 + COND + SAVE ON FUEL ; + +************************************************************ +* FILL MULTICOMPO * +************************************************************ + ECHO "fill multicompo at" Tbeg "days" ; + MCOMPO := COMPO: MCOMPO EDIBU BURNUP LIBRARY2 :: + EDIT 0 + STEP UP 'EDIBU_HOM' + SET <<Tbeg>> DAY ; + + EDIBU := DELETE: EDIBU ; + +********************************************************* +* BEGIN DEPLETION * +********************************************************* + IF istep maxstep < THEN + EVALUATE istepNext := istep 1 + ; + GREP: StepList :: GETVAL 'ListBU' <<istepNext>> >>BUend<< ; + EVALUATE Tend := BUend Norm_f2 / ; + + BURNUP LIBRARY2 := EVO: BURNUP LIBRARY2 FLUX TRACK_2L :: + EDIT 0 + NOEX DEPL <<Tbeg>> <<Tend>> DAY + POWR <<Norm_f2>> + EXPM 1.0E15 GLOB ; + + ENDIF ; + +********************************************************* +* END BURNUP LOOP * +********************************************************* +ENDWHILE ; +ECHO "test lumpSS completed" ; +END: ; |