diff options
Diffstat (limited to 'Donjon/data/testNAPhet.x2m')
| -rw-r--r-- | Donjon/data/testNAPhet.x2m | 541 |
1 files changed, 541 insertions, 0 deletions
diff --git a/Donjon/data/testNAPhet.x2m b/Donjon/data/testNAPhet.x2m new file mode 100644 index 0000000..20d3d47 --- /dev/null +++ b/Donjon/data/testNAPhet.x2m @@ -0,0 +1,541 @@ +*====================================================================== +* Name : testNAPhet.x2m +* Type : DONJON File +* +* Use : Input file to perform PPR +* Notes: +* 1) Homogenization can be homogeneous or heterogeneous +* 2) The pin-by-pin geometry is also created to be able to visualize it +*---------------------------------------------------------------------- +* Donjon: v5.0.1_ev346 modified (or more recent) +*---------------------------------------------------------------------- +* +* Author : R. Chambon (06/2014) +* +* modified : R. Chambon (09/2015) +* Note: The '_ACPOenrichU0' enrich multicompo can be obtained by running +* ----- twice the 'rep900EnrichCOMPOhet.x2m' input file once for +* homogeneous option and once for heterogeneous option. +* +SEQ_ASCII ACOMPOU :: FILE '_ACPOenrichU0' ; +SEQ_ASCII ACOMPOR :: FILE '_ACompo' ; +* +SEQ_ASCII AGEO1 :: FILE 'AGEO1' ; +SEQ_ASCII AGEOH AmapH AMatH :: FILE 'AGEOH' 'AmapH' 'AMatH' ; +SEQ_ASCII AMFH AMACROH AMACR :: FILE 'AMFH' 'AMACROH' 'AMACR' ; +SEQ_ASCII AFLU ATRK :: FILE 'AFLU_NAP' 'ATRK_NAP' ; +LINKED_LIST CpoU CpoR DONNEES Geo1 + GeoH MatexH FmapH MacroFH Macro2H + GeoP + CpoRefl Micro Macro + Track System Flux + FmapH2 FmapH3 PowerIter +; +PROCEDURE assertS ; +* +MODULE GEO: USPLIT: RESINI: MACINI: NCR: NAP: + TRIVAA: TRIVAT: FLUD: FLPOW: + END: DELETE: GREP: UTL: ABORT: ; +* +STRING Fuel := "UOX" ; ! UOX + +INTEGER npin maxr ; +REAL BURN1 BURN2 BURN3 := 0. 0. 0. ; +REAL CBOR := 600. ; +REAL Cote Lame mx1 mx2 mx3 mx4 mx5 mx6 Size wat ; +REAL Ptot := 9.0 ; ! MW +EVALUATE npin := 17 ; +* +EVALUATE Cote := 1.26 ; +EVALUATE Lame := 0.04 ; +EVALUATE wat := 10. ; +EVALUATE Size := npin I_TO_R Cote * Lame 2. * + ; +EVALUATE mx1 := Size -1.5 * wat - ; +EVALUATE mx2 := Size -1.5 * ; +EVALUATE mx3 := Size -0.5 * ; +EVALUATE mx4 := Size 0.5 * ; +EVALUATE mx5 := Size 1.5 * ; +EVALUATE mx6 := Size 1.5 * wat + ; + +* +STRING TypeCompo := "Homog" ; ! Homog : homogeneous properties + ! Heter : heterogeneous (3 mix) properties + ! PinByPin : PinByPin properties + +STRING DirAss := "Ass " ; ! the name used in the multicompo + ! Homog Ass + ! PinByPin Pin + ! Heter Het1 + +INTEGER ifx := 0 ; ! ifx value used in the enrich multicompo + ! Homog 0 + ! Heter 1 + + +STRING DirR := "EDI2B " ; +STRING AutoAss := "Manual" ; ! Definition of assembly in FmapH + ! "Auto" done by NAP + ! "Manual" by the user directly in the input + ! for homogeneous properties only +* +CpoU := ACOMPOU ; +CpoR := ACOMPOR ; + +!-------------------------- +!-------------------------- +* Test NAP: Geometry part +!-------------------------- +! Build geometry for homogeneous assembly +Geo1 := GEO: :: CAR3D 5 5 4 + EDIT 0 + X- VOID X+ VOID + Y- VOID Y+ VOID + Z- VOID Z+ VOID + + MIX +PLANE 1 + 0 4 4 4 0 + 4 4 4 4 4 + 4 4 4 4 4 + 4 4 4 4 4 + 0 4 4 4 0 +PLANE 2 + 0 4 4 4 0 + 4 3 2 3 4 + 4 2 1 2 4 + 4 3 2 3 4 + 0 4 4 4 0 +PLANE 3 SAME 2 +PLANE 4 SAME 1 + MESHX <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHY <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHZ 0. 10. 110. 210. 220. + SPLITX 2 1 1 1 2 + SPLITY 2 1 1 1 2 + SPLITZ 2 1 1 2 +; + +IF TypeCompo "Homog" = THEN + ECHO " Build geometry for homogeneous assembly " ; + GeoH := NAP: Geo1 CpoU :: + EDIT 100 + DIRGEO <<DirAss>> + MIXASS 3 1 2 3 + SPLITX-ASS 2 + SPLITY-ASS 2 + ; + ! Max reg=(2+2*3+2)^2*6=600 + EVALUATE maxr := 600 ; +ELSEIF TypeCompo "Heter" = THEN + ECHO " Build geometry for heterogeneous assembly " ; + GeoH := NAP: Geo1 CpoU :: + EDIT 100 + DIRGEO <<DirAss>> + MIXASS 3 1 2 3 + SPLITX-ASS 1 2 1 + SPLITY-ASS 1 2 1 + ; + ! Max reg=(2+4*3+2)^2*6=1536 + EVALUATE maxr := 1536 ; +ELSEIF TypeCompo "PinByPin" = THEN + ECHO " Build geometry for pin-by-pin assembly " ; + GeoH := NAP: Geo1 CpoU :: + EDIT 100 + DIRGEO <<DirAss>> + MIXASS 3 1 2 3 + ; + ! Max reg=(2+17*3+2)^2*6=18150 + EVALUATE maxr := 18150 ; +ENDIF ; + +AGEO1 := Geo1 ; +AGEOH := GeoH ; +!AGEOP := GeoP ; + +* Max reg=(2+4*3+2)^2*6=1536 +IF AutoAss "Auto" = THEN + GeoH MatexH := USPLIT: GeoH :: + NGRP 2 MAXR 1536 + NREFL 1 RMIX 4 + NFUEL ASBLY + ; +ELSEIF AutoAss "Manual" = THEN + GeoH MatexH := USPLIT: GeoH :: + NGRP 2 MAXR 1536 + NREFL 1 RMIX 4 + NFUEL 6 FMIX 1 5 ! Homog + 2 6 + 3 7 + ; +ENDIF ; +Track := TRIVAT: GeoH :: + MAXR <<maxr>> DUAL 3 3 ; +ATRK := Track ; + +!-------------------------- +* Test NAP: Geometry part embedded in RESINI: +!-------------------------- +IF TypeCompo "Homog" = AutoAss "Auto" = * THEN +ECHO " Build map for homogeneous assembly (auto)" ; + FmapH MatexH := RESINI: MatexH CpoU :: + EDIT 0 + ::: + SPLIT-NAP: + GEO: + CAR3D 5 5 4 + EDIT 10 + X- VOID X+ VOID + Y- VOID Y+ VOID + Z- VOID Z+ VOID + MIX + PLANE 1 + 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 + PLANE 2 + 0 0 0 0 0 + 0 3 2 3 0 + 0 2 1 2 0 + 0 3 2 3 0 + 0 0 0 0 0 + PLANE 3 SAME 2 + PLANE 4 SAME 1 + MESHX <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHY <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHZ 0. 10. 110. 210. 220. + ; + ::: NAP: + DIRGEO <<DirAss>> + MIXASS 3 1 2 3 + MAX-MIX-GEO 4 + ; + ASSEMBLY 9 3 3 + A-ZONE + ASBLY + AXNAME '1' '2' '3' + AYNAME 'A' 'B' 'C' + NXNAME '-' '1' '2' '3' '-' + NYNAME '-' 'A' 'B' 'C' '-' + NCOMB ASBLY + ADD-PARAM PNAME 'ppmBore' PARKEY 'ppmBore' GLOBAL + ; +ELSEIF TypeCompo "Homog" = AutoAss "Manual" = * THEN +ECHO " Build map for homogeneous assembly (manual)" ; + FmapH MatexH := RESINI: MatexH :: + EDIT 0 + ::: + GEO: + CAR3D 5 5 4 + EDIT 10 + X- VOID X+ VOID + Y- VOID Y+ VOID + Z- VOID Z+ VOID + MIX + PLANE 1 + 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 + PLANE 2 + 0 0 0 0 0 + 0 3 2 3 0 + 0 2 1 2 0 + 0 3 2 3 0 + 0 0 0 0 0 + PLANE 3 SAME 2 + PLANE 4 SAME 1 + MESHX <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHY <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHZ 0. 10. 110. 210. 220. + ; + ASSEMBLY 9 3 3 + A-ZONE + 1 2 3 + 4 5 6 + 7 8 9 + A-NX 3 3 3 + A-IBX 1 1 1 +! !--- + AXNAME '1' '2' '3' + AYNAME 'A' 'B' 'C' + NXNAME '-' '1' '2' '3' '-' + NYNAME '-' 'A' 'B' 'C' '-' + NCOMB ASBLY + ADD-PARAM PNAME 'ppmBore' PARKEY 'ppmBore' GLOBAL + ; +ELSEIF TypeCompo "Heter" = THEN + ECHO " Build map for heterogeneous assembly " ; + FmapH MatexH := RESINI: MatexH CpoU :: + EDIT 1000 + ::: + SPLIT-NAP: + GEO: + CAR3D 5 5 4 + EDIT 10 + X- VOID X+ VOID + Y- VOID Y+ VOID + Z- VOID Z+ VOID + MIX + PLANE 1 + 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 + PLANE 2 + 0 0 0 0 0 + 0 3 2 3 0 + 0 2 1 2 0 + 0 3 2 3 0 + 0 0 0 0 0 + PLANE 3 SAME 2 + PLANE 4 SAME 1 + MESHX <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHY <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHZ 0. 10. 110. 210. 220. + ; + ::: NAP: + DIRGEO <<DirAss>> + MIXASS 3 1 2 3 + MAX-MIX-GEO 4 + ; + ASSEMBLY 9 3 3 + A-ZONE + ASBLY + AXNAME '1' '2' '3' + AYNAME 'A' 'B' 'C' + NXNAME '-' '1' '1' '1' '2' '2' '2' '3' '3' '3' '-' + NYNAME '-' 'A' 'A' 'A' 'B' 'B' 'B' 'C' 'C' 'C' '-' + NCOMB ASBLY + ADD-PARAM PNAME 'ppmBore' PARKEY 'ppmBore' GLOBAL +; +ELSEIF TypeCompo "PinByPin" = THEN + ECHO " Build map for heterogeneous assembly " ; + FmapH MatexH := RESINI: MatexH CpoU :: + EDIT 0 + ::: + SPLIT-NAP: + GEO: + CAR3D 5 5 4 + EDIT 10 + X- VOID X+ VOID + Y- VOID Y+ VOID + Z- VOID Z+ VOID + MIX + PLANE 1 + 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 + PLANE 2 + 0 0 0 0 0 + 0 3 2 3 0 + 0 2 1 2 0 + 0 3 2 3 0 + 0 0 0 0 0 + PLANE 3 SAME 2 + PLANE 4 SAME 1 + MESHX <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHY <<mx1>> <<mx2>> <<mx3>> <<mx4>> <<mx5>> <<mx6>> + MESHZ 0. 10. 110. 210. 220. + ; + ::: NAP: + DIRGEO <<DirAss>> + MIXASS 3 1 2 3 + MAX-MIX-GEO 4 + ; + ASSEMBLY 9 3 3 + A-ZONE + ASBLY + AXNAME '1' '2' '3' + AYNAME 'A' 'B' 'C' + NXNAME '-' '1' '1' '1' '1' '1' '1' '1' '1' '1' + '1' '1' '1' '1' '1' '1' '1' '1' + '2' '2' '2' '2' '2' '2' '2' '2' '2' + '2' '2' '2' '2' '2' '2' '2' '2' + '3' '3' '3' '3' '3' '3' '3' '3' '3' + '3' '3' '3' '3' '3' '3' '3' '3' + '-' + NYNAME '-' 'A' 'A' 'A' 'A' 'A' 'A' 'A' 'A' 'A' + 'A' 'A' 'A' 'A' 'A' 'A' 'A' 'A' + 'B' 'B' 'B' 'B' 'B' 'B' 'B' 'B' 'B' + 'B' 'B' 'B' 'B' 'B' 'B' 'B' 'B' + 'C' 'C' 'C' 'C' 'C' 'C' 'C' 'C' 'C' + 'C' 'C' 'C' 'C' 'C' 'C' 'C' 'C' + '-' + NCOMB ASBLY + ADD-PARAM PNAME 'ppmBore' PARKEY 'ppmBore' GLOBAL + ; +ENDIF ; + +FmapH := RESINI: FmapH :: EDIT 2 + BTYPE INST-BURN + INST-BVAL ASBL + <<BURN3>> <<BURN2>> <<BURN3>> + <<BURN2>> <<BURN1>> <<BURN2>> + <<BURN3>> <<BURN2>> <<BURN3>> + + SET-PARAM 'ppmBore' 0.0 +; +AMatH := MatexH ; + +!-------------------------- +* Test NCR: assembly related keywords +!-------------------------- +! Compute cross-section for heterogeneous assembly +IF AutoAss "Auto" = THEN + MacroFH := NCR: CpoU FmapH :: + EDIT 0 + MACRO + LINEAR + TABLE CpoU <<DirAss>> 'burnup' + MIX 1 ASBLY INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + MIX 2 ASBLY INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + MIX 3 ASBLY INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + ; +ELSEIF AutoAss "Manual" = THEN ! for homogene only + MacroFH := NCR: CpoU FmapH :: + EDIT 100 + MACRO + LINEAR + TABLE CpoU <<DirAss>> 'burnup' +!* look back at mixture renumbering earlier, USPLIT: module +!* homogene + MIX 1 FROM 1 INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + MIX 2 FROM 1 INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + MIX 3 FROM 1 INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + MIX 5 FROM 1 INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + MIX 6 FROM 1 INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + MIX 7 FROM 1 INST-BURN + SET 'burnup' 0.0 + SET 'ppmBore' <<CBOR>> + SET 'TF' 800.0 + SET 'DCA' 0.659 + SET 'TCA' 600.0 + SET 'BARR' 0.0 + ENDMIX + ; +ENDIF ; +AMFH := MacroFH ; + +Micro := NCR: CpoR :: +* dummy reflector -- please use more realistic data + EDIT 0 + LINEAR NMIX 4 + COMPO CpoR <<DirR>> + MIX 4 SET 'burnup' 0.0 + SET 'ppmBore' 0.0 + ENDMIX +; +Macro := Micro :: STEP UP MACROLIB ; +AMACR := Macro ; +Micro := DELETE: Micro ; + +Macro2H MatexH := MACINI: MatexH Macro MacroFH :: +EDIT 0 +; +AMACROH := Macro2H ; + + +* Compute flux distribution +System := TRIVAA: Macro2H Track ; + +Flux := FLUD: System Track :: + EDIT 0 + ADI 4 + EXTE 500 0.00001 +; +AFLU := Flux ; + +REAL keff2 ; +GREP: Flux :: + GETVAL 'K-EFFECTIVE ' 1 >>keff2<< ; +ECHO "CBOR=" CBOR "K-effective 2 = " keff2 ; + +PowerIter FmapH := FLPOW: FmapH Flux Track MatexH :: + EDIT 5 PTOT <<Ptot>> PRINT DISTR POWER ; + +* Perform PPR +IF TypeCompo "Homog" = TypeCompo "Heter" = + THEN + FmapH := NAP: FmapH Track Flux MatexH MacroFH :: + EDIT 4 + PPR + NZASS 2 + METH GPPR <<ifx>> + ; +ENDIF ; + +AmapH := FmapH ; + +FmapH2 := FmapH :: + STEP UP ASSEMBLY + STEP AT 1 +; +FmapH2 := UTL: FmapH2 :: + DIR +; +assertS FmapH2 :: 'ASS-POWER' 1 6.361761E02 ; + +ECHO "test testNAPhet.x2m completed" ; +END: ; +QUIT . |