Initial commit from Polytechnique Montreal

1 files changed, 164 insertions, 0 deletions

diff --git a/Trivac/data/NodalTests_proc/prob5p3_nem.c2m b/Trivac/data/NodalTests_proc/prob5p3_nem.c2m
new file mode 100755
index 0000000..f895037
--- /dev/null
+++ b/Trivac/data/NodalTests_proc/prob5p3_nem.c2m
@@ -0,0 +1,164 @@
+*----
+*  TEST CASE prob5p3_nem
+*
+*  REF: A. Hebert, "Applied Reactor Physics", Presses Internationales
+*       Polytechnique, Problem 5.3 (2009).
+*
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST GEOM MACRO TRACK FLUX ;
+MODULE GEO: MAC: NSST: NSSF: DELETE: END: ;
+PROCEDURE assertS ;
+*----
+*  Macroscopic cross sections
+*----
+MACRO := MAC: ::
+ EDIT 2 NGRO 2 NMIX 2 NIFI 1
+ READ INPUT
+ MIX     1
+      DIFF  1.264E+00  0.9328E+00
+     TOTAL  8.154E-03  4.1000E-03
+    NUSIGF  0.000E+00  4.5620E-03
+       CHI  1.000E+00  0.000E+00
+  H-FACTOR  0.000E+00  4.5620E-03
+      SCAT  1 1 0.0 2 2 0.0   7.368E-03
+ MIX     2
+      DIFF  1.310E+00  0.8695E+00
+     TOTAL  1.018E-02  2.1170E-04
+      SCAT  1 1 0.0 2 2 0.0   1.0180E-02
+ ;
+*----
+*  Steady-state calculations with the Nodal Expansion Method
+*----
+GEOM := GEO: :: CAR1D 4
+           X- ZERO X+ ZERO
+           MIX  2 1 1 2
+           MESHX 0.0 40.0 350.0 660.0 700.0
+           SPLITX 1 4 4 1
+           ;
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9977308 ;
+TRACK FLUX := DELETE: TRACK FLUX ;
+*
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 CMFD ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9976518 ;
+GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
+*
+GEOM := GEO: :: CAR1D 4
+           X- VOID X+ VOID
+           MIX  2 1 1 2
+           MESHX 0.0 40.0 350.0 660.0 700.0
+           SPLITX 1 4 4 1
+           ;
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9977973 ;
+TRACK FLUX := DELETE: TRACK FLUX ;
+*
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 CMFD ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9977070 ;
+GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
+*
+GEOM := GEO: :: CAR1D 4
+           X- ALBE 0.85 X+ ALBE 0.95
+           MIX  2 1 1 2
+           MESHX 0.0 40.0 350.0 660.0 700.0
+           SPLITX 1 4 4 1
+           ;
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9988235 ;
+TRACK FLUX := DELETE: TRACK FLUX ;
+*
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 CMFD ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9986206 ;
+GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
+*
+GEOM := GEO: :: CAR1D 2
+           X- VOID X+ REFL
+           MIX  2 1
+           MESHX 0.0 40.0 350.0
+           SPLITX 1 4
+           ;
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9977973 ;
+TRACK FLUX := DELETE: TRACK FLUX ;
+*
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 CMFD ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9977079 ;
+GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
+*
+MACRO := MAC: MACRO ::
+  ALBP 2 0.8 0.9 0.7 0.75
+ ;
+GEOM := GEO: :: CAR1D 4
+           X- ALBE 1 X+ ALBE 2
+           MIX  2 1 1 2
+           MESHX 0.0 40.0 350.0 660.0 700.0
+           SPLITX 1 4 4 1
+           ;
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9983991 ;
+TRACK FLUX := DELETE: TRACK FLUX ;
+*
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 HYPE 2 CMFD ;
+FLUX := NSSF: TRACK MACRO :: EDIT 1 EXTE 1000 1.0E-7 ;
+assertS FLUX :: K-EFFECTIVE 1 0.9982761 ;
+GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
+*----
+*  Steady-state calculations with the Analytic Nodal Method
+*----
+GEOM := GEO: :: CAR1D 3
+           X- ZERO X+ ZERO
+           MIX  2 1 2
+           MESHX 0.0 40.0 660.0 700.0
+           ;
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 99 MAXR 20 ANM ;
+FLUX := NSSF: TRACK MACRO :: 
+    EDIT 1 NUPD 100 1.0E-7 EXTE 100 1.0E-5 LEAK quadratic ;
+assertS FLUX :: K-EFFECTIVE 1 0.99773398 ;
+GEOM TRACK FLUX := DELETE: GEOM TRACK FLUX ;
+*
+GEOM := GEO: :: CAR1D 3
+           X- ALBE 1 X+ ALBE 2
+           MIX  2 1 2
+           MESHX 0.0 40.0 660.0 700.0
+           ;
+TRACK := NSST: GEOM :: 
+    TITLE 'BENCHMARK ENE6103'
+    EDIT 2 MAXR 20 ANM ;
+FLUX := NSSF: TRACK MACRO ::
+    EDIT 1 NUPD 100 1.0E-7 EXTE 100 1.0E-5 LEAK quadratic ;
+assertS FLUX :: K-EFFECTIVE 1 0.9984057 ;
+ECHO "test prob5p3_nem completed" ;
+END: ;