Initial commit from Polytechnique Montreal

1 files changed, 196 insertions, 0 deletions

diff --git a/Dragon/data/TEST_GEO_hex_sect_tspc.x2m b/Dragon/data/TEST_GEO_hex_sect_tspc.x2m
new file mode 100644
index 0000000..60c56ad
--- /dev/null
+++ b/Dragon/data/TEST_GEO_hex_sect_tspc.x2m
@@ -0,0 +1,196 @@
+*-------------------------------------------------------------------------------------------------
+*  Name          : TEST_GEO_hex_sect_tspc.x2m                                                                
+*  Author        : G. Cirillo
+*  HEXAGONAL PIN MODEL TO TEST GLOW GEO
+*-------------------------------------------------------------------------------------------------
+*  Define STRUCTURES and MODULES used
+*-------------------------------------------------------------------------------------------------
+SEQ_ASCII TDT_LAT :: FILE './ALFRED_1cell.dat' ;
+SEQ_ASCII TDT_LAT_SECT :: FILE './ALFRED_1cell_sect.dat' ;
+SEQ_ASCII HexTEST.ps :: FILE './HexTEST.ps' ;
+SEQ_ASCII HexTEST2.ps :: FILE './HexTEST2.ps' ;
+SEQ_ASCII Lines_HEX.m :: FILE './Lines_HEX.m' ;
+LINKED_LIST 
+    MICRO MICRO2 TRACK TRACK2 CP FLUX OUT ;
+MODULE 
+    GEO: G2S: LIB: SALT: MCCGT: TLM: TONE: ASM: FLU: EDI: END: ;
+SEQ_BINARY
+    TRACK_f TRACK_f2 ;
+PROCEDURE assertS ;
+*
+*-------------------------------------------------------------------------------------------------
+*  Input Data 
+*-------------------------------------------------------------------------------------------------
+REAL temp_fuel := 1200.0 ;    (*in K*)
+REAL temp_clad := 600.0 ;     (*in K*)
+REAL temp_cool := 600.0 ;     (*in K*)
+REAL densUO2 := 10.3235 ;     (*in g/cm^3*)
+REAL densPuO2 := 10.7771 ;    (*in g/cm^3*)
+(*REAL densMOX := 10.41561 ;    in g/cm^3*)
+REAL densAIM1 := 7.972 ;      (*in g/cm^3*)
+REAL densLEAD := 10.503 ;     (*in g/cm^3*)
+REAL densHe := 0.0001663 ;    (*in g/cm^3*)
+
+*-------------------------------------------------------------------------------------------------
+*  Microlib: DRAGLIB ENDFB_VII 315 energy discr
+*-------------------------------------------------------------------------------------------------
+MICRO := LIB: ::
+    EDIT 2
+    NMIX 8         (*MAXIMUM OF MATERIAL MIXTURES*)
+    CTRA NONE
+    ANIS 2         (*LINEAR ANISOTROPY SCATTERING*)
+    MIXS LIB: DRAGON FIL: DLIB_315
+        MIX 1 <<temp_fuel>> <<densUO2>>
+            U235   = U235   1.762963E-03      1              O16    = O16    1.18235352E-01      1
+            U238   = U238   8.79681435E-01    1              (*O17    = O17    4.77255E-05*)
+            U234   = U234   5.81718E-06       1              (*O18    = O18    2.66717E-04*)
+            O17    = O17    3.144425E-04
+        MIX 2 <<temp_fuel>> <<densPuO2>>
+            Pu238  = Pu238  2.030129548E-02   1              Pu241  = Pu241  5.313709197E-02     1
+            Pu239  = Pu239  4.951133413E-01   1              Pu242  = Pu242  6.696932147E-02     1
+            Pu240  = Pu240  2.3501943E-01     1              Am241  = Am241  1.16169418E-02      1
+            O16    = O16    1.175300115E-01   1
+            O17    = O17    3.1256696E-04
+            (*O17    = O17    4.74407E-05*)
+            (*O18    = O18    2.6512626E-04*)
+        MIX 3 COMB 1 0.80204511 2 0.19795489
+        MIX 4 COMB 3 1.0
+        MIX 5 COMB 3 1.0
+        MIX 6 <<temp_clad>> <<densHe>> GAS
+            He4    = He4   1.0
+        MIX 7 <<temp_clad>> <<densAIM1>>
+            B10    = B10   1.1013297525043955E-05   2         Ni58   = Ni58   0.10415644320971433     2 
+            B11    = B11   4.898670247495604E-05    2         Ni60   = Ni60   0.04150274249380582     2
+            C0     = C0    9E-04                    2
+        (*C12    = C12   0.0008892061605771179*)              Ni61   = Ni61   0.0018342065822099006   2
+        (*C13    = C13   1.0793839422882127E-05*)             Ni62   = Ni62   0.005943974570405113    2
+            N14    = N14   0.0001494115819963943    2         Ni64   = Ni64   0.0015626331438648383   2
+            N15    = N15   5.884180036056514E-07    2         Cu63   = Cu63   0.00030815637842357554
+            O16    = O16   9.973468728094326E-05    2         Cu65   = Cu65   0.00014184362157642448
+            O17    = O17   2.653127191E-07
+        (*O17    = O17   4.02684437425691E-08*)               Zr90   = Zr90   0.00015211835252769947  2
+        (*O18    = O18   2.2504427531418012E-07*)             Zr91   = Zr91   3.3542662359302346E-05  2
+            Al27   = Al27  0.00015                  2         Zr92   = Zr92   5.183430439088447E-05   2
+            Si28   = Si28  0.007809248969041125     2         Zr94   = Zr94   5.367331302703024E-05   2
+            Si29   = Si29  0.0004106987748868437    2         Zr96   = Zr96   8.831367695083468E-06   2
+            Si30   = Si30  0.0002800522560720316    2         Nb93   = Nb93   0.00015                 2
+            P31    = P31   0.00045                            Mo92   = Mo92   0.0021047837070591823   2
+            S32    = S32   9.476861080827586E-05              Mo94   = Mo94   0.001348697757094504    2
+            S33    = S33   7.698819377540353E-07              Mo95   = Mo95   0.0023550696583545722   2
+            S34    = S34   4.445152400876023E-06              Mo96   = Mo96   0.00249980045270409     2
+            S36    = S36   1.6354853094087033E-08             Mo97   = Mo97   0.0014516393261295652   2
+            Ca40   = Ca40  0.00019332358625973815             Mo98   = Mo98   0.00371810585345211     2
+            Ca42   = Ca42  1.3547188215775078E-06   2         Mo100  = Mo100  0.001521903245205976    2
+            Ca43   = Ca43  2.8940716633277043E-07   2         Sn112  = Sn112  4.571964310655134E-07
+            Ca44   = Ca44  4.575631664190023E-06    2         Sn114  = Sn114  3.1663618573474405E-07
+            Ca46   = Ca46  9.172845337918173E-09              Sn115  = Sn115  1.6454847928465057E-07
+            Ca48   = Ca48  4.474832428236633E-07    2         Sn116  = Sn116  7.098010841901756E-06
+            Ti46   = Ti46  0.0003168038086310422    2         Sn117  = Sn117  3.7815425451416758E-06
+            Ti47   = Ti47  0.00029191129292780875   2         Sn118  = Sn118  1.2027521650336542E-05
+            Ti48   = Ti48  0.0029538020902329385    2         Sn119  = Sn119  4.3019900178835224E-06
+            Ti49   = Ti49  0.00022128762152676788   2         Sn120  = Sn120  1.6453584513264257E-05
+            Ti50   = Ti50  0.00021619518668144252   2         Sn122  = Sn122  2.377274090454982E-06
+            V50    = V50   7.353610074676385E-07              Sn124  = Sn124  3.021695244932352E-06
+            V51    = V51   0.0002992646389925323              Sb121  = Sb121  2.840390648290279E-05
+            Cr50   = Cr50  0.006051845602941554     2         Sb123  = Sb123  2.159609351709721E-05
+            Cr52   = Cr52  0.12136407540945154      2     (*Ta180  = Ta180 1.7915400047357467E-08*)
+            Cr53   = Cr53  0.014026700887079016     2     (*Ta181  = Ta181 0.00014998208459995264*)
+                                                          (*Ta181  = Ta181  1.5E-04*)
+            Cr54   = Cr54  0.00355737810052789      2     (*W180   = W180  3.5237265941784446E-07*)
+            Mn55   = Mn55  0.015                    2     (*W182   = W182  7.868114842886196E-05*)
+                                                              W182   = W182   7.903352109E-05
+            Fe54   = Fe54  0.03688469009223297                W183   = W183   4.2721807585340764E-05
+            Fe56   = Fe56  0.6004294481283934       2         W184   = W184   9.19745762621493E-05
+            Fe57   = Fe57  0.01411455230058187                W186   = W186   8.627009506423013E-05
+            Fe58   = Fe58  0.0019113094787918548    2
+            Co59   = Co59  0.0003
+        MIX 8 <<temp_cool>> <<densLEAD>>
+            Na23   = Na23  1E-05                              As75   = As75   5E-06
+            Mg24   = Mg24  7.7909641699443E-06                Ag107  = Ag107  5.137617589720804E-06
+            Mg25   = Mg25  1.0300414880798231E-06             Ag109  = Ag109  4.862382410279197E-06
+            Mg26   = Mg26  1.1789943419758772E-06             Sn112  = Sn112  4.571964310655134E-08
+            Ca40   = Ca40  9.666179312986907E-08              Sn114  = Sn114  3.166361857347441E-08
+            Ca42   = Ca42  6.773594107887539E-08    3         Sn115  = Sn115  1.645484792846506E-08
+            Ca43   = Ca43  1.4470358316638523E-08   3         Sn116  = Sn116  7.098010841901756E-07
+            Ca44   = Ca44  2.2878158320950115E-07   3         Sn117  = Sn117  3.781542545141676E-07
+            Ca46   = Ca46  4.586422668959087E-10              Sn118  = Sn118  1.2027521650336543E-06
+            Ca48   = Ca48  2.2374162141183165E-08   3         Sn119  = Sn119  4.3019900178835226E-07
+            Fe54   = Fe54  5.645558222706856E-07              Sn120  = Sn120  1.6453584513264259E-06
+            Fe56   = Fe56  9.190152877956246E-06    3         Sn122  = Sn122  2.3772740904549819E-07
+            Fe57   = Fe57  2.1603686136746366E-07             Sn124  = Sn124  3.021695244932352E-07
+            Fe58   = Fe58  2.9254438405605885E-08   3         Sb121  = Sb121  5.680781296580559E-06
+            Cu63   = Cu63  6.8479195205239015E-06             Sb123  = Sb123  4.319218703419443E-06       
+            Cu65   = Cu65  3.1520804794761E-06                Pb204  = Pb204  0.01377877075349251     3
+            (*Zn64   = Zn64  4.8080491426635065E-06*)         Pb206  = Pb206  0.2395190658912157      3
+            Zn66   = Zn66  2.796254092194871E-06              Pb207  = Pb207  0.22070984690476395     3
+            Zn67   = Zn67  4.135741306808834E-07              Pb208  = Pb208  0.5258423164505278      3
+            Zn68   = Zn68  1.916879606631653E-06              Bi209  = Bi209  6E-05
+            Zn70   = Zn70  6.524302782908711E-08                                          
+        ;
+
+*-------------------------------------------------------------------------------------------------
+*  Tracking + PostScript
+*-------------------------------------------------------------------------------------------------
+*---------------------------------------------------------
+*  PostScript
+*---------------------------------------------------------
+HexTEST.ps := G2S: TDT_LAT_SECT :: 
+    EDIT 2 
+    DRAWMIX ;
+
+HexTEST2.ps := G2S: TDT_LAT :: 
+    EDIT 2 
+    DRAWMIX ;
+*---------------------------------------------------------
+*  Tracking sectorized geometry for self-shielding calc
+*---------------------------------------------------------
+TRACK2 TRACK_f2 := SALT: TDT_LAT ::
+    TITL 'Tracking for self-shielding'
+    EDIT 2
+    ANIS 1
+    TSPC 12 50.0
+    ;
+
+*-------------------------------------------------------------------------------------------------
+*  Resonance self-shielding calculation 
+*-------------------------------------------------------------------------------------------------
+MICRO2 := TONE: MICRO TRACK2 TRACK_f2 :: EDIT 1 MXIT 3 TRAN PIJ ;
+
+*---------------------------------------------------------
+*  Tracking sectorized geometry for the flux calc
+*---------------------------------------------------------
+TRACK TRACK_f := SALT: TDT_LAT_SECT ::
+    TITL 'Tracking for flux calculation'
+    EDIT 2
+    ANIS 2
+    TSPC 12 50.0 LONG
+    ;
+
+TRACK := MCCGT: TRACK TRACK_f ::
+    EDIT 2 
+    GAUS 4 AAC 80 TMT STIS 1 EPSI 1E-5 MCU 10000
+    ;
+
+Lines_HEX.m := TLM: TRACK TRACK_f ::
+  EDIT 99
+  NTPO 1
+  DIRECTIONS NoPause DIR 0 ;
+*-------------------------------------------------------------------------------------------------
+*  Flux calculation 
+*-------------------------------------------------------------------------------------------------
+CP := ASM: MICRO2 TRACK TRACK_f :: ARM ;
+FLUX := FLU: CP MICRO2 TRACK TRACK_f :: TYPE K P1 ;
+
+*-------------------------------------------------------------------------------------------------
+*  Editing 
+*-------------------------------------------------------------------------------------------------
+OUT := EDI: MICRO2 TRACK FLUX TDT_LAT_SECT CP ::
+    EDIT 4 
+    MERG COMP 
+    COND
+    ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.341990 ;
+
+ECHO "test TEST_GEO_hex_sect_tspc completed" ;
+END: ;
+QUIT "LIST" .