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/twlup_proc | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Dragon/data/twlup_proc')
24 files changed, 3560 insertions, 0 deletions
diff --git a/Dragon/data/twlup_proc/Datos_pertur.c2m b/Dragon/data/twlup_proc/Datos_pertur.c2m new file mode 100644 index 0000000..a67fa1a --- /dev/null +++ b/Dragon/data/twlup_proc/Datos_pertur.c2m @@ -0,0 +1,101 @@ +*ENTRADA Datos_pertur.c2m +*---- +* Nombre : Datos_pertur.c2m +* Autor : R. Santi +* + +MODULE + END: ; + +LINKED_LIST Datos_pertur ; + +* Dimension + +REAL radio := 14.302 ; ! circular + +* Temperaturas, densidades y composiciones + +REAL PUR_MOL := 99.60 ; +REAL conc_B := 0.05 ; +REAL conc_Li := 0.26 ; + +REAL T_comb1 := 798.54 ; +REAL T_comb2 := 798.54 ; +REAL T_comb3 := 798.54 ; +REAL T_comb4 := 798.54 ; + +REAL Tvaina := 564.53 ; +REAL Dzry := 6.56 ; +REAL T_cool := 559.03 ; +REAL D_cool := 0.82430 ; +REAL Tcanal := 550.25 ; +REAL T_gap := 515.76 ; +REAL D_gap := 0.90467 ; +REAL T_foil := 485.20 ; +REAL T_mod := 448.66 ; +REAL D_mod := 0.99524 ; + +REAL D_comb := 10.07189 ; + +REAL W238 := 87.51763 ; +REAL W235 := 0.62707 ; +REAL W234 := 5.1127E-3 ; +REAL W239 := 1.0E-18 ; +REAL WO2 := 11.8502 ; +REAL WB := 1.13E-5 ; +REAL potencia := 21.8940148 ; + +* Temperaturas, densidades y composiciones perturbadas + +REAL PUR_MOL_P := 98.60 ; +REAL conc_B_p := 1.0 ; +REAL conc_Li_p := 1.0 ; +REAL T_comb1P := 960.15 ; +REAL T_comb2P := 960.15 ; +REAL T_comb3P := 960.15 ; +REAL T_comb4P := 960.15 ; +REAL T_modP := 460.65 ; +REAL T_coolP := 578.95 ; +REAL D_modP := 0.98144 ; +REAL D_coolP := 0.78430 ; +REAL D_coolV := 0.0013 ; + +* BUCKLING GEOMETRICO + +REAL buck_geom := 7.328175E-05 ; + +* SELECCION DE LA VARIABLE A PERTURBAR ("perturba"): +* a) Si la entrada es "da2_pertur_Qwims_T_Dens.inp" +* 1 Tcombustible +* 2 Trefrigerante +* 3 Dens refrigerante +* 4 Trefrig y Dens refrig +* 5 Tmoderador +* 6 Dens moderador +* 7 Tmoderador y Dens moderador +* 8 Vacio + +* b) Si la entrada es "da2_pertur_Qwims_D2O_B_Li.inp" +* 9 Pureza D2O +* 10 Boro en refrigerante y moderador +* 11 Litio en refrigerante y moderador + +INTEGER perturba := 3 ; ! opcion 1 a 8 + +* SELECCION DEL TIPO DE BUCKLING A USAR EN LA REFERENCIA ("buckling"): + +STRING buckling := "geometrico" ; ! opcion "critico" // "geometrico" + +:: <<radio>> <<PUR_MOL>> <<conc_B>> <<conc_Li>> + <<T_comb1>> <<T_comb2>> <<T_comb3>> <<T_comb4>> + <<Tvaina>> <<Dzry>> <<T_cool>> <<D_cool>> + <<Tcanal>> <<T_gap>> <<D_gap>> <<T_foil>> <<T_mod>> <<D_mod>> + <<D_comb>> <<W238>> <<W235>> <<W234>> <<W239>> <<WO2>> <<WB>> + <<potencia>> <<PUR_MOL_P>> <<conc_B_p>> <<conc_Li_p>> + <<T_comb1P>> <<T_comb2P>> <<T_comb3P>> <<T_comb4P>> + <<T_modP>> <<T_coolP>> <<D_modP>> <<D_coolP>> <<D_coolV>> + <<buck_geom>> <<perturba>> <<buckling>> + ; + +END: ; +QUIT . diff --git a/Dragon/data/twlup_proc/TCWU01.c2m b/Dragon/data/twlup_proc/TCWU01.c2m new file mode 100644 index 0000000..9f9f57b --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU01.c2m @@ -0,0 +1,115 @@ +*---- +* TEST CASE TCWU01 +* MOSTELLER BENCHMARK: 1-D ANNULAR CELL AND 2-D CARTESIAN CELL +* WIMSD4 69 GROUPS LIBRARY FILE iaea FROM WLUP +* +* 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 file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 600.0 O16 = '6016' 4.61309E-2 + U235 = '2235' 1.66078E-4 1 + U238 = '8238' 2.28994E-2 1 + MIX 2 600.0 + Zr91 = '91' 3.83243E-2 + MIX 3 600.0 + H1H2O = '3001' 4.42326E-2 O16H2O = '6016' 2.21163E-2 + BNat = '1011' 1.02133E-5 + ; +*---- +* 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 'TCWU01: 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.8276153 ; +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 'TCWU01: 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.8277465 ; +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 'TCWU01: 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.8287268 ; +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 'TCWU01: 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.8307597 ; +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 TCWU01 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU02.c2m b/Dragon/data/twlup_proc/TCWU02.c2m new file mode 100644 index 0000000..e6df980 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU02.c2m @@ -0,0 +1,159 @@ +*---- +* TEST CASE TCWU02 +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* WIMSD4 69 GROUPS LIBRARY FILE iaea FROM WLUP +* +* 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 file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 579.9 + H1H2O = '3001' 4.76690E-2 O16H2O = '6016' 2.38345E-2 + BNat = '1011' 2.38103E-5 + MIX 2 579.9 O16 = '6016' 3.06711E-4 + Cr52 = '52' 7.54987E-5 Fe56 = '2056' 1.47624E-4 + Zr91 = '91' 4.18621E-2 + MIX 3 579.9 + H1H2O = '3001' 4.65292E-2 O16H2O = '6016' 2.32646E-2 + Cr52 = '52' 4.79927E-5 + Fe56 = '2056' 4.45845E-5 Ni58 = '58' 1.13521E-4 + MoNat = '96' 4.03755E-6 + Al27 = '27' 2.35231E-6 + Mn55 = '55' 4.15901E-7 + BNat = '1011' 2.32761E-5 + Zr91 = '91' 8.92427E-4 + MIX 4 579.9 O16 = '6016' 2.87335E-4 + Cr52 = '52' 7.07291E-5 Fe56 = '2056' 1.38298E-4 + Zr91 = '91' 3.92175E-2 + MIX 5 579.9 + H1H2O = '3001' 4.71346E-2 O16H2O = '6016' 2.35673E-2 + Cr52 = '52' 2.24991E-5 + Fe56 = '2056' 2.09013E-5 Ni58 = '58' 5.32188E-5 + MoNat = '96' 1.89281E-6 + Al27 = '27' 1.10277E-6 + Mn55 = '55' 1.94976E-7 + BNat = '1011' 2.35598E-5 + Zr91 = '91' 4.18372E-4 + MIX 6 579.9 + H1H2O = '3001' 4.71676E-2 O16H2O = '6016' 2.35838E-2 + Cr52 = '52' 2.11122E-5 + Fe56 = '2056' 1.96130E-5 Ni58 = '58' 4.99383E-5 + MoNat = '96' 1.77614E-6 + Al27 = '27' 1.03479E-6 + Mn55 = '55' 1.82957E-7 + BNat = '1011' 2.35753E-5 + Zr91 = '91' 3.92583E-4 + MIX 7 579.9 + H1H2O = '3001' 4.72020E-2 O16H2O = '6016' 2.36010E-2 + Cr52 = '52' 1.96591E-5 + Fe56 = '2056' 1.82630E-5 Ni58 = '58' 4.65011E-5 + MoNat = '96' 1.65389E-6 + Al27 = '27' 9.63569E-7 + Mn55 = '55' 1.70365E-7 + BNat = '1011' 2.35914E-5 + Zr91 = '91' 3.65562E-4 + MIX 8 933.6 O16 = '6016' 4.49355E-2 + U235 = '2235' 7.39237E-4 1 + U238 = '8238' 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 +* Editing using SPH model for transport-diffusion +*---- +DISCR := SYBILT: ASSMB :: + TITLE 'TCWU02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 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 1.257190 ; + +OUT := EDI: LIBRARY DISCR CALC ASSMB :: + EDIT 3 UPS SAVE MICR RES MERGE CELL COND 4.0 + ; +SPHGEOM := OUT :: STEP UP 'MACRO-GEOM' ; +MTRACK := BIVACT: SPHGEOM :: PRIM 1 2 EDIT 0 ; +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 5.186729E-01 ; +assertS ISOT :: 'NWT0' 2 1.305927E-01 ; + +ECHO "test TCWU02 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU03.c2m b/Dragon/data/twlup_proc/TCWU03.c2m new file mode 100644 index 0000000..17d501f --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU03.c2m @@ -0,0 +1,150 @@ +*---- +* TEST CASE TCWU03 +* MULTICELL HEXAGONAL ASSEMBLY WITH POISON +* iaea WLUP Library +* +* 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 file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 11 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 579.9 + H1H2O = '3001' 4.76690E-2 O16H2O = '6016' 2.38345E-2 + BNat = '1011' 2.38103E-5 + MIX 2 579.9 O16 = '6016' 3.06711E-4 + Cr52 = '52' 7.54987E-5 Fe56 = '2056' 1.47624E-4 + Zr91 = '91' 4.18621E-2 + MIX 3 579.9 + H1H2O = '3001' 4.65292E-2 O16H2O = '6016' 2.32646E-2 + Cr52 = '52' 4.79927E-5 + Fe56 = '2056' 4.45845E-5 Ni58 = '58' 1.13521E-4 + MoNat = '96' 4.03755E-6 + Al27 = '27' 2.35231E-6 + Mn55 = '55' 4.15901E-7 + BNat = '1011' 2.32761E-5 + Zr91 = '91' 8.92427E-4 + MIX 4 933.6 O16 = '6016' 4.49355E-2 + U235 = '2235' 7.39237E-4 1 + U238 = '8238' 2.17285E-2 1 + MIX 5 579.9 In115 = '2115' 7.57464E-3 + Cd113 = '2113' 2.62493E-3 + Ag109 = '3109' 4.49188E-2 + MIX 6 579.9 Cr52 = '52' 1.52702E-2 + Fe56 = '2056' 5.57670E-2 Ni58 = '58' 7.51418E-3 + Mn55 = '55' 8.02943E-4 + MIX 7 579.9 + H1H2O = '3001' 3.06466E-2 O16H2O = '6016' 1.53233E-2 + Fe56 = '2056' 5.27485E-5 Cr52 = '52' 2.69769E-5 + BNat = '1011' 1.53077E-5 + Zr91 = '91' 1.49580E-2 + MIX 8 579.9 + H1H2O = '3001' 4.65292E-2 O16H2O = '6016' 2.32646E-2 + Cr52 = '52' 4.79927E-5 + Fe56 = '2056' 4.45845E-5 Ni58 = '58' 1.13521E-4 + MoNat = '96' 4.03755E-6 + Al27 = '27' 2.35231E-6 + Mn55 = '55' 4.15901E-7 + BNat = '1011' 2.32761E-5 + Zr91 = '91' 8.92427E-4 + MIX 9 579.9 O16 = '6016' 2.87335E-4 + Cr52 = '52' 7.07291E-5 Fe56 = '2056' 1.38298E-4 + Zr91 = '91' 3.92175E-2 + MIX 10 579.9 + H1H2O = '3001' 4.71346E-2 O16H2O = '6016' 2.35673E-2 + Cr52 = '52' 2.24991E-5 + Fe56 = '2056' 2.09013E-5 Ni58 = '58' 5.32188E-5 + MoNat = '96' 1.89281E-6 + Al27 = '27' 1.10277E-6 + Mn55 = '55' 1.94976E-7 BNat = '1011' 2.35598E-5 + Zr91 = '91' 4.18372E-4 + MIX 11 579.9 + H1H2O = '3001' 4.71676E-2 O16H2O = '6016' 2.35838E-2 + Cr52 = '52' 2.11122E-5 + Fe56 = '2056' 1.96130E-5 Ni58 = '58' 4.99383E-5 + MoNat = '96' 1.77614E-6 + Al27 = '27' 1.03479E-6 + Mn55 = '55' 1.82957E-7 + BNat = '1011' 2.35753E-5 + Zr91 = '91' 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 JPM +* Transport calculation SYBIL +* Flux calculation for B1 homogeneous leakage +* Editing using SPH model for transport-diffusion +*---- +DISCR := SYBILT: ASSMBH :: + TITLE 'TCWU03: MULTICELL HEXAGONAL ASSEMBLY WITH POISON' + MAXR 400 MAXZ 15000 QUA2 6 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.7116683 ; +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.442307E+00 ; +assertS ISOT :: 'NWT0' 2 1.477404E-01 ; + +ECHO "test TCWU03 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU04.c2m b/Dragon/data/twlup_proc/TCWU04.c2m new file mode 100644 index 0000000..9f44916 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU04.c2m @@ -0,0 +1,149 @@ +*---- +* TEST CASE TCWU04 +* iaea WLUP Library +* ANNULAR MOSTELLER BENCHMARK WITH BURNUP +* +* REF: R. Mosteller et al. Nucl. Sci. Eng. 107, 265 (1991) +* +*---- +* Define variables +*---- +INTEGER + istep := 1 ; +REAL + evobeg evoend ; +REAL + step2 step3 step4 step5 step6 step7 := + 1.0 27.1739 47.5544 67.9348 101.9022 135.8696 ; +*---- +* Define STRUCTURES and MODULES used +*---- +PARAMETER res :: ::: SEQ_ASCII res ; ; +LINKED_LIST + LIBRARY MOSTELAS MOSTELA TRACKS TRACK SYS FLUX BURNUP EDITION + DATABASE ISOT ; +MODULE + GEO: SYBILT: LIB: SHI: ASM: FLU: EVO: EDI: COMPO: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + DEPL LIB: WIMSD4 FIL: iaea + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 600.0 O16 = '6016' 4.61309E-2 + U235 = '2235' 1.66078E-4 1 + U238 = '8238' 2.28994E-2 1 + U236 = '236' 0.0 1 + Pu239 = '6239' 0.0 1 + MIX 2 600.0 + Zr91 = '91' 3.83243E-2 + MIX 3 600.0 + H1H2O = '3001' 4.42326E-2 O16H2O = '6016' 2.21163E-2 + BNat = '1011' 1.02133E-5 + ; +*---- +* 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 'TCWU04: MOSTELLER BENCHMARK WITH BURNUP' + EDIT 1 MAXR 3 ; +LIBRARY := SHI: LIBRARY TRACKS :: EDIT 0 NOLJ ; +TRACK := SYBILT: MOSTELA :: + TITLE 'TCWU04: MOSTELLER BENCHMARK WITH BURNUP' + EDIT 1 MAXR 4 ; +SYS := ASM: LIBRARY TRACK ; +FLUX := FLU: SYS LIBRARY TRACK :: + TYPE K B1 PNL BUCK 0.2948E-2 ; +EDITION := EDI: LIBRARY TRACK FLUX :: + EDIT 3 MICR RES 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 PNL 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 :: EDIT 0 NOLJ ; + SYS := DELETE: SYS ; + SYS := ASM: LIBRARY TRACK ; + FLUX := FLU: FLUX SYS LIBRARY TRACK :: + TYPE K B1 PNL 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 step6 step7 := + step3 step4 step5 step6 step7 step2 ; + EVALUATE istep := istep 1 + ; +ENDWHILE ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.7332095 ; + +res := DATABASE ; + +ISOT := DATABASE :: STEP UP default + STEP UP MIXTURES STEP AT 1 + STEP UP CALCULATIONS STEP AT 2 + STEP UP ISOTOPESLIST STEP AT 1 ; +assertS ISOT :: 'NWT0' 1 3.838717E+01 ; +assertS ISOT :: 'NWT0' 2 1.695042E+01 ; + +ECHO "test TCWU04 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU05.c2m b/Dragon/data/twlup_proc/TCWU05.c2m new file mode 100644 index 0000000..2c5d5fa --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU05.c2m @@ -0,0 +1,192 @@ +*---- +* TEST CASE TCWU05 +* CANDU-6 ANNULAR CELL +* iaea WLUP Library +* 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 file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<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 'TCWU05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + EDIT 0 MAXR 13 + ALLG BATCH 100 + TRAK TISO 5 10.0 SYMM 12 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; +VOLMATF INTLINF := EXCELT: CANDU6F :: + TITLE 'TCWU05: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + EDIT 0 MAXR 31 + ALLG BATCH 100 + 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.121034 ; +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 ; +assertS FLUX :: 'K-INFINITY' 1 0.9539985 ; +*---- +* 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.724794E0 ; +assertS ISOT :: 'NWT0' 2 2.582013E1 ; +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.327823E-2 ; +assertS ISOT :: 'NWT0' 2 3.186795E-1 ; + +INTLINF INTLINS := DELETE: INTLINF INTLINS ; +ECHO "test TCWU05 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU05Lib.c2m b/Dragon/data/twlup_proc/TCWU05Lib.c2m new file mode 100644 index 0000000..7988781 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU05Lib.c2m @@ -0,0 +1,62 @@ +*---- +* Procedure TCWU05Lib +* Create Library for test CASE TCWU05 +* Calling : +* LIBRARY := TCWU05Lib :: iedit ; +* with : +* LIBRARY = Linked list containing the result of LIB: for +* TCWU05 +* 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 file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + EDIT <<iedit>> + NMIX 17 CTRA WIMS + DEPL LIB: WIMSD4 FIL: iaea + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 560.66 0.81212 O16 = '6016' 7.99449E-1 + D2D2O = '3002' 1.99768E-1 H1H2O = '3001' 7.83774E-4 + MIX 2 560.66 6.57 Nb93 = '93' 2.50000 + BNat = '1011' 2.10000E-4 + Zr91 = '91' 9.75000E+1 + MIX 3 345.66 0.0014 He4 = '4' 1.00000E+2 + MIX 4 345.66 6.44 Fe56 = '2056' 1.60000E-1 + Ni58 = '58' 6.00000E-2 Cr52 = '52' 1.10000E-1 + BNat = '1011' 3.10000E-4 + Zr91 = '91' 9.97100E+1 + MIX 5 345.66 1.082885 O16 = '6016' 7.98895E-1 + D2D2O = '3002' 2.01016E-1 H1H2O = '3001' 8.96000E-5 + MIX 6 941.29 10.4375010 O16 = '6016' 1.18473E+1 + Xe135 = '4135' 0.0 + U235 = '2235' 6.27118E-1 1 + U238 = '8238' 8.75256E+1 1 + U236 = '236' 0.0 1 + Pu239 = '6239' 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 = '2056' 1.60000E-1 + Ni58 = '58' 6.00000E-2 Cr52 = '52' 1.10000E-1 + BNat = '1011' 3.10000E-4 + Zr91 = '91' 9.97100E+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/twlup_proc/TCWU06.c2m b/Dragon/data/twlup_proc/TCWU06.c2m new file mode 100644 index 0000000..23f5195 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU06.c2m @@ -0,0 +1,155 @@ +*---- +* TEST CASE TCWU06 +* CANDU-6 CARTESIAN CELL +* iaea WLUP Library +* 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 file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<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 'TCWU06: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 14 + ALLG BATCH 100 + 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 'TCWU06: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 32 + ALLG BATCH 100 + 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.120661 ; +*---- +* 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: WIMSD4 FIL: iaea + MIX 4 345.66 Fe56 = '2056' 6.19027E-2 + Cr52 = '52' 1.56659E-3 Ni58 = '58' 6.83337E-3 + Si29 = '29' 7.79072E-4 C12 = '2012' 1.46552E-4 + Mn55 = '55' 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.127695 ; +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.042922 ; +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 TCWU06 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU07.c2m b/Dragon/data/twlup_proc/TCWU07.c2m new file mode 100644 index 0000000..9f86d8a --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU07.c2m @@ -0,0 +1,142 @@ +*---- +* TEST CASE TCWU07 +* CANDU-6 CARTESIAN CELL +* iaea WLUP Library +* 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 file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<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 'TCWU07: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 14 + ALLG BATCH 100 + 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 'TCWU07: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 32 ANIS 2 + ALLG BATCH 100 + 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.120623 ; +EDITION := EDI: LIBRARY TRACK FLUX :: + 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.112290 ; +EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + 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.112264 ; +EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + 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.112270 ; +EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + 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 'TCWU07: 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 'TCWU07: 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.139288 ; +EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + 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.131289 ; +EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + 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.131540 ; +EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + 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.131549 ; +EDITION := EDI: EDITION LIBRARY TRACK FLUX :: + 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 TCWU07 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU08.c2m b/Dragon/data/twlup_proc/TCWU08.c2m new file mode 100644 index 0000000..c643461 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU08.c2m @@ -0,0 +1,121 @@ +*---- +* TEST CASE TCWW08 +* HOMOGENEOUS DEPLETION CASE +* iaea WLUP Library +* +* 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 file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 1 CTRA WIMS + DEPL LIB: WIMSD4 FIL: iaea + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 300.0 H1 = '3001' 2.00000E+1 + U235 = '2235' 1.0 1 + U236 = '8238' 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 'TCWW08: 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.871363 ; +EDITION := EDI: LIBRARY TRACK FLUX :: + 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 LIBRARY TRACK FLUX :: + 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.184694 ; +ECHO "test TCWU08 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU09.c2m b/Dragon/data/twlup_proc/TCWU09.c2m new file mode 100644 index 0000000..2d9d86f --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU09.c2m @@ -0,0 +1,72 @@ +*---- +* TEST CASE TCWU09 +* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL +* iaea WLUP Library +* 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 section from file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 600.0 O16 = '6016' 4.61309E-2 + U235 = '2235' 1.66078E-4 1 + U238 = '8238' 2.28994E-2 1 + MIX 2 600.0 + Zr91 = '91' 3.83243E-2 + MIX 3 600.0 + H1H2O = '3001' 4.42326E-2 O16H2O = '6016' 2.21163E-2 + BNat = '1011' 1.02133E-5 + ; +*---- +* 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 'TCWU09: 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.8276187 ; +OUT := EDI: LIBRARY VOLMAT FLUX :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ; +PIJ VOLMAT := DELETE: PIJ VOLMAT ; +VOLMAT := SYBILT: MOSTELV :: + TITLE 'TCWU09: 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.023486E-02 ; +OUT := EDI: OUT LIBRARY VOLMAT FLUX :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 SAVE ; +OUT FLUX PIJ LIBRARY VOLMAT := DELETE: + OUT FLUX PIJ LIBRARY VOLMAT ; +ECHO "test TCWU09 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU10.c2m b/Dragon/data/twlup_proc/TCWU10.c2m new file mode 100644 index 0000000..2635553 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU10.c2m @@ -0,0 +1,78 @@ +*---- +* TEST CASE TCWU10 +* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL +* iaea WLUP Library +* 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 section from file iaea format WIMSD4 +* Fixed source of 1.0E5 in group 6 +*---- +LIBRARY := LIB: :: + EDIT 0 NMIX 3 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 600.0 O16 = '6016' 4.61309E-2 + U235 = '2235' 1.66078E-4 1 + U238 = '8238' 2.28994E-2 1 + MIX 2 600.0 + Zr91 = '91' 3.83243E-2 + MIX 3 600.0 + H1H2O = '3001' 4.42326E-2 O16H2O = '6016' 2.21163E-2 + BNat = '1011' 1.02133E-5 + ; +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 + ; +*---- +* 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 'TCWU10: 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.8276187 ; +OUT := EDI: LIBRARY VOLMAT FLUX :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 ; +FLUX := DELETE: FLUX ; +FLUX := FLU: PIJ LIBRARY VOLMAT :: EDIT 99 + TYPE S EXTE 30 UNKT 1.0E-3 ; +assertV FLUX :: 'FLUX' (*GROUP*) 10 (*REGION*) 4 4.291752E+05 ; +OUT := EDI: OUT LIBRARY VOLMAT FLUX :: + EDIT 4 MERG MIX 1 2 3 COND 4.0 ; +OUT FLUX PIJ LIBRARY VOLMAT := DELETE: + OUT FLUX PIJ LIBRARY VOLMAT ; +ECHO "test TCWU10 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU11.c2m b/Dragon/data/twlup_proc/TCWU11.c2m new file mode 100644 index 0000000..ed0ea19 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU11.c2m @@ -0,0 +1,171 @@ +*---- +* TEST CASE TCWU11 +* CANDU-6 ANNULAR CELL +* iaea WLUP Library +* 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 file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<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 'TCWU11: FEW GROUP BURNUP / SELF-SHIELDING TRACKING' + EDIT 0 MAXR 13 + ALLG BATCH 100 + TRAK TISO 5 10.0 SYMM 12 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: EDIT 0 NOLJ ; +VOLMATF INTLINF := EXCELT: CANDU6F :: + TITLE 'TCWU11: FEW GROUP BURNUP / TRANSPORT TRACKING' + EDIT 0 MAXR 31 + ALLG BATCH 100 + 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.121035 ; +EDITION := EDI: LIBRARY VOLMATF FLUX 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 ; +*---- +* 69 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.9414081 ; +EDITION := EDI: EDITION LIBRARY VOLMATF FLUX 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.121052 ; +EDITION := EDI: LIBRARY VOLMATF FLUX :: + 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 LIBRARY VOLMATF FLUX :: + 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.075327 ; +res := EDITION ; +EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE: + EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ; +ECHO "test TCWU11 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU12.c2m b/Dragon/data/twlup_proc/TCWU12.c2m new file mode 100644 index 0000000..6259280 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU12.c2m @@ -0,0 +1,139 @@ +*---- +* TEST CASE TCWU12 +* iaea WLUP Library +* 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 941.29 0.72 10.437501 ; +REAL + WH1C WD2C WO16C DensCool Pres + WU235F WU238F WO16F ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY ISOT ; +MODULE + LIB: INFO: END: ; +PROCEDURE assertS ; +*---- +* Get Coolant properties +*---- +ECHO + "Case 1. Coolant density vs temperature " ; +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: WIMSD4 FIL: iaea + ISO: 3 '3001' '3002' '6016' + 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 ; +ECHO + "Case 2. Coolant density vs temperature and pressure " ; +EVALUATE + TempCool Purity Pres := + 366.72 100.0 20.8 ; +ECHO + "Input - Coolant temperature (K) " TempCool ; +ECHO + "Input - Coolant pressure (MPa) " Pres ; +ECHO + "Input - D2/(D2+H1) Weight % ratio in Coolant" Purity ; +INFO: :: + TMP: <<TempCool>> C + PRES: <<Pres>> MPa + PUR: <<Purity>> WGT% + CALC DENS PWATER >>DensCool<< + LIB: WIMSD4 FIL: iaea + ISO: 3 '3001' '3002' '6016' + 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 ; +ECHO + "Case 3. Fuel enrichment " ; +*---- +* 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: WIMSD4 FIL: iaea + ISO: 3 '2235' '8238' '6016' + 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 file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 5 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 <<TempCool>> <<DensCool>> O16 = '6016' <<WO16C>> + D2D2O = '3002' <<WD2C>> H1H2O = '3001' <<WH1C>> + MIX 2 <<TempFuel>> <<DensCool>> O16 = '6016' <<WO16C>> + D2D2O = '3002' <<WD2C>> H1H2O = '3001' <<WH1C>> + MIX 3 <<TempFuel>> <<DensFuel>> O16 = '6016' <<WO16F>> + U235 = '2235' <<WU235F>> 1 + U238 = '8238' <<WU238F>> 1 + MIX 4 COMB 2 0.5 0 0.5 + MIX 5 COMB 3 0.1 3 0.9 + ; +ISOT := LIBRARY :: STEP UP ISOTOPESLIST STEP AT 15 ; (*U238 0005*) +assertS ISOT :: NG 35 5.000504E-01 ; +ECHO "test TCWU12 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU13.c2m b/Dragon/data/twlup_proc/TCWU13.c2m new file mode 100644 index 0000000..7fc5388 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU13.c2m @@ -0,0 +1,126 @@ +*---- +* TEST CASE TCWU13 +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* WIMSD4 69 GROUPS LIBRARY FILE iaea FROM WLUP +* +* REF: none +* +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + ASSMB DISCR1 DISCR2 LIBRARY CP CALC OUT DATABASE ISOT SPHGEOM + MTRACK ; +SEQ_BINARY FILTRK ; +SEQ_ASCII + res ; +MODULE + LIB: GEO: SYBILT: NXT: BIVACT: SHI: ASM: FLU: EDI: COMPO: SPH: + DELETE: END: ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from file iaea format WIMSD4 +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 579.9 + H1H2O = '3001' 4.76690E-2 O16H2O = '6016' 2.38345E-2 + BNat = '1011' 2.38103E-5 + MIX 2 579.9 O16 = '6016' 3.06711E-4 + Cr52 = '52' 7.54987E-5 Fe56 = '2056' 1.47624E-4 + Zr91 = '91' 4.18621E-2 + MIX 3 579.9 + H1H2O = '3001' 4.65292E-2 O16H2O = '6016' 2.32646E-2 + Cr52 = '52' 4.79927E-5 + Fe56 = '2056' 4.45845E-5 Ni58 = '58' 1.13521E-4 + MoNat = '96' 4.03755E-6 + Al27 = '27' 2.35231E-6 + Mn55 = '55' 4.15901E-7 + BNat = '1011' 2.32761E-5 + Zr91 = '91' 8.92427E-4 + MIX 4 579.9 O16 = '6016' 2.87335E-4 + Cr52 = '52' 7.07291E-5 Fe56 = '2056' 1.38298E-4 + Zr91 = '91' 3.92175E-2 + MIX 5 579.9 + H1H2O = '3001' 4.71346E-2 O16H2O = '6016' 2.35673E-2 + Cr52 = '52' 2.24991E-5 + Fe56 = '2056' 2.09013E-5 Ni58 = '58' 5.32188E-5 + MoNat = '96' 1.89281E-6 + Al27 = '27' 1.10277E-6 + Mn55 = '55' 1.94976E-7 + BNat = '1011' 2.35598E-5 + Zr91 = '91' 4.18372E-4 + MIX 6 579.9 + H1H2O = '3001' 4.71676E-2 O16H2O = '6016' 2.35838E-2 + Cr52 = '52' 2.11122E-5 + Fe56 = '2056' 1.96130E-5 Ni58 = '58' 4.99383E-5 + MoNat = '96' 1.77614E-6 + Al27 = '27' 1.03479E-6 + Mn55 = '55' 1.82957E-7 + BNat = '1011' 2.35753E-5 + Zr91 = '91' 3.92583E-4 + MIX 7 579.9 + H1H2O = '3001' 4.72020E-2 O16H2O = '6016' 2.36010E-2 + Cr52 = '52' 1.96591E-5 + Fe56 = '2056' 1.82630E-5 Ni58 = '58' 4.65011E-5 + MoNat = '96' 1.65389E-6 + Al27 = '27' 9.63569E-7 + Mn55 = '55' 1.70365E-7 + BNat = '1011' 2.35914E-5 + Zr91 = '91' 3.65562E-4 + MIX 8 933.6 O16 = '6016' 4.49355E-2 + U235 = '2235' 7.39237E-4 1 + U238 = '8238' 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 + ::: 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 +* Editing using SPH model for transport-diffusion +*---- +DISCR1 := SYBILT: ASSMB :: + TITLE 'TCWU13: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; +DISCR2 FILTRK := NXT: ASSMB :: + TITLE 'TCWU13: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + ALLG BATCH 100 + TISO 10 20.0 ; +LIBRARY := SHI: LIBRARY DISCR1 :: EDIT 0 NOLJ ; + +CP := ASM: LIBRARY DISCR2 FILTRK ; +CALC := FLU: CP LIBRARY DISCR2 :: + TYPE B B1 PNL ; +assertS CALC :: 'K-INFINITY' 1 1.256567 ; + +OUT := EDI: LIBRARY DISCR1 CALC ASSMB :: + EDIT 3 UPS SAVE MICR RES MERGE CELL COND 4.0 + ; +ECHO "test TCWU13 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU14.c2m b/Dragon/data/twlup_proc/TCWU14.c2m new file mode 100644 index 0000000..3f9e542 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU14.c2m @@ -0,0 +1,121 @@ +*---- +* TEST CASE TCWU14 +* CANDU-6 ANNULAR CELL +* iaea WLUP Library +*---- +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX EDITION + DATABASE ISOT ; +SEQ_BINARY + INTLINS INTLINF ; +SEQ_ASCII + database ; +MODULE + GEO: NXT: SHI: ASM: FLU: EDI: DELETE: END: ; +PROCEDURE assertS ; +*---- +* Depletion data from file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<iedit>> ; +*---- +* Geometry CANDU6S : 13 regions annular cluster for self-shielding +* CANDU6F : 31 regions annular cluster for transport +*---- +CANDU6S := GEO: :: CARCEL 5 + X+ REFL X- REFL + Y+ REFL Y- REFL + MESHX -14.2875 14.2875 + 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: :: CAR2D 2 2 + X+ REFL X- REFL + Y+ REFL Y- REFL + CELL F1 F2 F3 F4 + ::: F1 := GEO: CARCEL 5 + MESHX -14.2875 14.2875 + 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 ; + ; + ::: F2 := GEO: CARCEL 5 + MESHX -14.2875 14.2875 + 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 ; + ; + ::: F3 := GEO: CARCEL 5 + MESHX -14.2875 14.2875 + 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 ; + ; + ::: F4 := GEO: CARCEL 5 + MESHX -14.2875 14.2875 + 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 ; + ; +; +*---- +* Self-Shielding calculation NXT +* Transport calculation NXT +* Flux calculation for keff +*---- +VOLMATS INTLINS := NXT: CANDU6S :: + TITLE 'TCWU14: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + ALLG BATCH 100 MERGMIX + EDIT 0 TISO 50 10.0 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: + EDIT 0 ; +VOLMATF INTLINF := NXT: CANDU6F :: + TITLE 'TCWU14: CANDU-6 ANNULAR POWER= 31.971 FUEL TEMP= 941.29' + ALLG BATCH 100 + EDIT 0 TISO 50 10.0 ; +PIJ := ASM: LIBRARY VOLMATF INTLINF ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-INFINITY' 1 1.129156 ; +EDITION := EDI: LIBRARY VOLMATF FLUX :: + COND 4.0 MERGE CELL MICR 1 Xe135 SAVE ON 'fuel' ; +INTLINF INTLINS := DELETE: INTLINF INTLINS ; +ECHO "test TCWU14 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU15.c2m b/Dragon/data/twlup_proc/TCWU15.c2m new file mode 100644 index 0000000..cf34a09 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU15.c2m @@ -0,0 +1,157 @@ +*---- +* TEST CASE TCWU05 +* CANDU-6 ANNULAR CELL +* iaea WLUP Library +* 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 + COMPO1 COMPO2 ; +SEQ_BINARY + INTLINS INTLINF ; +SEQ_ASCII + fuel mode ; +MODULE + GEO: NXT: SHI: ASM: FLU: EVO: EDI: CPO: DELETE: END: ; +*---- +* Depletion data from file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +PROCEDURE assertS ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<iedit>> ; +*---- +* Geometry CANDU6S : 13 regions annular cluster for self-shielding +* CANDU6F : 31 regions annular cluster for transport +*---- +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 +*---- +VOLMATS INTLINS := NXT: CANDU6S :: + TITLE 'TCWU05: CANDU-6 CARTESIAN POWER= 31.971 FUEL TEMP= 941.29' + EDIT 0 TRAK TISO 5 10.0 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: + EDIT 0 ; +VOLMATF INTLINF := NXT: CANDU6F :: + TITLE 'TCWU05: CANDU-6 CARTESIAN POWER= 31.971 FUEL TEMP= 941.29' + ALLG BATCH 100 + EDIT 0 TRAK TISO 5 10.0 ; +PIJ := ASM: LIBRARY VOLMATF INTLINF :: ; +FLUX := FLU: PIJ LIBRARY VOLMATF :: + TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.120613 ; +EDITION := EDI: FLUX LIBRARY VOLMATF :: + COND 4.0 MERGE MIX 0 0 0 0 1 0 0 0 0 0 SAVE ON 'EDITMOD' ; +EDITION := EDI: EDITION FLUX LIBRARY VOLMATF :: + COND 4.0 MERGE COMP MICR 1 Xe135 SAVE ; +*---- +* 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 ; + 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 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 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9537331 ; +*---- +* Save calculation results in CPO format file +*---- +COMPO1 := CPO: BURNUP EDITION :: + BURNUP REF-CASE EXTRACT Xe135 Xe135 NAME MIXTRXE ; +fuel := COMPO1 ; +COMPO2 := CPO: EDITION :: + STEP 'EDITMOD' NAME MIXTMOD ; +mode := COMPO2 ; +INTLINF INTLINS := DELETE: INTLINF INTLINS ; +ECHO "test TCWU15 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU17.c2m b/Dragon/data/twlup_proc/TCWU17.c2m new file mode 100644 index 0000000..4dd7de0 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU17.c2m @@ -0,0 +1,337 @@ +*---- +* Exemple of the use of HMIX for cell homogenization +* 2-D supercell with fuel clusters based on AECL supercell model +* for G2 SOR and MCA with fuel and reactivity devices parallel +* References +* PREPARED BY : G. Marleau on 2013/06/11 +* +*---- +* modules and data structures +*---- +SEQ_ASCII MACROLIBF ; +LINKED_LIST SORINS SORIN TRACK MicLib FLUX EDITION ; +XSM_FILE ASMPIJ ; +SEQ_ASCII FigReg.ps FigMix.ps FigHom.ps ; +SEQ_ASCII HomMix.txt HomHMix.txt ; +SEQ_BINARY Lines ; +MODULE GEO: EXCELT: EXCELL: SHI: ASM: LIB: FLU: + EDI: UTL: DELETE: FREE: END: INFO: CPO: MAC: + NXT: PSP: ; +PROCEDURE TCWU17Lib ; +PROCEDURE assertS ; +INTEGER iedit := 1 ; +MicLib := TCWU17Lib :: <<iedit>> ; +*---- +* DEFINE GEOMETRY FOR SUPERCELL CALCULATION +* SORINS : 2D self-shielding geometry with SHUT-OFF ROD & GT in +* for annular fuel. +* SORIN : 2D transport geometry with SHUT-OFF ROD & GT in +* for annular fuel. +*---- +SORINS := GEO: :: CAR2D 5 3 + EDIT 0 + X- REFL X+ REFL + Y- REFL Y+ REFL + CELL M MXL MX2 MXR M + MY FXYL BXY FXYR MY + M MXL MX2 MXR M + ::: M := GEO: CAR2D 1 1 + MESHX 0.0 7.14375 + MESHY 0.0 7.14375 + HMIX 0 + MIX 5 ; + ::: MXL := GEO: CAR2D 2 1 + MESHX -7.14375 0.0 7.14375 + MESHY 0.0 7.14375 + HMIX 0 1 + MIX 5 15 ; + ::: MX2 := GEO: CAR2D 1 1 + MESHX -7.14375 7.14375 + MESHY 0.0 7.14375 + HMIX 1 + MIX 15 ; + ::: MXR := GEO: CAR2D 2 1 + MESHX -7.14375 0.0 7.14375 + MESHY 0.0 7.14375 + HMIX 1 0 + MIX 15 5 ; + ::: MY := GEO: CAR2D 1 2 + MESHY -7.14375 0.0 7.14375 + MESHX 0.0 7.14375 + HMIX 0 0 + MIX 5 5 ; + ::: BXY := GEO: CARCEL 2 + MESHX -7.14375 7.14375 + MESHY -7.14375 7.14375 + RADIUS 0.0 6.380 6.530 + HMIX 1 1 1 + MIX 15 15 15 + CLUSTER ROD + ::: ROD := GEO: TUBE 4 + NPIN 1 RPIN 0.0 APIN 0.0 + RADIUS 0.0 5.4115 5.4877 5.5791 5.6553 + HMIX 1 1 1 1 + MIX 15 14 21 14 ; ; + ::: FXYL := GEO: CARCEL 5 2 1 + MESHX -7.14375 0.0 7.14375 + MESHY -7.14375 7.14375 + RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00 + HMIX 0 0 0 0 0 0 + 1 1 1 1 1 1 + MIX 1 2 3 4 5 5 + 11 12 13 14 15 15 + CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R + ::: ROD1 := GEO: TUBE 2 1 2 MIX 6 10 16 20 HMIX 0 0 1 1 + NPIN 1 RPIN 0.0000 APIN 0.0000 + MESHX -0.6540 0.6540 + MESHY -0.6540 0.0 0.6540 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2L := GEO: TUBE 2 MIX 7 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN 2.09439510 3.14159265 4.18879020 ; + ::: ROD2R := GEO: TUBE 2 MIX 17 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN -1.04719755 0.0000 1.04719755 ; + ::: ROD3L := GEO: TUBE 2 MIX 8 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -2.87979327 -2.35619449 -1.83259571 + 1.83259571 2.35619449 2.87979327 ; + ::: ROD3R := GEO: TUBE 2 MIX 18 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -1.30899694 -0.78539816 -0.26179939 + 0.26179939 0.78539816 1.30899694 ; + ::: ROD4L := GEO: TUBE 2 MIX 9 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN 1.74532925 2.09439510 2.44346095 + 2.79252680 3.14159265 3.49065850 + 3.83972435 4.18879020 4.53785606 ; + ::: ROD4R := GEO: TUBE 2 MIX 19 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN -1.39626340 -1.04719755 -0.69813170 + -0.34906585 0.0 0.34906585 + 0.69813170 1.04719755 1.39626340 ; + ; + ::: FXYR := GEO: CARCEL 5 2 1 + MESHX -7.14375 0.0 7.14375 + MESHY -7.14375 7.14375 + RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00 + HMIX 1 1 1 1 1 1 + 0 0 0 0 0 0 + MIX 11 12 13 14 15 15 + 1 2 3 4 5 5 + CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R + ::: ROD1 := GEO: TUBE 2 1 2 MIX 16 20 6 10 HMIX 1 1 0 0 + NPIN 1 RPIN 0.0000 APIN 0.0000 + MESHX -0.6540 0.6540 + MESHY -0.6540 0.0 0.6540 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2L := GEO: TUBE 2 MIX 17 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN 2.09439510 3.14159265 4.18879020 ; + ::: ROD2R := GEO: TUBE 2 MIX 7 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN -1.04719755 0.0000 1.04719755 ; + ::: ROD3L := GEO: TUBE 2 MIX 18 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -2.87979327 -2.35619449 -1.83259571 + 1.83259571 2.35619449 2.87979327 ; + ::: ROD3R := GEO: TUBE 2 MIX 8 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -1.30899694 -0.78539816 -0.26179939 + 0.26179939 0.78539816 1.30899694 ; + ::: ROD4L := GEO: TUBE 2 MIX 19 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN 1.74532925 2.09439510 2.44346095 + 2.79252680 3.14159265 3.49065850 + 3.83972435 4.18879020 4.53785606 ; + ::: ROD4R := GEO: TUBE 2 MIX 9 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN -1.39626340 -1.04719755 -0.69813170 + -0.34906585 0.0 0.34906585 + 0.69813170 1.04719755 1.39626340 ; + ; +; +SORIN := GEO: :: CAR2D 5 3 + EDIT 0 + X- REFL X+ REFL + Y- REFL Y+ REFL + CELL M MXL MX2 MXR M + MY FXYL BXY FXYR MY + M MXL MX2 MXR M + ::: M := GEO: CAR2D 1 1 + MESHX 0.0 7.14375 SPLITX 2 + MESHY 0.0 7.14375 SPLITY 2 + HMIX 0 + MIX 5 ; + ::: MXL := GEO: CAR2D 2 1 + MESHX -7.14375 0.0 7.14375 SPLITX 3 3 + MESHY 0.0 7.14375 SPLITY 3 + HMIX 0 1 + MIX 5 15 ; + ::: MX2 := GEO: CAR2D 1 1 + MESHX -7.14375 7.14375 SPLITX 6 + MESHY 0.0 7.14375 SPLITY 3 + HMIX 1 + MIX 15 ; + ::: MXR := GEO: CAR2D 2 1 + MESHX -7.14375 0.0 7.14375 SPLITX 3 3 + MESHY 0.0 7.14375 SPLITY 3 + HMIX 1 0 + MIX 15 5 ; + ::: MY := GEO: CAR2D 1 2 + MESHY -7.14375 0.0 7.14375 SPLITY 3 3 + MESHX 0.0 7.14375 SPLITX 3 + HMIX 0 0 + MIX 5 5 ; + ::: BXY := GEO: CARCEL 2 + MESHX -7.14375 7.14375 SPLITX 6 + MESHY -7.14375 7.14375 SPLITY 6 + RADIUS 0.0 6.380 6.530 + SPLITR 2 2 + HMIX 1 1 1 + MIX 15 15 15 + CLUSTER ROD + ::: ROD := GEO: TUBE 4 + NPIN 1 RPIN 0.0 APIN 0.0 + RADIUS 0.0 5.4115 5.4877 5.5791 5.6553 + SPLITR 1 1 1 1 + HMIX 1 1 1 1 + MIX 15 14 21 14 ; ; + ::: FXYL := GEO: CARCEL 5 2 1 + MESHX -7.14375 0.0 7.14375 SPLITX 3 3 + MESHY -7.14375 7.14375 SPLITY 6 + RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00 + HMIX 0 0 0 0 0 0 + 1 1 1 1 1 1 + MIX 1 2 3 4 5 5 + 11 12 13 14 15 15 + CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R + ::: ROD1 := GEO: TUBE 2 1 2 MIX 6 10 16 20 HMIX 0 0 1 1 + NPIN 1 RPIN 0.0000 APIN 0.0000 + MESHX -0.6540 0.6540 + MESHY -0.6540 0.0 0.6540 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2L := GEO: TUBE 2 MIX 7 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN 2.09439510 3.14159265 4.18879020 ; + ::: ROD2R := GEO: TUBE 2 MIX 17 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN -1.04719755 0.0000 1.04719755 ; + ::: ROD3L := GEO: TUBE 2 MIX 8 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -2.87979327 -2.35619449 -1.83259571 + 1.83259571 2.35619449 2.87979327 ; + ::: ROD3R := GEO: TUBE 2 MIX 18 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -1.30899694 -0.78539816 -0.26179939 + 0.26179939 0.78539816 1.30899694 ; + ::: ROD4L := GEO: TUBE 2 MIX 9 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN 1.74532925 2.09439510 2.44346095 + 2.79252680 3.14159265 3.49065850 + 3.83972435 4.18879020 4.53785606 ; + ::: ROD4R := GEO: TUBE 2 MIX 19 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN -1.39626340 -1.04719755 -0.69813170 + -0.34906585 0.0 0.34906585 + 0.69813170 1.04719755 1.39626340 ; + ; + ::: FXYR := GEO: CARCEL 5 2 1 + MESHX -7.14375 0.0 7.14375 SPLITX 3 3 + MESHY -7.14375 7.14375 SPLITY 6 + RADIUS 0.00000 5.16890 5.60320 6.44780 6.58750 7.00 + HMIX 1 1 1 1 1 1 + 0 0 0 0 0 0 + MIX 11 12 13 14 15 15 + 1 2 3 4 5 5 + CLUSTER ROD1 ROD2L ROD2R ROD3L ROD3R ROD4L ROD4R + ::: ROD1 := GEO: TUBE 2 1 2 MIX 16 20 6 10 HMIX 1 1 0 0 + NPIN 1 RPIN 0.0000 APIN 0.0000 + MESHX -0.6540 0.6540 + MESHY -0.6540 0.0 0.6540 + RADIUS 0.00000 0.6122 0.6540 ; + ::: ROD2L := GEO: TUBE 2 MIX 17 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN 2.09439510 3.14159265 4.18879020 ; + ::: ROD2R := GEO: TUBE 2 MIX 7 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 3 RPIN 1.4885 + APIN -1.04719755 0.0000 1.04719755 ; + ::: ROD3L := GEO: TUBE 2 MIX 18 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -2.87979327 -2.35619449 -1.83259571 + 1.83259571 2.35619449 2.87979327 ; + ::: ROD3R := GEO: TUBE 2 MIX 8 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 6 RPIN 2.8755 + APIN -1.30899694 -0.78539816 -0.26179939 + 0.26179939 0.78539816 1.30899694 ; + ::: ROD4L := GEO: TUBE 2 MIX 19 20 HMIX 1 1 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN 1.74532925 2.09439510 2.44346095 + 2.79252680 3.14159265 3.49065850 + 3.83972435 4.18879020 4.53785606 ; + ::: ROD4R := GEO: TUBE 2 MIX 9 10 HMIX 0 0 + RADIUS 0.00000 0.6122 0.6540 + NPIN 9 RPIN 4.3305 + APIN -1.39626340 -1.04719755 -0.69813170 + -0.34906585 0.0 0.34906585 + 0.69813170 1.04719755 1.39626340 ; + ; +; +*---- +* Rod and GT absent +*---- +TRACK Lines := NXT: SORINS :: EDIT 3 + ALLG BATCH 100 MERGMIX + TISO 40 30.0 ; +MicLib := SHI: MicLib TRACK Lines ; +TRACK Lines := DELETE: TRACK Lines ; +TRACK Lines := NXT: SORIN :: EDIT 3 + ALLG BATCH 100 + TISO 40 30.0 ; +FigReg.ps := PSP: TRACK :: TYPE REGI ; +FigMix.ps := PSP: TRACK :: TYPE MIXT ; +FigHom.ps := PSP: TRACK :: TYPE HMIX ; +ASMPIJ := ASM: MicLib TRACK Lines :: PIJ ; +FLUX := FLU: ASMPIJ MicLib TRACK :: TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.036408 ; +EDITION := EDI: FLUX MicLib TRACK :: + EDIT 3 + COND 0.625 + MERG HMIX + SAVE ON SORINHMIX ; +HomHMix.txt := EDITION ; +EDITION := DELETE: EDITION ; +EDITION := EDI: FLUX MicLib TRACK :: + EDIT 3 + COND 0.625 + MERG MIX 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 + SAVE ON SORINMIX ; +HomMix.txt := EDITION ; +EDITION := DELETE: EDITION ; +TRACK FLUX ASMPIJ Lines := DELETE: TRACK FLUX ASMPIJ Lines ; +ECHO "test TCWU17 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU17Lib.c2m b/Dragon/data/twlup_proc/TCWU17Lib.c2m new file mode 100644 index 0000000..730b54a --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU17Lib.c2m @@ -0,0 +1,85 @@ +*---- +* Create Library for test CASE TCWU17.x2m +* Calling : +* LIBRARY := TCWU17Lib.x2m :: iedit ; +* with : +* LIBRARY = Linked list containing the result of LIB: for +* TCWU17.x2m +* iprint = print level for LIB: module +*---- +* Define PARAMETERS,STRUCTURES and MODULES used +*---- +PARAMETER LIBRARY :: ::: LINKED_LIST LIBRARY ; ; +MODULE LIB: DELETE: END: ; +*---- +* Define and read LIB: EDIT option +INTEGER iedit ; +:: >>iedit<< ; +*---- +* Depletion data from file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +* All materials are duplicated for left and right cell identification +* during homogenization +*---- +LIBRARY := LIB: :: + EDIT <<iedit>> + NMIX 21 CTRA WIMS + DEPL LIB: WIMSD4 FIL: iaea + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 560.66 0.81212 O16 = '6016' 7.99449E-1 + D2D2O = '3002' 1.99768E-1 H1H2O = '3001' 7.83774E-4 + MIX 2 560.66 6.57 Nb93 = '93' 2.50000 + BNat = '1011' 2.10000E-4 + Zr91 = '91' 9.75000E+1 + MIX 3 345.66 0.0014 He4 = '4' 1.00000E+2 + MIX 4 345.66 6.44 Fe56 = '2056' 1.60000E-1 + Ni58 = '58' 6.00000E-2 Cr52 = '52' 1.10000E-1 + BNat = '1011' 3.10000E-4 + Zr91 = '91' 9.97100E+1 + MIX 5 345.66 1.082885 O16 = '6016' 7.98895E-1 + D2D2O = '3002' 2.01016E-1 H1H2O = '3001' 8.96000E-5 + MIX 6 941.29 10.4375010 O16 = '6016' 1.18473E+1 + Xe135 = '4135' 0.0 + U235 = '2235' 6.27118E-1 1 + U238 = '8238' 8.75256E+1 1 + U236 = '236' 0.0 1 + Pu239 = '6239' 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 = '2056' 1.60000E-1 + Ni58 = '58' 6.00000E-2 Cr52 = '52' 1.10000E-1 + BNat = '1011' 3.10000E-4 + Zr91 = '91' 9.97100E+1 + MIX 11 560.66 0.81212 O16 = '6016' 7.99449E-1 + D2D2O = '3002' 1.99768E-1 H1H2O = '3001' 7.83774E-4 + MIX 12 560.66 6.57 Nb93 = '93' 2.50000 + BNat = '1011' 2.10000E-4 + Zr91 = '91' 9.75000E+1 + MIX 13 345.66 0.0014 He4 = '4' 1.00000E+2 + MIX 14 345.66 6.44 Fe56 = '2056' 1.60000E-1 + Ni58 = '58' 6.00000E-2 Cr52 = '52' 1.10000E-1 + BNat = '1011' 3.10000E-4 + Zr91 = '91' 9.97100E+1 + MIX 15 345.66 1.082885 O16 = '6016' 7.98895E-1 + D2D2O = '3002' 2.01016E-1 H1H2O = '3001' 8.96000E-5 + MIX 16 941.29 10.4375010 O16 = '6016' 1.18473E+1 + Xe135 = '4135' 0.0 + U235 = '2235' 6.27118E-1 1 + U238 = '8238' 8.75256E+1 1 + U236 = '236' 0.0 1 + Pu239 = '6239' 0.0 1 + MIX 17 COMB 6 1.0 + MIX 18 COMB 6 1.0 + MIX 19 COMB 6 1.0 + MIX 20 560.66 6.44 Fe56 = '2056' 1.60000E-1 + Ni58 = '58' 6.00000E-2 Cr52 = '52' 1.10000E-1 + BNat = '1011' 3.10000E-4 + Zr91 = '91' 9.97100E+1 + MIX 21 345.66 Fe56 = '2056' 6.19027E-2 + Cr52 = '52' 1.56659E-3 Ni58 = '58' 6.83337E-3 + Si29 = '29' 7.79072E-4 C12 = '2012' 1.46552E-4 + Mn55 = '55' 1.25431E-3 + ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU18.c2m b/Dragon/data/twlup_proc/TCWU18.c2m new file mode 100644 index 0000000..6a015dc --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU18.c2m @@ -0,0 +1,94 @@ +*---- +* Exemple 3D NXT case with self-shielding +* +*---- +* modules and data structures +*---- +REAL VoidF := 0.0 ; +REAL DtCool := 0.653 ; +REAL TFuel := 900.0 ; +REAL TCool := 600.0 ; +REAL RFuel RClad := 0.608 0.628 ; +INTEGER IDcool := 4 ; +REAL HA1 HA2 HA3 := 20.0 10.0 50.0 ; +INTEGER NumZA1 NumZA2 NumZA3 := 2 1 10 ; +INTEGER IDA1 IDA2 IDA3 := 4 5 1 ; +EVALUATE DtCool := DtCool VoidF * ; ! Evaluate the coolant density +REAL pitch := 1.26 ; +INTEGER NumX NumRF NumRC := 1 -2 1 ; +LINKED_LIST BWRF MICLIB MICLIB2 LATTICE PIJ SYS FLUX EDITION TRKNAM ; +SEQ_ASCII res FigReg.ps micro ; +MODULE MAC: ASM: FLU: NXT: END: GEO: LIB: SHI: EDI: EXCELT: USS: ; +PROCEDURE assertS ; + +************************************************************************ +* Material Definitions * +************************************************************************ +MICLIB2 := LIB: :: + NMIX 5 CTRA WIMS SUBG +* CALENDF 4 (*CALENDF TYPE PROBABILITY TABLES*) + DEPL LIB: WIMSD4 FIL: iaea + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 <<TFuel>> 11.72 ! ThO2 11.72 g/cm3 + U233 = '9233' 8.82 1 ! UO2 10.3 g/cm3 + U235 = '2235' 0.0 1 ! UN 14.42 g/cm3 + Th232 = '2232' 79.06 1 + O16 = '6016' 12.12 + MIX 2 <<TFuel>> 6.55 ! Cladding + SN = '118' 1.45 + Zr4 = '91' 98.23 + Fe = '1054' 0.22 + Cr = '52' 0.1 + MIX 3 <<TCool>> <<DtCool>> ! H2O + HH = '3001' 20.0 + O16 = '6016' 80.0 + MIX 4 <<TCool>> <<DtCool>> ! D2O + HD = '3002' 20.0 + O16 = '6016' 80.0 + MIX 5 <<TFuel>> 11.72 ! ThO2 11.72 g/cm3 + U233 = '9233' 0.0 1 ! UO2 10.3 g/cm3 + U235 = '2235' 0.0 1 ! UN 14.42 g/cm3 + Th232 = '2232' 87.88 1 + O16 = '6016' 12.12 + ; + +************************************************************************ +* Geometry Definitions 3D Pin cell * +************************************************************************ +BWRF := GEO: :: CAR3D 1 1 5 + X- REFL X+ REFL + Y- REFL Y+ REFL + Z- REFL Z+ REFL + CELL A1 A2 A3 A2 A1 + ::: A1 := GEO: CARCELZ 2 1 + MESHX 0.0 <<pitch>> SPLITX <<NumX>> + MESHY 0.0 <<pitch>> SPLITY <<NumX>> + MESHZ 0.0 <<HA1>> SPLITZ <<NumZA1>> + RADIUS 0.0 <<RFuel>> <<RClad>> + SPLITR <<NumRF>> <<NumRC>> + MIX <<IDA1>> 2 <<IDcool>> ; + ::: A2 := GEO: A1 + MESHZ 0.0 <<HA2>> SPLITZ <<NumZA2>> + MIX <<IDA2>> 2 <<IDcool>> ; + ::: A3 := GEO: A1 + MESHZ 0.0 <<HA3>> SPLITZ <<NumZA3>> + MIX <<IDA3>> 2 <<IDcool>> ; + ; + +TRKNAM := NXT: BWRF :: + ALLG BATCH 100 + TISO 4 12.0 ; + +MICLIB := USS: MICLIB2 TRKNAM :: EDIT 2 ; + +SYS := ASM: MICLIB TRKNAM :: ; +FLUX := FLU: SYS MICLIB TRKNAM :: + TYPE K ; +EDITION := EDI: FLUX MICLIB TRKNAM :: + EDIT 0 + MERGE COMP + SAVE ; +assertS FLUX :: 'K-INFINITY' 1 1.090306 ; +ECHO "test TCWU18 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU31.c2m b/Dragon/data/twlup_proc/TCWU31.c2m new file mode 100644 index 0000000..35f8425 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU31.c2m @@ -0,0 +1,172 @@ +*---- +* TEST CASE TCWU31 +* CANDU-6 ANNULAR CELL +* iaea WLUP Library +* 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 file iaea format WIMSD4 +* Microscopic cross sections from file iaea format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<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 'TCWU31: 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 'TCWU31: 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.121035 ; +EDITION := EDI: LIBRARY VOLMATF FLUX 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 'TCWU31: 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.121044 ; +EDITION := EDI: LIBRARY VOLMATF FLUX :: + 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 LIBRARY VOLMATF FLUX :: + 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.075343 ; +res := EDITION ; +EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE: + EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ; +ECHO "test TCWU31 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU32.c2m b/Dragon/data/twlup_proc/TCWU32.c2m new file mode 100644 index 0000000..b4974bf --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU32.c2m @@ -0,0 +1,95 @@ +*---- +* TEST CASE TCWU32 +* VVER Assembly of central water hole with 6 identical fuel cells +* NXT self-shielding and flux +* Coarse mesh +*---- +MODULE + LIB: GEO: NXT: SHI: ASM: FLU: DELETE: END: ; +LINKED_LIST + LIBRARY GEOMS GEOMF TRACK PIJMAT FLUX ; +SEQ_BINARY + TRACK_LGT ; +PROCEDURE assertS ; +*---- +* Microscopic cross sections from file WNEALIB format WIMSD4 +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 4 + CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 1 579.9 + H1H2O = '3001' 4.76690E-2 O16H2O = '6016' 2.38345E-2 + BNat = '1011' 2.38103E-5 + MIX 2 609.0 O16 = '6016' 3.06711E-4 + Cr52 = '52' 7.54987E-5 Fe56 = '2056' 1.47624E-4 + Zr91 = '91' 4.18621E-2 + MIX 3 966.0 O16 = '6016' 4.49355E-2 + U235 = '2235' 7.39237E-4 1 + U238 = '8238' 2.17285E-2 1 + MIX 4 579.9 + H1H2O = '3001' 4.76690E-2 O16H2O = '6016' 2.38345E-2 + BNat = '1011' 2.38103E-5 + ; +*---- +* Geometries NXT +*---- +GEOMS := GEO: :: HEX 7 + EDIT 2 + HBC COMPLETE REFL + CELL CT CFI CFI CFI CFI CFI CFI + ::: CT := GEO: HEXTCEL 2 1 + SIDE 0.736121593 + MIX 4 2 1 + RADIUS 0.0 0.386 0.445 + ; + ::: CFI := GEO: HEXTCEL 2 1 + SIDE 0.736121593 + MIX 3 2 1 + RADIUS 0.0 0.386 0.445 + ; + ; +GEOMF := GEO: :: HEX 7 + EDIT 2 + HBC COMPLETE REFL + CELL CT CFI CFI CFI CFI CFI CFI + ::: CT := GEO: HEXTCEL 2 1 + SIDE 0.736121593 + SPLITR 3 1 + MIX 1 2 1 + RADIUS 0.0 0.386 0.445 + ; + ::: CFI := GEO: HEXTCEL 2 1 + SIDE 0.736121593 + SPLITR 3 1 + MIX 3 2 1 + RADIUS 0.0 0.386 0.445 + ; + ; +*---- +* Self-Shielding calculation +* Transport calculation +* Flux calculation +*---- +TRACK TRACK_LGT := NXT: GEOMS :: + EDIT 3 TITL 'VVER-NXT: Shielding' + ALLG BATCH 100 + TISO 50 20.0 ; +LIBRARY := SHI: LIBRARY TRACK TRACK_LGT ; +TRACK TRACK_LGT GEOMS := DELETE: TRACK TRACK_LGT GEOMS ; +TRACK TRACK_LGT := NXT: GEOMF :: + EDIT 3 TITL 'VVER-NXT: Flux' + ALLG BATCH 100 + TISO 50 100.0 ; +PIJMAT := ASM: LIBRARY TRACK TRACK_LGT :: + EDIT 1 ; +FLUX := FLU: PIJMAT LIBRARY TRACK :: + EDIT 1 TYPE K ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.266641 ; +TRACK TRACK_LGT GEOMF := DELETE: TRACK TRACK_LGT GEOMF ; +PIJMAT FLUX := DELETE: PIJMAT FLUX ; +LIBRARY := DELETE: LIBRARY ; +ECHO "test TCWU32 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU33.c2m b/Dragon/data/twlup_proc/TCWU33.c2m new file mode 100644 index 0000000..05fe0a3 --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU33.c2m @@ -0,0 +1,378 @@ +*---- +* TEST CASE TCWU33 +* CANDU cluster +*---- +*---- +* Definicion de Estructuras y Modulos usados +*---- +* +LINKED_LIST + A2S A2T TRACK_S TRACK_T SYS VectComp + Pij BURNUP Vector burnups Totburns + LIBRARY LIBRARY_S + LIBRARY_P + LIBRARY_N + EDITION0 EDITION1 EDITION2 EDITION3 + EDITION4 EDITION5 EDITION6 EDITION7 + FLUX FLUX1 FLUX2 +; + +MODULE + GEO: LIB: SHI: EXCELT: ASM: FLU: EDI: DELETE: END: EVO: PSP: + USS: UTL: GREP: COMPO: ; + +SEQ_BINARY + INTLIN_S INTLIN_T ; + +PROCEDURE Datos_pertur assertS ; + +SEQ_ASCII CEL37B2G_Bu :: FILE './CEL37B2G_Bu' ; +SEQ_ASCII QUEMADO :: FILE './QUEMADO' ; +SEQ_ASCII GRAF_A2S.PS :: FILE './GRAF_A2S.PS' ; +SEQ_ASCII GRAF_A2T.PS :: FILE './GRAF_A2T.PS' ; + +* XSref = no perturbada; XSPi = perturbada +SEQ_ASCII XS_DcoolP :: FILE './XS_DcoolP' ; +SEQ_ASCII XSref :: FILE './XSref' ; +SEQ_ASCII XS_combP :: FILE './XS_combP' ; +SEQ_ASCII XS_TrefP :: FILE './XS_TrefP' ; +SEQ_ASCII XS_DrefP :: FILE './XS_DrefP' ; +SEQ_ASCII XS_T_DrefP :: FILE './XS_T_DrefP' ; +SEQ_ASCII XS_TmodP :: FILE './XS_TmodP' ; +SEQ_ASCII XS_DmodP :: FILE './XS_DmodP' ; +SEQ_ASCII XS_T_DmodP :: FILE './XS_T_DmodP' ; +SEQ_ASCII XS_vacio :: FILE './XS_vacio' ; + +* +*---- +* Definición de variables +*---- + +REAL radio PUR_MOL conc_B conc_Li + T_comb1 T_comb2 T_comb3 T_comb4 + Tvaina Dzry T_cool D_cool Tcanal + T_gap D_gap T_foil T_mod D_mod + D_comb + W238 W235 W234 W239 WO2 WB potencia + PUR_MOL_P conc_B_p conc_Li_p + T_comb1P T_comb2P T_comb3P T_comb4P + T_modP T_coolP D_modP D_coolP D_coolV + buck_geom + ; + +INTEGER perturba ; +STRING buckling ; + +Datos_pertur :: + >>radio<< >>PUR_MOL<< >>conc_B<< >>conc_Li<< + >>T_comb1<< >>T_comb2<< >>T_comb3<< >>T_comb4<< + >>Tvaina<< >>Dzry<< >>T_cool<< >>D_cool<< + >>Tcanal<< >>T_gap<< >>D_gap<< >>T_foil<< + >>T_mod<< >>D_mod<< >>D_comb<< + >>W238<< >>W235<< >>W234<< >>W239<< >>WO2<< >>WB<< + >>potencia<< + >>PUR_MOL_P<< >>conc_B_p<< >>conc_Li_p<< + >>T_comb1P<< >>T_comb2P<< >>T_comb3P<< >>T_comb4P<< + >>T_modP<< >>T_coolP<< >>D_modP<< >>D_coolP<< >>D_coolV<< + >>buck_geom<< >>perturba<< >>buckling<< + ; + +ECHO "BUCKLING = " buckling ; +ECHO "PERTURBACION = " perturba ; +ECHO "GEOMETRICO = " buck_geom ; +*---- +* Geometry parameters. +*---- + +*--------------------------- +*Constantes de la Nuclear wallet card (Octubre 2011) +*--------------------------- +REAL JxMeV := 931.494061 ; !! uma a Mev +REAL EM_D := 13.136 ; !! exceso de masa D +REAL EM_H := 7.289 ; !! exceso de masa H +REAL EM_Li6 := 14.087 ; !! exceso de masa Li6 +REAL EM_Li7 := 14.907 ; !! exceso de masa Li7 +REAL AB_Li6 := 7.59 ; !! abundancia Li6 +REAL AB_Li7 := 92.41 ; !! abundancia Li7 +*------------------------ +REAL MA_D := EM_D JxMeV / 2.0 + ; !! masa atomica D +REAL MA_H := EM_H JxMeV / 1.0 + ; !! masa atomica H +REAL MA_O := 15.999 ; !! Masa atomica O +REAL MA_Li6 := EM_Li6 JxMeV / 6.0 + ; !! masa atomica Li6 +REAL MA_Li7 := EM_Li7 JxMeV / 7.0 + ; !! masa atomica Li7 +REAL MA_Li := MA_Li6 AB_Li6 * MA_Li7 AB_Li7 * + 100. / ; !! m at Li + +*------------ CALCULO DE LOS PORCENTAJES EN PESO ---------------- +REAL M_D := MA_D PUR_MOL * ; !! masa de D +REAL M_H := MA_H 100.0 PUR_MOL - * ; !! masa de H +REAL M_O := MA_O 50.0 * ; !! masa de O +REAL PP_D := 100.0 M_D M_D M_H + M_O + / * ; !! % en peso D +REAL PP_H := 100.0 M_H M_D M_H + M_O + / * ; !! % en peso H +REAL PP_O := 100.0 M_O M_D M_H + M_O + / * ; !! % en peso O +REAL PP_B := conc_B 10000.0 / ; !! % en peso B +REAL PP_Li6 := conc_Li MA_Li6 AB_Li6 * MA_Li 100. * / * 10000. / ; + !!% en peso Li6 +REAL PP_Li7 := conc_Li MA_Li7 AB_Li7 * MA_Li 100. * / * 10000. / ; + !!% en peso Li7 +ECHO "------------- % en peso calculados para D2O -------------------" ; +ECHO "% en peso de D = " PP_D ; +ECHO "% en peso de H = " PP_H ; +ECHO "% en peso de O = " PP_O ; +ECHO "% en peso de B = " PP_B ; +ECHO "% en peso de Li6 = " PP_Li6 ; +ECHO "% en peso de Li7 = " PP_Li7 ; +ECHO "---------------------------------------------------------------" ; + +LIBRARY := LIB: :: + EDIT 3 + NMIX 11 + CTRA WIMS + DEPL LIB: WIMSD4 FIL: iaea + MIXS LIB: WIMSD4 FIL: iaea + + MIX 1 <<T_comb1>> <<D_comb>> + U238 = '8238' <<W238>> 1 SHIB '8238.0' + U235 = '2235' <<W235>> 1 SHIB '2235.0' + U234 = '234' <<W234>> 1 SHIB '234.0' + Pu239 = '6239' <<W239>> 1 SHIB '6239.0' + O2 = '6016' <<WO2>> + '1011' <<WB>> + + MIX 2 COMB 1 1.0 + MIX 3 COMB 1 1.0 + MIX 4 COMB 1 1.0 + +***** vaina +MIX 5 <<Tvaina>> <<Dzry>> + '91' 98.5 + '2056' 0.2 + '52' 0.1 + '118' 1.4 + '6016' 0.12 + +***** refrigerante +MIX 6 <<T_cool>> <<D_cool>> +'6016' <<PP_O>> +'3001' <<PP_H>> +'3002' <<PP_D>> +'1011' <<PP_B>> + '6' <<PP_Li6>> + '7' <<PP_Li7>> + +***** tubo * Canal de Zry4 + MIX 7 <<Tcanal>> <<Dzry>> + '91' 98.5 + '2056' 0.2 + '52' 0.1 + '118' 1.4 + '6016' 0.12 + +***** gap +MIX 8 <<T_gap>> <<D_gap>> +'6016' <<PP_O>> +'3001' <<PP_H>> +'3002' <<PP_D>> +'1011' <<PP_B>> + '6' <<PP_Li6>> + '7' <<PP_Li7>> + +***** foil * Tubo de aislacion de Zry4 + MIX 9 <<T_foil>> <<Dzry>> + '91' 98.5 + '2056' 0.2 + '52' 0.1 + '118' 1.4 + '6016' 0.12 + +***** moderador +MIX 10 <<T_mod>> <<D_mod>> +'6016' <<PP_O>> +'3001' <<PP_H>> +'3002' <<PP_D>> +'1011' <<PP_B>> + '6' <<PP_Li6>> + '7' <<PP_Li7>> + +***** moderador del borde celda +MIX 11 <<T_mod>> <<D_mod>> +'6016' <<PP_O>> +'3001' <<PP_H>> +'3002' <<PP_D>> +'1011' <<PP_B>> + '6' <<PP_Li6>> + '7' <<PP_Li7>> +; + +*---- +* Geometry A2S : GEOMETRY FOR SELF-SHIELDING +* A2T : GEOMETRY FOR TRANSPORT +*---- +A2S := GEO: :: TUBE 5 + R+ REFL + RADIUS 0.0000 + 5.41850 + 5.59150 + 5.76000 + 5.80000 + <<radio>> + MIX 6 7 8 9 10 + CLUSTER CORONA_1 CORONA_2 CORONA_3 CORONA_4 + ::: CORONA_1 := GEO: TUBE 3 + RADIUS 0.00000 + 0.41578 + 0.58800 + 0.65200 + MIX 1 1 5 + NPIN 1 RPIN 0.0000 APIN 0.000000 ; + ::: CORONA_2 := GEO: CORONA_1 MIX 2 2 5 + NPIN 6 RPIN 1.5475 APIN 0.000000 ; + ::: CORONA_3 := GEO: CORONA_1 MIX 3 3 5 + NPIN 12 RPIN 3.0145 APIN 0.261799 ; + ::: CORONA_4 := GEO: CORONA_1 MIX 4 4 5 + NPIN 18 RPIN 4.4870 APIN 0.000000 ; + ; + +A2T := GEO: :: TUBE 11 + R+ REFL + RADIUS 0.00000 + 0.86213 + 1.54750 + 2.28100 + 3.01450 + 3.75075 + 4.48700 + 5.41850 + 5.59150 + 5.76000 + 5.80000 + <<radio>> + MIX 6 6 6 6 6 6 6 7 8 9 10 + SPLITR 1 1 1 1 1 1 1 1 1 1 20 + CLUSTER CORONA_1 CORONA_2 CORONA_3 CORONA_4 + ::: CORONA_1 := GEO: TUBE 3 + RADIUS 0.00000 + 0.41578 + 0.58800 + 0.65200 + MIX 1 1 5 + NPIN 1 RPIN 0.0000 APIN 0.000000 ; + ::: CORONA_2 := GEO: CORONA_1 MIX 2 2 5 + NPIN 6 RPIN 1.5475 APIN 0.000000 ; + ::: CORONA_3 := GEO: CORONA_1 MIX 3 3 5 + NPIN 12 RPIN 3.0145 APIN 0.261799 ; + ::: CORONA_4 := GEO: CORONA_1 MIX 4 4 5 + NPIN 18 RPIN 4.4870 APIN 0.000000 ; + ; + +*---- +* Graficos +*---- + GRAF_A2S.PS := PSP: A2S :: EDIT 2 FILL HSB TYPE MIXTURE ; + GRAF_A2T.PS := PSP: A2T :: EDIT 2 FILL HSB TYPE REGION ; + +*---- +* Tracking para Self-Shielding EXCELT +*---- + +TRACK_S INTLIN_S := EXCELT: A2S :: + TITLE 'CNA-UII: 37 B. UNAT - referencia' + EDIT 3 + MAXR 60 +ALLG BATCH 100 +TRAK TISO 18 20.0 + SYMM 12 + ; + +*---- +* Tracking for Transport: EXCELT +*---- +TRACK_T INTLIN_T := EXCELT: A2T :: + TITLE 'CNA-UII: 37 B. UNAT - evolución perturbada' + EDIT 3 + MAXR 60 +ALLG BATCH 100 +TRAK TISO 18 20.0 + SYMM 12 + ; + +LIBRARY_S := SHI: LIBRARY TRACK_S INTLIN_S :: + EDIT 3 + ; + +Pij := ASM: LIBRARY_S TRACK_T INTLIN_T :: + EDIT 2 + PIJ ! default + PNOR HELI ! default + ; + +IF buckling "critico" = THEN + FLUX := FLU: Pij LIBRARY_S TRACK_T :: + EDIT 2 + TYPE B B0TR BUCK 0.0 + ; +ELSEIF buckling "geometrico" = THEN + FLUX := FLU: Pij LIBRARY_S TRACK_T :: + EDIT 2 + TYPE K B0TR BUCK <<buck_geom>> + ; +ENDIF ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.078719 ; + +*---- +* Calculo de self-shielding y transporte para el +* caso de referencia a quemado Q=0 +*---- + +ECHO "PERTURBACION DE BIBLIOTECA DE REFERENCIA " ; + + LIBRARY_N := LIBRARY_S ; + LIBRARY_N := LIB: LIBRARY_N :: + EDIT 3 + CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 6 <<T_cool>> <<D_coolP>> + ; + +LIBRARY_N := SHI: LIBRARY LIBRARY_N TRACK_S INTLIN_S :: + EDIT 3 + ; + +Pij := DELETE: Pij ; + +Pij := ASM: LIBRARY_N TRACK_T INTLIN_T :: + EDIT 0 + PIJ + PNOR HELI + ; + +FLUX1 := FLUX ; +FLUX1 := FLU: FLUX1 Pij LIBRARY_N TRACK_T :: + TYPE K B0TR BUCK IDEM ; +assertS FLUX1 :: 'K-EFFECTIVE' 1 1.079210 ; + +*---- +* Condensacion de XS a 2 grupos y homogenizacion +*---- + +EDITION1 := EDI: FLUX LIBRARY_N TRACK_T :: + EDIT 3 + SAVE + COND 0.625 ! para 2 grupos + MERGE COMP + MICR ALL + ; + +*---- +* Se crean los archivos EDIT para extraer las XS +*---- +XS_DcoolP := EDITION1 ; + +*---- +* Se eliminan bibliotecas del directorio de trabajo +*---- +FLUX Pij TRACK_T INTLIN_T INTLIN_S := DELETE: +FLUX Pij TRACK_T INTLIN_T INTLIN_S ; + +END: ; + +QUIT "LIST" . diff --git a/Dragon/data/twlup_proc/TCWU41.c2m b/Dragon/data/twlup_proc/TCWU41.c2m new file mode 100644 index 0000000..e2c357c --- /dev/null +++ b/Dragon/data/twlup_proc/TCWU41.c2m @@ -0,0 +1,189 @@ +*--- +* Regression test 6 for MCCG: +*--- +* TEST CASE TCWU41 +* CANDU-6 CARTESIAN CELL +* WIMSD4 69 GROUPS LIBRARY FILE WNEALIB +* STAINLESS STELL RODS IN 3D SUPERCELL +*---- + +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY LIBRARY2 CANDU6F CANDU6S TRACK SYS FLUX EDITION0 EDITION + BCO BCI TEMP ; +SEQ_BINARY + INTLIN ; +MODULE + GEO: EXCELT: MCCGT: LIB: SHI: ASM: FLU: EDI: DELETE: END: ; +PROCEDURE assertS ; +INTEGER split1 splitz := 1 2 ; +INTEGER split2 := split1 2 * ; +INTEGER nang := 4 ; +REAL dens := 2.5 ; +INTEGER stis := 1 ; +*---- +* Microscopic cross sections from file WNEALIB format WIMSD4 +*---- +PROCEDURE TCWU05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWU05Lib :: <<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 +* Transport calculation +* Flux calculation +*---- +TRACK INTLIN := EXCELT: CANDU6S :: + TITLE 'TCWU41: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 14 + ALLG BATCH 100 + TRAK TISO 29 20.0 SYMM 4 ; +LIBRARY := SHI: LIBRARY TRACK INTLIN :: NOLJ + EDIT 0 ; +TRACK INTLIN := DELETE: TRACK INTLIN ; +TRACK INTLIN := EXCELT: CANDU6F :: + TITLE 'TCWU41: CANDU-6 CARTESIAN FUEL TEMP= 941.29' + EDIT 0 MAXR 32 + ALLG BATCH 100 + TRAK TISO 29 20.0 SYMM 4 ; +TRACK := MCCGT: TRACK INTLIN :: + EDIT 1 LCMD 2 + AAC 50 TMT EPSI 1E-5 + MAXI 1 KRYL 0 SCR 0 HDD 0.0 ; +SYS := ASM: LIBRARY TRACK INTLIN :: + EDIT 0 ARM ; +FLUX := FLU: SYS LIBRARY TRACK INTLIN :: + TYPE K THER 5 5E-5 EXTE 100 1E-5 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.120511 ; +*---- +* Microscopic cross sections from WIMSLIB for stainless steel +* MIX 1,2,3 from EDI: : fuel, structure material and moderator +* MIX 4 stainless steel rods +*---- +EDITION0 := EDI: LIBRARY TRACK FLUX :: + EDIT 0 MERGE MIX 1 2 2 2 3 1 1 1 1 1 SAVE ON SSRODS ; +EDITION := EDITION0 :: STEP UP SSRODS ; +LIBRARY2 := LIB: EDITION :: + EDIT 2 NMIX 4 CTRA WIMS + MIXS LIB: WIMSD4 FIL: iaea + MIX 4 345.66 Fe56 = '2056' 6.19027E-2 + Cr52 = '52' 1.56659E-3 Ni58 = '58' 6.83337E-3 + Si29 = '29' 7.79072E-4 C12 = '2012' 1.46552E-4 + Mn55 = '55' 1.25431E-3 + ; +EDITION0 EDITION TRACK INTLIN SYS FLUX := DELETE: + EDITION0 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 SPLITX <<split1>> + MESHY 0.0 7.14375 SPLITY <<split1>> + MESHZ -8.25500 8.25500 SPLITZ <<splitz>> ; + ::: MX := GEO: M MESHX -7.14375 +7.14375 SPLITX <<split2>> ; + ::: MXY := GEO: MX MESHY -7.14375 +7.14375 SPLITY <<split2>> ; + ::: BX := GEO: CARCELY 2 1 MIX 3 3 3 + MESHX -7.14375 7.14375 SPLITX <<split2>> + MESHY 0.0 7.14375 SPLITY <<split1>> + MESHZ -8.25500 8.25500 SPLITZ <<splitz>> + RADIUS 0.0 3.5100 3.8100 ; + ::: BXY := GEO: BX MESHY -7.14375 +7.14375 SPLITY <<split2>> ; + ::: FXY := GEO: CARCELZ 2 1 MIX 1 2 3 + MESHX -7.14375 7.14375 SPLITX <<split2>> + MESHY -7.14375 7.14375 SPLITY <<split2>> + MESHZ -8.25500 8.25500 SPLITZ <<splitz>> + RADIUS 0.0 5.16890 6.58750 ; + ; +BCI := GEO: BCO :: + ::: BX := GEO: BX MIX 3 4 3 ; + ::: BXY := GEO: BXY MIX 3 4 3 ; + ; +*---- +* Flux calculation +* Homogenized properties for rod out +*---- +TRACK INTLIN := EXCELT: BCO :: + EDIT 1 MAXR 2000 TRAK TISO <<nang>> <<dens>> ; +TEMP := TRACK ; +TRACK := MCCGT: TRACK INTLIN :: +EDIT 3 EPSI 1E-8 MAXI 100 +DIFC TMT STIS <<stis>> ; +SYS := ASM: LIBRARY2 TRACK INTLIN :: + EDIT 2 ARM ; +FLUX := FLU: SYS LIBRARY2 TRACK INTLIN :: + EDIT 1 TYPE K THER 5 5E-5 EXTE 10 1E-5 ; +!assertS FLUX :: 'K-EFFECTIVE' 1 1.118734 ; +TRACK := DELETE: TRACK ; +TRACK := TEMP ; +TRACK := MCCGT: TRACK INTLIN :: +EDIT 3 EPSI 1E-7 +AAC 50 TMT +SCR 0 +KRYL 0 MAXI 1 STIS <<stis>> ; +FLUX := FLU: FLUX + SYS LIBRARY2 TRACK INTLIN :: + EDIT 1 TYPE K + THER 5 5E-5 EXTE 100 1E-5 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.127860 ; +EDITION := EDI: LIBRARY2 TRACK FLUX :: + EDIT 2 MERG COMP COND 4.0 SAVE ON 'NOBC' ; +SYS TRACK INTLIN := DELETE: SYS TRACK INTLIN ; +*---- +* Flux calculation +* Homogenized properties for rod in +*---- +TRACK INTLIN := EXCELT: BCI :: + EDIT 0 MAXR 2000 TRAK TISO <<nang>> <<dens>> ; +TRACK := MCCGT: TRACK INTLIN :: +EDIT 2 EPSI 1E-7 +AAC 50 TMT +SCR 0 +KRYL 0 MAXI 1 STIS <<stis>> ; +SYS := ASM: LIBRARY2 TRACK INTLIN :: + EDIT 2 ARM ; +FLUX := FLU: FLUX +SYS LIBRARY2 TRACK INTLIN :: + EDIT 1 TYPE K + THER 5 5E-5 EXTE 100 1E-5 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.042392 ; +EDITION := EDI: EDITION LIBRARY2 TRACK FLUX :: + EDIT 2 MERG COMP COND 4.0 STAT DELS REFE 'NOBC' ; +*--- +* Memory Release +*--- +TRACK INTLIN SYS := DELETE: TRACK INTLIN SYS ; +ECHO "test TCWU41 completed" ; +END: ; +QUIT "LIST" . |