diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Donjon/data/Test_microdepl_NCR.x2m | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Donjon/data/Test_microdepl_NCR.x2m')
| -rw-r--r-- | Donjon/data/Test_microdepl_NCR.x2m | 347 |
1 files changed, 347 insertions, 0 deletions
diff --git a/Donjon/data/Test_microdepl_NCR.x2m b/Donjon/data/Test_microdepl_NCR.x2m new file mode 100644 index 0000000..acb80db --- /dev/null +++ b/Donjon/data/Test_microdepl_NCR.x2m @@ -0,0 +1,347 @@ +*---- +* Name : Test_microdepl_NCR.x2m +* Type : DONJON CLE2000 procedure +* Description : Test microdepletion capabilities with a simplified +* Xe-Sm depletion chain (8 particularized isotopes) +* 'Xe135' 'I135' 'Nd147' 'Pm147' 'Pm148' 'Pm148m' 'Pm149' +* 'Sm149' +*---- +* STRUCTURES, MODULES and PROCEDURES definitions +*---- +MODULE END: UTL: DELETE: NCR: GEO: USPLIT: TRIVAT: + RESINI: MAC: MACINI: TRIVAA: FLUD: FLPOW: EVO: ABORT: ; +LINKED_LIST compo Micro MacroF GeoCoeur GeoRes Matex Track Fmap + MacroRefl MacroT System Flux Power EvolFlux Burn + Geoevo Geoevo2 Trackevo Matex2 ; +SEQ_ASCII _compo :: FILE '_Compo_1BU.txt' ; +PROCEDURE assertS ; + +REAL dx := 21.61 ; +REAL dz := 20.0 ; +INTEGER Splitx Splitz := 2 1 ; + +REAL x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 x13 x14 + x15 x16 x17 x18 x19 := + 1. dx * 2. dx * 3. dx * 4. dx * 5. dx * 6. dx * + 7. dx * 8. dx * 9. dx * 10. dx * 11. dx * 12. dx * + 13. dx * 14. dx * 15. dx * 16. dx * 17. dx * 18. dx * 19. dx * ; + +REAL z1 z2 z3 z4 z5 z6 z7 z8 z9 z10 z11 z12 z13 z14 + z15 z16 z17 z18 z19 z20 z21 z22 z23 := + 1. dz * 2. dz * 3. dz * 4. dz * 5. dz * 6. dz * + 7. dz * 8. dz * 9. dz * 10. dz * 11. dz * 12. dz * + 13. dz * 14. dz * 15. dz * 16. dz * 17. dz * 18. dz * + 19. dz * 20. dz * 21. dz * 22. dz * 23. dz * ; + +GeoCoeur := GEO: :: CAR3D 19 19 23 + EDIT 0 + X+ VOID X- VOID + Y+ VOID Y- VOID + Z+ VOID Z- VOID + MIX + +PLANE 1 +0 0 0 0 0 2 2 2 2 2 2 2 2 2 0 0 0 0 0 +0 0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 0 +0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 +0 0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 0 +0 0 0 0 0 2 2 2 2 2 2 2 2 2 0 0 0 0 0 + +PLANE 2 +0 0 0 0 0 2 2 2 2 2 2 2 2 2 0 0 0 0 0 +0 0 0 2 2 2 1 1 1 1 1 1 1 2 2 2 0 0 0 +0 0 2 2 1 1 1 1 1 1 1 1 1 1 1 2 2 0 0 +0 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 0 +0 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 0 +2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 +2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 +2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 +2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 +2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 +2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 +2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 +2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 +2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 +0 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 0 +0 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 0 +0 0 2 2 1 1 1 1 1 1 1 1 1 1 1 2 2 0 0 +0 0 0 2 2 2 1 1 1 1 1 1 1 2 2 2 0 0 0 +0 0 0 0 0 2 2 2 2 2 2 2 2 2 0 0 0 0 0 + +PLANE 3 SAME 2 +PLANE 4 SAME 2 +PLANE 5 SAME 2 +PLANE 6 SAME 2 +PLANE 7 SAME 2 +PLANE 8 SAME 2 +PLANE 9 SAME 2 +PLANE 10 SAME 2 +PLANE 11 SAME 2 +PLANE 12 SAME 2 +PLANE 13 SAME 2 +PLANE 14 SAME 2 +PLANE 15 SAME 2 +PLANE 16 SAME 2 +PLANE 17 SAME 2 +PLANE 18 SAME 2 +PLANE 19 SAME 2 +PLANE 20 SAME 2 +PLANE 21 SAME 2 +PLANE 22 SAME 2 + +PLANE 23 +0 0 0 0 0 2 2 2 2 2 2 2 2 2 0 0 0 0 0 +0 0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 0 +0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 +0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 +0 0 0 2 2 2 2 2 2 2 2 2 2 2 2 2 0 0 0 +0 0 0 0 0 2 2 2 2 2 2 2 2 2 0 0 0 0 0 + + MESHX 0. <<x1>> <<x2>> <<x3>> <<x4>> <<x5>> <<x6>> <<x7>> + <<x8>> <<x9>> <<x10>> <<x11>> <<x12>> <<x13>> <<x14>> + <<x15>> <<x16>> <<x17>> <<x18>> <<x19>> + + MESHY 0. <<x1>> <<x2>> <<x3>> <<x4>> <<x5>> <<x6>> <<x7>> + <<x8>> <<x9>> <<x10>> <<x11>> <<x12>> <<x13>> <<x14>> + <<x15>> <<x16>> <<x17>> <<x18>> <<x19>> + + MESHZ 0.0 <<z1>> <<z2>> <<z3>> <<z4>> <<z5>> <<z6>> <<z7>> + <<z8>> <<z9>> <<z10>> <<z11>> <<z12>> <<z13>> <<z14>> + <<z15>> <<z16>> <<z17>> <<z18>> <<z19>> <<z20>> <<z21>> + <<z22>> <<z23>> + + SPLITX <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + SPLITY <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + <<Splitx>> <<Splitx>> <<Splitx>> <<Splitx>> + + SPLITZ <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> + <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> + <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> + <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> <<Splitz>> + <<Splitz>> <<Splitz>> <<Splitz>> +; + +GeoRes Matex := USPLIT: GeoCoeur :: NGRP 2 MAXR 1000000 + NREFL 1 RMIX 2 + NFUEL 1 FMIX 1 +; + +*-- +* Definition de la map de combustible et zones de combustion +*-- +Fmap Matex := RESINI: Matex :: + ::: GEO: CAR3D 19 19 3 + EDIT 0 + X- VOID X+ VOID + Y- VOID Y+ VOID + Z- VOID Z+ VOID +MIX +PLANE 1 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +PLANE 2 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +0 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0 0 0 +0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 +0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 +0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 +0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 +0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 +0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 +0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 +0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 +0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 +0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 +0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 +0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 +0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 +0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 +0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 +0 0 0 0 0 0 1 1 1 1 1 1 1 0 0 0 0 0 0 +0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +PLANE 3 SAME 1 +MESHX 0.0 <<x1>> <<x2>> <<x3>> <<x4>> <<x5>> <<x6>> <<x7>> + <<x8>> <<x9>> <<x10>> <<x11>> <<x12>> <<x13>> <<x14>> + <<x15>> <<x16>> <<x17>> <<x18>> <<x19>> +MESHY 0.0 <<x1>> <<x2>> <<x3>> <<x4>> <<x5>> <<x6>> <<x7>> + <<x8>> <<x9>> <<x10>> <<x11>> <<x12>> <<x13>> <<x14>> + <<x15>> <<x16>> <<x17>> <<x18>> <<x19>> +MESHZ 0.0 <<z1>> <<z22>> <<z23>> ; +NXNAME '-' '01' '02' '03' '04' '05' '06' '07' '08' '09' + '10' '11' '12' '13' '14' '15' '16' '17' '-' +NYNAME '-' 'A' 'B' 'C' 'D' 'E' 'F' 'G' 'H' 'J' 'K' 'L' + 'M' 'N' 'O' 'P' 'Q' 'R' '-' +NCOMB 241 B-ZONE + 1 2 3 4 5 6 7 + 8 9 10 11 12 13 14 15 16 17 18 + 19 20 21 22 23 24 25 26 27 28 29 30 31 + 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 + 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 + 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 + 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 + 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 + 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 + 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 + 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 + 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 + 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 + 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 + 211 212 213 214 215 216 217 218 219 220 221 222 223 + 224 225 226 227 228 229 230 231 232 233 234 + 235 236 237 238 239 240 241 +; + +Fmap := RESINI: Fmap :: EDIT 0 + + BTYPE INST-BURN + INST-BVAL CHAN + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 + 0.0 + ; + +MacroRefl := MAC: :: + EDIT 3 NGRO 2 NMIX 2 NIFI 0 + ANIS 1 CTRA NONE ! ENER 0.625 + READ INPUT + MIX 2 ! Reflecteur radial + TOTAL 7.161029E-03 1.066538E-01 + DIFF 0.851754657 0.686485503 + SCAT 2 2 (*2->1*) 0.0 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) 3.671925E-03 +; + +*--- +* Recover the fuel-map geometry to prepare the depletion calculation +*--- +Geoevo := Fmap :: STEP UP 'GEOMAP' ; +Geoevo2 Matex2 := USPLIT: Geoevo :: NGRP 2 MAXR 2000 + NFUEL 1 FMIX 1 ; +Trackevo := TRIVAT: Geoevo2 :: MAXR 2000 MCFD 1 ; + +*---- +* Convert the multicompo from ascii to linked list +*---- +compo := _compo ; + +*---- +* Multicompo interpolation +*---- +Micro := NCR: compo Fmap :: + EDIT 1 MICRO LINEAR + NMIX 241 TABLE compo 'Ass' 'BURN' + MIX 1 FROM 1 + SET 'TCOM' 6.5E+02 ! Celsius + SET 'DMOD' 0.7 ! g/cm3 + SET 'BARR' 1.0 ! Sans grappe + SET 'CBOR' 2400.0 ! ppm + MICRO ALL + ENDMIX + ; + +MacroF := Micro :: STEP UP 'MACROLIB' ; + +MacroT Matex := MACINI: Matex MacroRefl MacroF ; +Track := TRIVAT: GeoRes :: + MAXR 1000000 DUAL 2 3 ; +System := TRIVAA: MacroT Track ; + +UTL: Trackevo :: IMPR STATE-VECTOR * ; +UTL: MacroF :: IMPR STATE-VECTOR * ; +UTL: MacroT :: IMPR STATE-VECTOR * ; + +Flux := FLUD: System Track :: + EDIT 2 ADI 2 EXTE 1.0E-4 ; + +Power EvolFlux := FLPOW: Fmap Flux Track Matex :: + EDIT 2 + PTOT 4250.0 BUND + PRINT DISTR POWER ; + +Burn Micro := EVO: Micro Power :: + EDIT 1 RUNG FLUX_POW PIFI DEPL 4.483 DAY KEEP ; + +MacroF MacroT System Flux Power EvolFlux := DELETE: +MacroF MacroT System Flux Power EvolFlux ; + +*---- +* Xenon and Samarium are saturated +*---- +MacroF := Micro :: STEP UP 'MACROLIB' ; +MacroT Matex := MACINI: Matex MacroRefl MacroF ; +System := TRIVAA: MacroT Track ; +Flux := FLUD: System Track :: + EDIT 2 ADI 2 EXTE 1.0E-4 ; + +assertS Flux :: 'K-EFFECTIVE' 1 1.235925 ; + +ECHO "test Test_microdepl_NCR.x2m completed" ; +END: ; +QUIT "LIST" . |