diff options
Diffstat (limited to 'Dragon/data/tdraglib_proc/TCDR22.c2m')
| -rw-r--r-- | Dragon/data/tdraglib_proc/TCDR22.c2m | 208 |
1 files changed, 208 insertions, 0 deletions
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" . |