From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001 From: stainer_t Date: Mon, 8 Sep 2025 13:48:49 +0200 Subject: Initial commit from Polytechnique Montreal --- Dragon/data/twimsE_proc/TCWE01.c2m | 129 +++++++++++++ Dragon/data/twimsE_proc/TCWE02.c2m | 179 ++++++++++++++++++ Dragon/data/twimsE_proc/TCWE03.c2m | 172 +++++++++++++++++ Dragon/data/twimsE_proc/TCWE04.c2m | 151 +++++++++++++++ Dragon/data/twimsE_proc/TCWE05.c2m | 192 +++++++++++++++++++ Dragon/data/twimsE_proc/TCWE05Lib.c2m | 68 +++++++ Dragon/data/twimsE_proc/TCWE06.c2m | 155 ++++++++++++++++ Dragon/data/twimsE_proc/TCWE07.c2m | 142 ++++++++++++++ Dragon/data/twimsE_proc/TCWE08.c2m | 121 ++++++++++++ Dragon/data/twimsE_proc/TCWE09.c2m | 74 ++++++++ Dragon/data/twimsE_proc/TCWE10.c2m | 104 +++++++++++ Dragon/data/twimsE_proc/TCWE11.c2m | 171 +++++++++++++++++ Dragon/data/twimsE_proc/TCWE12.c2m | 107 +++++++++++ Dragon/data/twimsE_proc/TCWE13.c2m | 146 +++++++++++++++ Dragon/data/twimsE_proc/TCWE14.c2m | 121 ++++++++++++ Dragon/data/twimsE_proc/TCWE15.c2m | 159 ++++++++++++++++ Dragon/data/twimsE_proc/TCWE17.c2m | 337 ++++++++++++++++++++++++++++++++++ Dragon/data/twimsE_proc/TCWE17Lib.c2m | 97 ++++++++++ Dragon/data/twimsE_proc/TCWE31.c2m | 172 +++++++++++++++++ Dragon/data/twimsE_proc/TCWE32.c2m | 97 ++++++++++ 20 files changed, 2894 insertions(+) create mode 100644 Dragon/data/twimsE_proc/TCWE01.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE02.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE03.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE04.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE05.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE05Lib.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE06.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE07.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE08.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE09.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE10.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE11.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE12.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE13.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE14.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE15.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE17.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE17Lib.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE31.c2m create mode 100644 Dragon/data/twimsE_proc/TCWE32.c2m (limited to 'Dragon/data/twimsE_proc') diff --git a/Dragon/data/twimsE_proc/TCWE01.c2m b/Dragon/data/twimsE_proc/TCWE01.c2m new file mode 100644 index 0000000..a18e4da --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE01.c2m @@ -0,0 +1,129 @@ +*---- +* TEST CASE TCWE01 +* MOSTELLER BENCHMARK: 1-D ANNULAR CELL AND 2-D CARTESIAN CELL +* WIMSE 69 GROUPS LIBRARY FILE WIMS-E +* +* 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 ; +REAL P_Zr90 P_Zr91 P_Zr92 P_Zr94 P_Zr96 := + 51.45E-02 11.22E-02 17.15E-02 17.38E-02 2.80E-02 ; +REAL N_Zr0 := 3.83243E-2 ; +REAL N_Zr90 N_Zr91 N_Zr92 N_Zr94 N_Zr96 := + N_Zr0 P_Zr90 * N_Zr0 P_Zr91 * N_Zr0 P_Zr92 * N_Zr0 P_Zr94 * + N_Zr0 P_Zr96 * ; + +REAL P_B10 P_B11 := 0.199 0.801 ; +REAL N_B0 := 1.02133E-5 ; +REAL N_B10 N_B11 := N_B0 P_B10 * N_B0 P_B11 * ; +*---- +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 600.0 O16 = '8016' 4.61309E-2 + U235 = '92235' 1.66078E-4 1 + U238 = '92238' 2.28994E-2 1 + MIX 2 600.0 + Zr90 = '40090' <> + Zr91 = '40091' <> 2 + Zr92 = '40092' <> + Zr94 = '40094' <> + Zr96 = '40096' <> + MIX 3 600.0 + H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2 + B10 = '5010' <> B11 = '5011' <> + ; +*---- +* 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 'TCWE01: 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.8289803 ; +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 'TCWE01: 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.8290998 ; +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 'TCWE01: 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.8302385 ; +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 'TCWE01: 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.8325990 ; +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 TCWE01 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE02.c2m b/Dragon/data/twimsE_proc/TCWE02.c2m new file mode 100644 index 0000000..afb01e4 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE02.c2m @@ -0,0 +1,179 @@ +*---- +* TEST CASE TCWE02 +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* WIMSE 69 GROUPS LIBRARY FILE WIMS-E +* +* 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 WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 579.9 + H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2 + B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05 + MIX 2 579.9 O16 = '8016' 3.06711E-4 + Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4 + Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2 + Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03 + Zr96 = '40096' 1.17214E-03 + MIX 3 579.9 + H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2 + Cr52 = '24052' 4.79927E-5 + Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4 + Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07 + Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07 + Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07 + Al27 = '13027' 2.35231E-6 + Mn55 = '25055' 4.15901E-7 + B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05 + Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04 + Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04 + Zr96 = '40096' 2.49880E-05 + MIX 4 579.9 O16 = '8016' 2.87335E-4 + Cr52 = '24052' 7.07291E-5 Fe56 = '26056' 1.38298E-4 + Zr90 = '40090' 2.01774E-02 Zr91 = '40091' 4.40020E-03 2 + Zr92 = '40092' 6.72580E-03 Zr94 = '40094' 6.81600E-03 + Zr96 = '40096' 1.09809E-03 + MIX 5 579.9 + H1H2O = '1001' 4.71346E-2 O16H2O = '8016' 2.35673E-2 + Cr52 = '24052' 2.24991E-5 + Fe56 = '26056' 2.09013E-5 Ni58 = '28058' 5.32188E-5 + Mo92 = '42092' 2.77278E-07 Mo95 = '42095' 1.73892E-07 + Mo96 = '42096' 3.00446E-07 Mo97 = '42097' 1.81369E-07 + Mo98 = '42098' 4.59801E-07 Mo100 = '42100' 1.84435E-07 + Al27 = '13027' 1.10277E-6 + Mn55 = '25055' 1.94976E-7 + B10 = '5010' 4.68840E-06 B11 = '5011' 1.88714E-05 + Zr90 = '40090' 2.15252E-04 Zr91 = '40091' 4.69413E-05 + Zr92 = '40092' 7.17508E-05 Zr94 = '40094' 7.27131E-05 + Zr96 = '40096' 1.17144E-05 + MIX 6 579.9 + H1H2O = '1001' 4.71676E-2 O16H2O = '8016' 2.35838E-2 + Cr52 = '24052' 2.11122E-5 + Fe56 = '26056' 1.96130E-5 Ni58 = '28058' 4.99383E-5 + Mo92 = '42092' 2.60187E-07 Mo95 = '42095' 1.63174E-07 + Mo96 = '42096' 2.81927E-07 Mo97 = '42097' 1.70190E-07 + Mo98 = '42098' 4.31460E-07 Mo100 = '42100' 1.73067E-07 + Al27 = '13027' 1.03479E-6 + Mn55 = '25055' 1.82957E-7 + B10 = '5010' 4.69148E-06 B11 = '5011' 1.88838E-05 + Zr90 = '40090' 2.01984E-04 Zr91 = '40091' 4.40478E-05 + Zr92 = '40092' 6.73280E-05 Zr94 = '40094' 6.82309E-05 + Zr96 = '40096' 1.09923E-05 + MIX 7 579.9 + H1H2O = '1001' 4.72020E-2 O16H2O = '8016' 2.36010E-2 + Cr52 = '24052' 1.96591E-5 + Fe56 = '26056' 1.82630E-5 Ni58 = '28058' 4.65011E-5 + Mo92 = '42092' 2.42278E-07 Mo95 = '42095' 1.51943E-07 + Mo96 = '42096' 2.62522E-07 Mo97 = '42097' 1.58476E-07 + Mo98 = '42098' 4.01763E-07 Mo100 = '42100' 1.61155E-07 + Al27 = '13027' 9.63569E-7 + Mn55 = '25055' 1.70365E-7 + B10 = '5010' 4.69469E-06 B11 = '5011' 1.88967E-05 + Zr90 = '40090' 1.88082E-04 Zr91 = '40091' 4.10161E-05 + Zr92 = '40092' 6.26939E-05 Zr94 = '40094' 6.35347E-05 + Zr96 = '40096' 1.02357E-05 + MIX 8 933.6 O16 = '8016' 4.49355E-2 + U235 = '92235' 7.39237E-4 1 + U238 = '92238' 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 'TCWE02: 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.258811 ; + +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.182448E-01 ; +assertS ISOT :: 'NWT0' 2 1.297516E-01 ; + +ECHO "test TCWE02 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE03.c2m b/Dragon/data/twimsE_proc/TCWE03.c2m new file mode 100644 index 0000000..53d1dae --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE03.c2m @@ -0,0 +1,172 @@ +*---- +* TEST CASE TCWE03 +* MULTICELL HEXAGONAL ASSEMBLY WITH POISON +* WIMS-E 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 WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 11 CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 579.9 + H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2 + B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05 + MIX 2 579.9 O16 = '8016' 3.06711E-4 + Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4 + Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2 + Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03 + Zr96 = '40096' 1.17214E-03 + MIX 3 579.9 + H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2 + Cr52 = '24052' 4.79927E-5 + Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4 + Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07 + Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07 + Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07 + Mn55 = '25055' 4.15901E-7 + B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05 + Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04 + Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04 + Zr96 = '40096' 2.49880E-05 + MIX 4 933.6 O16 = '8016' 4.49355E-2 + U235 = '92235' 7.39237E-4 1 + U238 = '92238' 2.17285E-2 1 + MIX 5 579.9 In115 = '49115' 7.57464E-3 + Cd113 = '48113' 2.62493E-3 + Ag109 = '47109' 4.49188E-2 + MIX 6 579.9 Cr52 = '24052' 1.52702E-2 + Fe56 = '26056' 5.57670E-2 Ni58 = '28058' 7.51418E-3 + Mn55 = '25055' 8.02943E-4 + MIX 7 579.9 + H1H2O = '1001' 3.06466E-2 O16H2O = '8016' 1.53233E-2 + Fe56 = '26056' 5.27485E-5 Cr52 = '24052' 2.69769E-5 + B10 = '5010' 3.04623E-06 B11 = '5011' 1.22615E-05 + Zr90 = '40090' 7.69589E-03 Zr91 = '40091' 1.67829E-03 2 + Zr92 = '40092' 2.56530E-03 Zr94 = '40094' 2.59970E-03 + Zr96 = '40096' 4.18824E-04 + MIX 8 579.9 + H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2 + Cr52 = '24052' 4.79927E-5 + Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4 + Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07 + Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07 + Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07 + Al27 = '13027' 2.35231E-6 + Mn55 = '25055' 4.15901E-7 + B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05 + Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04 + Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04 + Zr96 = '40096' 2.49880E-05 + MIX 9 579.9 O16 = '8016' 2.87335E-4 + Cr52 = '24052' 7.07291E-5 Fe56 = '26056' 1.38298E-4 + Zr90 = '40090' 2.01774E-02 Zr91 = '40091' 4.40020E-03 2 + Zr92 = '40092' 6.72580E-03 Zr94 = '40094' 6.81600E-03 + Zr96 = '40096' 1.09809E-03 + MIX 10 579.9 + H1H2O = '1001' 4.71346E-2 O16H2O = '8016' 2.35673E-2 + Cr52 = '24052' 2.24991E-5 + Fe56 = '26056' 2.09013E-5 Ni58 = '28058' 5.32188E-5 + Mo92 = '42092' 2.77278E-07 Mo95 = '42095' 1.73892E-07 + Mo96 = '42096' 3.00446E-07 Mo97 = '42097' 1.81369E-07 + Mo98 = '42098' 4.59801E-07 Mo100 = '42100' 1.84435E-07 + Al27 = '13027' 1.10277E-6 + Mn55 = '25055' 1.94976E-7 + B10 = '5010' 4.68840E-06 B11 = '5011' 1.88714E-05 + Zr90 = '40090' 2.15252E-04 Zr91 = '40091' 4.69413E-05 + Zr92 = '40092' 7.17508E-05 Zr94 = '40094' 7.27131E-05 + Zr96 = '40096' 1.17144E-05 + MIX 11 579.9 + H1H2O = '1001' 4.71676E-2 O16H2O = '8016' 2.35838E-2 + Cr52 = '24052' 2.11122E-5 + Fe56 = '26056' 1.96130E-5 Ni58 = '28058' 4.99383E-5 + Mo92 = '42092' 2.60187E-07 Mo95 = '42095' 1.63174E-07 + Mo96 = '42096' 2.81927E-07 Mo97 = '42097' 1.70190E-07 + Mo98 = '42098' 4.31460E-07 Mo100 = '42100' 1.73067E-07 + Al27 = '13027' 1.03479E-6 + Mn55 = '25055' 1.82957E-7 + B10 = '5010' 4.69148E-06 B11 = '5011' 1.88838E-05 + Zr90 = '40090' 2.01984E-04 Zr91 = '40091' 4.40478E-05 + Zr92 = '40092' 6.73280E-05 Zr94 = '40094' 6.82309E-05 + Zr96 = '40096' 1.09923E-05 + ; +*---- +* 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 'TCWE03: 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.6724899 ; +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.483070E+00 ; +assertS ISOT :: 'NWT0' 2 1.400084E-01 ; + +ECHO "test TCWE03 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE04.c2m b/Dragon/data/twimsE_proc/TCWE04.c2m new file mode 100644 index 0000000..e9ab59f --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE04.c2m @@ -0,0 +1,151 @@ +*---- +* TEST CASE TCWE04 +* WIMS-E 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 := + 1.0 27.1739 67.9348 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + DEPL LIB: WIMSE FIL: WIMS-E + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 600.0 O16 = '8016' 4.61309E-2 + U235 = '92235' 1.66078E-4 1 + U238 = '92238' 2.28994E-2 1 + U236 = '92236' 0.0 1 + Pu239 = '94239' 0.0 1 + MIX 2 600.0 + Zr90 = '40090' 1.97179E-02 Zr91 = '40091' 4.29999E-03 2 + Zr92 = '40092' 6.57262E-03 Zr94 = '40094' 6.66076E-03 + Zr96 = '40096' 1.07308E-03 + MIX 3 600.0 + H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2 + B10 = '5010' 2.03245E-06 B11 = '5011' 8.18085E-06 + ; +*---- +* 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 'TCWE04: MOSTELLER BENCHMARK WITH BURNUP' + EDIT 1 MAXR 3 ; +LIBRARY := SHI: LIBRARY TRACKS :: EDIT 0 NOLJ ; +TRACK := SYBILT: MOSTELA :: + TITLE 'TCWE04: 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 <> DAY POWR 36.8 + DEPL <> <> DAY POWR 36.8 + SET <> DAY ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK :: + SAVE <> DAY POWR 36.8 + NOEX DEPL <> <> DAY POWR 36.8 + SET <> 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 <> DAY POWR 36.8 + DEPL <> <> DAY POWR 36.8 + SET <> 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 <> DAY POWR 36.8 ; + DATABASE := COMPO: DATABASE EDITION BURNUP :: + EDIT 3 + SET <> DAY + ; + EVALUATE step2 step3 step4 step5 := + step3 step4 step5 step2 ; + EVALUATE istep := istep 1 + ; +ENDWHILE ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.7397701 ; + +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.840256E+01 ; +assertS ISOT :: 'NWT0' 2 1.686708E+01 ; + +ECHO "test TCWE04 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE05.c2m b/Dragon/data/twimsE_proc/TCWE05.c2m new file mode 100644 index 0000000..93d8f95 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE05.c2m @@ -0,0 +1,192 @@ +*---- +* TEST CASE TCWE05 +* CANDU-6 ANNULAR CELL +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +PROCEDURE TCWE05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWE05Lib :: <> ; +*---- +* 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 'TCWE05: 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 'TCWE05: 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.124888 ; +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 <> <> DAY POWR <> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + NOEX DEPL <> <> DAY POWR <> ; + 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 <> DAY POWR <> ; + DATABASE := COMPO: DATABASE EDITION BURNUP LIBRARY :: + EDIT 3 + STEP UP * + SET <> 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.9690224 ; +*---- +* 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.663608E0 ; +assertS ISOT :: 'NWT0' 2 2.579294E1 ; +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.295185E-2 ; +assertS ISOT :: 'NWT0' 2 3.190355E-1 ; + +INTLINF INTLINS := DELETE: INTLINF INTLINS ; +ECHO "test TCWE05 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE05Lib.c2m b/Dragon/data/twimsE_proc/TCWE05Lib.c2m new file mode 100644 index 0000000..9114411 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE05Lib.c2m @@ -0,0 +1,68 @@ +*---- +* Procedure TCWE05Lib +* Create Library for test CASE TCWE05 +* Calling : +* LIBRARY := TCWE05Lib :: iedit ; +* with : +* LIBRARY = Linked list containing the result of LIB: for +* TCWE05 +* 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + EDIT <> + NMIX 17 CTRA WIMS + DEPL LIB: WIMSE FIL: WIMS-E + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 560.66 0.81212 O16 = '8016' 7.99449E-1 + D2D2O = '1002' 1.99768E-1 H1H2O = '1001' 7.83774E-4 + MIX 2 560.66 6.57 Nb93 = '41093' 2.50000 + B10 = '5010' 4.17900E-05 B11 = '5011' 1.68210E-04 + Zr90 = '40090' 5.01638E+01 Zr91 = '40091' 1.09395E+01 2 + Zr92 = '40092' 1.67213E+01 Zr94 = '40094' 1.69455E+01 + Zr96 = '40096' 2.73000E+00 + MIX 3 345.66 0.0014 He4 = '2004' 1.00000E+2 + MIX 4 345.66 6.44 Fe56 = '26056' 1.60000E-1 + Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1 + B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04 + Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2 + Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01 + Zr96 = '40096' 2.79188E+00 + MIX 5 345.66 1.082885 O16 = '8016' 7.98895E-1 + D2D2O = '1002' 2.01016E-1 H1H2O = '1001' 8.96000E-5 + MIX 6 941.29 10.4375010 O16 = '8016' 1.18473E+1 + Xe135 = '54135' 0.0 + U235 = '92235' 6.27118E-1 1 + U238 = '92238' 8.75256E+1 1 + U236 = '92236' 0.0 1 + Pu239 = '94239' 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 = '26056' 1.60000E-1 + Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1 + B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04 + Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2 + Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01 + Zr96 = '40096' 2.79188E+00 + 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/twimsE_proc/TCWE06.c2m b/Dragon/data/twimsE_proc/TCWE06.c2m new file mode 100644 index 0000000..287c539 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE06.c2m @@ -0,0 +1,155 @@ +*---- +* TEST CASE TCWE06 +* CANDU-6 CARTESIAN CELL +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +PROCEDURE TCWE05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWE05Lib :: <> ; +*---- +* 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 'TCWE06: 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 'TCWE06: 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.124544 ; +*---- +* 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: WIMSE FIL: WIMS-E + MIX 4 345.66 Fe56 = '26056' 6.19027E-2 + Cr52 = '24052' 1.56659E-3 Ni58 = '28058' 6.83337E-3 + Si29 = '14029' 7.79072E-4 C12 = '6000' 1.46552E-4 + Mn55 = '25055' 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.131641 ; +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.046070 ; +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 TCWE06 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE07.c2m b/Dragon/data/twimsE_proc/TCWE07.c2m new file mode 100644 index 0000000..afbf9c3 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE07.c2m @@ -0,0 +1,142 @@ +*---- +* TEST CASE TCWE07 +* CANDU-6 CARTESIAN CELL +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +PROCEDURE TCWE05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWE05Lib :: <> ; +*---- +* 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 'TCWE07: 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 'TCWE07: 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.124517 ; +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.116083 ; +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.116019 ; +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.116024 ; +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 'TCWE07: 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 'TCWE07: 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.144657 ; +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.136610 ; +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.136823 ; +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.136833 ; +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 TCWE07 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE08.c2m b/Dragon/data/twimsE_proc/TCWE08.c2m new file mode 100644 index 0000000..cb329ad --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE08.c2m @@ -0,0 +1,121 @@ +*---- +* TEST CASE TCWW08 +* HOMOGENEOUS DEPLETION CASE +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 1 CTRA WIMS + DEPL LIB: WIMSE FIL: WIMS-E + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 300.0 H1 = '1001' 2.00000E+1 + U235 = '92235' 1.0 1 + U236 = '92238' 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.843108 ; +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 <> DEPL <> <> DAY + POWR <> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX TRACK :: + EDIT <> NOEX DEPL <> <> DAY + POWR <> ; + 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.322846 ; +ECHO "test TCWE08 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE09.c2m b/Dragon/data/twimsE_proc/TCWE09.c2m new file mode 100644 index 0000000..a50e6d1 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE09.c2m @@ -0,0 +1,74 @@ +*---- +* TEST CASE TCWE09 +* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL +* WIMS-E 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 WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 3 CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 600.0 O16 = '8016' 4.61309E-2 + U235 = '92235' 1.66078E-4 1 + U238 = '92238' 2.28994E-2 1 + MIX 2 600.0 + Zr90 = '40090' 1.97179E-02 Zr91 = '40091' 4.29999E-03 2 + Zr92 = '40092' 6.57262E-03 Zr94 = '40094' 6.66076E-03 + Zr96 = '40096' 1.07308E-03 + MIX 3 600.0 + H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2 + B10 = '5010' 2.03245E-06 B11 = '5011' 8.18085E-06 + ; +*---- +* 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 'TCWE09: 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.8289741 ; +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 'TCWE09: 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 9.918292E-03 ; +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 TCWE09 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE10.c2m b/Dragon/data/twimsE_proc/TCWE10.c2m new file mode 100644 index 0000000..0005914 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE10.c2m @@ -0,0 +1,104 @@ +*---- +* TEST CASE TCWE10 +* MOSTELLER BENCHMARK FOR 1-D ANNULAR CELL +* WIMS-E 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 WIMS-E format WIMSE +* Fixed source of 1.0E5 in group 6 +*---- +LIBRARY := LIB: :: + EDIT 0 NMIX 3 CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 600.0 O16 = '8016' 4.61309E-2 + U235 = '92235' 1.66078E-4 1 + U238 = '92238' 2.28994E-2 1 + MIX 2 600.0 + Zr90 = '40090' 1.97179E-02 Zr91 = '40091' 4.29999E-03 2 + Zr92 = '40092' 6.57262E-03 Zr94 = '40094' 6.66076E-03 + Zr96 = '40096' 1.07308E-03 + MIX 3 600.0 + H1H2O = '1001' 4.42326E-2 O16H2O = '8016' 2.21163E-2 + B10 = '5010' 2.03245E-06 B11 = '5011' 8.18085E-06 + ; +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 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 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 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 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 'TCWE10: 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.8289741 ; +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 2.435309E+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 TCWE10 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE11.c2m b/Dragon/data/twimsE_proc/TCWE11.c2m new file mode 100644 index 0000000..23a73cf --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE11.c2m @@ -0,0 +1,171 @@ +*---- +* TEST CASE TCWE11 +* CANDU-6 ANNULAR CELL +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +PROCEDURE TCWE05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWE05Lib :: <> ; +*---- +* 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 'TCWE11: 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 'TCWE11: 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.124889 ; +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 <> <> DAY POWR <> ; +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.9580293 ; +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.124906 ; +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 <> <> DAY POWR <> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + EDIT 3 NOEX DEPL <> <> DAY POWR <> ; + 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.079448 ; +res := EDITION ; +EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE: + EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ; +ECHO "test TCWE11 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE12.c2m b/Dragon/data/twimsE_proc/TCWE12.c2m new file mode 100644 index 0000000..2ea0f32 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE12.c2m @@ -0,0 +1,107 @@ +*---- +* TEST CASE TCWE12 +* WIMS-E 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 + WU235F WU238F WO16F ; +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST + LIBRARY ISOT ; +MODULE + LIB: INFO: END: ; +PROCEDURE assertS ; +*---- +* Get Coolant properties +*---- +ECHO + "Input - Coolant temperature (K) " TempCool ; +ECHO + "Input - D2/(D2+H1) Weight % ratio in Coolant" Purity ; +INFO: :: + TMP: <> K + PUR: <> WGT% + CALC DENS WATER >>DensCool<< + LIB: WIMSE FIL: WIMS-E + ISO: 3 '1001' '1002' '8016' + 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 ; +*---- +* 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: <> WGT% + LIB: WIMSE FIL: WIMS-E + ISO: 3 '92235' '92238' '8016' + 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 WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 5 CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 <> <> O16 = '8016' <> + D2D2O = '1002' <> H1H2O = '1001' <> + MIX 2 <> <> O16 = '8016' <> + D2D2O = '1002' <> H1H2O = '1001' <> + MIX 3 <> <> O16 = '8016' <> + U235 = '92235' <> 1 + U238 = '92238' <> 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 50 4.918399E-01 ; +ECHO "test TCWE12 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE13.c2m b/Dragon/data/twimsE_proc/TCWE13.c2m new file mode 100644 index 0000000..21bf25c --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE13.c2m @@ -0,0 +1,146 @@ +*---- +* TEST CASE TCWE13 +* 17 X 17 PWR ASSEMBLY WITHOUT POISON +* WIMSE 69 GROUPS LIBRARY FILE WIMS-E +* +* 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 WIMS-E format WIMSE +*---- +LIBRARY := LIB: :: + NMIX 8 CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 579.9 + H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2 + B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05 + MIX 2 579.9 O16 = '8016' 3.06711E-4 + Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4 + Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2 + Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03 + Zr96 = '40096' 1.17214E-03 + MIX 3 579.9 + H1H2O = '1001' 4.65292E-2 O16H2O = '8016' 2.32646E-2 + Cr52 = '24052' 4.79927E-5 + Fe56 = '26056' 4.45845E-5 Ni58 = '28058' 1.13521E-4 + Mo92 = '42092' 5.91461E-07 Mo95 = '42095' 3.70930E-07 + Mo96 = '42096' 6.40880E-07 Mo97 = '42097' 3.86878E-07 + Mo98 = '42098' 9.80802E-07 Mo100 = '42100' 3.93419E-07 + Al27 = '13027' 2.35231E-6 + Mn55 = '25055' 4.15901E-7 + B10 = '5010' 4.63194E-06 B11 = '5011' 1.86442E-05 + Zr90 = '40090' 4.59154E-04 Zr91 = '40091' 1.00130E-04 + Zr92 = '40092' 1.53051E-04 Zr94 = '40094' 1.55104E-04 + Zr96 = '40096' 2.49880E-05 + MIX 4 579.9 O16 = '8016' 2.87335E-4 + Cr52 = '24052' 7.07291E-5 Fe56 = '26056' 1.38298E-4 + Zr90 = '40090' 2.01774E-02 Zr91 = '40091' 4.40020E-03 2 + Zr92 = '40092' 6.72580E-03 Zr94 = '40094' 6.81600E-03 + Zr96 = '40096' 1.09809E-03 + MIX 5 579.9 + H1H2O = '1001' 4.71346E-2 O16H2O = '8016' 2.35673E-2 + Cr52 = '24052' 2.24991E-5 + Fe56 = '26056' 2.09013E-5 Ni58 = '28058' 5.32188E-5 + Mo92 = '42092' 2.77278E-07 Mo95 = '42095' 1.73892E-07 + Mo96 = '42096' 3.00446E-07 Mo97 = '42097' 1.81369E-07 + Mo98 = '42098' 4.59801E-07 Mo100 = '42100' 1.84435E-07 + Al27 = '13027' 1.10277E-6 + Mn55 = '25055' 1.94976E-7 + B10 = '5010' 4.68840E-06 B11 = '5011' 1.88714E-05 + Zr90 = '40090' 2.15252E-04 Zr91 = '40091' 4.69413E-05 + Zr92 = '40092' 7.17508E-05 Zr94 = '40094' 7.27131E-05 + Zr96 = '40096' 1.17144E-05 + MIX 6 579.9 + H1H2O = '1001' 4.71676E-2 O16H2O = '8016' 2.35838E-2 + Cr52 = '24052' 2.11122E-5 + Fe56 = '26056' 1.96130E-5 Ni58 = '28058' 4.99383E-5 + Mo92 = '42092' 2.60187E-07 Mo95 = '42095' 1.63174E-07 + Mo96 = '42096' 2.81927E-07 Mo97 = '42097' 1.70190E-07 + Mo98 = '42098' 4.31460E-07 Mo100 = '42100' 1.73067E-07 + Al27 = '13027' 1.03479E-6 + Mn55 = '25055' 1.82957E-7 + B10 = '5010' 4.69148E-06 B11 = '5011' 1.88838E-05 + Zr90 = '40090' 2.01984E-04 Zr91 = '40091' 4.40478E-05 + Zr92 = '40092' 6.73280E-05 Zr94 = '40094' 6.82309E-05 + Zr96 = '40096' 1.09923E-05 + MIX 7 579.9 + H1H2O = '1001' 4.72020E-2 O16H2O = '8016' 2.36010E-2 + Cr52 = '24052' 1.96591E-5 + Fe56 = '26056' 1.82630E-5 Ni58 = '28058' 4.65011E-5 + Mo92 = '42092' 2.42278E-07 Mo95 = '42095' 1.51943E-07 + Mo96 = '42096' 2.62522E-07 Mo97 = '42097' 1.58476E-07 + Mo98 = '42098' 4.01763E-07 Mo100 = '42100' 1.61155E-07 + Al27 = '13027' 9.63569E-7 + Mn55 = '25055' 1.70365E-7 + B10 = '5010' 4.69469E-06 B11 = '5011' 1.88967E-05 + Zr90 = '40090' 1.88082E-04 Zr91 = '40091' 4.10161E-05 2 + Zr92 = '40092' 6.26939E-05 Zr94 = '40094' 6.35347E-05 + Zr96 = '40096' 1.02357E-05 + MIX 8 933.6 O16 = '8016' 4.49355E-2 + U235 = '92235' 7.39237E-4 1 + U238 = '92238' 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 'TCWE13: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON' + MAXR 400 QUA2 6 3 ; +DISCR2 FILTRK := NXT: ASSMB :: + TITLE 'TCWE13: 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.260388 ; + +OUT := EDI: LIBRARY DISCR1 CALC ASSMB :: + EDIT 3 UPS SAVE MICR RES MERGE CELL COND 4.0 + ; +ECHO "test TCWE13 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE14.c2m b/Dragon/data/twimsE_proc/TCWE14.c2m new file mode 100644 index 0000000..8f626ec --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE14.c2m @@ -0,0 +1,121 @@ +*---- +* TEST CASE TCWE14 +* CANDU-6 ANNULAR CELL +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +PROCEDURE TCWE05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWE05Lib :: <> ; +*---- +* 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 'TCWE14: 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 'TCWE14: 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.133684 ; +EDITION := EDI: LIBRARY VOLMATF FLUX :: + COND 4.0 MERGE CELL MICR 1 Xe135 SAVE ON 'fuel' ; +INTLINF INTLINS := DELETE: INTLINF INTLINS ; +ECHO "test TCWE14 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE15.c2m b/Dragon/data/twimsE_proc/TCWE15.c2m new file mode 100644 index 0000000..a2e54f4 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE15.c2m @@ -0,0 +1,159 @@ +*---- +* TEST CASE TCWE05 +* CANDU-6 ANNULAR CELL +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +PROCEDURE TCWE05Lib ; +PROCEDURE assertS ; +INTEGER iedit := 1 ; +LIBRARY := TCWE05Lib :: <> ; +*---- +* 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 'TCWE05: CANDU-6 CARTESIAN POWER= 31.971 FUEL TEMP= 941.29' + ALLG BATCH 100 MERGMIX + EDIT 0 TRAK TISO 5 10.0 ; +LIBRARY := SHI: LIBRARY VOLMATS INTLINS :: + EDIT 0 ; +VOLMATF INTLINF := NXT: CANDU6F :: + TITLE 'TCWE05: 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.124479 ; +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 <> <> DAY POWR <> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + NOEX DEPL <> <> DAY POWR <> ; + 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.9688265 ; +*---- +* 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 TCWE15 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE17.c2m b/Dragon/data/twimsE_proc/TCWE17.c2m new file mode 100644 index 0000000..1d0d955 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE17.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 TCWE17Lib ; +PROCEDURE assertS ; +INTEGER iedit := 1 ; +MicLib := TCWE17Lib :: <> ; +*---- +* 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.039269 ; +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 TCWE17 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE17Lib.c2m b/Dragon/data/twimsE_proc/TCWE17Lib.c2m new file mode 100644 index 0000000..e085080 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE17Lib.c2m @@ -0,0 +1,97 @@ +*---- +* Create Library for test CASE TCWE17.x2m +* Calling : +* LIBRARY := TCWE17Lib.x2m :: iedit ; +* with : +* LIBRARY = Linked list containing the result of LIB: for +* TCWE17.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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +* All materials are duplicated for left and right cell identification +* during homogenization +*---- +LIBRARY := LIB: :: + EDIT <> + NMIX 21 CTRA WIMS + DEPL LIB: WIMSE FIL: WIMS-E + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 560.66 0.81212 O16 = '8016' 7.99449E-1 + D2D2O = '1002' 1.99768E-1 H1H2O = '1001' 7.83774E-4 + MIX 2 560.66 6.57 Nb93 = '41093' 2.50000 + B10 = '5010' 4.17900E-05 B11 = '5011' 1.68210E-04 + Zr90 = '40090' 5.01638E+01 Zr91 = '40091' 1.09395E+01 2 + Zr92 = '40092' 1.67213E+01 Zr94 = '40094' 1.69455E+01 + Zr96 = '40096' 2.73000E+00 + MIX 3 345.66 0.0014 He4 = '2004' 1.00000E+2 + MIX 4 345.66 6.44 Fe56 = '26056' 1.60000E-1 + Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1 + B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04 + Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2 + Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01 + Zr96 = '40096' 2.79188E+00 + MIX 5 345.66 1.082885 O16 = '8016' 7.98895E-1 + D2D2O = '1002' 2.01016E-1 H1H2O = '1001' 8.96000E-5 + MIX 6 941.29 10.4375010 O16 = '8016' 1.18473E+1 + Xe135 = '54135' 0.0 + U235 = '92235' 6.27118E-1 1 + U238 = '92238' 8.75256E+1 1 + U236 = '92236' 0.0 1 + Pu239 = '94239' 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 = '26056' 1.60000E-1 + Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1 + B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04 + Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2 + Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01 + Zr96 = '40096' 2.79188E+00 + MIX 11 560.66 0.81212 O16 = '8016' 7.99449E-1 + D2D2O = '1002' 1.99768E-1 H1H2O = '1001' 7.83774E-4 + MIX 12 560.66 6.57 Nb93 = '41093' 2.50000 + B10 = '5010' 4.17900E-05 B11 = '5011' 1.68210E-04 + Zr90 = '40090' 5.01638E+01 Zr91 = '40091' 1.09395E+01 2 + Zr92 = '40092' 1.67213E+01 Zr94 = '40094' 1.69455E+01 + Zr96 = '40096' 2.73000E+00 + MIX 13 345.66 0.0014 He4 = '2004' 1.00000E+2 + MIX 14 345.66 6.44 Fe56 = '26056' 1.60000E-1 + Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1 + B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04 + Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2 + Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01 + Zr96 = '40096' 2.79188E+00 + MIX 15 345.66 1.082885 O16 = '8016' 7.98895E-1 + D2D2O = '1002' 2.01016E-1 H1H2O = '1001' 8.96000E-5 + MIX 16 941.29 10.4375010 O16 = '8016' 1.18473E+1 + Xe135 = '54135' 0.0 + U235 = '92235' 6.27118E-1 1 + U238 = '92238' 8.75256E+1 1 + U236 = '92236' 0.0 1 + Pu239 = '94239' 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 = '26056' 1.60000E-1 + Ni58 = '28058' 6.00000E-2 Cr52 = '24052' 1.10000E-1 + B10 = '5010' 6.16900E-05 B11 = '5011' 2.48310E-04 + Zr90 = '40090' 5.13008E+01 Zr91 = '40091' 1.11875E+01 2 + Zr92 = '40092' 1.71003E+01 Zr94 = '40094' 1.73296E+01 + Zr96 = '40096' 2.79188E+00 + MIX 21 345.66 Fe56 = '26056' 6.19027E-2 + Cr52 = '24052' 1.56659E-3 Ni58 = '28058' 6.83337E-3 + Si29 = '14029' 7.79072E-4 C12 = '6000' 1.46552E-4 + Mn55 = '25055' 1.25431E-3 + ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE31.c2m b/Dragon/data/twimsE_proc/TCWE31.c2m new file mode 100644 index 0000000..084b688 --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE31.c2m @@ -0,0 +1,172 @@ +*---- +* TEST CASE TCWE31 +* CANDU-6 ANNULAR CELL +* WIMS-E 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 WIMS-E format WIMSE +* Microscopic cross sections from file WIMS-E format WIMSE +*---- +PROCEDURE TCWE05Lib ; +INTEGER iedit := 1 ; +LIBRARY := TCWE05Lib :: <> ; +*---- +* 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 'TCWE31: 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 'TCWE31: 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.124889 ; +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 'TCWE31: 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.124861 ; +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 <> <> DAY POWR <> ; + ELSE + BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF :: + EDIT 3 NOEX DEPL <> <> DAY POWR <> ; + 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.079416 ; +res := EDITION ; +EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE: + EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ; +ECHO "test TCWE31 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/twimsE_proc/TCWE32.c2m b/Dragon/data/twimsE_proc/TCWE32.c2m new file mode 100644 index 0000000..68096ab --- /dev/null +++ b/Dragon/data/twimsE_proc/TCWE32.c2m @@ -0,0 +1,97 @@ +*---- +* TEST CASE TCWE32 +* 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 WIMSE +*---- +LIBRARY := LIB: :: + EDIT 1 + NMIX 4 + CTRA WIMS + MIXS LIB: WIMSE FIL: WIMS-E + MIX 1 579.9 + H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2 + BNat = '1011' 2.38103E-5 + MIX 2 609.0 O16 = '8016' 3.06711E-4 + Cr52 = '24052' 7.54987E-5 Fe56 = '26056' 1.47624E-4 + Zr90 = '40090' 2.15381E-02 Zr91 = '40091' 4.69693E-03 2 + Zr92 = '40092' 7.17935E-03 Zr94 = '40094' 7.27563E-03 + Zr96 = '40096' 1.17214E-03 + MIX 3 966.0 O16 = '8016' 4.49355E-2 + U235 = '92235' 7.39237E-4 1 + U238 = '92238' 2.17285E-2 1 + MIX 4 579.9 + H1H2O = '1001' 4.76690E-2 O16H2O = '8016' 2.38345E-2 + B10 = '5010' 4.73825E-06 B11 = '5011' 1.90721E-05 + ; +*---- +* 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.360705 ; +TRACK TRACK_LGT GEOMF := DELETE: TRACK TRACK_LGT GEOMF ; +PIJMAT FLUX := DELETE: PIJMAT FLUX ; +LIBRARY := DELETE: LIBRARY ; +ECHO "test TCWE32 completed" ; +END: ; +QUIT "LIST" . -- cgit v1.2.3