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/proc_LZC/MicWUC.c2m | 189 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 189 insertions(+) create mode 100644 Dragon/data/proc_LZC/MicWUC.c2m (limited to 'Dragon/data/proc_LZC/MicWUC.c2m') diff --git a/Dragon/data/proc_LZC/MicWUC.c2m b/Dragon/data/proc_LZC/MicWUC.c2m new file mode 100644 index 0000000..d24bbec --- /dev/null +++ b/Dragon/data/proc_LZC/MicWUC.c2m @@ -0,0 +1,189 @@ +*DECK MicWUC.c2m +*---------- +* Name : MicWUC.c2m +* Type : DRAGON procedure +* Use : Create Gentilly-2 fuel microlib +* with cross sections from the WLUP based 69 +* groups (LibName = IAEA69). +* Author : G. Marleau +* Date : 2006/03/01 +* +* Procedure description: +* MicroLib := MicWUC :: +* <> ; +* +* Procedure parameters: +* MicroLib : XSM_FILE containing the new microlib where: +* Mix 1 = coolant at full density +* Mix 2 = pressure tube +* Mix 3 = helium for gap and LZC +* Mix 4 = Zr for calandrial tube and adjuster +* Mix 5 = moderator +* Mix 6 = left central fuel +* Mix 7 = left 6 pin fuel +* Mix 8 = left 12 pin fuel +* Mix 9 = left outer fuel +* Mix 10 = right central fuel +* Mix 11 = right 6 pin fuel +* Mix 12 = right 12 pin fuel +* Mix 13 = right outer fuel +* Mix 14 = fuel sheating +* (see GeoLZC.c2m for mixture location in +* geometry) +* PrtLev : Print Level +* 0 -> Minimum print +* 1 -> Standard print +* 10 -> Debug +* +*---------- +* Define procedure parameters and options and read +* the procedure input data +*---- +PARAMETER MicroLib :: EDIT 0 ::: XSM_FILE MicroLib ; ; +INTEGER PrtLev ; +:: >>PrtLev<< ; +*---- +* Local modules used in this procedure +*---- +MODULE INFO: LIB: DELETE: ; +*---------- +* Define local fuel, coolant and moderator properties +* Fuel enrichment (%) = 0.71140; T = 941.26 K; RHO =10.4375 g/cc; +* Coolant D2 AT (%) = 99.222; T = 560.66 K; +* Moderator D2 AT (%) = 99.911; T = 345.66 K. +*---------- +REAL + Ftemp Fenri Fdens + Ctemp Cprty + Mtemp Mprty + Wprty := + 941.26 0.7114 10.4375 + 560.66 99.30 + 345.66 99.92 + 0.0 ; +REAL + FwgtU8 FwgtU5 FwgtO16 + CwgtH1 CwgtD2 CwgtO16 Cdens + MwgtH1 MwgtD2 MwgtO16 Mdens + WwgtH1 WwgtD2 WwgtO16 Wdens ; +REAL + BNZrNb BNZr BNZrII := + 0.00021 0.00031 0.00005962 ; +REAL + B10ZrNb B10Zr B10ZrII + B11ZrNb B11Zr B11ZrII := + 0.199 BNZrNb * 0.199 BNZr * 0.199 BNZrII * + 0.801 BNZrNb * 0.801 BNZr * 0.801 BNZrII * ; +*---------- +* Library and isotope names +*---------- +STRING LibName + IsoU5 IsoU8 IsoU6 IsoPu9 IsoO16 + IsoH1 IsoD2 IsoB10 IsoB11 IsoBn + IsoZr IsoFe IsoFeA IsoCr IsoNi + IsoNb IsoC IsoSi IsoMn IsoCu + IsoMo IsoAl IsoHe4 IsoN IsoTi + IsoCd IsoCo59 := + "IAEA69" + "2235" "8238" "236" "6239" "6016" + "3001" "3002" "10" "11" "1011" + "91" "2056" "2056" "52" "58" + "93" "2012" "29" "55" "3063" + "96" "27" "4" "14" "48" + "2113" "1059" ; +*---------- +* End region information +* endzrv = Volume of Zr in endregion +* endcoolv = Volume of coolant in endregion +* endgasv = Volume of gas (void) in endregion +* endtotv = Total volume of endregion +* endzrf = Volume fraction of Zr in endregion +* endcoolf = Volume fraction of coolant in endregion +* endgasf = Volume fraction of gas (void) in endregion +*---------- +REAL endzrv endcoolv endgasv := 42.40 77.43 8.65 ; +REAL endtotv := endzrv endcoolv endgasv + + ; +REAL endzrf endcoolf endgasf := + endzrv endtotv / endcoolv endtotv / endgasv endtotv / ; +*---------- +* Guide tube perforation +*---------- +REAL PerfM1 PerfM2 PerfM3 := 0.3566 0.02 0.005 ; +REAL PerfZ1 PerfZ2 PerfZ3 := + 1.0 PerfM1 - 1.0 PerfM2 - 1.0 PerfM3 - ; +*---------- +* COMPUTE WGT% FOR FUEL +* COMPUTE DENSITY AND WGT% FOR COOLANT AND MODERATOR +*---------- +INFO: :: + ENR: <> WGT% + LIB: WIMSD4 FIL: <> + ISO: 3 <> <> <> + CALC WGT% UO2 >>FwgtU5<< >>FwgtU8<< >>FwgtO16<< ; +INFO: :: + TMP: <> K + PUR: <> WGT% + CALC DENS WATER >>Cdens<< + LIB: WIMSD4 FIL: <> + ISO: 3 <> <> <> + CALC WGT% D2O >>CwgtH1<< >>CwgtD2<< >>CwgtO16<< ; +INFO: :: + TMP: <> K + PUR: <> WGT% + CALC DENS WATER >>Mdens<< + LIB: WIMSD4 FIL: <> + ISO: 3 <> <> <> + CALC WGT% D2O >>MwgtH1<< >>MwgtD2<< >>MwgtO16<< ; +INFO: :: + TMP: <> K + PUR: <> WGT% + CALC DENS WATER >>Wdens<< + LIB: WIMSD4 FIL: <> + ISO: 3 <> <> <> + CALC WGT% D2O >>WwgtH1<< >>WwgtD2<< >>WwgtO16<< ; +*---------- +* Generate Microscopic cross section library +*---------- +MicroLib := LIB: :: + EDIT <> + NMIX 14 CTRA WIMS + DEPL LIB: WIMSD4 FIL: <> + MIXS LIB: WIMSD4 FIL: <> + MIX 1 <> <> + H1 = <> <> D2 = <> <> + O16 = <> <> + MIX 2 <> 6.57 + ZR = <> 97.5 NB = <> 2.5 +* BN = <> <> + B10 = <> <> B11 = <> <> + MIX 3 <> 0.0014 + HE4 = <> 100.0 + MIX 4 <> 6.44 + ZR = <> 99.71 FE = <> 0.16 + CR = <> 0.11 NI = <> 0.06 +* BN = <> <> + B10 = <> <> B11 = <> <> + MIX 5 <> <> + H1 = <> <> D2 = <> <> + O16 = <> <> + MIX 6 <> <> + O16 = <> <> U5 = <> <> 1 + U8 = <> <> 1 U6 = <> 0.0 1 + PU9 = <> 0.0 1 + MIX 7 COMB 6 100.0 + MIX 8 COMB 6 100.0 + MIX 9 COMB 6 100.0 + MIX 10 <> <> + O16 = <> <> U5 = <> <> 2 + U8 = <> <> 2 U6 = <> 0.0 2 + PU9 = <> 0.0 2 + MIX 11 COMB 10 100.0 + MIX 12 COMB 10 100.0 + MIX 13 COMB 10 100.0 + MIX 14 <> 6.44 + ZR = <> 99.71 FE = <> 0.16 + CR = <> 0.11 NI = <> 0.06 +* BN = <> <> + B10 = <> <> B11 = <> <> + ; +QUIT . -- cgit v1.2.3