diff options
Diffstat (limited to 'Dragon/data/tdraglib_proc/TCDR33.c2m')
| -rw-r--r-- | Dragon/data/tdraglib_proc/TCDR33.c2m | 146 |
1 files changed, 146 insertions, 0 deletions
diff --git a/Dragon/data/tdraglib_proc/TCDR33.c2m b/Dragon/data/tdraglib_proc/TCDR33.c2m new file mode 100644 index 0000000..e7f8ed4 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR33.c2m @@ -0,0 +1,146 @@ +*---- +* TEST CASE TCDR33 +* HOMOGENEOUS DEPLETION CASE WITH USER-DEFINED CHAIN +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: None +* +*---- +* Define variables +* Burnup paremeters +* a) Power +* = 600.0 kw/kg for 0.0 to 1000.0 days +* = 0.0 kw/kg for 1000.0 to 2000.0 days +* b) Burnup time interval Delt +* = 10 days for 0 to 50 days +* = 50 days for 50 to 500 days +* = 100 days for 500 to 1000 days +* = 1000 days for 1000 to 2000 days +* c) Editing time Timec +* = 0.0, 50.0, 500.0, 1000.0 and 2000.0 days +* d) Burnup control time variables Timei, Timef, TotalTime +* Timei = initial time +* Timef = final time +* TotalTime = Final time reached +* d) Print variable Iprint +* = 1 reduced print +* = 3 full print +*---- +REAL + Power Delt Timec Timei Timef TotalTime := + 600.0 10.0 50.0 0.0 0.0 2000.0 ; +INTEGER + Iprint := 1 ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY HOM TRACK PIJ FLUX BURNUP EDITION ; +MODULE + GEO: SYBILT: LIB: SHI: ASM: FLU: EVO: EDI: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 1 CTRA WIMS + DEPL 100 CHAIN + Am243 NFTOT 222.2 NG + FROM NG 1.0000E+00 Pu242 NG 1.0000E+00 Am242 + Am242 DECAY 1.55780E-02 NFTOT NG + FROM NG 2.5000E-01 Am241 + Am241 NFTOT 217.1 NG + FROM DECAY 1.0000E+00 Pu241 + Pu242 NFTOT 217.0 NG N2N + FROM NG 1.0000E+00 Pu241 + Pu241 DECAY 1.52540E-01 NFTOT 214.0 NG N2N + FROM NG 1.0000E+00 Pu240 + Pu240 NFTOT 214.1 NG N2N + FROM NG 1.0000E+00 Np239 NG 1.0000E+00 Pu239 + Pu239 NFTOT 211.3 NG N2N + FROM DECAY 1.0000E+00 Np239 + Np239 DECAY 3.39940E+02 NFTOT NG N2N + FROM NG 1.0000E+00 U238 + U238 NFTOT 212.2 NG N2N + Np237 NFTOT 210.3 NG + FROM NG 1.0000E+00 U236 + U236 NFTOT 200.3 NG + FROM NG 1.0000E+00 U235 + U235 NFTOT 202.6 NG N2N + FROM NG 1.0000E+00 U234 + U234 NFTOT 203.4 NG + ENDCHAIN + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 300.0 H1 = H1 2.00000E+1 + U235 = U235 1.0 1 + U236 = U236 0.0 1 + ; +*---- +* Geometry HOM : Homogeneous geometry +*---- +HOM := GEO: :: HOMOGE + MIX 1 ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for keff +*---- +TRACK := SYBILT: HOM :: + TITLE 'TCDR33: HOMOGENEOUS BENCHMARK WITH BURNUP' ; +LIBRARY := SHI: LIBRARY TRACK :: EDIT 0 NOLJ ; +PIJ := ASM: LIBRARY TRACK ; +FLUX := FLU: PIJ LIBRARY TRACK :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.873188 ; +EDITION := EDI: FLUX LIBRARY TRACK :: + COND 4.0 MERGE COMP SAVE ; +*---- +* Burnup loop: for first step BURNUP is created +* while for other steps it is modified +*---- +WHILE Timei TotalTime < DO + EVALUATE Timef := Timei Delt + ; + IF Timef Timec = THEN + EVALUATE Iprint := 3 ; + ELSE + EVALUATE Iprint := 1 ; + ENDIF ; + IF Timei 0.0 = THEN + BURNUP LIBRARY := EVO: LIBRARY FLUX TRACK :: + EDIT <<Iprint>> DEPL <<Timei>> <<Timef>> DAY + POWR <<Power>> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK :: + EDIT <<Iprint>> NOEX DEPL <<Timei>> <<Timef>> DAY + POWR <<Power>> ; + ENDIF ; + LIBRARY := SHI: LIBRARY TRACK :: EDIT 0 NOLJ ; + PIJ := DELETE: PIJ ; + PIJ := ASM: LIBRARY TRACK ; + FLUX := FLU: FLUX PIJ LIBRARY TRACK :: + TYPE K ; + IF Iprint 3 = THEN + EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + SAVE ; + ENDIF ; +*---- +* change delta t for burnup, final time and power if required +*---- + IF Timef Timec = THEN + IF Timec 1000.0 = THEN + EVALUATE Power Delt Timec := 0.0 1000.0 2000.0 ; + ENDIF ; + IF Timec 500.0 = THEN + EVALUATE Delt Timec := 100.0 1000.0 ; + ENDIF ; + IF Timec 50.0 = THEN + EVALUATE Delt Timec := 50.0 500.0 ; + ENDIF ; + ENDIF ; + EVALUATE Timei := Timef ; +ENDWHILE ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.760241 ; +ECHO "test TCDR33 completed" ; +END: ; +QUIT "LIST" . |