From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001 From: stainer_t Date: Mon, 8 Sep 2025 13:48:49 +0200 Subject: Initial commit from Polytechnique Montreal --- Dragon/data/tdraglib_proc/TCDR20.c2m | 175 +++++++++++++++++++++++++++++++++++ 1 file changed, 175 insertions(+) create mode 100644 Dragon/data/tdraglib_proc/TCDR20.c2m (limited to 'Dragon/data/tdraglib_proc/TCDR20.c2m') 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 <> +GRA = 'C0_GR' 3.5353E-02 +C0 = 'C0' 5.0874E-02 +B10 = 'B10' <> +B11 = 'B11' <> +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 <> +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 <> +C0 = 'C0' 7.5811E-09 +N14 = 'N14' 3.9484E-05 +O16 = 'O16' 1.0608E-05 +*-- Reflector +MIX 4 <> +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 <> +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" . -- cgit v1.2.3