diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Dragon/data/tdraglib_proc | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Dragon/data/tdraglib_proc')
42 files changed, 6165 insertions, 0 deletions
diff --git a/Dragon/data/tdraglib_proc/EvolRef.c2m b/Dragon/data/tdraglib_proc/EvolRef.c2m new file mode 100644 index 0000000..acddbe2 --- /dev/null +++ b/Dragon/data/tdraglib_proc/EvolRef.c2m @@ -0,0 +1,120 @@ +*DECK EvolRef.c2m +*---- +* Name : EvolRef.c2m +* Type : Procedure +* Usage : Burnup calculation with DRAGON +* Author : G. Marleau +* Date : 2004/11/22 +* Updated : a) E. Varin 2005/04/03 +* b) E. Varin 2005/06/01 +* c) A. Hebert 2007/01/31 +* +* Procedure called as: +* +* CompoDS ConcIso Flux := EvoRef MicLib TrackingS TrackingF +* IntlineS IntlineF :: +* <<BExit>> <<Power>> <<NbSteps>> ; +* Output structures : +* CompoDS : Multiparameter compo data structure. +* ConcIso : Burnup data structure. +* Flux : Last flux computed. +* Input structures : +* MicLib : Microscopic cross section library. +* TrackingS : Tracking data structures for self shielding. +* TrackingF : Tracking data structures for flux solution. +* IntlineS : File with integration lines for self shielding. +* IntlineF : File with integration lines for flux solution. +* Input variables : +* BExit : Average exit burnup (MWd/T). +* Power : Burnup power (kw/kg) . +* NbSteps : Number of burnup steps. +* +*---- +* Definition of data structures for procedure +* +PARAMETER CompoDS ConcIso Flux MicLib + TrackingS TrackingF IntlineS IntlineF :: + ::: XSM_FILE CompoDS ConcIso ; + ::: LINKED_LIST MicLib Flux TrackingS TrackingF ; + ::: SEQ_BINARY IntlineS IntlineF ; ; +*---- +* Get input parameters +*---- +REAL Bexit Power ; +INTEGER NbSteps ; +:: >>Bexit<< >>Power<< >>NbSteps<< ; +ECHO "Specific burnup power = " Power "kw/kg" ; +*---- +* Local modules, data structures +*---- +MODULE SHI: ASM: FLU: EDI: EVO: COMPO: DELETE: ; +LINKED_LIST PIJ EditDS ; +REAL keff ; +INTEGER Istep ; +REAL NbPt R TimeO ; +REAL Delt Timef := 1.0 0.0 ; +REAL Tfinal := Bexit Power / ; +*---- +* Burnup steps set in arithmetic progression +*---- +EVALUATE NbPt := NbSteps I_TO_R ; +EVALUATE R := Tfinal LN NbPt 1. - / EXP ; +ECHO "Burnup from " 0.0 " to " Tfinal " days using " NbSteps " steps" ; +* +*---- +* First flux calculation +* a) resonance self-shielding +* b) compute Pij matrix +* c) solve for the flux +* d) define edit parameters and edit +* e) clean up Pij +*---- +MicLib := SHI: MicLib TrackingS IntlineS :: EDIT 0 NOLJ ; +PIJ := ASM: MicLib TrackingF IntlineF :: EDIT 0 ; +Flux := FLU: PIJ MicLib TrackingF :: TYPE B B1 PNL ; +EditDS := EDI: Flux MicLib TrackingF :: + COND 0.625 MERG COMP MICR 1 Xe135 SAVE ON 'fuel' ; +PIJ := DELETE: PIJ ; +*---------- +* Burnup loop +*---------- +EVALUATE Istep := 0 ; +WHILE Istep NbSteps < DO + ECHO " Burning for " Delt " until " Timef ; + IF Timef 0.0 = THEN + ConcIso MicLib := EVO: MicLib Flux TrackingF :: + DEPL <<Delt>> DAY POWR <<Power>> ; + ELSE + ConcIso MicLib := EVO: ConcIso MicLib Flux TrackingF :: + DEPL <<Delt>> DAY POWR <<Power>> ; + ENDIF ; + CompoDS := COMPO: CompoDS EditDS ConcIso MicLib :: + EDIT 3 + STEP UP * + SET <<Timef>> DAY + ; +*---- +* Flux calculation after each burnup steps +* a) resonance self-shielding +* b) compute Pij matrix +* c) solve for the flux +* d) edit +* e) clean up Pij +*---- + MicLib := SHI: MicLib TrackingS IntlineS :: EDIT 0 NOLJ ; + PIJ := ASM: MicLib TrackingF IntlineF :: EDIT 0 ; + Flux := FLU: Flux PIJ MicLib TrackingF :: TYPE B B1 PNL ; + EditDS := EDI: EditDS Flux MicLib TrackingF :: SAVE ON 'fuel' ; + PIJ := DELETE: PIJ ; +*---- +* Select next burnup interval +*---- + EVALUATE Timef := Timef Delt + ; + EVALUATE Delt := Timef R 1. - * ; + EVALUATE Istep := Istep 1 + ; +ENDWHILE ; +*---- +* Clean up +*---- +EditDS := DELETE: EditDS ; +QUIT . diff --git a/Dragon/data/tdraglib_proc/TCDR01.c2m b/Dragon/data/tdraglib_proc/TCDR01.c2m new file mode 100644 index 0000000..bf4081f --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR01.c2m @@ -0,0 +1,115 @@ +*---- +* TEST CASE TCDR01 +* MOSTELLER BENCHMARK: 1-D ANNULAR CELL AND 2-D CARTESIAN CELL +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: R. Mosteller et al. Nucl. Sci. Eng. 107, 265 (1991) +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + MOSTELA MOSTELC DISCR LIBRARY CP CALC OUT ; +SEQ_BINARY + TRKSPC ; +MODULE + LIB: GEO: SYBILT: EXCELT: SHI: ASM: FLU: EDI: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 600.0 O16 = O16 4.61309E-2 + U235 = U235 1.66078E-4 1 + U238 = U238 2.28994E-2 1 + MIX 2 600.0 + Zr91 = Zr0 3.83243E-2 + MIX 3 600.0 + H1H2O = H1_H2O 4.42326E-2 O16H2O = O16 2.21163E-2 + B10 = B10 2.02222E-6 B11 = B11 8.19107E-6 + ; +*---- +* Geometry MOSTELA : annular 3 region geometry +* MOSTELC : Cartesian 3 region geometry +*---- +MOSTELA := GEO: :: TUBE 3 + R+ REFL RADIUS 0.0 0.39306 0.45802 0.71206 SPLITR 2 1 1 + MIX 1 2 3 ; +MOSTELC := GEO: :: CARCEL 2 + X- REFL X+ REFL MESHX 0.0 1.26209 + Y- REFL Y+ REFL MESHY 0.0 1.26209 + RADIUS 0.0 0.39306 0.45802 SPLITR 2 1 + MIX 1 2 3 ; +*---- +* Case 1 -- annular +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for K no leakage +*---- +DISCR := SYBILT: MOSTELA :: + TITLE 'TCDR01: MOSTELLER BENCHMARK (SYBIL / SYBIL)' + MAXR 4 QUA1 5 ; +LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ; +CP := ASM: LIBRARY DISCR ; +CALC := FLU: CP LIBRARY DISCR :: + TYPE K ; +assertS CALC :: 'K-EFFECTIVE' 1 0.8237322 ; +OUT := EDI: LIBRARY DISCR CALC :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ; +DISCR CP := DELETE: DISCR CP ; +*---- +* Case 2 -- Cartesian +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for K no leakage +*---- +DISCR := SYBILT: MOSTELC :: + TITLE 'TCDR01: MOSTELLER BENCHMARK (SYBIL / SYBIL)' + MAXR 4 QUA1 5 QUA2 6 5 ; +LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ; +CP := ASM: LIBRARY DISCR ; +CALC := FLU: CALC CP LIBRARY DISCR :: + TYPE K ; +assertS CALC :: 'K-EFFECTIVE' 1 0.8239249 ; +OUT := EDI: OUT LIBRARY DISCR CALC :: + EDIT 1 MERG MIX 1 2 3 COND 4.0 STAT ALL REFE 1 ; +CALC DISCR CP := DELETE: CALC DISCR CP ; +*---- +* Case 3 -- annular +* Self-Shielding calculation EXCEL-ISO +* Transport calculation EXCEL-ISO +* Flux calculation for K no leakage +*---- +DISCR TRKSPC := EXCELT: MOSTELC :: + TITLE 'TCDR01: MOSTELLER BENCHMARK (EXCELL)' + MAXR 4 TRAK TISO 12 20.0 ; +LIBRARY := SHI: LIBRARY DISCR TRKSPC :: EDIT 0 NOLJ ; +CP := ASM: LIBRARY DISCR TRKSPC ; +CALC := FLU: CP LIBRARY DISCR :: + TYPE K ; +assertS CALC :: 'K-EFFECTIVE' 1 0.8247744 ; +OUT := EDI: OUT LIBRARY DISCR CALC :: + EDIT 1 MERG MIX 1 2 3 COND 4.0 STAT ALL REFE 1 ; +DISCR TRKSPC CP := DELETE: DISCR TRKSPC CP ; +*---- +* Case 4 -- Cartesian +* Self-Shielding calculation EXCEL-SPC +* Transport calculation EXCEL-SPC +* Flux calculation for K no leakage +*---- +DISCR TRKSPC := EXCELT: MOSTELC :: + TITLE 'TCDR01: MOSTELLER BENCHMARK (EXCELL)' + MAXR 4 TRAK TSPC 12 20.0 ; +LIBRARY := SHI: LIBRARY DISCR TRKSPC :: EDIT 0 NOLJ ; +CP := ASM: LIBRARY DISCR TRKSPC ; +CALC := FLU: CALC CP LIBRARY DISCR :: + TYPE K ; +assertS CALC :: 'K-EFFECTIVE' 1 0.8256814 ; +OUT := EDI: OUT LIBRARY DISCR CALC :: + EDIT 1 MERG MIX 1 2 3 COND 4.0 STAT ALL REFE 1 ; +DISCR TRKSPC CP := DELETE: DISCR TRKSPC CP ; +ECHO "test TCDR01 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR02.c2m b/Dragon/data/tdraglib_proc/TCDR02.c2m new file mode 100644 index 0000000..1c6bf0a --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR02.c2m @@ -0,0 +1,162 @@ +*---- +* TEST CASE TRDR02 +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + ASSMB DISCR LIBRARY CP CALC OUT DATABASE ISOT SPHGEOM MTRACK ; +SEQ_ASCII + res ; +MODULE + LIB: GEO: SYBILT: BIVACT: SHI: ASM: FLU: EDI: COMPO: SPH: DELETE: + END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 + H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2 + B10 = B10 4.71444E-6 B11 = B11 1.90959E-5 + MIX 2 579.9 O16 = O16 3.06711E-4 + Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4 + Zr91 = Zr0 4.18621E-2 + MIX 3 579.9 + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + MoNat = Mo95 4.03755E-6 + Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr0 8.92427E-4 + MIX 4 579.9 O16 = O16 2.87335E-4 + Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4 + Zr91 = Zr0 3.92175E-2 + MIX 5 579.9 + H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2 + Cr52 = Cr52 2.24991E-5 + Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5 + MoNat = Mo95 1.89281E-6 + Al27 = Al27 1.10277E-6 + Mn55 = Mn55 1.94976E-7 + B10 = B10 4.66484E-6 B11 = B11 1.88950E-5 + Zr91 = Zr0 4.18372E-4 + MIX 6 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + Cr52 = Cr52 2.11122E-5 + Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5 + MoNat = Mo95 1.77614E-6 + Al27 = Al27 1.03479E-6 + Mn55 = Mn55 1.82957E-7 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + Zr91 = Zr0 3.92583E-4 + MIX 7 579.9 + H1H2O = H1_H2O 4.72020E-2 O16H2O = O16 2.36010E-2 + Cr52 = Cr52 1.96591E-5 + Fe56 = Fe56 1.82630E-5 Ni58 = Ni58 4.65011E-5 + MoNat = Mo95 1.65389E-6 + Al27 = Al27 9.63569E-7 + Mn55 = Mn55 1.70365E-7 + B10 = B10 4.67110E-6 B10 = B10 1.89203E-5 + Zr91 = Zr0 3.65562E-4 + MIX 8 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + ; +*---- +* Geometry ASSMB : a 17 X 17 normal PWR assembly +* contains C1 : cell without fuel +* C2 : normal fuel cell +* C3 : peripheral cell +* C4 : corner cell +*---- +ASSMB := GEO: :: CAR2D 9 9 + X- DIAG X+ REFL Y- SYME Y+ DIAG + CELL C1 C2 C2 C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C2 C2 C2 C3 + C2 C2 C2 C2 C2 C2 C3 + C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C3 + C1 C2 C2 C3 + C2 C2 C3 + C2 C3 + C4 + MERGE 1 3 12 11 12 12 11 12 15 + 4 6 5 6 6 5 6 8 + 13 5 6 6 5 6 8 + 2 5 5 10 5 8 + 13 5 5 6 8 + 2 5 7 8 + 13 7 8 + 14 8 + 9 + TURN H H B H H B H H A + H G G H G G H A + A E E F E E A + H H F H H A + H E G H A + H H A A + H A A + A A + A + ::: C1 := GEO: CARCEL 2 + MESHX 0.0 1.26472 MESHY 0.0 1.26472 + RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ; + ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ; + ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ; + ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ; + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for B1 homogeneous leakage +*---- +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; +LIBRARY := SHI: LIBRARY DISCR :: LJ ; +CP := ASM: LIBRARY DISCR :: ARM ; +CALC := FLU: CP LIBRARY DISCR :: + TYPE B B1 PNL ; +assertS CALC :: 'K-INFINITY' 1 1.256444 ; + +*---- +* Editing using SPH model for transport-diffusion with Selengut +* normalization +*---- +OUT := EDI: LIBRARY DISCR CALC ASSMB :: + EDIT 4 UPS SAVE MICR RES MERGE CELL COND 4.0 JOUT + ; +SPHGEOM := OUT :: STEP UP 'MACRO-GEOM' ; +MTRACK := BIVACT: SPHGEOM :: PRIM 1 2 EDIT 0 ; +OUT := SPH: OUT MTRACK :: EDIT 2 SELE_ALB ; +SPHGEOM MTRACK := DELETE: SPHGEOM MTRACK ; + +DATABASE := COMPO: :: + EDIT 5 + COMM 'Multi-parameter reactor database' ENDC + INIT + ; +DATABASE := COMPO: DATABASE OUT :: + EDIT 3 + ; +res := DATABASE ; + +ISOT := DATABASE :: STEP UP default + STEP UP MIXTURES STEP AT 5 + STEP UP CALCULATIONS STEP AT 1 + STEP UP ISOTOPESLIST STEP AT 1 ; +assertS ISOT :: 'NWT0' 1 0.5218698 ; +assertS ISOT :: 'NWT0' 2 0.1277670 ; + +ECHO "test TCDR02 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR02_void.c2m b/Dragon/data/tdraglib_proc/TCDR02_void.c2m new file mode 100644 index 0000000..b8a5e79 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR02_void.c2m @@ -0,0 +1,140 @@ +*---- +* TEST CASE TRDR02_void +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* VOID BOUNDARY CONDITION +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + ASSMB DISCR LIBRARY CP CALC ; +MODULE + LIB: GEO: SYBILT: SHI: ASM: FLU: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 + H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2 + B10 = B10 4.71444E-6 B11 = B11 1.90959E-5 + MIX 2 579.9 O16 = O16 3.06711E-4 + Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4 + Zr91 = Zr0 4.18621E-2 + MIX 3 579.9 + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + MoNat = Mo95 4.03755E-6 + Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr0 8.92427E-4 + MIX 4 579.9 O16 = O16 2.87335E-4 + Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4 + Zr91 = Zr0 3.92175E-2 + MIX 5 579.9 + H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2 + Cr52 = Cr52 2.24991E-5 + Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5 + MoNat = Mo95 1.89281E-6 + Al27 = Al27 1.10277E-6 + Mn55 = Mn55 1.94976E-7 + B10 = B10 4.66484E-6 B11 = B11 1.88950E-5 + Zr91 = Zr0 4.18372E-4 + MIX 6 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + Cr52 = Cr52 2.11122E-5 + Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5 + MoNat = Mo95 1.77614E-6 + Al27 = Al27 1.03479E-6 + Mn55 = Mn55 1.82957E-7 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + Zr91 = Zr0 3.92583E-4 + MIX 7 579.9 + H1H2O = H1_H2O 4.72020E-2 O16H2O = O16 2.36010E-2 + Cr52 = Cr52 1.96591E-5 + Fe56 = Fe56 1.82630E-5 Ni58 = Ni58 4.65011E-5 + MoNat = Mo95 1.65389E-6 + Al27 = Al27 9.63569E-7 + Mn55 = Mn55 1.70365E-7 + B10 = B10 4.67110E-6 B10 = B10 1.89203E-5 + Zr91 = Zr0 3.65562E-4 + MIX 8 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + ; +*---- +* Geometry ASSMB : a 17 X 17 normal PWR assembly +* contains C1 : cell without fuel +* C2 : normal fuel cell +* C3 : peripheral cell +* C4 : corner cell +*---- +ASSMB := GEO: :: CAR2D 9 9 + X- DIAG X+ REFL Y- SYME Y+ DIAG + CELL C1 C2 C2 C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C2 C2 C2 C3 + C2 C2 C2 C2 C2 C2 C3 + C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C3 + C1 C2 C2 C3 + C2 C2 C3 + C2 C3 + C4 + MERGE 1 3 12 11 12 12 11 12 15 + 4 6 5 6 6 5 6 8 + 13 5 6 6 5 6 8 + 2 5 5 10 5 8 + 13 5 5 6 8 + 2 5 7 8 + 13 7 8 + 14 8 + 9 + TURN H H B H H B H H A + H G G H G G H A + A E E F E E A + H H F H H A + H E G H A + H H A A + H A A + A A + A + ::: C1 := GEO: CARCEL 2 + MESHX 0.0 1.26472 MESHY 0.0 1.26472 + RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ; + ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ; + ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ; + ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ; + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for B1 homogeneous leakage +*---- +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; +LIBRARY := SHI: LIBRARY DISCR :: LJ ; +DISCR := DELETE: DISCR ; + +ASSMB := GEO: ASSMB :: + X- DIAG X+ VOID Y- SYME Y+ DIAG +; +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; + +CP := ASM: LIBRARY DISCR :: ARM ; +CALC := FLU: CP LIBRARY DISCR :: + EDIT 1 TYPE K P1 HETE BUCK -0.017069 ; +assertS CALC :: 'K-EFFECTIVE' 1 1.000043 ; + +ECHO "test TCDR02_void completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR03.c2m b/Dragon/data/tdraglib_proc/TCDR03.c2m new file mode 100644 index 0000000..380eb42 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR03.c2m @@ -0,0 +1,153 @@ +*---- +* TEST CASE TCDR03 +* MULTICELL HEXAGONAL ASSEMBLY WITH POISON +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + ASSMBH DISCR LIBRARY CP CALC OUT DATABASE ISOT SPHGEOM MTRACK ; +SEQ_ASCII + res ; +MODULE + GEO: SYBILT: BIVACT: LIB: SHI: ASM: FLU: EDI: COMPO: SPH: DELETE: + END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 11 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 + H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2 + B10 = B10 4.71444E-6 B11 = B11 1.90959E-5 + MIX 2 579.9 O16 = O16 3.06711E-4 + Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4 + Zr90 = Zr90 2.15171E-2 Zr91 = Zr91 4.68856E-3 + Zr92 = Zr92 7.15842E-3 Zr94 = Zr94 7.32587E-3 + MIX 3 579.9 + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Nb93 = Nb93 7.07082E-6 Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + Co59 = Co59 1.07714E-7 Mo95 = Mo95 4.03755E-6 + Ti48 = Ti0 2.33305E-6 Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr91 8.92427E-4 + MIX 4 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + MIX 5 579.9 In115 = In115 7.57464E-3 + Cd113 = Cd113 2.62493E-3 + Ag109 = Ag109 4.49188E-2 + MIX 6 579.9 Cr52 = Cr52 1.52702E-2 + Fe56 = Fe56 5.57670E-2 Ni58 = Ni58 7.51418E-3 + Mn55 = Mn55 8.02943E-4 + MIX 7 579.9 + H1H2O = H1_H2O 3.06466E-2 O16H2O = O16 1.53233E-2 + Fe56 = Fe56 5.27485E-5 Cr52 = Cr52 2.69769E-5 + B10 = B10 3.03092E-6 B11 = B11 1.22768E-5 + Zr90 = Zr90 7.68841E-3 Zr91 = Zr91 1.67530E-3 + Zr92 = Zr92 2.55782E-3 Zr94 = Zr94 2.61765E-3 + MIX 8 579.9 + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Nb93 = Nb93 7.07082E-6 Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + Co59 = Co59 1.07714E-7 Mo95 = Mo95 4.03755E-6 + Ti48 = Ti0 2.33305E-6 Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr91 8.92427E-4 + MIX 9 579.9 O16 = O16 2.87335E-4 + Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4 + Zr90 = Zr90 2.01578E-2 Zr91 = Zr91 4.39236E-3 + Zr92 = Zr92 6.70619E-3 Zr94 = Zr94 6.86306E-3 + MIX 10 579.9 + H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2 + Nb93 = Nb93 3.31482E-6 Cr52 = Cr52 2.24991E-5 + Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5 + Co59 = Co59 5.04968E-8 Mo95 = Mo95 1.89281E-6 + Ti48 = Ti0 1.09374E-6 Al27 = Al27 1.10277E-6 + Mn55 = Mn55 1.94976E-7 B10 = B10 4.66484E-6 + B11 = B11 1.88950E-5 Zr91 = Zr91 4.18372E-4 + MIX 11 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + Nb93 = Nb93 3.11049E-6 Cr52 = Cr52 2.11122E-5 + Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5 + Co59 = Co59 4.73842E-8 Mo95 = Mo95 1.77614E-6 + Ti48 = Ti0 1.02632E-6 Al27 = Al27 1.03479E-6 + Mn55 = Mn55 1.82957E-7 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + Zr91 = Zr91 3.92583E-4 + ; +*---- +* Geometry ASSMBH : hexagonal assembly with poison +* contains C1 : cell without fuel +* C2 : poison cell +* C3 : normal fuel cell +* C4 : peripheral cell +*---- +ASSMBH := GEO: :: HEX 36 + HBC S30 REFL + CELL C1 C3 C3 C3 C3 C3 C2 C3 C3 C3 C2 C3 C3 C3 C3 C3 C3 C2 + C3 C3 C2 C3 C3 C3 C3 C3 C3 C3 C3 C3 C4 C4 C4 C4 C4 C4 + TURN A A A A A A A A B D A I A C F J B A + F A A E E A A E A A A A A A A A A A + MERGE 1 2 3 4 5 4 6 7 8 7 9 8 10 7 7 4 7 11 + 12 13 14 15 12 16 17 12 16 18 18 19 20 21 21 22 22 23 + ::: C1 := GEO: HEXCEL 2 + SIDE 0.707297 RADIUS 0.0 0.412282 0.475917 + MIX 1 2 3 ; + ::: C2 := GEO: HEXCEL 5 + SIDE 0.707297 RADIUS 0.0 0.25057 0.354359 0.436 0.486 0.6125 + MIX 5 5 5 6 7 8 ; + ::: C3 := GEO: C1 MIX 4 9 10 ; + ::: C4 := GEO: C3 MIX 4 9 11 ; + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for B1 homogeneous leakage +* Editing using SPH model for transport-diffusion +*---- +DISCR := SYBILT: ASSMBH :: + TITLE 'TCDR03: MULTICELL HEXAGONAL ASSEMBLY WITH POISON' + MAXR 400 MAXZ 15000 QUA2 4 3 ; +LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ; +CP := ASM: LIBRARY DISCR ; +CALC := FLU: CP LIBRARY DISCR :: + TYPE B B1 PNL ; +assertS CALC :: 'K-INFINITY' 1 0.6800977 ; +OUT := EDI: LIBRARY DISCR CALC ASSMBH :: + EDIT 3 UPS SAVE MICR RES MERGE CELL COND 4.0 + ; +SPHGEOM := OUT :: STEP UP 'MACRO-GEOM' ; +SPHGEOM := GEO: SPHGEOM :: SPLITH 0 ; +MTRACK := BIVACT: SPHGEOM :: PRIM 1 2 ; +OUT := SPH: OUT MTRACK ; +SPHGEOM MTRACK := DELETE: SPHGEOM MTRACK ; + +DATABASE := COMPO: :: + EDIT 5 + COMM 'Multi-parameter reactor database' ENDC + INIT + ; +DATABASE := COMPO: DATABASE OUT :: + EDIT 3 + ; +res := DATABASE ; + +ISOT := DATABASE :: STEP UP default + STEP UP MIXTURES STEP AT 5 + STEP UP CALCULATIONS STEP AT 1 + STEP UP ISOTOPESLIST STEP AT 1 ; +assertS ISOT :: 'NWT0' 1 1.504198E+00 ; +assertS ISOT :: 'NWT0' 2 1.369299E-01 ; + +ECHO "test TCDR03 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR04.c2m b/Dragon/data/tdraglib_proc/TCDR04.c2m new file mode 100644 index 0000000..c98134b --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR04.c2m @@ -0,0 +1,158 @@ +*---- +* TEST CASE TCDR04 +* ANNULAR MOSTELLER BENCHMARK WITH BURNUP +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: R. Mosteller et al. Nucl. Sci. Eng. 107, 265 (1991) +* +*---- +* Define variables +*---- +INTEGER + istep := 1 ; +REAL + evobeg evoend ; +REAL + step2 step3 step4 step5 := + 1.0 27.1739 67.9348 135.8696 ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY MOSTELAS MOSTELA TRACKS TRACK SYS FLUX BURNUP + EDITION DATABASE ISOT ; +SEQ_ASCII + res ; +MODULE + GEO: SYBILT: LIB: SHI: ASM: FLU: EVO: EDI: COMPO: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + DEPL LIB: DRAGON FIL: DLIB_J2 + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 600.0 O16 = O16 4.61309E-2 + U235 = U235 1.66078E-4 1 + U238 = U238 2.28994E-2 1 + U236 = U236 0.0 1 + Pu239 = Pu239 0.0 1 + MIX 2 600.0 + Zr90 = Zr90 1.96987E-2 Zr91 = Zr91 4.29232E-3 + Zr92 = Zr92 6.55346E-3 Zr94 = Zr94 6.70675E-3 + MIX 3 600.0 + H1H2O = H1_H2O 4.42326E-2 O16H2O = O16 2.21163E-2 + B10 = B10 2.02223E-6 B11 = B11 8.19107E-6 + ; +*---- +* Geometry MOSTELAS : 3 regions annular cell for self-shielding +* MOSTELA : 4 regions annular cell for transport +*---- +MOSTELAS := GEO: :: TUBE 3 + R+ REFL RADIUS 0.0 0.39306 0.45802 0.71206 + MIX 1 2 3 ; +MOSTELA := GEO: MOSTELAS :: + SPLITR 2 1 1 ; +*---- +* Create the reactor database +*---- +DATABASE := COMPO: :: + EDIT 5 + COMM 'Multi-parameter reactor database' ENDC + PARA 'BURN' IRRA + PARA 'FLUB' FLUB + INIT + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for keff with imposed buckling +* using B1 homogeneous leakage model +*---- +TRACKS := SYBILT: MOSTELAS :: + TITLE 'TCDR04: MOSTELLER BENCHMARK WITH BURNUP' + EDIT 1 MAXR 3 ; +LIBRARY := SHI: LIBRARY TRACKS :: EDIT 0 NOLJ ; +TRACK := SYBILT: MOSTELA :: + TITLE 'TCDR04: MOSTELLER BENCHMARK WITH BURNUP' + EDIT 1 MAXR 4 ; +SYS := ASM: LIBRARY TRACK ; +FLUX := FLU: SYS LIBRARY TRACK :: + TYPE K B1 SIGS BUCK 0.2948E-2 ; +assertS FLUX :: 'K-INFINITY' 1 0.8247780 ; +EDITION := EDI: LIBRARY TRACK FLUX :: + EDIT 3 MICR 25 Th232 U233 U234 U235 U236 U238 Np237 Np238 Pu238 Pu239 + Pu240 Pu241 Pu242 Am241 Am242 Am243 Cm242 Cm243 Cm244 + Cm245 Cf252 I135 Xe135 Pm149 Sm149 + MERG COMP COND 4.0 SAVE ; +*---- +* Burnup loop: for first step BURNUP is created +* while for other steps it is modified +* two burnup per step: +* 1) get a first approximation of final composition followed +* by a transport calculation +* 2) use approximation for final flux distribution to get a +* better approximation for final composition +*---- +EVALUATE evoend := 0.0 ; +WHILE evoend step2 < DO + EVALUATE evobeg := evoend ; + EVALUATE evoend := step2 ; + IF istep 1 = THEN + BURNUP LIBRARY := EVO: LIBRARY FLUX TRACK :: + SAVE <<evobeg>> DAY POWR 36.8 + DEPL <<evobeg>> <<evoend>> DAY POWR 36.8 + SET <<evoend>> DAY ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK :: + SAVE <<evobeg>> DAY POWR 36.8 + NOEX DEPL <<evobeg>> <<evoend>> DAY POWR 36.8 + SET <<evoend>> DAY ; + ENDIF ; + LIBRARY := SHI: LIBRARY TRACKS :: EDIT 0 NOLJ ; + SYS := DELETE: SYS ; + SYS := ASM: LIBRARY TRACK ; + FLUX := FLU: FLUX SYS LIBRARY TRACK :: + TYPE K B1 SIGS BUCK 0.2948E-2 ; + + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK :: + SAVE <<evoend>> DAY POWR 36.8 + DEPL <<evobeg>> <<evoend>> DAY POWR 36.8 + SET <<evoend>> DAY ; + LIBRARY := SHI: LIBRARY TRACKS ; + SYS := DELETE: SYS ; + SYS := ASM: LIBRARY TRACK ; + FLUX := FLU: FLUX SYS LIBRARY TRACK :: + TYPE K B1 SIGS BUCK 0.2948E-2 ; + + EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + EDIT 3 SAVE ; + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK :: + SAVE <<evoend>> DAY POWR 36.8 ; + DATABASE := COMPO: DATABASE EDITION BURNUP :: + EDIT 3 + SET <<evoend>> DAY + ; + EVALUATE step2 step3 step4 step5 := + step3 step4 step5 step2 ; + EVALUATE istep := istep 1 + ; +ENDWHILE ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.7645292 ; + +res := DATABASE ; + +*---- +* Set isotope *MAC*RES on DATABASE +*---- +ISOT := DATABASE :: STEP UP default + STEP UP MIXTURES STEP AT 1 + STEP UP CALCULATIONS STEP AT 2 + STEP UP ISOTOPESLIST STEP AT 26 ; +assertS ISOT :: 'NWT0' 1 37.66850 ; +assertS ISOT :: 'NWT0' 2 16.65745 ; + +ECHO "test TCDR04 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR05.c2m b/Dragon/data/tdraglib_proc/TCDR05.c2m new file mode 100644 index 0000000..7f8d6cb --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR05.c2m @@ -0,0 +1,187 @@ +*---- +* TEST CASE TCDR05 +* CANDU-6 ANNULAR CELL +* 172-GROUP JEF2.2 DRAGLIB +* POWER (KW) = 615.00000 +* BURN POWER (KW/KG) = 31.97130 +* URANIUM MASS = 19.23600 +* UO2 REAL DENSITY = 10.59300 +* UO2 EFF DENSITY = 10.43750 +* UO2 TEMPERATURE = 941.28998 +* ENRICHMENT = 0.71140 +* COOLANT D2 AT % = 99.222 +* MODERATOR D2 AT % = 99.911 +* NUMBER OF DAYS = 50 +* +*---- +* Define variables and initialize +* Burnup paremeters +* a) Power +* = 31.9713 kw/kg for 0.0 to 300.0 days +* b) Burnup time interval Delt +* = 1 day for 0 to 1 day +* = 4 days for 1 to 5 days +* = 5 days for 5 to 10 days +* = 10 days for 10 to 50 days +* = 20 days for 50 to 150 days +* = 50 days for 150 to 300 days +* c) Days with burnup interval changes +* = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days +* d) Burnup control time variables Timei, Timef +* Timei = initial time +* Timef = final time +*---- +REAL + Power Delt Timec Timei Timef := + 31.9713 1.0 1.0 0.0 0.0 ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION + DATABASE ISOT ; +SEQ_BINARY + INTLINS INTLINF ; +SEQ_ASCII + database ; +MODULE + GEO: EXCELT: SHI: ASM: FLU: EVO: EDI: COMPO: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +PROCEDURE TCDR05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCDR05Lib :: <<iedit>> ; +*---- +* Geometry CANDU6S : 13 regions annular cluster for self-shielding +* CANDU6F : 31 regions annular cluster for transport +*---- +CANDU6S := GEO: :: TUBE 5 + R+ REFL RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 16.12171 + MIX 1 2 3 4 5 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ; + ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ; + ; +CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10 + ::: ROD1 := GEO: ROD1 SPLITR 2 1 ; + ::: ROD2 := GEO: ROD2 SPLITR 2 1 ; + ::: ROD3 := GEO: ROD3 SPLITR 2 1 ; + ::: ROD4 := GEO: ROD4 SPLITR 2 1 ; + ; +*---- +* Create the reactor database +*---- +DATABASE := COMPO: :: + EDIT 5 + STEP UP 'moderator' + COMM 'Multi-parameter reactor database for moderator' ENDC + INIT + STEP UP 'fuel' + COMM 'Multi-parameter reactor database for fuel' ENDC + PARA 'BURN' IRRA + PARA 'FLUB' FLUB + INIT + ; +*---- +* Self-Shielding calculation EXCEL +* Transport calculation EXCEL +* Flux calculation for keff +*---- +VOLMATS INTLINS := EXCELT: CANDU6S :: + TITLE 'TCDR05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + EDIT 0 MAXR 13 TRAK TISO 5 10.0 SYMM 12 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; +VOLMATF INTLINF := EXCELT: CANDU6F :: + TITLE 'TCDR05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-INFINITY' 1 1.117541 ; +EDITION := EDI: LIBRARY VOLMATF FLUX :: + COND 4.0 MERGE MIX 0 0 0 0 1 0 0 0 0 0 MICR RES SAVE ON 'moderator' ; +DATABASE := COMPO: DATABASE EDITION :: + EDIT 3 + STEP UP * + ; +EDITION := EDI: EDITION LIBRARY VOLMATF FLUX :: + COND 4.0 MERGE COMP MICR 1 Xe135 SAVE ON 'fuel' ; +*---- +* Burnup loop: for first step BURNUP is created +* while for other steps it is modified +*---- +WHILE Timei Timec < DO + EVALUATE Timef := Timei Delt + ; + IF Timei 0.0 = THEN + BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF :: + DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ENDIF ; + LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; + PIJ := DELETE: PIJ ; + PIJ := ASM: LIBRARY VOLMATF INTLINF ; + FLUX := FLU: FLUX PIJ LIBRARY VOLMATF :: + TYPE K ; + EDITION := EDI: EDITION LIBRARY VOLMATF FLUX :: + SAVE ON 'fuel' ; + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + SAVE <<Timef>> DAY POWR <<Power>> ; + DATABASE := COMPO: DATABASE EDITION BURNUP LIBRARY :: + EDIT 3 + STEP UP * + SET <<Timef>> DAY + ; +*---- +* change delta t for burnup and final time if required +*---- + IF Timef Timec = THEN + IF Timec 150.0 = THEN + EVALUATE Delt Timec := 50.0 300.0 ; + ENDIF ; + IF Timec 50.0 = THEN + EVALUATE Delt Timec := 20.0 150.0 ; + ENDIF ; + 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 ; +*---- +* Export and access the database +*---- +database := DATABASE ; + +ISOT := DATABASE :: STEP UP fuel + STEP UP MIXTURES STEP AT 1 + STEP UP CALCULATIONS STEP AT 5 + STEP UP ISOTOPESLIST STEP AT 1 ; +assertS ISOT :: 'NWT0' 1 9.761326E0 ; +assertS ISOT :: 'NWT0' 2 2.651622E1 ; +ISOT := DELETE: ISOT ; + +ISOT := DATABASE :: STEP UP moderator + STEP UP MIXTURES STEP AT 1 + STEP UP CALCULATIONS STEP AT 1 + STEP UP ISOTOPESLIST STEP AT 1 ; +assertS ISOT :: 'NWT0' 1 9.095968E-2 ; +assertS ISOT :: 'NWT0' 2 3.246063E-1 ; + +INTLINF INTLINS := DELETE: INTLINF INTLINS ; +ECHO "test TCDR05 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR05Lib.c2m b/Dragon/data/tdraglib_proc/TCDR05Lib.c2m new file mode 100644 index 0000000..bd132b3 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR05Lib.c2m @@ -0,0 +1,65 @@ +*---- +* Procedure TCDR05Lib +* Create Library for test CASE TCDR05 +* Calling : +* LIBRARY := TCDR05Lib :: iedit ; +* with : +* LIBRARY = Linked list containing the result of LIB: for +* TCDR05 +* iprint = print level for LIB: module +*---- +* Define PARAMETERS,STRUCTURES and MODULES used +*---- +PARAMETER LIBRARY :: ::: LINKED_LIST LIBRARY ; ; +MODULE LIB: END: ; +*---- +* Define and read LIB: EDIT option +INTEGER iedit ; +:: >>iedit<< ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + EDIT <<iedit>> + NMIX 17 CTRA WIMS + DEPL LIB: DRAGON FIL: DLIB_J2 + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 560.66 0.81212 O16 = O16 7.99449E-1 + D2D2O = H2_D2O 1.99768E-1 H1H2O = H1_H2O 7.83774E-4 + MIX 2 560.66 6.57 Nb93 = Nb93 2.50000 + B10 = B10 4.15800E-5 B11 = B11 1.68420E-4 + Zr90 = Zr90 5.01150E+1 Zr91 = Zr91 1.09200E+1 + Zr92 = Zr92 1.66725E+1 Zr94 = Zr94 1.70625E+1 + MIX 3 345.66 0.0014 He4 = He4 1.00000E+2 + MIX 4 345.66 6.44 Fe56 = Fe56 1.60000E-1 + Ni58 = Ni58 6.00000E-2 Cr52 = Cr52 1.10000E-1 + B10 = B10 6.13800E-5 B11 = B11 2.48620E-4 + Zr90 = Zr90 5.12509E+1 Zr91 = Zr91 1.11675E+1 + Zr92 = Zr92 1.70504E+1 Zr94 = Zr94 1.74493E+1 + MIX 5 345.66 1.082885 O16 = O16 7.98895E-1 + D2D2O = H2_D2O 2.01016E-1 H1H2O = H1_H2O 8.96000E-5 + MIX 6 941.29 10.4375010 O16 = O16 1.18473E+1 + Xe135 = Xe135 0.0 + U235 = U235 6.27118E-1 1 + U238 = U238 8.75256E+1 1 + U236 = U236 0.0 1 + Pu239 = Pu239 0.0 1 + MIX 7 COMB 6 1.0 + MIX 8 COMB 6 1.0 + MIX 9 COMB 6 1.0 + MIX 10 560.66 6.44 Fe56 = Fe56 1.60000E-1 + Ni58 = Ni58 6.00000E-2 Cr52 = Cr52 1.10000E-1 + B10 = B10 6.13800E-5 B11 = B11 2.48620E-4 + Zr90 = Zr90 5.12509E+1 Zr91 = Zr91 1.11675E+1 + Zr92 = Zr92 1.70504E+1 Zr94 = Zr94 1.74493E+1 + MIX 11 COMB 10 1.0 + MIX 12 COMB 10 1.0 + MIX 13 COMB 10 1.0 + MIX 14 COMB 1 1.0 + MIX 15 COMB 1 1.0 + MIX 16 COMB 1 1.0 + MIX 17 COMB 1 1.0 + ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR06.c2m b/Dragon/data/tdraglib_proc/TCDR06.c2m new file mode 100644 index 0000000..bae3d04 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR06.c2m @@ -0,0 +1,151 @@ +*---- +* TEST CASE TCDR06 +* CANDU-6 CARTESIAN CELL +* 172-GROUP JEF2.2 DRAGLIB +* STAINLESS STELL RODS IN 3D SUPERCELL +* +* REF: R. Roy et al. Ann. Nucl. Energy 21, 115 (1994) +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY LIBRARY2 CANDU6F CANDU6S TRACK SYS FLUX EDITION BCO BCI + SSRODS ; +SEQ_BINARY + INTLIN ; +MODULE + GEO: EXCELT: LIB: SHI: ASM: FLU: EDI: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +PROCEDURE TCDR05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCDR05Lib :: <<iedit>> ; +*---- +* CELL CALCULATION +* Geometry CANDU6S : 14 regions Cartesian cluster for self-shielding +* CANDU6F : 32 regions Cartesian cluster for transport +* BCO : 48 regions 3D Cartesian geometry +* BCI : 48 regions 3D Cartesian geometry +*---- +CANDU6S := GEO: :: CARCEL 5 + X+ REFL X- REFL MESHX -14.2875 14.2875 + Y+ REFL Y- REFL MESHY -14.2875 14.2875 + RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00 + MIX 1 2 3 4 5 5 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ; + ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ; + ; +CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10 + ::: ROD1 := GEO: ROD1 SPLITR 2 1 ; + ::: ROD2 := GEO: ROD2 SPLITR 2 1 ; + ::: ROD3 := GEO: ROD3 SPLITR 2 1 ; + ::: ROD4 := GEO: ROD4 SPLITR 2 1 ; + ; +*---- +* Self-Shielding calculation EXCEL +* Transport calculation EXCEL +* Flux calculation for keff +*---- +TRACK INTLIN := EXCELT: CANDU6S :: + TITLE 'TCDR06: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 14 TRAK TISO 29 20.0 SYMM 4 ; +LIBRARY := SHI: LIBRARY TRACK INTLIN :: EDIT 0 NOLJ ; +TRACK INTLIN := DELETE: TRACK INTLIN ; +TRACK INTLIN := EXCELT: CANDU6F :: + TITLE 'TCDR06: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 32 TRAK TISO 29 20.0 SYMM 4 ; +SYS := ASM: LIBRARY TRACK INTLIN :: + EDIT 0 ; +FLUX := FLU: SYS LIBRARY TRACK :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.117081 ; +*---- +* Microscopic cross sections from WIMSLIB for stainless steel +* MIX 1,2,3 from EDI: : fuel, structure material and moderator +* MIX 4 stainless steel rods +*---- +EDITION := EDI: LIBRARY TRACK FLUX :: + EDIT 0 MERGE MIX 1 2 2 2 3 1 1 1 1 1 SAVE ON SSRODS ; +SSRODS := EDITION :: STEP UP SSRODS STEP UP MACROLIB ; +LIBRARY2 := LIB: SSRODS :: + EDIT 0 NMIX 4 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 4 345.66 Fe56 = Fe56 6.19027E-2 + Cr52 = Cr52 1.56659E-3 Ni58 = Ni58 6.83337E-3 + Si29 = Si0 7.79072E-4 C12 = C0 1.46552E-4 + Mn55 = Mn55 1.25431E-3 + ; +EDITION TRACK INTLIN SYS FLUX := DELETE: + EDITION TRACK INTLIN SYS FLUX ; +*---- +* SUPERCELL CALCULATION +* Geometry BCO : 27 regions 3D Cartesian geometry with rods out +* BCI : 27 regions 3D Cartesian geometry with rods in +*---- +BCO := GEO: :: CAR3D 3 2 2 + X- REFL X+ SYME Y- REFL Y+ SYME Z- REFL Z+ SYME + CELL M MX MX MX FXY MXY M MX BX MX FXY BXY + TURN A A A F A A A A A F A A + ::: M := GEO: CAR3D 1 1 1 MIX 3 + MESHX 0.0 7.14375 + MESHY 0.0 7.14375 + MESHZ -8.25500 8.25500 SPLITZ 2 ; + ::: MX := GEO: M MESHX -7.14375 +7.14375 SPLITX 2 ; + ::: MXY := GEO: MX MESHY -7.14375 +7.14375 SPLITY 2 ; + ::: BX := GEO: CARCELY 2 1 MIX 3 3 3 + MESHX -7.14375 7.14375 SPLITX 2 + MESHY 0.0 7.14375 + MESHZ -8.25500 8.25500 SPLITZ 2 + RADIUS 0.0 3.5100 3.8100 ; + ::: BXY := GEO: BX MESHY -7.14375 +7.14375 SPLITY 2 ; + ::: FXY := GEO: CARCELZ 2 1 MIX 1 2 3 + MESHX -7.14375 7.14375 SPLITX 2 + MESHY -7.14375 7.14375 SPLITY 2 + MESHZ -8.25500 8.25500 SPLITZ 2 + RADIUS 0.0 5.16890 6.58750 ; + ; +BCI := GEO: BCO :: + ::: BX := GEO: BX MIX 3 4 3 ; + ::: BXY := GEO: BXY MIX 3 4 3 ; + ; +*---- +* Transport calculation EXCEL +* Flux calculation for keff +* Homogenized properties for rod out +*---- +TRACK INTLIN := EXCELT: BCO :: + EDIT 0 MAXR 40 TRAK TISO 2 1.0 ; +SYS := ASM: LIBRARY2 TRACK INTLIN :: + EDIT 0 ; +FLUX := FLU: SYS LIBRARY2 TRACK :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.127179 ; +EDITION := EDI: LIBRARY2 TRACK FLUX :: + EDIT 2 MERG COMP COND 4.0 SAVE ON 'NOBC' ; +SYS TRACK INTLIN := DELETE: SYS TRACK INTLIN ; +*---- +* Transport calculation EXCEL +* Flux calculation for keff +* Homogenized properties for rod in +*---- +TRACK INTLIN := EXCELT: BCI :: + EDIT 0 MAXR 40 TRAK TISO 2 1.0 ; +SYS := ASM: LIBRARY2 TRACK INTLIN :: + EDIT 0 ; +FLUX := FLU: FLUX SYS LIBRARY2 TRACK :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.041988 ; +EDITION := EDI: EDITION LIBRARY2 TRACK FLUX :: + EDIT 2 MERG COMP COND 4.0 STAT DELS REFE 'NOBC' ; +TRACK INTLIN SYS := DELETE: TRACK INTLIN SYS ; +ECHO "test TCDR06 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR07.c2m b/Dragon/data/tdraglib_proc/TCDR07.c2m new file mode 100644 index 0000000..be6360a --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR07.c2m @@ -0,0 +1,138 @@ +*---- +* TEST CASE TCDR07 +* CANDU-6 CARTESIAN CELL +* 172-GROUP JEF2.2 DRAGLIB +* TEST VARIOUS LEAKAGE OPTIONS +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY CANDU6S CANDU6T CANDU6SV CANDU6TV TRACK + SYS FLUX EDITION ; +MODULE + GEO: EXCELT: LIB: SHI: ASM: FLU: EDI: DELETE: END: ; +SEQ_BINARY + INTLIN ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +PROCEDURE TCDR05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCDR05Lib :: <<iedit>> ; +*---- +* Geometry CANDU6S : GEOMETRY FOR SELF-SHIELDING (NO VOID) +* CANDU6F : GEOMETRY FOR TRANSPORT (NO VOID) +* CANDU6FV: GEOMETRY FOR TRANSPORT (COOLANT VOID) +* CANDU6FV: GEOMETRY FOR TRANSPORT (COOLANT VOID) +*---- +CANDU6S := GEO: :: CARCEL 5 + X+ REFL X- REFL MESHX -14.2875 14.2875 + Y+ REFL Y- REFL MESHY -14.2875 14.2875 + RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 14.00 + MIX 1 2 3 4 5 5 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ; + ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ; + ; +CANDU6T := GEO: CANDU6S :: SPLITR 6 1 1 1 10 + ::: ROD1 := GEO: ROD1 SPLITR 2 1 ; + ::: ROD2 := GEO: ROD2 SPLITR 2 1 ; + ::: ROD3 := GEO: ROD3 SPLITR 2 1 ; + ::: ROD4 := GEO: ROD4 SPLITR 2 1 ; + ; +CANDU6SV := GEO: CANDU6S :: MIX 0 2 3 4 5 5 ; +CANDU6TV := GEO: CANDU6SV :: SPLITR 6 1 1 1 10 + ::: ROD1 := GEO: ROD1 SPLITR 2 1 ; + ::: ROD2 := GEO: ROD2 SPLITR 2 1 ; + ::: ROD3 := GEO: ROD3 SPLITR 2 1 ; + ::: ROD4 := GEO: ROD4 SPLITR 2 1 ; + ; +*---- +* CASE WITH NO VOID +* Self-Shielding calculation EXCEL +* Transport calculation EXCEL +* Flux TYPE K AND B WITH VARIOUS LEAKAGE OPTIONS +*---- +TRACK INTLIN := EXCELT: CANDU6S :: + TITLE 'TCDR07: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 14 TRAK TISO 7 20.0 SYMM 4 ; +LIBRARY := SHI: LIBRARY TRACK INTLIN :: EDIT 0 NOLJ ; +TRACK INTLIN := DELETE: TRACK INTLIN ; +TRACK INTLIN := EXCELT: CANDU6T :: + TITLE 'TCDR07: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 32 ANIS 2 TRAK TISO 7 20.0 SYMM 4 ; +SYS := ASM: LIBRARY TRACK INTLIN :: + EDIT 0 PIJK ; +FLUX := FLU: SYS LIBRARY TRACK :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.117062 ; +EDITION := EDI: FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX := FLU: FLUX SYS LIBRARY TRACK :: + TYPE B B1 PNL ; +assertS FLUX :: 'K-INFINITY' 1 1.108985 ; +EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX := FLU: FLUX SYS LIBRARY TRACK :: + TYPE B B1 TIBERE ; +assertS FLUX :: 'K-INFINITY' 1 1.108995 ; +EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX SYS := DELETE: FLUX SYS ; +SYS := ASM: LIBRARY TRACK INTLIN :: + EDIT 0 PIJ ECCO ; +FLUX := FLU: SYS LIBRARY TRACK :: + TYPE B B1 ECCO ; +assertS FLUX :: 'K-INFINITY' 1 1.108996 ; +EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX SYS TRACK INTLIN := DELETE: FLUX SYS TRACK INTLIN ; +*---- +* CASE WITH COOLANT VOIDED +* Self-Shielding calculation EXCEL +* Transport calculation EXCEL +* Flux TYPE K AND B WITH VARIOUS LEAKAGE OPTIONS +*---- +TRACK INTLIN := EXCELT: CANDU6SV :: + TITLE 'TCDR07: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 14 TRAK TISO 7 20.0 SYMM 4 ; +LIBRARY := SHI: LIBRARY TRACK INTLIN :: EDIT 0 NOLJ ; +TRACK INTLIN := DELETE: TRACK INTLIN ; +TRACK INTLIN := EXCELT: CANDU6TV :: + TITLE 'TCDR07: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 32 ANIS 2 TRAK TISO 7 20.0 SYMM 4 ; +SYS := ASM: LIBRARY TRACK INTLIN :: + EDIT 0 PIJK ; +FLUX := FLU: SYS LIBRARY TRACK :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.140569 ; +EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX := FLU: FLUX SYS LIBRARY TRACK :: + TYPE B B1 PNL ; +assertS FLUX :: 'K-INFINITY' 1 1.132941 ; +EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX := FLU: FLUX SYS LIBRARY TRACK :: + TYPE B B1 TIBERE ; +assertS FLUX :: 'K-INFINITY' 1 1.133257 ; +EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX SYS := DELETE: FLUX SYS ; +SYS := ASM: LIBRARY TRACK INTLIN :: + EDIT 0 PIJ ECCO ; +FLUX := FLU: SYS LIBRARY TRACK :: + TYPE B B1 ECCO ; +assertS FLUX :: 'K-INFINITY' 1 1.133257 ; +EDITION := EDI: EDITION FLUX LIBRARY TRACK :: + EDIT 3 SAVE COND 4.0 TAKE REGI 1 4 7 10 16 24 ; +FLUX SYS TRACK INTLIN := DELETE: FLUX SYS TRACK INTLIN ; +ECHO "test TCDR07 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR08.c2m b/Dragon/data/tdraglib_proc/TCDR08.c2m new file mode 100644 index 0000000..311ac4c --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR08.c2m @@ -0,0 +1,124 @@ +*---- +* TEST CASE TCDR08 +* HOMOGENEOUS DEPLETION CASE +* 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 LIB: DRAGON FIL: DLIB_J2 + 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 'TCDR08: 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 + IF Power 0.0 = THEN + ECHO "Cooling stage...." ; + ENDIF ; + 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.217538 ; +ECHO "test TCDR08 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR09.c2m b/Dragon/data/tdraglib_proc/TCDR09.c2m new file mode 100644 index 0000000..9af763b --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR09.c2m @@ -0,0 +1,72 @@ +*---- +* TEST CASE TCDR09 +* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL +* 172-GROUP JEF2.2 DRAGLIB +* REFLECTIVE AND VOID BC +* +* REF: None +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + MOSTELA MOSTELV VOLMAT LIBRARY PIJ FLUX OUT ; +MODULE + LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 600.0 O16 = O16 4.61309E-2 + U235 = U235 1.66078E-4 1 + U238 = U238 2.28994E-2 1 + MIX 2 600.0 + Zr91 = Zr91 3.83243E-2 + MIX 3 600.0 + H1H2O = H1_H2O 4.42326E-2 O16H2O = O16 2.21163E-2 + B10 = B10 2.02223E-6 B11 = B11 8.19107E-6 + ; +*---- +* Geometry MOSTELA : Annular cell with reflective BC +* MOSTELV : Annular cell with void BC +*---- +MOSTELA := GEO: :: TUBE 3 + RADIUS 0.0 0.39306 0.45802 0.71206 + SPLITR 2 1 1 + MIX 1 2 3 + R+ REFL ; +MOSTELV := GEO: MOSTELA :: + R+ VOID ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for keff +*---- +VOLMAT := SYBILT: MOSTELA :: + TITLE 'TCDR09: SYBIL TRACK MOSTELLER BENCHMARK REFLECTIVE BC ' + MAXR 4 QUA1 5 ; +LIBRARY := SHI: LIBRARY VOLMAT :: EDIT 0 NOLJ ; +PIJ := ASM: LIBRARY VOLMAT ; +FLUX := FLU: PIJ LIBRARY VOLMAT :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.7670202 ; +OUT := EDI: FLUX LIBRARY VOLMAT :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ; +PIJ VOLMAT := DELETE: PIJ VOLMAT ; +VOLMAT := SYBILT: MOSTELV :: + TITLE 'TCDR09: SYBIL TRACK MOSTELLER BENCHMARK VOID BC ' + MAXR 4 QUA1 5 ; +PIJ := ASM: LIBRARY VOLMAT ; +FLUX := FLU: FLUX PIJ LIBRARY VOLMAT :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.048279E-02 ; +OUT := EDI: OUT FLUX LIBRARY VOLMAT :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ; +OUT FLUX PIJ LIBRARY VOLMAT := DELETE: + OUT FLUX PIJ LIBRARY VOLMAT ; +ECHO "test TCDR09 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR10.c2m b/Dragon/data/tdraglib_proc/TCDR10.c2m new file mode 100644 index 0000000..d34a4cf --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR10.c2m @@ -0,0 +1,89 @@ +*---- +* TEST CASE TCDR10 +* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL +* 172-GROUP JEF2.2 DRAGLIB +* FIXED SOURCE PROBLEM IN MULTIPLICATIVE MEDIA +* +* REF: None +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + MOSTELA VOLMAT LIBRARY PIJ FLUX OUT ; +MODULE + LIB: MAC: GEO: SYBILT: SHI: ASM: FLU: EDI: DELETE: END: ; +PROCEDURE assertS assertV ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +* Fixed source of 1.0E5 in group 6 +*---- +LIBRARY := LIB: :: + EDIT 0 NMIX 3 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 600.0 O16 = O16 4.61309E-2 + U235 = U235 1.66078E-4 1 + U238 = U238 2.28994E-2 1 + MIX 2 600.0 + Zr91 = Zr91 3.83243E-2 + MIX 3 600.0 + H1H2O = H1_H2O 4.42326E-2 O16H2O = O16 2.21163E-2 + B10 = B10 2.02223E-6 B11 = B11 8.19107E-6 + ; +LIBRARY := MAC: LIBRARY :: + EDIT 0 + READ INPUT + MIX 3 FIXE + 0.0 0.0 0.0 0.0 0.0 1.0E+5 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 + ; +*---- +* Geometry MOSTELA : Annular cell with reflective BC +*---- +MOSTELA := GEO: :: TUBE 3 + RADIUS 0.0 0.39306 0.45802 0.71206 + SPLITR 2 1 1 + MIX 1 2 3 + R+ REFL ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for keff +*---- +VOLMAT := SYBILT: MOSTELA :: + TITLE 'TCDR10: SYBIL TRACK MOSTELLER BENCHMARK' + MAXR 4 QUA1 5 ; +LIBRARY := SHI: LIBRARY VOLMAT :: EDIT 0 NOLJ ; +PIJ := ASM: LIBRARY VOLMAT ; +FLUX := FLU: PIJ LIBRARY VOLMAT :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.7670202 ; +OUT := EDI: FLUX LIBRARY VOLMAT :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 ; +FLUX := DELETE: FLUX ; +FLUX := FLU: PIJ LIBRARY VOLMAT :: + TYPE S EXTE 30 UNKT 1.0E-3 ; +assertV FLUX :: 'FLUX' (*GROUP*) 10 (*REGION*) 4 2.656012E+05 ; +OUT := EDI: OUT FLUX LIBRARY VOLMAT :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 ; +OUT FLUX PIJ LIBRARY VOLMAT := DELETE: + OUT FLUX PIJ LIBRARY VOLMAT ; +ECHO "test TCDR10 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR11.c2m b/Dragon/data/tdraglib_proc/TCDR11.c2m new file mode 100644 index 0000000..86c654e --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR11.c2m @@ -0,0 +1,167 @@ +*---- +* TEST CASE TCDR11 +* CANDU-6 ANNULAR CELL +* 172-GROUP JEF2.2 DRAGLIB +* TWO GROUP BURNUP +* POWER (KW) = 615.00000 +* BURN POWER (KW/KG) = 31.97130 +* URANIUM MASS = 19.23600 +* UO2 REAL DENSITY = 10.59300 +* UO2 EFF DENSITY = 10.43750 +* UO2 TEMPERATURE = 941.28998 +* ENRICHMENT = 0.71140 +* COOLANT D2 AT % = 99.222 +* MODERATOR D2 AT % = 99.911 +* NUMBER OF DAYS = 50 +* +*---- +* Define variables +* Burnup paremeters +* a) Power +* = 31.9713 kw/kg for 0.0 to 300.0 days +* b) 69 Groups Burnup time interval Delt +* = 300 day for 0 to 300 day +* c) 2 Groups Burnup time interval Delt +* = 1 day for 0 to 1 day +* = 4 days for 1 to 5 days +* = 5 days for 5 to 10 days +* = 10 days for 10 to 50 days +* = 20 days for 50 to 150 days +* = 50 days for 150 to 300 days +* c) Days with burnup interval changes +* = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days +* d) Burnup control time variables Timei, Timef +* Timei = initial time +* Timef = final time +*---- +REAL + Power Delt Timec Timei Timef := + 31.9713 1.0 300.0 0.0 0.0 ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY ; +LINKED_LIST + CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION ; +SEQ_BINARY + INTLINS INTLINF ; +SEQ_ASCII + res ; +MODULE + GEO: EXCELT: LIB: SHI: ASM: FLU: EVO: EDI: SPH: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +PROCEDURE TCDR05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCDR05Lib :: <<iedit>> ; +*---- +* Geometry CANDU6S : 13 regions annular cluster for self-shielding +* CANDU6F : 31 regions annular cluster for transport +*---- +CANDU6S := GEO: :: TUBE 5 + R+ REFL RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 16.12171 + MIX 1 2 3 4 5 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ; + ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ; + ; +CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10 + ::: ROD1 := GEO: ROD1 SPLITR 2 1 ; + ::: ROD2 := GEO: ROD2 SPLITR 2 1 ; + ::: ROD3 := GEO: ROD3 SPLITR 2 1 ; + ::: ROD4 := GEO: ROD4 SPLITR 2 1 ; + ; +*---- +* Self-Shielding calculation EXCEL +* Transport calculation EXCEL +* Flux calculation for keff +*---- +VOLMATS INTLINS := EXCELT: CANDU6S :: + TITLE 'TCDR11: FEW GROUP BURNUP / SELF-SHIELDING TRACKING' + EDIT 0 MAXR 13 TRAK TISO 5 10.0 SYMM 12 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; +VOLMATF INTLINF := EXCELT: CANDU6F :: + TITLE 'TCDR11: FEW GROUP BURNUP / TRANSPORT TRACKING' + EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.117541 ; +EDITION := EDI: FLUX LIBRARY VOLMATF CANDU6F :: + MERG REGI 6 6 10 7 7 10 1 1 8 8 10 1 + 1 9 9 10 1 1 2 3 4 5 5 5 + 5 5 5 5 5 5 5 + COND 4.0 MICR ALL SAVE + MGEO CANDU6F + ; +EDITION := SPH: EDITION VOLMATF INTLINF ; +*---- +* 172 group Burnup +*---- +BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF :: + EDIT 3 EXPM 200.0 DEPL <<Timei>> <<Timec>> DAY POWR <<Power>> ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; +PIJ := DELETE: PIJ ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: FLUX PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9515619 ; +EDITION := EDI: EDITION FLUX LIBRARY VOLMATF CANDU6F :: + MGEO CANDU6F + ; +EDITION := SPH: EDITION VOLMATF INTLINF ; +BURNUP FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S := DELETE: + BURNUP FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S ; +*---- +* 2 group Burnup +*---- +LIBRARY := EDITION :: STEP UP 'REF-CASE0001' ; +EDITION := DELETE: EDITION ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.117564 ; +EDITION := EDI: FLUX LIBRARY VOLMATF :: + EDIT 1 MERGE COMP MICR ALL SAVE ; +EVALUATE Timec := 1.0 ; +WHILE Timei Timec < DO + EVALUATE Timef := Timei Delt + ; + IF Timei 0.0 = THEN + BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF :: + EDIT 3 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + EDIT 3 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ENDIF ; + PIJ := DELETE: PIJ ; + PIJ := ASM: LIBRARY VOLMATF INTLINF ; + FLUX := FLU: FLUX PIJ LIBRARY VOLMATF :: + TYPE K ; + EDITION := EDI: EDITION FLUX LIBRARY VOLMATF :: + SAVE ; +*---- +* change delta t for burnup and final time if required +*---- + IF Timef Timec = THEN + 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 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.075418 ; +res := EDITION ; +EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE: + EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ; +ECHO "test TCDR11 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR12.c2m b/Dragon/data/tdraglib_proc/TCDR12.c2m new file mode 100644 index 0000000..5a1162e --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR12.c2m @@ -0,0 +1,106 @@ +*---- +* TEST CASE TCDR12 +* 172-GROUP JEF2.2 DRAGLIB +* GENERATE A LIBRARY USING INFO AND OTHER OPTIONS +* +* REF: None +* +*---- +* Define variables and initialize +* Coolant properties +* a) Input +* TempCool = Coolant temperature (K) +* Purity = D2/(D2+H1) Weight % ratio in Coolant +* b) Output +* DensCool = Coolant Density (g/cm**3) +* WH1C = H1 Weight % in Coolant +* WD2C = D2 Weight % in Coolant +* WO16C = O16 Weight % in Coolant +* Fuel properties +* a) Input +* TempFuel = Fuel temperature (K) +* Enrichment = U235/(U235+U238) Weight % ratio in Fuel +* DensFuel = Fuel Density (g/cm**3) +* b) Output +* WU235F = U235 Weight % in Fuel +* WU238F = U238 Weight % in Fuel +* WO16F = O16 Weight % in Fuel +*---- +REAL + TempCool Purity TempFuel Enrichment DensFuel := + 560.66 99.95 560.66 0.72 10.437501 ; +REAL + WH1C WD2C WO16C DensCool + WU235F WU238F WO16F ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY ISOT ; +MODULE + LIB: INFO: END: ; +PROCEDURE assertS ; +*---- +* Get Coolant properties +*---- +ECHO + "Input - Coolant temperature (K) " TempCool ; +ECHO + "Input - D2/(D2+H1) Weight % ratio in Coolant" Purity ; +INFO: :: + TMP: <<TempCool>> K + PUR: <<Purity>> WGT% + CALC DENS WATER >>DensCool<< + LIB: DRAGON FIL: DLIB_J2 + ISO: 3 H1_H2O H2_D2O O16 + CALC WGT% D2O >>WH1C<< >>WD2C<< >>WO16C<< + ; +ECHO + "Output - Coolant Density (g/cm**3)" DensCool ; +ECHO + "Output - H1 Weight % in Coolant " WH1C ; +ECHO + "Output - D2 Weight % in Coolant " WD2C ; +ECHO + "Output - O16 Weight % in Coolant " WO16C ; +*---- +* Get Fuel properties +*---- +ECHO + "Input - Fuel temperature (K) " TempFuel ; +ECHO + "Input - U235/(U235+U238) Weight % ratio in Fuel" Enrichment ; +ECHO + "Input - Fuel Density (g/cm**3) " DensFuel ; +INFO: :: + ENR: <<Enrichment>> WGT% + LIB: DRAGON FIL: DLIB_J2 + ISO: 3 U235 U238 O16 + CALC WGT% UO2 >>WU235F<< >>WU238F<< >>WO16F<< + ; +ECHO + "Output - U235 Weight % in Fuel" WU235F ; +ECHO + "Output - U238 Weight % in Fuel" WU238F ; +ECHO + "Output - O16 Weight % in Fuel " WO16F ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 <<TempCool>> <<DensCool>> O16 = O16 <<WO16C>> + D2D2O = H2_D2O <<WD2C>> H1H2O = H1_H2O <<WH1C>> + MIX 2 <<TempFuel>> <<DensFuel>> O16 = O16 <<WO16F>> + U235 = U235 <<WU235F>> 1 + U238 = U238 <<WU238F>> 1 + MIX 3 COMB 1 0.5 0 0.5 + MIX 4 COMB 1 0.1 2 0.9 + ; + +ISOT := LIBRARY :: STEP UP ISOTOPESLIST STEP AT 14 ; (*U238 0004*) +assertS ISOT :: NG 35 2.219075E-01 ; +ECHO "test TCDR12 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR13.c2m b/Dragon/data/tdraglib_proc/TCDR13.c2m new file mode 100644 index 0000000..d524cd5 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR13.c2m @@ -0,0 +1,132 @@ +*---- +* TEST CASE TCDR13 +* CANDU-6 ANNULAR CELL +* 172-GROUP JEF2.2 DRAGLIB +* POWER (KW) = 615.00000 +* BURN POWER (KW/KG) = 31.97130 +* URANIUM MASS = 19.23600 +* UO2 REAL DENSITY = 10.59300 +* UO2 EFF DENSITY = 10.43750 +* UO2 TEMPERATURE = 941.28998 +* ENRICHMENT = 0.71140 +* COOLANT D2 AT % = 99.222 +* MODERATOR D2 AT % = 99.911 +* NUMBER OF DAYS = 50 +* +*---- +* Select burnup options +*---- +REAL Power BExit := 31.9713 9609.42 ; +INTEGER NbSteps := 16 ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX EDITION ISOT ; + XSM_FILE BURNUP DATABASE ; +SEQ_BINARY + INTLINS INTLINF ; +SEQ_ASCII + database.txt ; +MODULE + GEO: EXCELT: SHI: ASM: FLU: EDI: COMPO: DELETE: END: ; +PROCEDURE EvolRef ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +PROCEDURE TCDR05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCDR05Lib :: <<iedit>> ; +*---- +* Geometry CANDU6S : 13 regions annular cluster for self-shielding +* CANDU6F : 31 regions annular cluster for transport +*---- +CANDU6S := GEO: :: TUBE 5 + R+ REFL RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 16.12171 + MIX 1 2 3 4 5 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ; + ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ; + ; +CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10 + ::: ROD1 := GEO: ROD1 SPLITR 2 1 ; + ::: ROD2 := GEO: ROD2 SPLITR 2 1 ; + ::: ROD3 := GEO: ROD3 SPLITR 2 1 ; + ::: ROD4 := GEO: ROD4 SPLITR 2 1 ; + ; +*---- +* Create the reactor database +*---- +DATABASE := COMPO: :: + EDIT 5 + STEP UP 'moderator' + COMM 'Multi-parameter reactor database for moderator' ENDC + INIT + STEP UP 'fuel' + COMM 'Multi-parameter reactor database for fuel' ENDC + PARA 'BURN' IRRA + PARA 'FLUB' FLUB + INIT + ; +*---- +* Self-Shielding calculation EXCEL +* Transport calculation EXCEL +*---- +VOLMATS INTLINS := EXCELT: CANDU6S :: + TITLE 'TCDR05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + EDIT 0 MAXR 13 TRAK TISO 5 10.0 SYMM 12 ; +VOLMATF INTLINF := EXCELT: CANDU6F :: + TITLE 'TCDR05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ; +*---- +* Reflector calculation +*---- +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-INFINITY' 1 1.117541 ; +EDITION := EDI: FLUX LIBRARY VOLMATF :: + COND 4.0 MERGE MIX 0 0 0 0 1 0 0 0 0 0 MICR RES SAVE ON 'moderator' ; +DATABASE := COMPO: DATABASE EDITION :: + EDIT 3 + STEP UP * + ; +EDITION := EDI: EDITION FLUX LIBRARY VOLMATF :: + COND 4.0 MERGE COMP MICR 1 Xe135 SAVE ON 'fuel' ; +PIJ FLUX EDITION := DELETE: PIJ FLUX EDITION ; +*---- +* Burnup calculation using procedure EvolRef +*---- +DATABASE BURNUP FLUX LIBRARY := EvolRef LIBRARY VOLMATS + VOLMATF INTLINS INTLINF :: + <<BExit>> <<Power>> <<NbSteps>> ; +*---- +* Export and access the database +*---- +database.txt := DATABASE ; + +ISOT := DATABASE :: STEP UP fuel + STEP UP MIXTURES STEP AT 1 + STEP UP CALCULATIONS STEP AT 5 + STEP UP ISOTOPESLIST STEP AT 1 ; +assertS ISOT :: 'NWT0' 1 1.149079E1 ; +assertS ISOT :: 'NWT0' 2 2.586806E1 ; +ISOT := DELETE: ISOT ; + +ISOT := DATABASE :: STEP UP moderator + STEP UP MIXTURES STEP AT 1 + STEP UP CALCULATIONS STEP AT 1 + STEP UP ISOTOPESLIST STEP AT 1 ; +assertS ISOT :: 'NWT0' 1 9.095950E-2 ; +assertS ISOT :: 'NWT0' 2 3.246065E-1 ; + +INTLINF INTLINS := DELETE: INTLINF INTLINS ; +ECHO "test TCDR13 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR14.c2m b/Dragon/data/tdraglib_proc/TCDR14.c2m new file mode 100644 index 0000000..768e2f7 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR14.c2m @@ -0,0 +1,293 @@ +*---- +* TEST CASE TCDR14 +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* DEPLETION WITH CPO AND MULTICOMPO DATABASE +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +PARAMETER saCOMPO :: ::: SEQ_ASCII saCOMPO ; ; +LINKED_LIST + ASSMB DISCR LIBRARY CP CALC OUT COMPO MCOMPO ISOT ASSMB_HOM BURNUP + RECORD MTRACK ; +MODULE + LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: CPO: COMPO: DELETE: END: EVO: + BIVACT: SPH: ; +PROCEDURE assertS ; +INTEGER ISTEP ; +STRING HSTEP ; +REAL + Power Delt Timec Timei Timef := + 31.9713 1.0 300.0 0.0 0.0 ; +*---- +* MULTICOMPO +*---- +MCOMPO := COMPO: :: + EDIT 3 + MAXCAL 3 + COMM 'generated for TINST test case by TCDR14' ENDC + PARA 'BURN' IRRA + INIT + ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + DEPL LIB: DRAGON FIL: DLIB_J2 + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 NOEV + H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2 + B10 = B10 4.71444E-6 B11 = B11 1.90959E-5 + MIX 2 579.9 NOEV O16 = O16 3.06711E-4 + Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4 + Zr91 = Zr0 4.18621E-2 + MIX 3 579.9 NOEV + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + MoNat = Mo95 4.03755E-6 + Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr0 8.92427E-4 + MIX 4 579.9 NOEV O16 = O16 2.87335E-4 + Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4 + Zr91 = Zr0 3.92175E-2 + MIX 5 579.9 NOEV + H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2 + Cr52 = Cr52 2.24991E-5 + Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5 + MoNat = Mo95 1.89281E-6 + Al27 = Al27 1.10277E-6 + Mn55 = Mn55 1.94976E-7 + B10 = B10 4.66484E-6 B11 = B11 1.88950E-5 + Zr91 = Zr0 4.18372E-4 + MIX 6 579.9 NOEV + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + Cr52 = Cr52 2.11122E-5 + Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5 + MoNat = Mo95 1.77614E-6 + Al27 = Al27 1.03479E-6 + Mn55 = Mn55 1.82957E-7 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + Zr91 = Zr0 3.92583E-4 + MIX 7 579.9 NOEV + H1H2O = H1_H2O 4.72020E-2 O16H2O = O16 2.36010E-2 + Cr52 = Cr52 1.96591E-5 + Fe56 = Fe56 1.82630E-5 Ni58 = Ni58 4.65011E-5 + MoNat = Mo95 1.65389E-6 + Al27 = Al27 9.63569E-7 + Mn55 = Mn55 1.70365E-7 + B10 = B10 4.67110E-6 B10 = B10 1.89203E-5 + Zr91 = Zr0 3.65562E-4 + MIX 8 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + ; +*---- +* Geometry ASSMB : a 17 X 17 normal PWR assembly +* contains C1 : cell without fuel +* C2 : normal fuel cell +* C3 : peripheral cell +* C4 : corner cell +*---- +ASSMB := GEO: :: CAR2D 9 9 + X- DIAG X+ REFL Y- SYME Y+ DIAG + CELL C1 C2 C2 C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C2 C2 C2 C3 + C2 C2 C2 C2 C2 C2 C3 + C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C3 + C1 C2 C2 C3 + C2 C2 C3 + C2 C3 + C4 + MERGE 1 3 12 11 12 12 11 12 15 + 4 6 5 6 6 5 6 8 + 13 5 6 6 5 6 8 + 2 5 5 10 5 8 + 13 5 5 6 8 + 2 5 7 8 + 13 7 8 + 14 8 + 9 + TURN H H B H H B H H A + H G G H G G H A + A E E F E E A + H H F H H A + H E G H A + H H A A + H A A + A A + A + ::: C1 := GEO: CARCEL 2 + MESHX 0.0 1.26472 MESHY 0.0 1.26472 + RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ; + ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ; + ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ; + ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ; + ; +ASSMB_HOM := GEO: :: CAR2D 9 9 + X- DIAG X+ REFL Y- SYME Y+ DIAG + MIX 1 2 2 1 2 2 1 2 2 + 2 2 2 2 2 2 2 2 + 2 2 2 2 2 2 2 + 1 2 2 1 2 2 + 2 2 2 2 2 + 1 2 2 2 + 2 2 2 + 2 2 + 2 + MESHX 0.0 1.26472 2.52944 3.79416 5.05888 6.3236 7.58832 8.85304 + 10.11776 11.43248 + MESHY 0.0 1.26472 2.52944 3.79416 5.05888 6.3236 7.58832 8.85304 + 10.11776 11.43248 + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* OUT calculation for B1 homogeneous leakage +* Editing using SPH model for transport-diffusion +*---- +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR14: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 ; +LIBRARY := SHI: LIBRARY DISCR :: LJ ; +CP := ASM: LIBRARY DISCR ; +CALC := FLU: CP LIBRARY DISCR :: + TYPE B B1 PNL ; +assertS CALC :: 'K-INFINITY' 1 1.256372 ; + +OUT := EDI: CALC LIBRARY DISCR ASSMB_HOM :: + EDIT 3 UPS SAVE MICR ALL MERGE MIX 1 1 1 2 2 2 2 2 + COND 4.0 MGEO ASSMB_HOM + ; +MTRACK := BIVACT: ASSMB_HOM :: DUAL 1 1 EDIT 0 ; +OUT := SPH: OUT MTRACK ; +MTRACK := DELETE: MTRACK ; +BURNUP LIBRARY := EVO: LIBRARY CALC DISCR :: + EDIT 0 + SAVE <<Timei>> DAY POWR <<Power>> +; +MCOMPO := COMPO: MCOMPO OUT BURNUP :: + EDIT 3 + MACRO + SET <<Timei>> DAY + ; +*---- +* 172 group Burnup +*---- +EVALUATE Timec := 1.0 ; +WHILE Timei Timec < DO + EVALUATE Timef := Timei Delt + ; + BURNUP LIBRARY := EVO: BURNUP LIBRARY CALC DISCR :: + EDIT 3 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + CP := DELETE: CP ; + CP := ASM: LIBRARY DISCR ; + CALC := FLU: CALC CP LIBRARY DISCR :: + TYPE B B1 PNL ; + OUT := EDI: OUT CALC LIBRARY DISCR ASSMB_HOM :: + EDIT 3 UPS SAVE MICR ALL MERGE MIX 1 1 1 2 2 2 2 2 + COND 4.0 MGEO ASSMB_HOM + ; + MTRACK := BIVACT: ASSMB_HOM :: DUAL 1 1 EDIT 0 ; + OUT := SPH: OUT MTRACK ; + MTRACK := DELETE: MTRACK ; + BURNUP LIBRARY := EVO: BURNUP LIBRARY CALC DISCR :: + EDIT 0 + SAVE <<Timef>> DAY POWR <<Power>> + ; + MCOMPO := COMPO: MCOMPO OUT BURNUP :: + EDIT 3 + MACRO + SET <<Timef>> DAY + ; +*---- +* change delta t for burnup and final time if required +*---- + IF Timef Timec = THEN + 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 ; +assertS CALC :: 'K-INFINITY' 1 1.211175 ; + +COMPO := CPO: OUT BURNUP :: + BURNUP REF-CASE NAME COMPO ; + +RECORD := COMPO :: STEP UP 'COMPO 2' ; +assertS RECORD :: 'BURNUP' 4 3.197132E+02 ; +*---- +* 2 group Burnup initialization +*---- +EVALUATE Delt := 1.0 ; +EVALUATE Timei := 0.0 ; +CALC CP DISCR LIBRARY BURNUP := DELETE: CALC CP DISCR LIBRARY BURNUP ; +EVALUATE HSTEP := "REF-CASE0001" ; +ECHO "-----> Perform micro-depletion at hstep=" HSTEP ; +LIBRARY := OUT :: STEP UP <<HSTEP>> ; +DISCR := BIVACT: ASSMB_HOM :: + DUAL 1 1 EDIT 0 ; +CP := ASM: LIBRARY DISCR :: ARM ; +CALC := FLU: CP LIBRARY DISCR :: + TYPE K LKRD IDEM ; +assertS CALC :: 'K-INFINITY' 1 1.256390 ; +*---- +* 2 group Burnup +*---- +EVALUATE Timec := 1.0 ; +EVALUATE ISTEP := 1 ; +WHILE Timei Timec < DO + ECHO "-----> Perform micro-depletion at hstep=" HSTEP ; + EVALUATE Timef := Timei Delt + ; + IF Timei 0.0 = THEN + BURNUP LIBRARY := EVO: LIBRARY CALC DISCR :: + EDIT 3 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY CALC DISCR :: + EDIT 3 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ENDIF ; + EVALUATE ISTEP := ISTEP 1 + ; + IF ISTEP 10 < THEN + EVALUATE HSTEP := "REF-CASE000" ISTEP I_TO_S + ; + ELSEIF ISTEP 100 < THEN + EVALUATE HSTEP := "REF-CASE00" ISTEP I_TO_S + ; + ELSEIF ISTEP 10000 < THEN + EVALUATE HSTEP := "REF-CASE0" ISTEP I_TO_S + ; + ENDIF ; + LIBRARY := DELETE: LIBRARY ; + LIBRARY := OUT :: STEP UP <<HSTEP>> ; + BURNUP LIBRARY := EVO: BURNUP LIBRARY DISCR :: + EDIT 2 NOSA SET <<Timef>> DAY ; + CP := DELETE: CP ; + CP := ASM: LIBRARY DISCR :: ARM ; + CALC := FLU: CALC CP LIBRARY DISCR :: + TYPE K LKRD IDEM ; +*---- +* change delta t for burnup and final time if required +*---- + IF Timef Timec = THEN + 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 ; +saCOMPO := MCOMPO ; +assertS CALC :: 'K-INFINITY' 1 1.212845 ; + +ECHO "test TCDR14 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR15.c2m b/Dragon/data/tdraglib_proc/TCDR15.c2m new file mode 100644 index 0000000..7b70a14 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR15.c2m @@ -0,0 +1,200 @@ +*-------------------------------------------------------------------- +* TEST CASE TCDR15 +* HTTR Fuel +* Block ID:f393320 +* JAVIER ORTENSI +* Jef2.2 evaluation +* XMAS172 Group Structure +* Description: models a fuel block in infinite medium with spatial +* homogenization of the burnable poison region (BP x-section) +* and the rest of the assembly (background x-section). +*-------------------------------------------------------------------- +* Define STRUCTURES and MODULES used +*-------------------------------------------------------------------- +LINKED_LIST LATTIC LIBRARY LIBRARY2 OUT1G OUT2G FLUX TRACK SYS ; +SEQ_ASCII EDIT.out ; +SEQ_BINARY TRACK_FIL ; +PROCEDURE assertS ; + +MODULE GEO: EXCELT: LIB: ASM: FLU: EDI: USS: END: DELETE: ; +*-------------------------------------------------------------------- +* Define variables used +*-------------------------------------------------------------------- +LIBRARY := LIB: :: + NMIX 14 CTRA WIMS + SUBG + ANIS 2 + MIXS LIB: DRAGON FIL: DLIB_J2 +* Fuel Kernel + MIX 1 300.0 + B10 = 'B10' 1.1010E-07 + B11 = 'B11' 4.4317E-07 + O16 = 'O16' 4.8150E-02 + O17 = 'O17' 1.8304E-05 + U235 = 'U235' 8.0367E-04 1 IRSET 0.0 81 + U238 = 'U238' 2.3280E-02 1 IRSET 0.0 81 +* Porous C buffer + MIX 2 300.0 + C12 = 'C0_GR' 5.8060E-02 +* IPyC + MIX 3 300.0 + C12 = 'C0_GR' 9.3909E-02 +* SiC coating layer + MIX 4 300.0 + C12 = 'C0' 4.8151E-02 + SI28 = 'Si0' 4.4410E-02 + SI29 = 'Si0' 2.2487E-03 + SI30 = 'Si0' 1.4927E-03 +* OPyC + MIX 5 300.0 + C12 = 'C0_GR' 9.3007E-02 +* Graphite matrix (Compact) + MIX 6 300.0 + C12 = 'C0_GR' 8.5085E-02 + B10 = 'B10' 2.4543E-08 + B11 = 'B11' 9.8790E-08 +* Graphite sleeve + MIX 7 300.0 + C12 = 'C0_GR' 8.8745E-02 + B10 = 'B10' 7.1811E-09 + B11 = 'B11' 2.8905E-08 +* Helium Coolant + MIX 8 300.0 + HE3 = 'He3' 1.3230E-09 + HE4 = 'He4' 9.6573E-04 +* Block Graphite inside fuel cell + MIX 9 300.0 + C12 = 'C0_GR' 8.8745E-02 + B10 = 'B10' 7.8482E-09 + B11 = 'B11' 3.1590E-08 +* Block Graphite around fuel handling hole (rho adjusted) + MIX 10 300.0 + C12 = 'C0_GR' 6.7791E-02 + B10 = 'B10' 5.9951E-09 + B11 = 'B11' 2.4131E-08 + HE3 = 'He3' 3.1239E-10 + HE4 = 'He4' 2.2802E-04 +* Block Graphite outside fuel cell (rho adjusted) + MIX 11 300.0 + C12 = 'C0_GR' 9.0845E-02 + B10 = 'B10' 8.0339E-09 + B11 = 'B11' 3.2338E-08 +* Burnable Poison + MIX 12 300.0 + C12 = 'C0_GR' 8.7821E-02 + B10 = 'B10' 3.5588E-04 + B11 = 'B11' 1.4071E-03 +* Helium Coolant + MIX 13 300.0 + HE3 = 'He3' 1.3230E-09 + HE4 = 'He4' 9.6573E-04 +* Block Graphite around BP + MIX 14 300.0 + C12 = 'C0_GR' 8.8745E-02 + B10 = 'B10' 7.8482E-09 + B11 = 'B11' 3.1590E-08 +; +*-------------------------------------------------------------------- +* FUEL BLOCK SECTION +*-------------------------------------------------------------------- +* Define mixture for the fuel block +*-------------------------------------------------------------------- +LATTIC := GEO: :: HEX 61 + HBC COMPLETE REFL + CELL +* Ring 0 (center cell) +FH +* Ring 1 (6 cells) +F F F F F F +* Ring 2 (12 cells) +F F F F F F F F F F F F +* Ring 3 (18 cells) +BPE F F F F F BP F F F F F BP F F F F F +* Ring 4 (24 cells) +C C C C C C C C C C C C C C C C C C C C C C C C +* Spherical microstruture model + BIHET SPHE 1 1 + (* NS *) 5 + (* RS *) 0.0 0.0297 0.0359 0.0389 0.0418 0.0464 + (* milie *) 15 + (* mixdil *) 6 + (* loop over NMISTR *) + (* fract *) 0.3040 + (* mixgr *) 1 2 3 4 5 +* Side length based on the pin pitch +* Fuel cell + ::: F := GEO: HEXCEL 5 + SIDE 2.973354 + RADIUS 0.0 0.4995 1.3000 1.3130 1.7000 1.9800 + MIX 8 15 8 7 8 9 ; +* Fuel handling hole + ::: FH := GEO: HEXCEL 0 + SIDE 2.973354 + MIX 10 ; +* Burnable Poison Cell (Empty) +* To be homogenized with fuel (use mat 8&9) + ::: BPE := GEO: HEXCEL 1 + SIDE 2.973354 + RADIUS 0.0 0.6950 + MIX 8 9 ; +* Burnable Poison Cell + ::: BP := GEO: HEXCEL 2 + SIDE 2.973354 + RADIUS 0.0 0.6950 0.7500 + MIX 12 13 14 ; +* graphite cell + ::: C := GEO: HEXCEL 0 + SIDE 2.973354 + MIX 11 ; +; +*-------------------------------------------------------------------- +* Use SYBILT, USS, ASM and FLU for flux solution +*-------------------------------------------------------------------- +*---- +* Tracking calculation +*---- +TRACK TRACK_FIL := EXCELT: LATTIC :: + EDIT 0 + MAXR 400 + ANIS 2 + TRAK TISO 8 20.0 +; +*--- +* Self Shielding calculation +*--- +LIBRARY2 := USS: LIBRARY TRACK TRACK_FIL :: + EDIT 1 +; +*--- +* Multigroup Flux Calculation +*--- + SYS := ASM: LIBRARY2 TRACK TRACK_FIL :: + EDIT 1 ; + FLUX := FLU: SYS LIBRARY2 TRACK TRACK_FIL :: + EDIT 1 + TYPE B B1 PNL + EXTE 100 1E-5 THER 5 1E-5 ; + OUT1G := EDI: FLUX LIBRARY2 TRACK LATTIC :: + EDIT 3 + MERG MIX 1 1 1 1 1 1 1 1 1 1 1 0 0 0 + COND 7.4080E06 3.6790E06 6.7210E05 1.1110E05 1.9310E04 3.3550E03 + 1.5850E03 7.4850E02 2.7540E02 1.3010E02 6.1440E01 2.9000E01 + 1.3700E01 8.3200E00 5.0400E00 2.3800E00 1.2900E00 6.5000E-01 + 3.5000E-01 2.0000E-01 1.2000E-01 8.0000E-02 5.0000E-02 2.0000E-02 + 1.0000E-02 + MICR ISOTXS ASCII ALL + SAVE ; + OUT1G := EDI: OUT1G FLUX LIBRARY2 TRACK :: + EDIT 3 + MERG MIX 0 0 0 0 0 0 0 0 0 0 0 1 1 1 + COND 7.4080E06 3.6790E06 6.7210E05 1.1110E05 1.9310E04 3.3550E03 + 1.5850E03 7.4850E02 2.7540E02 1.3010E02 6.1440E01 2.9000E01 + 1.3700E01 8.3200E00 5.0400E00 2.3800E00 1.2900E00 6.5000E-01 + 3.5000E-01 2.0000E-01 1.2000E-01 8.0000E-02 5.0000E-02 2.0000E-02 + 1.0000E-02 + MICR ISOTXS ASCII ALL + SAVE ; +assertS FLUX :: K-INFINITY 1 1.145036 ; +ECHO "test TCDR15 completed" ; +EDIT.out := OUT1G ; +END: ; diff --git a/Dragon/data/tdraglib_proc/TCDR16.c2m b/Dragon/data/tdraglib_proc/TCDR16.c2m new file mode 100644 index 0000000..e228900 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR16.c2m @@ -0,0 +1,152 @@ +*-------------------------------------------------------------------- +* TEST CASE TCDR16 +* DB-VHTR Fuel Block Lattice calculation +* Fresh fuel TRU loading +* Author : V. Descotes +* Core corner block(6,1), no BP, no CR. +*-------------------------------------------------------------------- +* +*-------------------------------------------------------------------- +* Define STRUCTURES and MODULES used +*-------------------------------------------------------------------- +LINKED_LIST LIBRARY LIBRARY2 GEOM TRACK SYS FLUX ; +SEQ_BINARY TRACK_FIL ; +MODULE GEO: LIB: ASM: FLU: USS: EXCELT: MCCGT: END: ; +PROCEDURE assertS ; + +*-------------------------------------------------------- +* GEOMETRY +*-------------------------------------------------------- +* ---- +* F = Fuel cell +* C = coolant hole cell (large) +* K = coolant hole cell (small) +* G = Full density graphite cell +* J = Reduced density graphite cell +* ---- +GEOM := GEO: :: HEX 1 + EDIT 1 + HBC COMPLETE REFL + CELL F0 + ::: F0 := GEO: HEXCEL 2 + SIDE 1.08542 + RADIUS 0.0 0.6225 0.6350 + SECT -1 + MIX 9 9 9 9 9 9 + 7 7 7 7 7 7 + 8 8 8 8 8 8 ; +* +*--------------------------------------------------------- +* SPHERICAL MICRO-STRUCTURES * +* NMISTR = NUMBER OF MICRO STRUCTURE TYPES IN REGION +* NMILG = NUMBER OF MICRO STRUCTURE REGIONS +* NS = ARRAY OF SUB REGIONS IN THE MICRO STRUCTURES; LEN=NMILG +* RS = RADIUS OF MICROS STRUCTURES; LEN=NS(I); I=1,NMISTR +* milie = COMPOSITION OF EACH shell OF MICRO STRUCTURE; LEN=NMISTR; +* !!!!!!!! milie #'s ARE > THAN MIX #'S !!!!!!!! +* mixdil = BASE COMPOSITION OF EACH REGION; LEN=NMILG +* fract = MICROSTRUCTURE TYPE VOLUME FRACTION IN REGION LEN=NMILG +* mixgr = LIBRARY MIXTURES FOR EACH MICROSTRUCTURE SHELL; LEN=NS(I) + +* NMISTR, NMILG + BIHET SPHE 1 1 + (* NS *) 5 + (* RS *) 0.0 0.0100 0.0220 0.0260 0.0295 0.0335 + (* milie *) 9 + (* mixdil *) 1 + (* loop over NMISTR *) + (* fract *) 0.175016 + (* mixgr *) 2 3 4 5 6 +; + +*-------------------------------------------------------------------- +* Use EXCELT-MCCGT, USS, ASM and FLU for flux solution +*-------------------------------------------------------------------- + +*** +* Tracking EXCELT-MCCGT +*** +TRACK TRACK_FIL := EXCELT: GEOM :: + EDIT 1 + MAXR 50000 + ANIS 1 + ALLG + QUAB 5 HEBE + PISO + TRAK TISO 4 10.0 +; + +TRACK := MCCGT: TRACK TRACK_FIL GEOM :: + EDIT 0 LCMD 3 + EPSI 5E-5 + MAXI 100 KRYL 10 SCR 1 AAC 1 TMT +; + +*-------------------------------------------------------------------- +* CREATION OF THE LIBRARY +*-------------------------------------------------------------------- +LIBRARY := LIB: :: + EDIT 0 + NMIX 8 CTRA APOL ANIS 2 + SUBG + MIXS LIB: DRAGON FIL: DLIB_J2 + +*---- +* Mixtures for fuel cell F0 +*---- +* Graphite matrix + MIX 1 293.6 + GRAB = 'C0_GR' 8.524E-02 +* Fuel Kernel + MIX 2 293.6 + O16 = 'O16' 4.436E-02 + NP27 = 'Np237' 1.522E-03 1 + PU38 = 'Pu238' 6.464E-04 1 + PU39 = 'Pu239' 1.099E-02 1 + PU40 = 'Pu240' 5.084E-03 1 + PU41 = 'Pu241' 1.937E-03 1 + PU42 = 'Pu242' 1.074E-03 1 + AM41 = 'Am241' 6.163E-04 1 + AM42 = 'Am242m' 4.385E-06 1 + AM43 = 'Am243' 3.057E-04 1 +* Porous C buffer + MIX 3 293.6 + GRAB = 'C0_GR' 5.265E-02 +* IPyC + MIX 4 293.6 + GRAB = 'C0_GR' 9.526E-02 +* SiC coating layer + MIX 5 293.6 + SI0 = 'Si0' 4.402E-02 + C12 = 'C0' 4.772E-02 +* OPyC + MIX 6 293.6 + GRAB = 'C0_GR' 9.526E-02 +* Helium Coolant surrounding fuel + MIX 7 293.6 + HE4 = 'He4' 7.0000E-04 +* Block Graphite (full density surrounding fuel) + MIX 8 293.6 + GRAB = 'C0_GR' 8.774E-02 +; + +*--- +* Self Shielding calculation +*--- +LIBRARY2 := USS: LIBRARY TRACK TRACK_FIL :: + EDIT 0 PASS 2 ARM ; + +*--- +* Multigroup Flux Calculation +*--- + SYS := ASM: LIBRARY2 TRACK TRACK_FIL :: + EDIT 0 ARM PNOR HELI ; + + FLUX := FLU: SYS LIBRARY2 TRACK TRACK_FIL :: + EDIT 1 + TYPE K P1 EXTE 300 1E-5 THER 10 1E-5 ; + +assertS FLUX :: 'K-EFFECTIVE' 1 1.144677 ; +ECHO "test TCDR16 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR17.c2m b/Dragon/data/tdraglib_proc/TCDR17.c2m new file mode 100644 index 0000000..ce2751c --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR17.c2m @@ -0,0 +1,206 @@ +*---- +* TEST CASE TCDR17 +* Test editing capability and create a database for a fuel ring +* depletion +* 172-GROUP JEF2.2 DRAGLIB +*---- +* Define STRUCTURES and MODULES used +*---- +PARAMETER saCOMPO2 :: ::: SEQ_ASCII saCOMPO2 ; ; +LINKED_LIST + GEOM DISCR LIBRARY CP CALC OUT BURNUP DATABASE ; +MODULE + LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: DELETE: + EVO: COMPO: END: ; +REAL + Power Delt Timec Timei Timef := + 31.9713 1.0 300.0 0.0 0.0 ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +* - Mix 1 : fuel 1 +* - Mix 2 : fuel 2 +* - Mix 3 : water +*---- +LIBRARY := LIB: :: + EDIT 0 + NMIX 3 CTRA WIMS + ANIS 2 + DEPL LIB: DRAGON FIL: DLIB_J2 + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 600.0 + O16 = O16 4.61309E-2 1 + U235 = U235 1.66078E-4 1 + U238 = U238 2.28994E-2 1 + Np239 = Np239 0.0 1 + Pu239 = Pu239 0.0 1 + Pu240 = Pu240 0.0 1 + Pu241 = Pu241 0.0 1 + Rh105 = Rh105 0.0 1 + Cd113 = Cd113 0.0 1 + I135 = I135 0.0 1 + Xe135 = Xe135 0.0 1 + Nd145 = Nd145 0.0 1 + Nd146 = Nd146 0.0 1 + Nd147 = Nd147 0.0 1 + Nd148 = Nd148 0.0 1 + Pm147 = Pm147 0.0 1 + Pm148 = Pm148 0.0 1 + Pm148m = Pm148m 0.0 1 + Pm149 = Pm149 0.0 1 + Pm151 = Pm151 0.0 1 + Sm147 = Sm147 0.0 1 + Sm149 = Sm149 0.0 1 + Sm150 = Sm150 0.0 1 + Sm151 = Sm151 0.0 1 + Sm152 = Sm152 0.0 1 + Eu153 = Eu153 0.0 1 + Eu154 = Eu154 0.0 1 + Eu155 = Eu155 0.0 1 + Eu156 = Eu156 0.0 1 + Gd157 = Gd157 0.0 1 + MIX 2 600.0 + O16 = O16 4.61309E-2 1 + U235 = U235 3.2E-4 1 + U238 = U238 2.28994E-2 1 + Np239 = Np239 0.0 1 + Pu239 = Pu239 0.0 1 + Pu240 = Pu240 0.0 1 + Pu241 = Pu241 0.0 1 + Rh105 = Rh105 0.0 1 + Cd113 = Cd113 0.0 1 + I135 = I135 0.0 1 + Xe135 = Xe135 0.0 1 + Nd145 = Nd145 0.0 1 + Nd146 = Nd146 0.0 1 + Nd147 = Nd147 0.0 1 + Nd148 = Nd148 0.0 1 + Pm147 = Pm147 0.0 1 + Pm148 = Pm148 0.0 1 + Pm148m = Pm148m 0.0 1 + Pm149 = Pm149 0.0 1 + Pm151 = Pm151 0.0 1 + Sm147 = Sm147 0.0 1 + Sm149 = Sm149 0.0 1 + Sm150 = Sm150 0.0 1 + Sm151 = Sm151 0.0 1 + Sm152 = Sm152 0.0 1 + Eu153 = Eu153 0.0 1 + Eu154 = Eu154 0.0 1 + Eu155 = Eu155 0.0 1 + Eu156 = Eu156 0.0 1 + Gd157 = Gd157 0.0 1 + I129 = I129 0.0 1 + MIX 3 600.0 + H1H2O = H1_H2O 4.42326E-2 O16H2O = O16 2.21163E-2 + B10 = B10 2.02222E-6 B11 = B11 8.19107E-6 ; +*---- +* 2D Geometry +*---- +GEOM := GEO: :: CARCEL 2 + X- REFL X+ REFL MESHX 0.0 1.26209 + Y- REFL Y+ REFL MESHY 0.0 1.26209 + RADIUS 0.0 0.39306 0.45802 SPLITR 2 1 + MIX 1 2 3 ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +*---- +DISCR := SYBILT: GEOM :: + MAXR 4 QUA1 5 QUA2 6 5 ; +LIBRARY := SHI: LIBRARY DISCR :: EDIT 0 NOLJ ; + +CP := ASM: LIBRARY DISCR ; +CALC := FLU: CP LIBRARY DISCR :: EDIT 0 + TYPE K B0 PNL ; + +*---- +* During the homogeneization, we keep a single mixture with +* all the isotopes and the homogeneized cross sections (ALLX) +*---- +OUT := EDI: CALC LIBRARY DISCR :: + EDIT 1 UPS SAVE ON 'fuel' MERGE COMP MICR ALLX 29 U235 U238 Np239 + Pu239 Pu240 Pu241 Rh105 Cd113 I135 Xe135 Nd145 Nd146 Nd147 Nd148 + Pm147 Pm148 Pm148m Pm149 Pm151 Sm147 Sm149 Sm150 Sm151 Sm152 + Eu153 Eu154 Eu155 Eu156 Gd157 + COND 4.0 ; +*--- +* Creation of databases for other DONJON and DRAGON calculations +*--- +DATABASE := COMPO: :: + EDIT 5 + STEP UP 'fuel' + COMM 'Multi-parameter reactor database for fuel' ENDC + PARA 'BURN' IRRA + PARA 'FLUB' FLUB + INIT ; +*--- +* Register the isotopic properties for the initial burnup +*--- +BURNUP LIBRARY := EVO: LIBRARY CALC DISCR :: + SAVE 0.0 DAY POWR <<Power>> ; +assertS CALC :: 'K-EFFECTIVE' 1 9.225872E-01 ; +DATABASE := COMPO: DATABASE OUT BURNUP LIBRARY :: + EDIT 0 + ALLX + STEP UP * + SET 0.0 DAY ; +*---- +* Perform depletion +*---- +EVALUATE Timec := 1.0 ; +EVALUATE Delt := 1.0 ; +EVALUATE Timei := 0.0 ; + +WHILE Timei Timec < DO + EVALUATE Timef := Timei Delt + ; + IF Timei 0.0 = THEN + BURNUP LIBRARY := EVO: BURNUP LIBRARY CALC DISCR :: + EDIT 0 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY CALC DISCR :: + EDIT 0 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ENDIF ; + + CP := DELETE: CP ; + CP := ASM: LIBRARY DISCR :: ; + CALC := FLU: CALC CP LIBRARY DISCR :: TYPE K B0 PNL ; + + BURNUP LIBRARY := EVO: BURNUP LIBRARY CALC DISCR :: + SAVE <<Timef>> DAY POWR <<Power>> ; + + OUT := EDI: OUT CALC LIBRARY DISCR :: + EDIT 1 UPS SAVE ON 'fuel' MERGE COMP MICR ALLX 29 U235 U238 Np239 + Pu239 Pu240 Pu241 Rh105 Cd113 I135 Xe135 Nd145 Nd146 Nd147 Nd148 + Pm147 Pm148 Pm148m Pm149 Pm151 Sm147 Sm149 Sm150 Sm151 Sm152 + Eu153 Eu154 Eu155 Eu156 Gd157 + COND 4.0 ; + + DATABASE := COMPO: DATABASE OUT BURNUP LIBRARY :: + EDIT 0 + ALLX + STEP UP * + SET <<Timef>> DAY ; +*---- +* change delta t for burnup and final time if required +*---- + IF Timef Timec = THEN + IF Timec 10.0 = THEN + EVALUATE Delt Timec := 10.0 200.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 ; + +saCOMPO2 := DATABASE ; + +assertS CALC :: 'K-EFFECTIVE' 1 9.392296E-01 ; +ECHO "test TCDR17 completed" ; +END: ; +QUIT . diff --git a/Dragon/data/tdraglib_proc/TCDR18.c2m b/Dragon/data/tdraglib_proc/TCDR18.c2m new file mode 100644 index 0000000..e3dba38 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR18.c2m @@ -0,0 +1,83 @@ +*---- +* TEST CASE TCDR18 +* SOURCE PROBLEM WITH FISSION +* 172-GROUP JEF2.2 DRAGLIB +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY LIBRARY2 GEOMETRY TRACK SYS FLUX EDITION GROUP ; +SEQ_BINARY + YTH ; +*---- +MODULE + GEO: EXCELT: MAC: LIB: SHI: ASM: FLU: EDI: GREP: END: ; +PROCEDURE assertS ; +*--- +LIBRARY := LIB: :: + NMIX 2 CTRA WIMS + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 COMB 0 1.0 + MIX 2 300.0 + U235 = U235 8.8365E-05 1 + U238 = U238 1.2133E-02 1 + O16 = O16 4.1252E-02 + H1H2O = H1_H2O 3.3615E-02 + ; +GEOMETRY := GEO: :: CAR2D 10 5 + X- REFL X+ REFL Y- REFL Y+ REFL + MESHX 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 + MESHY 0.0 0.2 0.4 0.6 0.8 1.0 + MIX 1 2 2 2 2 2 2 2 2 2 + 1 2 2 2 2 2 2 2 2 2 + 1 2 2 2 2 2 2 2 2 2 + 1 2 2 2 2 2 2 2 2 2 + 1 2 2 2 2 2 2 2 2 2 ; +*--- +TRACK YTH := EXCELT: GEOMETRY :: + MAXR 50 + TRAK TISO 12 20.0 + ; +LIBRARY2 := SHI: LIBRARY TRACK YTH :: + NOLJ + ; +LIBRARY2 := MAC: LIBRARY2 :: + READ INPUT + MIX 1 FIXE + 0.0 0.0 0.0 0.0 10.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 + ; +*--- +SYS := ASM: LIBRARY2 TRACK YTH :: + PIJ ; +FLUX := FLU: SYS LIBRARY2 TRACK YTH :: + TYPE S ; +*-- +EDITION := EDI: LIBRARY2 FLUX TRACK :: + EDIT 2 + MERG COMP + COND NONE + SAVE ON 'EDITMOD' ; +*-- +GROUP := EDITION :: STEP UP EDITMOD STEP UP MACROLIB STEP UP GROUP + STEP AT 1 ; +assertS GROUP :: 'NTOT0' 1 0.1447909 ; +ECHO "test TCDR18 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR19.c2m b/Dragon/data/tdraglib_proc/TCDR19.c2m new file mode 100644 index 0000000..26266db --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR19.c2m @@ -0,0 +1,158 @@ +*-------------------------------------------------------------------- +* TEST CASE TCDR19 +* VHTRC Fuel Cell B-2 +* JAVIER ORTENSI +* CP METHOD +* Material Temp 300 K +* ENDFB-VII +* XMAS172 Group Structure +* Based on design data from +* NEA/NSC VHTRC-GCR-EXP-001 +*-------------------------------------------------------------------- +* Define STRUCTURES and MODULES used +*-------------------------------------------------------------------- +LINKED_LIST + GEOMETRY LIBRARY LIBRARY2 SYS TRACK FLUX ; +SEQ_BINARY TRACK_FIL ; +MODULE + GEO: EXCELT: MCCGT: LIB: ASM: FLU: EDI: END: USS: ; +PROCEDURE assertS ; +*---- +* Mix number assignments +* 1 = fuel kernel material (UO2) +* 2 = 1st coating +* 3 = 2nd coating +* 4 = compact matrix +* 5 = graphite sheath +* 6 = air +* 7 = block graphite +* 8 = fuel compact mixture +*---- +LIBRARY := LIB: :: + NMIX 7 + CTRA NONE + SUBG + MIXS LIB: DRAGON FIL: DLIB_J2 +* Fuel Kernel + MIX 1 300. + U234 = 'U234' 3.19004E-06 1 + U235 = 'U235' 4.70582E-04 1 + U236 = 'U236' 3.74868E-07 1 + U238 = 'U238' 2.27638E-02 1 + O16 = 'O16' 4.65455E-02 + B10 = 'B10' 1.84810E-08 + B11 = 'B11' 7.43882E-08 +* 1st coating + MIX 2 300. + C12 = 'C0_GR' 5.96649E-02 + B10 = 'B10' 2.11059E-09 + B11 = 'B11' 8.49539E-09 +* 2nd coating + MIX 3 300. + C12 = 'C0_GR' 9.37591E-02 + B10 = 'B10' 3.31664E-09 + B11 = 'B11' 1.33499E-08 +* Graphite matrix (Compact) + MIX 4 300. + C12 = 'C0_GR' 8.65629E-02 + H1 = 'H1' 3.23272E-05 + O16 = 'O16' 1.61635E-05 + B10 = 'B10' 3.06284E-09 + B11 = 'B11' 1.23284E-08 +* Graphite sheath + MIX 5 300. + C12 = 'C0_GR' 8.39319E-02 + B10 = 'B10' 1.29894E-09 + B11 = 'B11' 5.22841E-09 + H1 = 'H1' 1.23109E-05 + O16 = 'O16' 8.77944E-06 + N14 = 'N14' 9.77961E-06 +* air + MIX 6 300. + O16 = 'O16' 1.029640E-5 + N14 = 'N14' 3.837430E-5 +* Block graphite + MIX 7 300. + C12 = 'C0_GR' 8.34807E-02 + B10 = 'B10' 1.29196E-09 + B11 = 'B11' 5.20030E-09 + H1 = 'H1' 1.22447E-05 + O16 = 'O16' 8.78735E-06 + N14 = 'N14' 9.93241E-06 +; +*-------------------------------------------------------------------- +* Geometry: GEOMETRY - Infinite hexagonal lattice of fuel pin +*-------------------------------------------------------------------- +GEOMETRY := GEO: :: HEX 1 + HBC COMPLETE REFL + CELL B2 + +* SPHERICAL MICRO-STRUCTURES * +* NMISTR = NUMBER OF MICRO STRUCTURE TYPES IN REGION +* NMILG = NUMBER OF MICRO STRUCTURE REGIONS +* NS = ARRAY OF SUB REGIONS IN THE MICRO STRUCTURES; LEN=NMILG +* RS = RADIUS OF MICROS STRUCTURES; LEN=NS(I); I=1,NMISTR +* milie = COMPOSITION OF EACH shell OF MICRO STRUCTURE; LEN=NMISTR; +* !!!!!!!! milie #'s ARE > THAN MIX #'S !!!!!!!! +* mixdil = BASE COMPOSITION OF EACH REGION; LEN=NMILG +* fract = MICROSTRUCTURE TYPE VOLUME FRACTION IN REGION LEN=NMILG +* mixgr = LIBRARY MIXTURES FOR EACH MICROSTRUCTURE SHELL; LEN=NS(I) + +* NMISTR, NMILG + BIHET SPHE 1 1 + (* NS *) 3 + (* RS *) 0.0 .0301 .0380 .04590 + (* milie *) 8 + (* mixdil *) 4 + (* loop over NMISTR *) + (* fract *) 0.2977 + (* mixgr *) 1 2 3 + + ::: B2 := GEO: HEXCEL 5 + SIDE 3.75278 + RADIUS 0.0 0.8975 1.7925 1.8250 2.3400 2.3600 + MIX 6 8 6 5 6 7 ; +; +*---- +* +* Use EXCELT, MCCGT, USS, ASM and FLU for flux solution +*---- +TRACK TRACK_FIL := EXCELT: GEOMETRY :: + EDIT 1 + TITLE 'VHTRC B2 FUEL PIN WITH REPRESENTATIVE GRAPHITE ' + MAXR 500000 + ANIS 1 + ALLG + QUAB 5 SAPO + PISO + TRAK TISO 8 20.0 +; +TRACK := MCCGT: TRACK TRACK_FIL GEOMETRY :: + EDIT 1 LCMD 3 + EPSI 5E-5 + MAXI 100 KRYL 10 SCR 1 AAC 1 TMT +; +*--- +* Self Shielding calculation +*--- +LIBRARY2 := USS: LIBRARY TRACK TRACK_FIL :: + EDIT 2 + PASS 2 + ARM +; +*--- +* Multigroup Flux Calculation +*--- + SYS := ASM: LIBRARY2 TRACK TRACK_FIL :: + EDIT 1 + ARM + ; + FLUX := FLU: SYS LIBRARY2 TRACK TRACK_FIL :: + EDIT 1 + TYPE K + EXTE 500 1E-5 + THER 10 1E-5 ; +assertS FLUX :: 'K-INFINITY' 1 1.208958 ; +ECHO "test TCDR19 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR20.c2m b/Dragon/data/tdraglib_proc/TCDR20.c2m new file mode 100644 index 0000000..772c061 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR20.c2m @@ -0,0 +1,175 @@ +*---- +* +* TEST CASE TCDR20 +* 1D SN CASE WITH SPH/P5 EQUIVALENCE +* 172-GROUP JEF2.2 DRAGLIB +* +* Author: J. Ortensi INL +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY LIBRARY2 SYS GEOM SPHGEOM MACRO EDIT + TRACK SPHTRACK FLUX ; +MODULE LIB: GEO: SNT: USS: BIVACT: ASM: FLU: EDI: DELETE: + SPH: END: ; +SEQ_ASCII MAC_P3SPH :: FILE 'MAC_P3SPH.txt' ; +REAL Tfuel NB10 NB11 := 300.0 1.4495E-07 5.8343E-07 ; +PROCEDURE assertS ; +*---- +LIBRARY := LIB: :: +EDIT 0 +NMIX 5 +DEPL LIB: DRAGON FIL: DLIB_J2 +CTRA NONE +ANIS 2 +SUBG +MIXS LIB: DRAGON FIL: DLIB_J2 +*-- Fuel +MIX 1 <<Tfuel>> +GRA = 'C0_GR' 3.5353E-02 +C0 = 'C0' 5.0874E-02 +B10 = 'B10' <<NB10>> +B11 = 'B11' <<NB11>> +U235 = 'U235' 8.6849E-06 1 IRSET 0.0 81 +U238 = 'U238' 6.2967E-07 1 IRSET 0.0 81 +Fe54 = 'Fe54' 6.5286E-07 +Fe56 = 'Fe56' 1.0239E-05 +Fe57 = 'Fe57' 2.3659E-07 +Fe58 = 'Fe58' 3.1248E-08 +O16 = 'O16' 1.8623E-05 +*-- Zircaloy Clad +MIX 2 <<Tfuel>> +Sn115 = 'Sn115' 3.1545E-07 +Sn116 = 'Sn116' 1.3490E-05 +Sn117 = 'Sn117' 7.1255E-06 +Sn118 = 'Sn118' 2.2471E-05 +Sn119 = 'Sn119' 7.9698E-06 +Sn120 = 'Sn120' 3.0228E-05 +Sn122 = 'Sn122' 4.2957E-06 +Sn124 = 'Sn124' 5.3720E-06 +Zr0 = 'Zr0' 4.2049E-02 2 IRSET 0.0 81 +Fe54 = 'Fe54' 1.1659E-05 +Fe56 = 'Fe56' 1.8285E-04 +Fe57 = 'Fe57' 4.2252E-06 +Fe58 = 'Fe58' 5.5804E-07 +B10 = 'B10' 1.3751E-07 +B11 = 'B11' 5.9006E-07 +Cr50 = 'Cr50' 9.8706E-08 +Cr52 = 'Cr52' 1.9013E-06 +Cr53 = 'Cr53' 2.1556E-07 +Cr54 = 'Cr54' 5.3551E-08 +Ni58 = 'Ni58' 2.2803E-07 +Ni60 = 'Ni60' 8.7821E-08 +Ni61 = 'Ni61' 3.8183E-09 +Ni62 = 'Ni62' 2.2207E-08 +Ni64 = 'Ni64' 3.1149E-09 +Cd110 = 'Cd110' 8.7398E-10 +Cd111 = 'Cd111' 8.9567E-10 +Cd112 = 'Cd112' 1.6885E-09 +Cd113 = 'Cd113' 8.5508E-10 +Cd114 = 'Cd114' 2.0104E-09 +Cd116 = 'Cd116' 5.2411E-10 +*-- Air +MIX 3 <<Tfuel>> +C0 = 'C0' 7.5811E-09 +N14 = 'N14' 3.9484E-05 +O16 = 'O16' 1.0608E-05 +*-- Reflector +MIX 4 <<Tfuel>> +GRA = 'C0_GR' 4.936E-02 +C0 = 'C0' 3.430E-02 +B10 = 'B10' 3.517E-08 +B11 = 'B11' 1.509E-07 +Fe54 = 'Fe54' 1.054E-06 +Fe56 = 'Fe56' 1.652E-05 +Fe57 = 'Fe57' 3.818E-07 +Fe58 = 'Fe58' 5.043E-08 +*-- Aluminum +MIX 5 <<Tfuel>> +Al27 = 'Al27' 5.9477E-02 +Fe54 = 'Fe54' 2.5591E-05 +Fe56 = 'Fe56' 4.0136E-04 +Fe57 = 'Fe57' 9.2739E-06 +Fe58 = 'Fe58' 1.2249E-06 +; + +************************************************************************ +* build and track geometry +************************************************************************ +* Fuel Zr void Fuel ref Al void ref +* DX 26.913 1.402 5.732 26.913 26.642 2.789 4.886 26.642 +************************************************************************ + +*-- Single block model +GEOM := GEO: :: CAR1D 8 +X- REFL X+ REFL +MESHX 0.0 26.913 28.315 34.047 60.960 87.602 90.391 95.278 121.920 +MIX 1 2 3 1 4 5 3 4 +SPLITX 20 1 5 20 20 2 4 20 ; + +*-- Single block model +SPHGEOM := GEO: :: CAR1D 8 +X- REFL X+ REFL +MESHX 0.0 26.913 28.315 34.047 60.960 87.602 90.391 95.278 121.920 +MIX 1 2 3 1 4 5 3 4 +SPLITX 8 1 1 8 8 1 1 8 ; + +TRACK := SNT: GEOM :: + EDIT 0 + MAXR 15000 + DIAM 0 + SN 12 + SCAT 2 +; + +SPHTRACK := BIVACT: SPHGEOM :: +EDIT 2 +SPN 5 +* SPH with SCAT > 1 not implemented for PN +SCAT 1 +MAXR 500 +DUAL (*IELEM=*) 1 (*ICOL=*) 1 ; + +************************************************************************ + +LIBRARY2 := USS: LIBRARY TRACK :: + EDIT 0 + PASS 2 + MAXST 100 + ARM ; + +SYS := ASM: LIBRARY2 TRACK :: + EDIT 0 ARM ; + +FLUX := FLU: SYS LIBRARY2 TRACK :: + EDIT 1 + TYPE K P1 ; +SYS := DELETE: SYS ; + +EDIT := EDI: FLUX LIBRARY2 TRACK GEOM :: + EDIT 1 + MERG MIX 1 2 3 4 5 + COND 2.E+05 9.5E+02 2.3E+00 6.0E-01 3.0E-01 + P0W MICRO RES + SAVE ; +FLUX := DELETE: FLUX ; + +MAC_P3SPH := EDIT :: STEP UP 'REF-CASE0001' STEP UP MACROLIB ; + +EDIT := SPH: EDIT SPHTRACK :: + EDIT 2 + MAXNB 200 + ITER 400 1E-5 ; +* +************************** +* Verification calculation +************************** +* +MACRO := EDIT :: STEP UP 'REF-CASE0001' STEP UP MACROLIB ; +SYS := ASM: MACRO SPHTRACK :: ARM ; +FLUX := FLU: MACRO SPHTRACK SYS :: TYPE K EXTE 100 5.0E-8 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.503346 ; +ECHO "test TCDR20 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR21.c2m b/Dragon/data/tdraglib_proc/TCDR21.c2m new file mode 100644 index 0000000..ab526c2 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR21.c2m @@ -0,0 +1,145 @@ +*---- +* TEST CASE TRDR21 +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* DIFFUSION COEFFICIENTS COMPUTED WITH P0TR EQUATION +* 172-GROUP JEF2.2 DRAGLIB +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + ASSMB DISCR LIBRARY CP FLUX EDITION MACDIR + HOMXS_EDI HOMXS_MIC HOM_GEOM HOM_DISCR HOM_CP FU_FLUX ; +MODULE + LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: SPH: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA OLDW + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 + H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2 + B10 = B10 4.71444E-6 B11 = B11 1.90959E-5 + MIX 2 579.9 O16 = O16 3.06711E-4 + Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4 + Zr91 = Zr0 4.18621E-2 + MIX 3 579.9 + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + MoNat = Mo95 4.03755E-6 + Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr0 8.92427E-4 + MIX 4 579.9 O16 = O16 2.87335E-4 + Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4 + Zr91 = Zr0 3.92175E-2 + MIX 5 579.9 + H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2 + Cr52 = Cr52 2.24991E-5 + Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5 + MoNat = Mo95 1.89281E-6 + Al27 = Al27 1.10277E-6 + Mn55 = Mn55 1.94976E-7 + B10 = B10 4.66484E-6 B11 = B11 1.88950E-5 + Zr91 = Zr0 4.18372E-4 + MIX 6 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + Cr52 = Cr52 2.11122E-5 + Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5 + MoNat = Mo95 1.77614E-6 + Al27 = Al27 1.03479E-6 + Mn55 = Mn55 1.82957E-7 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + Zr91 = Zr0 3.92583E-4 + MIX 7 579.9 + H1H2O = H1_H2O 4.72020E-2 O16H2O = O16 2.36010E-2 + Cr52 = Cr52 1.96591E-5 + Fe56 = Fe56 1.82630E-5 Ni58 = Ni58 4.65011E-5 + MoNat = Mo95 1.65389E-6 + Al27 = Al27 9.63569E-7 + Mn55 = Mn55 1.70365E-7 + B10 = B10 4.67110E-6 B10 = B10 1.89203E-5 + Zr91 = Zr0 3.65562E-4 + MIX 8 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + ; +*---- +* Geometry ASSMB : a 17 X 17 normal PWR assembly +* contains C1 : cell without fuel +* C2 : normal fuel cell +* C3 : peripheral cell +* C4 : corner cell +*---- +ASSMB := GEO: :: CAR2D 9 9 + X- DIAG X+ REFL Y- SYME Y+ DIAG + CELL C1 C2 C2 C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C2 C2 C2 C3 + C2 C2 C2 C2 C2 C2 C3 + C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C3 + C1 C2 C2 C3 + C2 C2 C3 + C2 C3 + C4 + MERGE 1 3 12 11 12 12 11 12 15 + 4 6 5 6 6 5 6 8 + 13 5 6 6 5 6 8 + 2 5 5 10 5 8 + 13 5 5 6 8 + 2 5 7 8 + 13 7 8 + 14 8 + 9 + TURN H H B H H B H H A + H G G H G G H A + A E E F E E A + H H F H H A + H E G H A + H H A A + H A A + A A + A + ::: C1 := GEO: CARCEL 2 + MESHX 0.0 1.26472 MESHY 0.0 1.26472 + RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ; + ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ; + ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ; + ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ; + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for keff no leakage +* Diffusion coefficients with P0TR equation +*---- +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; +LIBRARY := SHI: LIBRARY DISCR :: LJ ; +CP := ASM: LIBRARY DISCR :: ARM ; +FLUX := FLU: CP LIBRARY DISCR :: EDIT 2 TYPE K P0TR PNL ; +EDITION := EDI: LIBRARY DISCR FLUX ASSMB :: + EDIT 4 SAVE MICR RES MERGE COMP COND 0.625 JOUT ; +EDITION := SPH: EDITION :: EDIT 2 ALBS SELE_ALB ; +* +MACDIR := EDITION :: STEP UP 'REF-CASE0001' STEP UP MACROLIB + STEP UP GROUP STEP AT 1 ; +assertS MACDIR :: 'NTOT0' 1 5.363227E-01 ; +assertS MACDIR :: 'NUSIGF' 1 6.825784E-03 ; +assertS MACDIR :: 'DIFF' 1 1.532271E+00 ; +assertS MACDIR :: 'SIGW00' 1 5.097417E-01 ; +MACDIR := DELETE: MACDIR ; +MACDIR := EDITION :: STEP UP 'REF-CASE0001' STEP UP MACROLIB + STEP UP GROUP STEP AT 2 ; +assertS MACDIR :: 'NTOT0' 1 1.338058E+00 ; +assertS MACDIR :: 'NUSIGF' 1 1.413793E-01 ; +assertS MACDIR :: 'DIFF' 1 3.821711E-01 ; +assertS MACDIR :: 'SIGW00' 1 1.247216E+00 ; +MACDIR := DELETE: MACDIR ; +ECHO "test TCDR21 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR22.c2m b/Dragon/data/tdraglib_proc/TCDR22.c2m new file mode 100644 index 0000000..8983c9f --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR22.c2m @@ -0,0 +1,208 @@ +*---- +* TEST CASE TRDR22 +* 17 X 17 PWR VOIDED ASSEMBLY +* 172-GROUP JEF2.2 DRAGLIB +* Author : Alain Hebert (2019) +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + ASSMB DISCR LIBRARY CP CALC OUT SPHGEOM MTRACK MACRO MACRO2 SYS FLUX + EDIT ; +MODULE + LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: FLUD: MAC: DELETE: UTL: BIVACT: + BIVACA: SPH: OUT: ERROR: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS ANIS 2 + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 + H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2 + B10 = B10 4.71444E-4 B11 = B11 1.90959E-3 + MIX 2 579.9 O16 = O16 3.06711E-4 + Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4 + Zr91 = Zr0 4.18621E-2 + MIX 3 579.9 + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + MoNat = Mo95 4.03755E-6 + Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr0 8.92427E-4 + MIX 4 579.9 O16 = O16 2.87335E-4 + Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4 + Zr91 = Zr0 3.92175E-2 + MIX 5 579.9 + H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2 + Cr52 = Cr52 2.24991E-5 + Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5 + MoNat = Mo95 1.89281E-6 + Al27 = Al27 1.10277E-6 + Mn55 = Mn55 1.94976E-7 + B10 = B10 4.66484E-6 B11 = B11 1.88950E-5 + Zr91 = Zr0 4.18372E-4 + MIX 6 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + Cr52 = Cr52 2.11122E-5 + Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5 + MoNat = Mo95 1.77614E-6 + Al27 = Al27 1.03479E-6 + Mn55 = Mn55 1.82957E-7 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + Zr91 = Zr0 3.92583E-4 + MIX 7 579.9 + H1H2O = H1_H2O 4.72020E-2 O16H2O = O16 2.36010E-2 + Cr52 = Cr52 1.96591E-5 + Fe56 = Fe56 1.82630E-5 Ni58 = Ni58 4.65011E-5 + MoNat = Mo95 1.65389E-6 + Al27 = Al27 9.63569E-7 + Mn55 = Mn55 1.70365E-7 + B10 = B10 4.67110E-6 B10 = B10 1.89203E-5 + Zr91 = Zr0 3.65562E-4 + MIX 8 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + ; +*---- +* Geometry ASSMB : a 17 X 17 normal PWR assembly +* contains C1 : cell without fuel +* C2 : normal fuel cell +* C3 : peripheral cell +* C4 : corner cell +*---- +ASSMB := GEO: :: CAR2D 9 9 + X- DIAG X+ REFL Y- SYME Y+ DIAG + CELL C1 C2 C2 C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C2 C2 C2 C3 + C2 C2 C2 C2 C2 C2 C3 + C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C3 + C1 C2 C2 C3 + C2 C2 C3 + C2 C3 + C4 + MERGE 1 3 12 11 12 12 11 12 15 + 4 6 5 6 6 5 6 8 + 13 5 6 6 5 6 8 + 2 5 5 10 5 8 + 13 5 5 6 8 + 2 5 7 8 + 13 7 8 + 14 8 + 9 + TURN H H B H H B H H A + H G G H G G H A + A E E F E E A + H H F H H A + H E G H A + H H A A + H A A + A A + A + ::: C1 := GEO: CARCEL 2 + MESHX 0.0 1.26472 MESHY 0.0 1.26472 + RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ; + ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ; + ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ; + ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ; + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for P1 homogeneous leakage +*---- +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; +LIBRARY := SHI: LIBRARY DISCR :: LJ ; +ASSMB DISCR := DELETE: ASSMB DISCR ; +*---- +* New lattice geometry with VOID boundary condition +*---- +ASSMB := GEO: :: CAR2D 9 9 + X- DIAG X+ VOID Y- SYME Y+ DIAG + CELL C1 C2 C2 C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C2 C2 C2 C3 + C2 C2 C2 C2 C2 C2 C3 + C1 C2 C2 C1 C2 C3 + C2 C2 C2 C2 C3 + C1 C2 C2 C3 + C2 C2 C3 + C2 C3 + C4 + ::: C1 := GEO: CARCEL 2 + MESHX 0.0 1.26472 MESHY 0.0 1.26472 + RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ; + ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ; + ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ; + ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ; + ; +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; +CP := ASM: LIBRARY DISCR :: ARM ; +CALC := FLU: CP LIBRARY DISCR :: + TYPE K P1 HETE ; +assertS CALC :: 'K-EFFECTIVE' 1 0.4551677 ; + +*---- +* Editing using SPH model for transport-diffusion with flux-volume +* normalization +*---- +OUT := EDI: LIBRARY DISCR CALC ASSMB :: + EDIT 4 UPS SAVE MICR RES MERGE CELL COND 4.0 + ; + +*---- +* TRANSPORT-DIFFUSION SPH EQUIVALENCE IN NON-FUNDAMENTAL MODE +*---- +SPHGEOM := GEO: :: CAR2D 9 9 + X- DIAG X+ ALBE 1 Y- SYME Y+ DIAG + MIX 1 2 3 4 5 6 7 8 9 + 10 11 12 13 14 15 16 17 + 18 19 20 21 22 23 24 + 25 26 27 28 29 30 + 31 32 33 34 35 + 36 37 38 39 + 40 41 42 + 43 44 + 45 + MESHX 0.0 1.26471996 2.52943993 3.79415989 5.05887985 + 6.32359982 7.58831978 8.85303974 10.1177597 11.4324799 + MESHY 0.0 1.26471996 2.52943993 3.79415989 5.05887985 + 6.32359982 7.58831978 8.85303974 10.1177597 11.4324799 +; +MTRACK := BIVACT: SPHGEOM :: DUAL (*IELEM=*) 1 (*ICOL=*) 1 ; +MACRO := OUT :: STEP UP 'REF-CASE0001' STEP UP 'MACROLIB' ; +*---- +* Set physical albedos +*---- +MACRO := MAC: MACRO :: ALBP 1 0.0 0.0 ; +*---- +* Perform SPH equivalence +*---- +MACRO2 := SPH: MACRO MTRACK :: EDIT 2 ITER 1.0E-7 ; +*---- +* Verification calculation +*---- +SYS := BIVACA: MACRO2 MTRACK ; +FLUX := FLUD: SYS MTRACK :: EDIT 2 ADI 20 ACCE 5 3 ; +EDIT := OUT: FLUX SPHGEOM MACRO2 MTRACK :: + EDIT 0 INTG MIX ; + +ECHO "Error statistics after correction" ; +ERROR: MACRO EDIT :: HREA 'NFTOT' ; + +assertS FLUX :: K-EFFECTIVE 1 0.4551687 ; + +ECHO "test TCDR22 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR31.c2m b/Dragon/data/tdraglib_proc/TCDR31.c2m new file mode 100644 index 0000000..23530bd --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR31.c2m @@ -0,0 +1,172 @@ +*---- +* TEST CASE TCDR31 +* CANDU-6 ANNULAR CELL +* 172-GROUP JEF2.2 DRAGLIB +* MULTI-PARAMETER COMPO ACCESS FOR MACRO-DEPLETION +* TWO GROUP BURNUP +* POWER (KW) = 615.00000 +* BURN POWER (KW/KG) = 31.97130 +* URANIUM MASS = 19.23600 +* UO2 REAL DENSITY = 10.59300 +* UO2 EFF DENSITY = 10.43750 +* UO2 TEMPERATURE = 941.28998 +* ENRICHMENT = 0.71140 +* COOLANT D2 AT % = 99.222 +* MODERATOR D2 AT % = 99.911 +* NUMBER OF DAYS = 50 +* +*---- +* Define variables +* Burnup paremeters +* a) Power +* = 31.9713 kw/kg for 0.0 to 300.0 days +* b) 69 Groups Burnup time interval Delt +* = 300 day for 0 to 300 day +* c) 2 Groups Burnup time interval Delt +* = 1 day for 0 to 1 day +* = 4 days for 1 to 5 days +* = 5 days for 5 to 10 days +* = 10 days for 10 to 50 days +* = 20 days for 50 to 150 days +* = 50 days for 150 to 300 days +* c) Days with burnup interval changes +* = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days +* d) Burnup control time variables Timei, Timef +* Timei = initial time +* Timef = final time +*---- +REAL + Power Delt Timec Timei Timef := + 31.9713 1.0 300.0 0.0 0.0 ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY ; +LINKED_LIST + CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION DATABASE ; +SEQ_BINARY + INTLINS INTLINF ; +SEQ_ASCII + res ; +MODULE + GEO: EXCELT: LIB: SHI: ASM: FLU: EVO: EDI: COMPO: SPH: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from Jef 2.2 Draglib +* Microscopic cross sections from Jef 2.2 Draglib +*---- +PROCEDURE TCDR05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCDR05Lib :: <<iedit>> ; +*---- +* Geometry CANDU6S : 13 regions annular cluster for self-shielding +* CANDU6F : 31 regions annular cluster for transport +*---- +CANDU6S := GEO: :: TUBE 5 + R+ REFL RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 16.12171 + MIX 1 2 3 4 5 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN 1 RPIN 0.0000 APIN 0.0000 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2 := GEO: ROD1 MIX 7 10 NPIN 6 RPIN 1.4885 APIN 0.0000 ; + ::: ROD3 := GEO: ROD1 MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ; + ; +CANDU6F := GEO: CANDU6S :: SPLITR 6 1 1 1 10 + ::: ROD1 := GEO: ROD1 SPLITR 2 1 ; + ::: ROD2 := GEO: ROD2 SPLITR 2 1 ; + ::: ROD3 := GEO: ROD3 SPLITR 2 1 ; + ::: ROD4 := GEO: ROD4 SPLITR 2 1 ; + ; +*---- +* Self-Shielding calculation EXCEL +* Transport calculation EXCEL +* Flux calculation for keff +*---- +VOLMATS INTLINS := EXCELT: CANDU6S :: + TITLE 'TCDR31: FEW GROUP BURNUP / SELF-SHIELDING TRACKING' + EDIT 0 MAXR 13 TRAK TISO 5 10.0 SYMM 12 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; +VOLMATF INTLINF := EXCELT: CANDU6F :: + TITLE 'TCDR31: FEW GROUP BURNUP / TRANSPORT TRACKING' + EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.117541 ; +EDITION := EDI: FLUX LIBRARY VOLMATF CANDU6F :: + MERG REGI 6 6 10 7 7 10 1 1 8 8 10 1 + 1 9 9 10 1 1 2 3 4 5 5 5 + 5 5 5 5 5 5 5 + COND 4.0 MICR ALL SAVE + MGEO CANDU6F + ; +EDITION := SPH: EDITION VOLMATF INTLINF ; +FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S VOLMATF CANDU6F +INTLINF := DELETE: FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S +VOLMATF CANDU6F INTLINF ; +*---- +* 2 group Burnup (macro-depletion) +*---- +DATABASE := COMPO: :: + EDIT 5 + COMM 'Multi-parameter reactor database' ENDC + INIT + ; +DATABASE := COMPO: DATABASE EDITION :: + EDIT 3 + ; +EDITION := DELETE: EDITION ; +* +* COMPO INTERPOLATION +LIBRARY := SPH: DATABASE :: + EDIT 1 + STEP UP default STEP AT 1 + MICRO OFF + ; +CANDU6F := DATABASE :: STEP UP default STEP UP 'GEOMETRIES' STEP AT 1 ; +VOLMATF INTLINF := EXCELT: CANDU6F :: + TITLE 'TCDR31: FEW GROUP BURNUP / TRANSPORT TRACKING' + EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.117542 ; +EDITION := EDI: FLUX LIBRARY VOLMATF :: + EDIT 1 MERGE COMP MICR ALL SAVE ; +EVALUATE Timec := 1.0 ; +WHILE Timei Timec < DO + EVALUATE Timef := Timei Delt + ; + IF Timei 0.0 = THEN + BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF :: + EDIT 3 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + EDIT 3 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ; + ENDIF ; + PIJ := DELETE: PIJ ; + PIJ := ASM: LIBRARY VOLMATF INTLINF ; + FLUX := FLU: FLUX PIJ LIBRARY VOLMATF :: + TYPE K ; + EDITION := EDI: EDITION FLUX LIBRARY VOLMATF :: + SAVE ; +*---- +* change delta t for burnup and final time if required +*---- + IF Timef Timec = THEN + 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 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.075418 ; +res := EDITION ; +EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE: + EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ; +ECHO "test TCDR31 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR32.c2m b/Dragon/data/tdraglib_proc/TCDR32.c2m new file mode 100644 index 0000000..c0653b0 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR32.c2m @@ -0,0 +1,207 @@ +*---- +* +* TEST CASE TCDR32 +* FIRST MOX ROWLAND'S BENCHMARK +* DISTRIBUTED SELF-SHIELDING +* 172-GROUP JEF2.2 DRAGLIB +* +* Author: A. Hebert +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ; +MODULE LIB: GEO: SYBILT: SHI: USS: ASM: FLU: DELETE: END: ; +INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1 + MODE1 MODE2 := + 1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *) +REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ; +PROCEDURE assertS ; +* +EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.259307 +0.343031 0.366715 0.388960 0.3996186 0.41 ; +GEOM := GEO: :: TUBE 9 + R+ REFL + RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>> + <<RCOMB5>> <<RCOMB6>> 0.475 0.6045552 0.710879 + MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>> + <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>> + ; +TRACK := SYBILT: GEOM :: + TITLE 'FIRST MOX ROWLAND S BENCHMARK' + MAXR 20 QUA1 5 ; +*---- +* SHI LEVEL 1 +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 + ADED 4 NELAS N4N N2N N3N +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX <<COMB0101>> 300.0 (*COMB0101*) + U235 = U235 5.105E-5 1 IRSET 0.0 81 + U238 = U238 2.037E-2 1 IRSET 0.0 81 + PU238 = Pu238 4.669E-5 1 IRSET 0.0 81 + PU239 = Pu239 1.465E-3 1 IRSET 0.0 81 + PU240 = Pu240 5.691E-4 1 IRSET 0.0 81 + PU241 = Pu241 2.713E-4 1 IRSET 0.0 81 + PU242 = Pu242 1.413E-4 1 IRSET 0.0 81 + AM241 = Am241 3.028E-5 1 IRSET 0.0 81 + O16 = O16 4.588E-2 + MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) + MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) + MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) + MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) + MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*) + MIX <<GAIN1>> 300.0 (*GAIN1*) + ZR = Zr0 3.88E-2 2 IRSET 0.0 81 + MIX <<MODE1>> 293.0 (*MODE1*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + MIX <<MODE2>> 293.0 (*MODE2*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + ; +LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 45 LEVEL 1 MXIT 3 ; +CP := ASM: LIBRARY2 TRACK :: PIJ ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.212449 ; +LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ; +*---- +* SHI LEVEL 2 +*---- +LIBRARY2 := SHI: LIBRARY TRACK :: EDIT 1 LJ GRMIN 45 LEVEL 2 MXIT 3 ; +CP := ASM: LIBRARY2 TRACK :: PIJ ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.217413 ; +LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ; +*---- +* USS +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 + ADED 4 NELAS N4N N2N N3N + SUBG (*HELIOS TYPE PROBABILITY TABLES*) +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX <<COMB0101>> 300.0 (*COMB0101*) + U235 = U235 5.105E-5 1 IRSET 0.0 81 + U238 = U238 2.037E-2 1 IRSET 0.0 81 + PU238 = Pu238 4.669E-5 1 IRSET 0.0 81 + PU239 = Pu239 1.465E-3 1 IRSET 0.0 81 + PU240 = Pu240 5.691E-4 1 IRSET 0.0 81 + PU241 = Pu241 2.713E-4 1 IRSET 0.0 81 + PU242 = Pu242 1.413E-4 1 IRSET 0.0 81 + AM241 = Am241 3.028E-5 1 IRSET 0.0 81 + O16 = O16 4.588E-2 + MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) + MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) + MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) + MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) + MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*) + MIX <<GAIN1>> 300.0 (*GAIN1*) + ZR = Zr0 3.88E-2 2 IRSET 0.0 81 + MIX <<MODE1>> 293.0 (*MODE1*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + MIX <<MODE2>> 293.0 (*MODE2*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + ; +LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 45 ; +CP := ASM: LIBRARY2 TRACK :: PIJ ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.216461 ; +LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ; +*---- +* Ribon extended +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 + ADED 4 NELAS N4N N2N N3N + PTSL (*CALENDF TYPE PROBABILITY TABLES*) +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX <<COMB0101>> 300.0 (*COMB0101*) + U235 = U235 5.105E-5 1 IRSET PTSL 1 + U238 = U238 2.037E-2 1 IRSET PTSL 1 + PU238 = Pu238 4.669E-5 1 IRSET PTSL 1 + PU239 = Pu239 1.465E-3 1 IRSET PTSL 1 + PU240 = Pu240 5.691E-4 1 IRSET PTSL 1 + PU241 = Pu241 2.713E-4 1 IRSET PTSL 1 + PU242 = Pu242 1.413E-4 1 IRSET PTSL 1 + AM241 = Am241 3.028E-5 1 IRSET PTSL 1 + O16 = O16 4.588E-2 + MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) + MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) + MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) + MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) + MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*) + MIX <<GAIN1>> 300.0 (*GAIN1*) + ZR = Zr0 3.88E-2 2 IRSET PTSL 1 + MIX <<MODE1>> 293.0 (*MODE1*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + MIX <<MODE2>> 293.0 (*MODE2*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + ; +LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 3 GRMIN 45 ; +CP := ASM: LIBRARY2 TRACK :: PIJ ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.216642 ; +LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ; +*---- +* Ribon extended with correlation +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 + ADED 4 NELAS N4N N2N N3N + PTSL (*CALENDF TYPE PROBABILITY TABLES*) +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX <<COMB0101>> 300.0 (*COMB0101*) + PU240 = Pu240 5.691E-4 CORR 1 IRSET PTSL 1 + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX <<COMB0101>> 300.0 (*COMB0101*) + U235 = U235 5.105E-5 1 IRSET PTSL 1 + U238 = U238 2.037E-2 CORR 1 IRSET PTSL 1 + PU238 = Pu238 4.669E-5 1 IRSET PTSL 1 + PU239 = Pu239 1.465E-3 1 IRSET PTSL 1 + PU241 = Pu241 2.713E-4 1 IRSET PTSL 1 + PU242 = Pu242 1.413E-4 1 IRSET PTSL 1 + AM241 = Am241 3.028E-5 1 IRSET PTSL 1 + O16 = O16 4.588E-2 + MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) + MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) + MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) + MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) + MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*) + MIX <<GAIN1>> 300.0 (*GAIN1*) + ZR = Zr0 3.88E-2 2 IRSET PTSL 1 + MIX <<MODE1>> 293.0 (*MODE1*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + MIX <<MODE2>> 293.0 (*MODE2*) + H1 = H1_H2O 4.744E-2 + O16 = O16 2.372E-2 + ; +LIBRARY2 := USS: LIBRARY TRACK :: EDIT 1 TRAN PASS 4 GRMIN 45 ; +CP := ASM: LIBRARY2 TRACK :: PIJ ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.218970 ; +LIBRARY LIBRARY2 CP FLUX := DELETE: LIBRARY LIBRARY2 CP FLUX ; +ECHO "test TCDR32 completed" ; +END: ; 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" . diff --git a/Dragon/data/tdraglib_proc/TCDR34.c2m b/Dragon/data/tdraglib_proc/TCDR34.c2m new file mode 100644 index 0000000..062008f --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR34.c2m @@ -0,0 +1,65 @@ +*---- +* TEST CASE TCDR34 +* PWR UOX pin cell cylindrical / 293K / uss-6c / cylindrical +* heterogeneous leakage / isotropic streaming +* gamma(B,Sigma)=1 +* critical B square +*---- +LINKED_LIST LIBRARY LIBRARY2 C3 TRACK SYS FLUX EDIT MACDIR ; +MODULE GEO: SYBILT: USS: MAC: LIB: ASM: FLU: EDI: END: ; +INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1 + MODE1 MODE2 := 1 2 3 4 5 6 7 8 9 ; (* 6 RINGS *) +REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ; +PROCEDURE assertS ; +* +LIBRARY := LIB: :: + EDIT 1 + NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 + SUBG (*COMPUTE SUBGROUPS*) +* + MIXS LIB: DRAGON FIL: DLIB_J2 (*JEF2.2*) + MIX <<COMB0101>> 293.0 (*COMB0101*) + O16 = O16 4.6624E-2 + U235 = U235 7.0803E-4 1 IRSET 0.0 81 + U238 = U238 2.2604E-2 1 IRSET 0.0 81 + MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) + MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) + MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) + MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) + MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*) + MIX <<GAIN1>> 293.0 (*GAIN1*) + Zr0 = Zr0 4.3241E-2 2 IRSET 0.0 81 + MIX <<MODE1>> 293.0 (*MODE1*) + H1 = H1_H2O 6.6988E-2 + O16 = O16 3.3494E-2 + MIX <<MODE2>> 293.0 (*MODE2*) + H1 = H1_H2O 6.6988E-2 + O16 = O16 3.3494E-2 + ; +EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 +0.334664 0.3577709 0.3794733 0.3898718 0.40 ; +C3 := GEO: :: TUBE 9 + R+ REFL + RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>> + <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275 + MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>> + <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>> + ; +TRACK := SYBILT: C3 :: + EDIT 1 + MAXR 20 (*MAXIMUM OF 9 REGIONS*) + QUA1 5 TITLE 'SIMPLE WIGNER-SEITZ CELL' ; +LIBRARY2 := USS: LIBRARY TRACK :: EDIT 0 TRAN PASS 2 GRMIN 45 ; +SYS := ASM: LIBRARY2 TRACK :: PIJ ECCO ; +FLUX := FLU: SYS LIBRARY2 TRACK :: TYPE B P1 ECCO ; +assertS FLUX :: 'K-INFINITY' 1 1.378569 ; +EDIT := EDI: FLUX LIBRARY2 TRACK :: EDIT 3 + COND 75 99 MERGE COMP + SAVE ON 'EDITCDAT 1' ; +MACDIR := EDIT :: STEP UP 'EDITCDAT 1' STEP UP MACROLIB + STEP UP GROUP STEP AT 2 ; +assertS MACDIR :: 'DIFF' 1 0.64423 ; +ECHO "test TCDR34 completed" ; +END: ; diff --git a/Dragon/data/tdraglib_proc/TCDR35.c2m b/Dragon/data/tdraglib_proc/TCDR35.c2m new file mode 100644 index 0000000..7156fc0 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR35.c2m @@ -0,0 +1,79 @@ +*---- +* +* TEST CASE TCDR35 WITH SN SELF-SHIELDING AND FLUX SOLUTIONS +* FIRST UO2 ROWLAND'S BENCHMARK +* DISTRIBUTED SELF-SHIELDING +* 172-GROUP JEF2.2 DRAGLIB +* +* Author: A. Hebert +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY LIBRARY2 GEOM1 GEOM2 TRACK SYS FLUX ; +MODULE GEO: SNT: USS: MAC: LIB: ASM: FLU: DELETE: END: ; +INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1 + MODE1 MODE2 := + 1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *) +REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ; +PROCEDURE assertS ; +* +LIBRARY := LIB: :: + EDIT 3 + NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) + ANIS 2 + CTRA APOL + PTSL (*COMPUTE SUBGROUPS*) +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX <<COMB0101>> 293.0 (*COMB0101*) + O16 = O16 4.6624E-2 + U235 = U235 7.0803E-4 CORR 1 IRSET PTSL 1 + U238 = U238 2.2604E-2 CORR 1 IRSET PTSL 1 + MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) + MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) + MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) + MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) + MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*) + MIX <<GAIN1>> 293.0 (*GAIN1*) + Zr0 = Zr0 4.3241E-2 2 IRSET PTSL 1 + MIX <<MODE1>> 293.0 (*MODE1*) + H1 = H1_H2O 6.6988E-2 + O16 = O16 3.3494E-2 + MIX <<MODE2>> 293.0 (*MODE2*) + H1 = H1_H2O 6.6988E-2 + O16 = O16 3.3494E-2 + ; +EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 +0.334664 0.3577709 0.3794733 0.3898718 0.40 ; +GEOM1 := GEO: :: TUBE 9 + R+ REFL + RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>> + <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275 + MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>> + <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>> + SPLITR 1 1 1 1 1 1 1 -2 -2 + ; +GEOM2 := GEO: :: TUBE 9 + R+ REFL + RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>> + <<RCOMB5>> <<RCOMB6>> 0.45 0.5748331 0.6770275 + MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>> + <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>> + SPLITR 1 1 1 1 1 1 1 -8 -8 + ; +TRACK := SNT: GEOM1 :: + EDIT 1 + MAXR 50 (*MAXIMUM OF 50 REGIONS*) + SN 16 SCAT 1 QUAD 10 ; +LIBRARY2 := USS: LIBRARY TRACK :: EDIT 0 ARM PASS 4 GRMIN 45 ; +TRACK := DELETE: TRACK ; +TRACK := SNT: GEOM2 :: + EDIT 1 + MAXR 50 (*MAXIMUM OF 50 REGIONS*) + SN 16 SCAT 1 QUAD 10 ; +SYS := ASM: LIBRARY2 TRACK :: ARM ; +FLUX := FLU: SYS LIBRARY2 TRACK :: EDIT 1 TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.387591 ; +ECHO "test TCDR35 completed" ; +END: ; diff --git a/Dragon/data/tdraglib_proc/TCDR36.c2m b/Dragon/data/tdraglib_proc/TCDR36.c2m new file mode 100644 index 0000000..8a67373 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR36.c2m @@ -0,0 +1,337 @@ +*---- +* TEST CASE TCDR36 +* MULTICELL CARTESIAN ASSEMBLY WITH WATER HOLES +* 172-GROUP JEF2.2 DRAGLIB +* SPH with BIVAC / multicompo creation / verification calculation +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST GEOM_F1 TRACK_F1 GEOM_SS TRACK_SS LIBRARY AUTO_LIB + SYS FLUX1 BURN EDITION DATABASE MACROLIB SPHGEOM MTRACK ; +MODULE DELETE: END: LIB: GEO: SYBILT: USS: ASM: FLU: EVO: EDI: + COMPO: BIVACT: GREP: SPH: ABORT: ; +REAL b2 ; +REAL VALUE REFVALUE DELTA ; +PROCEDURE assertS ; +* +STRING LibName := "DLIB_J2" ; +INTEGER F11 F12 F13 F14 F15 := 1 2 3 4 5 ; +INTEGER GAP CLAD MODE := 6 7 8 ; +* +REAL COTE := 1.26 ; +REAL RF1 RF2 RF3 RF4 RF5 := 0.1767 0.2499 0.3060 0.3534 0.3951 ; +REAL RG RC := 0.4010 0.4583 ; +REAL Tmode Tclad Tfuel Tgap := 580.0 600.0 900.0 900.0 ; +REAL Keff ; +* +************************************************************************ +* Zero burnup Calculation -- Zero burnup Calculation -- Zero burnup Calc +************************************************************************ +* +* Microscopic XS Library *********************************************** +* +LIBRARY := LIB: :: + EDIT 0 + NMIX <<MODE>> + CTRA APOL + ANIS 2 + SUBG + DEPL LIB: DRAGON FIL: <<LibName>> + MIXS LIB: DRAGON FIL: <<LibName>> +* + MIX <<F11>> <<Tfuel>> (*F11*) + U234 0.0 1 IRSET 0.0 81 + U235 9.71492E-04 1 IRSET 0.0 81 + U236 0.0 1 IRSET 0.0 81 + U238 2.18794E-02 1 IRSET 0.0 81 + Np237 0.0 1 IRSET 0.0 81 + Pu238 0.0 1 IRSET 0.0 81 + Pu239 0.0 1 IRSET 0.0 81 + Pu240 0.0 1 IRSET 0.0 81 + Pu241 0.0 1 IRSET 0.0 81 + Pu242 0.0 1 IRSET 0.0 81 + Am241 0.0 1 IRSET 0.0 81 + O16 4.57018E-02 + MIX <<F12>> COMB <<F11>> 1.0 (*F12*) + MIX <<F13>> COMB <<F11>> 1.0 (*F13*) + MIX <<F14>> COMB <<F11>> 1.0 (*F14*) + MIX <<F15>> COMB <<F11>> 1.0 (*F15*) +* + MIX <<GAP>> <<Tgap>> (*GAP*) + O16 3.76661E-05 +* + MIX <<CLAD>> <<Tclad>> (*CLAD*) + N14 2.15407E-05 + Cr50 1.78447E-06 + Cr52 3.58179E-05 + Cr53 4.13810E-06 + Cr54 1.04843E-06 + Fe54 2.87539E-06 + Fe56 4.71718E-05 + Fe57 1.09881E-06 + Fe58 1.48904E-07 + Zr0 4.20495E-02 2 IRSET 0.0 81 + Sn115 2.11926E-06 + Sn116 9.11285E-05 + Sn117 4.85505E-05 + Sn118 1.54514E-04 + Sn119 5.52294E-05 + Sn120 2.11413E-04 + Sn122 3.05046E-05 + Sn124 3.87890E-05 +* + MIX <<MODE>> <<Tmode>> (*MODE*) + H1_H2O 4.76055E-02 NOEV + B10 7.89221E-06 NOEV + B11 3.17673E-05 NOEV + O16 2.38129E-02 NOEV + ; +* +* Self-Shielding ******************************************************* +*---- +* Geometry +*---- +GEOM_SS := GEO: :: CAR2D 3 3 + X- DIAG X+ REFL + Y- REFL Y+ DIAG + CELL C1 C2 C1 + C2 C2 + C1 + MERGE 1 2 1 + 2 3 + 1 + ::: C1 := GEO: CARCEL 0 (*WATER HOLE*) + MESHX 0.0 <<COTE>> + MESHY 0.0 <<COTE>> + MIX <<MODE>> ; + ::: C2 := GEO: CARCEL 7 (*STANDARD CELL*) + X- REFL X+ REFL + Y- REFL Y+ REFL + MESHX 0.0 <<COTE>> + MESHY 0.0 <<COTE>> + RADIUS 0.0 <<RF1>> <<RF2>> <<RF3>> <<RF4>> <<RF5>> + <<RG>> <<RC>> + MIX <<F11>> <<F12>> <<F13>> <<F14>> <<F15>> + <<GAP>> <<CLAD>> <<MODE>> ; + ; +*---- +* Tracking +*---- +TRACK_SS := SYBILT: GEOM_SS :: + EDIT 0 + MAXR 400 + MAXZ 100000 + DP01 + QUA2 20 3 + ; +*---- +* Self-Shielding +*---- +AUTO_LIB := USS: LIBRARY TRACK_SS :: + EDIT 0 GRMIN 45 PASS 2 + CALC REGI W1 U234 ALL + REGI W1 U235 ALL + REGI W1 U236 ALL + REGI W1 Np237 ALL + REGI W1 Pu238 ALL + REGI W1 Pu239 ALL + REGI W1 Pu240 ALL + REGI W1 Pu241 ALL + REGI W1 Pu242 ALL + REGI W1 Am241 ALL + REGI W1 Zr0 ALL + ENDC + ; +*---- +* Create the reactor database +*---- +DATABASE := COMPO: :: + EDIT 5 + STEP UP 'fuel' + COMM 'Multi-parameter reactor database for fuel' ENDC + PARA 'BURN' IRRA + PARA 'FLUB' FLUB + INIT + ; +*---- +* Geometry +*---- +GEOM_F1 := GEO: :: CAR2D 3 3 + X- DIAG X+ REFL + Y- REFL Y+ DIAG + CELL C1 C2 C1 + C2 C2 + C1 + MERGE 1 2 3 + 2 4 + 5 + ::: C1 := GEO: CARCEL 0 (*WATER HOLE*) + MESHX 0.0 <<COTE>> + MESHY 0.0 <<COTE>> + MIX <<MODE>> ; + ::: C2 := GEO: CARCEL 7 (*STANDARD CELL*) + MESHX 0.0 <<COTE>> + MESHY 0.0 <<COTE>> + RADIUS 0.0 <<RF1>> <<RF2>> <<RF3>> <<RF4>> <<RF5>> + <<RG>> <<RC>> + MIX <<F11>> <<F12>> <<F13>> <<F14>> <<F15>> + <<GAP>> <<CLAD>> <<MODE>> ; + ; +*---- +* Tracking +*---- +TRACK_F1 := SYBILT: GEOM_F1 :: + EDIT 0 + MAXR 400 + MAXZ 100000 + DP01 + QUA2 20 3 + ; +*---- +* Flux Calculation +*---- +SYS := ASM: AUTO_LIB TRACK_F1 :: ; +FLUX1 := FLU: AUTO_LIB TRACK_F1 SYS :: TYPE B B1 SIGS ; +SYS := DELETE: SYS ; +assertS FLUX1 :: 'K-INFINITY' 1 1.217577 ; +* +************************************************************************ +* burnup loop -- burnup loop -- burnup loop -- burnup loop -- burnup loo +************************************************************************ +* +REAL Power Delt ; +INTEGER i f := 2 1 ; +REAL + Timei Timef Delb Burnc Burni Burnf Burns := + 0.0 0.0 75.0 150.0 0.0 0.0 0.0 ; +EVALUATE Power := 37.87 ; +*---- +* Fill multicompo +*---- +BURN AUTO_LIB := EVO: AUTO_LIB FLUX1 TRACK_F1 :: + EDIT 0 + SAVE <<Timef>> DAY POWR <<Power>> +; +EDITION := EDI: FLUX1 AUTO_LIB TRACK_F1 GEOM_F1 :: + COND 0.625 MERGE CELL MICR RES SAVE ON 'fuel' + ; + +DATABASE := COMPO: DATABASE EDITION BURN :: +EDIT 3 +STEP UP * SET <<Timef>> DAY +; +* +WHILE Burni Burnc < DO +* + EVALUATE Burnf := Burni Delb + ; + EVALUATE Delt := Delb Power / ; + EVALUATE Timef := Timei Delt + ; +* + BURN AUTO_LIB := EVO: BURN AUTO_LIB FLUX1 TRACK_F1 :: + EDIT 1 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> + EXPM 1.0E15 GLOB ; +* + AUTO_LIB := USS: LIBRARY AUTO_LIB TRACK_SS :: + EDIT 0 GRMIN 45 PASS 2 + ; +*---- +* Flux Calculation +*---- + SYS := ASM: AUTO_LIB TRACK_F1 ; + FLUX1 := FLU: FLUX1 AUTO_LIB TRACK_F1 SYS :: TYPE B B1 SIGS ; + SYS := DELETE: SYS ; +*---- +* Fill multicompo +*---- + BURN AUTO_LIB := EVO: BURN AUTO_LIB FLUX1 TRACK_F1 :: + EDIT 0 + SAVE <<Timef>> DAY POWR <<Power>> + ; + EDITION := EDI: EDITION FLUX1 AUTO_LIB TRACK_F1 GEOM_F1 :: + COND 0.625 MERGE CELL SAVE ON 'fuel' + ; + + DATABASE := COMPO: DATABASE EDITION BURN :: + EDIT 3 + STEP UP * SET <<Timef>> DAY + ; +* + IF Burnf Burnc = THEN + IF Burnc 1000.0 = THEN + EVALUATE Delb Burnc := 1000.0 2000.0 ; + GREP: FLUX1 :: GETVAL 'B2 B1HOM' 1 1 1 >>b2<< ; + ECHO "buckling at 1000 MWd/tonne=" b2 ; + ENDIF ; + IF Burnc 500.0 = THEN + EVALUATE Delb Burnc := 500.0 1000.0 ; + ENDIF ; + IF Burnc 150.0 = THEN + EVALUATE Delb Burnc := 350.0 500.0 ; + ENDIF ; + ENDIF ; + EVALUATE Burni := Burnf ; + EVALUATE Timei := Timef ; +* +ENDWHILE ; +assertS FLUX1 :: 'K-INFINITY' 1 1.161285 ; +* +LIBRARY GEOM_SS TRACK_SS AUTO_LIB BURN := DELETE: +LIBRARY GEOM_SS TRACK_SS AUTO_LIB BURN ; +GEOM_F1 TRACK_F1 FLUX1 := DELETE: GEOM_F1 TRACK_F1 FLUX1 ; +* +******************************************** +* SPH equivalences +******************************************** +* +SPHGEOM := DATABASE :: STEP UP fuel STEP UP 'GEOMETRIES' STEP AT 1 ; +MTRACK := BIVACT: SPHGEOM :: DUAL (*IELEM=*) 1 (*ICOL=*) 2 ; +* +INTEGER istep := 0 ; +WHILE istep 6 < DO + EVALUATE istep := istep 1 + ; + DATABASE := SPH: DATABASE MTRACK :: + EDIT 1 + STEP UP fuel STEP AT <<istep>> + ; +ENDWHILE ; +* +******************************************** +* Verification calculation at 1000 MWd/tonne +******************************************** +* +MACROLIB := SPH: DATABASE :: + EDIT 5 + STEP UP fuel STEP AT 5 + MACRO LEAK <<b2>> OFF + ; +GEOM_F1 := DATABASE :: STEP UP fuel STEP UP 'GEOMETRIES' STEP AT 1 ; +TRACK_F1 := BIVACT: GEOM_F1 :: + TITLE 'VERIFICATION CALCULATION' + EDIT 2 MAXR 50 DUAL (*IELEM=*) 1 (*ICOL=*) 2 ; +SYS := ASM: MACROLIB TRACK_F1 :: ARM ; +FLUX1 := FLU: MACROLIB TRACK_F1 SYS :: TYPE K EXTE 100 5.0E-8 ; +* +* Non regression assert +* +EVALUATE REFVALUE := 1.0 ; +* +GREP: FLUX1 :: GETVAL K-EFFECTIVE 1 >>VALUE<< ; +EVALUATE DELTA := VALUE REFVALUE - REFVALUE / ABS ; +IF DELTA 5.0E-4 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "DELTA=" DELTA ; + PRINT "REFERENCE=" REFVALUE " CALCULATED=" VALUE ; + ABORT: ; +ENDIF ; +GEOM_F1 TRACK_F1 SYS FLUX1 := DELETE: GEOM_F1 TRACK_F1 SYS FLUX1 ; +* +ECHO "test TCDR36 completed" ; +END: ; +QUIT . diff --git a/Dragon/data/tdraglib_proc/TCDR37.c2m b/Dragon/data/tdraglib_proc/TCDR37.c2m new file mode 100644 index 0000000..8f0422a --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR37.c2m @@ -0,0 +1,164 @@ +*-------------------------------------------------------------------- +* TEST CASE TCDR37 +* 172-GROUP JEF2.2 DRAGLIB +* DB-VHTR Fuel Block Lattice calculation with NXT and BIHET +*-------------------------------------------------------------------- +* +*-------------------------------------------------------------------- +* Define STRUCTURES and MODULES used +*-------------------------------------------------------------------- +LINKED_LIST LIBRARY LIBRARY2 GEOM TRACK SYS FLUX + ausgabe EDIT DATABASE sphgeom mtrack ; +SEQ_BINARY TRACK_FIL mtrack_fil ; +SEQ_ASCII Fig_mix.ps Fig_reg.ps ; +MODULE GEO: LIB: ASM: FLU: SHI: NXT: END: PSP: ; +PROCEDURE assertS ; + +*-------------------------------------------------------------------- +* CREATION OF THE LIBRARY +*-------------------------------------------------------------------- +LIBRARY := LIB: :: + EDIT 0 + NMIX 8 CTRA APOL ANIS 2 + MIXS LIB: DRAGON FIL: DLIB_J2 + +*---- +* Mixtures for fuel cell F0 +*---- +* Graphite matrix + MIX 1 293.6 + GRAB = 'C0_GR' 8.524E-02 +* Fuel Kernel + MIX 2 293.6 + O16 = 'O16' 4.436E-02 + NP27 = 'Np237' 1.522E-03 1 + PU38 = 'Pu238' 6.464E-04 + PU39 = 'Pu239' 1.099E-02 1 + PU40 = 'Pu240' 5.084E-03 1 + PU41 = 'Pu241' 1.937E-03 1 + PU42 = 'Pu242' 1.074E-03 1 +* Porous C buffer + MIX 3 293.6 + GRAB = 'C0_GR' 5.265E-02 +* IPyC + MIX 4 293.6 + GRAB = 'C0_GR' 9.526E-02 +* SiC coating layer + MIX 5 293.6 + SI0 = 'Si0' 4.402E-02 + C12 = 'C0' 4.772E-02 +* OPyC + MIX 6 293.6 + GRAB = 'C0_GR' 9.526E-02 +* Helium Coolant surrounding fuel + MIX 7 293.6 + HE4 = 'He4' 7.0000E-04 +* Block Graphite (full density surrounding fuel) + MIX 8 293.6 + GRAB = 'C0_GR' 8.774E-02 +; + +*-------------------------------------------------------- +* GEOMETRY +*-------------------------------------------------------- +* ---- +* F = Fuel cell +* C = coolant hole cell (large) +* K = coolant hole cell (small) +* G = Full density graphite cell +* J = Reduced density graphite cell +* ---- +*orig+ +GEOM := GEO: :: HEX 7 + EDIT 1 + HBC COMPLETE REFL + CELL F0 F1 F1 F1 F1 F1 F1 + ::: F0 := GEO: HEXT 1 + SIDE 1.08542 + MIX 8 8 8 8 8 8 + CLUSTER FF0 + ::: FF0 := GEO: TUBE 2 NPIN 1 RPIN 0.0 APIN 0.0 + RADIUS 0.0 0.6 0.63 + MIX 9 7 ; + ; + ::: F1 := GEO: HEXT 1 + SIDE 1.08542 + MIX 8 8 8 8 8 8 + CLUSTER FF1 + ::: FF1 := GEO: TUBE 2 NPIN 1 RPIN 0.0 APIN 0.0 + RADIUS 0.0 0.6 0.63 + MIX 8 7 ; + ; +*--------------------------------------------------------- +* SPHERICAL MICRO-STRUCTURES * +* NMISTR = NUMBER OF MICRO STRUCTURE TYPES IN REGION +* NMILG = NUMBER OF MICRO STRUCTURE REGIONS +* NS = ARRAY OF SUB REGIONS IN THE MICRO STRUCTURES; LEN=NMILG +* RS = RADIUS OF MICROS STRUCTURES; LEN=NS(I); I=1,NMISTR +* milie = COMPOSITION OF EACH shell OF MICRO STRUCTURE; LEN=NMISTR; +* !!!!!!!! milie #'s ARE > THAN MIX #'S !!!!!!!! +* mixdil = BASE COMPOSITION OF EACH REGION; LEN=NMILG +* fract = MICROSTRUCTURE TYPE VOLUME FRACTION IN REGION LEN=NMILG +* mixgr = LIBRARY MIXTURES FOR EACH MICROSTRUCTURE SHELL; LEN=NS(I) +* NMISTR, NMILG + BIHET SPHE 1 1 + (* NS *) 5 + (* RS *) 0.0 0.0100 0.0220 0.0260 0.0295 0.0335 + (* milie *) 9 + (* mixdil *) 1 + (* loop over NMISTR *) + (* fract *) 0.02 + (* mixgr *) 2 3 4 5 6 +; + +*-------------------------------------------------------------------- +* Use NXT, SHI, ASM and FLU for flux solution +*-------------------------------------------------------------------- + +*** +* Tracking NXT +*** +TRACK TRACK_FIL := NXT: GEOM :: + EDIT 1 + MAXR 50000 + SLSI 0.65 QUAB -100 + TRAK TISO 4 10.0 +; + +*--- +* Self Shielding calculation +*--- +LIBRARY := SHI: LIBRARY TRACK TRACK_FIL :: + EDIT 0 +; + +*--- +* Multigroup Flux Calculation +*--- + SYS := ASM: LIBRARY TRACK TRACK_FIL :: + EDIT 0 PNOR HELI ; + + FLUX := FLU: SYS LIBRARY TRACK TRACK_FIL :: + EDIT 1 + TYPE K P1 EXTE 300 1E-5 THER 10 1E-5 ; + assertS FLUX :: K-EFFECTIVE 1 1.083298 ; + +*-------------------------------------------------------- +* Graphik +*-------------------------------------------------------- +* +Fig_reg.ps := PSP: TRACK :: + EDIT 5 + FILL CMYK NOCONTOUR + TYPE REGION +; + +Fig_mix.ps := PSP: TRACK :: + EDIT 5 + FILL CMYK NOCONTOUR + TYPE MIXTURE +; + +ECHO "test TCDR37 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR38.c2m b/Dragon/data/tdraglib_proc/TCDR38.c2m new file mode 100644 index 0000000..bf4eeed --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR38.c2m @@ -0,0 +1,136 @@ +*---- +* +* TEST CASE TCDR38 +* DISTRIBUTED SELF-SHIELDING +* 172-GROUP JEF2.2 DRAGLIB +* CARTESIAN CANDU6 CLUSTER WITH NXT +* +* Author: G. Marleau +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY LIBRARY2 CANDU6S CANDU6F TRACKS TRACKF SYS FLUX ; +SEQ_BINARY + INTLINS INTLINF ; +MODULE + GEO: NXT: LIB: USS: ASM: FLU: PSP: GREP: ABORT: END: ; +SEQ_ASCII Geos.ps Geof.ps :: FILE './Geos.ps' './Geof.ps' ; +REAL Kinf REFVALUE DELTA ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + EDIT 0 + NMIX 22 CTRA WIMS + SUBG +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 560.66 0.81212 O16 = O16 7.99449E-1 + D2D2O = H2_D2O 1.99768E-1 H1H2O = H1_H2O 7.83774E-4 + MIX 2 560.66 6.57 Nb93 = Nb93 2.50000 + B10 = B10 4.15800E-5 B11 = B11 1.68420E-4 + Zr90 = Zr90 5.01150E+1 Zr91 = Zr91 1.09200E+1 + Zr92 = Zr92 1.66725E+1 Zr94 = Zr94 1.70625E+1 + MIX 3 345.66 0.0014 He4 = He4 1.00000E+2 + MIX 4 345.66 6.44 Fe56 = Fe56 1.60000E-1 + Ni58 = Ni58 6.00000E-2 Cr52 = Cr52 1.10000E-1 + B10 = B10 6.13800E-5 B11 = B11 2.48620E-4 + Zr90 = Zr90 5.12509E+1 Zr91 = Zr91 1.11675E+1 + Zr92 = Zr92 1.70504E+1 Zr94 = Zr94 1.74493E+1 + MIX 5 345.66 1.082885 O16 = O16 7.98895E-1 + D2D2O = H2_D2O 2.01016E-1 H1H2O = H1_H2O 8.96000E-5 + MIX 6 941.29 10.4375010 O16 = O16 1.18473E+1 + U235 = U235 6.27118E-1 1 IRSET 0.0 81 + U238 = U238 8.75256E+1 1 IRSET 0.0 81 + MIX 7 COMB 6 1.0 MIX 8 COMB 6 1.0 MIX 9 COMB 6 1.0 + MIX 10 COMB 6 1.0 MIX 11 COMB 6 1.0 MIX 12 COMB 6 1.0 + MIX 13 COMB 6 1.0 MIX 14 COMB 6 1.0 MIX 15 COMB 6 1.0 + MIX 16 COMB 6 1.0 MIX 17 COMB 6 1.0 MIX 18 COMB 6 1.0 + MIX 19 COMB 6 1.0 MIX 20 COMB 6 1.0 MIX 21 COMB 6 1.0 + MIX 22 560.66 6.44 Fe56 = Fe56 1.60000E-1 + Ni58 = Ni58 6.00000E-2 Cr52 = Cr52 1.10000E-1 + B10 = B10 6.13800E-5 B11 = B11 2.48620E-4 + Zr90 = Zr90 5.12509E+1 Zr91 = Zr91 1.11675E+1 + Zr92 = Zr92 1.70504E+1 Zr94 = Zr94 1.74493E+1 + ; +*---- +* Geometry CANDU6S +*---- +REAL Rcomb4 := 0.6122 ; (*fuel radius*) +REAL Rcomb1 := 0.5 SQRT Rcomb4 * ; +REAL Rcomb2 := 0.8 SQRT Rcomb4 * ; +REAL Rcomb3 := 0.95 SQRT Rcomb4 * ; +REAL Rcool1 := 1.4885 2.8755 + 2. / ; +REAL Rcool2 := 2.8755 4.3305 + 2. / ; +CANDU6S := GEO: :: CARCEL 6 + X+ ALBE 1.0 X- ALBE 1.0 MESHX -14.2875 14.2875 + Y+ ALBE 1.0 Y- ALBE 1.0 MESHY -14.2875 14.2875 + RADIUS 0.0 <<Rcool1>> <<Rcool2>> 5.1689 5.6032 6.4478 6.5875 + MIX 1 1 1 2 3 4 5 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 5 MIX 6 7 8 9 22 NPIN 1 RPIN 0.0 APIN 0.0 + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> 0.654 ; + ::: ROD2 := GEO: ROD1 MIX 10 11 12 13 22 NPIN 6 RPIN 1.4885 APIN 0.0 ; + ::: ROD3 := GEO: ROD1 MIX 14 15 16 17 22 NPIN 12 RPIN 2.8755 + APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 18 19 20 21 22 NPIN 18 RPIN 4.3305 APIN 0.0 ; +; +TRACKS INTLINS := NXT: CANDU6S :: + EDIT 0 TISO 5 10.0 ; +Geos.ps := PSP: TRACKS ; +*---- +* Geometry CANDU6F +*---- +REAL xyp := 14.0 2.0 SQRT / ; +REAL xym := xyp -1. * ; +CANDU6F := GEO: :: CARCEL 7 3 3 + X+ REFL X- REFL MESHX -14.2875 <<xym>> <<xyp>> 14.2875 + Y+ REFL Y- REFL MESHY -14.2875 <<xym>> <<xyp>> 14.2875 + RADIUS 0.0 <<Rcool1>> <<Rcool2>> 5.16890 5.60320 6.44780 + 6.58750 14.0 + MIX 1 1 1 2 3 4 5 5 1 1 1 2 3 4 5 5 1 1 1 2 3 4 5 5 + 1 1 1 2 3 4 5 5 1 1 1 2 3 4 5 5 1 1 1 2 3 4 5 5 + 1 1 1 2 3 4 5 5 1 1 1 2 3 4 5 5 1 1 1 2 3 4 5 5 + SPLITX 5 1 5 + SPLITY 5 1 5 + SPLITR 3 2 2 1 1 1 10 + CLUSTER ROD1 ROD2 ROD3 ROD4 + ::: ROD1 := GEO: TUBE 5 MIX 6 7 8 9 22 NPIN 1 RPIN 0.0 APIN 0.0 + RADIUS 0.0 <<Rcomb1>> <<Rcomb2>> <<Rcomb3>> <<Rcomb4>> 0.654 ; + ::: ROD2 := GEO: ROD1 MIX 10 11 12 13 22 NPIN 6 RPIN 1.4885 APIN 0.0 ; + ::: ROD3 := GEO: ROD1 MIX 14 15 16 17 22 NPIN 12 RPIN 2.8755 + APIN 0.261799 ; + ::: ROD4 := GEO: ROD1 MIX 18 19 20 21 22 NPIN 18 RPIN 4.3305 APIN 0.0 ; +; +TRACKF INTLINF := NXT: CANDU6F :: + EDIT 0 TISO 5 10.0 ; +Geof.ps := PSP: TRACKF ; +*---- +* Self-Shielding calculation NXT +* Transport calculation NXT +* Flux calculation for critical buckling +*---- +LIBRARY2 := USS: LIBRARY TRACKS INTLINS :: EDIT 1 PASS 2 ; +SYS := ASM: LIBRARY2 TRACKF INTLINF :: EDIT 0 ; +FLUX := FLU: SYS LIBRARY2 TRACKF :: + EDIT 2 TYPE B B1 SIGS ; + +GREP: FLUX :: GETVAL 'K-INFINITY' 1 >>Kinf<< ; +ECHO "Kinf=" Kinf ; +EVALUATE REFVALUE := 1.110458 ; +EVALUATE DELTA := Kinf REFVALUE - REFVALUE / ABS ; +IF DELTA 1.0E-3 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFVALUE " CALCULATED=" Kinf ; + ABORT: ; +ENDIF ; + +ECHO "test TCDR38 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR41.c2m b/Dragon/data/tdraglib_proc/TCDR41.c2m new file mode 100644 index 0000000..7c4e5a1 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR41.c2m @@ -0,0 +1,169 @@ +*--- +* Regression test 2 for MCCG: and MOCC: +*--- +* TEST CASE TCDR41 +* ROWLANDS PINCELL BENCHMARK 1 +* DRAGLIB JEF2.2 XMAS 172 GROUPS LIBRARY FILE DLIB_J2 +*--- + +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY LIBRARY2 CELL TRACK SYS FLUX EDIT TEMP ; +MODULE GEO: EXCELT: MCCGT: LIB: USS: ASM: FLU: EDI: DELETE: END: ; +SEQ_BINARY TRACK_FIL ; +CHARACTER PolarAng := "LCMD" ; +REAL densH := 0.066988 ; +REAL densO := 0.033414 ; +REAL temp1 := 293.0 ; +REAL temp2 := 293.0 ; +REAL temp3 := 293.0 ; +REAL Buckl := 0.0 ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections +*---- +LIBRARY := LIB: :: +EDIT 1 NMIX 9 PTSL + CTRA APOL + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 <<temp1>> (*fuel*) + U238 = U238 0.022604 1 IRSET PTSL 1 + U235 = U235 7.0803E-4 1 IRSET PTSL 1 + O16 = O16 0.046624 NOEV + MIX 2 COMB 1 1.0 + MIX 3 COMB 1 1.0 + MIX 4 COMB 1 1.0 + MIX 5 COMB 1 1.0 + MIX 6 COMB 1 1.0 + MIX 7 <<temp2>> (*clad*) + ZR = Zr0 0.043241 2 IRSET PTSL NONE NOEV + MIX 8 <<temp3>> (*water1*) + H1 = H1_H2O <<densH>> NOEV + O16 = O16 <<densO>> NOEV + MIX 9 <<temp3>> (*water2*) + H1 = H1_H2O <<densH>> NOEV + O16 = O16 <<densO>> NOEV + ; +*---- +* Geometry CELL : a single UO2 pincell with 6 annuli in fuel +*---- +CELL := GEO: :: +CARCEL 8 EDIT 1 + X- REFL X+ REFL Y- REFL Y+ REFL + MESHX 0.0 1.2 + MESHY 0.0 1.2 + RADIUS 0.0 0.2529822 0.334664 0.3577709 0.3794733 0.3898718 0.4 + 0.45 0.5748331 + MIX 1 2 3 4 5 6 7 8 9 ; +*---- +* MCCG (isotropic BC) +*---- +*---- +* Tracking calculation +* Self-Shielding calculation +*---- +TRACK TRACK_FIL := EXCELT: CELL :: + EDIT 1 ANIS 1 + MAXR 40 + TITLE 'UO2 Rowlands benchmark 1 with MCCG' + TRAK TISO 12 12.0 NOSY ; +TEMP := TRACK ; +TRACK := MCCGT: TRACK TRACK_FIL :: + EDIT 1 <<PolarAng>> 2 + AAC 1 SCR 0 EPSI 1E-5 STIS 1 + MAXI 100 KRYL 30 HDD 0.0 ; +LIBRARY2 := USS: LIBRARY TRACK TRACK_FIL :: + EDIT 1 ARM + CALC REGI W1 U235 ALL + REGI W1 ZR ALL + REGI W1 U238 1 + REGI W2 U238 2 + REGI W3 U238 3 + REGI W4 U238 4 + REGI W5 U238 5 + REGI W6 U238 6 + ENDC ; +*--- +* Multigroup Flux Calculation +* Edition +*-- +TRACK := DELETE: TRACK ; +TRACK := TEMP ; +TRACK := MCCGT: TRACK TRACK_FIL :: + EDIT 1 <<PolarAng>> 2 + AAC 80 TMT SCR 0 EPSI 1E-5 STIS 1 + MAXI 1 KRYL 0 HDD 0.0 ; +SYS := ASM: LIBRARY2 TRACK TRACK_FIL :: + EDIT 2 ARM ; +FLUX := FLU: LIBRARY2 TRACK SYS TRACK_FIL :: + EDIT 1 TYPE K B0 SIGS BUCK <<Buckl>> ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.392635 ; +EDIT := EDI: FLUX LIBRARY2 TRACK CELL :: +EDIT 2 +COND 27 59 95 ; +*--- +* Memory Release +*--- +TRACK LIBRARY2 SYS FLUX TRACK_FIL EDIT TEMP := DELETE: +TRACK LIBRARY2 SYS FLUX TRACK_FIL EDIT TEMP ; +*---- +* MOCC (cyclic tracking) +*---- +*---- +* Tracking calculation +* Self-Shielding calculation +*---- +TRACK TRACK_FIL := EXCELT: CELL :: + EDIT 1 ANIS 1 + MAXR 40 + TITLE 'UO2 Rowlands benchmark 1 with MOCC' + TRAK TSPC 7 12.0 NOSY ; +TEMP := TRACK ; +TRACK := MCCGT: TRACK TRACK_FIL :: + EDIT 1 <<PolarAng>> 2 + AAC 1 TMT SCR 0 EPSI 1E-5 + MAXI 100 KRYL 30 HDD 100000.0 ; +LIBRARY2 := USS: LIBRARY TRACK TRACK_FIL :: + EDIT 1 ARM + CALC REGI W1 U235 ALL + REGI W1 ZR ALL + REGI W1 U238 1 + REGI W2 U238 2 + REGI W3 U238 3 + REGI W4 U238 4 + REGI W5 U238 5 + REGI W6 U238 6 + ENDC ; +*--- +* Multigroup Flux Calculation +* Edition +*-- +!TRACK TRACK_FIL := DELETE: TRACK TRACK_FIL ; +!TRACK TRACK_FIL := EXCELT: CELL :: +! EDIT 1 ANIS 1 +! MAXR 40 +! TITLE 'UO2 Rowlands benchmark 1 with MOCC' +! TRAK TSPC 7 12.0 NOSY ; +TRACK := DELETE: TRACK ; +TRACK := TEMP ; +TRACK := MCCGT: TRACK TRACK_FIL :: + EDIT 1 <<PolarAng>> 2 + AAC 80 TMT SCR 0 EPSI 1E-5 + MAXI 1 KRYL 0 HDD 100000.0 ; +SYS := ASM: LIBRARY2 TRACK TRACK_FIL :: + EDIT 2 ARM ; +FLUX := FLU: LIBRARY2 TRACK SYS TRACK_FIL :: + EDIT 1 TYPE K B0 SIGS BUCK <<Buckl>> ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.393138 ; +EDIT := EDI: FLUX LIBRARY2 TRACK CELL :: +EDIT 2 +COND 27 59 95 ; +*--- +* Memory Release +*--- +CELL TRACK LIBRARY LIBRARY2 SYS FLUX TRACK_FIL EDIT := DELETE: +CELL TRACK LIBRARY LIBRARY2 SYS FLUX TRACK_FIL EDIT ; +ECHO "test TCDR41 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR42.c2m b/Dragon/data/tdraglib_proc/TCDR42.c2m new file mode 100644 index 0000000..2023f36 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR42.c2m @@ -0,0 +1,129 @@ +*---- +* +* TEST CASE TCDR42 +* TEST SHI: WITH FLUX-CURRENT EURYDICE ITERATIONS +* DISTRIBUTED SELF-SHIELDING +* 172-GROUP JEF2.2 DRAGLIB +* +* Author: A. Hebert +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST GEOM TRACK LIBRARY LIBRARY2 CP FLUX ; +MODULE LIB: GEO: SYBILT: SHI: ASM: FLU: DELETE: END: ; +INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 GAIN1 MODE1 := + 1 2 3 4 5 6 7 ; (* DISTRIBUTED SELF-SHIELDING *) +REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 ; +PROCEDURE assertS ; +* +EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 := 0.2529822 0.334664 +0.3577709 0.3794733 0.40 ; +* +GEOM := GEO: :: CAR2D 5 5 (*ASSEMBLAGE 10 X 10*) + X- DIAG X+ REFL + Y- REFL Y+ DIAG + CELL C1 C3 C2 C3 C4 + C3 C3 C3 C4 + C2 C3 C4 + C3 C4 + C5 + MERGE 1 2 3 4 5 + 6 7 8 9 + 10 11 9 + 12 9 + 13 + ::: C1 := GEO: CARCEL 1 (*CELLULE SANS COMBUSTIBLE*) + MESHX 0.0 1.262082 + MESHY 0.0 1.262082 + RADIUS 0.0 4.12660E-01 + MIX <<MODE1>> <<MODE1>> + ; + ::: C2 := GEO: CARCEL 3 (*CELLULE EMPOISONNEE*) + MESHX 0.0 1.262082 + MESHY 0.0 1.262082 + RADIUS 0.0 3.25296E-01 4.60039E-01 5.63430E-01 + MIX <<GAIN1>> <<GAIN1>> <<GAIN1>> <<MODE1>> + ; + ::: C3 := GEO: CARCEL 7 (*CELLULE NORMALE*) + MESHX 0.0 1.262082 + MESHY 0.0 1.262082 + RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>> + <<RCOMB5>> 0.45 0.5748331 + MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>> + <<COMB0105>> <<GAIN1>> <<MODE1>> <<MODE1>> + ; + ::: C4 := GEO: C3 (*CELLULE PERIPHERIQUE*) + MESHX 0.0 1.322082 + ; + ::: C5 := GEO: C4 (*CELLULE DE COIN*) + MESHY 0.0 1.322082 + ; + ; +TRACK := SYBILT: GEOM :: + MAXR 100 MAXZ 100000 + TITLE 'TCDR42: MULTICELL 10 X 10 PWR BENCHMARK WITH POISON' + DP00 ; +*---- +* MAKE MICROLIB +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 + ADED 4 NELAS N4N N2N N3N +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX <<COMB0101>> 900.0 (*COMB0101*) + O16 = O16 4.6624E-2 + U235 = U235 7.0803E-4 1 IRSET 0.0 81 + U238 = U238 2.2604E-2 1 IRSET 0.0 81 + MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) + MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) + MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) + MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) + MIX <<GAIN1>> 600.0 (*GAIN1*) + Zr0 = Zr0 4.3241E-2 2 IRSET 0.0 81 + MIX <<MODE1>> 574.0 (*MODE1*) + H1 = H1_H2O 4.6892E-2 + O16 = O16 2.3446E-2 +; +*---- +* SHI LEVEL 0 WITH PIJ +*---- +LIBRARY2 := SHI: LIBRARY TRACK :: + EDIT 1 LJ GRMIN 45 LEVEL 0 MXIT 3 PIJ ; +CP := ASM: LIBRARY2 TRACK :: PIJ ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.342790 ; +LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ; +*---- +* SHI LEVEL 0 WITH ARM (FLUX-CURRENT EURYDICE ITERATIONS) +*---- +LIBRARY2 := SHI: LIBRARY TRACK :: + EDIT 1 LJ GRMIN 45 LEVEL 0 MXIT 3 ARM ; +CP := ASM: LIBRARY2 TRACK :: ARM ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.342790 ; +LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ; +*---- +* SHI LEVEL 1 WITH PIJ +*---- +LIBRARY2 := SHI: LIBRARY TRACK :: + EDIT 1 LJ GRMIN 45 LEVEL 1 MXIT 3 PIJ ; +CP := ASM: LIBRARY2 TRACK :: PIJ ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.345739 ; +LIBRARY2 CP FLUX := DELETE: LIBRARY2 CP FLUX ; +*---- +* SHI LEVEL 1 WITH ARM (FLUX-CURRENT EURYDICE ITERATIONS) +*---- +LIBRARY2 := SHI: LIBRARY TRACK :: + EDIT 1 LJ GRMIN 45 LEVEL 1 MXIT 3 ARM ; +CP := ASM: LIBRARY2 TRACK :: ARM ; +FLUX := FLU: CP LIBRARY2 TRACK :: TYPE K ; +assertS FLUX :: K-EFFECTIVE 1 1.345739 ; +* +ECHO "test TCDR42 completed" ; +END: ; diff --git a/Dragon/data/tdraglib_proc/TCDR43.c2m b/Dragon/data/tdraglib_proc/TCDR43.c2m new file mode 100644 index 0000000..01c97fd --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR43.c2m @@ -0,0 +1,277 @@ +*---- +* +* TEST CASE TCDR43 +* DRAGON DATA FILE FOR MULTI-CELL ICSBEP PROBLEMS. +* LEU-COMP-THERM-041-c01 BENCHMARK. +* DISTRIBUTED SELF-SHIELDING +* 172-GROUP JEF2.2 DRAGLIB +* +* Author: N. Martin +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST GEOMETRY MACRO CP1 TRACK1 FLUX1 + TRACK1_F TRACK1_SS GEOMETRY_F GEOMETRY_SS + ASSMB AUTO_LIB LIBRARY1 EDIT1 SPHGEOM MTRACK ; + +MODULE LIB: GEO: SYBILT: USS: ASM: FLU: EDI: M2T: DELETE: GREP: SPH: + END: ; +SEQ_ASCII APOTR ; +REAL kiso ; +PROCEDURE assertS ; + +GEOMETRY_F := GEO: :: CAR2D 8 8 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + + CELL C1 C1 C1 C1 C1 C1 C1 C3 + C2 C1 C1 C2 C1 C1 C3 + C1 C1 C1 C1 C1 C3 + C2 C1 C1 C1 C3 + C1 C1 C1 C3 + C1 C1 C3 + C1 C3 + C4 + + MERGE 1 2 3 4 5 6 7 8 + 9 10 11 12 13 14 15 + 16 17 18 19 20 21 + 22 23 24 25 26 + 27 28 29 30 + 31 32 33 + 34 35 + 36 + + +!Fuel Pin and clad + ::: C1 := GEO: CARCEL 7 + MESHX 0.0 1.26 + MESHY 0.0 1.26 + SPLITX 4 + SPLITY 4 + RADIUS 0.0 0.2529822 0.334664 0.3577709 + 0.3794733 0.3898718 0.4015 0.4715 + MIX 1 2 3 4 5 6 7 8 ; + +!Water hole + + ::: C2 := GEO: CARCEL 2 + MESHX 0.0 1.26 MESHY 0.0 1.26 + SPLITX 4 + SPLITY 4 + RADIUS 0.0 0.41 0.47 + MIX 9 10 11 ; + +!Water rectangular cell + ::: C3 := GEO: CARCEL 0 + MESHX 0.0 2.24 + MESHY 0.0 1.26 + SPLITX 8 + SPLITY 4 + MIX 12 ; + +!Water square cell + ::: C4 := GEO: CARCEL 0 + MESHX 0.0 2.24 + MESHY 0.0 2.24 + SPLITX 8 + SPLITY 8 + MIX 13 ; +; + +GEOMETRY_SS := GEO: :: CAR2D 8 8 + X- DIAG X+ REFL + Y- SYME Y+ DIAG + + CELL C1 C1 C1 C1 C1 C1 C1 C3 + C2 C1 C1 C2 C1 C1 C3 + C1 C1 C1 C1 C1 C3 + C2 C1 C1 C1 C3 + C1 C1 C1 C3 + C1 C1 C3 + C1 C3 + C4 + + MERGE 1 1 1 1 1 1 1 3 + 2 1 1 2 1 1 3 + 1 1 1 1 1 3 + 2 1 1 1 3 + 1 1 1 3 + 1 1 3 + 1 3 + 4 + +!Fuel Pin and clad + ::: C1 := GEO: CARCEL 7 + MESHX 0.0 1.26 + MESHY 0.0 1.26 + RADIUS 0.0 0.2529822 0.334664 0.3577709 + 0.3794733 0.3898718 0.4015 0.4715 + MIX 1 2 3 4 5 6 7 8 ; + +!Water hole + + ::: C2 := GEO: CARCEL 2 + MESHX 0.0 1.26 MESHY 0.0 1.26 + RADIUS 0.0 0.41 0.47 + MIX 9 10 11 ; + +!Water rectangular cell + ::: C3 := GEO: CARCEL 0 + MESHX 0.0 2.24 MESHY 0.0 1.26 + MIX 12 ; + +!Water square cell + ::: C4 := GEO: C3 + MESHY 0.0 2.24 + MIX 13 ; +; + +LIBRARY1 := LIB: :: + EDIT 1 + NMIX 13 + CTRA APOL + ANIS 2 + ADED 4 NELAS N4N N2N N3N + SUBG +MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 293.0 + U234 = U234 4.55680E-06 1 IRSET 0.0 81 + U235 = U235 6.92410E-04 1 IRSET 0.0 81 + U236 = U236 5.46810E-06 1 IRSET 0.0 81 + U238 = U238 2.20820E-02 1 IRSET 0.0 81 + O16 = O16 4.55680E-02 + MIX 2 COMB 1 1.0 + MIX 3 COMB 1 1.0 + MIX 4 COMB 1 1.0 + MIX 5 COMB 1 1.0 + MIX 6 COMB 1 1.0 + + MIX 7 293.0 + Fe54 = Fe54 7.84700E-06 + Fe56 = Fe56 1.21988E-04 + Fe57 = Fe57 2.79300E-06 + Fe58 = Fe58 3.72400E-07 + Cr50 = Cr50 3.41082E-06 + Cr52 = Cr52 6.57752E-05 + Cr53 = Cr53 7.45750E-06 + Cr54 = Cr54 1.85652E-06 + Zr0 = Zr0 3.67800E-02 2 IRSET 0.0 81 + Sn115 = Sn115 1.54440E-06 + Sn116 = Sn116 6.23337E-05 + Sn117 = Sn117 3.29472E-05 + Sn118 = Sn118 1.03904E-04 + Sn119 = Sn119 3.68082E-05 + Sn120 = Sn120 1.39811E-04 + Sn122 = Sn122 1.98627E-05 + Sn124 = Sn124 2.48391E-05 + + MIX 8 293.0 + H1_H2O = H1_H2O 6.67340E-02 + O16 = O16 3.33670E-02 + MIX 9 293.0 + H1_H2O = H1_H2O 6.67340E-02 + O16 = O16 3.33670E-02 + MIX 11 293.0 + H1_H2O = H1_H2O 6.67340E-02 + O16 = O16 3.33670E-02 + + MIX 12 293.0 + H1_H2O = H1_H2O 6.67340E-02 + O16 = O16 3.33670E-02 + + MIX 13 293.0 + H1_H2O = H1_H2O 6.67340E-02 + O16 = O16 3.33670E-02 + + MIX 10 293.0 + Fe54 = Fe54 8.85000E-06 + Fe56 = Fe56 1.37580E-04 + Fe57 = Fe57 3.15000E-06 + Fe58 = Fe58 4.20000E-07 + Cr50 = Cr50 3.85836E-06 + Cr52 = Cr52 7.44055E-05 + Cr53 = Cr53 8.43600E-06 + Cr54 = Cr54 2.10012E-06 + Zr0 = Zr0 4.15800E-02 3 IRSET 0.0 81 + Sn115 = Sn115 1.74600E-06 + Sn116 = Sn116 7.04705E-05 + Sn117 = Sn117 3.72480E-05 + Sn118 = Sn118 1.17467E-04 + Sn119 = Sn119 4.16130E-05 + Sn120 = Sn120 1.58062E-04 + Sn122 = Sn122 2.24555E-05 + Sn124 = Sn124 2.80815E-05 + ; + +TRACK1_SS := SYBILT: GEOMETRY_SS :: MAXR 50000 MAXZ 50000 +DP01 ; + +TRACK1_F := SYBILT: GEOMETRY_F :: MAXR 50000 MAXZ 50000 +DP01 ; + +AUTO_LIB := USS: LIBRARY1 TRACK1_SS :: + CALC REGI W1 U235 ALL + REGI W1 U234 ALL + REGI W1 U236 ALL + REGI W1 Zr0 ALL + + REGI W1 U238 1 + REGI W2 U238 2 + REGI W3 U238 3 + REGI W4 U238 4 + REGI W5 U238 5 + REGI W6 U238 6 + ENDC + + ; +* +CP1 := ASM: AUTO_LIB TRACK1_F :: PIJ ; +FLUX1 := FLU: CP1 AUTO_LIB TRACK1_F :: + TYPE K ; + +GREP: FLUX1 :: GETVAL 'K-INFINITY' 1 1 1 >>kiso<< ; +ECHO "K-inf is" kiso ; + +SPHGEOM := GEO: :: CAR2D 8 8 (*HOMOGENIZED ASSEMBLY 15 X 15*) + X- DIAG X+ REFL + Y- SYME Y+ DIAG + + MIX 1 1 1 1 1 1 1 2 + 1 1 1 1 1 1 2 + 1 1 1 1 1 2 + 1 1 1 1 2 + 1 1 1 2 + 1 1 2 + 1 2 + 2 + + MESHX 0.0 1.26 2.52 3.78 5.04 6.3 7.56 8.82 11.06 + MESHY 0.0 1.26 2.52 3.78 5.04 6.3 7.56 8.82 11.06 + ; + +EDIT1 := EDI: FLUX1 AUTO_LIB TRACK1_F SPHGEOM :: + EDIT 1 + MERG MIX 1 1 1 1 1 1 1 1 1 1 1 2 2 + SAVE ON 'EDITCDAT 1' + ; +MTRACK := SYBILT: SPHGEOM :: MAXJ 500 QUA2 20 3 DP01 ; +EDIT1 := SPH: EDIT1 MTRACK :: STD ARM ; +MTRACK := DELETE: MTRACK ; + +MACRO := EDIT1 :: STEP UP 'EDITCDAT 1' STEP UP 'MACROLIB' ; +APOTR := M2T: MACRO :: NOMA PN 2 TRAN +MIX Fuel FROM 1 ENDMIX +MIX Coolant FROM 2 ENDMIX +; +* +* VERIFICATION CALCULATION +TRACK1_F CP1 FLUX1 := DELETE: TRACK1_F CP1 FLUX1 ; +TRACK1_F := SYBILT: SPHGEOM :: MAXJ 500 QUA2 20 3 DP01 ; +CP1 := ASM: MACRO TRACK1_F :: PIJ ; +FLUX1 := FLU: MACRO TRACK1_F CP1 :: TYPE K ; + +assertS FLUX1 :: K-EFFECTIVE 1 1.264118 ; +ECHO "test TCDR43 completed" ; +END: ; diff --git a/Dragon/data/tdraglib_proc/TCDR44.c2m b/Dragon/data/tdraglib_proc/TCDR44.c2m new file mode 100644 index 0000000..2da37e5 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR44.c2m @@ -0,0 +1,98 @@ +*---- +* +* TEST CASE TCDR44 WITH MONTE CARLO MODULE +* FIRST UO2 ROWLAND'S BENCHMARK +* DISTRIBUTED SELF-SHIELDING +* 172-GROUP JEF2.2 DRAGLIB +* +* Author: A. Hebert +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST GEOM TRACKSS TRACK LIBRARY LIBRARY2 FLUX ; +MODULE LIB: GEO: NXT: USS: FLU: GREP: MC: ABORT: END: ; +SEQ_BINARY TRKFIL ; +INTEGER COMB0101 COMB0102 COMB0103 COMB0104 COMB0105 COMB0106 GAIN1 +MODE1 MODE2 := +1 2 3 4 5 6 7 8 9 ; (* DISTRIBUTED SELF-SHIELDING *) +REAL RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 ; +REAL VALUE REFVALUE DELTA ; +* +EVALUATE RCOMB1 RCOMB2 RCOMB3 RCOMB4 RCOMB5 RCOMB6 := 0.2529822 +0.334664 0.3577709 0.3794733 0.3898718 0.40 ; + +GEOM := GEO: :: CAR2D 1 1 +X- REFL X+ REFL +Y- REFL Y+ REFL + CELL C1 + ::: C1 := GEO: CARCEL 8 + MESHX 0.0 1.20 + MESHY 0.0 1.20 + RADIUS 0.0 <<RCOMB1>> <<RCOMB2>> <<RCOMB3>> <<RCOMB4>> + <<RCOMB5>> <<RCOMB6>> 0.45 0.5 + MIX <<COMB0101>> <<COMB0102>> <<COMB0103>> <<COMB0104>> + <<COMB0105>> <<COMB0106>> <<GAIN1>> <<MODE1>> <<MODE2>> +; +; + +TRACKSS TRKFIL := NXT: GEOM :: + TITLE 'FIRST UO2 ROWLAND S BENCHMARK (USS)' + EDIT 0 TISO 20 5.0 + ; + +TRACK := NXT: GEOM :: + TITLE 'FIRST UO2 ROWLAND S BENCHMARK (MC)' + EDIT 0 MC ; + +LIBRARY := LIB: :: +EDIT 0 +NMIX 9 (*MAXIMUM OF MATERIAL MIXTURES*) +ANIS 2 +ADED 4 NELAS N4N N2N N3N +CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) +SUBG (*HELIOS TYPE PROBABILITY TABLES*) +* +MIXS LIB: DRAGON FIL: DLIB_J2 +MIX <<COMB0101>> 293.0 (*COMB0101*) + O16 = O16 4.6624E-2 + U235 = U235 7.0803E-4 1 IRSET 0.0 81 + U238 = U238 2.2604E-2 1 IRSET 0.0 81 +MIX <<COMB0102>> COMB <<COMB0101>> 1.0 (*COMB0102*) +MIX <<COMB0103>> COMB <<COMB0101>> 1.0 (*COMB0103*) +MIX <<COMB0104>> COMB <<COMB0101>> 1.0 (*COMB0104*) +MIX <<COMB0105>> COMB <<COMB0101>> 1.0 (*COMB0105*) +MIX <<COMB0106>> COMB <<COMB0101>> 1.0 (*COMB0106*) +MIX <<GAIN1>> 293.0 (*GAIN1*) + Zr0 = Zr0 4.3241E-2 2 IRSET 0.0 81 +MIX <<MODE1>> 293.0 (*MODE1*) + H1 = H1_H2O 6.6988E-2 + O16 = O16 3.3494E-2 +MIX <<MODE2>> 293.0 (*MODE2*) + H1 = H1_H2O 6.6988E-2 + O16 = O16 3.3494E-2 +; + +LIBRARY2 := USS: LIBRARY TRACKSS TRKFIL :: + EDIT 0 TRAN PASS 2 GRMIN 45 MAXST 50 ; + +FLUX := MC: TRACK LIBRARY2 :: + EDIT 2 KCODE 1000 25 50 SEED 23412 ; +* +* Non regression assert +* +EVALUATE REFVALUE := 1.395400 ; +* +GREP: FLUX :: GETVAL K-EFFECTIVE 1 >>VALUE<< ; +EVALUATE DELTA := VALUE REFVALUE - REFVALUE / ABS ; +IF DELTA 1.0E-2 < THEN + PRINT "TEST SUCCESSFUL; DELTA=" DELTA ; +ELSE + PRINT "------------" ; + PRINT "TEST FAILURE" ; + PRINT "------------" ; + PRINT "REFERENCE=" REFVALUE " CALCULATED=" VALUE ; + ABORT: ; +ENDIF ; +ECHO "test TCDR44 completed" ; +END: ; diff --git a/Dragon/data/tdraglib_proc/TCDR45.c2m b/Dragon/data/tdraglib_proc/TCDR45.c2m new file mode 100644 index 0000000..46710e8 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR45.c2m @@ -0,0 +1,62 @@ +*---- +* TEST CASE TRDR45 +* MIXTURES DEFINED WITH COMB INSTRUCTION +* 172-GROUP JEF2.2 DRAGLIB +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY ; +MODULE LIB: END: ; +PROCEDURE assertS ; +*---- +* Concentrations given to LIB: +*---- +REAL Ni58 := 3.9171E-3 ; +REAL O16 := 2.3954E-2 ; +REAL H1_H2O := O16 2.0 * ; +REAL B10 := 4.7892E-6 ; +REAL B11 := 1.9157E-5 ; + +*---- +* Using COMB instruction to define MIX 3 +*---- +LIBRARY := LIB: :: + NMIX 3 + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 + Ni58 = Ni58 <<Ni58>> + MIX 2 579.9 + O16 = O16 <<O16>> + H1H2O = H1_H2O <<H1_H2O>> + B10 = B10 <<B10>> + B11 = B11 <<B11>> + MIX 3 + COMB 1 0.01 + 2 0.99 +; +*---- +* Out-of-LIB: dilution +*---- +REAL Ni58dil := Ni58 0.01 * ; +REAL O16dil := O16 0.99 * ; +REAL H1_H2Odil := H1_H2O 0.99 * ; +REAL B10dil := B10 0.99 * ; +REAL B11dil := B11 0.99 * ; +* +assertS LIBRARY :: 'ISOTOPESDENS' 1 <<Ni58>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 2 <<O16>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 3 <<H1_H2O>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 4 <<B10>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 5 <<B11>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 6 <<Ni58dil>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 7 <<O16dil>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 8 <<H1_H2Odil>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 9 <<B10dil>> ; +assertS LIBRARY :: 'ISOTOPESDENS' 10 <<B11dil>> ; +* +ECHO "test TCDR45 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR46.c2m b/Dragon/data/tdraglib_proc/TCDR46.c2m new file mode 100644 index 0000000..9d36bdc --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR46.c2m @@ -0,0 +1,111 @@ +*---- +* TEST CASE TRDR46 +* 15 X 15 PWR ASSEMBLY WITHOUT POISON +* 172-GROUP JEF2.2 DRAGLIB +* SPH EQUIVALENCE WITH SELENGUT NORMALIZATION +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY LIBRARY_SS CP FLUX DISCR DISCR_SS EDITION + CELLS EIGHTH ASSMB ASSMB_SS MACDIR ; +MODULE + LIB: GEO: SYBILT: USS: ASM: FLU: EDI: SPH: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 103 CTRA APOL SUBG + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 10 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + MIX 20 579.9 + Zr90 = Zr0 4.18621E-2 2 + MIX 100 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + MIX 101 COMB 100 1.0 + MIX 102 COMB 100 1.0 + MIX 103 COMB 100 1.0 + ; +*---- +* Geometry : a 15 X 15 normal PWR assembly +* contains T : guide thimble tube +* C1 : normal fuel cell +* C2 : peripheral cell +* C3 : corner cell +*---- +REAL Rfuel := 0.4646 ; +REAL Rfuel1 := Rfuel .50 SQRT * ; +REAL Rfuel2 := Rfuel .80 SQRT * ; +REAL Rfuel3 := Rfuel .95 SQRT * ; +CELLS := GEO: :: + ::: T := GEO: CARCEL 2 + MESHX 0. 1.43 + MESHY 0. 1.43 + RADIUS 0. 0.55 0.57 + MIX 10 20 10 ; + ::: C1 := GEO: CARCEL 5 + MESHX 0. 1.43 + MESHY 0. 1.43 + RADIUS 0. <<Rfuel1>> <<Rfuel2>> <<Rfuel3>> <<Rfuel>> 0.53594 + MIX 100 101 102 103 20 10 ; + ::: C2 := GEO: C1 MESHX 0. 1.48 MIX 100 101 102 103 20 10 ; + ::: C3 := GEO: C2 MESHY 0. 1.48 MIX 100 101 102 103 20 10 ; +; +EIGHTH := GEO: CELLS :: CAR2D 8 8 +X- DIAG X+ REFL +Y- SYME Y+ DIAG ; +ASSMB := GEO: EIGHTH :: +CELL T C1 C1 C1 T C1 C1 C2 + C1 C1 C1 C1 C1 C1 C2 + C1 C1 C1 T C1 C2 + T C1 C1 C1 C2 + C1 C1 C1 C2 + T C1 C2 + C1 C2 + C3 ; +ASSMB_SS := GEO: ASSMB :: +MERGE 1 2 2 2 1 2 2 4 + 2 2 2 2 2 2 4 + 2 2 2 1 2 4 + 1 2 2 2 4 + 2 2 2 4 + 1 2 4 + 2 4 + 3 ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for B1 homogeneous leakage +*---- +DISCR := SYBILT: ASSMB :: MAXZ 1000000 MAXR 1000000 + QUA2 3 3 GAUSS DP01 ; +DISCR_SS := SYBILT: ASSMB_SS :: MAXZ 1000000 MAXR 1000000 + QUA2 3 3 GAUSS DP01 ; +LIBRARY_SS := USS: LIBRARY DISCR_SS ; +CP := ASM: LIBRARY_SS DISCR :: ARM ; +FLUX := FLU: CP LIBRARY_SS DISCR :: EDIT 2 TYPE B B1 SIGS ; +assertS FLUX :: 'K-INFINITY' 1 1.279648 ; +EDITION := EDI: LIBRARY_SS DISCR FLUX ASSMB :: + EDIT 4 SAVE MICR RES MERGE COMP COND 0.625 JOUT ; +EDITION := SPH: EDITION :: EDIT 9 ALBS SELE_ALB ; +* +MACDIR := EDITION :: STEP UP 'REF-CASE0001' STEP UP MACROLIB + STEP UP GROUP STEP AT 1 ; +assertS MACDIR :: 'NTOT0' 1 5.305629E-01 ; +assertS MACDIR :: 'NUSIGF' 1 6.686473E-03 ; +assertS MACDIR :: 'DIFF' 1 1.400524E+00 ; +assertS MACDIR :: 'SIGW00' 1 5.050910E-01 ; +MACDIR := DELETE: MACDIR ; +MACDIR := EDITION :: STEP UP 'REF-CASE0001' STEP UP MACROLIB + STEP UP GROUP STEP AT 2 ; +assertS MACDIR :: 'NTOT0' 1 1.359428E+00 ; +assertS MACDIR :: 'NUSIGF' 1 1.365678E-01 ; +assertS MACDIR :: 'DIFF' 1 4.055927E-01 ; +assertS MACDIR :: 'SIGW00' 1 1.271332E+00 ; +MACDIR := DELETE: MACDIR ; +ECHO "test TCDR46 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR47.c2m b/Dragon/data/tdraglib_proc/TCDR47.c2m new file mode 100644 index 0000000..1c9f41d --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR47.c2m @@ -0,0 +1,147 @@ +*---- +* TEST CASE TRDR47 +* 15 X 15 PWR ASSEMBLY WITHOUT POISON +* SPH EQUIVALENCE WITH SELENGUT NORMALIZATION +* 172-GROUP JEF2.2 DRAGLIB +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST LIBRARY LIBRARY_SS LIBRARY2 CP FLUX DISCR DISCR_SS EDITION + EIGHTH EIGHTH_SS MACDIR ; +MODULE + LIB: GEO: SYBILT: USS: ASM: FLU: EDI: SPH: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + EDIT 2 NMIX 6 CTRA APOL SUBG + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + MIX 2 COMB 1 1.0 + MIX 3 COMB 1 1.0 + MIX 4 COMB 1 1.0 + MIX 5 579.9 + Zr90 = Zr0 4.18621E-2 2 + MIX 6 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + ; +*---- +* Geometry : a 15 X 15 normal PWR assembly +* contains T : guide thimble tube +* C1 : normal fuel cell +* C2 : peripheral cell +* C3 : corner cell +*---- +REAL Rfuel := 0.4646 ; +REAL Rfuel1 := Rfuel .50 SQRT * ; +REAL Rfuel2 := Rfuel .80 SQRT * ; +REAL Rfuel3 := Rfuel .95 SQRT * ; + +EIGHTH_SS := GEO: :: CAR2D 8 8 +X- DIAG X+ REFL +Y- SYME Y+ DIAG + ::: T := GEO: CARCEL 2 + MESHX 0. 1.43 + MESHY 0. 1.43 + RADIUS 0. 0.55 0.57 + MIX 6 5 6 ; + ::: C1 := GEO: CARCEL 5 + MESHX 0. 1.43 + MESHY 0. 1.43 + RADIUS 0. <<Rfuel1>> <<Rfuel2>> <<Rfuel3>> <<Rfuel>> 0.53594 + MIX 1 2 3 4 5 6 ; + ::: C2 := GEO: C1 MESHX 0. 1.48 MIX 1 2 3 4 5 6 ; + ::: C3 := GEO: C2 MESHY 0. 1.48 MIX 1 2 3 4 5 6 ; +CELL T C1 C1 C1 T C1 C1 C2 + C1 C1 C1 C1 C1 C1 C2 + C1 C1 C1 T C1 C2 + T C1 C1 C1 C2 + C1 C1 C1 C2 + T C1 C2 + C1 C2 + C3 +MERGE 1 2 2 2 1 2 2 4 + 2 2 2 2 2 2 4 + 2 2 2 1 2 4 + 1 2 2 2 4 + 2 2 2 4 + 1 2 4 + 2 4 + 3 +; +*---- +* Self-Shielding calculation SYBIL +*---- +DISCR_SS := SYBILT: EIGHTH_SS :: MAXZ 1000000 MAXR 1000000 + QUA2 3 3 GAUSS DP01 ; +LIBRARY_SS := USS: LIBRARY DISCR_SS ; +*---- +* Permute mixtures before flux calculation +*---- +LIBRARY2 := LIB: LIBRARY_SS :: + EDIT 2 NMIX 6 + CATL + MIX 1 6 + MIX 2 5 + MIX 3 4 + MIX 4 3 + MIX 5 2 + MIX 6 1 + ; +*---- +* Transport calculation SYBIL +* Flux calculation for B1 homogeneous leakage +*---- +EIGHTH := GEO: :: CAR2D 8 8 +X- DIAG X+ REFL +Y- SYME Y+ DIAG + ::: T := GEO: CARCEL 2 + MESHX 0. 1.43 + MESHY 0. 1.43 + RADIUS 0. 0.55 0.57 + MIX 1 2 1 ; + ::: C1 := GEO: CARCEL 5 + MESHX 0. 1.43 + MESHY 0. 1.43 + RADIUS 0. <<Rfuel1>> <<Rfuel2>> <<Rfuel3>> <<Rfuel>> 0.53594 + MIX 6 5 4 3 2 1 ; + ::: C2 := GEO: C1 MESHX 0. 1.48 MIX 6 5 4 3 2 1 ; + ::: C3 := GEO: C2 MESHY 0. 1.48 MIX 6 5 4 3 2 1 ; +CELL T C1 C1 C1 T C1 C1 C2 + C1 C1 C1 C1 C1 C1 C2 + C1 C1 C1 T C1 C2 + T C1 C1 C1 C2 + C1 C1 C1 C2 + T C1 C2 + C1 C2 + C3 +; +DISCR := SYBILT: EIGHTH :: MAXZ 1000000 MAXR 1000000 + QUA2 3 3 GAUSS DP01 ; +CP := ASM: LIBRARY2 DISCR :: ARM ; +FLUX := FLU: CP LIBRARY2 DISCR :: EDIT 1 TYPE B B1 SIGS ; +EDITION := EDI: LIBRARY2 DISCR FLUX EIGHTH :: + EDIT 1 SAVE MICR RES MERGE COMP COND 0.625 JOUT ; +EDITION := SPH: EDITION :: EDIT 2 ALBS SELE_ALB ; +* +MACDIR := EDITION :: STEP UP 'REF-CASE0001' STEP UP MACROLIB + STEP UP GROUP STEP AT 1 ; +assertS MACDIR :: 'NTOT0' 1 5.305629E-01 ; +assertS MACDIR :: 'NUSIGF' 1 6.686473E-03 ; +assertS MACDIR :: 'DIFF' 1 1.400524E+00 ; +assertS MACDIR :: 'SIGW00' 1 5.050910E-01 ; +MACDIR := DELETE: MACDIR ; +MACDIR := EDITION :: STEP UP 'REF-CASE0001' STEP UP MACROLIB + STEP UP GROUP STEP AT 2 ; +assertS MACDIR :: 'NTOT0' 1 1.359428E+00 ; +assertS MACDIR :: 'NUSIGF' 1 1.365678E-01 ; +assertS MACDIR :: 'DIFF' 1 4.055927E-01 ; +assertS MACDIR :: 'SIGW00' 1 1.271332E+00 ; +MACDIR := DELETE: MACDIR ; +ECHO "test TCDR47 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/tdraglib_proc/TCDR48.c2m b/Dragon/data/tdraglib_proc/TCDR48.c2m new file mode 100644 index 0000000..10dda98 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR48.c2m @@ -0,0 +1,43 @@ +*---- +* +* TEST TCDR48 +* CRITICAL FLUX CALCULATION IN PURE WATER +* 172-GROUP JEF2.2 DRAGLIB +* +* Author: A. Hebert +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST GEOM TRACK LIBRARY ASB FLUX ; +MODULE LIB: GEO: SYBILT: ASM: FLU: END: ; +PROCEDURE assertS ; +* +GEOM := GEO: :: HOMOGE + MIX 1 ; +TRACK := SYBILT: GEOM :: + MAXR 20 MAXZ 15000 QUA2 20 3 ; +*---- +* LIB +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 2 (*MAXIMUM OF MATERIAL MIXTURES*) + CTRA APOL (*APOLLO TYPE TRANSPORT CORRECTION*) + ANIS 2 (*PROCESS SCATTERING DATA TO P1 ORDER*) +* + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 293.0 + H1 = H1_H2O 4.6892E-2 + O16 = O16 2.3446E-2 + ; +*---- +* FLUX +*---- +ASB := ASM: LIBRARY TRACK :: PIJ ; +FLUX := FLU: LIBRARY TRACK ASB :: EDIT 2 TYPE L P1 SIGS BUCK -1.0E-2 ; +assertS FLUX :: 'B2 B1HOM' 1 -6.071696E-3 ; +FLUX := FLU: FLUX LIBRARY TRACK ASB :: TYPE L RHS SIGS ; +assertS FLUX :: 'B2 B1HOM' 1 -6.071696E-3 ; +ECHO "test TCDR48 completed" ; +END: ; diff --git a/Dragon/data/tdraglib_proc/TCDR49.c2m b/Dragon/data/tdraglib_proc/TCDR49.c2m new file mode 100644 index 0000000..9086450 --- /dev/null +++ b/Dragon/data/tdraglib_proc/TCDR49.c2m @@ -0,0 +1,124 @@ +*---- +* TEST CASE TRDR49 +* 5 X 5 PWR ASSEMBLY WITH MULTICELL SURFACIC APPROXIMATION +* 172-GROUP JEF2.2 DRAGLIB +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + ASSMB DISCR LIBRARY LIBRARY2 CP CALC ; +MODULE + LIB: GEO: G2S: SYBILT: SALT: USS: ASM: FLU: DELETE: END: ; +PROCEDURE assertS ; +SEQ_ASCII Fic_SAL :: FILE './Fic_SAL_TCDR49.tdt' ; +SEQ_ASCII Fic_EPS.eps :: FILE './Fic_EPS_TCDR49.eps' ; +SEQ_BINARY DISCR_FIL ; +*---- +* Microscopic cross sections from Jef 2.2 Draglib +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + SUBG + MIXS LIB: DRAGON FIL: DLIB_J2 + MIX 1 579.9 + H1H2O = H1_H2O 4.76690E-2 O16H2O = O16 2.38345E-2 + B10 = B10 4.71444E-6 B11 = B11 1.90959E-5 + MIX 2 579.9 O16 = O16 3.06711E-4 + Cr52 = Cr52 7.54987E-5 Fe56 = Fe56 1.47624E-4 + Zr91 = Zr0 4.18621E-2 + MIX 3 579.9 + H1H2O = H1_H2O 4.65292E-2 O16H2O = O16 2.32646E-2 + Cr52 = Cr52 4.79927E-5 + Fe56 = Fe56 4.45845E-5 Ni58 = Ni58 1.13521E-4 + MoNat = Mo95 4.03755E-6 + Al27 = Al27 2.35231E-6 + Mn55 = Mn55 4.15901E-7 + B10 = B10 4.60867E-6 B11 = B11 1.86674E-5 + Zr91 = Zr0 8.92427E-4 + MIX 4 579.9 O16 = O16 2.87335E-4 + Cr52 = Cr52 7.07291E-5 Fe56 = Fe56 1.38298E-4 + Zr91 = Zr0 3.92175E-2 + MIX 5 579.9 + H1H2O = H1_H2O 4.71346E-2 O16H2O = O16 2.35673E-2 + Cr52 = Cr52 2.24991E-5 + Fe56 = Fe56 2.09013E-5 Ni58 = Ni58 5.32188E-5 + MoNat = Mo95 1.89281E-6 + Al27 = Al27 1.10277E-6 + Mn55 = Mn55 1.94976E-7 + B10 = B10 4.66484E-6 B11 = B11 1.88950E-5 + Zr91 = Zr0 4.18372E-4 + MIX 6 579.9 + H1H2O = H1_H2O 4.71676E-2 O16H2O = O16 2.35838E-2 + Cr52 = Cr52 2.11122E-5 + Fe56 = Fe56 1.96130E-5 Ni58 = Ni58 4.99383E-5 + MoNat = Mo95 1.77614E-6 + Al27 = Al27 1.03479E-6 + Mn55 = Mn55 1.82957E-7 + B10 = B10 4.66791E-6 B11 = B11 1.89074E-5 + Zr91 = Zr0 3.92583E-4 + MIX 7 579.9 + H1H2O = H1_H2O 4.72020E-2 O16H2O = O16 2.36010E-2 + Cr52 = Cr52 1.96591E-5 + Fe56 = Fe56 1.82630E-5 Ni58 = Ni58 4.65011E-5 + MoNat = Mo95 1.65389E-6 + Al27 = Al27 9.63569E-7 + Mn55 = Mn55 1.70365E-7 + B10 = B10 4.67110E-6 B10 = B10 1.89203E-5 + Zr91 = Zr0 3.65562E-4 + MIX 8 933.6 O16 = O16 4.49355E-2 + U235 = U235 7.39237E-4 1 + U238 = U238 2.17285E-2 1 + ; +*---- +* Geometry ASSMB : a 17 X 17 normal PWR assembly +* contains C1 : cell without fuel +* C2 : normal fuel cell +* C3 : peripheral cell +* C4 : corner cell +*---- +ASSMB := GEO: :: CAR2D 3 3 + X- DIAG X+ ALBE 1.0 Y- SYME Y+ DIAG + CELL C1 C2 C3 + C2 C3 + C4 + ::: C1 := GEO: CARCEL 2 + MESHX 0.0 1.26472 MESHY 0.0 1.26472 + RADIUS 0.0 0.572435 0.613142 MIX 1 2 3 ; + ::: C2 := GEO: C1 RADIUS 0.0 0.412660 0.474364 MIX 8 4 5 ; + ::: C3 := GEO: C2 MESHX 0.0 1.31472 MIX 8 4 6 ; + ::: C4 := GEO: C3 MESHY 0.0 1.31472 MIX 8 4 7 ; + ; +*---- +* Self-Shielding calculation SYBIL +* Transport calculation SYBIL +* Flux calculation for B1 homogeneous leakage +*---- +DISCR := SYBILT: ASSMB :: + TITLE 'TCDR49: 5 X 5 MULTICELL 5x5 CASE WITH SYBIL' + MAXR 100 QUA2 6 3 ; +LIBRARY2 := USS: LIBRARY DISCR :: EDIT 2 PASS 2 ARM ; +CP := ASM: LIBRARY2 DISCR :: ARM ; +CALC := FLU: CP LIBRARY2 DISCR :: TYPE K B1 SIGS ; +assertS CALC :: 'K-INFINITY' 1 1.233258 ; +DISCR LIBRARY2 CP CALC := DELETE: DISCR LIBRARY2 CP CALC ; +*---- +* Self-Shielding calculation SALT WITH MULTICELL SURFACIC +* Transport calculation SALT WITH MULTICELL SURFACIC +* Flux calculation for B1 homogeneous leakage +*---- +Fic_SAL Fic_EPS.eps := G2S: ASSMB :: MACRO ; + +DISCR DISCR_FIL := SALT: Fic_SAL :: + EDIT 3 + TITLE "TCDR49: MULTICELL SURFACIC 5x5 CASE" + TISO 8 40.0 + IC EPSJ 1.0E-5 + ; +LIBRARY2 := USS: LIBRARY DISCR DISCR_FIL :: EDIT 2 PASS 2 ARM ; +CP := ASM: LIBRARY2 DISCR DISCR_FIL :: EDIT 1 ARM ; +CALC := FLU: CP LIBRARY2 DISCR :: TYPE K B1 SIGS ; +assertS CALC :: 'K-INFINITY' 1 1.232764 ; +ECHO "test TCDR49 completed" ; +END: ; +QUIT "LIST" . |