*---- * TEST CASE TDCM33 * MACROSCOPIC CROSS SECTIONS * PULSE REACTOR WITH MOC *---- * Define STRUCTURES and MODULES used *---- LINKED_LIST GEOM TRACK MACRO PIJ FLUX EDITION ; MODULE GEO: G2S: SALT: MCCGT: MAC: ASM: FLU: EDI: GREP: END: ; SEQ_BINARY TRACK_FIL ; SEQ_ASCII Fic_SAL :: FILE './Fic_SAL.tdt' ; SEQ_ASCII Fic_SAL.eps :: FILE './Fic_SAL.eps' ; REAL Kinf ; INTEGER an2d := 6 ; REAL densur := 15.0 ; PROCEDURE assertS ; * * evaluate pin positions in triangular assemblies REAL x y := 12.6 13.91 ; REAL rpin1 := x x * y y * + SQRT ; REAL apin1 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin2 := x x * y y * + SQRT ; REAL apin2 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin3 := x x * y y * + SQRT ; REAL apin3 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin4 := x x * y y * + SQRT ; REAL apin4 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin5 := x x * y y * + SQRT ; REAL apin5 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin6 := x x * y y * + SQRT ; REAL apin6 := y x / ARCTAN ; EVALUATE x y := 11.34 13.91 ; REAL rpin7 := x x * y y * + SQRT ; REAL apin7 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin8 := x x * y y * + SQRT ; REAL apin8 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin9 := x x * y y * + SQRT ; REAL apin9 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin10 := x x * y y * + SQRT ; REAL apin10 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin11 := x x * y y * + SQRT ; REAL apin11 := y x / ARCTAN ; EVALUATE x y := 10.08 13.91 ; REAL rpin12 := x x * y y * + SQRT ; REAL apin12 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin13 := x x * y y * + SQRT ; REAL apin13 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin14 := x x * y y * + SQRT ; REAL apin14 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin15 := x x * y y * + SQRT ; REAL apin15 := y x / ARCTAN ; EVALUATE x y := 8.82 13.91 ; REAL rpin16 := x x * y y * + SQRT ; REAL apin16 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin17 := x x * y y * + SQRT ; REAL apin17 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin18 := x x * y y * + SQRT ; REAL apin18 := y x / ARCTAN ; EVALUATE x y := 7.56 13.91 ; REAL rpin19 := x x * y y * + SQRT ; REAL apin19 := y x / ARCTAN ; EVALUATE y := y 1.26 - ; REAL rpin20 := x x * y y * + SQRT ; REAL apin20 := y x / ARCTAN ; EVALUATE x y := 6.3 13.91 ; REAL rpin21 := x x * y y * + SQRT ; REAL apin21 := y x / ARCTAN ; * REAL apin22 := $Pi_R apin1 + CHS ; REAL apin23 := $Pi_R apin2 + CHS ; REAL apin24 := $Pi_R apin3 + CHS ; REAL apin25 := $Pi_R apin4 + CHS ; REAL apin26 := $Pi_R apin5 + CHS ; REAL apin27 := $Pi_R apin6 + CHS ; REAL apin28 := $Pi_R apin7 + CHS ; REAL apin29 := $Pi_R apin8 + CHS ; REAL apin30 := $Pi_R apin9 + CHS ; REAL apin31 := $Pi_R apin10 + CHS ; REAL apin32 := $Pi_R apin11 + CHS ; REAL apin33 := $Pi_R apin12 + CHS ; REAL apin34 := $Pi_R apin13 + CHS ; REAL apin35 := $Pi_R apin14 + CHS ; REAL apin36 := $Pi_R apin15 + CHS ; REAL apin37 := $Pi_R apin16 + CHS ; REAL apin38 := $Pi_R apin17 + CHS ; REAL apin39 := $Pi_R apin18 + CHS ; REAL apin40 := $Pi_R apin19 + CHS ; REAL apin41 := $Pi_R apin20 + CHS ; REAL apin42 := $Pi_R apin21 + CHS ; * REAL apin43 := $Pi_R apin1 + ; REAL apin44 := $Pi_R apin2 + ; REAL apin45 := $Pi_R apin3 + ; REAL apin46 := $Pi_R apin4 + ; REAL apin47 := $Pi_R apin5 + ; REAL apin48 := $Pi_R apin6 + ; REAL apin49 := $Pi_R apin7 + ; REAL apin50 := $Pi_R apin8 + ; REAL apin51 := $Pi_R apin9 + ; REAL apin52 := $Pi_R apin10 + ; REAL apin53 := $Pi_R apin11 + ; REAL apin54 := $Pi_R apin12 + ; REAL apin55 := $Pi_R apin13 + ; REAL apin56 := $Pi_R apin14 + ; REAL apin57 := $Pi_R apin15 + ; REAL apin58 := $Pi_R apin16 + ; REAL apin59 := $Pi_R apin17 + ; REAL apin60 := $Pi_R apin18 + ; REAL apin61 := $Pi_R apin19 + ; REAL apin62 := $Pi_R apin20 + ; REAL apin63 := $Pi_R apin21 + ; * REAL apin64 := apin1 CHS ; REAL apin65 := apin2 CHS ; REAL apin66 := apin3 CHS ; REAL apin67 := apin4 CHS ; REAL apin68 := apin5 CHS ; REAL apin69 := apin6 CHS ; REAL apin70 := apin7 CHS ; REAL apin71 := apin8 CHS ; REAL apin72 := apin9 CHS ; REAL apin73 := apin10 CHS ; REAL apin74 := apin11 CHS ; REAL apin75 := apin12 CHS ; REAL apin76 := apin13 CHS ; REAL apin77 := apin14 CHS ; REAL apin78 := apin15 CHS ; REAL apin79 := apin16 CHS ; REAL apin80 := apin17 CHS ; REAL apin81 := apin18 CHS ; REAL apin82 := apin19 CHS ; REAL apin83 := apin20 CHS ; REAL apin84 := apin21 CHS ; * REAL Rtest4 := 0.4083 ; REAL Rtest1 := 0.5 SQRT Rtest4 * ; REAL Rtest2 := 0.8 SQRT Rtest4 * ; REAL Rtest3 := 0.95 SQRT Rtest4 * ; REAL Rcomb4 := 0.46 ; REAL Rcomb1 := 0.5 SQRT Rcomb4 * ; REAL Rcomb2 := 0.8 SQRT Rcomb4 * ; REAL Rcomb3 := 0.95 SQRT Rcomb4 * ; * GEOM := GEO: :: CARCEL 19 X- VOID X+ VOID Y- VOID Y+ VOID RADIUS 0.0 <> <> <> <> 0.475 0.75 1.55 2.4 3.92 4.1 6.1 7.7 8.47 12.5 12.9 13.75 14.9 16.2 17.55 MESHX 0.0 26.46 MESHY 0.0 29.08 SPLITX 6 SPLITY 6 MIX 7 8 9 10 11 12 13 14 15 16 17 18 19 7 20 7 7 7 7 7 CLUSTER PIN1 PIN2 PIN3 PIN4 PIN5 PIN6 PIN7 PIN8 PIN9 PIN10 PIN11 PIN12 PIN13 PIN14 PIN15 PIN16 PIN17 PIN18 PIN19 PIN20 PIN21 PIN22 PIN23 PIN24 PIN25 PIN26 PIN27 PIN28 PIN29 PIN30 PIN31 PIN32 PIN33 PIN34 PIN35 PIN36 PIN37 PIN38 PIN39 PIN40 PIN41 PIN42 PIN43 PIN44 PIN45 PIN46 PIN47 PIN48 PIN49 PIN50 PIN51 PIN52 PIN53 PIN54 PIN55 PIN56 PIN57 PIN58 PIN59 PIN60 PIN61 PIN62 PIN63 PIN64 PIN65 PIN66 PIN67 PIN68 PIN69 PIN70 PIN71 PIN72 PIN73 PIN74 PIN75 PIN76 PIN77 PIN78 PIN79 PIN80 PIN81 PIN82 PIN83 PIN84 ::: PIN1 := GEO: TUBE 5 RADIUS 0.0 <> <> <> <> 0.5 MIX 1 2 3 4 5 NPIN 1 RPIN <> APIN <> ; ::: PIN2 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN3 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN4 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN5 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN6 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN7 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN8 := GEO: TUBE 1 RADIUS 0.0 0.5 MIX 21 NPIN 1 RPIN <> APIN <> ; ::: PIN9 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN10 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN11 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN12 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN13 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN14 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN15 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN16 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN17 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN18 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN19 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN20 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN21 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; * ::: PIN22 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN23 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN24 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN25 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN26 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN27 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN28 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN29 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN30 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN31 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN32 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN33 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN34 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN35 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN36 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN37 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN38 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN39 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN40 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN41 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN42 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; * ::: PIN43 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN44 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN45 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN46 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN47 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN48 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN49 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN50 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN51 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN52 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN53 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN54 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN55 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN56 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN57 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN58 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN59 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN60 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN61 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN62 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN63 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; * ::: PIN64 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN65 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN66 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN67 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN68 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN69 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN70 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN71 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN72 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN73 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN74 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN75 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN76 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN77 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN78 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN79 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN80 := GEO: PIN8 NPIN 1 RPIN <> APIN <> ; ::: PIN81 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN82 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN83 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ::: PIN84 := GEO: PIN1 NPIN 1 RPIN <> APIN <> ; ; Fic_SAL Fic_SAL.eps := G2S: GEOM ; *---- * Macroscopic XS *---- MACRO := MAC: :: EDIT 2 NGRO 1 NMIX 21 NIFI 1 READ INPUT MIX 1 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 2 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 3 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 4 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 5 (*CLAD*) TOTAL 0.4029 SCAT 1 1 0.4000 MIX 6 (*H2O*) TOTAL 0.3683 SCAT 1 1 0.3661 MIX 7 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 8 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 9 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 10 (*FUEL*) TOTAL 0.3652 SCAT 1 1 0.3234 NUSIGF 0.089 CHI 1.0 MIX 11 (*CLAD*) TOTAL 0.4029 SCAT 1 1 0.4000 MIX 12 (*H2O*) TOTAL 0.3683 SCAT 1 1 0.3661 MIX 13 (*TUBE*) TOTAL 0.6029 SCAT 1 1 0.6000 MIX 14 (*GAS*) TOTAL 0.0650 SCAT 1 1 0.0600 MIX 15 (*PRESSURE TUBE*) TOTAL 0.4029 SCAT 1 1 0.4000 MIX 16 (*H2O*) TOTAL 0.3683 SCAT 1 1 0.3661 MIX 17 (*BETA TUBE*) TOTAL 0.4029 SCAT 1 1 0.4000 MIX 18 (*VOID*) TOTAL 0.006 SCAT 1 1 0.0060 MIX 19 (*ALPHA TUBE*) TOTAL 0.4029 SCAT 1 1 0.4000 MIX 20 (*NICKEL TUBE*) TOTAL 0.4029 SCAT 1 1 0.4000 MIX 21 (*BETA TUBE*) TOTAL 0.4029 SCAT 1 1 0.4000 ; *---- * Perform tracking *---- TRACK TRACK_FIL := SALT: Fic_SAL :: EDIT 1 TISO <> <> LONG ; TRACK := MCCGT: TRACK TRACK_FIL :: CACB 4 AAC 80 TMT EPSI 1E-5 MCU 100000 ; *---- * Compute flux *---- PIJ := ASM: MACRO TRACK TRACK_FIL :: ARM ; FLUX := FLU: PIJ MACRO TRACK TRACK_FIL :: TYPE K ; GREP: FLUX :: GETVAL 'K-INFINITY ' 1 1 1 >>Kinf<< ; ECHO "@@@ Kinf=" Kinf ; EDITION := EDI: FLUX MACRO TRACK Fic_SAL :: EDIT 3 MERG G2S 36 RECT 0.0 4.41 0.0 4.846666667 RECT 0.0 4.41 4.846666667 9.693333333 RECT 0.0 4.41 9.693333333 14.54 RECT 0.0 4.41 14.54 19.386666667 RECT 0.0 4.41 19.386666667 24.233333333 RECT 0.0 4.41 24.233333333 29.08 RECT 4.41 8.82 0.0 4.846666667 RECT 4.41 8.82 4.846666667 9.693333333 RECT 4.41 8.82 9.693333333 14.54 RECT 4.41 8.82 14.54 19.386666667 RECT 4.41 8.82 19.386666667 24.233333333 RECT 4.41 8.82 24.233333333 29.08 RECT 8.82 13.23 0.0 4.846666667 RECT 8.82 13.23 4.846666667 9.693333333 RECT 8.82 13.23 9.693333333 14.54 RECT 8.82 13.23 14.54 19.386666667 RECT 8.82 13.23 19.386666667 24.233333333 RECT 8.82 13.23 24.233333333 29.08 RECT 13.23 17.64 0.0 4.846666667 RECT 13.23 17.64 4.846666667 9.693333333 RECT 13.23 17.64 9.693333333 14.54 RECT 13.23 17.64 14.54 19.386666667 RECT 13.23 17.64 19.386666667 24.233333333 RECT 13.23 17.64 24.233333333 29.08 RECT 17.64 22.05 0.0 4.846666667 RECT 17.64 22.05 4.846666667 9.693333333 RECT 17.64 22.05 9.693333333 14.54 RECT 17.64 22.05 14.54 19.386666667 RECT 17.64 22.05 19.386666667 24.233333333 RECT 17.64 22.05 24.233333333 29.08 RECT 22.05 26.46 0.0 4.846666667 RECT 22.05 26.46 4.846666667 9.693333333 RECT 22.05 26.46 9.693333333 14.54 RECT 22.05 26.46 14.54 19.386666667 RECT 22.05 26.46 19.386666667 24.233333333 RECT 22.05 26.46 24.233333333 29.08 SAVE ON TEST ; assertS FLUX :: 'K-EFFECTIVE' 1 1.036426 ; ECHO "test TDCM33 completed" ; END: ; QUIT "LIST" .