Initial commit from Polytechnique Montreal

1 files changed, 211 insertions, 0 deletions

diff --git a/Dragon/data/twinfrith_proc/TCWW31.c2m b/Dragon/data/twinfrith_proc/TCWW31.c2m
new file mode 100644
index 0000000..e0af5e2
--- /dev/null
+++ b/Dragon/data/twinfrith_proc/TCWW31.c2m
@@ -0,0 +1,211 @@
+*----
+*  TEST CASE TCWW31
+*  CANDU-6 ANNULAR CELL
+*  WIMS-AECL 69 GROUPS LIBRARY FILE WIMSLIB
+*  MULTI-PARAMETER COMPO ACCESS FOR MACRO-DEPLETION
+*  TWO GROUP BURNUP
+*  POWER (KW)         =  615.00000
+*  BURN POWER (KW/KG) =   31.97130
+*  URANIUM MASS       =   19.23600
+*  UO2 REAL DENSITY   =   10.59300
+*  UO2 EFF DENSITY    =   10.43750
+*  UO2 TEMPERATURE    =  941.28998
+*  ENRICHMENT         =    0.71140
+*  COOLANT D2 AT %    =   99.222
+*  MODERATOR D2 AT %  =   99.911
+*  NUMBER OF DAYS     =   50
+*
+*  REF: R. Roy et al.  Ann. Nucl. Energy 21, 115 (1994)
+*
+*----
+*  Define variables
+*     Burnup paremeters
+*     a) Power
+*        =  31.9713 kw/kg for    0.0 to 300.0 days
+*     b) 69 Groups Burnup time interval Delt
+*        =   300 day  for    0 to  300 day
+*     c) 2 Groups Burnup time interval Delt
+*        =     1 day  for    0 to    1 day
+*        =     4 days for    1 to    5 days
+*        =     5 days for    5 to   10 days
+*        =    10 days for   10 to   50 days
+*        =    20 days for   50 to  150 days
+*        =    50 days for  150 to  300 days
+*     c) Days with burnup interval changes
+*        = 1.0, 5.0, 10.0, 50.0, 150.0 and 300.0 days
+*     d) Burnup control time variables Timei, Timef, TotalTime
+*        Timei = initial time
+*        Timef = final time
+*        TotalTime = Final time reached
+*----
+REAL
+  Power    Delt   Timec   Timei Timef TotalTime :=
+  31.9713  1.0    1.0     0.0   0.0   300.0 ;
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+  LIBRARY CANDU6S CANDU6F VOLMATS VOLMATF PIJ FLUX BURNUP EDITION
+  DATABASE ;
+SEQ_BINARY
+  INTLINS INTLINF ;
+SEQ_ASCII
+  res ;
+MODULE
+  GEO: EXCELT: LIB: SHI: ASM: FLU: EVO: EDI: COMPO: SPH: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+*  Depletion data from file WIMSLIB format WIMS-AECL
+*  Microscopic cross sections from file WIMSLIB format WIMS-AECL
+*----
+LIBRARY := LIB: ::
+  NMIX 10 CTRA WIMS
+  DEPL LIB: WIMSAECL FIL: WIMSLIB
+  MIXS LIB: WIMSAECL FIL: WIMSLIB
+  MIX  1  560.66   0.81212           O16      = O16      7.99449E-1
+    D2D2O    = D2D2O    1.99768E-1   H1H2O    = H1H2O    7.83774E-4
+  MIX  2  560.66   6.57              Nb93     = NB93     2.50000
+    BNat     = B-NAT    2.10000E-4
+    Zr91     = ZR91     9.75000E+1
+  MIX  3  345.66   0.0014            He4      = HE4      1.00000E+2
+  MIX  4  345.66   6.44              Fe56     = FE56A    1.60000E-1
+    Ni58     = NI58     6.00000E-2   Cr52     = CR52     1.10000E-1
+    BNat     = B-NAT    3.10000E-4
+    Zr91     = ZR91     9.97100E+1
+  MIX  5  345.66   1.082885          O16      = O16      7.98895E-1
+    D2D2O    = D2D2O    2.01016E-1   H1H2O    = H1H2O    8.96000E-5
+  MIX  6  941.29  10.4375010         O16      = O16      1.18473E+1
+    Xe135    = XE135    0.0
+    U235     = U235     6.27118E-1 1
+    U238     = U238     8.75256E+1 1
+    U236     = U236     0.0        1
+    Pu239    = PU239    0.0        1
+  MIX  7  COMB 6 1.0
+  MIX  8  COMB 6 1.0
+  MIX  9  COMB 6 1.0
+  MIX 10   560.66   6.44             Fe56     = FE56A    1.60000E-1
+    Ni58     = NI58     6.00000E-2   Cr52     = CR52     1.10000E-1
+    BNat     = B-NAT    3.10000E-4
+    Zr91     = ZR91     9.97100E+1
+  ;
+*----
+*  Geometry CANDU6S : 13 regions annular cluster for self-shielding
+*           CANDU6F : 31 regions annular cluster for transport
+*----
+CANDU6S := GEO: ::  TUBE 5
+  R+ REFL RADIUS  0.00000 5.16890 5.60320  6.44780 6.58750 16.12171
+  MIX  1 2 3 4 5
+  CLUSTER ROD1 ROD2 ROD3 ROD4
+  ::: ROD1 := GEO: TUBE 2 MIX 6 10 NPIN  1 RPIN 0.0000 APIN 0.0000
+    RADIUS 0.00000 0.6122 0.6540 ;
+  ::: ROD2 := GEO: ROD1   MIX 7 10 NPIN  6 RPIN 1.4885 APIN 0.0000 ;
+  ::: ROD3 := GEO: ROD1   MIX 8 10 NPIN 12 RPIN 2.8755 APIN 0.261799 ;
+  ::: ROD4 := GEO: ROD1   MIX 9 10 NPIN 18 RPIN 4.3305 APIN 0.0 ;
+  ;
+CANDU6F := GEO: CANDU6S :: SPLITR  6 1 1 1 10
+  ::: ROD1 := GEO: ROD1 SPLITR 2 1 ;
+  ::: ROD2 := GEO: ROD2 SPLITR 2 1 ;
+  ::: ROD3 := GEO: ROD3 SPLITR 2 1 ;
+  ::: ROD4 := GEO: ROD4 SPLITR 2 1 ;
+  ;
+*----
+*  Self-Shielding calculation EXCEL
+*  Transport calculation      EXCEL
+*  Flux calculation for keff
+*----
+VOLMATS INTLINS := EXCELT: CANDU6S ::
+  TITLE 'TCWW31: FEW GROUP BURNUP / SELF-SHIELDING TRACKING'
+  EDIT 0 MAXR 13 TRAK TISO 5 10.0 SYMM 12 ;
+LIBRARY := SHI: LIBRARY VOLMATS INTLINS ::
+  NOLJ EDIT 0 ;
+VOLMATF INTLINF := EXCELT: CANDU6F ::
+  TITLE 'TCWW31: FEW GROUP BURNUP / TRANSPORT TRACKING'
+  EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+  TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.117653 ;
+EDITION := EDI: FLUX LIBRARY VOLMATF CANDU6F ::
+  MERG REGI  6  6 10  7  7 10  1  1  8  8 10  1
+             1  9  9 10  1  1  2  3  4  5  5  5
+             5  5  5  5  5  5  5
+  COND 4.0 MICR ALL SAVE
+  MGEO CANDU6F
+  ;
+EDITION := SPH: EDITION VOLMATF INTLINF ;
+FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S VOLMATF CANDU6F
+INTLINF := DELETE: FLUX PIJ LIBRARY INTLINS VOLMATS CANDU6S
+VOLMATF CANDU6F INTLINF ;
+*----
+*  2 group Burnup (macro-depletion)
+*----
+DATABASE := COMPO: ::
+  EDIT 5
+  COMM  'Multi-parameter reactor database' ENDC
+  INIT
+  ;
+DATABASE := COMPO: DATABASE EDITION ::
+  EDIT 3
+  ;
+EDITION := DELETE: EDITION ;
+*
+* COMPO INTERPOLATION
+LIBRARY := SPH: DATABASE ::
+   EDIT 1
+   STEP UP default STEP AT 1
+   MICRO OFF
+   ;
+CANDU6F := DATABASE :: STEP UP default STEP UP 'GEOMETRIES' STEP AT 1 ;
+VOLMATF INTLINF := EXCELT: CANDU6F ::
+  TITLE 'TCWW31: FEW GROUP BURNUP / TRANSPORT TRACKING'
+  EDIT 0 MAXR 31 TRAK TISO 5 10.0 SYMM 12 ;
+PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+FLUX := FLU: PIJ LIBRARY VOLMATF ::
+  TYPE K ;
+assertS FLUX :: 'K-EFFECTIVE' 1 1.117663 ;
+EDITION := EDI: FLUX LIBRARY VOLMATF ::
+  EDIT 1 MERGE COMP MICR  ALL SAVE ;
+WHILE Timei TotalTime < DO
+  EVALUATE Timef := Timei Delt + ;
+  IF Timei 0.0 = THEN
+    BURNUP LIBRARY := EVO: LIBRARY FLUX VOLMATF ::
+      EDIT 3 DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+  ELSE
+    BURNUP LIBRARY := EVO: BURNUP LIBRARY FLUX VOLMATF ::
+      EDIT 3 NOEX DEPL <<Timei>> <<Timef>> DAY POWR <<Power>> ;
+  ENDIF ;
+  PIJ := DELETE: PIJ ;
+  PIJ := ASM: LIBRARY VOLMATF INTLINF ;
+  FLUX := FLU: FLUX PIJ LIBRARY VOLMATF ::
+    TYPE K ;
+  EDITION := EDI: EDITION FLUX LIBRARY VOLMATF ::
+    SAVE  ;
+*----
+*  change delta t for burnup and final time if required
+*----
+  IF Timef Timec = THEN
+    IF Timec 150.0 = THEN
+      EVALUATE Delt Timec := 50.0 300.0 ;
+    ENDIF ;
+    IF Timec 50.0 = THEN
+      EVALUATE Delt Timec := 20.0 150.0 ;
+    ENDIF ;
+    IF Timec 10.0 = THEN
+      EVALUATE Delt Timec := 10.0 50.0 ;
+    ENDIF ;
+    IF Timec 5.0 = THEN
+      EVALUATE Delt Timec := 5.0 10.0 ;
+    ENDIF ;
+    IF Timec 1.0 = THEN
+      EVALUATE Delt Timec := 4.0 5.0 ;
+    ENDIF ;
+  ENDIF ;
+  EVALUATE Timei := Timef ;
+ENDWHILE ;
+assertS FLUX :: 'K-EFFECTIVE' 1 0.9505190 ;
+res := EDITION ;
+EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F := DELETE:
+  EDITION BURNUP FLUX PIJ LIBRARY INTLINF VOLMATF CANDU6F ;
+ECHO "test TCWW31 completed" ;
+END: ; 
+QUIT "LIST" .