*---- * 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 LIBRARY2 CP CALC OUT SPHGEOM MTRACK MACRO MACRO2 SYS FLUX EDIT OPTIM ; MODULE LIB: GEO: SYBILT: USS: ASM: FLU: EDI: FLUD: MAC: DELETE: UTL: BIVACT: BIVACA: OUT: ERROR: END: ; PROCEDURE SPH_FixedPt ; PROCEDURE assertS ; INTEGER IterEmax := 50 ; *---- * Microscopic cross sections from Jef 2.2 Draglib *---- LIBRARY := LIB: :: NMIX 8 CTRA WIMS ANIS 2 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-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 ; LIBRARY2 := USS: LIBRARY DISCR :: PASS 2 ; 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: LIBRARY2 DISCR :: ARM ; CALC := FLU: CP LIBRARY2 DISCR :: TYPE K P1 HETE ; assertS CALC :: 'K-EFFECTIVE' 1 0.4565791 ; *---- * Editing using SPH model for transport-diffusion with flux-volume * normalization *---- OUT := EDI: LIBRARY2 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 fixed-point SPH equivalence *---- MACRO2 OPTIM := SPH_FixedPt MACRO SPHGEOM MTRACK :: 0.0 10.0 1.0E-7 <> ; ECHO "control variables at convergence" ; UTL: OPTIM :: IMPR 'VAR-VALUE' * ; *---- * Verification calculation *---- CALC := DELETE: CALC ; SYS := BIVACA: MACRO2 MTRACK ; CALC := FLUD: SYS MTRACK :: EDIT 2 ADI 5 ACCE 5 3 ; assertS CALC :: K-EFFECTIVE 1 0.4565792 ; ECHO "test TCDR22 completed" ; END: ; QUIT "LIST" .