*DECK CreaMix.c2m *====================================================== * Name: CreaMix.c2m * Type: DRAGON procedure * Use: Create isotopic mixtures depending on * assembly type. * Author: A. Bruneau 04.2015, 09.2015 * [CS-SI | IRSN/PSN-EXP/SNC/LNR] * Mod: R. Nguyen Van Ho 11.2016 * [URANUS | IRSN/PSN-EXP/SNC/LNR] * Ajout selecteur donnees nucleaires *------------------------------------------------------ * Input Required: * > Molar Masses * > Metal Mixes * > Gestion (STR) AssType (STR) Teneur_I (INT) ThermaExpans (INT) * Output Provided: * < Library * < Moderator concentrations * < Power Conversion Factor (PowMWT, REAL) *====================================================== PARAMETER Library ConcMode FeNat CrNat SnNat NiNat BNat ZrNat2 TiNat2 SiNat2 KNat2 AgNat CdNat SNat InNat Zirc4 Incon SS304 M5 AIC GdNat Gd2O3 BP MolarMasses :: ::: LINKED_LIST FeNat CrNat SnNat NiNat BNat ZrNat2 TiNat2 SiNat2 KNat2 AgNat CdNat SNat InNat Zirc4 Incon SS304 M5 AIC GdNat Gd2O3 BP MolarMasses ; ; PROCEDURE AdaptIso ThermaExpans DeplChain ; MODULE LIB: GREP: ABORT: UTL: END: ; LINKED_LIST ListeIsot ListeAutop ; LINKED_LIST ZrNat TiNat SiNat KNat ; ZrNat := ZrNat2 ; TiNat := TiNat2 ; SiNat := SiNat2 ; KNat := KNat2 ; *-------------------------------------------------------------------- * Constant parameters *-------------------------------------------------------------------- * Avogadro de 1974. Ref : History and progress in the accurate * determination of the Avogadro constant, P. Becker, 2001. * L'impact de l'utilisation d'un Avogadro plus 'recent' est environ * d'un pcm. REAL avogadro := 6.022094E+23 ; REAL avo_conc := avogadro 1.E-24 * ; *-------------------------------------------------------------------- * General assembly parameters *-------------------------------------------------------------------- REAL cbor dmod tmod tgain tcom tclad ; STRING Gestion AssType ; STRING NuclData ; ! Nom de la bibliotheque de XS micro JEFF ou ENDF STRING LibType ; ! Type de bibliotheque (ENDFBVII0...) STRING LibGenre ; STRING Model := "TIH" ; ! Model for thermal expansion INTEGER LibId ; ! Type de bibliotheque de XS micro INTEGER Teneur_I ; ! in 100 * % (ex: 370 for 3.7%) REAL Main_enri ; ! = Teneur_I / 100 [%] REAL PowMWth ; ! Core Thermal Power in MW REAL PowMWT ; ! Power Conversion Factor in MW/T REAL Hactiv ; ! Active height in cm REAL NbAss ; ! Nb of assemblies REAL ReNorm ; ! Re-normalization factor (for isotopic ratios) REAL H2OTube GrPreTube ; ! water & Grid proportions for tube REAL H2OCntr GrPreCntr ; ! water & Grid prop. for central cells REAL H2OSide GrPreSide ; ! water & Grid prop. for side cells REAL H2OCorn GrPreCorn ; ! water & Grid prop. for corner cells INTEGER ThermExp ; STRING CTRA ; STRING Grid ; INTEGER iU5 iU8 ; ! SHIB inrs for U INTEGER iPu9 iPu0 iPu1 iPu2 ; ! SHIB inrs for Pu INTEGER iAg7 iAg9 iIn ; ! SHIB inrs for Ag & In (AIC) INTEGER iZrG iZrT ; ! SHIB inrs Zr in Clad & Tube INTEGER iGd4 iGd5 iGd6 iGd7 iGd8 iGd0 ; ! SHIB inrs Gd INTEGER iSi8 iSi9 iSi0 ; INTEGER iAl7 iK9 iK0 iK1 iNa3 ; ! SHIB inrs for Pyrex *-------------------------------------------------------------------- * Specific assembly parameters : UOX *-------------------------------------------------------------------- REAL enriU ; ! U-235 enrichment in regular fuel pin REAL MUO2 ; ! UO2-Molar Mass REAL MUPuAm ; ! Average 'U' Molar Mass REAL dU ; ! Heavy metal density [g/cm**3] *--- 'isotopic enrichment' (final isotopic fractions) --------------- REAL U4enri U5enri U6enri U8enri Pu8enri Pu9enri Pu0enri Pu1enri Pu2enri Am1enri ; *--- fuel concentrations in at/barn/cm ------------------------------ REAL U_tot ; ! Total U concentration [at/barn/cm] REAL U4fuel U5fuel U6fuel U8fuel Pu8fuel Pu9fuel Pu0fuel Pu1fuel Pu2fuel Am1fuel O_fuel := 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 ; *-------------------------------------------------------------------- * Specific assembly parameters : MOX *-------------------------------------------------------------------- REAL enriPu ; ! Fuel-average Pu vector enrichment REAL f_corn ; ! Pu enrichment fraction for corner pins REAL f_side ; ! Pu enrichment fraction for side pins REAL f_cntr ; ! Pu enrichment fraction for all remaining pins REAL MUtot ; ! Average 'U' Molar Mass REAL MPutot ; ! Average 'Pu' Molar Mass REAL M_corn ; ! Average 'U+Pu' Molar Mass for corner pins REAL M_side ; ! Average 'U+Pu' Molar Mass for side pins REAL M_cntr ; ! Average 'U+Pu' Molar Mass for center pins REAL MMOXcorn ; ! corner MOX molar mass REAL MMOXside ; ! side MOX molar mass REAL MMOXcntr ; ! center MOX molar mass REAL MOX_corn ; ! corner MOX concentration [molec./barn/cm] REAL MOX_side ; ! side MOX concentration [molec./barn/cm] REAL MOX_cntr ; ! center MOX concentration [molec./barn/cm] REAL f_cornPin ; ! corner pins fraction REAL f_sidePin ; ! side pins fraction REAL f_cntrPin ; ! center pins fraction REAL dUcorn ; ! corner heavy metal density [g/cm**3] REAL dUside ; ! side heavy metal density [g/cm**3] REAL dUcntr ; ! center heavy metal density [g/cm**3] *--- 'isotopic enrichment' (final isotopic fractions) --------------- * isotopic ratios & enrichments for corner pins REAL Ue_corn ; ! Total U fraction in corner pins REAL U4e_corn U5e_corn U6e_corn U8e_corn ; REAL Pue_corn ; ! Total Pu fraction in corner pins REAL Pu8e_corn Pu9e_corn Pu0e_corn Pu1e_corn Pu2e_corn Am1e_corn ; * isotopic ratios & enrichments for side pins REAL Ue_side ; ! Total U fraction in side pins REAL U4e_side U5e_side U6e_side U8e_side ; REAL Pue_side ; ! Total Pu fraction in side pins REAL Pu8e_side Pu9e_side Pu0e_side Pu1e_side Pu2e_side Am1e_side ; * isotopic ratios & enrichments for center pins REAL Ue_cntr ; ! Total U fraction in center pins REAL U4e_cntr U5e_cntr U6e_cntr U8e_cntr ; REAL Pue_cntr ; ! Total Pu fraction in center pins REAL Pu8e_cntr Pu9e_cntr Pu0e_cntr Pu1e_cntr Pu2e_cntr Am1e_cntr ; *--- corner fuel concentrations in at/barn/cm ----------------------- REAL UPu_corn_tot ; ! Total corner U+Pu concentration [at/barn/cm] REAL U4cornfuel U5cornfuel U6cornfuel U8cornfuel Pu8cornfuel Pu9cornfuel Pu0cornfuel Pu1cornfuel Pu2cornfuel Am1cornfuel O_corn_fuel := 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. ; *--- side fuel concentrations in at/barn/cm ------------------------- REAL UPu_side_tot ; ! Total side U+Pu concentration [at/barn/cm] REAL U4sidefuel U5sidefuel U6sidefuel U8sidefuel Pu8sidefuel Pu9sidefuel Pu0sidefuel Pu1sidefuel Pu2sidefuel Am1sidefuel O_side_fuel := 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. ; *--- center fuel concentrations in at/barn/cm ----------------------- REAL UPu_cntr_tot ; ! Total center U+Pu concentration [at/barn/cm] REAL U4cntrfuel U5cntrfuel U6cntrfuel U8cntrfuel Pu8cntrfuel Pu9cntrfuel Pu0cntrfuel Pu1cntrfuel Pu2cntrfuel Am1cntrfuel O_cntr_fuel := 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. ; *-------------------------------------------------------------------- * Specific assembly parameters : UOX + Gd (UGd) *-------------------------------------------------------------------- REAL enriUGd ; ! U-235 enrichment in Gd pins REAL enriGd ; ! Gadolinia (Gd2O3) massic enrichment REAL Gd2O3f ; ! Gd2O3 fraction in UO2 + Gd2O3 pins (= enriGd / 100) REAL UO2Gdf ; ! UO2 fraction un UO2 + Gd2O3 pins (= 1 - Gd2O3f) REAL MGd2O3 ; ! Gd2O3 Molar Mass REAL MUO2Gd ; ! UO2 for UOX+Gd Molar Mass REAL MUPuGd ; ! Average 'U' Molar Mass in Gd pins REAL M_UO2Gd2O3 ; ! Average 'UO2+Gd2O3' Molar Mass REAL f_GadoPin ; ! Gado'ed pins fraction (out of 264 pins) REAL f_NoGdPin ; ! Regular UOX (no Gd) pins fraction (out of 264 pins) REAL dUGado ; ! Gd pins Heavy metal density [g/cm**3] *--- Gd2O3 isotopic fractions --------------------------------------- REAL Gd2O3_G4f Gd2O3_G5f Gd2O3_G6f Gd2O3_G7f Gd2O3_G8f Gd2O3_G0f Gd2O3_Of ; *--- 'isotopic enrichment' (final isotopic fractions) --------------- REAL U4eGd U5eGd U6eGd U8eGd Pu9eGd Pu0eGd Pu1eGd Pu2eGd OeUO2Gd ; *--- fuel concentrations in at/barn/cm ------------------------------ REAL UO2Gd2O3 ; ! Total U+Pu+Gd concentration [at/barn/cm REAL U4Gd_fuel U5Gd_fuel U6Gd_fuel U8Gd_fuel Pu9Gd_fuel Pu0Gd_fuel Pu1Gd_fuel Pu2Gd_fuel Gd154_fuel Gd155_fuel Gd156_fuel Gd157_fuel Gd158_fuel Gd160_fuel O_UGd_fuel := 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. ; * All other SS304 densities based on dSS with dilution coefficients REAL dRodSS ; ! Stainless Steel (SS304) control rod density REAL dCladAIC ; ! Stainless Steel (SS304) Clad dens. (AIC Rods) REAL dCladSS ; ! Stainless Steel (SS304) Clad dens. (Steel Rods) REAL dCladB4C ; ! Stainless Steel (SS304) Clad dens. (B4C Rods) *------------------------------------------------------ * Get NuclData (nom de la bibliotheque) * LibType (nom generique de la bibli, ex. ENDFVII1) * Gestion (TIH, BEAV) * AssType (UOX, MOX, UGd) * Teneur (370, 865, etc...) '370' => 3.70% * Densités d'UO2, de MOX et de combustible gadolinié * Enrichissement de l'U dans UGd * Teneur en Gd2O3 dans les crayons empoisonnés * Nombre de crayons empoisonnés *------------------------------------------------------ :: >>NuclData<< >>LibType<< >>Gestion<< >>AssType<< >>Teneur_I<< >>ThermExp<< >>CTRA<< >>Grid<< ; EVALUATE Main_enri := Teneur_I I_TO_R 100. / ; *------------------------------------------------------ * Molar masses (from LINKED_LIST MolarMasses) *------------------------------------------------------ REAL al27 b10 b10b b11 b11b co59 cr50 cr52 cr53 cr54 fe54 fe56 fe57 fe58 h2o mn55 mo95 nb93 ni58 ni60 ni61 ni62 ni64 o16 pu238 pu239 pu240 pu241 pu242 am241 sinat si28 si29 si30 sn112 sn114 sn115 sn116 sn117 sn118 sn119 sn120 sn122 sn124 tinat ti46 ti47 ti48 ti49 ti50 u234 u235 u236 u238 zrnat zr90 zr91 zr92 zr93 zr94 zr95 zr96 cnat ag107 ag109 cd106 cd108 cd110 cd111 cd112 cd113 cd114 cd116 in115 gd154 gd155 gd156 gd157 gd158 gd160 knat k39 k40 k41 na23 h1 s32 s33 s34 in113 ; *------------------------------------------------------ ! Isotopes non conflictuels GREP: MolarMasses :: GETVAL 'Am241' 1 >>am241<< ; GREP: MolarMasses :: GETVAL 'Al27 ' 1 >>al27<< ; GREP: MolarMasses :: GETVAL 'B10 ' 1 >>b10<< ; GREP: MolarMasses :: GETVAL 'B10B ' 1 >>b10b<< ; GREP: MolarMasses :: GETVAL 'B11 ' 1 >>b11<< ; GREP: MolarMasses :: GETVAL 'B11B ' 1 >>b11b<< ; GREP: MolarMasses :: GETVAL 'Co59 ' 1 >>co59<< ; GREP: MolarMasses :: GETVAL 'Cr50 ' 1 >>cr50<< ; GREP: MolarMasses :: GETVAL 'Cr52 ' 1 >>cr52<< ; GREP: MolarMasses :: GETVAL 'Cr53 ' 1 >>cr53<< ; GREP: MolarMasses :: GETVAL 'Cr54 ' 1 >>cr54<< ; GREP: MolarMasses :: GETVAL 'Fe54 ' 1 >>fe54<< ; GREP: MolarMasses :: GETVAL 'Fe56 ' 1 >>fe56<< ; GREP: MolarMasses :: GETVAL 'Fe57 ' 1 >>fe57<< ; GREP: MolarMasses :: GETVAL 'Fe58 ' 1 >>fe58<< ; GREP: MolarMasses :: GETVAL 'Mn55 ' 1 >>mn55<< ; GREP: MolarMasses :: GETVAL 'Mo95 ' 1 >>mo95<< ; GREP: MolarMasses :: GETVAL 'Nb93 ' 1 >>nb93<< ; GREP: MolarMasses :: GETVAL 'Ni58 ' 1 >>ni58<< ; GREP: MolarMasses :: GETVAL 'Ni60 ' 1 >>ni60<< ; GREP: MolarMasses :: GETVAL 'Ni61 ' 1 >>ni61<< ; GREP: MolarMasses :: GETVAL 'Ni62 ' 1 >>ni62<< ; GREP: MolarMasses :: GETVAL 'Ni64 ' 1 >>ni64<< ; GREP: MolarMasses :: GETVAL 'O16 ' 1 >>o16<< ; GREP: MolarMasses :: GETVAL 'Pu238' 1 >>pu238<< ; GREP: MolarMasses :: GETVAL 'Pu239' 1 >>pu239<< ; GREP: MolarMasses :: GETVAL 'Pu240' 1 >>pu240<< ; GREP: MolarMasses :: GETVAL 'Pu241' 1 >>pu241<< ; GREP: MolarMasses :: GETVAL 'Pu242' 1 >>pu242<< ; GREP: MolarMasses :: GETVAL 'Sn115' 1 >>sn115<< ; GREP: MolarMasses :: GETVAL 'Sn116' 1 >>sn116<< ; GREP: MolarMasses :: GETVAL 'Sn117' 1 >>sn117<< ; GREP: MolarMasses :: GETVAL 'Sn118' 1 >>sn118<< ; GREP: MolarMasses :: GETVAL 'Sn119' 1 >>sn119<< ; GREP: MolarMasses :: GETVAL 'Sn120' 1 >>sn120<< ; GREP: MolarMasses :: GETVAL 'Sn122' 1 >>sn122<< ; GREP: MolarMasses :: GETVAL 'Sn124' 1 >>sn124<< ; GREP: MolarMasses :: GETVAL 'U234 ' 1 >>u234<< ; GREP: MolarMasses :: GETVAL 'U235 ' 1 >>u235<< ; GREP: MolarMasses :: GETVAL 'U236 ' 1 >>u236<< ; GREP: MolarMasses :: GETVAL 'U238 ' 1 >>u238<< ; GREP: MolarMasses :: GETVAL 'CNat ' 1 >>cnat<< ; GREP: MolarMasses :: GETVAL 'Ag107' 1 >>ag107<< ; GREP: MolarMasses :: GETVAL 'Ag109' 1 >>ag109<< ; GREP: MolarMasses :: GETVAL 'Cd110' 1 >>cd110<< ; GREP: MolarMasses :: GETVAL 'Cd111' 1 >>cd111<< ; GREP: MolarMasses :: GETVAL 'Cd112' 1 >>cd112<< ; GREP: MolarMasses :: GETVAL 'Cd113' 1 >>cd113<< ; GREP: MolarMasses :: GETVAL 'Cd114' 1 >>cd114<< ; GREP: MolarMasses :: GETVAL 'Cd116' 1 >>cd116<< ; GREP: MolarMasses :: GETVAL 'Gd154' 1 >>gd154<< ; GREP: MolarMasses :: GETVAL 'Gd155' 1 >>gd155<< ; GREP: MolarMasses :: GETVAL 'Gd156' 1 >>gd156<< ; GREP: MolarMasses :: GETVAL 'Gd157' 1 >>gd157<< ; GREP: MolarMasses :: GETVAL 'Gd158' 1 >>gd158<< ; GREP: MolarMasses :: GETVAL 'Gd160' 1 >>gd160<< ; GREP: MolarMasses :: GETVAL 'Na23' 1 >>na23<< ; * Isotopes conflictuels * Pour les isotopes conflictuels (absents de certaines bibli.) IF LibType 'ENDFVI8' = THEN GREP: MolarMasses :: GETVAL 'H1 ' 1 >>h1<< ; GREP: MolarMasses :: GETVAL 'SiNat' 1 >>sinat<< ; GREP: MolarMasses :: GETVAL 'TiNat' 1 >>tinat<< ; GREP: MolarMasses :: GETVAL 'ZrNat' 1 >>zrnat<< ; ! ENDF6 connait a la fois ZrNat et ses differents isotopes GREP: MolarMasses :: GETVAL 'InNat' 1 >>in115<< ; GREP: MolarMasses :: GETVAL 'Knat' 1 >>knat<< ; GREP: MolarMasses :: GETVAL 'S32 ' 1 >>s32<< ; ELSEIF LibType 'JEF22' = THEN GREP: MolarMasses :: GETVAL 'H1 ' 1 >>h1<< ; GREP: MolarMasses :: GETVAL 'SiNat' 1 >>sinat<< ; GREP: MolarMasses :: GETVAL 'TiNat' 1 >>tinat<< ; GREP: MolarMasses :: GETVAL 'ZrNat' 1 >>zrnat<< ; GREP: MolarMasses :: GETVAL 'In115' 1 >>in115<< ; GREP: MolarMasses :: GETVAL 'Knat' 1 >>knat<< ; GREP: MolarMasses :: GETVAL 'S32 ' 1 >>s32<< ; GREP: MolarMasses :: GETVAL 'S33 ' 1 >>s33<< ; GREP: MolarMasses :: GETVAL 'S34 ' 1 >>s34<< ; ELSEIF LibType 'JEFF31' = LibType 'JEFF311' = + LibType 'JEFF312' = + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN GREP: MolarMasses :: GETVAL 'H1 ' 1 >>h1<< ; GREP: MolarMasses :: GETVAL 'Si28 ' 1 >>si28<< ; GREP: MolarMasses :: GETVAL 'Si29 ' 1 >>si28<< ; GREP: MolarMasses :: GETVAL 'Si30 ' 1 >>si28<< ; GREP: MolarMasses :: GETVAL 'Ti46 ' 1 >>ti46<< ; GREP: MolarMasses :: GETVAL 'Ti47 ' 1 >>ti47<< ; GREP: MolarMasses :: GETVAL 'Ti48 ' 1 >>ti48<< ; GREP: MolarMasses :: GETVAL 'Ti49 ' 1 >>ti49<< ; GREP: MolarMasses :: GETVAL 'Ti50 ' 1 >>ti50<< ; GREP: MolarMasses :: GETVAL 'Zr90 ' 1 >>zr90<< ; GREP: MolarMasses :: GETVAL 'Zr91 ' 1 >>zr91<< ; GREP: MolarMasses :: GETVAL 'Zr92 ' 1 >>zr92<< ; GREP: MolarMasses :: GETVAL 'Zr93 ' 1 >>zr93<< ; GREP: MolarMasses :: GETVAL 'Zr94 ' 1 >>zr94<< ; GREP: MolarMasses :: GETVAL 'Zr95 ' 1 >>zr95<< ; GREP: MolarMasses :: GETVAL 'Zr96 ' 1 >>zr96<< ; GREP: MolarMasses :: GETVAL 'Cd106' 1 >>cd106<< ; GREP: MolarMasses :: GETVAL 'Cd108' 1 >>cd108<< ; GREP: MolarMasses :: GETVAL 'In115' 1 >>in115<< ; GREP: MolarMasses :: GETVAL 'Sn112' 1 >>sn112<< ; GREP: MolarMasses :: GETVAL 'Sn114' 1 >>sn114<< ; GREP: MolarMasses :: GETVAL 'K39 ' 1 >>k39<< ; GREP: MolarMasses :: GETVAL 'K40 ' 1 >>k40<< ; GREP: MolarMasses :: GETVAL 'K41 ' 1 >>k41<< ; GREP: MolarMasses :: GETVAL 'S32 ' 1 >>s32<< ; GREP: MolarMasses :: GETVAL 'S33 ' 1 >>s33<< ; GREP: MolarMasses :: GETVAL 'S34 ' 1 >>s34<< ; ELSE ! CEA99v9p5, CEA514T2, CLA99CEA93 GREP: MolarMasses :: GETVAL 'H1 ' 1 >>h1<< ; GREP: MolarMasses :: GETVAL 'SiNat' 1 >>sinat<< ; GREP: MolarMasses :: GETVAL 'TiNat' 1 >>tinat<< ; GREP: MolarMasses :: GETVAL 'ZrNat' 1 >>zrnat<< ; GREP: MolarMasses :: GETVAL 'In115' 1 >>in115<< ; GREP: MolarMasses :: GETVAL 'Knat' 1 >>knat<< ; GREP: MolarMasses :: GETVAL 'S32 ' 1 >>s32<< ; GREP: MolarMasses :: GETVAL 'S33 ' 1 >>s33<< ; GREP: MolarMasses :: GETVAL 'S34 ' 1 >>s34<< ; ENDIF ; *------------------------------------------------------ * Recuperation des noms des isotopes *------------------------------------------------------ STRING NomH2O ; STRING NomH1 ; STRING NomO16 ; STRING NomAm241 ; STRING NomAl27 ; STRING NomB10 NomB11 ; STRING NomCo59 NomCr50 NomCr52 NomCr53 NomCr54 ; STRING NomFe54 NomFe56 NomFe57 NomFe58 ; STRING NomMn55 ; STRING NomMo95 ; STRING NomNb93 ; STRING NomNi58 NomNi60 NomNi61 NomNi62 NomNi64 ; STRING NomPu238 NomPu239 NomPu240 NomPu241 NomPu242 ; STRING NomSiNat NomSi28 NomSi29 NomSi30 ; STRING NomSn112 NomSn114 NomSn115 NomSn116 NomSn117 ; STRING NomSn118 NomSn119 NomSn120 NomSn122 NomSn124 ; STRING NomTiNat NomTi46 NomTi47 NomTi48 NomTi49 NomTi50 ; STRING NomU234 NomU235 NomU236 NomU238 ; STRING NomZrNat ; STRING NomZr90 NomZr91 NomZr92 NomZr93 NomZr94 NomZr95 NomZr96 ; STRING NomCNat ; STRING NomAg107 NomAg109 ; STRING NomCd106 NomCd108 NomCd110 NomCd111 NomCd112 NomCd113 NomCd114 ; STRING NomCd116 ; STRING NomIn115 NomInNat ; STRING NomGd154 NomGd155 NomGd156 NomGd157 NomGd158 NomGd160 ; STRING NomKnat NomK39 NomK40 NomK41 NomNa23 ; STRING NomS32 NomS33 NomS34 NomN14 ; STRING NomIn113 ; ListeIsot ListeAutop := AdaptIso :: <> <> >>LibId<< >>LibGenre<< ; UTL: ListeIsot :: DUMP ; GREP: ListeIsot :: GETVAL 'H2O ' <> NVAL 8 >>NomH2O<< ; GREP: ListeIsot :: GETVAL 'H1 ' <> NVAL 8 >>NomH1<< ; GREP: ListeIsot :: GETVAL 'O16 ' <> NVAL 8 >>NomO16<< ; GREP: ListeIsot :: GETVAL 'Am241' <> NVAL 8 >>NomAm241<< ; GREP: ListeIsot :: GETVAL 'Al27 ' <> NVAL 8 >>NomAl27<< ; GREP: ListeIsot :: GETVAL 'B10 ' <> NVAL 8 >>NomB10<< ; GREP: ListeIsot :: GETVAL 'B11 ' <> NVAL 8 >>NomB11<< ; GREP: ListeIsot :: GETVAL 'Co59 ' <> NVAL 8 >>NomCo59<< ; GREP: ListeIsot :: GETVAL 'Cr50 ' <> NVAL 8 >>NomCr50<< ; GREP: ListeIsot :: GETVAL 'Cr52 ' <> NVAL 8 >>NomCr52<< ; GREP: ListeIsot :: GETVAL 'Cr53 ' <> NVAL 8 >>NomCr53<< ; GREP: ListeIsot :: GETVAL 'Cr54 ' <> NVAL 8 >>NomCr54<< ; GREP: ListeIsot :: GETVAL 'Fe54 ' <> NVAL 8 >>NomFe54<< ; GREP: ListeIsot :: GETVAL 'Fe56 ' <> NVAL 8 >>NomFe56<< ; GREP: ListeIsot :: GETVAL 'Fe57 ' <> NVAL 8 >>NomFe57<< ; GREP: ListeIsot :: GETVAL 'Fe58 ' <> NVAL 8 >>NomFe58<< ; GREP: ListeIsot :: GETVAL 'Mn55 ' <> NVAL 8 >>NomMn55<< ; GREP: ListeIsot :: GETVAL 'Mo95 ' <> NVAL 8 >>NomMo95<< ; GREP: ListeIsot :: GETVAL 'Nb93 ' <> NVAL 8 >>NomNb93<< ; GREP: ListeIsot :: GETVAL 'Ni58 ' <> NVAL 8 >>NomNi58<< ; GREP: ListeIsot :: GETVAL 'Ni60 ' <> NVAL 8 >>NomNi60<< ; GREP: ListeIsot :: GETVAL 'Ni61 ' <> NVAL 8 >>NomNi61<< ; GREP: ListeIsot :: GETVAL 'Ni62 ' <> NVAL 8 >>NomNi62<< ; GREP: ListeIsot :: GETVAL 'Ni64 ' <> NVAL 8 >>NomNi64<< ; GREP: ListeIsot :: GETVAL 'Pu238' <> NVAL 8 >>NomPu238<< ; GREP: ListeIsot :: GETVAL 'Pu239' <> NVAL 8 >>NomPu239<< ; GREP: ListeIsot :: GETVAL 'Pu240' <> NVAL 8 >>NomPu240<< ; GREP: ListeIsot :: GETVAL 'Pu241' <> NVAL 8 >>NomPu241<< ; GREP: ListeIsot :: GETVAL 'Pu242' <> NVAL 8 >>NomPu242<< ; GREP: ListeIsot :: GETVAL 'SiNat' <> NVAL 8 >>NomSiNat<< ; GREP: ListeIsot :: GETVAL 'Si28 ' <> NVAL 8 >>NomSi28<< ; GREP: ListeIsot :: GETVAL 'Si29 ' <> NVAL 8 >>NomSi29<< ; GREP: ListeIsot :: GETVAL 'Si30 ' <> NVAL 8 >>NomSi30<< ; GREP: ListeIsot :: GETVAL 'Sn112' <> NVAL 8 >>NomSn112<< ; GREP: ListeIsot :: GETVAL 'Sn114' <> NVAL 8 >>NomSn114<< ; GREP: ListeIsot :: GETVAL 'Sn115' <> NVAL 8 >>NomSn115<< ; GREP: ListeIsot :: GETVAL 'Sn116' <> NVAL 8 >>NomSn116<< ; GREP: ListeIsot :: GETVAL 'Sn117' <> NVAL 8 >>NomSn117<< ; GREP: ListeIsot :: GETVAL 'Sn118' <> NVAL 8 >>NomSn118<< ; GREP: ListeIsot :: GETVAL 'Sn119' <> NVAL 8 >>NomSn119<< ; GREP: ListeIsot :: GETVAL 'Sn120' <> NVAL 8 >>NomSn120<< ; GREP: ListeIsot :: GETVAL 'Sn122' <> NVAL 8 >>NomSn122<< ; GREP: ListeIsot :: GETVAL 'Sn124' <> NVAL 8 >>NomSn124<< ; GREP: ListeIsot :: GETVAL 'TiNat' <> NVAL 8 >>NomTiNat<< ; GREP: ListeIsot :: GETVAL 'Ti46 ' <> NVAL 8 >>NomTi46<< ; GREP: ListeIsot :: GETVAL 'Ti47 ' <> NVAL 8 >>NomTi47<< ; GREP: ListeIsot :: GETVAL 'Ti48 ' <> NVAL 8 >>NomTi48<< ; GREP: ListeIsot :: GETVAL 'Ti49 ' <> NVAL 8 >>NomTi49<< ; GREP: ListeIsot :: GETVAL 'Ti50 ' <> NVAL 8 >>NomTi50<< ; GREP: ListeIsot :: GETVAL 'U234 ' <> NVAL 8 >>NomU234<< ; GREP: ListeIsot :: GETVAL 'U235 ' <> NVAL 8 >>NomU235<< ; GREP: ListeIsot :: GETVAL 'U236 ' <> NVAL 8 >>NomU236<< ; GREP: ListeIsot :: GETVAL 'U238 ' <> NVAL 8 >>NomU238<< ; GREP: ListeIsot :: GETVAL 'ZrNat' <> NVAL 8 >>NomZrNat<< ; GREP: ListeIsot :: GETVAL 'Zr90 ' <> NVAL 8 >>NomZr90<< ; GREP: ListeIsot :: GETVAL 'Zr91 ' <> NVAL 8 >>NomZr91<< ; GREP: ListeIsot :: GETVAL 'Zr92 ' <> NVAL 8 >>NomZr92<< ; GREP: ListeIsot :: GETVAL 'Zr93 ' <> NVAL 8 >>NomZr93<< ; GREP: ListeIsot :: GETVAL 'Zr94 ' <> NVAL 8 >>NomZr94<< ; GREP: ListeIsot :: GETVAL 'Zr95 ' <> NVAL 8 >>NomZr95<< ; GREP: ListeIsot :: GETVAL 'Zr96 ' <> NVAL 8 >>NomZr96<< ; GREP: ListeIsot :: GETVAL 'CNat ' <> NVAL 8 >>NomCNat<< ; GREP: ListeIsot :: GETVAL 'Ag107' <> NVAL 8 >>NomAg107<< ; GREP: ListeIsot :: GETVAL 'Ag109' <> NVAL 8 >>NomAg109<< ; GREP: ListeIsot :: GETVAL 'Cd106' <> NVAL 8 >>NomCd106<< ; GREP: ListeIsot :: GETVAL 'Cd108' <> NVAL 8 >>NomCd108<< ; GREP: ListeIsot :: GETVAL 'Cd110' <> NVAL 8 >>NomCd110<< ; GREP: ListeIsot :: GETVAL 'Cd111' <> NVAL 8 >>NomCd111<< ; GREP: ListeIsot :: GETVAL 'Cd112' <> NVAL 8 >>NomCd112<< ; GREP: ListeIsot :: GETVAL 'Cd113' <> NVAL 8 >>NomCd113<< ; GREP: ListeIsot :: GETVAL 'Cd114' <> NVAL 8 >>NomCd114<< ; GREP: ListeIsot :: GETVAL 'Cd116' <> NVAL 8 >>NomCd116<< ; GREP: ListeIsot :: GETVAL 'In115' <> NVAL 8 >>NomIn115<< ; GREP: ListeIsot :: GETVAL 'In113' <> NVAL 8 >>NomIn113<< ; GREP: ListeIsot :: GETVAL 'Gd154' <> NVAL 8 >>NomGd154<< ; GREP: ListeIsot :: GETVAL 'Gd155' <> NVAL 8 >>NomGd155<< ; GREP: ListeIsot :: GETVAL 'Gd156' <> NVAL 8 >>NomGd156<< ; GREP: ListeIsot :: GETVAL 'Gd157' <> NVAL 8 >>NomGd157<< ; GREP: ListeIsot :: GETVAL 'Gd158' <> NVAL 8 >>NomGd158<< ; GREP: ListeIsot :: GETVAL 'Gd160' <> NVAL 8 >>NomGd160<< ; GREP: ListeIsot :: GETVAL 'Knat ' <> NVAL 8 >>NomKnat<< ; GREP: ListeIsot :: GETVAL 'K39 ' <> NVAL 8 >>NomK39<< ; GREP: ListeIsot :: GETVAL 'K40 ' <> NVAL 8 >>NomK40<< ; GREP: ListeIsot :: GETVAL 'K41 ' <> NVAL 8 >>NomK41<< ; GREP: ListeIsot :: GETVAL 'Na23 ' <> NVAL 8 >>NomNa23<< ; GREP: ListeIsot :: GETVAL 'S32 ' <> NVAL 8 >>NomS32<< ; GREP: ListeIsot :: GETVAL 'S33 ' <> NVAL 8 >>NomS33<< ; GREP: ListeIsot :: GETVAL 'S34 ' <> NVAL 8 >>NomS34<< ; GREP: ListeIsot :: GETVAL 'N14 ' <> NVAL 8 >>NomN14<< ; GREP: ListeIsot :: GETVAL 'InNat' <> NVAL 8 >>NomInNat<< ; * Les noms des isotopes autoproteges STRING NomSSZrNat NomSSZr90 NomSSZr91 NomSSZr92 NomSSZr93 NomSSZr94 ; STRING NomSSZr95 NomSSZr96 ; STRING NomSSAg7 NomSSAg9 ; STRING NomSSIn115 ; STRING NomSSGd4 NomSSGd5 NomSSGd6 NomSSGd7 NomSSGd8 ; STRING NomSSGd0 ; STRING NomSSU235 NomSSU236 NomSSU238 ; STRING NomSSPu8 NomSSPu9 NomSSPu0 NomSSPu1 NomSSPu2 ; STRING NomSSAm241 ; STRING NomSSInNat ; STRING NomSSIn113 ; GREP: ListeIsot :: GETVAL 'Ag107' <> NVAL 8 >>NomSSAg7<< ; GREP: ListeIsot :: GETVAL 'Ag109' <> NVAL 8 >>NomSSAg9<< ; GREP: ListeIsot :: GETVAL 'In115' <> NVAL 8 >>NomSSIn115<< ; GREP: ListeIsot :: GETVAL 'Gd154' <> NVAL 8 >>NomSSGd4<< ; GREP: ListeIsot :: GETVAL 'Gd155' <> NVAL 8 >>NomSSGd5<< ; GREP: ListeIsot :: GETVAL 'Gd156' <> NVAL 8 >>NomSSGd6<< ; GREP: ListeIsot :: GETVAL 'Gd157' <> NVAL 8 >>NomSSGd7<< ; GREP: ListeIsot :: GETVAL 'Gd158' <> NVAL 8 >>NomSSGd8<< ; GREP: ListeIsot :: GETVAL 'Gd160' <> NVAL 8 >>NomSSGd0<< ; GREP: ListeIsot :: GETVAL 'U235 ' <> NVAL 8 >>NomSSU235<< ; GREP: ListeIsot :: GETVAL 'U236 ' <> NVAL 8 >>NomSSU236<< ; GREP: ListeIsot :: GETVAL 'U238 ' <> NVAL 8 >>NomSSU238<< ; GREP: ListeIsot :: GETVAL 'Pu238' <> NVAL 8 >>NomSSPu8<< ; GREP: ListeIsot :: GETVAL 'Pu239' <> NVAL 8 >>NomSSPu9<< ; GREP: ListeIsot :: GETVAL 'Pu240' <> NVAL 8 >>NomSSPu0<< ; GREP: ListeIsot :: GETVAL 'Pu241' <> NVAL 8 >>NomSSPu1<< ; GREP: ListeIsot :: GETVAL 'Pu242' <> NVAL 8 >>NomSSPu2<< ; GREP: ListeIsot :: GETVAL 'Am241' <> NVAL 8 >>NomSSAm241<< ; IF LibType 'ENDFVII1' = LibType 'ENDFVII0' = LibType 'JEFF31' = LibType 'JEFF311' = LibType 'JEFF312' = + + + + THEN GREP: ListeIsot :: GETVAL 'Zr90' <> NVAL 8 >>NomSSZr90<< ; GREP: ListeIsot :: GETVAL 'Zr91' <> NVAL 8 >>NomSSZr91<< ; GREP: ListeIsot :: GETVAL 'Zr92' <> NVAL 8 >>NomSSZr92<< ; GREP: ListeIsot :: GETVAL 'Zr94' <> NVAL 8 >>NomSSZr94<< ; GREP: ListeIsot :: GETVAL 'Zr96' <> NVAL 8 >>NomSSZr96<< ; ELSEIF LibType 'ENDFVI8' = LibType 'JEF22' = + THEN GREP: ListeIsot :: GETVAL 'ZrNat' <> NVAL 8 >>NomSSZrNat<< ; IF LibType 'ENDFVI8' = THEN GREP: ListeIsot :: GETVAL 'InNat' <> NVAL 8 >>NomSSInNat<< ; ENDIF ; ENDIF ; *------------------------------------------------------ REAL rFuel dUO2 TempF Alpha Beta ; IF Gestion "TIH" = THEN * Specifications TIHANGE CEA-N-2092 ! Thermal dilatation (p.153) EVALUATE TempF := 286.0 (* C *) 9.0 * 5.0 / 32.0 + ; ! C to F EVALUATE Beta := 1.0 -3.31925E-4 + 4.87957E-6 TempF * + 1.37101E-11 TempF 2.0 ** * + 1.6161E-13 TempF 3.0 ** * + ; EVALUATE rFuel := .464693 Beta * ; ! (p. 158) *------------------------------------------------------------------- * Theoric density = 10.96 g/cm3 * Porosity = 5.0 % (x0.95) -> 10.412 (p.154) * Thermal Dilatation Effects -> 10.338 (p.153/154) * Free Rod Space Modelisation, Grain boundary, Others Effects * --> Dilution factor = 0.9881 (p.159) * Final density at 286 deg[C] = 10.215 g/cm3 *------------------------------------------------------------------- EVALUATE dUO2 := 10.96 0.95 * Beta 3.0 ** / 0.9881 * ; ELSEIF Gestion "BEAV" = THEN IF Teneur_I 160 = THEN EVALUATE dUO2 := 10.31341 ; ELSEIF Teneur_I 240 = THEN EVALUATE dUO2 := 10.29748 ; ELSEIF Teneur_I 310 = THEN EVALUATE dUO2 := 10.30166 ; ELSE ECHO "BEAVRS dispose d'enrichissements specifiques. Les seuls" ; ECHO "enrichissements disponibles sont 1,6, 2,4 et 3,10" ; ECHO "pourcents. Specifier l'un de ces trois enrichissements." ; ECHO "Actuellement :" ; ECHO Main_enri ; ABORT: ; ENDIF ; EVALUATE rFuel := .39218 ; ! [cm] ELSE EVALUATE rFuel := .4095 ; ! [cm] EVALUATE dUO2 := 10.07 ; ! UO2-density ENDIF ; REAL dMOX := 10.02 ; REAL dGado := 9.83 ; REAL dInc dZr4 ; REAL U5FC := 1E-9 ; *-------------------------------------------------------------------- * Assembly structure *-------------------------------------------------------------------- REAL dTube dClad dAICr dSS dGrid dB4C dBP ; * Specifications TIHANGE CEA-N-2092 EVALUATE dInc := 8.089 ; ! In-718 density EVALUATE dZr4 := 6.52 ; IF Gestion "BEAV" = THEN EVALUATE dInc := 8.2 ; ! In-718 density EVALUATE dZr4 := 6.55 ; ENDIF ; IF Gestion "TIH" = THEN * Specifications TIHANGE CEA-N-2092 ! AgInCd Control Rod density EVALUATE dAICr := 10.17 ; ! 20 C EVALUATE TempF := 286.0 (* C *) 9.0 * 5.0 / 32.0 + ; ! C to F EVALUATE Alpha := 1.25E-5 ; EVALUATE dAICr := dAICr TempF 68.0 - Alpha * 1.0 + 3.0 ** / ; EVALUATE dSS := 7.79 ; ! Stainless Steel (SS304) density EVALUATE dTube := dZr4 ; ! Guide Tube density EVALUATE dClad := dZr4 ; ! Fuel Cladding density EVALUATE dBP := 2.27 ; ! BoroSilicate Poison Rod (BP) ELSEIF Gestion "BEAV" = THEN EVALUATE dAICr := 10.16 ; ! 20 C EVALUATE dSS := 8.03 ; ! Stainless Steel (SS304) density EVALUATE dTube := dZr4 ; ! Guide Tube density EVALUATE dClad := dZr4 ; ! Fuel Cladding density EVALUATE dBP := 2.26 ; ! BoroSilicate Poison Rod (BP) ELSE ECHO "Error --- Gestion:" Gestion "not recognized." ; ABORT: ; ENDIF ; EVALUATE dB4C := 1.79 ; EVALUATE dGrid := 6.691985 ; *------------------------------------------------------ * Gestion Branching *------------------------------------------------------ IF Gestion "TIH" = THEN EVALUATE PowMWth := 2652. ; ! [MW] (CEA-N-2092) EVALUATE cbor := 600. ; ! [ppm] EVALUATE Hactiv := 368. ; ! [cm] EVALUATE NbAss := 157. ; ! [900MWe] EVALUATE dmod := .716403 ; ! [g.cm-3] EVALUATE tmod := 305. ; ! [degC] EVALUATE tcom := 650. ; ! [degC] EVALUATE tgain := 335. ; ! [degC] ELSEIF Gestion "BEAV" = THEN EVALUATE PowMWth := 3411. ; ! [MW] (SLB 1100MWe) EVALUATE Hactiv := 365.76 ; ! [cm] EVALUATE NbAss := 193. ; ! [1300MWe] EVALUATE cbor := 975. ; ! [ppm] EVALUATE dmod := .7398625 ; ! [g.cm-3] EVALUATE tmod := 293.33 ; ! [degC] EVALUATE tcom := 360. ; ! [degC] EVALUATE tgain := 335. ; ! [degC] ELSE ! Input Error ECHO "Error --- Gestion:" Gestion "not recognized." ; ABORT: ; ENDIF ; *------------------------------------------------------ * Thermal Expansion of densities (If requested) *------------------------------------------------------ IF ThermExp 1 = THEN EVALUATE tclad := 0.12 tcom * 0.88 tmod * + ; EVALUATE Model := "TIH" ; * EVALUATE tclad := tgain ; ThermaExpans :: <> "UO2" <> "Density" <> "NONE" 0. 0. 0. >>dUO2<< ; ThermaExpans :: <> "AIC" <> "Density" <> "NONE" 0. 0. 0. >>dAICr<< ; ThermaExpans :: <> "SS304" <> "Density" <> "NONE" 0. 0. 0. >>dTube<< ; ThermaExpans :: <> "Pyrex" <> "Density" <> "NONE" 0. 0. 0. >>dBP<< ; ThermaExpans :: <> "Zr4" <> "Density" <> "NONE" 0. 0. 0. >>dClad<< ; ThermaExpans :: <> "In718" <> "Density" <> "NONE" 0. 0. 0. >>dInc<< ; ThermaExpans :: <> "SS304" <> "Density" <> "NONE" 0. 0. 0. >>dSS<< ; ECHO "---------------THERMAL EXPANSION OF DENSITIES------------" ; ECHO " dUO2 = " dUO2 " at " tcom "C." ; ECHO " dAICr = " dAICr " at " tmod "C." ; ECHO " dSS = " dSS " at " tmod "C." ; ECHO " dClad = " dClad " at " tclad "C." ; ECHO " dTube = " dTube " at " tcom "C." ; ECHO " dBP = " dBP " at " tcom "C." ; ECHO " dInc = " dInc " at " tcom "C." ; ENDIF ; *------ GAINC ------ *---- made of Zircalloy4 (with fuel-clad gap) ---- EVALUATE dClad := dClad 0.53677 2.0 ** 0.474956 2.0 ** - 0.53677 2.0 ** 0.465793 2.0 ** - / * ; *------ GRILLE ------------------- *--------------------------------- * Mass Ratios Grille precalculee : * Inconel .131668 * Zirc4 .868332 *--------------------------------- REAL InGridf := .131668 ; REAL Z4Gridf := .868332 ; REAL SSGridf := 0.0 ; * SS rods EVALUATE dRodSS := dSS .5525475 * ; ! * AIC IF Gestion "TIH" = THEN EVALUATE dCladAIC := dSS 0.5601362 2.0 ** 0.5110126 2.0 ** - 0.5601362 2.0 ** 0.5078476 2.0 ** - / * ; ELSEIF Gestion "BEAV" = THEN EVALUATE dCladAIC := dSS 0.48387 2.0 ** 0.43688 2.0 ** - 0.48387 2.0 ** 0.43310 2.0 ** - / * ; ELSE ! Input Error ECHO "Error --- Gestion:" Gestion "not recognized." ; ABORT: ; ENDIF ; * EVALUATE dCladSS := dSS .87575664 * ; EVALUATE dCladB4C := dSS .8868781 * ; * PYREX REAL dCladIntBP := dSS 0.301759 2.0 ** 0.2851718 2.0 ** - 0.31115 2.0 ** 0.2851718 2.0 ** - / * ; ! 0.56013513697 = 0,55753*(1+9,8E-6*(546,8-70)) REAL dCladExtBP := dSS 0.5601362 2.0 ** 0.5110126 2.0 ** - 0.5601362 2.0 ** 0.49784 2.0 ** - / * ; IF Gestion "BEAV" = THEN EVALUATE dCladIntBP := dSS 0.2305 2.0 ** 0.2140 2.0 ** - 0.2413 2.0 ** 0.2140 2.0 ** - / * ; EVALUATE dCladExtBP := dSS 0.4839 2.0 ** 0.4369 2.0 ** - 0.4839 2.0 ** 0.4267 2.0 ** - / * ; ENDIF ; *------------------------------------------------------ * Convert values (degC -> K ; ppm -> value) *------------------------------------------------------ EVALUATE tmod := tmod 273.15 + ; ! user temp in degC, LIB in K EVALUATE tcom := tcom 273.15 + ; ! user temp in degC, LIB in K EVALUATE tgain := tgain 273.15 + ; ! user temp in degC, LIB in K EVALUATE cbor := cbor 1.E-06 * ; ! *------------------------------------------------------ * Grid Dilution in Moderator for every cell *------------------------------------------------------ * Le mot-clef COMB interdit 0.0, mais autorise une quantite * negligeable. Faut-il faire une demande de dev pour qu'il * accepte une quantite negligeable ? REAL SurfEauCntr SurfEauSide SurfEauCorn ; REAL ss304Cntr IncCntr Zr4Cntr ss304Side IncSide Zr4Side ss304Corn IncCorn Zr4Corn ss304Tube IncTube Zr4Tube ; REAL DemiLameH2O SurfEauTube := 0.0 0.0 ; REAL VtotZr4 VtotInc VtotSS Vtot Zf Incf SSf := 0.0 .0 .0 .0 .0 .0 .0 ; REAL MtotZr4 MtotInc ; IF Gestion 'TIH' = THEN ! p155 et 159 EVALUATE SurfEauCntr := 1.43538098 2.0 ** $Pi_R 0.53677 2.0 ** * - ; EVALUATE SurfEauSide := 1.43538098 0.07343 2.0 / * ; EVALUATE SurfEauCorn := SurfEauSide 2.0 * 0.07343 2.0 / 0.07343 2.0 / * + ; * EVALUATE ss304Tube := 2.205 100.0 / ; EVALUATE IncTube := 1.162 100.0 / ; EVALUATE Zr4Tube := 1.0E-20 ; * EVALUATE ss304Cntr := 1.0E-20 ; EVALUATE IncCntr := 0.5484 100.0 / ; EVALUATE Zr4Cntr := 1.0E-20 ; * EVALUATE ss304Side := 1.0E-20 ; EVALUATE IncSide := 0.5484 100.0 / SurfEauCntr SurfEauCntr SurfEauSide + / * 3.041 100.0 / SurfEauSide SurfEauCntr SurfEauSide + / * + ; EVALUATE Zr4Side := 1.0E-20 ; * EVALUATE ss304Corn := 1.0E-20 ; EVALUATE IncCorn := 0.5484 100.0 / SurfEauCntr SurfEauCntr SurfEauCorn + / * 3.041 100.0 / SurfEauCorn SurfEauCntr SurfEauCorn + / * + ; EVALUATE Zr4Corn := 1.0E-20 ; * EVALUATE GrPreTube := 1.0E-20 ; EVALUATE GrPreCntr := 1.0E-20 ; EVALUATE GrPreSide := 1.0E-20 ; EVALUATE GrPreCorn := 1.0E-20 ; * EVALUATE H2OTube := 1.0 ss304Tube - IncTube - Zr4Tube - ; EVALUATE H2OCntr := 1.0 ss304Cntr - IncCntr - Zr4Cntr - ; EVALUATE H2OSide := 1.0 ss304Side - IncSide - Zr4Side - ; EVALUATE H2OCorn := 1.0 ss304Corn - IncCorn - Zr4Corn - ; ELSEIF Gestion "BEAV" = THEN ! Le calcul des rapport volumique est caclule dans le fichier Dilution_grille.xslx * EVALUATE ss304Tube := 1.0E-20 ; EVALUATE IncTube := 1.00254E-03 ; EVALUATE Zr4Tube := 2.08644E-02 ; * EVALUATE ss304Cntr := 1.0E-20 ; EVALUATE IncCntr := 4.835E-04 ; EVALUATE Zr4Cntr := 1.006E-02 ; * EVALUATE ss304Side := 4.06144E-04 ; EVALUATE IncSide := 4.56781E-04 ; EVALUATE Zr4Side := 1.28293E-02 ; * EVALUATE ss304Corn := 7.76842E-04 ; EVALUATE IncCorn := 4.32097E-04 ; EVALUATE Zr4Corn := 1.53486E-02 ; * EVALUATE GrPreTube := 1.0E-20 ; EVALUATE GrPreCntr := 1.0E-20 ; EVALUATE GrPreSide := 1.0E-20 ; EVALUATE GrPreCorn := 1.0E-20 ; * EVALUATE H2OTube := 1.0 ss304Tube - IncTube - Zr4Tube - ; EVALUATE H2OCntr := 1.0 ss304Cntr - IncCntr - Zr4Cntr - ; EVALUATE H2OSide := 1.0 ss304Side - IncSide - Zr4Side - ; EVALUATE H2OCorn := 1.0 ss304Corn - IncCorn - Zr4Corn - ; ELSE ! Input Error ECHO "Error --- Gestion:" Gestion "not recognized." ; ABORT: ; ENDIF ; ECHO "CreaMix: GrPreTube=" GrPreTube ; *------------------------------------------------------ * Evaluate Total Fuel Volume *------------------------------------------------------ REAL VolPin := $Pi_R rFuel * rFuel * Hactiv * ; REAL VolAss ; IF Gestion "TIH" = THEN EVALUATE VolAss := VolPin 204. * ; ! 204 fuel pins ELSE EVALUATE VolAss := VolPin 264. * ; ! 264 fuel pins ENDIF ; REAL TotFuelVol := VolAss NbAss * ; *------------------------------------------------------ * AssType Branching *------------------------------------------------------ IF AssType "UOX" = THEN EVALUATE iU8 iPu9 iPu0 iPu1 iPu2 iU5 iAg9 iAg7 iIn iZrG iZrT := 1 2 3 4 5 6 7 8 9 10 11 ; EVALUATE iGd4 iGd5 iGd6 iGd7 iGd8 iGd0 := 0 0 0 0 0 0 ; EVALUATE enriU := Main_enri ; IF Gestion "TIH" = THEN IF Teneur_I 195 = THEN EVALUATE U4enri := 0.019 100.0 / ; EVALUATE U6enri := 0.012 100.0 / ; ELSEIF Teneur_I 255 = THEN EVALUATE U4enri := 0.025 100.0 / ; EVALUATE U6enri := 0.016 100.0 / ; ELSEIF Teneur_I 310 = THEN EVALUATE U4enri := 0.031 100.0 / ; EVALUATE U6enri := 0.019 100.0 / ; ELSE ECHO "Tihange dispose d'enrichissements specifiques (uranium " ; ECHO "de retraitement re-enrichi). Les seuls enrichissements " ; ECHO "disponibles sont 1,95, 2,55 et 3,10 pourcents. Specifier" ; ECHO " l'un de ces trois enrichissements. Actuellement :" ; ECHO Main_enri ; ABORT: ; ENDIF ; ELSEIF Gestion "BEAV" = THEN IF Teneur_I 160 = THEN EVALUATE U4enri := 0.01288048 100.0 / ; EVALUATE U6enri := 0.0 ; EVALUATE enriU := 1.61 ; ELSEIF Teneur_I 240 = THEN EVALUATE U4enri := 0.01919944 100.0 / ; EVALUATE U6enri := 0.0 ; ELSEIF Teneur_I 310 = THEN EVALUATE U4enri := 0.02481768 100.0 / ; EVALUATE U6enri := 0.0 ; ELSE ECHO "BEAVRS dispose d'enrichissements specifiques. " ; ECHO "Les seuls enrichissements " ; ECHO "disponibles sont 1,60, 2,40 et 3,10 pourcents. Specifier" ; ECHO " l'un de ces trois enrichissements. Actuellement :" ; ECHO Main_enri ; ABORT: ; ENDIF ; ELSE ! Input Error ECHO "Error --- Gestion:" Gestion "not recognized." ; ABORT: ; ENDIF ; EVALUATE U5enri := enriU 100. / ; ! X.XX % => 0.0XXX EVALUATE U8enri := 1. U6enri - U5enri - U4enri - ; EVALUATE Pu8enri := 1.E-15 ; EVALUATE Pu9enri := 1.E-15 ; EVALUATE Pu0enri := 1.E-15 ; EVALUATE Pu1enri := 1.E-15 ; EVALUATE Pu2enri := 1.E-15 ; EVALUATE Am1enri := 1.E-15 ; *---- massic -> isotopic ratio ---- EVALUATE U4enri := U4enri avo_conc * u234 / ; EVALUATE U5enri := U5enri avo_conc * u235 / ; EVALUATE U6enri := U6enri avo_conc * u236 / ; EVALUATE U8enri := U8enri avo_conc * u238 / ; EVALUATE Pu8enri := Pu8enri avo_conc * pu238 / ; EVALUATE Pu9enri := Pu9enri avo_conc * pu239 / ; EVALUATE Pu0enri := Pu0enri avo_conc * pu240 / ; EVALUATE Pu1enri := Pu1enri avo_conc * pu241 / ; EVALUATE Pu2enri := Pu2enri avo_conc * pu242 / ; EVALUATE Am1enri := Am1enri avo_conc * am241 / ; EVALUATE ReNorm := U4enri U5enri + U6enri + U8enri + Pu8enri + Pu9enri + Pu0enri + Pu1enri + Pu2enri + Am1enri + ; EVALUATE U4enri := U4enri ReNorm / ; EVALUATE U5enri := U5enri ReNorm / ; EVALUATE U6enri := U6enri ReNorm / ; EVALUATE U8enri := U8enri ReNorm / ; EVALUATE Pu8enri := Pu8enri ReNorm / ; EVALUATE Pu9enri := Pu9enri ReNorm / ; EVALUATE Pu0enri := Pu0enri ReNorm / ; EVALUATE Pu1enri := Pu1enri ReNorm / ; EVALUATE Pu2enri := Pu2enri ReNorm / ; EVALUATE Am1enri := Am1enri ReNorm / ; *---- average U, UO2 molar masses ---- EVALUATE MUPuAm := U4enri u234 * U5enri u235 * + U6enri u236 * + U8enri u238 * + Pu8enri pu238 * + Pu9enri pu239 * + Pu0enri pu240 * + Pu1enri pu241 * + Pu2enri pu242 * + Am1enri am241 * + ; EVALUATE MUO2 := MUPuAm 2. o16 * + ; *---- U density; U & 0 concentrations ---- EVALUATE U_tot := dUO2 avo_conc * MUO2 / ; ! conc_U = conc_UO2 EVALUATE O_fuel := 2. U_tot * ; ! conc_O = 2 x conc_UO2 EVALUATE dU := U_tot MUPuAm * avo_conc / ; ! U density in g/cm**3 *---------------------------------------------------------------------- * Final U+Pu+Am Concentrations in atm/cm/barn *---------------------------------------------------------------------- EVALUATE U4fuel := U_tot U4enri * ; EVALUATE U5fuel := U_tot U5enri * ; EVALUATE U6fuel := U_tot U6enri * ; EVALUATE U8fuel := U_tot U8enri * ; EVALUATE Pu8fuel := U_tot Pu8enri * ; EVALUATE Pu9fuel := U_tot Pu9enri * ; EVALUATE Pu0fuel := U_tot Pu0enri * ; EVALUATE Pu1fuel := U_tot Pu1enri * ; EVALUATE Pu2fuel := U_tot Pu2enri * ; EVALUATE Am1fuel := U_tot Am1enri * ; ECHO "U4fuel" U4fuel ; ECHO "U5fuel" U5fuel ; ECHO "U6fuel" U6fuel ; ECHO "U8fuel" U8fuel ; ECHO "Pu8fuel" Pu8fuel ; ECHO "Pu9fuel" Pu9fuel ; ECHO "Pu0fuel" Pu0fuel ; ECHO "Pu1fuel" Pu1fuel ; ECHO "Pu2fuel" Pu2fuel ; ECHO "Am1fuel" Am1fuel ; *---------------------------------------------------------------------- * Power conversion factor for UOX: * PowMWT = PowMWth / ( dU * TotFuelVol * 1.E-6 ) * (MW/T) (MW) (g/cm**3) (cm**3) (T/g) *---------------------------------------------------------------------- EVALUATE PowMWT := PowMWth dU TotFuelVol * 1.E-6 * / ; *---------------------------------------------------------------------- ELSEIF AssType "MOX" = THEN EVALUATE iU8 iPu9 iPu0 iPu1 iPu2 iU5 iAg9 iAg7 iIn iZrG iZrT := 1 2 3 4 5 6 7 8 9 10 11 ; EVALUATE iGd4 iGd5 iGd6 iGd7 iGd8 iGd0 := 0 0 0 0 0 0 ; *---- enrichment in isotopic ratio ---- *---- U enrichment ---- EVALUATE U4enri := U4enri avo_conc * u234 / ; EVALUATE U5enri := U5enri avo_conc * u235 / ; EVALUATE U6enri := U6enri avo_conc * u236 / ; EVALUATE U8enri := U8enri avo_conc * u238 / ; EVALUATE ReNorm := U4enri U5enri + U6enri + U8enri + ; EVALUATE U4enri := U4enri ReNorm / ; EVALUATE U5enri := U5enri ReNorm / ; EVALUATE U6enri := U6enri ReNorm / ; EVALUATE U8enri := U8enri ReNorm / ; *---- average Utot ---- EVALUATE MUtot := U4enri u234 * U5enri u235 * + U6enri u236 * + U8enri u238 * + ; *---- Pu enrichment ---- EVALUATE Pu8enri := Pu8enri avo_conc * pu238 / ; EVALUATE Pu9enri := Pu9enri avo_conc * pu239 / ; EVALUATE Pu0enri := Pu0enri avo_conc * pu240 / ; EVALUATE Pu1enri := Pu1enri avo_conc * pu241 / ; EVALUATE Pu2enri := Pu2enri avo_conc * pu242 / ; EVALUATE Am1enri := Am1enri avo_conc * am241 / ; EVALUATE ReNorm := Pu8enri Pu9enri + Pu0enri + Pu1enri + Pu2enri + Am1enri + ; EVALUATE Pu8enri := Pu8enri ReNorm / ; EVALUATE Pu9enri := Pu9enri ReNorm / ; EVALUATE Pu0enri := Pu0enri ReNorm / ; EVALUATE Pu1enri := Pu1enri ReNorm / ; EVALUATE Pu2enri := Pu2enri ReNorm / ; EVALUATE Am1enri := Am1enri ReNorm / ; *---- average Putot ---- EVALUATE MPutot := Pu8enri pu238 * Pu9enri pu239 * + Pu0enri pu240 * + Pu1enri pu241 * + Pu2enri pu242 * + Am1enri am241 * + ; *-------- CORNER -------- EVALUATE Pue_corn := f_corn ; EVALUATE Ue_corn := 1. Pue_corn - ; EVALUATE Pue_corn := Pue_corn avo_conc * MPutot / ; EVALUATE Ue_corn := Ue_corn avo_conc * MUtot / ; EVALUATE ReNorm := Pue_corn Ue_corn + ; EVALUATE Pue_corn := Pue_corn ReNorm / ; EVALUATE Ue_corn := Ue_corn ReNorm / ; EVALUATE U4e_corn := U4enri Ue_corn * ; EVALUATE U5e_corn := U5enri Ue_corn * ; EVALUATE U6e_corn := U6enri Ue_corn * ; EVALUATE U8e_corn := U8enri Ue_corn * ; EVALUATE Pu8e_corn := Pu8enri Pue_corn * ; EVALUATE Pu9e_corn := Pu9enri Pue_corn * ; EVALUATE Pu0e_corn := Pu0enri Pue_corn * ; EVALUATE Pu1e_corn := Pu1enri Pue_corn * ; EVALUATE Pu2e_corn := Pu2enri Pue_corn * ; EVALUATE Am1e_corn := Am1enri Pue_corn * ; *-------- SIDE -------- EVALUATE Pue_side := f_side ; EVALUATE Ue_side := 1. Pue_side - ; EVALUATE Pue_side := Pue_side avo_conc * MPutot / ; EVALUATE Ue_side := Ue_side avo_conc * MUtot / ; EVALUATE ReNorm := Pue_side Ue_side + ; EVALUATE Pue_side := Pue_side ReNorm / ; EVALUATE Ue_side := Ue_side ReNorm / ; EVALUATE U4e_side := U4enri Ue_side * ; EVALUATE U5e_side := U5enri Ue_side * ; EVALUATE U6e_side := U6enri Ue_side * ; EVALUATE U8e_side := U8enri Ue_side * ; EVALUATE Pu8e_side := Pu8enri Pue_side * ; EVALUATE Pu9e_side := Pu9enri Pue_side * ; EVALUATE Pu0e_side := Pu0enri Pue_side * ; EVALUATE Pu1e_side := Pu1enri Pue_side * ; EVALUATE Pu2e_side := Pu2enri Pue_side * ; EVALUATE Am1e_side := Am1enri Pue_side * ; *-------- CENTER -------- EVALUATE Pue_cntr := f_cntr ; EVALUATE Ue_cntr := 1. Pue_cntr - ; EVALUATE Pue_cntr := Pue_cntr avo_conc * MPutot / ; EVALUATE Ue_cntr := Ue_cntr avo_conc * MUtot / ; EVALUATE ReNorm := Pue_cntr Ue_cntr + ; EVALUATE Pue_cntr := Pue_cntr ReNorm / ; EVALUATE Ue_cntr := Ue_cntr ReNorm / ; EVALUATE U4e_cntr := U4enri Ue_cntr * ; EVALUATE U5e_cntr := U5enri Ue_cntr * ; EVALUATE U6e_cntr := U6enri Ue_cntr * ; EVALUATE U8e_cntr := U8enri Ue_cntr * ; EVALUATE Pu8e_cntr := Pu8enri Pue_cntr * ; EVALUATE Pu9e_cntr := Pu9enri Pue_cntr * ; EVALUATE Pu0e_cntr := Pu0enri Pue_cntr * ; EVALUATE Pu1e_cntr := Pu1enri Pue_cntr * ; EVALUATE Pu2e_cntr := Pu2enri Pue_cntr * ; EVALUATE Am1e_cntr := Am1enri Pue_cntr * ; *---- average UPuAm, MOX molar masses ---- EVALUATE M_corn := Pue_corn MPutot * Ue_corn MUtot * + ; EVALUATE M_side := Pue_side MPutot * Ue_side MUtot * + ; EVALUATE M_cntr := Pue_cntr MPutot * Ue_cntr MUtot * + ; EVALUATE MMOXcorn := M_corn 2. o16 * + ; EVALUATE MMOXside := M_side 2. o16 * + ; EVALUATE MMOXcntr := M_cntr 2. o16 * + ; *---- UPu density; MOX, UPu & 0 concentrations ---- EVALUATE MOX_corn := dMOX avo_conc * MMOXcorn / ; EVALUATE MOX_side := dMOX avo_conc * MMOXside / ; EVALUATE MOX_cntr := dMOX avo_conc * MMOXcntr / ; EVALUATE UPu_corn_tot := MOX_corn ; ! conc_UPu = conc_MOX EVALUATE UPu_side_tot := MOX_side ; ! EVALUATE UPu_cntr_tot := MOX_cntr ; ! EVALUATE O_corn_fuel := 2. MOX_corn * ; ! conc_O = 2 x conc_MOX EVALUATE O_side_fuel := 2. MOX_side * ; ! EVALUATE O_cntr_fuel := 2. MOX_cntr * ; ! EVALUATE dUcorn := UPu_corn_tot M_corn * avo_conc / ; EVALUATE dUside := UPu_side_tot M_side * avo_conc / ; EVALUATE dUcntr := UPu_cntr_tot M_cntr * avo_conc / ; *---------------------------------------------------------------------- * Final U+Pu+Am Concentrations in atm/cm/barn *--- Corner ----------------------------------------------------------- EVALUATE U4cornfuel := UPu_corn_tot U4e_corn * ; EVALUATE U5cornfuel := UPu_corn_tot U5e_corn * ; EVALUATE U6cornfuel := UPu_corn_tot U6e_corn * ; EVALUATE U8cornfuel := UPu_corn_tot U8e_corn * ; EVALUATE Pu8cornfuel := UPu_corn_tot Pu8e_corn * ; EVALUATE Pu9cornfuel := UPu_corn_tot Pu9e_corn * ; EVALUATE Pu0cornfuel := UPu_corn_tot Pu0e_corn * ; EVALUATE Pu1cornfuel := UPu_corn_tot Pu1e_corn * ; EVALUATE Pu2cornfuel := UPu_corn_tot Pu2e_corn * ; EVALUATE Am1cornfuel := UPu_corn_tot Am1e_corn * ; *--- Side ------------------------------------------------------------- EVALUATE U4sidefuel := UPu_side_tot U4e_side * ; EVALUATE U5sidefuel := UPu_side_tot U5e_side * ; EVALUATE U6sidefuel := UPu_side_tot U6e_side * ; EVALUATE U8sidefuel := UPu_side_tot U8e_side * ; EVALUATE Pu8sidefuel := UPu_side_tot Pu8e_side * ; EVALUATE Pu9sidefuel := UPu_side_tot Pu9e_side * ; EVALUATE Pu0sidefuel := UPu_side_tot Pu0e_side * ; EVALUATE Pu1sidefuel := UPu_side_tot Pu1e_side * ; EVALUATE Pu2sidefuel := UPu_side_tot Pu2e_side * ; EVALUATE Am1sidefuel := UPu_side_tot Am1e_side * ; *--- Center ----------------------------------------------------------- EVALUATE U4cntrfuel := UPu_cntr_tot U4e_cntr * ; EVALUATE U5cntrfuel := UPu_cntr_tot U5e_cntr * ; EVALUATE U6cntrfuel := UPu_cntr_tot U6e_cntr * ; EVALUATE U8cntrfuel := UPu_cntr_tot U8e_cntr * ; EVALUATE Pu8cntrfuel := UPu_cntr_tot Pu8e_cntr * ; EVALUATE Pu9cntrfuel := UPu_cntr_tot Pu9e_cntr * ; EVALUATE Pu0cntrfuel := UPu_cntr_tot Pu0e_cntr * ; EVALUATE Pu1cntrfuel := UPu_cntr_tot Pu1e_cntr * ; EVALUATE Pu2cntrfuel := UPu_cntr_tot Pu2e_cntr * ; EVALUATE Am1cntrfuel := UPu_cntr_tot Am1e_cntr * ; *---------------------------------------------------------------------- * Power conversion factor for MOX: * PowMWT = PowMWth / [ ( dUcorn * TotFuelVol * f_cornPin ) + * (MW/T) (MW) ( dUside * TotFuelVol * f_sidePin ) + * ( dUcntr * TotFuelVol * f_cntrPin ) ] * 1.E-6 * (g/cm**3) (cm**3) (unitless) (T/g) *---------------------------------------------------------------------- EVALUATE PowMWT := PowMWth dUcorn TotFuelVol * f_cornPin * dUside TotFuelVol * f_sidePin * + dUcntr TotFuelVol * f_cntrPin * + 1.E-6 * / ; *---------------------------------------------------------------------- ELSEIF AssType "UGd" = AssType "UGd8" = + THEN EVALUATE iU8 iGd5 iGd7 iGd4 iGd6 iGd8 iGd0 iPu9 iPu0 iPu1 iPu2 iU5 iAg9 iAg7 iIn iZrG iZrT := 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 ; EVALUATE enriU := Main_enri ; EVALUATE enriUGd := 3.0 ; EVALUATE enriGd := 8. ; EVALUATE f_GadoPin := 12. 264. / ; ! 12 pins with Gd out of 264 tot EVALUATE f_NoGdPin := 252. 264. / ; ! 252 regular pins out of 264 tot *---------------------------------------------------------------------- * GADOLINIUM FUEL *---------------------------------------------------------------------- EVALUATE U4eGd := 1.E-15 ; EVALUATE U5eGd := enriUGd 100. / ; EVALUATE U6eGd := 1.E-15 ; EVALUATE Pu9eGd := 1.E-15 ; EVALUATE Pu0eGd := 1.E-15 ; EVALUATE Pu1eGd := 1.E-15 ; EVALUATE Pu2eGd := 1.E-15 ; EVALUATE U8eGd := 1. U4eGd - U5eGd - U6eGd - Pu9eGd - Pu0eGd - Pu1eGd - Pu2eGd - ; *---- enrichment in isotopic ratio ---- EVALUATE U4eGd := U4eGd avo_conc * u234 / ; EVALUATE U5eGd := U5eGd avo_conc * u235 / ; EVALUATE U6eGd := U6eGd avo_conc * u236 / ; EVALUATE U8eGd := U8eGd avo_conc * u238 / ; EVALUATE Pu9eGd := Pu9eGd avo_conc * pu239 / ; EVALUATE Pu0eGd := Pu0eGd avo_conc * pu240 / ; EVALUATE Pu1eGd := Pu1eGd avo_conc * pu241 / ; EVALUATE Pu2eGd := Pu2eGd avo_conc * pu242 / ; EVALUATE ReNorm := U4eGd U5eGd + U6eGd + U8eGd + Pu9eGd + Pu0eGd + Pu1eGd + Pu2eGd + ; EVALUATE U4eGd := U4eGd ReNorm / ; EVALUATE U5eGd := U5eGd ReNorm / ; EVALUATE U6eGd := U6eGd ReNorm / ; EVALUATE U8eGd := U8eGd ReNorm / ; EVALUATE Pu9eGd := Pu9eGd ReNorm / ; EVALUATE Pu0eGd := Pu0eGd ReNorm / ; EVALUATE Pu1eGd := Pu1eGd ReNorm / ; EVALUATE Pu2eGd := Pu2eGd ReNorm / ; *---- average U, UO2 molar masses for fuel with Gd ---- EVALUATE MUPuGd := U4eGd u234 * U5eGd u235 * + U6eGd u236 * + U8eGd u238 * + Pu9eGd pu239 * + Pu0eGd pu240 * + Pu1eGd pu241 * + Pu2eGd pu242 * + ; EVALUATE MUO2Gd := MUPuGd 2. o16 * + ; *---------------------------------------------------------------------- * Gd2O3 and UO2 for UO2 + Gd2O3 mix * Massic proportions: Gd2O3 (enriGd / 100) * UO2 for Gd (1 - enriGd / 100) *---------------------------------------------------------------------- EVALUATE Gd2O3f := enriGd 100. / ; EVALUATE UO2Gdf := 1. Gd2O3f - ; GREP: Gd2O3 :: GETVAL 'AvgMolM ' 1 >>MGd2O3<< ; GREP: Gd2O3 :: GETVAL 'GadOxG4frIso' 1 >>Gd2O3_G4f<< ; GREP: Gd2O3 :: GETVAL 'GadOxG5frIso' 1 >>Gd2O3_G5f<< ; GREP: Gd2O3 :: GETVAL 'GadOxG6frIso' 1 >>Gd2O3_G6f<< ; GREP: Gd2O3 :: GETVAL 'GadOxG7frIso' 1 >>Gd2O3_G7f<< ; GREP: Gd2O3 :: GETVAL 'GadOxG8frIso' 1 >>Gd2O3_G8f<< ; GREP: Gd2O3 :: GETVAL 'GadOxG0frIso' 1 >>Gd2O3_G0f<< ; GREP: Gd2O3 :: GETVAL 'GadOxOfrIso ' 1 >>Gd2O3_Of<< ; *---- massic to isotopic ratios ---- EVALUATE Gd2O3f := Gd2O3f avo_conc * MGd2O3 / ; EVALUATE UO2Gdf := UO2Gdf avo_conc * MUO2Gd / ; EVALUATE ReNorm := Gd2O3f UO2Gdf + ; EVALUATE Gd2O3f := Gd2O3f ReNorm / ; EVALUATE UO2Gdf := UO2Gdf ReNorm / ; *---- average UO2+Gd2O3 molar mass ---- EVALUATE M_UO2Gd2O3 := Gd2O3f MGd2O3 * UO2Gdf MUO2Gd * + ; *---- U density; U & O total concentrations ---- * conc_U = conc_UO2+Gd2O3 * conc_O = 3 x conc_Gd2O3 + 2 x conc_UO2(Gd) EVALUATE UO2Gd2O3 := dGado avo_conc * M_UO2Gd2O3 / ; EVALUATE dUGado := UO2Gd2O3 MUO2Gd * avo_conc / ; EVALUATE U4eGd := U4eGd UO2Gdf * ; EVALUATE U5eGd := U5eGd UO2Gdf * ; EVALUATE U6eGd := U6eGd UO2Gdf * ; EVALUATE U8eGd := U8eGd UO2Gdf * ; EVALUATE Pu9eGd := Pu9eGd UO2Gdf * ; EVALUATE Pu0eGd := Pu0eGd UO2Gdf * ; EVALUATE Pu1eGd := Pu1eGd UO2Gdf * ; EVALUATE Pu2eGd := Pu2eGd UO2Gdf * ; EVALUATE OeUO2Gd := 2. UO2Gdf * ; EVALUATE Gd2O3_G4f := Gd2O3_G4f Gd2O3f * ; EVALUATE Gd2O3_G5f := Gd2O3_G5f Gd2O3f * ; EVALUATE Gd2O3_G6f := Gd2O3_G6f Gd2O3f * ; EVALUATE Gd2O3_G7f := Gd2O3_G7f Gd2O3f * ; EVALUATE Gd2O3_G8f := Gd2O3_G8f Gd2O3f * ; EVALUATE Gd2O3_G0f := Gd2O3_G0f Gd2O3f * ; EVALUATE Gd2O3_Of := Gd2O3_Of Gd2O3f * ; *---------------------------------------------------------------------- * REGULAR UOX FUEL *---------------------------------------------------------------------- EVALUATE enriU := Main_enri ; *EVALUATE U4enri := 1.E-15 ; *EVALUATE U4enri := enriU 100. / 0.01003 * 0.00006 - ; EVALUATE U4enri := 0.0335 100. / ; EVALUATE U5enri := enriU 100. / ; EVALUATE U6enri := 0.0061 100. / ; *EVALUATE U6enri := 1.E-15 ; EVALUATE U8enri := 1. U6enri - U5enri - U4enri - ; EVALUATE Pu8enri := 1.E-15 ; EVALUATE Pu9enri := 1.E-15 ; EVALUATE Pu0enri := 1.E-15 ; EVALUATE Pu1enri := 1.E-15 ; EVALUATE Pu2enri := 1.E-15 ; EVALUATE Am1enri := 1.E-15 ; *---- enrichment in isotopic ratio ---- EVALUATE U4enri := U4enri avo_conc * u234 / ; EVALUATE U5enri := U5enri avo_conc * u235 / ; EVALUATE U6enri := U6enri avo_conc * u236 / ; EVALUATE U8enri := U8enri avo_conc * u238 / ; EVALUATE Pu8enri := Pu8enri avo_conc * pu238 / ; EVALUATE Pu9enri := Pu9enri avo_conc * pu239 / ; EVALUATE Pu0enri := Pu0enri avo_conc * pu240 / ; EVALUATE Pu1enri := Pu1enri avo_conc * pu241 / ; EVALUATE Pu2enri := Pu2enri avo_conc * pu242 / ; EVALUATE Am1enri := Am1enri avo_conc * am241 / ; EVALUATE ReNorm := U4enri U5enri + U6enri + U8enri + Pu8enri + Pu9enri + Pu0enri + Pu1enri + Pu2enri + Am1enri + ; EVALUATE U4enri := U4enri ReNorm / ; EVALUATE U5enri := U5enri ReNorm / ; EVALUATE U6enri := U6enri ReNorm / ; EVALUATE U8enri := U8enri ReNorm / ; EVALUATE Pu8enri := Pu8enri ReNorm / ; EVALUATE Pu9enri := Pu9enri ReNorm / ; EVALUATE Pu0enri := Pu0enri ReNorm / ; EVALUATE Pu1enri := Pu1enri ReNorm / ; EVALUATE Pu2enri := Pu2enri ReNorm / ; EVALUATE Am1enri := Am1enri ReNorm / ; *---- average U, UO2 molar masses ---- EVALUATE MUPuAm := U4enri u234 * U5enri u235 * + U6enri u236 * + U8enri u238 * + Pu8enri pu238 * + Pu9enri pu239 * + Pu0enri pu240 * + Pu1enri pu241 * + Pu2enri pu242 * + Am1enri am241 * + ; EVALUATE MUO2 := MUPuAm 2. o16 * + ; *---- U density; U & 0 concentrations ---- EVALUATE U_tot := dUO2 avo_conc * MUO2 / ; ! conc_u = conc_uo2 EVALUATE O_fuel := 2. U_tot * ; ! conc_o = 2 x conc_uo2 EVALUATE dU := U_tot MUPuAm * avo_conc / ; ! U density in g/cm**3 *---------------------------------------------------------------------- * Final U+Pu+Am Concentrations in atm/cm/barn *--- Regular UOX ------------------------------------------------------ EVALUATE U4fuel := U_tot U4enri * ; EVALUATE U5fuel := U_tot U5enri * ; EVALUATE U6fuel := U_tot U6enri * ; EVALUATE U8fuel := U_tot U8enri * ; EVALUATE Pu8fuel := U_tot Pu8enri * ; EVALUATE Pu9fuel := U_tot Pu9enri * ; EVALUATE Pu0fuel := U_tot Pu0enri * ; EVALUATE Pu1fuel := U_tot Pu1enri * ; EVALUATE Pu2fuel := U_tot Pu2enri * ; EVALUATE Am1fuel := U_tot Am1enri * ; *--- UOX with Gd ----------------------------------------------------- EVALUATE U4Gd_fuel := U4eGd UO2Gd2O3 * ; EVALUATE U5Gd_fuel := U5eGd UO2Gd2O3 * ; EVALUATE U6Gd_fuel := U6eGd UO2Gd2O3 * ; EVALUATE U8Gd_fuel := U8eGd UO2Gd2O3 * ; EVALUATE Pu9Gd_fuel := Pu9eGd UO2Gd2O3 * ; EVALUATE Pu0Gd_fuel := Pu0eGd UO2Gd2O3 * ; EVALUATE Pu1Gd_fuel := Pu1eGd UO2Gd2O3 * ; EVALUATE Pu2Gd_fuel := Pu2eGd UO2Gd2O3 * ; EVALUATE Gd154_fuel := Gd2O3_G4f UO2Gd2O3 * ; EVALUATE Gd155_fuel := Gd2O3_G5f UO2Gd2O3 * ; EVALUATE Gd156_fuel := Gd2O3_G6f UO2Gd2O3 * ; EVALUATE Gd157_fuel := Gd2O3_G7f UO2Gd2O3 * ; EVALUATE Gd158_fuel := Gd2O3_G8f UO2Gd2O3 * ; EVALUATE Gd160_fuel := Gd2O3_G0f UO2Gd2O3 * ; EVALUATE O_UGd_fuel := OeUO2Gd Gd2O3_Of + UO2Gd2O3 * ; *---------------------------------------------------------------------- * Power conversion factor for UOX+Gd: * PowMWT = PowMWth / [ ( dU * TotFuelVol * f_NoGdPin ) + * (MW/T) (MW) ( dUGado * TotFuelVol * f_GadoPin ) ] * 1.E-6 * (g/cm**3) (cm**3) (unitless) (T/g) *---------------------------------------------------------------------- EVALUATE PowMWT := PowMWth dU TotFuelVol * f_NoGdPin * dUGado TotFuelVol * f_GadoPin * + 1.E-6 * / ; *---------------------------------------------------------------------- ELSE ! Input Error ECHO "Error --- AssType:" AssType "not recognized." ; ABORT: ; ENDIF ; *---------------------------------------------------------------------- :: <> ; *------ MODETI ------ REAL MBNat B10f B11f ; REAL BTot ; REAL H_mti O_mti H2O_mti B10_mti B11_mti ; GREP: BNat :: GETVAL 'AvgMolM ' 1 >>MBNat<< ; GREP: BNat :: GETVAL 'B10frIso' 1 >>B10f<< ; GREP: BNat :: GETVAL 'B11frIso' 1 >>B11f<< ; EVALUATE BTot := dmod cbor * avo_conc * MBNat / ; REAL MH2O := 18.01528 ; EVALUATE H2O_mti := 0. ; EVALUATE H_mti := dmod avo_conc * MH2O / 2. * ; EVALUATE O_mti := dmod avo_conc * MH2O / ; EVALUATE B10_mti := BTot B10f * ; EVALUATE B11_mti := BTot B11f * ; *ECHO "CONC MODETI" ; *ECHO "H2O_mti" H2O_mti ; *ECHO "H_mti " H_mti ; *ECHO "O_mti " O_mti ; *ECHO "B10_mti" B10_mti ; *ECHO "B11_mti" B11_mti ; *------ TUBE ------ *---- Made of Zircalloy4 ---- REAL MZ4 ; REAL ZrNatZ4f Zr90Z4f Zr91Z4f Zr92Z4f Zr93Z4f Zr94Z4f Zr95Z4f Zr96Z4f ; REAL Fe54Z4f Fe56Z4f Fe57Z4f Fe58Z4f ; REAL Cr50Z4f Cr52Z4f Cr53Z4f Cr54Z4f ; REAL Sn112Z4f Sn114Z4f Sn115Z4f Sn116Z4f Sn117Z4f Sn118Z4f Sn119Z4f Sn120Z4f Sn122Z4f Sn124Z4f ; REAL O16Z4f ; *---- Isotopic Ratios ---- GREP: Zirc4 :: GETVAL 'AvgMolM ' 1 >>MZ4<< ; GREP: Zirc4 :: GETVAL 'O16Z4frIso ' 1 >>O16Z4f<< ; GREP: Zirc4 :: GETVAL 'Fe54Z4frIso ' 1 >>Fe54Z4f<< ; GREP: Zirc4 :: GETVAL 'Fe56Z4frIso ' 1 >>Fe56Z4f<< ; GREP: Zirc4 :: GETVAL 'Fe57Z4frIso ' 1 >>Fe57Z4f<< ; GREP: Zirc4 :: GETVAL 'Fe58Z4frIso ' 1 >>Fe58Z4f<< ; GREP: Zirc4 :: GETVAL 'Cr50Z4frIso ' 1 >>Cr50Z4f<< ; GREP: Zirc4 :: GETVAL 'Cr52Z4frIso ' 1 >>Cr52Z4f<< ; GREP: Zirc4 :: GETVAL 'Cr53Z4frIso ' 1 >>Cr53Z4f<< ; GREP: Zirc4 :: GETVAL 'Cr54Z4frIso ' 1 >>Cr54Z4f<< ; GREP: Zirc4 :: GETVAL 'ZrNatfrIso ' 1 >>ZrNatZ4f<< ; GREP: Zirc4 :: GETVAL 'Zr90Z4frIso ' 1 >>Zr90Z4f<< ; GREP: Zirc4 :: GETVAL 'Zr91Z4frIso ' 1 >>Zr91Z4f<< ; GREP: Zirc4 :: GETVAL 'Zr92Z4frIso ' 1 >>Zr92Z4f<< ; GREP: Zirc4 :: GETVAL 'Zr93Z4frIso ' 1 >>Zr93Z4f<< ; GREP: Zirc4 :: GETVAL 'Zr94Z4frIso ' 1 >>Zr94Z4f<< ; GREP: Zirc4 :: GETVAL 'Zr95Z4frIso ' 1 >>Zr95Z4f<< ; GREP: Zirc4 :: GETVAL 'Zr96Z4frIso ' 1 >>Zr96Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn112Z4frIso' 1 >>Sn112Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn114Z4frIso' 1 >>Sn114Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn115Z4frIso' 1 >>Sn115Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn116Z4frIso' 1 >>Sn116Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn117Z4frIso' 1 >>Sn117Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn118Z4frIso' 1 >>Sn118Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn119Z4frIso' 1 >>Sn119Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn120Z4frIso' 1 >>Sn120Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn122Z4frIso' 1 >>Sn122Z4f<< ; GREP: Zirc4 :: GETVAL 'Sn124Z4frIso' 1 >>Sn124Z4f<< ; *---- Atomic concentrations ---- REAL Tubconc := dTube avo_conc * MZ4 / ; REAL Fe54T := Fe54Z4f Tubconc * ; REAL Fe56T := Fe56Z4f Tubconc * ; REAL Fe57T := Fe57Z4f Tubconc * ; REAL Fe58T := Fe58Z4f Tubconc * ; REAL Cr50T := Cr50Z4f Tubconc * ; REAL Cr52T := Cr52Z4f Tubconc * ; REAL Cr53T := Cr53Z4f Tubconc * ; REAL Cr54T := Cr54Z4f Tubconc * ; REAL O16T := O16Z4f Tubconc * ; REAL ZrNatT := ZrNatZ4f Tubconc * ; REAL Zr90T := Zr90Z4f Tubconc * ; REAL Zr91T := Zr91Z4f Tubconc * ; REAL Zr92T := Zr92Z4f Tubconc * ; REAL Zr93T := Zr93Z4f Tubconc * ; REAL Zr94T := Zr94Z4f Tubconc * ; REAL Zr95T := Zr95Z4f Tubconc * ; REAL Zr96T := Zr96Z4f Tubconc * ; REAL Sn112T := Sn112Z4f Tubconc * ; REAL Sn114T := Sn114Z4f Tubconc * ; REAL Sn115T := Sn115Z4f Tubconc * ; REAL Sn116T := Sn116Z4f Tubconc * ; REAL Sn117T := Sn117Z4f Tubconc * ; REAL Sn118T := Sn118Z4f Tubconc * ; REAL Sn119T := Sn119Z4f Tubconc * ; REAL Sn120T := Sn120Z4f Tubconc * ; REAL Sn122T := Sn122Z4f Tubconc * ; REAL Sn124T := Sn124Z4f Tubconc * ; *ECHO "CONC TUBE" ; *ECHO "Fe54T " Fe54T ; *ECHO "Fe56T " Fe56T ; *ECHO "Fe57T " Fe57T ; *ECHO "Fe58T " Fe58T ; *ECHO "Cr50T " Cr50T ; *ECHO "Cr52T " Cr52T ; *ECHO "Cr53T " Cr53T ; *ECHO "Cr54T " Cr54T ; *ECHO "O16T " O16T ; *ECHO "ZrNatT " ZrNatT ; *ECHO "Zr90T " Zr90T ; *ECHO "Zr91T " Zr91T ; *ECHO "Zr92T " Zr92T ; *ECHO "Zr93T " Zr93T ; *ECHO "Zr94T " Zr94T ; *ECHO "Zr95T " Zr95T ; *ECHO "Zr96T " Zr96T ; *ECHO "Sn112T " Sn112T ; *ECHO "Sn114T " Sn114T ; *ECHO "Sn115T " Sn115T ; *ECHO "Sn116T " Sn116T ; *ECHO "Sn117T " Sn117T ; *ECHO "Sn118T " Sn118T ; *ECHO "Sn119T " Sn119T ; *ECHO "Sn120T " Sn120T ; *ECHO "Sn122T " Sn122T ; *ECHO "Sn124T " Sn124T ; *------ GAINC ------ *---- made of Zircalloy4 (with fuel-clad gap) ---- REAL Gaiconc := dClad avo_conc * MZ4 / ; REAL Fe54Cl := Fe54Z4f Gaiconc * ; REAL Fe56Cl := Fe56Z4f Gaiconc * ; REAL Fe57Cl := Fe57Z4f Gaiconc * ; REAL Fe58Cl := Fe58Z4f Gaiconc * ; REAL Cr50Cl := Cr50Z4f Gaiconc * ; REAL Cr52Cl := Cr52Z4f Gaiconc * ; REAL Cr53Cl := Cr53Z4f Gaiconc * ; REAL Cr54Cl := Cr54Z4f Gaiconc * ; REAL O16Cl := O16Z4f Gaiconc * ; REAL ZrNatCl := ZrNatZ4f Gaiconc * ; REAL Zr90Cl := Zr90Z4f Gaiconc * ; REAL Zr91Cl := Zr91Z4f Gaiconc * ; REAL Zr92Cl := Zr92Z4f Gaiconc * ; REAL Zr93Cl := Zr93Z4f Gaiconc * ; REAL Zr94Cl := Zr94Z4f Gaiconc * ; REAL Zr95Cl := Zr95Z4f Gaiconc * ; REAL Zr96Cl := Zr96Z4f Gaiconc * ; REAL Sn112Cl := Sn112Z4f Gaiconc * ; REAL Sn114Cl := Sn114Z4f Gaiconc * ; REAL Sn115Cl := Sn115Z4f Gaiconc * ; REAL Sn116Cl := Sn116Z4f Gaiconc * ; REAL Sn117Cl := Sn117Z4f Gaiconc * ; REAL Sn118Cl := Sn118Z4f Gaiconc * ; REAL Sn119Cl := Sn119Z4f Gaiconc * ; REAL Sn120Cl := Sn120Z4f Gaiconc * ; REAL Sn122Cl := Sn122Z4f Gaiconc * ; REAL Sn124Cl := Sn124Z4f Gaiconc * ; *ECHO "CONC GAINC" ; *ECHO "Fe54Cl " Fe54Cl ; *ECHO "Fe56Cl " Fe56Cl ; *ECHO "Fe57Cl " Fe57Cl ; *ECHO "Fe58Cl " Fe58Cl ; *ECHO "Cr50Cl " Cr50Cl ; *ECHO "Cr52Cl " Cr52Cl ; *ECHO "Cr53Cl " Cr53Cl ; *ECHO "Cr54Cl " Cr54Cl ; *ECHO "O16Cl " O16Cl ; *ECHO "ZrNatCl " ZrNatCl ; *ECHO "Zr90Cl " Zr90Cl ; *ECHO "Zr91Cl " Zr91Cl ; *ECHO "Zr92Cl " Zr92Cl ; *ECHO "Zr93Cl " Zr93Cl ; *ECHO "Zr94Cl " Zr94Cl ; *ECHO "Zr95Cl " Zr95Cl ; *ECHO "Zr96Cl " Zr96Cl ; *ECHO "Sn112Cl " Sn112Cl ; *ECHO "Sn114Cl " Sn114Cl ; *ECHO "Sn115Cl " Sn115Cl ; *ECHO "Sn116Cl " Sn116Cl ; *ECHO "Sn117Cl " Sn117Cl ; *ECHO "Sn118Cl " Sn118Cl ; *ECHO "Sn119Cl " Sn119Cl ; *ECHO "Sn120Cl " Sn120Cl ; *ECHO "Sn122Cl " Sn122Cl ; *ECHO "Sn124Cl " Sn124Cl ; *------ GAINC ------ *---- made of M5 (with fuel-clad gap) ---- REAL MM5 ; REAL Fe54M5f Fe56M5f Fe57M5f Fe58M5f ; REAL O16M5f ZrM5f Nb93M5f ; REAL Zr90M5f Zr91M5f Zr92M5f Zr93M5f Zr94M5f Zr95M5f Zr96M5f ; REAL S32M5f S33M5f S34M5f ; *---- Isotopic Ratios ---- GREP: M5 :: GETVAL 'AvgMolM ' 1 >>MM5<< ; GREP: M5 :: GETVAL 'Fe54M5frIso' 1 >>Fe54M5f<< ; GREP: M5 :: GETVAL 'Fe56M5frIso' 1 >>Fe56M5f<< ; GREP: M5 :: GETVAL 'Fe57M5frIso' 1 >>Fe57M5f<< ; GREP: M5 :: GETVAL 'Fe58M5frIso' 1 >>Fe58M5f<< ; GREP: M5 :: GETVAL 'O16M5frIso ' 1 >>O16M5f<< ; GREP: M5 :: GETVAL 'ZrM5frIso ' 1 >>ZrM5f<< ; GREP: M5 :: GETVAL 'Zr90M5frIso' 1 >>Zr90M5f<< ; GREP: M5 :: GETVAL 'Zr91M5frIso' 1 >>Zr91M5f<< ; GREP: M5 :: GETVAL 'Zr92M5frIso' 1 >>Zr92M5f<< ; GREP: M5 :: GETVAL 'Zr93M5frIso' 1 >>Zr93M5f<< ; GREP: M5 :: GETVAL 'Zr94M5frIso' 1 >>Zr94M5f<< ; GREP: M5 :: GETVAL 'Zr95M5frIso' 1 >>Zr95M5f<< ; GREP: M5 :: GETVAL 'Zr96M5frIso' 1 >>Zr96M5f<< ; GREP: M5 :: GETVAL 'Nb93M5frIso' 1 >>Nb93M5f<< ; GREP: M5 :: GETVAL 'S32M5frIso ' 1 >>S32M5f<< ; GREP: M5 :: GETVAL 'S33M5frIso ' 1 >>S33M5f<< ; GREP: M5 :: GETVAL 'S34M5frIso ' 1 >>S34M5f<< ; REAL GaiconcM5 := dClad avo_conc * MM5 / ; REAL Fe54M5Cl := Fe54M5f GaiconcM5 * ; REAL Fe56M5Cl := Fe56M5f GaiconcM5 * ; REAL Fe57M5Cl := Fe57M5f GaiconcM5 * ; REAL Fe58M5Cl := Fe58M5f GaiconcM5 * ; REAL O16M5Cl := O16M5f GaiconcM5 * ; REAL ZrM5Cl := ZrM5f GaiconcM5 * ; REAL Zr90M5Cl := Zr90M5f GaiconcM5 * ; REAL Zr91M5Cl := Zr91M5f GaiconcM5 * ; REAL Zr92M5Cl := Zr92M5f GaiconcM5 * ; REAL Zr93M5Cl := Zr93M5f GaiconcM5 * ; REAL Zr94M5Cl := Zr94M5f GaiconcM5 * ; REAL Zr95M5Cl := Zr95M5f GaiconcM5 * ; REAL Zr96M5Cl := Zr96M5f GaiconcM5 * ; REAL Nb93M5Cl := Nb93M5f GaiconcM5 * ; REAL S32M5Cl := S32M5f GaiconcM5 * ; REAL S33M5Cl := S33M5f GaiconcM5 * ; REAL S34M5Cl := S34M5f GaiconcM5 * ; *ECHO "CONC GAINC" ; *ECHO "Fe54M5Cl " Fe54M5Cl ; *ECHO "Fe56M5Cl " Fe56M5Cl ; *ECHO "Fe57M5Cl " Fe57M5Cl ; *ECHO "Fe58M5Cl " Fe58M5Cl ; *ECHO "O16M5Cl " O16M5Cl ; *ECHO "ZrM5Cl " ZrM5Cl ; *ECHO "Zr90M5Cl " Zr90M5Cl ; *ECHO "Zr91M5Cl " Zr91M5Cl ; *ECHO "Zr92M5Cl " Zr92M5Cl ; *ECHO "Zr93M5Cl " Zr93M5Cl ; *ECHO "Zr94M5Cl " Zr94M5Cl ; *ECHO "Zr95M5Cl " Zr95M5Cl ; *ECHO "Zr96M5Cl " Zr96M5Cl ; *ECHO "Nb93M5Cl " Nb93M5Cl ; *ECHO "S32M5Cl " S32M5Cl ; *ECHO "S33M5Cl " S33M5Cl ; *ECHO "S34M5Cl " S34M5Cl ; *------ GRILLE ------ REAL MSS304 ; REAL Fe54SSf Fe56SSf Fe57SSf Fe58SSf ; REAL Ni58SSf Ni60SSf Ni61SSf Ni62SSf Ni64SSf ; REAL Cr50SSf Cr52SSf Cr53SSf Cr54SSf ; REAL MnSSf SiSSf Si28SSf Si29SSf Si30SSf ; GREP: SS304 :: GETVAL 'AvgMolM ' 1 >>MSS304<< ; GREP: SS304 :: GETVAL 'Fe54SSfrIso ' 1 >>Fe54SSf<< ; GREP: SS304 :: GETVAL 'Fe56SSfrIso ' 1 >>Fe56SSf<< ; GREP: SS304 :: GETVAL 'Fe57SSfrIso ' 1 >>Fe57SSf<< ; GREP: SS304 :: GETVAL 'Fe58SSfrIso ' 1 >>Fe58SSf<< ; GREP: SS304 :: GETVAL 'Ni58SSfrIso ' 1 >>Ni58SSf<< ; GREP: SS304 :: GETVAL 'Ni60SSfrIso ' 1 >>Ni60SSf<< ; GREP: SS304 :: GETVAL 'Ni61SSfrIso ' 1 >>Ni61SSf<< ; GREP: SS304 :: GETVAL 'Ni62SSfrIso ' 1 >>Ni62SSf<< ; GREP: SS304 :: GETVAL 'Ni64SSfrIso ' 1 >>Ni64SSf<< ; GREP: SS304 :: GETVAL 'Cr50SSfrIso ' 1 >>Cr50SSf<< ; GREP: SS304 :: GETVAL 'Cr52SSfrIso ' 1 >>Cr52SSf<< ; GREP: SS304 :: GETVAL 'Cr53SSfrIso ' 1 >>Cr53SSf<< ; GREP: SS304 :: GETVAL 'Cr54SSfrIso ' 1 >>Cr54SSf<< ; GREP: SS304 :: GETVAL 'Mn55SSfrIso ' 1 >>MnSSf<< ; GREP: SS304 :: GETVAL 'SiNatSSfrIso' 1 >>SiSSf<< ; GREP: SS304 :: GETVAL 'Si28SSfrIso ' 1 >>Si28SSf<< ; GREP: SS304 :: GETVAL 'Si29SSfrIso ' 1 >>Si29SSf<< ; GREP: SS304 :: GETVAL 'Si30SSfrIso ' 1 >>Si30SSf<< ; REAL MInc ; REAL Ni58Inf Ni60Inf Ni61Inf Ni62Inf Ni64Inf ; REAL Cr50Inf Cr52Inf Cr53Inf Cr54Inf ; REAL Fe54Inf Fe56Inf Fe57Inf Fe58Inf ; REAL B10BInf B11BInf ; REAL TiNatInf Ti46Inf Ti47Inf Ti48Inf Ti49Inf Ti50Inf ; REAL SiNatInf Si28Inf Si29Inf Si30Inf ; REAL Al27Inf Mo95Inf Co59Inf Mn55Inf Nb93Inf ; GREP: Incon :: GETVAL 'AvgMolM ' 1 >>MInc<< ; GREP: Incon :: GETVAL 'Ni58InfrIso' 1 >>Ni58Inf<< ; GREP: Incon :: GETVAL 'Ni60InfrIso' 1 >>Ni60Inf<< ; GREP: Incon :: GETVAL 'Ni61InfrIso' 1 >>Ni61Inf<< ; GREP: Incon :: GETVAL 'Ni62InfrIso' 1 >>Ni62Inf<< ; GREP: Incon :: GETVAL 'Ni64InfrIso' 1 >>Ni64Inf<< ; GREP: Incon :: GETVAL 'Cr50InfrIso' 1 >>Cr50Inf<< ; GREP: Incon :: GETVAL 'Cr52InfrIso' 1 >>Cr52Inf<< ; GREP: Incon :: GETVAL 'Cr53InfrIso' 1 >>Cr53Inf<< ; GREP: Incon :: GETVAL 'Cr54InfrIso' 1 >>Cr54Inf<< ; GREP: Incon :: GETVAL 'Fe54InfrIso' 1 >>Fe54Inf<< ; GREP: Incon :: GETVAL 'Fe56InfrIso' 1 >>Fe56Inf<< ; GREP: Incon :: GETVAL 'Fe57InfrIso' 1 >>Fe57Inf<< ; GREP: Incon :: GETVAL 'Fe58InfrIso' 1 >>Fe58Inf<< ; GREP: Incon :: GETVAL 'B10InfrIso ' 1 >>B10BInf<< ; GREP: Incon :: GETVAL 'B11InfrIso ' 1 >>B11BInf<< ; GREP: Incon :: GETVAL 'Al27InfrIso' 1 >>Al27Inf<< ; GREP: Incon :: GETVAL 'Mo95InfrIso' 1 >>Mo95Inf<< ; GREP: Incon :: GETVAL 'TiNatInfrIso' 1 >>TiNatInf<< ; GREP: Incon :: GETVAL 'Ti46InfrIso ' 1 >>Ti46Inf<< ; GREP: Incon :: GETVAL 'Ti47InfrIso ' 1 >>Ti47Inf<< ; GREP: Incon :: GETVAL 'Ti48InfrIso ' 1 >>Ti48Inf<< ; GREP: Incon :: GETVAL 'Ti49InfrIso ' 1 >>Ti49Inf<< ; GREP: Incon :: GETVAL 'Ti50InfrIso ' 1 >>Ti50Inf<< ; GREP: Incon :: GETVAL 'SiNatInfrIso' 1 >>SiNatInf<< ; GREP: Incon :: GETVAL 'Si28InfrIso ' 1 >>Si28Inf<< ; GREP: Incon :: GETVAL 'Si29InfrIso ' 1 >>Si29Inf<< ; GREP: Incon :: GETVAL 'Si30InfrIso ' 1 >>Si30Inf<< ; GREP: Incon :: GETVAL 'Co59InfrIso' 1 >>Co59Inf<< ; GREP: Incon :: GETVAL 'Mn55InfrIso' 1 >>Mn55Inf<< ; GREP: Incon :: GETVAL 'Nb93InfrIso' 1 >>Nb93Inf<< ; *---- Grille Isotopic mass ratios ---- EVALUATE InGridf := InGridf avo_conc * MInc / ; EVALUATE Z4Gridf := Z4Gridf avo_conc * MZ4 / ; EVALUATE SSGridf := SSGridf avo_conc * MSS304 / ; EVALUATE ReNorm := InGridf Z4Gridf SSGridf + + ; EVALUATE InGridf := InGridf ReNorm / ; EVALUATE Z4Gridf := Z4Gridf ReNorm / ; EVALUATE SSGridf := SSGridf ReNorm / ; *---- Grille Atomic Concentration ---- REAL MGrid := InGridf MInc * Z4Gridf MZ4 * + SSGridf MSS304 * + ; REAL Griconc := dGrid avo_conc * MGrid / ; REAL Ni58Gr := Ni58Inf InGridf * Ni58SSf SSGridf * + Griconc * ; REAL Ni60Gr := Ni60Inf InGridf * Ni60SSf SSGridf * + Griconc * ; REAL Ni61Gr := Ni61Inf InGridf * Ni61SSf SSGridf * + Griconc * ; REAL Ni62Gr := Ni62Inf InGridf * Ni62SSf SSGridf * + Griconc * ; REAL Ni64Gr := Ni64Inf InGridf * Ni64SSf SSGridf * + Griconc * ; REAL Cr50Gr := Cr50Inf InGridf * Cr50Z4f Z4Gridf * + Cr50SSf SSGridf * + Griconc * ; REAL Cr52Gr := Cr52Inf InGridf * Cr52Z4f Z4Gridf * + Cr52SSf SSGridf * + Griconc * ; REAL Cr53Gr := Cr53Inf InGridf * Cr53Z4f Z4Gridf * + Cr53SSf SSGridf * + Griconc * ; REAL Cr54Gr := Cr54Inf InGridf * Cr54Z4f Z4Gridf * + Cr54SSf SSGridf * + Griconc * ; REAL Fe54Gr := Fe54Inf InGridf * Fe54Z4f Z4Gridf * + Fe54SSf SSGridf * + Griconc * ; REAL Fe56Gr := Fe56Inf InGridf * Fe56Z4f Z4Gridf * + Fe56SSf SSGridf * + Griconc * ; REAL Fe57Gr := Fe57Inf InGridf * Fe57Z4f Z4Gridf * + Fe57SSf SSGridf * + Griconc * ; REAL Fe58Gr := Fe58Inf InGridf * Fe58Z4f Z4Gridf * + Fe58SSf SSGridf * + Griconc * ; REAL B10BGr := B10BInf InGridf * Griconc * ; REAL B11BGr := B11BInf InGridf * Griconc * ; REAL Al27Gr := Al27Inf InGridf * Griconc * ; REAL Mo95Gr := Mo95Inf InGridf * Griconc * ; REAL TiNatGr := TiNatInf InGridf * Griconc * ; REAL Ti46Gr := Ti46Inf InGridf * Griconc * ; REAL Ti47Gr := Ti47Inf InGridf * Griconc * ; REAL Ti48Gr := Ti48Inf InGridf * Griconc * ; REAL Ti49Gr := Ti49Inf InGridf * Griconc * ; REAL Ti50Gr := Ti50Inf InGridf * Griconc * ; REAL SiNatGr := SiNatInf InGridf * SiSSf SSGridf * + Griconc * ; REAL Si28Gr := Si28Inf InGridf * Si28SSf SSGridf * + Griconc * ; REAL Si29Gr := Si29Inf InGridf * Si29SSf SSGridf * + Griconc * ; REAL Si30Gr := Si30Inf InGridf * Si30SSf SSGridf * + Griconc * ; REAL Co59Gr := Co59Inf InGridf * Griconc * ; REAL Mn55Gr := Mn55Inf InGridf * MnSSf SSGridf * + Griconc * ; REAL Nb93Gr := Nb93Inf InGridf * Griconc * ; REAL Sn112Gr := Sn112Z4f Z4Gridf * Griconc * ; REAL Sn114Gr := Sn114Z4f Z4Gridf * Griconc * ; REAL Sn115Gr := Sn115Z4f Z4Gridf * Griconc * ; REAL Sn116Gr := Sn116Z4f Z4Gridf * Griconc * ; REAL Sn117Gr := Sn117Z4f Z4Gridf * Griconc * ; REAL Sn118Gr := Sn118Z4f Z4Gridf * Griconc * ; REAL Sn119Gr := Sn119Z4f Z4Gridf * Griconc * ; REAL Sn120Gr := Sn120Z4f Z4Gridf * Griconc * ; REAL Sn122Gr := Sn122Z4f Z4Gridf * Griconc * ; REAL Sn124Gr := Sn124Z4f Z4Gridf * Griconc * ; REAL O16Gr := O16Z4f Z4Gridf * Griconc * ; REAL ZrNatGr := ZrNatZ4f Z4Gridf * Griconc * ; REAL Zr90Gr := Zr90Z4f Z4Gridf * Griconc * ; REAL Zr91Gr := Zr91Z4f Z4Gridf * Griconc * ; REAL Zr92Gr := Zr92Z4f Z4Gridf * Griconc * ; REAL Zr93Gr := Zr93Z4f Z4Gridf * Griconc * ; REAL Zr94Gr := Zr94Z4f Z4Gridf * Griconc * ; REAL Zr95Gr := Zr95Z4f Z4Gridf * Griconc * ; REAL Zr96Gr := Zr96Z4f Z4Gridf * Griconc * ; *ECHO "CONC GRILLE" ; *ECHO "Ni58Gr " Ni58Gr ; *ECHO "Ni60Gr " Ni60Gr ; *ECHO "Ni61Gr " Ni61Gr ; *ECHO "Ni62Gr " Ni62Gr ; *ECHO "Ni64Gr " Ni64Gr ; *ECHO "Cr50Gr " Cr50Gr ; *ECHO "Cr52Gr " Cr52Gr ; *ECHO "Cr53Gr " Cr53Gr ; *ECHO "Cr54Gr " Cr54Gr ; *ECHO "Fe54Gr " Fe54Gr ; *ECHO "Fe56Gr " Fe56Gr ; *ECHO "Fe57Gr " Fe57Gr ; *ECHO "Fe58Gr " Fe58Gr ; *ECHO "B10BGr " B10BGr ; *ECHO "B11BGr " B11BGr ; *ECHO "Al27Gr " Al27Gr ; *ECHO "Mo95Gr " Mo95Gr ; *ECHO "TiNatGr " TiNatGr ; *ECHO "Ti46Gr " Ti46Gr ; *ECHO "Ti47Gr " Ti47Gr ; *ECHO "Ti48Gr " Ti48Gr ; *ECHO "Ti49Gr " Ti49Gr ; *ECHO "Ti50Gr " Ti50Gr ; *ECHO "SiNatGr " SiNatGr ; *ECHO "Si28Gr " Si28Gr ; *ECHO "Si29Gr " Si29Gr ; *ECHO "Si30Gr " Si30Gr ; *ECHO "Co59Gr " Co59Gr ; *ECHO "Mn55Gr " Mn55Gr ; *ECHO "Nb93Gr " Nb93Gr ; *ECHO "Sn112Gr" Sn112Gr ; *ECHO "Sn114Gr" Sn114Gr ; *ECHO "Sn115Gr" Sn115Gr ; *ECHO "Sn116Gr" Sn116Gr ; *ECHO "Sn117Gr" Sn117Gr ; *ECHO "Sn118Gr" Sn118Gr ; *ECHO "Sn119Gr" Sn119Gr ; *ECHO "Sn120Gr" Sn120Gr ; *ECHO "Sn122Gr" Sn122Gr ; *ECHO "Sn124Gr" Sn124Gr ; *ECHO "O16Gr " O16Gr ; *ECHO "ZrNatGr " ZrNatGr ; *ECHO "Zr90Gr " Zr90Gr ; *ECHO "Zr91Gr " Zr91Gr ; *ECHO "Zr92Gr " Zr92Gr ; *ECHO "Zr93Gr " Zr93Gr ; *ECHO "Zr94Gr " Zr94Gr ; *ECHO "Zr95Gr " Zr95Gr ; *ECHO "Zr96Gr " Zr96Gr ; *------ Composants differencies des grilles : Zr4 ------ REAL GrConcZr := dZr4 avo_conc * MZ4 / ; REAL Fe54Z4 := Fe54Z4f GrConcZr * ; REAL Fe56Z4 := Fe56Z4f GrConcZr * ; REAL Fe57Z4 := Fe57Z4f GrConcZr * ; REAL Fe58Z4 := Fe58Z4f GrConcZr * ; REAL Cr50Z4 := Cr50Z4f GrConcZr * ; REAL Cr52Z4 := Cr52Z4f GrConcZr * ; REAL Cr53Z4 := Cr53Z4f GrConcZr * ; REAL Cr54Z4 := Cr54Z4f GrConcZr * ; REAL O16Z4 := O16Z4f GrConcZr * ; REAL ZrNatZ4 := ZrNatZ4f GrConcZr * ; REAL Zr90Z4 := Zr90Z4f GrConcZr * ; REAL Zr91Z4 := Zr91Z4f GrConcZr * ; REAL Zr92Z4 := Zr92Z4f GrConcZr * ; REAL Zr93Z4 := Zr93Z4f GrConcZr * ; REAL Zr94Z4 := Zr94Z4f GrConcZr * ; REAL Zr95Z4 := Zr95Z4f GrConcZr * ; REAL Zr96Z4 := Zr96Z4f GrConcZr * ; REAL Sn112Z4 := Sn112Z4f GrConcZr * ; REAL Sn114Z4 := Sn114Z4f GrConcZr * ; REAL Sn115Z4 := Sn115Z4f GrConcZr * ; REAL Sn116Z4 := Sn116Z4f GrConcZr * ; REAL Sn117Z4 := Sn117Z4f GrConcZr * ; REAL Sn118Z4 := Sn118Z4f GrConcZr * ; REAL Sn119Z4 := Sn119Z4f GrConcZr * ; REAL Sn120Z4 := Sn120Z4f GrConcZr * ; REAL Sn122Z4 := Sn122Z4f GrConcZr * ; REAL Sn124Z4 := Sn124Z4f GrConcZr * ; *------ Composants differencies des grilles : Inconel ------ REAL GrConcInc := dInc avo_conc * MInc / ; REAL Ni58Inc := Ni58Inf GrConcInc * ; REAL Ni60Inc := Ni60Inf GrConcInc * ; REAL Ni61Inc := Ni61Inf GrConcInc * ; REAL Ni62Inc := Ni62Inf GrConcInc * ; REAL Ni64Inc := Ni64Inf GrConcInc * ; REAL Cr50Inc := Cr50Inf GrConcInc * ; REAL Cr52Inc := Cr52Inf GrConcInc * ; REAL Cr53Inc := Cr53Inf GrConcInc * ; REAL Cr54Inc := Cr54Inf GrConcInc * ; REAL Fe54Inc := Fe54Inf GrConcInc * ; REAL Fe56Inc := Fe56Inf GrConcInc * ; REAL Fe57Inc := Fe57Inf GrConcInc * ; REAL Fe58Inc := Fe58Inf GrConcInc * ; REAL B10BInc := B10BInf GrConcInc * ; REAL B11BInc := B11BInf GrConcInc * ; REAL Al27Inc := Al27Inf GrConcInc * ; REAL Mo95Inc := Mo95Inf GrConcInc * ; REAL TiNatInc := TiNatInf GrConcInc * ; REAL Ti46Inc := Ti46Inf GrConcInc * ; REAL Ti47Inc := Ti47Inf GrConcInc * ; REAL Ti48Inc := Ti48Inf GrConcInc * ; REAL Ti49Inc := Ti49Inf GrConcInc * ; REAL Ti50Inc := Ti50Inf GrConcInc * ; REAL SiNatInc := SiNatInf GrConcInc * ; REAL Si28Inc := Si28Inf GrConcInc * ; REAL Si29Inc := Si29Inf GrConcInc * ; REAL Si30Inc := Si30Inf GrConcInc * ; REAL Co59Inc := Co59Inf GrConcInc * ; REAL Mn55Inc := Mn55Inf GrConcInc * ; REAL Nb93Inc := Nb93Inf GrConcInc * ; *------ GAINA, GAINAA, ABS_ACIER ------ *---- made of SS304 Stainless Steel (with rod-clad gap) ---- REAL SSconc := dSS avo_conc * MSS304 / ; REAL SiSS := SiSSf SSconc * ; REAL Si28SS := Si28SSf SSconc * ; REAL Si29SS := Si29SSf SSconc * ; REAL Si30SS := Si30SSf SSconc * ; REAL Fe54SS := Fe54SSf SSconc * ; REAL Fe56SS := Fe56SSf SSconc * ; REAL Fe57SS := Fe57SSf SSconc * ; REAL Fe58SS := Fe58SSf SSconc * ; REAL Ni58SS := Ni58SSf SSconc * ; REAL Ni60SS := Ni60SSf SSconc * ; REAL Ni61SS := Ni61SSf SSconc * ; REAL Ni62SS := Ni62SSf SSconc * ; REAL Ni64SS := Ni64SSf SSconc * ; REAL Cr50SS := Cr50SSf SSconc * ; REAL Cr52SS := Cr52SSf SSconc * ; REAL Cr53SS := Cr53SSf SSconc * ; REAL Cr54SS := Cr54SSf SSconc * ; REAL MnSS := MnSSf SSconc * ; *----- ASB_ACIER: Steel Rods ------ REAL RodSSconc := dRodSS avo_conc * MSS304 / ; REAL SiSSRo := SiSSf RodSSconc * ; REAL Si28SSRo := Si28SSf RodSSconc * ; REAL Si29SSRo := Si29SSf RodSSconc * ; REAL Si30SSRo := Si30SSf RodSSconc * ; REAL Fe54SSRo := Fe54SSf RodSSconc * ; REAL Fe56SSRo := Fe56SSf RodSSconc * ; REAL Fe57SSRo := Fe57SSf RodSSconc * ; REAL Fe58SSRo := Fe58SSf RodSSconc * ; REAL Ni58SSRo := Ni58SSf RodSSconc * ; REAL Ni60SSRo := Ni60SSf RodSSconc * ; REAL Ni61SSRo := Ni61SSf RodSSconc * ; REAL Ni62SSRo := Ni62SSf RodSSconc * ; REAL Ni64SSRo := Ni64SSf RodSSconc * ; REAL Cr50SSRo := Cr50SSf RodSSconc * ; REAL Cr52SSRo := Cr52SSf RodSSconc * ; REAL Cr53SSRo := Cr53SSf RodSSconc * ; REAL Cr54SSRo := Cr54SSf RodSSconc * ; REAL MnSSRo := MnSSf RodSSconc * ; *ECHO "----ABS_ACIER----" ; *ECHO "SiSSRo" SiSSRo ; *ECHO "Si28SSRo " Si28SSRo ; *ECHO "Si29SSRo " Si29SSRo ; *ECHO "Si30SSRo " Si30SSRo ; *ECHO "Fe54SSRo" Fe54SSRo ; *ECHO "Fe56SSRo" Fe56SSRo ; *ECHO "Fe57SSRo" Fe57SSRo ; *ECHO "Fe58SSRo" Fe58SSRo ; *ECHO "Ni58SSRo" Ni58SSRo ; *ECHO "Ni60SSRo" Ni60SSRo ; *ECHO "Ni61SSRo" Ni61SSRo ; *ECHO "Ni62SSRo" Ni62SSRo ; *ECHO "Ni64SSRo" Ni64SSRo ; *ECHO "Cr50SSRo" Cr50SSRo ; *ECHO "Cr52SSRo" Cr52SSRo ; *ECHO "Cr53SSRo" Cr53SSRo ; *ECHO "Cr54SSRo" Cr54SSRo ; *ECHO "MnSSRo " MnSSRo ; *----- GAINA & GAINAA: clad-rod gap for AIC Rods ------ REAL CladAICconc := dCladAIC avo_conc * MSS304 / ; REAL SiCladAIC := SiSSf CladAICconc * ; REAL Si28CladAIC := Si28SSf CladAICconc * ; REAL Si29CladAIC := Si29SSf CladAICconc * ; REAL Si30CladAIC := Si30SSf CladAICconc * ; REAL Fe54CladAIC := Fe54SSf CladAICconc * ; REAL Fe56CladAIC := Fe56SSf CladAICconc * ; REAL Fe57CladAIC := Fe57SSf CladAICconc * ; REAL Fe58CladAIC := Fe58SSf CladAICconc * ; REAL Ni58CladAIC := Ni58SSf CladAICconc * ; REAL Ni60CladAIC := Ni60SSf CladAICconc * ; REAL Ni61CladAIC := Ni61SSf CladAICconc * ; REAL Ni62CladAIC := Ni62SSf CladAICconc * ; REAL Ni64CladAIC := Ni64SSf CladAICconc * ; REAL Cr50CladAIC := Cr50SSf CladAICconc * ; REAL Cr52CladAIC := Cr52SSf CladAICconc * ; REAL Cr53CladAIC := Cr53SSf CladAICconc * ; REAL Cr54CladAIC := Cr54SSf CladAICconc * ; REAL MnCladAIC := MnSSf CladAICconc * ; *ECHO "----GAINA_AIC----" ; *ECHO "SiCladAIC" SiCladAIC ; *ECHO "Si28CladAIC " Si28CladAIC ; *ECHO "Si29CladAIC " Si29CladAIC ; *ECHO "Si30CladAIC " Si30CladAIC ; *ECHO "Fe54CladAIC" Fe54CladAIC ; *ECHO "Fe56CladAIC" Fe56CladAIC ; *ECHO "Fe57CladAIC" Fe57CladAIC ; *ECHO "Fe58CladAIC" Fe58CladAIC ; *ECHO "Ni58CladAIC" Ni58CladAIC ; *ECHO "Ni60CladAIC" Ni60CladAIC ; *ECHO "Ni61CladAIC" Ni61CladAIC ; *ECHO "Ni62CladAIC" Ni62CladAIC ; *ECHO "Ni64CladAIC" Ni64CladAIC ; *ECHO "Cr50CladAIC" Cr50CladAIC ; *ECHO "Cr52CladAIC" Cr52CladAIC ; *ECHO "Cr53CladAIC" Cr53CladAIC ; *ECHO "Cr54CladAIC" Cr54CladAIC ; *ECHO "MnCladAIC " MnCladAIC ; *----- GAINA & GAINAA: clad-rod gap for Steel Rods ------ REAL CladSSconc := dCladSS avo_conc * MSS304 / ; REAL SiCladSS := SiSSf CladSSconc * ; REAL Si28CladSS := Si28SSf CladSSconc * ; REAL Si29CladSS := Si29SSf CladSSconc * ; REAL Si30CladSS := Si30SSf CladSSconc * ; REAL Fe54CladSS := Fe54SSf CladSSconc * ; REAL Fe56CladSS := Fe56SSf CladSSconc * ; REAL Fe57CladSS := Fe57SSf CladSSconc * ; REAL Fe58CladSS := Fe58SSf CladSSconc * ; REAL Ni58CladSS := Ni58SSf CladSSconc * ; REAL Ni60CladSS := Ni60SSf CladSSconc * ; REAL Ni61CladSS := Ni61SSf CladSSconc * ; REAL Ni62CladSS := Ni62SSf CladSSconc * ; REAL Ni64CladSS := Ni64SSf CladSSconc * ; REAL Cr50CladSS := Cr50SSf CladSSconc * ; REAL Cr52CladSS := Cr52SSf CladSSconc * ; REAL Cr53CladSS := Cr53SSf CladSSconc * ; REAL Cr54CladSS := Cr54SSf CladSSconc * ; REAL MnCladSS := MnSSf CladSSconc * ; *ECHO "----GAINA_ABS_ACIER----" ; *ECHO "SiCladSS" SiCladSS ; *ECHO "Si28CladSS " Si28CladSS ; *ECHO "Si29CladSS " Si29CladSS ; *ECHO "Si30CladSS " Si30CladSS ; *ECHO "Fe54CladSS" Fe54CladSS ; *ECHO "Fe56CladSS" Fe56CladSS ; *ECHO "Fe57CladSS" Fe57CladSS ; *ECHO "Fe58CladSS" Fe58CladSS ; *ECHO "Ni58CladSS" Ni58CladSS ; *ECHO "Ni60CladSS" Ni60CladSS ; *ECHO "Ni61CladSS" Ni61CladSS ; *ECHO "Ni62CladSS" Ni62CladSS ; *ECHO "Ni64CladSS" Ni64CladSS ; *ECHO "Cr50CladSS" Cr50CladSS ; *ECHO "Cr52CladSS" Cr52CladSS ; *ECHO "Cr53CladSS" Cr53CladSS ; *ECHO "Cr54CladSS" Cr54CladSS ; *ECHO "MnCladSS " MnCladSS ; *----- GAINA & GAINAA: clad-rod gap for B4C Rods ------ REAL CladB4Cconc := dCladB4C avo_conc * MSS304 / ; REAL SiCladB4C := SiSSf CladB4Cconc * ; REAL Si28CladB4C := Si28SSf CladB4Cconc * ; REAL Si29CladB4C := Si29SSf CladB4Cconc * ; REAL Si30CladB4C := Si30SSf CladB4Cconc * ; REAL Fe54CladB4C := Fe54SSf CladB4Cconc * ; REAL Fe56CladB4C := Fe56SSf CladB4Cconc * ; REAL Fe57CladB4C := Fe57SSf CladB4Cconc * ; REAL Fe58CladB4C := Fe58SSf CladB4Cconc * ; REAL Ni58CladB4C := Ni58SSf CladB4Cconc * ; REAL Ni60CladB4C := Ni60SSf CladB4Cconc * ; REAL Ni61CladB4C := Ni61SSf CladB4Cconc * ; REAL Ni62CladB4C := Ni62SSf CladB4Cconc * ; REAL Ni64CladB4C := Ni64SSf CladB4Cconc * ; REAL Cr50CladB4C := Cr50SSf CladB4Cconc * ; REAL Cr52CladB4C := Cr52SSf CladB4Cconc * ; REAL Cr53CladB4C := Cr53SSf CladB4Cconc * ; REAL Cr54CladB4C := Cr54SSf CladB4Cconc * ; REAL MnCladB4C := MnSSf CladB4Cconc * ; *ECHO "----GAINA_B4C----" ; *ECHO "SiCladB4C" SiCladB4C ; *ECHO "Si28CladB4C " Si28CladB4C ; *ECHO "Si29CladB4C " Si29CladB4C ; *ECHO "Si30CladB4C " Si30CladB4C ; *ECHO "Fe54CladB4C" Fe54CladB4C ; *ECHO "Fe56CladB4C" Fe56CladB4C ; *ECHO "Fe57CladB4C" Fe57CladB4C ; *ECHO "Fe58CladB4C" Fe58CladB4C ; *ECHO "Ni58CladB4C" Ni58CladB4C ; *ECHO "Ni60CladB4C" Ni60CladB4C ; *ECHO "Ni61CladB4C" Ni61CladB4C ; *ECHO "Ni62CladB4C" Ni62CladB4C ; *ECHO "Ni64CladB4C" Ni64CladB4C ; *ECHO "Cr50CladB4C" Cr50CladB4C ; *ECHO "Cr52CladB4C" Cr52CladB4C ; *ECHO "Cr53CladB4C" Cr53CladB4C ; *ECHO "Cr54CladB4C" Cr54CladB4C ; *ECHO "MnCladB4C " MnCladB4C ; *==================================================================== * TIHANGE GAINE (BP Interior Cladding) *==================================================================== REAL CladIntconc := dCladIntBP avo_conc * MSS304 / ; *---- REAL Fe54CladInt := Fe54SSf CladIntconc * ; REAL Fe56CladInt := Fe56SSf CladIntconc * ; REAL Fe57CladInt := Fe57SSf CladIntconc * ; REAL Fe58CladInt := Fe58SSf CladIntconc * ; REAL Ni58CladInt := Ni58SSf CladIntconc * ; REAL Ni60CladInt := Ni60SSf CladIntconc * ; REAL Ni61CladInt := Ni61SSf CladIntconc * ; REAL Ni62CladInt := Ni62SSf CladIntconc * ; REAL Ni64CladInt := Ni64SSf CladIntconc * ; REAL Cr50CladInt := Cr50SSf CladIntconc * ; REAL Cr52CladInt := Cr52SSf CladIntconc * ; REAL Cr53CladInt := Cr53SSf CladIntconc * ; REAL Cr54CladInt := Cr54SSf CladIntconc * ; REAL MnCladInt := MnSSf CladIntconc * ; *---------------------------- *ECHO "----GAININT ----" ; *ECHO "Fe54CladInt" Fe54CladInt ; *ECHO "Fe56CladInt" Fe56CladInt ; *ECHO "Fe57CladInt" Fe57CladInt ; *ECHO "Fe58CladInt" Fe58CladInt ; *ECHO "Ni58CladInt" Ni58CladInt ; *ECHO "Ni60CladInt" Ni60CladInt ; *ECHO "Ni61CladInt" Ni61CladInt ; *ECHO "Ni62CladInt" Ni62CladInt ; *ECHO "Ni64CladInt" Ni64CladInt ; *ECHO "Cr50CladInt" Cr50CladInt ; *ECHO "Cr52CladInt" Cr52CladInt ; *ECHO "Cr53CladInt" Cr53CladInt ; *ECHO "Cr54CladInt" Cr54CladInt ; *ECHO "MnCladInt " MnCladInt ; *---------------------------- *==================================================================== * TIHANGE GAINE (BP Exterior Cladding) *==================================================================== REAL CladExtconc := dCladExtBP avo_conc * MSS304 / ; *---- *---- REAL Fe54CladExt := Fe54SSf CladExtconc * ; REAL Fe56CladExt := Fe56SSf CladExtconc * ; REAL Fe57CladExt := Fe57SSf CladExtconc * ; REAL Fe58CladExt := Fe58SSf CladExtconc * ; REAL Ni58CladExt := Ni58SSf CladExtconc * ; REAL Ni60CladExt := Ni60SSf CladExtconc * ; REAL Ni61CladExt := Ni61SSf CladExtconc * ; REAL Ni62CladExt := Ni62SSf CladExtconc * ; REAL Ni64CladExt := Ni64SSf CladExtconc * ; REAL Cr50CladExt := Cr50SSf CladExtconc * ; REAL Cr52CladExt := Cr52SSf CladExtconc * ; REAL Cr53CladExt := Cr53SSf CladExtconc * ; REAL Cr54CladExt := Cr54SSf CladExtconc * ; REAL MnCladExt := MnSSf CladExtconc * ; *---------------------------- *ECHO "----GAINEXT ----" ; *ECHO "Fe54CladExt" Fe54CladExt ; *ECHO "Fe56CladExt" Fe56CladExt ; *ECHO "Fe57CladExt" Fe57CladExt ; *ECHO "Fe58CladExt" Fe58CladExt ; *ECHO "Ni58CladExt" Ni58CladExt ; *ECHO "Ni60CladExt" Ni60CladExt ; *ECHO "Ni61CladExt" Ni61CladExt ; *ECHO "Ni62CladExt" Ni62CladExt ; *ECHO "Ni64CladExt" Ni64CladExt ; *ECHO "Cr50CladExt" Cr50CladExt ; *ECHO "Cr52CladExt" Cr52CladExt ; *ECHO "Cr53CladExt" Cr53CladExt ; *ECHO "Cr54CladExt" Cr54CladExt ; *ECHO "MnCladExt " MnCladExt ; *---------------------------- *------ AIC ------ *---- made of Ag+In+Cd ---- REAL MAIC ; REAL Ag107AICf Ag109AICf ; REAL Cd106AICf Cd108AICf Cd110AICf Cd111AICf Cd112AICf Cd113AICf ; REAL Cd114AICf Cd116AICf ; *REAL InAICf ; REAL In115AICf In113AICf ; GREP: AIC :: GETVAL 'AvgMolM ' 1 >>MAIC<< ; GREP: AIC :: GETVAL 'Ag107AICfIso' 1 >>Ag107AICf<< ; GREP: AIC :: GETVAL 'Ag109AICfIso' 1 >>Ag109AICf<< ; GREP: AIC :: GETVAL 'Cd106AICfIso' 1 >>Cd106AICf<< ; GREP: AIC :: GETVAL 'Cd108AICfIso' 1 >>Cd108AICf<< ; GREP: AIC :: GETVAL 'Cd110AICfIso' 1 >>Cd110AICf<< ; GREP: AIC :: GETVAL 'Cd111AICfIso' 1 >>Cd111AICf<< ; GREP: AIC :: GETVAL 'Cd112AICfIso' 1 >>Cd112AICf<< ; GREP: AIC :: GETVAL 'Cd113AICfIso' 1 >>Cd113AICf<< ; GREP: AIC :: GETVAL 'Cd114AICfIso' 1 >>Cd114AICf<< ; GREP: AIC :: GETVAL 'Cd116AICfIso' 1 >>Cd116AICf<< ; *GREP: AIC :: GETVAL 'In115AICfIso' 1 >>InAICf<< ; GREP: AIC :: GETVAL 'In115AICfIso' 1 >>In115AICf<< ; GREP: AIC :: GETVAL 'In113AICfIso' 1 >>In113AICf<< ; REAL AICconc := dAICr avo_conc * MAIC / ; REAL Ag107AIC := Ag107AICf AICconc * ; REAL Ag109AIC := Ag109AICf AICconc * ; REAL Cd106AIC := Cd106AICf AICconc * ; REAL Cd108AIC := Cd108AICf AICconc * ; REAL Cd110AIC := Cd110AICf AICconc * ; REAL Cd111AIC := Cd111AICf AICconc * ; REAL Cd112AIC := Cd112AICf AICconc * ; REAL Cd113AIC := Cd113AICf AICconc * ; REAL Cd114AIC := Cd114AICf AICconc * ; REAL Cd116AIC := Cd116AICf AICconc * ; *REAL InAIC := InAICf AICconc * ; REAL In115AIC := In115AICf AICconc * ; REAL In113AIC := In113AICf AICconc * ; *ECHO "Ag107AIC" Ag107AIC ; *ECHO "Ag109AIC" Ag109AIC ; *ECHO "Cd106AIC" Cd106AIC ; *ECHO "Cd108AIC" Cd108AIC ; *ECHO "Cd110AIC" Cd110AIC ; *ECHO "Cd111AIC" Cd111AIC ; *ECHO "Cd112AIC" Cd112AIC ; *ECHO "Cd113AIC" Cd113AIC ; *ECHO "Cd114AIC" Cd114AIC ; *ECHO "Cd116AIC" Cd116AIC ; *ECHO "In115AIC" In115AIC ; *ECHO "In113AIC" In113AIC ; *------ BNATB ------ *---- Natural Boron for structure ---- REAL MBNatB B10Bf B11Bf ; GREP: BNat :: GETVAL 'AvgMolM ' 1 >>MBNatB<< ; GREP: BNat :: GETVAL 'B10frIso' 1 >>B10Bf<< ; GREP: BNat :: GETVAL 'B11frIso' 1 >>B11Bf<< ; *------ B4C Control Rods ------ * Isotopic Composition : BNatB .8 * CNat .2 REAL MB4C ; REAL BNatBrodf := .8 ; REAL CNatrodf := .2 ; REAL B10Brodf := BNatBrodf B10Bf * ; REAL B11Brodf := BNatBrodf B11Bf * ; EVALUATE MB4C := MBNatB BNatBrodf * cnat CNatrodf * + ; REAL B4Cconc := dB4C avo_conc * MB4C / ; REAL B10Brod := B10Brodf B4Cconc * ; REAL B11Brod := B11Brodf B4Cconc * ; REAL Crod := CNatrodf B4Cconc * ; *ECHO "----B4C----" ; *ECHO "B10Brod" B10Brod ; *ECHO "B11Brod" B11Brod ; *ECHO "Crod " Crod ; *-------------------------------------------------------------------- * BP_T : Tihange Burnable Poison *-------------------------------------------------------------------- REAL MBP ; REAL B10BPf B11BPf ; REAL O16BPf SiNatBPf Al27BPf KNatBPf Na23BPf ; REAL Si28BPf Si29BPf Si30BPf K39BPf K40BPf K41BPf ; * GREP: BP :: GETVAL 'AvgMolM ' 1 >>MBP<< ; GREP: BP :: GETVAL 'B10BPfrIso ' 1 >>B10BPf<< ; GREP: BP :: GETVAL 'B11BPfrIso ' 1 >>B11BPf<< ; GREP: BP :: GETVAL 'O16BPfrIso ' 1 >>O16BPf<< ; GREP: BP :: GETVAL 'SiNatBPfrIso' 1 >>SiNatBPf<< ; GREP: BP :: GETVAL 'Si28BPfrIso ' 1 >>Si28BPf<< ; GREP: BP :: GETVAL 'Si29BPfrIso ' 1 >>Si29BPf<< ; GREP: BP :: GETVAL 'Si30BPfrIso ' 1 >>Si30BPf<< ; GREP: BP :: GETVAL 'Al27BPfrIso ' 1 >>Al27BPf<< ; GREP: BP :: GETVAL 'KNatBPfrIso ' 1 >>KNatBPf<< ; GREP: BP :: GETVAL 'K39BPfrIso ' 1 >>K39BPf<< ; GREP: BP :: GETVAL 'K40BPfrIso ' 1 >>K40BPf<< ; GREP: BP :: GETVAL 'K41BPfrIso ' 1 >>K41BPf<< ; GREP: BP :: GETVAL 'Na23BPfrIso ' 1 >>Na23BPf<< ; *-------------------------------------------------------------------- REAL BP_Tconc := dBP avo_conc * MBP / ; * REAL B10BP := B10BPf BP_Tconc * ; REAL B11BP := B11BPf BP_Tconc * ; REAL O16BP := O16BPf BP_Tconc * ; REAL SiNatBP := SiNatBPf BP_Tconc * ; REAL Si28BP := Si28BPf BP_Tconc * ; REAL Si29BP := Si29BPf BP_Tconc * ; REAL Si30BP := Si30BPf BP_Tconc * ; REAL Al27BP := Al27BPf BP_Tconc * ; REAL KNatBP := KNatBPf BP_Tconc * ; REAL K39BP := K39BPf BP_Tconc * ; REAL K40BP := K40BPf BP_Tconc * ; REAL K41BP := K41BPf BP_Tconc * ; REAL Na23BP := Na23BPf BP_Tconc * ; *---------------------------- *ECHO "----BP_T----" ; *ECHO "MBP " MBP ; *ECHO "BP_Tconc " BP_Tconc ; *ECHO "B10BP " B10BP ; *ECHO "B11BP " B11BP ; *ECHO "O16BP " O16BP ; *ECHO "SiNatBP " SiNatBP ; *ECHO "Si28BP " Si28BP ; *ECHO "Si29BP " Si29BP ; *ECHO "Si30BP " Si30BP ; *ECHO "Al27BP " Al27BP ; *ECHO "KNatBP " KNatBP ; *ECHO "K39BP " K39BP ; *ECHO "K40BP " K40BP ; *ECHO "K41BP " K41BP ; *ECHO "Na23BP " Na23BP ; *---------------------------- *------------------------------------------------------ * inrs order ('NAMISO = NISO <> inrs SHIB') * adding Mixes: ABSO, ABS_ACIER, GAINA+GAINAA * WARNING: Last (highest) Mix number MUST be equal to Nmix *------------------------------------------------------ *------------------------------------------------------ * Materiaux de base: 1 -> 9 * Moderateurs: 10 -> 19 * Structures: Tubes 20 * Barres 21 -> 29 * Gaines 30 -> 39 * Combustibles: UOX 100 -> 109 * MOX 110 -> 139 * UGd 140 -> Nmix *------------------------------------------------------ * Materiaux de base * 1 - SS304 * 2 - Inconel * 3 - Zircalloy-4 * 4 - Grilles precalculees * 5 - Eau boree * 6 - Isotopes de barre virtuelle * 7 - Eau boree tube instrumentation * Moderateurs * 10 - MODE [moderateur cellule centrale] * 11 - MODEL [moderateur cellule en peripherie (avec lame d'eau)] * 12 - MODEC [moderateur cellule en coin (double lame d'eau)] * 13 - MODETI [moderateur cellule tube (interieur du tube)+Isotopes de barres AIC a 0] * 14 - MODETE [moderateur cellule tube (exterieur du tube)] * 15 - MODETII [moderateur cellule tube (interieur du tube d'instrumentation)] * Structures: Tubes * 20 - TUBE * Structures: Barres * 21 - AIC [Barres AIC] * 22 - ACIER [Barres acier] * 23 - B4C [Barres B4C] * Structures: Gaines * 30 - GAINC [Gaine combustible UOX/MOX/UGd] * 31 - GAINA [Gaine Barres AIC] * 32 - GAINS [Gaine Barres ACIER] * 33 - GAINB [Gaine Barres B4C] * 34 - GAINE BP INT ACIER SS304 * 35 - GAINE BP EXT ACIER SS304 * 36 - GAINC M5 [Gaine combustible] * Poison consommable (Pyrex) * 60 - PYREX * Combustibles UOX * 100 - COMB [UOX 1ere couronne (centre a Rcomb*0.5^1/2) <=> 0 a r1] * 101 - COMB [UOX 2e couronne] * 102 - COMB [UOX 3e couronne] * 103 - COMB [UOX 4e couronne (r3 a Rcomb)] * Combustibles MOX * 110 - COMB1 [MOX 'faible' 1ere couronne] * 111 - COMB1 [MOX 'faible' 2e couronne] * 112 - COMB1 [MOX 'faible' 3e couronne] * 113 - COMB1 [MOX 'faible' 4e couronne] * 120 - COMB2 [MOX 'moyen' 1ere couronne] * 121 - COMB2 [MOX 'moyen' 2e couronne] * 122 - COMB2 [MOX 'moyen' 3e couronne] * 123 - COMB2 [MOX 'moyen' 4e couronne] * 130 - COMB3 [MOX 'fort' 1ere couronne] * 131 - COMB3 [MOX 'fort' 2e couronne] * 132 - COMB3 [MOX 'fort' 3e couronne] * 133 - COMB3 [MOX 'fort' 4e couronne] * Combustibles UGd * 140 - COMB [UOX+Gd 1ere couronne (centre a Rcomb*0.2^1/2)] * 141 - COMB [UOX+Gd 2e couronne] * 142 - COMB [UOX+Gd 3e couronne] * 143 - COMB [UOX+Gd 4e couronne] * 144 - COMB [UOX+Gd 5e couronne] * Nmix - COMB [UOX+Gd 6e couronne] --- Nmix = 145 *------------------------------------------------------ *------------------------------------------------------ * Depletion Chains *------------------------------------------------------ INTEGER Nmix := 145 ; ! Highest mix number must be equal to Nmix ! Procedure de definition de la chaine de decroissance STRING Autop := "SUBG" ; Library := DeplChain :: <> <> <> <> <> <> ; Library := LIB: Library :: EDIT 1 CTRA <> * Attention: si plusieurs appels a LIB: doivent etre faits, il faut * necessairement specifier pour chaque appel CTRA <>. Un bug * (non resolu actuellement) est present qui entraine des ecarts * lorsque cette condition n'est pas specifiee a chaque appel. MIXS LIB: <> FIL: <> MIX 1 <> NOEV (* SS304 *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Mn55 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> MIX 2 <> NOEV (* Inconel *) Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> B10 = <> <> B11 = <> <> Al27 = <> <> Mo95 = <> <> TiNat = <> <> ! Ti46 = <> <> ! Ti47 = <> <> ! Ti48 = <> <> ! Ti49 = <> <> ! Ti50 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> Co59 = <> <> Mn55 = <> <> Nb93 = <> <> MIX 3 <> NOEV (* Zr4 *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> O16 = <> <> ZrNat = <> <> ! Zr90 = <> <> ! Zr91 = <> <> ! Zr92 = <> <> ! Zr93 = <> <> ! Zr94 = <> <> ! Zr95 = <> <> ! Zr96 = <> <> ! Sn112 = <> <> ! Sn114 = <> <> Sn115 = <> <> Sn116 = <> <> Sn117 = <> <> Sn118 = <> <> Sn119 = <> <> Sn120 = <> <> Sn122 = <> <> Sn124 = <> <> MIX 4 <> NOEV (* GRILLE PRECALCULEE *) Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> B10 = <> <> B11 = <> <> Al27 = <> <> Mo95 = <> <> TiNat = <> <> ! Ti46 = <> <> ! Ti47 = <> <> ! Ti48 = <> <> ! Ti49 = <> <> ! Ti50 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> Co59 = <> <> Mn55 = <> <> Nb93 = <> <> ! Sn112 = <> <> ! Sn114 = <> <> Sn115 = <> <> Sn116 = <> <> Sn117 = <> <> Sn118 = <> <> Sn119 = <> <> Sn120 = <> <> Sn122 = <> <> Sn124 = <> <> O16 = <> <> ZrNat = <> <> ! Zr90 = <> <> ! Zr91 = <> <> ! Zr92 = <> <> ! Zr93 = <> <> ! Zr94 = <> <> ! Zr95 = <> <> ! Zr96 = <> <> MIX 5 <> NOEV ! (* Eau normale *) ! H2O = <> <> H1_H2O = <> <> O16 = <> <> B10 = <> <> B11 = <> <> MIX 6 <> NOEV *--- Isotopes de Barres pour Modeti a 0 --------------- * Utilité à réévaluer sur un calcul de Saphyb ! Ag107 = <> 0. Ag109 = <> 0. ! Cd106 = <> 0. ! Cd108 = <> 0. Cd110 = <> 0. Cd111 = <> 0. Cd112 = <> 0. Cd113 = <> 0. Cd114 = <> 0. Cd116 = <> 0. In115 = <> 0. ! In113 = <> 0. B10 = <> 0. B11 = <> 0. CNat = <> 0. MIX 7 <> NOEV ! (* Eau normale pour Tube Instr *) ! H2O = <> <> H1_H2O = <> <> O16 = <> <> B10 = <> <> B11 = <> <> U5DET = <> <> *---- Fin Isotopes de Barres pour Modeti a 0 ---------- MIX 10 (* MODE *) COMB 1 <> 2 <> 3 <> 4 <> 5 <> MIX 11 (* MODEL *) COMB 1 <> 2 <> 3 <> 4 <> 5 <> MIX 12 (* MODEC *) COMB 1 <> 2 <> 3 <> 4 <> 5 <> MIX 13 (* MODETI *) COMB 5 1.0 6 1.0E-20 MIX 14 (* MODETE *) COMB 1 <> 2 <> 3 <> 4 <> 5 <> MIX 15 (* MODETI POUR TUBE INSTR*) COMB 6 1.0E-20 7 1.0 *--- TUBES -------------------------------------------- MIX 20 <> NOEV (* TUBE *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> O16 = <> <> ZrNat = <> <> <> SHIB <> ! Zr90 = <> <> <> SHIB <> ! Zr91 = <> <> <> SHIB <> ! Zr92 = <> <> <> SHIB <> ! Zr93 = <> <> ! Zr94 = <> <> <> SHIB <> ! Zr95 = <> <> ! Zr96 = <> <> <> SHIB <> ! Sn112 = <> <> ! Sn114 = <> <> Sn115 = <> <> Sn116 = <> <> Sn117 = <> <> Sn118 = <> <> Sn119 = <> <> Sn120 = <> <> Sn122 = <> <> Sn124 = <> <> *--- BARRES ------------------------------------------- MIX 21 <> NOEV (* ABSO AIC *) ! Ag107 = <> <> <> SHIB <> Ag109 = <> <> <> SHIB <> ! Cd106 = <> <> ! Cd108 = <> <> Cd110 = <> <> Cd111 = <> <> Cd112 = <> <> Cd113 = <> <> Cd114 = <> <> Cd116 = <> <> In115 = <> <> <> SHIB <> ! In113 = <> <> <> SHIB <> MIX 22 <> NOEV (* ABS_ACIER *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Mn55 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> MIX 23 <> NOEV (* ABSO B4C *) B10B = <> <> B11B = <> <> CNat = <> <> *--- GAINES ------------------------------------------- MIX 30 <> NOEV (* GAINE COMB *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> O16 = <> <> ZrNat = <> <> <> SHIB <> ! Zr90 = <> <> <> SHIB <> ! Zr91 = <> <> <> SHIB <> ! Zr92 = <> <> <> SHIB <> ! Zr93 = <> <> ! Zr94 = <> <> <> SHIB <> ! Zr95 = <> <> ! Zr96 = <> <> <> SHIB <> ! Sn112 = <> <> ! Sn114 = <> <> Sn115 = <> <> Sn116 = <> <> Sn117 = <> <> Sn118 = <> <> Sn119 = <> <> Sn120 = <> <> Sn122 = <> <> Sn124 = <> <> MIX 31 <> NOEV (* GAINA & GAINAA for ABSO AIC *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Mn55 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> MIX 32 <> NOEV (* GAINA & GAINAA for ABS_ACIER *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Mn55 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> MIX 33 <> NOEV (* GAINA & GAINAA for ABSO B4C *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Mn55 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> MIX 34 <> NOEV (* GAINE BP INT ACIER SS304 *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Mn55 = <> <> MIX 35 <> NOEV (* GAINE BP EXT ACIER SS304 *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> Ni58 = <> <> Ni60 = <> <> Ni61 = <> <> Ni62 = <> <> Ni64 = <> <> Cr50 = <> <> Cr52 = <> <> Cr53 = <> <> Cr54 = <> <> Mn55 = <> <> MIX 36 <> NOEV (* GAINE M5 *) Fe54 = <> <> Fe56 = <> <> Fe57 = <> <> Fe58 = <> <> O16 = <> <> ZrNat = <> <> <> SHIB <> ! Zr90 = <> <> <> SHIB <> ! Zr91 = <> <> <> SHIB <> ! Zr92 = <> <> <> SHIB <> ! Zr93 = <> <> ! Zr94 = <> <> <> SHIB <> ! Zr95 = <> <> ! Zr96 = <> <> <> SHIB <> Nb93 = <> <> S32 = <> <> S33 = <> <> S34 = <> <> MIX 59 <> NOEV (* Air, origine CRISTAL *) N14 = <> 4.19850E-5 O16 = <> 1.12630E-5 MIX 60 <> (* BP Pyrex *) B10 = <> <> B11 = <> <> O16 = <> <> SiNat = <> <> ! Si28 = <> <> ! Si29 = <> <> ! Si30 = <> <> Al27 = <> <> KNat = <> <> ! K39 = <> <> ! K40 = <> <> ! K41 = <> <> Na23 = <> <> *--- COMB UOX ----------------------------------------- MIX 100 <> (* COMB *) U234 = <> <> U235 = <> <> <> SHIB <> U236 = <> <> U238 = <> <> <> SHIB <> PU238 = <> <> PU239 = <> <> <> SHIB <> PU240 = <> <> <> SHIB <> PU241 = <> <> <> SHIB <> PU242 = <> <> <> SHIB <> AM241 = <> <> O16 = <> <> MIX 101 COMB 100 1. MIX 102 COMB 100 1. MIX 103 COMB 100 1. *--- COMB MOX ----------------------------------------- MIX 110 <> (* COMB_P [LOW PU] *) U234 = <> <> U235 = <> <> <> SHIB <> U236 = <> <> U238 = <> <> <> SHIB <> PU238 = <> <> PU239 = <> <> <> SHIB <> PU240 = <> <> <> SHIB <> PU241 = <> <> <> SHIB <> PU242 = <> <> <> SHIB <> AM241 = <> <> O16 = <> <> MIX 111 COMB 110 1. MIX 112 COMB 110 1. MIX 113 COMB 110 1. MIX 120 <> (* COMB_I [AVG PU] *) U234 = <> <> U235 = <> <> <> SHIB <> U236 = <> <> U238 = <> <> <> SHIB <> PU238 = <> <> PU239 = <> <> <> SHIB <> PU240 = <> <> <> SHIB <> PU241 = <> <> <> SHIB <> PU242 = <> <> <> SHIB <> AM241 = <> <> O16 = <> <> MIX 121 COMB 120 1. MIX 122 COMB 120 1. MIX 123 COMB 120 1. MIX 130 <> (* COMB_C [HIGH PU] *) U234 = <> <> U235 = <> <> <> SHIB <> U236 = <> <> U238 = <> <> <> SHIB <> PU238 = <> <> PU239 = <> <> <> SHIB <> PU240 = <> <> <> SHIB <> PU241 = <> <> <> SHIB <> PU242 = <> <> <> SHIB <> AM241 = <> <> O16 = <> <> MIX 131 COMB 130 1. MIX 132 COMB 130 1. MIX 133 COMB 130 1. *--- COMB UGd ----------------------------------------- MIX 140 <> (* COMB with Gd *) Gd154 = <> <> <> SHIB <> Gd155 = <> <> <> SHIB <> Gd156 = <> <> <> SHIB <> Gd157 = <> <> <> SHIB <> Gd158 = <> <> <> SHIB <> Gd160 = <> <> <> SHIB <> U234 = <> <> U235 = <> <> <> SHIB <> U236 = <> <> U238 = <> <> <> SHIB <> PU238 = <> 0. PU239 = <> <> <> SHIB <> PU240 = <> <> <> SHIB <> PU241 = <> <> <> SHIB <> PU242 = <> <> <> SHIB <> AM241 = <> 0. O16 = <> <> MIX 141 COMB 140 1. MIX 142 COMB 140 1. MIX 143 COMB 140 1. MIX 144 COMB 140 1. ! Highest Mix number MUST be equal to NMIX MIX <> COMB 140 1. ; *------------------------------------------------------ * Moderator for every cell, all based on MODETI (internal guide tube): * MODETE (external guide tube) proportion * MODE (central fuel pins) proportion * MODEL (side fuel pins) proportion * MODEC (corner fuel pin) proportion *------------------------------------------------------ * 'MODETE' (MODErateur Tube Externe) *------------------------------------------------------ REAL H2O_mte := H2OTube H2O_mti * ; REAL H_mte := H2OTube H_mti * ; REAL O_mte := H2OTube O_mti * ; REAL B10_mte := H2OTube B10_mti * ; REAL B11_mte := H2OTube B11_mti * ; REAL B10B_mte := GrPreTube B10BGr * ; REAL B11B_mte := GrPreTube B11BGr * ; *ECHO "MODETE CONC" ; *ECHO "H2O_mte " H2O_mte ; *ECHO "H_mte " H_mte ; *ECHO "O_mte " O_mte ; *ECHO "B10_mte " B10_mte ; *ECHO "B11_mte " B11_mte ; *ECHO "B10B_mte " B10B_mte ; *ECHO "B11B_mte " B11B_mte ; *------------------------------------------------------ * 'MODE' (MODErateur cellule centrale) *------------------------------------------------------ REAL H2O_mod := H2OCntr H2O_mti * ; REAL H_mod := H2OCntr H_mti * ; REAL O_mod := H2OCntr O_mti * ; REAL B10_mod := H2OCntr B10_mti * ; REAL B11_mod := H2OCntr B11_mti * ; REAL B10B_mod := GrPreCntr B10BGr * ; REAL B11B_mod := GrPreCntr B11BGr * ; *ECHO "MODE CONC" ; *ECHO "H2O_mod " H2O_mod ; *ECHO "H_mod " H_mod ; *ECHO "O_mod " O_mod ; *ECHO "B10_mod " B10_mod ; *ECHO "B11_mod " B11_mod ; *ECHO "B10B_mod " B10B_mod ; *ECHO "B11B_mod " B11B_mod ; *------------------------------------------------------ * 'MODEL' (MODErateur cellule Laterale) *------------------------------------------------------ REAL H2O_mdl := H2OSide H2O_mti * ; REAL H_mdl := H2OSide H_mti * ; REAL O_mdl := H2OSide O_mti * ; REAL B10_mdl := H2OSide B10_mti * ; REAL B11_mdl := H2OSide B11_mti * ; REAL B10B_mdl := GrPreSide B10BGr * ; REAL B11B_mdl := GrPreSide B11BGr * ; *ECHO "MODEL CONC" ; *ECHO "H2O_mdl " H2O_mdl ; *ECHO "H_mdl " H_mdl ; *ECHO "O_mdl " O_mdl ; *ECHO "B10_mdl " B10_mdl ; *ECHO "B11_mdl " B11_mdl ; *ECHO "B10B_mdl " B10B_mdl ; *ECHO "B11B_mdl " B11B_mdl ; *------------------------------------------------------ * 'MODEC' (MODErateur cellule Coin) *------------------------------------------------------ REAL H2O_mdc := H2OCorn H2O_mti * ; REAL H_mdc := H2OCorn H_mti * ; REAL O_mdc := H2OCorn O_mti * ; REAL B10_mdc := H2OCorn B10_mti * ; REAL B11_mdc := H2OCorn B11_mti * ; REAL B10B_mdc := GrPreCorn B10BGr * ; REAL B11B_mdc := GrPreCorn B11BGr * ; *ECHO "MODEC CONC" ; *ECHO "H2O_mdc " H2O_mdc ; *ECHO "H_mdc " H_mdc ; *ECHO "O_mdc " O_mdc ; *ECHO "B10_mdc " B10_mdc ; *ECHO "B11_mdc " B11_mdc ; *ECHO "B10B_mdc " B10B_mdc ; *ECHO "B11B_mdc " B11B_mdc ; *------------------------------------------------------ * Somme des B10/B11 des grilles et de ceux de l'eau, * pour ConcMode (variation de CB lors du MakeSap). * Il serait mieux de supprimer cette 'astuce' et * d'isoler les calculs de H2O, B10 et B11 dans une * procedure appelee par CreaMix et MakeSap1BU. * NB : d'ici la, en realite, le calcul realise dans * MakeSap1BU n'est pas tout a fait exact (auparavant, * il ne l'etait pas non plus suite a un bug ...). * L'astuce ne marche que si cette variation lineaire * ne s'applique que sur le B10 et B11 de l'eau. En * pratique, le B10 et le B11 des grilles est trop peu * present pour avoir un impact. *------------------------------------------------------ EVALUATE B10_mod := B10_mod B10B_mod + ; EVALUATE B11_mod := B11_mod B11B_mod + ; EVALUATE B10_mdl := B10_mdl B10B_mdl + ; EVALUATE B11_mdl := B11_mdl B11B_mdl + ; EVALUATE B10_mdc := B10_mdc B10B_mdc + ; EVALUATE B11_mdc := B11_mdc B11B_mdc + ; EVALUATE B10_mte := B10_mte B10B_mte + ; EVALUATE B11_mte := B11_mte B11B_mte + ; *------------------------------------------------------ * Moderator Concentrations - for CBOR DMOD variation *------------------------------------------------------ ConcMode := UTL: :: CREA 'MODETI' 4 = <> <> <> <> ; ConcMode := UTL: ConcMode :: CREA 'MODETE' 4 = <> <> <> <> ; ConcMode := UTL: ConcMode :: CREA 'MODE ' 4 = <> <> <> <> ; ConcMode := UTL: ConcMode :: CREA 'MODEL ' 4 = <> <> <> <> ; ConcMode := UTL: ConcMode :: CREA 'MODEC ' 4 = <> <> <> <> ; *------------------------------------------------------ END: ; QUIT "LIST" .