diff options
Diffstat (limited to 'Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m')
| -rw-r--r-- | Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m | 1548 |
1 files changed, 1548 insertions, 0 deletions
diff --git a/Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m b/Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m new file mode 100644 index 0000000..607e3a6 --- /dev/null +++ b/Dragon/data/ErmBeavrsPwrRefl_proc/MetalMix.c2m @@ -0,0 +1,1548 @@ +*DECK MetalMix.c2m +*====================================================== +* Name: MetalMix.c2m +* Type: DRAGON procedure +* Use: Create isotopic mixtures of metals. +* Author: A. Bruneau 10.2013 +* [CS-SI | IRSN/PSN-EXP/SNC/LNR] +* Mod: A. Bruneau 12.2013 +* [CS-SI | IRSN/PSN-EXP/SNC/LNR] +* Added AIC with AgNat, CdNat & In for AIC rod description. +* Mod: A. Bruneau 02.2014 +* [CS-SI | IRSN/PSN-EXP/SNC/LNR] +* Added GdNat with Gd for Gd2O3 +* Mod: G. Tixier 09.2016 +* [CS-SI | IRSN/PSN-EXP/SNC/LNR] +* Ajouts S32, S33 et S34 pour materiau M5 +*------------------------------------------------------ +* Procedure called as: +*------------------------------------------------------ +* Input Required: +* > Molar Masses +* Output Provided: +* < Molar Masses +* < Metals Concentration +*====================================================== +PARAMETER FeNat CrNat SnNat NiNat BNat + ZrNat TiNat SiNat KNat + AgNat CdNat SNat InNat Zirc4 Incon SS304 + M5 AIC GdNat Gd2O3 BP MolarMasses :: + ::: LINKED_LIST + FeNat CrNat SnNat NiNat BNat + ZrNat TiNat SiNat KNat + AgNat CdNat SNat InNat Zirc4 Incon SS304 + M5 AIC GdNat Gd2O3 BP MolarMasses ; ; +MODULE GREP: UTL: END: ABORT: ; +STRING LibType ; ! ENDFVII1 JEFF31 JEFF311 JEFF312... +STRING Gestion ; +REAL avogadro := 6.022094E+23 ; +REAL avo_conc := avogadro 1.E-24 * ; +REAL ReNorm ; ! Renormalization Factor +REAL SumFr ; ! Fraction Sum (needed when &IMPU) +:: >>LibType<< >>Gestion<< ; + +*------------------------------------------------------ +REAL am241 ; +REAL al27 ; +REAL b10 ; +REAL b10b ; +REAL b11 ; +REAL b11b ; +REAL co59 ; +REAL cr50 ; +REAL cr52 ; +REAL cr53 ; +REAL cr54 ; +REAL fe54 ; +REAL fe56 ; +REAL fe57 ; +REAL fe58 ; +REAL h2o h1 ; +REAL mn55 ; +REAL mo95 ; +REAL nb93 ; +REAL ni58 ; +REAL ni60 ; +REAL ni61 ; +REAL ni62 ; +REAL ni64 ; +REAL o16 ; +REAL pu238 ; +REAL pu239 ; +REAL pu240 ; +REAL pu241 ; +REAL pu242 ; +REAL siNat si28 si29 si30 ; +REAL sn112 ; +REAL sn114 ; +REAL sn115 ; +REAL sn116 ; +REAL sn117 ; +REAL sn118 ; +REAL sn119 ; +REAL sn120 ; +REAL sn122 ; +REAL sn124 ; +REAL tiNat ti46 ti47 ti48 ti49 ti50 ; +REAL u234 ; +REAL u235 ; +REAL u236 ; +REAL u238 ; +REAL zrNat zr90 zr91 zr92 zr93 zr94 zr95 zr96 ; +REAL cnat ; +*---- added for AIC +REAL ag107 ; +REAL ag109 ; +REAL cd106 ; +REAL cd108 ; +REAL cd110 ; +REAL cd111 ; +REAL cd112 ; +REAL cd113 ; +REAL cd114 ; +REAL cd116 ; +REAL in115 in113 ; +*---- added for Gd2O3 +REAL gd154 ; +REAL gd155 ; +REAL gd156 ; +REAL gd157 ; +REAL gd158 ; +REAL gd160 ; +*---- added for Pyrex +REAL kNat k39 k40 k41 ; +REAL na23 ; +*---- added for M5 +REAL s32 ; +REAL s33 ; +REAL s34 ; + +*------------------------------------------------------ +*------------------------------------------------------ +* Recuperation des fractions massiques moyennes +*------------------------------------------------------ + +* Pour les isotopes non conflictuels + +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 'CNat ' 1 >>cnat<< ; +GREP: MolarMasses :: GETVAL 'O16 ' 1 >>o16<< ; +GREP: MolarMasses :: GETVAL 'Al27 ' 1 >>al27<< ; +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 'Mn55 ' 1 >>mn55<< ; +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 'Co59 ' 1 >>co59<< ; +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 'Nb93 ' 1 >>nb93<< ; +GREP: MolarMasses :: GETVAL 'Mo95 ' 1 >>mo95<< ; +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 '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 '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 'U234 ' 1 >>u234<< ; +GREP: MolarMasses :: GETVAL 'U235 ' 1 >>u235<< ; +GREP: MolarMasses :: GETVAL 'U236 ' 1 >>u236<< ; +GREP: MolarMasses :: GETVAL 'U238 ' 1 >>u238<< ; +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 'Am241' 1 >>am241<< ; +GREP: MolarMasses :: GETVAL 'Na23 ' 1 >>na23<< ; + +* Pour les isotopes conflictuels (absent 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<< ; + GREP: MolarMasses :: GETVAL 'InNat' 1 >>in115<< ; + ! In115 non defini dans ENDFVI8 + 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 >>si29<< ; + GREP: MolarMasses :: GETVAL 'Si30 ' 1 >>si30<< ; + 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 ; + +*------------------------------------------------------ +* ZrNat : Natural Zirconium +* Isotopic Composition : Zr90 .5145 +* Zr91 .1123 +* Zr92 .1716 +* Zr93 . +* Zr94 .1738 +* Zr95 . +* Zr96 .0281 +*------------------------------------------------------ +REAL Zr90f := .0 ; +REAL Zr91f := .0 ; +REAL Zr92f := .0 ; +REAL Zr93f := .0 ; +REAL Zr94f := .0 ; +REAL Zr95f := .0 ; +REAL Zr96f := .0 ; +REAL MZrNat ; ! Average Weight Ratio + +IF LibType 'JEFF31' = LibType 'JEFF311' = + + LibType 'JEFF312' = + + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN + + EVALUATE Zr90f := .5145 ; + EVALUATE Zr91f := .1122 ; + EVALUATE Zr92f := .1715 ; + EVALUATE Zr93f := 0. ; + EVALUATE Zr94f := .1738 ; + EVALUATE Zr95f := 0. ; + EVALUATE Zr96f := .0280 ; + + *---- ZrNat Average Weight Ratio ---- + EVALUATE MZrNat := zr90 Zr90f * zr91 Zr91f * + + zr92 Zr92f * + zr93 Zr93f * + + zr94 Zr94f * + zr95 Zr95f * + + zr96 Zr96f * + ; + + ECHO "Masse molaire du zirconium " MZrNat ; + +ELSE ! JEF22 + + EVALUATE MZrNat := zrNat ; + +ENDIF ; + +*---- ZrNat Table Creation ---- +ZrNat := UTL: :: CREA 'AvgMolM ' 1 = <<MZrNat>> ; +ZrNat := UTL: ZrNat :: CREA 'Zr90frIso' 1 = <<Zr90f>> ; +ZrNat := UTL: ZrNat :: CREA 'Zr91frIso' 1 = <<Zr91f>> ; +ZrNat := UTL: ZrNat :: CREA 'Zr92frIso' 1 = <<Zr92f>> ; +ZrNat := UTL: ZrNat :: CREA 'Zr93frIso' 1 = <<Zr93f>> ; +ZrNat := UTL: ZrNat :: CREA 'Zr94frIso' 1 = <<Zr94f>> ; +ZrNat := UTL: ZrNat :: CREA 'Zr95frIso' 1 = <<Zr95f>> ; +ZrNat := UTL: ZrNat :: CREA 'Zr96frIso' 1 = <<Zr96f>> ; + +*------------------------------------------------------ +* FeNat : Natural Iron +* Isotopic Composition : Fe54 .0580 +* Fe56 .9172 +* Fe57 .0220 +* Fe58 .0028 +*------------------------------------------------------ +REAL Fe54f := .0580 ; +REAL Fe56f := .9172 ; +REAL Fe57f := .0220 ; +REAL Fe58f := .0028 ; + +*---- FeNat Average Weight Ratio ---- +REAL MFeNat ; ! Average Weight Ratio +EVALUATE MFeNat := fe54 Fe54f * fe56 Fe56f * + ; +EVALUATE MFeNat := MFeNat fe57 Fe57f * + fe58 Fe58f * + ; + +*---- FeNat Table Creation ---- +FeNat := UTL: :: CREA 'AvgMolM ' 1 = <<MFeNat>> ; +FeNat := UTL: FeNat :: CREA 'Fe54frIso' 1 = <<Fe54f>> ; +FeNat := UTL: FeNat :: CREA 'Fe56frIso' 1 = <<Fe56f>> ; +FeNat := UTL: FeNat :: CREA 'Fe57frIso' 1 = <<Fe57f>> ; +FeNat := UTL: FeNat :: CREA 'Fe58frIso' 1 = <<Fe58f>> ; + +*------------------------------------------------------ +* TiNat : Natural Titanium +* Isotopic Composition : Ti46 .0825 +* Ti47 .0744 +* Ti48 .7372 +* Ti49 .0541 +* Ti50 .0518 +*------------------------------------------------------ +REAL Ti46f := .0 ; +REAL Ti47f := .0 ; +REAL Ti48f := .0 ; +REAL Ti49f := .0 ; +REAL Ti50f := .0 ; +REAL MTiNat ; ! Average Weight Ratio + +IF LibType 'JEFF31' = LibType 'JEFF311' = + + LibType 'JEFF312' = + + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN + + EVALUATE Ti46f := .0825 ; + EVALUATE Ti47f := .0744 ; + EVALUATE Ti48f := .7372 ; + EVALUATE Ti49f := .0541 ; + EVALUATE Ti50f := .0518 ; + + *---- TiNat Average Weight Ratio ---- + EVALUATE MTiNat := ti46 Ti46f * ti47 Ti47f * + + ti48 Ti48f * + + ti49 Ti49f * + ti50 Ti50f * + ; + + ECHO "Masse molaire du titane " MTiNat ; + +ELSE ! JEF22 ENDFVI8 + + EVALUATE MTiNat := tiNat ; + +ENDIF ; + +*---- TiNat Table Creation ---- +TiNat := UTL: :: CREA 'AvgMolM ' 1 = <<MTiNat>> ; +TiNat := UTL: TiNat :: CREA 'Ti46frIso' 1 = <<Ti46f>> ; +TiNat := UTL: TiNat :: CREA 'Ti47frIso' 1 = <<Ti47f>> ; +TiNat := UTL: TiNat :: CREA 'Ti48frIso' 1 = <<Ti48f>> ; +TiNat := UTL: TiNat :: CREA 'Ti49frIso' 1 = <<Ti49f>> ; +TiNat := UTL: TiNat :: CREA 'Ti50frIso' 1 = <<Ti50f>> ; + +*------------------------------------------------------ +* SiNat : Natural Silicium +* Isotopic Composition : Si28 .922298 +* Si29 .046833 +* Si30 .030873 +*------------------------------------------------------ +REAL Si28f := .0 ; +REAL Si29f := .0 ; +REAL Si30f := .0 ; +REAL MSiNat ; ! Average Weight Ratio + + +IF LibType 'JEFF31' = LibType 'JEFF311' = + + LibType 'JEFF312' = + + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN + + EVALUATE Si28f := .9222 ; + EVALUATE Si29f := .0468 ; + EVALUATE Si30f := .0309 ; + + *---- SiNat Average Weight Ratio ---- + EVALUATE MSiNat := si28 Si28f * si29 Si29f * + + si30 Si30f * + ; + + ECHO "Masse molaire du silicium " MSiNat ; + +ELSE ! JEF22 ENDFVI8 + + EVALUATE MSiNat := siNat ; + +ENDIF ; + +*---- SiNat Table Creation ---- +SiNat := UTL: :: CREA 'AvgMolM ' 1 = <<MSiNat>> ; +SiNat := UTL: SiNat :: CREA 'Si28frIso' 1 = <<Si28f>> ; +SiNat := UTL: SiNat :: CREA 'Si29frIso' 1 = <<Si29f>> ; +SiNat := UTL: SiNat :: CREA 'Si30frIso' 1 = <<Si30f>> ; + +*------------------------------------------------------ +* CrNat : Natural Chromium +* Isotopic Composition : Cr50 .0435 +* Cr52 .8379 +* Cr53 .0950 +* Cr54 .0236 +*------------------------------------------------------ +REAL Cr50f := .0435 ; +REAL Cr52f := .8379 ; +REAL Cr53f := .0950 ; +REAL Cr54f := .0236 ; + +*---- CrNat Average Weight Ratio ---- +REAL MCrNat ; ! Average Weight Ratio +EVALUATE MCrNat := cr50 Cr50f * cr52 Cr52f * + ; +EVALUATE MCrNat := MCrNat cr53 Cr53f * + cr54 Cr54f * + ; + +*---- CrNat Table Creation ---- +CrNat := UTL: :: CREA 'AvgMolM ' 1 = <<MCrNat>> ; +CrNat := UTL: CrNat :: CREA 'Cr50frIso' 1 = <<Cr50f>> ; +CrNat := UTL: CrNat :: CREA 'Cr52frIso' 1 = <<Cr52f>> ; +CrNat := UTL: CrNat :: CREA 'Cr53frIso' 1 = <<Cr53f>> ; +CrNat := UTL: CrNat :: CREA 'Cr54frIso' 1 = <<Cr54f>> ; + +*------------------------------------------------------ +* SnNat : Natural Tin +* Isotopic Composition : Sn112 .010 Sauf JEF22 et ENDFVI8 +* Sn114 .007 Sauf JEF22 et ENDFVI8 +* Sn115 .004 +* Sn116 .147 / Sn116 .164 (JEF22, ENDFVI8) +* Sn117 .077 +* Sn118 .243 +* Sn119 .086 +* Sn120 .324 +* Sn122 .046 +* Sn124 .056 +* Pour Sn112 et Sn114 avec JEF22 ou ENDFVI8, on reporte les +* fractions atomiques sur Sn116 (premier isotope pair defini +* dans ces bibliotheques) +*------------------------------------------------------ +REAL Sn112f := .0 ; +REAL Sn114f := .0 ; +REAL Sn115f := .0 ; +REAL Sn116f := .0 ; +REAL Sn117f := .0 ; +REAL Sn118f := .0 ; +REAL Sn119f := .0 ; +REAL Sn120f := .0 ; +REAL Sn122f := .0 ; +REAL Sn124f := .0 ; +REAL MSnNat ; ! Average Weight Ratio + +IF LibType 'JEFF31' = LibType 'JEFF311' = + + LibType 'JEFF312' = + + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN + + EVALUATE Sn112f := .010 ; + EVALUATE Sn114f := .007 ; + EVALUATE Sn115f := .004 ; + EVALUATE Sn116f := .147 ; + EVALUATE Sn117f := .077 ; + EVALUATE Sn118f := .243 ; + EVALUATE Sn119f := .086 ; + EVALUATE Sn120f := .324 ; + EVALUATE Sn122f := .046 ; + EVALUATE Sn124f := .056 ; + + *---- SnNat Average Weight Ratio ---- + EVALUATE MSnNat := sn112 Sn112f * sn114 Sn114f * + ; + EVALUATE MSnNat := MSnNat sn115 Sn115f * + sn116 Sn116f * + ; + EVALUATE MSnNat := MSnNat sn117 Sn117f * + sn118 Sn118f * + ; + EVALUATE MSnNat := MSnNat sn119 Sn119f * + sn120 Sn120f * + ; + EVALUATE MSnNat := MSnNat sn122 Sn122f * + sn124 Sn124f * + ; + +ELSE ! JEF22, ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93 + + EVALUATE Sn115f := .004 ; + EVALUATE Sn116f := .164 ; + EVALUATE Sn117f := .077 ; + EVALUATE Sn118f := .243 ; + EVALUATE Sn119f := .086 ; + EVALUATE Sn120f := .324 ; + EVALUATE Sn122f := .046 ; + EVALUATE Sn124f := .056 ; + + *---- SnNat Average Weight Ratio ---- + EVALUATE MSnNat := sn115 Sn115f * sn116 Sn116f * + ; + EVALUATE MSnNat := MSnNat sn117 Sn117f * + sn118 Sn118f * + ; + EVALUATE MSnNat := MSnNat sn119 Sn119f * + sn120 Sn120f * + ; + EVALUATE MSnNat := MSnNat sn122 Sn122f * + sn124 Sn124f * + ; + +ENDIF ; + +*---- SnNat Table Creation ---- +SnNat := UTL: :: CREA 'AvgMolM ' 1 = <<MSnNat>> ; +SnNat := UTL: SnNat :: CREA 'Sn112frIso' 1 = <<Sn112f>> ; +SnNat := UTL: SnNat :: CREA 'Sn114frIso' 1 = <<Sn114f>> ; +SnNat := UTL: SnNat :: CREA 'Sn115frIso' 1 = <<Sn115f>> ; +SnNat := UTL: SnNat :: CREA 'Sn116frIso' 1 = <<Sn116f>> ; +SnNat := UTL: SnNat :: CREA 'Sn117frIso' 1 = <<Sn117f>> ; +SnNat := UTL: SnNat :: CREA 'Sn118frIso' 1 = <<Sn118f>> ; +SnNat := UTL: SnNat :: CREA 'Sn119frIso' 1 = <<Sn119f>> ; +SnNat := UTL: SnNat :: CREA 'Sn120frIso' 1 = <<Sn120f>> ; +SnNat := UTL: SnNat :: CREA 'Sn122frIso' 1 = <<Sn122f>> ; +SnNat := UTL: SnNat :: CREA 'Sn124frIso' 1 = <<Sn124f>> ; + +*------------------------------------------------------ +* NiNat : Natural Nickel +* Isotopic Composition : Ni58 .6827 +* Ni60 .2610 +* Ni61 .0113 +* Ni62 .0359 +* Ni64 .0091 +*------------------------------------------------------ +REAL Ni58f := .6827 ; +REAL Ni60f := .2610 ; +REAL Ni61f := .0113 ; +REAL Ni62f := .0359 ; +REAL Ni64f := .0091 ; + +*---- NiNat Average Weight Ratio ---- +REAL MNiNat ; ! Average Weight Ratio +EVALUATE MNiNat := ni58 Ni58f * ni60 Ni60f * + + ni61 Ni61f * + ni62 Ni62f * + + ni64 Ni64f * + ; + +*---- NiNat Table Creation ---- +NiNat := UTL: :: CREA 'AvgMolM ' 1 = <<MNiNat>> ; +NiNat := UTL: NiNat :: CREA 'Ni58frIso' 1 = <<Ni58f>> ; +NiNat := UTL: NiNat :: CREA 'Ni60frIso' 1 = <<Ni60f>> ; +NiNat := UTL: NiNat :: CREA 'Ni61frIso' 1 = <<Ni61f>> ; +NiNat := UTL: NiNat :: CREA 'Ni62frIso' 1 = <<Ni62f>> ; +NiNat := UTL: NiNat :: CREA 'Ni64frIso' 1 = <<Ni64f>> ; + +*------------------------------------------------------ +* BNat : Natural Boron +* Isotopic Composition : B10 .2 +* B11 .8 +*------------------------------------------------------ +REAL B10f := .2 ; +REAL B11f := .8 ; +IF +Gestion "BEAV" = THEN +EVALUATE B10f := 0.199 ; +EVALUATE B11f := 0.801 ; +ENDIF ; + +*---- BNat Average Weight Ratio ---- +REAL MBNat ; ! Average Weight Ratio +EVALUATE MBNat := b10 B10f * b11 B11f * + ; +IF Gestion "BEAV" = THEN +EVALUATE MBNat := 10.8135 ; +ENDIF ; + +*---- BNat Table Creation ---- +BNat := UTL: :: CREA 'AvgMolM ' 1 = <<MBNat>> ; +BNat := UTL: BNat :: CREA 'B10frIso' 1 = <<B10f>> ; +BNat := UTL: BNat :: CREA 'B11frIso' 1 = <<B11f>> ; + +*------------------------------------------------------ +* AgNat : Natural Silver +* Isotopic Composition : Ag107 .5183 +* Ag109 .4817 +*------------------------------------------------------ +REAL Ag107f := .5183 ; +REAL Ag109f := .4817 ; + +*---- AgNat Average Weight Ratio ---- +REAL MAgNat ; +EVALUATE MAgNat := ag107 Ag107f * ag109 Ag109f * + ; + +*---- AgNat Table Creation ---- +AgNat := UTL: :: CREA 'AvgMolM ' 1 = <<MAgNat>> ; +AgNat := UTL: AgNat :: CREA 'Ag107frIso ' 1 = <<Ag107f>> ; +AgNat := UTL: AgNat :: CREA 'Ag109frIso ' 1 = <<Ag109f>> ; + +*------------------------------------------------------ +* CdNat : Natural Cadmium +* Isotopic Composition : Cd106 .0125 * +* Cd108 .0089 ** +* Cd110 .1251 +* Cd111 .1281 +* Cd112 .2413 ** +* Cd113 .1222 +* Cd114 .2872 +* Cd116 .0747 +* +* * : Cd106 is not available on JEF22 and ENDFVI8. It is replaced +* by Cd112. +* For JEF22 and ENDFVI8 : +* ==> Cd106 .0000 +* ==> Cd108 .0000 +* ==> Cd112 .2627 +* ** : Cd108 is not available on JEF22 and ENDFVI8. +* It is replaced by Cd112 which has relatively close enough +* neutronic properties +* +* ==> Cd108 .0000 +* ==> Cd112 .2502 +*------------------------------------------------------ +REAL Cd106f := .0 ; +REAL Cd108f := .0 ; +REAL Cd110f := .0 ; +REAL Cd111f := .0 ; +REAL Cd112f := .0 ; +REAL Cd113f := .0 ; +REAL Cd114f := .0 ; +REAL Cd116f := .0 ; +REAL MCdNat ; + + +IF LibType 'JEFF31' = LibType 'JEFF311' = + + LibType 'JEFF312' = + + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN + + EVALUATE Cd106f := .0125 ; + EVALUATE Cd108f := .0089 ; + EVALUATE Cd110f := .1251 ; + EVALUATE Cd111f := .1281 ; + EVALUATE Cd112f := .2413 ; + EVALUATE Cd113f := .1222 ; + EVALUATE Cd114f := .2872 ; + EVALUATE Cd116f := .0747 ; + +*---- CdNat Average Weight Ratio ---- + EVALUATE MCdNat := cd106 Cd106f * cd108 Cd108f * + ; + EVALUATE MCdNat := MCdNat cd110 Cd110f * + cd111 Cd111f * + ; + EVALUATE MCdNat := MCdNat cd112 Cd112f * + cd113 Cd113f * + ; + EVALUATE MCdNat := MCdNat cd114 Cd114f * + cd116 Cd116f * + ; + +ELSE ! JEF22, ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93 + + EVALUATE Cd110f := .1251 ; + EVALUATE Cd111f := .1281 ; + EVALUATE Cd112f := .2627 ; + EVALUATE Cd113f := .1222 ; + EVALUATE Cd114f := .2872 ; + EVALUATE Cd116f := .0747 ; + +*---- CdNat Average Weight Ratio ---- + EVALUATE MCdNat := cd110 Cd110f * ; + EVALUATE MCdNat := MCdNat cd111 Cd111f * + cd112 Cd112f * + ; + EVALUATE MCdNat := MCdNat cd113 Cd113f * + cd114 Cd114f * + ; + EVALUATE MCdNat := MCdNat cd116 Cd116f * + ; + +ENDIF ; + +*---- CdNat Table Creation ---- +CdNat := UTL: :: CREA 'AvgMolM ' 1 = <<MCdNat>> ; +CdNat := UTL: CdNat :: CREA 'Cd106frIso ' 1 = <<Cd106f>> ; +CdNat := UTL: CdNat :: CREA 'Cd108frIso ' 1 = <<Cd108f>> ; +CdNat := UTL: CdNat :: CREA 'Cd110frIso ' 1 = <<Cd110f>> ; +CdNat := UTL: CdNat :: CREA 'Cd111frIso ' 1 = <<Cd111f>> ; +CdNat := UTL: CdNat :: CREA 'Cd112frIso ' 1 = <<Cd112f>> ; +CdNat := UTL: CdNat :: CREA 'Cd113frIso ' 1 = <<Cd113f>> ; +CdNat := UTL: CdNat :: CREA 'Cd114frIso ' 1 = <<Cd114f>> ; +CdNat := UTL: CdNat :: CREA 'Cd116frIso ' 1 = <<Cd116f>> ; + +*------------------------------------------------------ +* KNat : Natural Potassium +* Reference : www.wikipedia.org/wiki/Potassium +* Isotopic Composition : K39 .9326 +* K40 .0001167 +* K41 .0673 +*------------------------------------------------------ +REAL K39f := .0 ; +REAL K40f := .0 ; +REAL K41f := .0 ; +REAL MKNat ; ! Average Weight Ratio + + +IF LibType 'JEFF31' = LibType 'JEFF311' = + + LibType 'JEFF312' = + + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN + + EVALUATE K39f := .9326 ; + EVALUATE K40f := .0001167 ; + EVALUATE K41f := .0673 ; + +*---- KNat Average Weight Ratio ---- + EVALUATE MKNat := k39 K39f * k40 K40f * + k41 K41f * + ; + + ECHO "Masse molaire du potassium " MKNat ; + +ELSE ! JEF22, ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93 + + EVALUATE MKNat := kNat ; + +ENDIF ; + +*---- KNat Table Creation ---- +KNat := UTL: :: CREA 'AvgMolM ' 1 = <<MKNat>> ; +KNat := UTL: KNat :: CREA 'K39frIso' 1 = <<K39f>> ; +KNat := UTL: KNat :: CREA 'K40frIso' 1 = <<K40f>> ; +KNat := UTL: KNat :: CREA 'K41frIso' 1 = <<K41f>> ; + + +*------------------------------------------------------ +* SNat : Natural Sulfur +* Reference : www.wikipedia.org/wiki/Soufre +* Isotopic Composition : S32 0.9502 +* S33 0.0076 +* S34 0.0422 +* NB: pas de S33 et S34 dans ENDFBVI8 +* Pour S34 dans ENDFVI8, on reporte les +* fractions atomiques sur S32 (premier isotope pair defini +* dans ces bibliotheques) +*------------------------------------------------------ + +REAL S32f := .0 ; +REAL S33f := .0 ; +REAL S34f := .0 ; + +REAL MSNat ; ! Average Weight Ratio +ECHO "libtype:" LibType ; +IF LibType 'JEFF31' = LibType 'JEFF311' = + + LibType 'JEFF312' = + LibType 'JEFF22' = + + LibType 'ENDFVII0' = + LibType 'ENDFVII1' = + THEN + + EVALUATE S32f := .9502 ; + EVALUATE S33f := .0076 ; + EVALUATE S34f := .0422 ; + +*---- SNat Average Weight Ratio ---- + EVALUATE MSNat := S32f s32 * S33f s33 * + ; + EVALUATE MSNat := MSNat S34f s34 * + ; + +ELSE ! ENDFVI8, CEA99v9p5, CEA514T2, CLA99CEA93 + + EVALUATE S32f := .9924 ; + + *---- SnNat Average Weight Ratio ---- + EVALUATE MSNat := s32 S32f * ; + +ENDIF ; + +*---- SNat Table Creation ---- +SNat := UTL: :: CREA 'AvgMolM ' 1 = <<MSNat>> ; +SNat := UTL: SNat :: CREA 'S32frIso' 1 = <<S32f>> ; +SNat := UTL: SNat :: CREA 'S33frIso' 1 = <<S33f>> ; +SNat := UTL: SNat :: CREA 'S34frIso' 1 = <<S34f>> ; + +*------------------------------------------------------ +* InNat : Natural Indium +* Reference : www.webelements.com/indium/isotopes.html +* Isotopic Composition : In115 0.9571 +* In113 0.0429 +* NB : In113 non disponible. On considere que +* l'InNat correspond a 100% de In115 +*------------------------------------------------------ +REAL In113f := 0. ; +REAL In115f := 1. ; +REAL MInNat ; ! Average Weight Ratio + +*---- SNat Average Weight Ratio ---- +EVALUATE MInNat := In115f in115 * ; + +*---- InNat Table Creation ---- +InNat := UTL: :: CREA 'AvgMolM ' 1 = <<MInNat>> ; +InNat := UTL: InNat :: CREA 'In115frIso' 1 = <<In115f>> ; +InNat := UTL: InNat :: CREA 'In113frIso' 1 = <<In113f>> ; + +*------------------------------------------------------ +* Zirc4 : Zircalloy4 +* Massic Composition : FeNat .002 +* CrNat .001 +* 016 .001 +* ZrNat .981 +* SnNat .015 +*------------------------------------------------------ +REAL FeZ4f := .002 ; +REAL CrZ4f := .001 ; +REAL O16Z4f := .001 ; +REAL ZrZ4f := .981 ; +REAL SnZ4f := .015 ; + +REAL 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 ; + +*---- Mass Ratios -> Isotopic Ratios ---- +EVALUATE FeZ4f := FeZ4f avo_conc * MFeNat / ; +EVALUATE CrZ4f := CrZ4f avo_conc * MCrNat / ; +EVALUATE O16Z4f := O16Z4f avo_conc * o16 / ; +EVALUATE ZrZ4f := ZrZ4f avo_conc * MZrNat / ; +EVALUATE SnZ4f := SnZ4f avo_conc * MSnNat / ; + +EVALUATE ReNorm := FeZ4f CrZ4f + O16Z4f + ZrZ4f + SnZ4f + ; + +EVALUATE FeZ4f := FeZ4f ReNorm / ; +EVALUATE CrZ4f := CrZ4f ReNorm / ; +EVALUATE O16Z4f := O16Z4f ReNorm / ; +EVALUATE ZrZ4f := ZrZ4f ReNorm / ; +EVALUATE SnZ4f := SnZ4f ReNorm / ; + +*---- Zirc4 Average Weight Ratio --- +REAL MZirc4 ; ! Average Weight Ratio +EVALUATE MZirc4 := FeZ4f MFeNat * CrZ4f MCrNat * + + O16Z4f o16 * + ZrZ4f MZrNat * + + SnZ4f MSnNat * + ; + +ECHO "Masse molaire du zircaloy " MZirc4 ; + +*---- Isotopic Ratios for all Zr, Fe, Cr & Sn ---- +* Rq : si les isotopes n'existent pas ce n'est pas grave, +* leur valeur sera egale a 0. dans tous les calculs +* et leur appel sera commente dans l'appel a LIB de +* la procedure CreaMix.c2m. Leur presence n'a aucune +* influence sur les calculs du fait de leur valeur = 0. +*------------------------------------------------- +EVALUATE Zr90Z4f := ZrZ4f Zr90f * ; +EVALUATE Zr91Z4f := ZrZ4f Zr91f * ; +EVALUATE Zr92Z4f := ZrZ4f Zr92f * ; +EVALUATE Zr93Z4f := ZrZ4f Zr93f * ; +EVALUATE Zr94Z4f := ZrZ4f Zr94f * ; +EVALUATE Zr95Z4f := ZrZ4f Zr95f * ; +EVALUATE Zr96Z4f := ZrZ4f Zr96f * ; + +EVALUATE Fe54Z4f := FeZ4f Fe54f * ; +EVALUATE Fe56Z4f := FeZ4f Fe56f * ; +EVALUATE Fe57Z4f := FeZ4f Fe57f * ; +EVALUATE Fe58Z4f := FeZ4f Fe58f * ; + +EVALUATE Cr50Z4f := CrZ4f Cr50f * ; +EVALUATE Cr52Z4f := CrZ4f Cr52f * ; +EVALUATE Cr53Z4f := CrZ4f Cr53f * ; +EVALUATE Cr54Z4f := CrZ4f Cr54f * ; + +EVALUATE Sn112Z4f := SnZ4f Sn112f * ; +EVALUATE Sn114Z4f := SnZ4f Sn114f * ; +EVALUATE Sn115Z4f := SnZ4f Sn115f * ; +EVALUATE Sn116Z4f := SnZ4f Sn116f * ; +EVALUATE Sn117Z4f := SnZ4f Sn117f * ; +EVALUATE Sn118Z4f := SnZ4f Sn118f * ; +EVALUATE Sn119Z4f := SnZ4f Sn119f * ; +EVALUATE Sn120Z4f := SnZ4f Sn120f * ; +EVALUATE Sn122Z4f := SnZ4f Sn122f * ; +EVALUATE Sn124Z4f := SnZ4f Sn124f * ; + +*---- Zirc4 Table Creation ---- + +Zirc4 := UTL: :: CREA 'AvgMolM ' 1 = <<MZirc4>> ; +Zirc4 := UTL: Zirc4 :: CREA 'O16Z4frIso ' 1 = <<O16Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Fe54Z4frIso ' 1 = <<Fe54Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Fe56Z4frIso ' 1 = <<Fe56Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Fe57Z4frIso ' 1 = <<Fe57Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Fe58Z4frIso ' 1 = <<Fe58Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Cr50Z4frIso ' 1 = <<Cr50Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Cr52Z4frIso ' 1 = <<Cr52Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Cr53Z4frIso ' 1 = <<Cr53Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Cr54Z4frIso ' 1 = <<Cr54Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'ZrNatfrIso ' 1 = <<ZrZ4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Zr90Z4frIso ' 1 = <<Zr90Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Zr91Z4frIso ' 1 = <<Zr91Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Zr92Z4frIso ' 1 = <<Zr92Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Zr93Z4frIso ' 1 = <<Zr93Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Zr94Z4frIso ' 1 = <<Zr94Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Zr95Z4frIso ' 1 = <<Zr95Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Zr96Z4frIso ' 1 = <<Zr96Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn112Z4frIso' 1 = <<Sn112Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn114Z4frIso' 1 = <<Sn114Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn115Z4frIso' 1 = <<Sn115Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn116Z4frIso' 1 = <<Sn116Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn117Z4frIso' 1 = <<Sn117Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn118Z4frIso' 1 = <<Sn118Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn119Z4frIso' 1 = <<Sn119Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn120Z4frIso' 1 = <<Sn120Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn122Z4frIso' 1 = <<Sn122Z4f>> ; +Zirc4 := UTL: Zirc4 :: CREA 'Sn124Z4frIso' 1 = <<Sn124Z4f>> ; + +*------------------------------------------------------ +* Incon : Inconel - AFA 2G +* Massic Composition : NiNat .5250 +* CrNat .1900 +* FeNat .1823 +* BNatB .00003 +* Al27 .0050 +* Mo95 .0305 +* TiNat .0088 +* SiNat .0018 +* Co59 .0005 +* Mn55 .0018 +* Nb93 .05175 +* (&IMPU .00252) +*------------------------------------------------------ +REAL NiIncf := .5250 ; +REAL CrIncf := .1900 ; +REAL FeIncf := .1823 ; +REAL BIncf := .00003 ; +REAL AlIncf := .0050 ; +REAL MoIncf := .0305 ; +REAL TiIncf := .0088 ; +REAL SiIncf := .0018 ; +REAL CoIncf := .0005 ; +REAL MnIncf := .0018 ; +REAL NbIncf := .05175 ; + +IF Gestion "BEAV" = THEN +EVALUATE NiIncf := .5119 ; +EVALUATE CrIncf := .1896 ; +EVALUATE FeIncf := .2863 ; +EVALUATE BIncf := 1E-15 ; +EVALUATE AlIncf := 1E-15 ; +EVALUATE MoIncf := 1E-15 ; +EVALUATE TiIncf := 1E-15 ; +EVALUATE SiIncf := .0035 ; +EVALUATE CoIncf := .0005 ; +EVALUATE MnIncf := .0087 ; +EVALUATE NbIncf := 1E-15 ; +ENDIF ; + +REAL Ni58Inf Ni60Inf Ni61Inf Ni62Inf Ni64Inf ; +REAL Cr50Inf Cr52Inf Cr53Inf Cr54Inf ; +REAL Fe54Inf Fe56Inf Fe57Inf Fe58Inf ; +REAL Ti46Inf Ti47Inf Ti48Inf Ti49Inf Ti50Inf ; +REAL Si28Inf Si29Inf Si30Inf ; +REAL B10Inf B11Inf ; + +*---- Fraction Sum ---- +EVALUATE SumFr := NiIncf CrIncf + FeIncf + BIncf + AlIncf + MoIncf + ; +EVALUATE SumFr := SumFr TiIncf + SiIncf + CoIncf + MnIncf + NbIncf + ; + +*---- Mass Ratios -> Isotopic Ratios ---- +EVALUATE NiIncf := NiIncf avo_conc * MNiNat / ; +EVALUATE CrIncf := CrIncf avo_conc * MCrNat / ; +EVALUATE FeIncf := FeIncf avo_conc * MFeNat / ; +EVALUATE BIncf := BIncf avo_conc * MBNat / ; +EVALUATE AlIncf := AlIncf avo_conc * al27 / ; +EVALUATE MoIncf := MoIncf avo_conc * mo95 / ; +EVALUATE TiIncf := TiIncf avo_conc * MTiNat / ; +EVALUATE SiIncf := SiIncf avo_conc * MSiNat / ; +EVALUATE CoIncf := CoIncf avo_conc * co59 / ; +EVALUATE MnIncf := MnIncf avo_conc * mn55 / ; +EVALUATE NbIncf := NbIncf avo_conc * nb93 / ; + +EVALUATE ReNorm := NiIncf CrIncf + FeIncf + BIncf + AlIncf + MoIncf + ; +EVALUATE ReNorm := ReNorm TiIncf + SiIncf + CoIncf + MnIncf + NbIncf + ; +EVALUATE ReNorm := ReNorm SumFr / ; + +EVALUATE NiIncf := NiIncf ReNorm / ; +EVALUATE CrIncf := CrIncf ReNorm / ; +EVALUATE FeIncf := FeIncf ReNorm / ; +EVALUATE BIncf := BIncf ReNorm / ; +EVALUATE AlIncf := AlIncf ReNorm / ; +EVALUATE MoIncf := MoIncf ReNorm / ; +EVALUATE TiIncf := TiIncf ReNorm / ; +EVALUATE SiIncf := SiIncf ReNorm / ; +EVALUATE CoIncf := CoIncf ReNorm / ; +EVALUATE MnIncf := MnIncf ReNorm / ; +EVALUATE NbIncf := NbIncf ReNorm / ; + +*---- Isotopic Ratios for all Ni, Cr, Fe & B ---- +EVALUATE Ni58Inf := Ni58f NiIncf * ; +EVALUATE Ni60Inf := Ni60f NiIncf * ; +EVALUATE Ni61Inf := Ni61f NiIncf * ; +EVALUATE Ni62Inf := Ni62f NiIncf * ; +EVALUATE Ni64Inf := Ni64f NiIncf * ; + +EVALUATE Ti46Inf := Ti46f TiIncf * ; +EVALUATE Ti47Inf := Ti47f TiIncf * ; +EVALUATE Ti48Inf := Ti48f TiIncf * ; +EVALUATE Ti49Inf := Ti49f TiIncf * ; +EVALUATE Ti50Inf := Ti50f TiIncf * ; + +EVALUATE Si28Inf := Si28f SiIncf * ; +EVALUATE Si29Inf := Si29f SiIncf * ; +EVALUATE Si30Inf := Si30f SiIncf * ; + +EVALUATE Cr50Inf := Cr50f CrIncf * ; +EVALUATE Cr52Inf := Cr52f CrIncf * ; +EVALUATE Cr53Inf := Cr53f CrIncf * ; +EVALUATE Cr54Inf := Cr54f CrIncf * ; + +EVALUATE Fe54Inf := Fe54f FeIncf * ; +EVALUATE Fe56Inf := Fe56f FeIncf * ; +EVALUATE Fe57Inf := Fe57f FeIncf * ; +EVALUATE Fe58Inf := Fe58f FeIncf * ; + +EVALUATE B10Inf := B10f BIncf * ; +EVALUATE B11Inf := B11f BIncf * ; + +*---- Inconel Average Weight Ratio --- +REAL MIncon ; ! Average Weight Ratio +EVALUATE MIncon := NiIncf MNiNat * CrIncf MCrNat * + + FeIncf MFeNat * + BIncf MBNat * + + AlIncf al27 * + MoIncf mo95 * + + TiIncf MTiNat * + SiIncf MSiNat * + + CoIncf co59 * + MnIncf mn55 * + + NbIncf nb93 * + ; + +ECHO "Masse molaire de l'inconel " MIncon ; + +*---- Due to &IMPU, fractions are summing to "1-&IMPU". +* These fractions are stored, as they represent the correct proportion. +* MIncon however has to be renormalized so that fractions are summing to 1. ---- +EVALUATE MIncon := MIncon SumFr / ; + +*---- Incon Table Creation ---- + +Incon := UTL: :: CREA 'AvgMolM ' 1 = <<MIncon>> ; +Incon := UTL: Incon :: CREA 'B10InfrIso ' 1 = <<B10Inf>> ; +Incon := UTL: Incon :: CREA 'B11InfrIso ' 1 = <<B11Inf>> ; +Incon := UTL: Incon :: CREA 'Al27InfrIso ' 1 = <<AlIncf>> ; +Incon := UTL: Incon :: CREA 'SiNatInfrIso' 1 = <<SiIncf>> ; +Incon := UTL: Incon :: CREA 'Si28InfrIso ' 1 = <<Si28Inf>> ; +Incon := UTL: Incon :: CREA 'Si29InfrIso ' 1 = <<Si29Inf>> ; +Incon := UTL: Incon :: CREA 'Si30InfrIso ' 1 = <<Si30Inf>> ; +Incon := UTL: Incon :: CREA 'TiNatInfrIso' 1 = <<TiIncf>> ; +Incon := UTL: Incon :: CREA 'Ti46InfrIso ' 1 = <<Ti46Inf>> ; +Incon := UTL: Incon :: CREA 'Ti47InfrIso ' 1 = <<Ti47Inf>> ; +Incon := UTL: Incon :: CREA 'Ti48InfrIso ' 1 = <<Ti48Inf>> ; +Incon := UTL: Incon :: CREA 'Ti49InfrIso ' 1 = <<Ti49Inf>> ; +Incon := UTL: Incon :: CREA 'Ti50InfrIso ' 1 = <<Ti50Inf>> ; +Incon := UTL: Incon :: CREA 'Cr50InfrIso ' 1 = <<Cr50Inf>> ; +Incon := UTL: Incon :: CREA 'Cr52InfrIso ' 1 = <<Cr52Inf>> ; +Incon := UTL: Incon :: CREA 'Cr53InfrIso ' 1 = <<Cr53Inf>> ; +Incon := UTL: Incon :: CREA 'Cr54InfrIso ' 1 = <<Cr54Inf>> ; +Incon := UTL: Incon :: CREA 'Mn55InfrIso ' 1 = <<MnIncf>> ; +Incon := UTL: Incon :: CREA 'Fe54InfrIso ' 1 = <<Fe54Inf>> ; +Incon := UTL: Incon :: CREA 'Fe56InfrIso ' 1 = <<Fe56Inf>> ; +Incon := UTL: Incon :: CREA 'Fe57InfrIso ' 1 = <<Fe57Inf>> ; +Incon := UTL: Incon :: CREA 'Fe58InfrIso ' 1 = <<Fe58Inf>> ; +Incon := UTL: Incon :: CREA 'Co59InfrIso ' 1 = <<CoIncf>> ; +Incon := UTL: Incon :: CREA 'Ni58InfrIso ' 1 = <<Ni58Inf>> ; +Incon := UTL: Incon :: CREA 'Ni60InfrIso ' 1 = <<Ni60Inf>> ; +Incon := UTL: Incon :: CREA 'Ni61InfrIso ' 1 = <<Ni61Inf>> ; +Incon := UTL: Incon :: CREA 'Ni62InfrIso ' 1 = <<Ni62Inf>> ; +Incon := UTL: Incon :: CREA 'Ni64InfrIso ' 1 = <<Ni64Inf>> ; +Incon := UTL: Incon :: CREA 'Nb93InfrIso ' 1 = <<NbIncf>> ; +Incon := UTL: Incon :: CREA 'Mo95InfrIso ' 1 = <<MoIncf>> ; + +*------------------------------------------------------ +* SS304 : control rod clad -- Stainless Steel +* Massic Composition : FeNat .706 +* NiNat .100 +* CrNat .180 +* Mn55 .010 +* SiNat .004 +*------------------------------------------------------ +REAL FeSSf := .706 ; +REAL NiSSf := .100 ; +REAL CrSSf := .180 ; +REAL MnSSf := .010 ; +REAL SiSSf := .004 ; + +IF Gestion "BEAV" = THEN + EVALUATE FeSSf := .684 ; + EVALUATE NiSSf := .100 ; + EVALUATE CrSSf := .190 ; + EVALUATE MnSSf := .020 ; + EVALUATE SiSSf := .006 ; +ENDIF ; + +REAL Fe54SSf Fe56SSf Fe57SSf Fe58SSf ; +REAL Ni58SSf Ni60SSf Ni61SSf Ni62SSf Ni64SSf ; +REAL Cr50SSf Cr52SSf Cr53SSf Cr54SSf ; +REAL Si28SSf Si29SSf Si30SSf ; + +*---- Mass Ratios -> Isotopic Ratios ---- +EVALUATE FeSSf := FeSSf avo_conc * MFeNat / ; +EVALUATE NiSSf := NiSSf avo_conc * MNiNat / ; +EVALUATE CrSSf := CrSSf avo_conc * MCrNat / ; +EVALUATE MnSSf := MnSSf avo_conc * mn55 / ; +EVALUATE SiSSf := SiSSf avo_conc * MSiNat / ; + +EVALUATE ReNorm := FeSSf NiSSf + CrSSf + MnSSf + SiSSf + ; + +EVALUATE FeSSf := FeSSf ReNorm / ; +EVALUATE NiSSf := NiSSf ReNorm / ; +EVALUATE CrSSf := CrSSf ReNorm / ; +EVALUATE MnSSf := MnSSf ReNorm / ; +EVALUATE SiSSf := SiSSf ReNorm / ; + +*---- SS304 Average Weight Ratio --- +REAL MSS304 ; ! Average Weight Ratio +EVALUATE MSS304 := FeSSf MFeNat * NiSSf MNiNat * + + CrSSf MCrNat * + MnSSf mn55 * + + SiSSf MSiNat * + ; + +ECHO "Masse molaire du SS304 " MSS304 ; + +*---- Isotopic Ratios for all Fe, Ni & Cr ---- +EVALUATE Fe54SSf := FeSSf Fe54f * ; +EVALUATE Fe56SSf := FeSSf Fe56f * ; +EVALUATE Fe57SSf := FeSSf Fe57f * ; +EVALUATE Fe58SSf := FeSSf Fe58f * ; + +EVALUATE Ni58SSf := NiSSf Ni58f * ; +EVALUATE Ni60SSf := NiSSf Ni60f * ; +EVALUATE Ni61SSf := NiSSf Ni61f * ; +EVALUATE Ni62SSf := NiSSf Ni62f * ; +EVALUATE Ni64SSf := NiSSf Ni64f * ; + +EVALUATE Cr50SSf := CrSSf Cr50f * ; +EVALUATE Cr52SSf := CrSSf Cr52f * ; +EVALUATE Cr53SSf := CrSSf Cr53f * ; +EVALUATE Cr54SSf := CrSSf Cr54f * ; + +EVALUATE Si28SSf := SiSSf Si28f * ; +EVALUATE Si29SSf := SiSSf Si29f * ; +EVALUATE Si30SSf := SiSSf Si30f * ; + +*---- SS304 Table Creation ---- + +* Isotopes non conflictuels +SS304 := UTL: :: CREA 'AvgMolM ' 1 = <<MSS304>> ; +SS304 := UTL: SS304 :: CREA 'SiNatSSfrIso' 1 = <<SiSSf>> ; +SS304 := UTL: SS304 :: CREA 'Si28SSfrIso ' 1 = <<Si28SSf>> ; +SS304 := UTL: SS304 :: CREA 'Si29SSfrIso ' 1 = <<Si29SSf>> ; +SS304 := UTL: SS304 :: CREA 'Si30SSfrIso ' 1 = <<Si30SSf>> ; +SS304 := UTL: SS304 :: CREA 'Fe54SSfrIso ' 1 = <<Fe54SSf>> ; +SS304 := UTL: SS304 :: CREA 'Fe56SSfrIso ' 1 = <<Fe56SSf>> ; +SS304 := UTL: SS304 :: CREA 'Fe57SSfrIso ' 1 = <<Fe57SSf>> ; +SS304 := UTL: SS304 :: CREA 'Fe58SSfrIso ' 1 = <<Fe58SSf>> ; +SS304 := UTL: SS304 :: CREA 'Ni58SSfrIso ' 1 = <<Ni58SSf>> ; +SS304 := UTL: SS304 :: CREA 'Ni60SSfrIso ' 1 = <<Ni60SSf>> ; +SS304 := UTL: SS304 :: CREA 'Ni61SSfrIso ' 1 = <<Ni61SSf>> ; +SS304 := UTL: SS304 :: CREA 'Ni62SSfrIso ' 1 = <<Ni62SSf>> ; +SS304 := UTL: SS304 :: CREA 'Ni64SSfrIso ' 1 = <<Ni64SSf>> ; +SS304 := UTL: SS304 :: CREA 'Cr50SSfrIso ' 1 = <<Cr50SSf>> ; +SS304 := UTL: SS304 :: CREA 'Cr52SSfrIso ' 1 = <<Cr52SSf>> ; +SS304 := UTL: SS304 :: CREA 'Cr53SSfrIso ' 1 = <<Cr53SSf>> ; +SS304 := UTL: SS304 :: CREA 'Cr54SSfrIso ' 1 = <<Cr54SSf>> ; +SS304 := UTL: SS304 :: CREA 'Mn55SSfrIso ' 1 = <<MnSSf>> ; + +*------------------------------------------------------ +* M5 +* Ref: C. Lemaignan "Zirconium Alloys: Properties +* and Characteristics" +* Massic Composition : FeNat .00050 +* O16 .00135 +* ZrNat .98813 +* Nb93 .01000 +* SNat .00002 +*------------------------------------------------------ +REAL FeM5f := 0.00050 ; +REAL O16M5f := 0.00135 ; +REAL ZrM5f := 0.98813 ; +REAL Nb93M5f := 0.01000 ; +REAL SM5f := 0.00002 ; + +REAL Zr90M5f Zr91M5f Zr92M5f Zr93M5f Zr94M5f Zr95M5f Zr96M5f ; +REAL Fe54M5f Fe56M5f Fe57M5f Fe58M5f ; +REAL S32M5f S33M5f S34M5f ; + +*---- Mass Ratios -> Isotopic Ratios ---- +EVALUATE FeM5f := FeM5f avo_conc * MFeNat / ; +EVALUATE O16M5f := O16M5f avo_conc * o16 / ; +EVALUATE ZrM5f := ZrM5f avo_conc * MZrNat / ; +EVALUATE Nb93M5f := Nb93M5f avo_conc * nb93 / ; +EVALUATE SM5f := SM5f avo_conc * MSNat / ; + +EVALUATE ReNorm := FeM5f O16M5f + ZrM5f + Nb93M5f + SM5f + ; + +EVALUATE FeM5f := FeM5f ReNorm / ; +EVALUATE O16M5f := O16M5f ReNorm / ; +EVALUATE ZrM5f := ZrM5f ReNorm / ; +EVALUATE Nb93M5f := Nb93M5f ReNorm / ; +EVALUATE SM5f := SM5f ReNorm / ; + +*---- M5 Average Weight Ratio ---- +REAL MM5 ; ! Average Weight Ratio +EVALUATE MM5 := FeM5f MFeNat * O16M5f o16 * + ; +EVALUATE MM5 := MM5 ZrM5f MZrNat * + Nb93M5f nb93 * + ; +EVALUATE MM5 := MM5 SM5f MSNat * + ; + +*---- Isotopic Ratios for all Zr & Fe & S ---- +EVALUATE Zr90M5f := ZrM5f Zr90f * ; +EVALUATE Zr91M5f := ZrM5f Zr91f * ; +EVALUATE Zr92M5f := ZrM5f Zr92f * ; +EVALUATE Zr93M5f := ZrM5f Zr93f * ; +EVALUATE Zr94M5f := ZrM5f Zr94f * ; +EVALUATE Zr95M5f := ZrM5f Zr95f * ; +EVALUATE Zr96M5f := ZrM5f Zr96f * ; + +EVALUATE Fe54M5f := FeM5f Fe54f * ; +EVALUATE Fe56M5f := FeM5f Fe56f * ; +EVALUATE Fe57M5f := FeM5f Fe57f * ; +EVALUATE Fe58M5f := FeM5f Fe58f * ; + +EVALUATE S32M5f := SM5f S32f * ; +EVALUATE S33M5f := SM5f S33f * ; +EVALUATE S34M5f := SM5f S34f * ; + +*---- M5 Table Creation ---- +M5 := UTL: :: CREA 'AvgMolM ' 1 = <<MM5>> ; +M5 := UTL: M5 :: CREA 'Fe54M5frIso' 1 = <<Fe54M5f>> ; +M5 := UTL: M5 :: CREA 'Fe56M5frIso' 1 = <<Fe56M5f>> ; +M5 := UTL: M5 :: CREA 'Fe57M5frIso' 1 = <<Fe57M5f>> ; +M5 := UTL: M5 :: CREA 'Fe58M5frIso' 1 = <<Fe58M5f>> ; +M5 := UTL: M5 :: CREA 'O16M5frIso ' 1 = <<O16M5f>> ; +M5 := UTL: M5 :: CREA 'ZrM5frIso ' 1 = <<ZrM5f>> ; +M5 := UTL: M5 :: CREA 'Zr90M5frIso' 1 = <<Zr90M5f>> ; +M5 := UTL: M5 :: CREA 'Zr91M5frIso' 1 = <<Zr91M5f>> ; +M5 := UTL: M5 :: CREA 'Zr92M5frIso' 1 = <<Zr92M5f>> ; +M5 := UTL: M5 :: CREA 'Zr93M5frIso' 1 = <<Zr93M5f>> ; +M5 := UTL: M5 :: CREA 'Zr94M5frIso' 1 = <<Zr94M5f>> ; +M5 := UTL: M5 :: CREA 'Zr95M5frIso' 1 = <<Zr95M5f>> ; +M5 := UTL: M5 :: CREA 'Zr96M5frIso' 1 = <<Zr96M5f>> ; +M5 := UTL: M5 :: CREA 'Nb93M5frIso' 1 = <<Nb93M5f>> ; +M5 := UTL: M5 :: CREA 'S32M5frIso ' 1 = <<S32M5f>> ; +M5 := UTL: M5 :: CREA 'S33M5frIso ' 1 = <<S33M5f>> ; +M5 := UTL: M5 :: CREA 'S34M5frIso ' 1 = <<S34M5f>> ; + +*------------------------------------------------------ +* AIC : control rod neutron absorber -- AgInCd +* Massic Composition : AgNat .80611 +* CdNat .04909 +* In115 (1-(Ag+Cd)) +* Dans le cas de ENDFVI8 : In115 est en realite InNat +*------------------------------------------------------ +REAL AgAICf := .80611 ; +REAL CdAICf := .04909 ; +REAL InAICf := 1. AgAICf - CdAICf - ; + +REAL Ag107AICf Ag109AICf ; +REAL Cd106AICf Cd108AICf Cd110AICf Cd111AICf Cd112AICf Cd113AICf ; +REAL Cd114AICf Cd116AICf ; +REAL In115AICf In113AICf ; + +*---- Mass Ratios -> Isotopic Ratios ---- +EVALUATE AgAICf := AgAICf avo_conc * MAgNat / ; +EVALUATE CdAICf := CdAICf avo_conc * MCdNat / ; +*EVALUATE InAICf := InAICf avo_conc * in115 / ; +EVALUATE InAICf := InAICf avo_conc * MInNat / ; + +EVALUATE ReNorm := AgAICf CdAICf + InAICf + ; + +EVALUATE AgAICf := AgAICf ReNorm / ; +EVALUATE CdAICf := CdAICf ReNorm / ; +EVALUATE InAICf := InAICf ReNorm / ; + +*---- AIC Average Weight Ratio --- +REAL MAIC ; ! Average Weight Ratio +EVALUATE MAIC := AgAICf MAgNat * CdAICf MCdNat * + ; +*EVALUATE MAIC := MAIC InAICf in115 * + ; +EVALUATE MAIC := MAIC InAICf MInNat * + ; + +*---- Isotopic Ratios for all Ag & Cd ---- +EVALUATE Ag107AICf := AgAICf Ag107f * ; +EVALUATE Ag109AICf := AgAICf Ag109f * ; + +EVALUATE Cd106AICf := CdAICf Cd106f * ; +EVALUATE Cd108AICf := CdAICf Cd108f * ; +EVALUATE Cd110AICf := CdAICf Cd110f * ; +EVALUATE Cd111AICf := CdAICf Cd111f * ; +EVALUATE Cd112AICf := CdAICf Cd112f * ; +EVALUATE Cd113AICf := CdAICf Cd113f * ; +EVALUATE Cd114AICf := CdAICf Cd114f * ; +EVALUATE Cd116AICf := CdAICf Cd116f * ; + +EVALUATE In115AICf := InAICf In115f * ; +EVALUATE In113AICf := InAICf In113f * ; + +*---- AIC Table Creation ---- +AIC := UTL: :: CREA 'AvgMolM ' 1 = <<MAIC>> ; +AIC := UTL: AIC :: CREA 'Ag107AICfIso' 1 = <<Ag107AICf>> ; +AIC := UTL: AIC :: CREA 'Ag109AICfIso' 1 = <<Ag109AICf>> ; +AIC := UTL: AIC :: CREA 'Cd106AICfIso' 1 = <<Cd106AICf>> ; +AIC := UTL: AIC :: CREA 'Cd108AICfIso' 1 = <<Cd108AICf>> ; +AIC := UTL: AIC :: CREA 'Cd110AICfIso' 1 = <<Cd110AICf>> ; +AIC := UTL: AIC :: CREA 'Cd111AICfIso' 1 = <<Cd111AICf>> ; +AIC := UTL: AIC :: CREA 'Cd112AICfIso' 1 = <<Cd112AICf>> ; +AIC := UTL: AIC :: CREA 'Cd113AICfIso' 1 = <<Cd113AICf>> ; +AIC := UTL: AIC :: CREA 'Cd114AICfIso' 1 = <<Cd114AICf>> ; +AIC := UTL: AIC :: CREA 'Cd116AICfIso' 1 = <<Cd116AICf>> ; +AIC := UTL: AIC :: CREA 'In115AICfIso' 1 = <<In115AICf>> ; +AIC := UTL: AIC :: CREA 'In113AICfIso' 1 = <<In113AICf>> ; + +*------------------------------------------------------ +* GdNat : Natural Gadolinium +* Isotopic Composition : Gd152 .002 *** +* Gd154 .021 *** +* Gd155 .148 +* Gd156 .206 +* Gd157 .157 +* Gd158 .248 +* Gd160 .218 +* +* *** : Gd152 is not available. It is replaced by Gd154 +* which has relatively close enough neutronic properties +* +* => Gd152 .000 +* => Gd154 .023 +*------------------------------------------------------ +*REAL Gd152f := .000 ; +REAL Gd154f := .023 ; +REAL Gd155f := .148 ; +REAL Gd156f := .206 ; +REAL Gd157f := .157 ; +REAL Gd158f := .248 ; +REAL Gd160f := .218 ; + +*---- GdNat Average Weight Ratio ---- +REAL MGdNat ; ! Average Weight Ratio +EVALUATE MGdNat := gd154 Gd154f * gd155 Gd155f * + ; +EVALUATE MGdNat := MGdNat gd156 Gd156f * + gd157 Gd157f * + ; +EVALUATE MGdNat := MGdNat gd158 Gd158f * + gd160 Gd160f * + ; + +*---- GdNat Table Creation ---- +GdNat := UTL: :: CREA 'AvgMolM ' 1 = <<MGdNat>> ; +GdNat := UTL: GdNat :: CREA 'Gd154frIso' 1 = <<Gd154f>> ; +GdNat := UTL: GdNat :: CREA 'Gd155frIso' 1 = <<Gd155f>> ; +GdNat := UTL: GdNat :: CREA 'Gd156frIso' 1 = <<Gd156f>> ; +GdNat := UTL: GdNat :: CREA 'Gd157frIso' 1 = <<Gd157f>> ; +GdNat := UTL: GdNat :: CREA 'Gd158frIso' 1 = <<Gd158f>> ; +GdNat := UTL: GdNat :: CREA 'Gd160frIso' 1 = <<Gd160f>> ; + +*------------------------------------------------------ +* Gd2O3 : Gadolinium Oxyde +* Isotopic Composition : GdNat 2 +* O16 3 +*------------------------------------------------------ +REAL GadOxGdf := 2. ; +REAL GadOxOf := 3. ; + +REAL GadOxG4f GadOxG5f GadOxG6f GadOxG7f GadOxG8f GadOxG0f ; + +EVALUATE GadOxG4f := GadOxGdf Gd154f * ; +EVALUATE GadOxG5f := GadOxGdf Gd155f * ; +EVALUATE GadOxG6f := GadOxGdf Gd156f * ; +EVALUATE GadOxG7f := GadOxGdf Gd157f * ; +EVALUATE GadOxG8f := GadOxGdf Gd158f * ; +EVALUATE GadOxG0f := GadOxGdf Gd160f * ; + +REAL MGd2O3 ; ! Average Weight Ratio +EVALUATE MGd2O3 := MGdNat GadOxGdf * o16 GadOxOf * + ; + +*---- Gd2O3 Table Creation ---- +Gd2O3 := UTL: :: CREA 'AvgMolM ' 1 = <<MGd2O3>> ; +Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG4frIso' 1 = <<GadOxG4f>> ; +Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG5frIso' 1 = <<GadOxG5f>> ; +Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG6frIso' 1 = <<GadOxG6f>> ; +Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG7frIso' 1 = <<GadOxG7f>> ; +Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG8frIso' 1 = <<GadOxG8f>> ; +Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxG0frIso' 1 = <<GadOxG0f>> ; +Gd2O3 := UTL: Gd2O3 :: CREA 'GadOxOfrIso ' 1 = <<GadOxOf>> ; + +*------------------------------------------------------ +* B2O3 : Solid White Glass Diboron Trioxide # B +* Quantity of BNat in B2O3 : BNatBq = 2. +* Quantity of O16 in B2O3 : O16Bq = 3. +*------------------------------------------------------ +REAL BNatBq := 2. ; +REAL O16Bq := 3. ; +*---- B2O3 Average Weight Ratio ---- +REAL MB2O3 ; ! Average Weight Ratio +EVALUATE MB2O3 := MBNat BNatBq * o16 O16Bq * + ; +*------------------------------------------------------ +* SiO2 : Silica Crystalline # S +* Quantity of SiNat in SiO2 : SiNatSq = 1. +* Quantity of O16 in SiO2 : O16Sq = 2. +*------------------------------------------------------ +REAL SiNatSq := 1. ; +REAL O16Sq := 2. ; +*---- SiO2 Average Weight Ratio ---- +REAL MSiO2 ; ! Average Weight Ratio +EVALUATE MSiO2 := MSiNat SiNatSq * o16 O16Sq * + ; +*------------------------------------------------------ +* Al2O3 : Aluminium Oxide # A +* Quantity of Aluminium in Al2O3 : Al27Aq = 2. +* Quantity of O16 in Al2O3 : O16Aq = 3. +*------------------------------------------------------ +REAL Al27Aq := 2. ; +REAL O16Aq := 3. ; +*---- Al2O3 Average Weight Ratio ---- +REAL MAl2O3 ; ! Average Weight Ratio +EVALUATE MAl2O3 := al27 Al27Aq * o16 O16Aq * + ; +*------------------------------------------------------ +* K2O : Potassium Oxide # K +* Quantity of KNat in K2O : KNatKq = 2. +* Quantity of O16 in K2O : O16Kq = 1. +*------------------------------------------------------ +REAL KNatKq := 2. ; +REAL O16Kq := 1. ; +*---- K2O Average Weight Ratio ---- +REAL MK2O ; ! Average Weight Ratio +EVALUATE MK2O := MKNat KNatKq * o16 O16Kq * + ; +*------------------------------------------------------ +* Na2O : Sodium oxide # N +* Quantity of Na23 in Na2O : Na23Nq = 2. +* Quantity of O16 in Na2O : O16Nq = 1. +*------------------------------------------------------ +REAL Na23Nq := 2. ; +REAL O16Nq := 1. ; +*---- Na2O Average Weight Ratio ---- +REAL MNa2O ; ! Average Weight Ratio +EVALUATE MNa2O := na23 Na23Nq * o16 O16Nq * + ; +*------------------------------------------------------ + +*------------------------------------------------------ +* BP : Burnable Poison Borosilicate Glass +* Tihange Specifications CEA-N-2092 +* Massic Composition : B2O3 .181 +* SiO2 .648 +* Al2O3 .076 +* K2O .032 +* Na2O .031 +* The sum of massic composition %m is equal to 0.968 !! +*------------------------------------------------------ +* Massic Fraction : %m +*------------------------------------------------------ +REAL B2O3BPf ; +REAL SiO2BPf ; +REAL Al2O3BPf ; +REAL K2OBPf ; +REAL Na2OBPf ; + +IF Gestion "BEAV" = THEN +* Weight Fraction : Onat 0.5481 +* AlNat 0.0344 +* B10 0.0071 +* B11 0.0317 +* Sinat 0.3787 + EVALUATE B2O3BPf := .125 ; + EVALUATE SiO2BPf := .810 ; ! .806 + EVALUATE Al2O3BPf := .065 ; ! .024 + EVALUATE K2OBPf := .0 ; ! .022 + EVALUATE Na2OBPf := .0 ; ! .021 +ELSE + EVALUATE B2O3BPf := .181 ; + EVALUATE SiO2BPf := .648 ; + EVALUATE Al2O3BPf := .076 ; + EVALUATE K2OBPf := .032 ; + EVALUATE Na2OBPf := .031 ; +ENDIF ; +EVALUATE SumFr := B2O3BPf SiO2BPf Al2O3BPf K2OBPf Na2OBPf + + + + ; + +* +*----------------------------------------------------- +* Material Massic To Material Fraction +* +* %m_i : Massic fraction of material i +* %n_i : Material fraction of material i +* +* %n_i = [%m_i / M_i] / [ Sum ( %m_k / M_k ) ] +*----------------------------------------------------- +* +EVALUATE B2O3BPf := B2O3BPf MB2O3 / ; +EVALUATE SiO2BPf := SiO2BPf MSiO2 / ; +EVALUATE Al2O3BPf := Al2O3BPf MAl2O3 / ; +EVALUATE K2OBPf := K2OBPf MK2O / ; +EVALUATE Na2OBPf := Na2OBPf MNa2O / ; +* +EVALUATE ReNorm := B2O3BPf SiO2BPf Al2O3BPf K2OBPf Na2OBPf + + + + ; +* +EVALUATE B2O3BPf := B2O3BPf ReNorm / ; +EVALUATE SiO2BPf := SiO2BPf ReNorm / ; +EVALUATE Al2O3BPf := Al2O3BPf ReNorm / ; +EVALUATE K2OBPf := K2OBPf ReNorm / ; +EVALUATE Na2OBPf := Na2OBPf ReNorm / ; + +* +*----------------------------------------------------- +* Molar Mass of BP Calculation : MBP +*----------------------------------------------------- +REAL MBP := 0.0 ; ! BP Average Weight Ratio +EVALUATE MBP := B2O3BPf MB2O3 * SiO2BPf MSiO2 * + ; +EVALUATE MBP := MBP Al2O3BPf MAl2O3 * + K2OBPf MK2O * + ; +EVALUATE MBP := MBP Na2OBPf MNa2O * + ; +* +*---------------------------------------------------- +* AMF : Ajusted Material Fraction +* because Sum of %m is 0.968 and not 1 +* --> AMF = Material Fraction (%n_i) * SumFr +*---------------------------------------------------- +EVALUATE B2O3BPf := B2O3BPf SumFr * ; +EVALUATE SiO2BPf := SiO2BPf SumFr * ; +EVALUATE Al2O3BPf := Al2O3BPf SumFr * ; +EVALUATE K2OBPf := K2OBPf SumFr * ; +EVALUATE Na2OBPf := Na2OBPf SumFr * ; +*----------------------------------------------------- +* Isotopic Ratio Calculation +*----------------------------------------------------- +* %f = AMF x Qem x Fie +* Qem = Quantity of Element in Material +* Fie = Fraction of Isotope In Element +* +* example : B10 isotopic fraction +* We find 0,20% of B10 in element BNat -> Fie = 0.2 +* There are 2 BNat in B203 -> Qem = 2. +* The molar fraction of B2O3 is %n -> B2O3BPf +* +* --> B10BPf = 0.2 * 2. * B2O3BPf +*----------------------------------------------------- +*------ Isotopic Ratios for all isotopes in BP ------- +REAL B10BPf ; +REAL B11BPf ; +REAL O16BPf ; +REAL SiNatBPf ; +REAL Al27BPf ; +REAL KNatBPf ; +REAL Na23BPf ; + +* +EVALUATE B10BPf := B10f BNatBq B2O3BPf * * ; +EVALUATE B11BPf := B11f BNatBq B2O3BPf * * ; +EVALUATE O16BPf := O16Bq B2O3BPf * O16Sq SiO2BPf * + ; +EVALUATE O16BPf := O16BPf O16Aq Al2O3BPf * + ; +EVALUATE O16BPf := O16BPf O16Kq K2OBPf * + ; +EVALUATE O16BPf := O16BPf O16Nq Na2OBPf * + ; +EVALUATE Na23BPf := Na23Nq Na2OBPf * ; +EVALUATE Al27BPf := Al27Aq Al2O3BPf * ; +EVALUATE KNatBPf := KNatKq K2OBPf * ; +REAL K39BPf := KNatKq K39f * K2OBPf * ; +REAL K40BPf := KNatKq K40f * K2OBPf * ; +REAL K41BPf := KNatKq K41f * K2OBPf * ; +EVALUATE SiNatBPf := SiNatSq SiO2BPf * ; +REAL Si28BPf := SiNatSq Si28f * SiO2BPf * ; +REAL Si29BPf := SiNatSq Si29f * SiO2BPf * ; +REAL Si30BPf := SiNatSq Si30f * SiO2BPf * ; +* +*---- BP Table Creation ---- +BP := UTL: :: CREA 'AvgMolM ' 1 = <<MBP>> ; +BP := UTL: BP :: CREA 'B10BPfrIso ' 1 = <<B10BPf>> ; +BP := UTL: BP :: CREA 'B11BPfrIso ' 1 = <<B11BPf>> ; +BP := UTL: BP :: CREA 'O16BPfrIso ' 1 = <<O16BPf>> ; +BP := UTL: BP :: CREA 'SiNatBPfrIso' 1 = <<SiNatBPf>> ; +BP := UTL: BP :: CREA 'Si28BPfrIso' 1 = <<Si28BPf>> ; +BP := UTL: BP :: CREA 'Si29BPfrIso' 1 = <<Si29BPf>> ; +BP := UTL: BP :: CREA 'Si30BPfrIso' 1 = <<Si30BPf>> ; +BP := UTL: BP :: CREA 'Al27BPfrIso ' 1 = <<Al27BPf>> ; +BP := UTL: BP :: CREA 'KNatBPfrIso ' 1 = <<KNatBPf>> ; +BP := UTL: BP :: CREA 'K39BPfrIso ' 1 = <<K39BPf>> ; +BP := UTL: BP :: CREA 'K40BPfrIso ' 1 = <<K40BPf>> ; +BP := UTL: BP :: CREA 'K41BPfrIso ' 1 = <<K41BPf>> ; +BP := UTL: BP :: CREA 'Na23BPfrIso ' 1 = <<Na23BPf>> ; + +END: ; +QUIT "LIST" . + |