From 68c03c65cfcfd66a3b7c324cfa8770133538f5b0 Mon Sep 17 00:00:00 2001 From: Alain Hebert Date: Tue, 23 Sep 2025 11:38:23 -0400 Subject: #5: Corrections and improvements to the Baff-Refl reflector model in Donjon --- Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m | 8 +++++--- Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m | 8 ++++---- 2 files changed, 9 insertions(+), 7 deletions(-) (limited to 'Donjon/data/Reflector_beavrs_APEX_DF-RT_proc') diff --git a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m index 48db04e..9b23d60 100755 --- a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m +++ b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/CreaDilut.c2m @@ -36,6 +36,7 @@ REAL avo_conc := avogadro 1.E-24 * ; REAL dSS := 7.792815 ; ! Densite de l'acier inoxydable (SS304) REAL dZr4 := 6.49012 ; ! Densite du zircaloy 4 REAL dInc := 8.303 ; ! Densite inconel (a 20 degres C...) +INTEGER nbZones ; ! Nombre de zones dans le reflecteur INTEGER iMil ; ! Numero de la zone reflecteur REAL TMil ; ! Temperature du moderateur en degres C (suppose = T de tout le milieu) REAL dZone ; ! Densite de la zone consideree @@ -49,7 +50,7 @@ INTEGER LibId ; *-------------------------------------------------------------------- * Fractions volumiques des differents milieux presents dans la zone *-------------------------------------------------------------------- -:: >>dMod<< >>iMil<< >>TMil<< >>CB<< >>fvMod<< >>fvZr4<< >>fvInc<< >>fvSS304<< +:: >>dMod<< >>nbZones<< >>iMil<< >>TMil<< >>CB<< >>fvMod<< >>fvZr4<< >>fvInc<< >>fvSS304<< >>NuclData<< >>LibType<< ; *------------------------------------------------------ @@ -423,9 +424,10 @@ ECHO "Nb93 " Nb93 ; * Definition du melange correspondant dans la librairie *-------------------------------------------------------------------- ECHO "CreaDilut: call LIB" ; +INTEGER nmix := nbZones 1 + ; IF iMil 1 = THEN LibDilut := LIB: :: - MXIS 300 NMIX 8 + MXIS 300 NMIX <> MIXS LIB: <> FIL: <> * Milieu generique issu de la dilution de tous les constituants MIX 1 <> @@ -465,7 +467,7 @@ IF iMil 1 = THEN ; ELSE LibDilut := LIB: LibDilut :: - MXIS 300 NMIX 8 + MXIS 300 NMIX <> MIXS LIB: <> FIL: <> * Milieu generique issu de la dilution de tous les constituants MIX <> <> diff --git a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m index 6a2f2cc..8d55fb4 100755 --- a/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m +++ b/Donjon/data/Reflector_beavrs_APEX_DF-RT_proc/REFL_APEX.c2m @@ -7,7 +7,7 @@ * CALL : * * ApexRefl := REFL_APEX Geom FVol Param :: <> <> * * <> <> <> <> * -* <> <> <> ; * +* <> <> <> <> ; * * * * Input objects: * * Geom : macro-geometry LCM object * @@ -25,6 +25,7 @@ * OptSph : 'SPH'/'NOSP' transform DF into ADF/keep DF * * OptNoal : 'ALBE'/'NOAL' compute albedo/force right albedo to 0 * * OptSS : 'TONE'/'NOSS' do/don't activate self-shielding of Fe56 and Cr52 * +* nbZones : number of zones in the reflector * * * * Output objects: * * ApexRefl : reflector APEX file * @@ -38,6 +39,7 @@ STRING nameFeed nameRefl Palier htype ; INTEGER nlf ; :: >>nlf<< ; STRING LibType NuclData OptSph OptNoal OptSS ; :: >>LibType<< >>NuclData<< >>OptSph<< >>OptNoal<< >>OptSS<< ; +INTEGER nbZones ; :: >>nbZones<< ; PROCEDURE GetMolarMass MetalMix CreaDilut assertS ; LINKED_LIST DONNEES Geom_refl LibDilut MACRO LibFuel MolarMasses FeNat CrNat SnNat @@ -166,8 +168,6 @@ TRACK_refl := SNT: Geom_refl :: MAXR 100000 SN 16 SCAT <> ; -INTEGER nbZones := 7 ; - REAL dModBas := 0.753 ; ! Densite du moderateur en entree du coeur REAL TMilBas := 286.4 ; ! Temp du reflecteur bas (soit Tmod,entree) EVALUATE TMilBas := TMilBas 273.15 + ; ! Conversion en Kelvin @@ -251,7 +251,7 @@ WHILE IndexCB iBoreMax < DO ! Define MIX 1 in reflector macrolib LibDilut := CreaDilut LibDilut MolarMasses BNat Zirc4 Incon SS304 :: - <> <> <> <> <> <> <> + <> <> <> <> <> <> <> <> <> <> <> ; EVALUATE iZone := iZone 1 + ; -- cgit v1.2.3