Initial commit from Polytechnique Montreal

1 files changed, 208 insertions, 0 deletions

diff --git a/Dragon/data/tdraglib_proc/TCDR22.c2m b/Dragon/data/tdraglib_proc/TCDR22.c2m
new file mode 100644
index 0000000..8983c9f
--- /dev/null
+++ b/Dragon/data/tdraglib_proc/TCDR22.c2m
@@ -0,0 +1,208 @@
+*----
+*  TEST CASE TRDR22
+*  17 X 17 PWR VOIDED ASSEMBLY
+*  172-GROUP JEF2.2 DRAGLIB
+*  Author        : Alain Hebert (2019)
+*
+*  REF: none
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+  ASSMB DISCR LIBRARY CP CALC OUT SPHGEOM MTRACK MACRO MACRO2 SYS FLUX
+  EDIT ;
+MODULE
+  LIB: GEO: SYBILT: SHI: ASM: FLU: EDI: FLUD: MAC: DELETE: UTL: BIVACT:
+  BIVACA: SPH: OUT: ERROR: END: ;
+PROCEDURE assertS ;
+*----
+*  Microscopic cross sections from Jef 2.2 Draglib
+*----
+LIBRARY := LIB: ::
+  NMIX 8 CTRA WIMS ANIS 2
+  MIXS LIB: DRAGON FIL: DLIB_J2
+  MIX 1 579.9
+    H1H2O    = H1_H2O   4.76690E-2   O16H2O   = O16      2.38345E-2
+    B10      = B10      4.71444E-4   B11      = B11      1.90959E-3
+  MIX 2 579.9                        O16      = O16      3.06711E-4
+    Cr52     = Cr52     7.54987E-5   Fe56     = Fe56     1.47624E-4
+    Zr91     = Zr0      4.18621E-2
+  MIX 3 579.9
+    H1H2O    = H1_H2O   4.65292E-2   O16H2O   = O16      2.32646E-2
+                                     Cr52     = Cr52     4.79927E-5
+    Fe56     = Fe56     4.45845E-5   Ni58     = Ni58     1.13521E-4
+                                     MoNat    = Mo95     4.03755E-6
+                                     Al27     = Al27     2.35231E-6
+    Mn55     = Mn55     4.15901E-7
+    B10      = B10      4.60867E-6   B11      = B11      1.86674E-5
+    Zr91     = Zr0      8.92427E-4
+  MIX 4 579.9                        O16      = O16      2.87335E-4
+    Cr52     = Cr52     7.07291E-5   Fe56     = Fe56     1.38298E-4
+    Zr91     = Zr0      3.92175E-2
+  MIX 5 579.9
+    H1H2O    = H1_H2O   4.71346E-2   O16H2O   = O16      2.35673E-2
+                                     Cr52     = Cr52     2.24991E-5
+    Fe56     = Fe56     2.09013E-5   Ni58     = Ni58     5.32188E-5
+                                     MoNat    = Mo95     1.89281E-6
+                                     Al27     = Al27     1.10277E-6
+    Mn55     = Mn55     1.94976E-7
+    B10      = B10      4.66484E-6   B11      = B11      1.88950E-5
+    Zr91     = Zr0      4.18372E-4
+  MIX 6 579.9
+    H1H2O    = H1_H2O   4.71676E-2   O16H2O   = O16      2.35838E-2
+                                     Cr52     = Cr52     2.11122E-5
+    Fe56     = Fe56     1.96130E-5   Ni58     = Ni58     4.99383E-5
+                                     MoNat    = Mo95     1.77614E-6
+                                     Al27     = Al27     1.03479E-6
+    Mn55     = Mn55     1.82957E-7
+    B10      = B10      4.66791E-6   B11      = B11      1.89074E-5
+    Zr91     = Zr0      3.92583E-4
+  MIX 7 579.9
+    H1H2O    = H1_H2O   4.72020E-2   O16H2O   = O16      2.36010E-2
+                                     Cr52     = Cr52     1.96591E-5
+    Fe56     = Fe56     1.82630E-5   Ni58     = Ni58     4.65011E-5
+                                     MoNat    = Mo95     1.65389E-6
+                                     Al27     = Al27     9.63569E-7
+    Mn55     = Mn55     1.70365E-7
+    B10      = B10      4.67110E-6   B10      = B10      1.89203E-5
+    Zr91     = Zr0      3.65562E-4
+  MIX 8 933.6                        O16      = O16      4.49355E-2
+    U235     = U235     7.39237E-4 1 
+    U238     = U238     2.17285E-2 1 
+  ;
+*----
+*  Geometry ASSMB : a 17 X 17 normal PWR assembly
+*  contains  C1 : cell without fuel
+*            C2 : normal fuel cell
+*            C3 : peripheral cell
+*            C4 : corner cell
+*----
+ASSMB := GEO: :: CAR2D 9 9
+  X- DIAG X+ REFL Y- SYME Y+ DIAG
+  CELL C1 C2 C2 C1 C2 C2 C1 C2 C3
+          C2 C2 C2 C2 C2 C2 C2 C3
+             C2 C2 C2 C2 C2 C2 C3
+                C1 C2 C2 C1 C2 C3
+                   C2 C2 C2 C2 C3
+                      C1 C2 C2 C3
+                         C2 C2 C3
+                            C2 C3
+                               C4
+  MERGE 1  3 12 11 12 12 11 12 15
+           4  6  5  6  6  5  6  8
+             13  5  6  6  5  6  8
+                 2  5  5 10  5  8
+                   13  5  5  6  8
+                       2  5  7  8
+                         13  7  8
+                            14  8
+                                9
+  TURN  H  H  B  H  H  B  H  H  A
+           H  G  G  H  G  G  H  A
+              A  E  E  F  E  E  A
+                 H  H  F  H  H  A
+                    H  E  G  H  A
+                       H  H  A  A
+                          H  A  A
+                             A  A
+                                A
+  ::: C1 := GEO: CARCEL 2
+    MESHX 0.0 1.26472 MESHY 0.0 1.26472
+    RADIUS 0.0 0.572435 0.613142  MIX 1 2 3 ;
+  ::: C2 := GEO: C1  RADIUS 0.0 0.412660 0.474364  MIX 8 4 5 ;
+  ::: C3 := GEO: C2  MESHX 0.0 1.31472 MIX 8 4 6 ;
+  ::: C4 := GEO: C3  MESHY 0.0 1.31472  MIX 8 4 7 ;
+  ;
+*----
+*  Self-Shielding calculation SYBIL
+*  Transport calculation      SYBIL
+*  Flux calculation for P1 homogeneous leakage
+*----
+DISCR := SYBILT: ASSMB ::
+  TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON'
+  MAXR 400 QUA2 6 3 ;
+LIBRARY := SHI: LIBRARY DISCR :: LJ ;
+ASSMB DISCR := DELETE: ASSMB DISCR ;
+*----
+*  New lattice geometry with VOID boundary condition
+*----
+ASSMB := GEO: :: CAR2D 9 9
+  X- DIAG X+ VOID Y- SYME Y+ DIAG
+  CELL C1 C2 C2 C1 C2 C2 C1 C2 C3
+          C2 C2 C2 C2 C2 C2 C2 C3
+             C2 C2 C2 C2 C2 C2 C3
+                C1 C2 C2 C1 C2 C3
+                   C2 C2 C2 C2 C3
+                      C1 C2 C2 C3
+                         C2 C2 C3
+                            C2 C3
+                               C4
+  ::: C1 := GEO: CARCEL 2
+    MESHX 0.0 1.26472 MESHY 0.0 1.26472
+    RADIUS 0.0 0.572435 0.613142  MIX 1 2 3 ;
+  ::: C2 := GEO: C1  RADIUS 0.0 0.412660 0.474364  MIX 8 4 5 ;
+  ::: C3 := GEO: C2  MESHX 0.0 1.31472 MIX 8 4 6 ;
+  ::: C4 := GEO: C3  MESHY 0.0 1.31472  MIX 8 4 7 ;
+  ;
+DISCR := SYBILT: ASSMB ::
+  TITLE 'TCDR02: 17 X 17 MULTICELL PWR BENCHMARK WITHOUT POISON'
+  MAXR 400 QUA2 6 3 ;
+CP := ASM: LIBRARY DISCR :: ARM ;
+CALC := FLU: CP LIBRARY DISCR ::
+  TYPE K P1 HETE ;
+assertS CALC :: 'K-EFFECTIVE' 1 0.4551677 ;
+
+*----
+*  Editing using SPH model for transport-diffusion with flux-volume
+*  normalization
+*----
+OUT := EDI: LIBRARY DISCR CALC ASSMB ::
+  EDIT 4 UPS SAVE MICR RES MERGE CELL COND 4.0
+  ;
+
+*----
+* TRANSPORT-DIFFUSION SPH EQUIVALENCE IN NON-FUNDAMENTAL MODE
+*----
+SPHGEOM := GEO: :: CAR2D 9 9
+  X- DIAG X+ ALBE 1 Y- SYME Y+ DIAG
+  MIX   1  2  3  4  5  6  7  8  9
+          10 11 12 13 14 15 16 17
+             18 19 20 21 22 23 24
+                25 26 27 28 29 30
+                   31 32 33 34 35
+                      36 37 38 39
+                         40 41 42
+                            43 44
+                               45
+  MESHX 0.0 1.26471996 2.52943993 3.79415989 5.05887985
+        6.32359982 7.58831978 8.85303974 10.1177597 11.4324799
+  MESHY 0.0 1.26471996 2.52943993 3.79415989 5.05887985
+        6.32359982 7.58831978 8.85303974 10.1177597 11.4324799
+;
+MTRACK := BIVACT: SPHGEOM :: DUAL (*IELEM=*) 1 (*ICOL=*) 1 ;
+MACRO := OUT :: STEP UP 'REF-CASE0001' STEP UP 'MACROLIB' ;
+*----
+* Set physical albedos
+*----
+MACRO := MAC: MACRO :: ALBP 1 0.0 0.0 ;
+*----
+* Perform SPH equivalence
+*----
+MACRO2 := SPH: MACRO MTRACK :: EDIT 2 ITER 1.0E-7 ;
+*----
+* Verification calculation
+*----
+SYS := BIVACA: MACRO2 MTRACK ;
+FLUX := FLUD: SYS MTRACK :: EDIT 2 ADI 20 ACCE 5 3 ;
+EDIT := OUT: FLUX SPHGEOM MACRO2 MTRACK ::
+      EDIT 0 INTG MIX ;
+
+ECHO "Error statistics after correction" ;
+ERROR: MACRO EDIT :: HREA 'NFTOT' ;
+
+assertS FLUX :: K-EFFECTIVE 1 0.4551687 ;
+
+ECHO "test TCDR22 completed" ;
+END: ; 
+QUIT "LIST" .