diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Trivac/data/DIFtst_proc | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Trivac/data/DIFtst_proc')
| -rwxr-xr-x | Trivac/data/DIFtst_proc/iaea2d.c2m | 125 | ||||
| -rwxr-xr-x | Trivac/data/DIFtst_proc/iaea3d.c2m | 157 | ||||
| -rwxr-xr-x | Trivac/data/DIFtst_proc/iaea_hexa.c2m | 176 | ||||
| -rwxr-xr-x | Trivac/data/DIFtst_proc/iaea_hexb.c2m | 72 | ||||
| -rwxr-xr-x | Trivac/data/DIFtst_proc/monju_diff.c2m | 348 | ||||
| -rwxr-xr-x | Trivac/data/DIFtst_proc/pertdiff.c2m | 243 | ||||
| -rwxr-xr-x | Trivac/data/DIFtst_proc/vv1k3d.c2m | 173 |
7 files changed, 1294 insertions, 0 deletions
diff --git a/Trivac/data/DIFtst_proc/iaea2d.c2m b/Trivac/data/DIFtst_proc/iaea2d.c2m new file mode 100755 index 0000000..9dc01ea --- /dev/null +++ b/Trivac/data/DIFtst_proc/iaea2d.c2m @@ -0,0 +1,125 @@ +*---- +* TEST CASE iaea2d +* IAEA 2D BENCHMARK IN DIFFUSION THEORY +* MACROLIB-DEFINED CROSS SECTIONS +*---- +* Define STRUCTURES and MODULES used +*---- +LINKED_LIST IAEA MACRO TRACK SYSTEM FLUX EDIT ; +MODULE GEO: MAC: BIVACT: BIVACA: TRIVAT: TRIVAA: FLUD: OUT: DELETE: + END: ; +PROCEDURE assertS ; +* +IAEA := GEO: :: CAR2D 9 9 + EDIT 2 + X- DIAG X+ VOID + Y- SYME Y+ DIAG + MIX 3 2 2 2 3 2 2 1 4 + 2 2 2 2 2 2 1 4 + 2 2 2 2 1 1 4 + 2 2 2 1 4 4 + 3 1 1 4 0 + 1 4 4 0 + 4 0 0 + 0 0 + 0 + MESHX 0.0 20.0 40.0 60.0 80.0 100.0 120.0 140.0 160.0 180.0 + ; +MACRO := MAC: :: + EDIT 2 NGRO 2 NMIX 4 NIFI 1 + READ INPUT + MIX 1 + DIFF 1.500E+00 4.0000E-01 + TOTAL 3.012E-02 8.0032E-02 + NUSIGF 0.000E+00 1.3500E-01 + H-FACTOR 0.000E+00 1.3500E-01 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 2 + DIFF 1.500E+00 4.0000E-01 + TOTAL 3.012E-02 8.5032E-02 + NUSIGF 0.000E+00 1.3500E-01 + H-FACTOR 0.000E+00 1.3500E-01 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 3 + DIFF 1.500E+00 4.00000E-01 + TOTAL 3.012E-02 1.30032E-01 + NUSIGF 0.000E+00 1.35000E-01 + H-FACTOR 0.000E+00 1.35000E-01 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 4 + DIFF 2.000E+00 3.0000E-01 + TOTAL 4.016E-02 1.0024E-02 + SCAT 1 1 0.0 2 2 0.0 0.4E-01 + ; +*---- +* BIVAC +*---- +TRACK := BIVACT: IAEA :: + TITLE 'IAEA-2D BENCHMARK' + EDIT 2 MAXR 81 DUAL 2 1 ; +SYSTEM := BIVACA: MACRO TRACK ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.029094 ; +EDIT := OUT: FLUX TRACK MACRO IAEA :: + EDIT 2 INTG + 1 2 3 4 5 6 7 8 0 + 9 10 11 12 13 14 15 0 + 16 17 18 19 20 21 0 + 22 23 24 25 0 0 + 26 27 28 0 0 + 29 0 0 0 + 0 0 0 + 0 0 + 0 + ; +TRACK SYSTEM FLUX EDIT := DELETE: TRACK SYSTEM FLUX EDIT ; +* +TRACK := BIVACT: IAEA :: + TITLE 'IAEA-2D BENCHMARK' + EDIT 2 MAXR 81 PRIM 2 2 ; +SYSTEM := BIVACA: MACRO TRACK ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.032310 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +*---- +* TRIVAC +*---- +TRACK := TRIVAT: IAEA :: + TITLE 'IAEA-2D BENCHMARK' + EDIT 2 MAXR 81 DUAL 2 1 ; +SYSTEM := TRIVAA: MACRO TRACK ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADJ ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.029094 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +TRACK := TRIVAT: IAEA :: + TITLE 'IAEA-2D BENCHMARK' + EDIT 2 MAXR 81 DUAL 2 2 ; +SYSTEM := TRIVAA: MACRO TRACK ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.028690 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +TRACK := TRIVAT: IAEA :: + TITLE 'IAEA-2D BENCHMARK' + EDIT 2 MAXR 81 PRIM 2 ; +SYSTEM := TRIVAA: MACRO TRACK ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.032310 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +TRACK := TRIVAT: IAEA :: + TITLE 'IAEA-2D BENCHMARK' + EDIT 2 MAXR 81 MCFD 2 ; +SYSTEM := TRIVAA: MACRO TRACK ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.028690 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +ECHO "test iaea2d completed" ; +END: ; diff --git a/Trivac/data/DIFtst_proc/iaea3d.c2m b/Trivac/data/DIFtst_proc/iaea3d.c2m new file mode 100755 index 0000000..b88125d --- /dev/null +++ b/Trivac/data/DIFtst_proc/iaea3d.c2m @@ -0,0 +1,157 @@ +LINKED_LIST IAEA3D MACRO TRACK SYSTEM FLUX EDIT REF ; +MODULE GEO: MAC: TRIVAT: TRIVAA: FLUD: OUT: DELETE: END: ; +PROCEDURE assertS ; +* +IAEA3D := GEO: :: CAR3D 9 9 4 + EDIT 2 + X- DIAG X+ VOID + Y- SYME Y+ DIAG + Z- VOID Z+ VOID + MESHX 0.0 20.0 40.0 60.0 80.0 100.0 120.0 140.0 160.0 180.0 + MESHZ 0.0 20.0 280.0 360.0 380.0 + SPLITZ 1 2 1 1 + ! PLANE NB 1 + MIX 4 4 4 4 4 4 4 4 4 + 4 4 4 4 4 4 4 4 + 4 4 4 4 4 4 4 + 4 4 4 4 4 4 + 4 4 4 4 0 + 4 4 4 0 + 4 0 0 + 0 0 + 0 + ! PLANE NB 2 + 3 2 2 2 3 2 2 1 4 + 2 2 2 2 2 2 1 4 + 2 2 2 2 1 1 4 + 2 2 2 1 4 4 + 3 1 1 4 0 + 1 4 4 0 + 4 0 0 + 0 0 + 0 + ! PLANE NB 3 + 3 2 2 2 3 2 2 1 4 + 2 2 2 2 2 2 1 4 + 3 2 2 2 1 1 4 + 2 2 2 1 4 4 + 3 1 1 4 0 + 1 4 4 0 + 4 0 0 + 0 0 + 0 + ! PLANE NB 4 + 5 4 4 4 5 4 4 4 4 + 4 4 4 4 4 4 4 4 + 5 4 4 4 4 4 4 + 4 4 4 4 4 4 + 5 4 4 4 0 + 4 4 4 0 + 4 0 0 + 0 0 + 0 + ; +MACRO := MAC: :: + EDIT 2 NGRO 2 NMIX 5 NIFI 1 + READ INPUT + MIX 1 + DIFFX 1.500E+00 4.0000E-01 + TOTAL 3.000E-02 8.0000E-02 + NUSIGF 0.000E+00 1.3500E-01 + H-FACTOR 0.000E+00 1.3500E-01 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 2 + DIFFX 1.500E+00 4.0000E-01 + TOTAL 3.000E-02 8.5000E-02 + NUSIGF 0.000E+00 1.3500E-01 + H-FACTOR 0.000E+00 1.3500E-01 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 3 + DIFFX 1.500E+00 4.00000E-01 + TOTAL 3.000E-02 1.30000E-01 + NUSIGF 0.000E+00 1.35000E-01 + H-FACTOR 0.000E+00 1.35000E-01 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 4 + DIFFX 2.000E+00 3.0000E-01 + TOTAL 4.000E-02 1.0000E-02 + SCAT 1 1 0.0 2 2 0.0 0.4E-01 + MIX 5 + DIFFX 2.000E+00 3.0000E-01 + TOTAL 4.000E-02 5.5000E-02 + SCAT 1 1 0.0 2 2 0.0 0.4E-01 + ; +TRACK := TRIVAT: IAEA3D :: + TITLE 'TEST IAEA 3D' + EDIT 5 MAXR 405 DUAL 3 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 5 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.028981 ; +EDIT := OUT: FLUX TRACK MACRO IAEA3D :: + EDIT 2 INTG + (*PLANE NB 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 0 0 0 0 0 + 0 0 0 0 0 + 0 0 0 0 + 0 0 0 + 0 0 + 0 + (*PLANE NB 2*) + 1 2 3 4 5 6 7 8 0 + 9 10 11 12 13 14 15 0 + 16 17 18 19 20 21 0 + 22 23 24 25 0 0 + 26 27 28 0 0 + 29 0 0 0 + 0 0 0 + 0 0 + 0 + (*PLANE NB 3*) + 30 31 32 33 34 35 36 37 0 + 38 39 40 41 42 43 44 0 + 45 46 47 48 49 50 0 + 51 52 53 54 0 0 + 55 56 57 0 0 + 58 0 0 0 + 0 0 0 + 0 0 + 0 + (*PLANE NB 4*) + 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 0 0 0 0 0 + 0 0 0 0 0 + 0 0 0 0 + 0 0 0 + 0 0 + 0 + ; +TRACK SYSTEM FLUX EDIT := DELETE: TRACK SYSTEM FLUX EDIT ; +* +TRACK := TRIVAT: IAEA3D :: + TITLE 'TEST IAEA 3D' + EDIT 5 MAXR 405 PRIM 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 5 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.029308 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +TRACK := TRIVAT: IAEA3D :: + TITLE 'TEST IAEA 3D' + EDIT 5 MAXR 405 MCFD 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 5 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.028842 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +ECHO "test iaea3d completed" ; +END: ; diff --git a/Trivac/data/DIFtst_proc/iaea_hexa.c2m b/Trivac/data/DIFtst_proc/iaea_hexa.c2m new file mode 100755 index 0000000..443cf1d --- /dev/null +++ b/Trivac/data/DIFtst_proc/iaea_hexa.c2m @@ -0,0 +1,176 @@ +LINKED_LIST HEX MACRO TRACK SYSTEM FLUX ; +MODULE GEO: MAC: BIVACT: BIVACA: TRIVAT: TRIVAA: FLUD: DELETE: END: ; +PROCEDURE assertS ; +* +HEX := GEO: :: HEX 16 + EDIT 2 + HBC S30 ALBE 0.0 + SIDE 11.5470054 + MIX + 3 + 2 + 2 3 + 2 2 + 3 2 2 + 2 2 2 + 1 1 1 1 + ; +MACRO := MAC: :: + EDIT 2 NGRO 2 NMIX 3 NIFI 1 + READ INPUT + MIX 1 + DIFF 1.5E+00 4.0E-01 + TOTAL 3.0E-02 8.0E-02 + NUSIGF 0.0E+00 1.35E-01 + H-FACTOR 0.0E+00 5.6E-02 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 2 + DIFF 1.5E+00 4.0E-01 + TOTAL 3.0E-02 8.5E-02 + NUSIGF 0.0E+00 1.35E-01 + H-FACTOR 0.0E+00 5.6E-02 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 3 + DIFF 1.5E+00 4.0E-01 + TOTAL 3.0E-02 1.3E-01 + NUSIGF 0.0E+00 1.35E-01 + H-FACTOR 0.0E+00 5.6E-02 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + ; +*---- +* BIVAC +*---- +HEX := GEO: HEX :: + SPLITH 0 SPLITL 0 ; +TRACK := BIVACT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 500 MCFD ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9917505 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 0 SPLITL 0 ; +TRACK := BIVACT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 500 PRIM (*IELEM=*) 1 (*ICOL=*) 2 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9610498 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 1 SPLITL 0 ; +TRACK := BIVACT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 500 MCFD ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9822278 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 1 SPLITL 0 ; +TRACK := BIVACT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 500 PRIM (*IELEM=*) 1 (*ICOL=*) 2 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9709085 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 0 SPLITL 1 ; +TRACK := BIVACT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 5000 DUAL (*IELEM=*) 2 (*ICOL=*) 3 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9782231 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 0 SPLITL 2 ; +TRACK := BIVACT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 5000 DUAL (*IELEM=*) 1 (*ICOL=*) 3 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9775723 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +*---- +* TRIVAC +*---- +HEX := GEO: HEX :: + SPLITH 0 SPLITL 0 ; +TRACK := TRIVAT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 500 MCFD (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9917505 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 0 SPLITL 0 ; +TRACK := TRIVAT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 2 MAXR 500 PRIM (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9610475 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 1 SPLITL 0 ; +TRACK := TRIVAT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 500 MCFD (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9822259 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 1 SPLITL 0 ; +TRACK := TRIVAT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 500 PRIM (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9709070 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX := GEO: HEX :: + SPLITH 0 SPLITL 2 ; +TRACK := TRIVAT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE A (NO REFLECTOR).' + EDIT 5 MAXR 5000 DUAL (*IELEM=*) 2 (*ICOL=*) 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 4 ; +assertS FLUX :: 'K-EFFECTIVE' 1 0.9780799 ; +ECHO "test iaea_hexa completed" ; +END: ; diff --git a/Trivac/data/DIFtst_proc/iaea_hexb.c2m b/Trivac/data/DIFtst_proc/iaea_hexb.c2m new file mode 100755 index 0000000..47a42f9 --- /dev/null +++ b/Trivac/data/DIFtst_proc/iaea_hexb.c2m @@ -0,0 +1,72 @@ +LINKED_LIST HEX MACRO TRACK SYSTEM FLUX EDIT ; +MODULE GEO: MAC: TRIVAT: TRIVAA: FLUD: OUT: DELETE: END: ; +PROCEDURE assertS ; +* +HEX := GEO: :: HEX 20 + EDIT 2 + HBC S30 ALBE 0.0 + SIDE 11.5470054 + MIX + 3 + 2 + 2 3 + 2 2 + 3 2 2 + 2 2 2 + 1 1 1 1 + 4 4 4 4 + ; +MACRO := MAC: :: + EDIT 2 NGRO 2 NMIX 4 NIFI 1 + READ INPUT + MIX 1 + DIFF 1.5E+00 4.0E-01 + TOTAL 3.0E-02 8.0E-02 + NUSIGF 0.0E+00 1.35E-01 + H-FACTOR 0.0E+00 5.6E-02 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 2 + DIFF 1.5E+00 4.0E-01 + TOTAL 3.0E-02 8.5E-02 + NUSIGF 0.0E+00 1.35E-01 + H-FACTOR 0.0E+00 5.6E-02 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 3 + DIFF 1.5E+00 4.0E-01 + TOTAL 3.0E-02 1.3E-01 + NUSIGF 0.0E+00 1.35E-01 + H-FACTOR 0.0E+00 5.6E-02 + SCAT 1 1 0.0 2 2 0.0 0.2E-01 + MIX 4 + DIFF 1.5E+00 4.0E-01 + TOTAL 4.0E-02 1.0E-02 + SCAT 1 1 0.0 2 2 0.0 0.4E-01 + ; +HEX := GEO: HEX :: + SPLITH 0 SPLITL 0 ; +TRACK := TRIVAT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE B (WITH REFLECTOR).' + EDIT 5 MAXR 500 MCFD (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +EDIT := OUT: FLUX HEX MACRO TRACK :: + EDIT 2 INTG IN ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.009611 ; +TRACK SYSTEM FLUX EDIT := DELETE: TRACK SYSTEM FLUX EDIT ; +* +HEX := GEO: HEX :: + SPLITH 0 SPLITL 2 ; +TRACK := TRIVAT: HEX :: + TITLE 'TEST HEXAGONAL IAEA BENCHMARK CASE B (WITH REFLECTOR).' + EDIT 5 MAXR 5000 DUAL (*IELEM=*) 2 (*ICOL=*) 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + ADI 2 EDIT 2 ; +EDIT := OUT: FLUX HEX MACRO TRACK :: + EDIT 2 INTG IN ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.005512 ; +ECHO "test iaea_hexb completed" ; +END: ; diff --git a/Trivac/data/DIFtst_proc/monju_diff.c2m b/Trivac/data/DIFtst_proc/monju_diff.c2m new file mode 100755 index 0000000..45b8d1d --- /dev/null +++ b/Trivac/data/DIFtst_proc/monju_diff.c2m @@ -0,0 +1,348 @@ +LINKED_LIST HEX2D HEX3D MACRO TRACK SYSTEM FLUX ; +MODULE GEO: MAC: BIVACT: BIVACA: TRIVAT: TRIVAA: FLUD: DELETE: END: ; +PROCEDURE assertS ; +* +MACRO := MAC: :: + EDIT 2 NGRO 3 NMIX 5 NIFI 1 + READ INPUT + MIX 1 (* ACTIVE INNER CORE *) + DIFF 2.540E+00 1.724E+00 1.264E+00 + TOTAL 3.09865E-02 9.490E-03 7.333E-03 + NUSIGF 1.235E-02 5.225E-03 7.684E-03 + H-FACTOR 1.235E-02 5.225E-03 7.684E-03 + SCAT 1 1 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) 2.544E-02 + 3 3 (*3->3*) 0.0 (*2->3*) 6.551E-03 (*1->3*) 5.625E-04 + MIX 2 (* ACTIVE OUTER CORE *) + DIFF 2.548E+00 1.725E+00 1.269E+00 + TOTAL 3.12138E-02 9.875E-03 8.099E-03 + NUSIGF 1.467E-02 6.955E-03 9.986E-03 + H-FACTOR 1.467E-02 6.955E-03 9.986E-03 + SCAT 1 1 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) 2.497E-02 + 3 3 (*3->3*) 0.0 (*2->3*) 6.341E-03 (*1->3*) 5.548E-04 + MIX 3 (* RADIAL/AXIAL BLANKET *) + DIFF 2.173E+00 1.439E+00 1.026E+00 + TOTAL 3.79308E-02 1.1843E-02 7.611E-03 + NUSIGF 8.631E-03 5.995E-04 1.381E-03 + H-FACTOR 8.631E-03 5.995E-04 1.381E-03 + SCAT 1 1 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) 3.288E-02 + 3 3 (*3->3*) 0.0 (*2->3*) 1.000E-02 (*1->3*) 7.468E-04 + MIX 4 (* CONTROL ROD *) + DIFF 2.500E+00 1.681E+00 1.269E+00 + TOTAL 2.32803E-02 1.2727E-02 1.497E-02 + SCAT 1 1 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) 2.185E-02 + 3 3 (*3->3*) 0.0 (*2->3*) 9.379E-03 (*1->3*) 2.163E-04 + MIX 5 (* SODIUM CHANNEL *) + DIFF 4.805E+00 3.262E+00 2.431E+00 + TOTAL 1.152508E-02 3.64874E-03 3.072E-04 + SCAT 1 1 (*1->1*) 0.0 + 2 2 (*2->2*) 0.0 (*1->2*) 1.130E-02 + 3 3 (*3->3*) 0.0 (*2->3*) 3.571E-03 (*1->3*) 6.718E-05 + ; +* +HEX2D := GEO: :: HEX 133 + EDIT 2 + HBC R120 VOID + SIDE 6.67417 + MIX + 4 + 1 1 + 1 1 1 1 + 1 4 1 1 4 1 + 1 1 1 1 1 1 1 1 + 4 1 1 1 1 4 1 1 1 1 + 1 1 1 4 1 1 1 1 1 4 1 1 + 2 2 2 2 2 2 2 2 2 2 2 2 2 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 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 0 3 3 3 3 3 3 3 3 3 3 0 + ; +* +HEX3D := GEO: :: HEXZ 133 4 + EDIT 2 + HBC R120 VOID + SIDE 6.67417 + Z- VOID Z+ VOID + MESHZ 0.0 30.0 79.0 123.0 158.0 + SPLITZ 1 2 2 1 + MIX + (* UPPER BLANKET *) + 4 + 3 3 + 3 3 3 3 + 3 4 3 3 4 3 + 3 3 3 3 3 3 3 3 + 4 3 3 3 3 4 3 3 3 3 + 3 3 3 4 3 3 3 3 3 4 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 0 3 3 3 3 3 3 3 3 3 3 0 + (* UPPER INNER/OUTER CORE *) + 4 + 1 1 + 1 1 1 1 + 1 4 1 1 4 1 + 1 1 1 1 1 1 1 1 + 4 1 1 1 1 4 1 1 1 1 + 1 1 1 4 1 1 1 1 1 4 1 1 + 2 2 2 2 2 2 2 2 2 2 2 2 2 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 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 0 3 3 3 3 3 3 3 3 3 3 0 + (* LOWER INNER/OUTER CORE *) + 4 + 1 1 + 1 1 1 1 + 1 4 1 1 4 1 + 1 1 1 1 1 1 1 1 + 5 1 1 1 1 5 1 1 1 1 + 1 1 1 5 1 1 1 1 1 5 1 1 + 2 2 2 2 2 2 2 2 2 2 2 2 2 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 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 0 3 3 3 3 3 3 3 3 3 3 0 + (* LOWER BLANKET *) + 4 + 3 3 + 3 3 3 3 + 3 4 3 3 4 3 + 3 3 3 3 3 3 3 3 + 5 3 3 3 3 5 3 3 3 3 + 3 3 3 5 3 3 3 3 3 5 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 + 3 3 3 3 3 3 3 3 3 3 0 3 3 3 3 3 3 3 3 3 3 0 + ; +*---- +* Mesh-corner finite differences +*---- +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 0 ; +TRACK := BIVACT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 PRIM (*IELEM=*) 1 (*ICOL=*) 2 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.145389 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 1 SPLITL 0 ; +TRACK := BIVACT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 PRIM (*IELEM=*) 1 (*ICOL=*) 2 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.151643 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 0 ; +TRACK := TRIVAT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 PRIM (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.145392 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 1 SPLITL 0 ; +TRACK := TRIVAT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 PRIM (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.151643 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX3D := GEO: HEX3D :: + SPLITH 0 SPLITL 0 ; +TRACK := TRIVAT: HEX3D :: + TITLE 'TEST 3D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 PRIM (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.007556 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX3D := GEO: HEX3D :: + SPLITH 1 SPLITL 0 ; +TRACK := TRIVAT: HEX3D :: + TITLE 'TEST 3D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 PRIM (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.012467 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +*---- +* Mesh-centered finite differences +*---- +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 0 ; +TRACK := BIVACT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 MCFD ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.167697 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 1 SPLITL 0 ; +TRACK := BIVACT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 MCFD ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.156221 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 0 ; +TRACK := TRIVAT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 MCFD (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.167697 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 1 SPLITL 0 ; +TRACK := TRIVAT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 MCFD (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.156212 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX3D := GEO: HEX3D :: + SPLITH 0 SPLITL 0 ; +TRACK := TRIVAT: HEX3D :: + TITLE 'TEST 3D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 MCFD (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.066331 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX3D := GEO: HEX3D :: + SPLITH 1 SPLITL 0 ; +TRACK := TRIVAT: HEX3D :: + TITLE 'TEST 3D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 5000 MCFD (*IELEM=*) 1 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.057418 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +*---- +* Thomas-Raviart-Schneider +*---- +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 1 ; +TRACK := BIVACT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 8000 DUAL (*IELEM=*) 1 (*ICOL=*) 3 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.152478 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 2 ; +TRACK := BIVACT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 8000 DUAL (*IELEM=*) 1 (*ICOL=*) 3 ; +SYSTEM := BIVACA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.153644 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 2 ; +TRACK := TRIVAT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 8000 DUAL (*IELEM=*) 1 (*ICOL=*) 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.153649 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX2D := GEO: HEX2D :: + SPLITH 0 SPLITL 1 ; +TRACK := TRIVAT: HEX2D :: + TITLE 'TEST 2D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 8000 DUAL (*IELEM=*) 2 (*ICOL=*) 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.154039 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX3D := GEO: HEX3D :: + SPLITH 0 SPLITL 1 ; +TRACK := TRIVAT: HEX3D :: + TITLE 'TEST 3D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 8000 DUAL (*IELEM=*) 1 (*ICOL=*) 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.029397 ; +TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ; +* +HEX3D := GEO: HEX3D :: + SPLITH 0 SPLITL 2 ; +TRACK := TRIVAT: HEX3D :: + TITLE 'TEST 3D FBR MONJU (3 GROUPS).' + EDIT 2 MAXR 30000 DUAL (*IELEM=*) 1 (*ICOL=*) 3 ; +SYSTEM := TRIVAA: MACRO TRACK :: + EDIT 2 ; +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADI 6 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.030225 ; +ECHO "test monju_diff completed" ; +END: ; diff --git a/Trivac/data/DIFtst_proc/pertdiff.c2m b/Trivac/data/DIFtst_proc/pertdiff.c2m new file mode 100755 index 0000000..2196541 --- /dev/null +++ b/Trivac/data/DIFtst_proc/pertdiff.c2m @@ -0,0 +1,243 @@ +*----- +* GPT TEST pertdiff +*----- +LINKED_LIST IAEA MACRO_REF MACRO_0 DMACRO TRACK SYSTEM_REF FLUX DSYSTEM + DSOUR DASOUR DFLUX SYSTEM_0 FLUX2 EDIT ADFLUX MACRO ; +MODULE GEO: MAC: TRIVAT: TRIVAA: FLUD: DELETE: DELTA: GPTFLU: OUT: + END: ADD: UTL: STAT: ; +PROCEDURE assertS ; +* +IAEA := GEO: :: CAR2D 9 9 + EDIT 2 + X- DIAG X+ VOID + Y- SYME Y+ DIAG + MIX 3 2 2 2 3 2 2 1 4 + 2 2 2 2 2 2 1 4 + 2 2 2 2 1 1 4 + 2 2 2 1 4 4 + 3 1 1 4 0 + 1 4 4 0 + 4 0 0 + 0 0 + 0 + MESHX 0.0 20.0 40.0 60.0 80.0 100.0 120.0 140.0 160.0 180.0 + ; +* +* REFERENCE CASE: +MACRO_REF := MAC: :: + EDIT 2 NGRO 2 NMIX 4 NIFI 1 + READ INPUT + MIX 1 + DIFF 1.900000E+00 4.400000E-01 + TOTAL 3.012000E-02 8.303201E-02 + NUSIGF 0.000000E+00 1.650000E-01 + H-FACTOR 0.000000E+00 1.650000E-01 + SCAT 1 1 0.0 2 2 0.0 0.200000E-01 + MIX 2 + DIFF 1.500000E+00 4.000000E-01 + TOTAL 3.012000E-02 8.503199E-02 + NUSIGF 0.000000E+00 1.350000E-01 + H-FACTOR 0.000000E+00 1.350000E-01 + SCAT 1 1 0.0 2 2 0.0 0.200000E-01 + MIX 3 + DIFF 1.500000E+00 4.000000E-01 + TOTAL 2.912000E-02 1.260320E-01 + NUSIGF 0.000000E+00 1.350000E-01 + H-FACTOR 0.000000E+00 1.350000E-01 + SCAT 1 1 0.0 2 2 0.0 0.200000E-01 + MIX 4 + DIFF 2.000000E+00 3.000000E-01 + TOTAL 4.016000E-02 1.002400E-02 + SCAT 1 1 0.0 2 2 0.0 0.400000E-01 + ; +* +* UNPERTURBED CASE: +MACRO_0 := MAC: :: + EDIT 2 NGRO 2 NMIX 4 NIFI 1 + READ INPUT + MIX 1 + DIFF 1.500000E+00 4.000000E-01 + TOTAL 3.012000E-02 8.003199E-02 + NUSIGF 0.000000E+00 1.350000E-01 + H-FACTOR 0.000000E+00 1.350000E-01 + SCAT 1 1 0.0 2 2 0.0 0.200000E-01 + MIX 2 + DIFF 1.500000E+00 4.000000E-01 + TOTAL 3.012000E-02 8.503199E-02 + NUSIGF 0.000000E+00 1.350000E-01 + H-FACTOR 0.000000E+00 1.350000E-01 + SCAT 1 1 0.0 2 2 0.0 0.200000E-01 + MIX 3 + DIFF 1.500000E+00 4.000000E-01 + TOTAL 3.012000E-02 1.300320E-01 + NUSIGF 0.000000E+00 1.350000E-01 + H-FACTOR 0.000000E+00 1.350000E-01 + SCAT 1 1 0.0 2 2 0.0 0.200000E-01 + MIX 4 + DIFF 2.000000E+00 3.000000E-01 + TOTAL 4.016000E-02 1.002400E-02 + SCAT 1 1 0.0 2 2 0.0 0.400000E-01 + ; +* +* PERTURBATION: +DMACRO := MAC: :: + EDIT 2 NGRO 2 NMIX 4 NIFI 1 + READ INPUT + MIX 1 + DIFF 4.000000E-01 4.000000E-02 + TOTAL 0.000000E+00 3.000000E-03 + NUSIGF 0.000000E+00 3.000000E-02 + H-FACTOR 0.000000E+00 3.000000E-02 + SCAT 1 1 0.0 2 2 0.0 1.0E-10 + MIX 2 + SCAT 1 1 0.0 2 2 0.0 1.0E-10 + MIX 3 + TOTAL -9.999999E-04 -4.000001E-03 + SCAT 1 1 0.0 2 2 0.0 1.0E-10 + MIX 4 + SCAT 1 1 0.0 2 2 0.0 1.0E-10 + ; +*---- +* MCFD1 case, Livolant iteration +*---- +TRACK := TRIVAT: IAEA :: + TITLE 'MODIFIED TEST IAEA-2D (ANL VERSION)' + EDIT 5 MAXR 1156 + MCFD 1 ; +SYSTEM_REF := TRIVAA: MACRO_REF TRACK :: EDIT 5 ; +FLUX := FLUD: SYSTEM_REF TRACK :: EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.113392 ; +FLUX := DELETE: FLUX ; +* +SYSTEM_0 := TRIVAA: MACRO_0 TRACK :: EDIT 5 ; +FLUX := FLUD: SYSTEM_0 TRACK :: EDIT 2 ADJ ; +DSYSTEM := TRIVAA: MACRO_0 TRACK DMACRO :: EDIT 5 PERT ; +DSOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ; +DFLUX := GPTFLU: DSOUR FLUX SYSTEM_0 TRACK :: EDIT 2 EXPLICIT + FROM-TO 1 1 ; +assertS DFLUX :: 'K-EFFECTIVE' 1 1.032713 ; +DASOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ADJ ; +ADFLUX := GPTFLU: DASOUR FLUX SYSTEM_0 TRACK :: EDIT 2 IMPLICIT + FROM-TO 1 1 ; +assertS ADFLUX :: 'K-EFFECTIVE' 1 1.032713 ; +* Reset the perturbation flag of DSYSTEM to 0 so that it can be added. +DSYSTEM := UTL: DSYSTEM :: CREA STATE-VECTOR 9 9 = 0 ; +SYSTEM_0 := ADD: SYSTEM_0 DSYSTEM ; +MACRO := ADD: MACRO_0 DMACRO ; +* +STAT: SYSTEM_REF SYSTEM_0 ; +STAT: MACRO_REF MACRO ; +FLUX2 := FLUD: SYSTEM_0 TRACK MACRO :: EDIT 2 ; +assertS FLUX2 :: 'K-EFFECTIVE' 1 1.113392 ; +EDIT := OUT: FLUX2 TRACK MACRO IAEA :: + EDIT 2 INTG + 1 2 3 4 5 6 7 8 0 + 9 10 11 12 13 14 15 0 + 16 17 18 19 20 21 0 + 22 23 24 25 0 0 + 26 27 28 0 0 + 29 0 0 0 + 0 0 0 + 0 0 + 0 + ; +TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2 +EDIT ADFLUX MACRO +:= DELETE: +TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2 +EDIT ADFLUX MACRO ; +*---- +* RAVIART-THOMAS DUAL 1 2 case, Livolant iteration +*---- +TRACK := TRIVAT: IAEA :: + TITLE 'MODIFIED TEST IAEA-2D (ANL VERSION)' + EDIT 5 MAXR 1156 + DUAL 1 2 ; +SYSTEM_REF := TRIVAA: MACRO_REF TRACK :: EDIT 5 ; +FLUX := FLUD: SYSTEM_REF TRACK :: EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.113392 ; +FLUX := DELETE: FLUX ; +* +SYSTEM_0 := TRIVAA: MACRO_0 TRACK :: EDIT 5 ; +FLUX := FLUD: SYSTEM_0 TRACK :: EDIT 2 ADJ ; +DSYSTEM := TRIVAA: MACRO_0 TRACK DMACRO :: EDIT 5 PERT ; +DSOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ; +DFLUX := GPTFLU: DSOUR FLUX SYSTEM_0 TRACK :: EDIT 2 EXPLICIT + FROM-TO 1 1 ; +assertS DFLUX :: 'K-EFFECTIVE' 1 1.032713 ; +DASOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ADJ ; +ADFLUX := GPTFLU: DASOUR FLUX SYSTEM_0 TRACK :: EDIT 2 IMPLICIT + FROM-TO 1 1 ; +assertS ADFLUX :: 'K-EFFECTIVE' 1 1.032713 ; +* Reset the perturbation flag of DSYSTEM to 0 so that it can be added. +DSYSTEM := UTL: DSYSTEM :: CREA STATE-VECTOR 9 9 = 0 ; +SYSTEM_0 := ADD: SYSTEM_0 DSYSTEM ; +MACRO := ADD: MACRO_0 DMACRO ; +* +STAT: SYSTEM_REF SYSTEM_0 ; +STAT: MACRO_REF MACRO ; +FLUX2 := FLUD: SYSTEM_0 TRACK MACRO :: EDIT 2 ; +assertS FLUX2 :: 'K-EFFECTIVE' 1 1.113392 ; +EDIT := OUT: FLUX2 TRACK MACRO IAEA :: + EDIT 2 INTG + 1 2 3 4 5 6 7 8 0 + 9 10 11 12 13 14 15 0 + 16 17 18 19 20 21 0 + 22 23 24 25 0 0 + 26 27 28 0 0 + 29 0 0 0 + 0 0 0 + 0 0 + 0 + ; +TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2 +EDIT ADFLUX MACRO +:= DELETE: +TRACK SYSTEM_REF FLUX DSYSTEM DSOUR DASOUR DFLUX SYSTEM_0 FLUX2 +EDIT ADFLUX MACRO ; +*---- +* RAVIART-THOMAS DUAL 1 2 case, GMRES iteration +*---- +TRACK := TRIVAT: IAEA :: + TITLE 'MODIFIED TEST IAEA-2D (ANL VERSION)' + EDIT 5 MAXR 1156 + DUAL 1 2 ; +SYSTEM_REF := TRIVAA: MACRO_REF TRACK :: EDIT 5 ; +FLUX := FLUD: SYSTEM_REF TRACK :: EDIT 2 ; +assertS FLUX :: 'K-EFFECTIVE' 1 1.113392 ; +FLUX := DELETE: FLUX ; +* +SYSTEM_0 := TRIVAA: MACRO_0 TRACK :: EDIT 5 ; +FLUX := FLUD: SYSTEM_0 TRACK :: EDIT 2 ADJ ; +DSYSTEM := TRIVAA: MACRO_0 TRACK DMACRO :: EDIT 5 PERT ; +DSOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ; +DFLUX := GPTFLU: DSOUR FLUX SYSTEM_0 TRACK :: EDIT 2 GMRES 10 EXPLICIT + FROM-TO 1 1 ; +assertS DFLUX :: 'K-EFFECTIVE' 1 1.032713 ; +DASOUR := DELTA: FLUX SYSTEM_0 DSYSTEM TRACK :: EDIT 2 ADJ ; +ADFLUX := GPTFLU: DASOUR FLUX SYSTEM_0 TRACK :: EDIT 2 GMRES 10 IMPLICIT + FROM-TO 1 1 ; +assertS ADFLUX :: 'K-EFFECTIVE' 1 1.032713 ; +* Reset the perturbation flag of DSYSTEM to 0 so that it can be added. +DSYSTEM := UTL: DSYSTEM :: CREA STATE-VECTOR 9 9 = 0 ; +SYSTEM_0 := ADD: SYSTEM_0 DSYSTEM ; +MACRO := ADD: MACRO_0 DMACRO ; +* +STAT: SYSTEM_REF SYSTEM_0 ; +STAT: MACRO_REF MACRO ; +FLUX2 := FLUD: SYSTEM_0 TRACK MACRO :: EDIT 2 ; +assertS FLUX2 :: 'K-EFFECTIVE' 1 1.113392 ; +EDIT := OUT: FLUX2 TRACK MACRO IAEA :: + EDIT 2 INTG + 1 2 3 4 5 6 7 8 0 + 9 10 11 12 13 14 15 0 + 16 17 18 19 20 21 0 + 22 23 24 25 0 0 + 26 27 28 0 0 + 29 0 0 0 + 0 0 0 + 0 0 + 0 + ; +ECHO "test pertdiff completed" ; +END: ; diff --git a/Trivac/data/DIFtst_proc/vv1k3d.c2m b/Trivac/data/DIFtst_proc/vv1k3d.c2m new file mode 100755 index 0000000..683ada5 --- /dev/null +++ b/Trivac/data/DIFtst_proc/vv1k3d.c2m @@ -0,0 +1,173 @@ +*------------------------ +* Benchmark VV1K3D +* +* Reference: A. Bernal, A. Hebert, J. E. Roman, R. Miro and G. Verdu, +* "A Krylov-Schur solution of the eigenvalue problem for the neutron +* diffusion equation discretized with the Raviart-Thomas method," J. of +* Nuclear Science and Technology, Vol. 54, No. 10, 1085-1094 (2017). +*------------------------ +LINKED_LIST GEOM XSEC TRACK SYSTEM FLUX RESU ; +MODULE GEO: MAC: TRIVAT: TRIVAA: FLUD: OUT: END: ABORT: ; +PROCEDURE assertS ; + +GEOM := GEO: :: + HEXZ 169 10 + EDIT 2 + Z- ZERO Z+ ZERO HBC COMPLETE ZERO + MESHZ 0. 20. 40. 60. 80. 100. 120. 140. 160. 180. 200. + SIDE 13.6255 + SPLITL 1 + MIX + ! PLANE 1 + 2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 2 + 2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 3 + 2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 4 + 2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 5 + 2 1 1 1 1 1 1 1 4 1 4 1 4 1 4 1 4 1 4 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 6 + 2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 7 + 2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 8 + 2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 9 + 2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ! PLANE 10 + 2 1 1 1 1 1 1 1 3 1 3 1 3 1 3 1 3 1 3 3 1 1 + 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3 2 1 4 1 4 + 1 2 3 1 4 1 2 3 1 3 2 1 4 1 4 1 2 3 1 3 2 1 4 + 5 4 1 3 1 3 5 3 1 3 1 4 5 3 1 3 1 4 5 4 1 3 1 3 + 5 3 1 3 1 4 5 4 1 3 1 3 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 + 5 5 5 5 5 5 5 5 + ; + +XSEC := MAC: :: + EDIT 2 + NGRO 2 NIFI 1 NMIX 5 + READ INPUT + MIX 1 + TOTAL 2.48836E-02 6.73049E-2 + NUSIGF 4.81619E-3 8.46154E-2 + DIFF 1.38320 3.86277E-1 + CHI 1. 0. + SCAT 1 1 0. + 1 1 1.64977E-2 + MIX 2 + TOTAL 2.62865E-02 8.10328E-2 + NUSIGF 4.66953E-3 8.52264E-2 + DIFF 1.38299 3.89403E-1 + CHI 1. 0. + SCAT 1 1 0. + 1 1 1.47315E-2 + MIX 3 + TOTAL 2.45662E-02 8.44801E-2 + NUSIGF 6.04889E-3 1.19428E-1 + DIFF 1.39522 3.86225E-1 + CHI 1. 0. + SCAT 1 1 0. + 1 1 1.56219E-2 + MIX 4 + TOTAL 2.60117E-02 9.89671E-2 + NUSIGF 5.91507E-3 1.20497E-1 + DIFF 1.39446 3.87723E-1 + CHI 1. 0. + SCAT 1 1 0. + 1 1 1.40185E-2 + MIX 5 + TOTAL 2.46141E-02 8.93878E-2 + NUSIGF 6.40256E-3 1.29281E-1 + DIFF 1.39506 3.84492E-1 + CHI 1. 0. + SCAT 1 1 0. + 1 1 1.54981E-2 + ; + +TRACK := TRIVAT: GEOM :: + EDIT 2 + TITL 'VV1K3D' + MAXR 10000 DUAL 2 1 + ; + +SYSTEM := TRIVAA: XSEC TRACK :: EDIT 3 ; + +FLUX := FLUD: SYSTEM TRACK :: + EDIT 2 ADJ IRAM 3 6 EXTE 30 1.0E-8 + ; + +RESU := OUT: FLUX SYSTEM TRACK XSEC GEOM :: + EDIT 2 + MODE 1 POWR 1.0 INTG MIX + ; + +assertS FLUX :: 'K-EFFECTIVE' 1 1.005450 ; +ECHO "test vv1k3d completed" ; +END: ; |