diff options
Diffstat (limited to 'Donjon/data/VHTRc.x2m')
| -rw-r--r-- | Donjon/data/VHTRc.x2m | 423 |
1 files changed, 423 insertions, 0 deletions
diff --git a/Donjon/data/VHTRc.x2m b/Donjon/data/VHTRc.x2m new file mode 100644 index 0000000..bdc1150 --- /dev/null +++ b/Donjon/data/VHTRc.x2m @@ -0,0 +1,423 @@ +*********************************************************** +* * +* Input file : core10.x2m * +* Purpose : Full-core calculation for a VHTR * +* Author(s) : Vincent Descotes (2010) * +* * +*********************************************************** +PROCEDURE Pgeohex assertS ; +MODULE DELETE: END: LIB: NCR: MACINI: USPLIT: TRIVAT: TRIVAA: + FLUD: FLPOW: RESINI: GREP: ; +LINKED_LIST GEOM MATEX MACRO1 MACRO2 CPO1 TRACK SYSTEM FLUX POWER + FLMAP MACFL ; +REAL Keff mbp mcp ; +*-- +* compo files: +*-- +SEQ_ASCII Multicompo ; + + +*---------------------------------------------------------- +* +* FULL-CORE CALCULATION +* ====================== +***** +* Geometry construction: +***** + +GEOM := Pgeohex ; + +***** +* Reactor material index: +***** + +GEOM MATEX := USPLIT: GEOM :: EDIT 0 NGRP 2 MAXR 728 + NREFL 1 RMIX 3 + NFUEL 2 FMIX 1 2 ; + +****** +* Fuel Map construction +****** + +FLMAP MATEX := RESINI: MATEX :: + EDIT 0 + ::: GEO: HEXZ 91 8 + Z- REFL Z+ REFL + HBC COMPLETE REFL + SIDE 1.0854185 + MIX +*---- +* 1st floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +*---- +*2nd floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +*---- +*3rd floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +*---- +*4th floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +*---- +*5th floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +*---- +*6th floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +*---- +*7th floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +*---- +*8th floor +*---- + +* Ring 0 (center cell) +2 +* Ring 1 +1 1 1 1 1 1 +* Ring 2 +1 2 1 2 1 2 +1 2 1 2 1 2 +* Ring 3 +2 1 1 2 1 1 +2 1 1 2 1 1 +2 1 1 2 1 1 +* Ring 4 +1 1 2 1 1 1 +2 1 1 1 2 1 +1 1 2 1 1 1 +2 1 1 1 2 1 +* Ring 5 : reflector, not in fuel-map. +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 +0 0 0 0 0 0 + +MESHZ 0.0 0.8 1.6 2.4 3.2 4.0 4.8 5.6 6.4 +; + +NCOMB ALL + +EDIT 0 +BTYPE INST-BURN +INST-BVAL BUND + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 +*10 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 +*20 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 +*30 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 +*40 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 +*50 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 +*60 + 60.0 60.0 60.0 60.0 60.0 60.0 60.0 60.0 +*61 = ncomb + +ADD-PARAM PNAME Temp PARKEY Temp GLOBAL +SET-PARAM Temp 740.0 ; + +***** +* Macrolib for the fuel +***** + +*---- +* Mix number assignments in Multicompo +* 1 = fuel material (TRISO) +* 2 = helium +* 3 = Fe56 (reflector) +* Temp in Multicompo : 700.0 and 750.0 K +*---- + +CPO1 := Multicompo ; + +*--- +*Fuel-map macrolib +*--- +MACFL := NCR: CPO1 FLMAP :: + EDIT 0 MACRO CUBIC + TABLE CPO1 default Combustion + MIX 1 FROM 00000001 + INST-BURN + ENDMIX + MIX 2 FROM 00000002 + INST-BURN + ENDMIX + ; +*--- +*Reflector macrolib +*--- +MACRO1 := NCR: CPO1 :: + EDIT 0 MACRO CUBIC + NMIX 3 COMPO CPO1 default + MIX 3 FROM 00000003 + SET 'Combustion' 0.0 + SET 'Temp' 700.0 + ENDMIX + ; + +*--- +* Extended macrolib (=complete because no device) +*--- +MACRO2 MATEX := MACINI: MATEX MACRO1 MACFL :: EDIT 0 ; + +MACFL := DELETE: MACFL ; + +***** +*Tracking, finite-differences method +***** + +TRACK := TRIVAT: GEOM :: EDIT 0 + MAXR 5000 MCFD 1 ; + +***** +* Numerical solution +***** + +SYSTEM := TRIVAA: MACRO2 TRACK :: EDIT 0 ; + +MACRO2 := DELETE: MACRO2 ; + +FLUX := FLUD: SYSTEM TRACK :: EDIT 0 + ACCE 3 3 EXTE 1000 1E-04 ADI 6 ; +assertS FLUX :: "K-EFFECTIVE" 1 1.119211 ; + +POWER := FLPOW: FLMAP FLUX TRACK MATEX :: + EDIT 0 PTOT 600.0 ; + +GREP: POWER :: GETVAL "PMAX-BUND" 1 >>mbp<< ; +assertS POWER :: "PMAX-BUND" 1 1.782570E+03 ; + +GREP: POWER :: GETVAL "PMAX-CHAN" 1 >>mcp<< ; +assertS POWER :: "PMAX-CHAN" 1 1.426045E+04 ; + +GREP: POWER :: GETVAL "K-EFFECTIVE" 1 >>Keff<< ; +assertS POWER :: "K-EFFECTIVE" 1 1.119211 ; + +END: ; +QUIT "LIST" . |