diff options
Diffstat (limited to 'Donjon/data/Fessenheim_proc/MarguetRefl.c2m')
| -rwxr-xr-x | Donjon/data/Fessenheim_proc/MarguetRefl.c2m | 111 |
1 files changed, 111 insertions, 0 deletions
diff --git a/Donjon/data/Fessenheim_proc/MarguetRefl.c2m b/Donjon/data/Fessenheim_proc/MarguetRefl.c2m new file mode 100755 index 0000000..477e952 --- /dev/null +++ b/Donjon/data/Fessenheim_proc/MarguetRefl.c2m @@ -0,0 +1,111 @@ +***************************************************************** +* * +* Procedure : MarguetRefl.c2m * +* Purpose : For a given boron concentration, obtain radial, * +* bottom and top reflector cross sections and * +* diffusion coefficients, as written in Serge * +* Marguet's book (see ref. below), with * +* D1 adjustement or without. * +* Author : V. Salino * +* Date : 09/2021 * +* * +* CALL : * +* MacroRefl := MarguetRefl :: <<CB>> <<TypeRefl>> * +* * +* TypeRefl can be either "MargCpyAdjus" or "MargCpy" * +* * +* Source : "La physique des réacteurs nucléaires", * +* S. Marguet, 2nd édition, p.864. * +* * +***************************************************************** +*---- +* Definition STRUCTURES, MODULES et PROCEDURES +*---- +PARAMETER MacroRefl :: ::: LINKED_LIST MacroRefl ; ; +MODULE MAC: ABORT: END: ; + +REAL CB ; +STRING TypeRefl ; +:: >>CB<< >>TypeRefl<< ; + +REAL D1h D1b D2h D2b := 1.30 1.30 0.40 0.40 ; +REAL D1r ; +IF TypeRefl "MargCpyAdjus" = THEN + EVALUATE D1r := 1.78 ; +ELSEIF TypeRefl "MargCpy" = THEN + EVALUATE D1r := 1.30 ; +ELSE + ECHO "Unknown reflector modelling" ; + ABORT: ; +ENDIF ; +REAL D2r := 0.40 ; +REAL SigmaRalenth SigmaRalentb SigmaRalentr := + 0.0218 0.0210 0.00474 ; + +* Boron dependency +REAL densh densb densr := 0.667 0.754 0.714 ; +REAL Voliqh Voliqb Voliqr := 0.61796 0.66922 0.74273 ; + +REAL SigmaBoAbs1h := 38.82 ; +REAL SigmaBoAbs1b := 2.28 ; +REAL SigmaBoAbs1r := 15.11 ; + +REAL SigmaBoAbs2h := 412.23 ; +REAL SigmaBoAbs2b := 344.97 ; +REAL SigmaBoAbs2r := 63.30 ; + +REAL b10 := 10.01294 ; +REAL b11 := 11.00928 ; +REAL B10f := .2 ; +REAL B11f := .8 ; +REAL MBNat := b10 B10f * b11 B11f * + ; + +REAL avogadro := 6.022094E+23 ; +REAL avo_conc := avogadro 1.E-24 * ; +REAL BTot := CB 1.0E-6 * avo_conc * MBNat / ; + +REAL SigmaAbs1h := 0.0119 SigmaBoAbs1h BTot * densh * Voliqh * + ; +REAL SigmaAbs1b := 0.0157 SigmaBoAbs1b BTot * densb * Voliqb * + ; +REAL SigmaAbs1r := 0.0097 SigmaBoAbs1r BTot * densr * Voliqr * + ; + +REAL SigmaAbs2h := 0.0176 SigmaBoAbs2h BTot * densh * Voliqh * + ; +REAL SigmaAbs2b := 0.0105 SigmaBoAbs2b BTot * densb * Voliqb * + ; +REAL SigmaAbs2r := 0.0831 SigmaBoAbs2r BTot * densr * Voliqr * + ; + +* p.863 +EVALUATE SigmaRalenth := SigmaRalenth SigmaAbs1h 0.0119 - - ; +EVALUATE SigmaRalentb := SigmaRalentb SigmaAbs1b 0.0157 - - ; +EVALUATE SigmaRalentr := SigmaRalentr SigmaAbs1r 0.0097 - - ; + +REAL Total1h := SigmaAbs1h SigmaRalenth + ; +REAL Total2h := SigmaAbs2h ; +REAL Total1b := SigmaAbs1b SigmaRalentb + ; +REAL Total2b := SigmaAbs2b ; +REAL Total1r := SigmaAbs1r SigmaRalentr + ; +REAL Total2r := SigmaAbs2r ; + +MacroRefl := MAC: :: + EDIT 3 NGRO 2 NMIX 3 NIFI 0 ANIS 1 CTRA NONE NADF 1 + ENER 2.0E7 0.625 1.0E-3 VOLUME 1.0 1.0 1.0 READ INPUT + MIX 1 ! Radial + FLUX-INTG 1.0 1.0 + TOTAL <<Total1r>> <<Total2r>> + DIFF <<D1r>> <<D2r>> + SCAT 2 2 (*2->1*) 0.0 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) <<SigmaRalentr>> + MIX 2 ! Bottom + FLUX-INTG 1.0 1.0 + TOTAL <<Total1b>> <<Total2b>> + DIFF <<D1b>> <<D2b>> + SCAT 2 2 (*2->1*) 0.0 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) <<SigmaRalentb>> + MIX 3 ! Top + FLUX-INTG 1.0 1.0 + TOTAL <<Total1h>> <<Total2h>> + DIFF <<D1h>> <<D2h>> + SCAT 2 2 (*2->1*) 0.0 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) <<SigmaRalenth>> + ; + +END: ; +QUIT "LIST" . |