blob: 3a0afa665a73e39fad5421d50ef7f4cfac7e522e (
plain)
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
|
*----
* TEST CASE TCM49
* NSE 2D CARTESIAN PROBLEM MODIFIED FOR 3D SN REGRESSION TEST
* ONE-GROUP MACROLIB-DEFINED CROSS SECTIONS
*
* Ref: A. Hebert, "The Search for Superconvergence in Spherical
* Harmonics Approximations", Nucl. Sci. Eng, 2006.
*----
* Define STRUCTURES and MODULES used
*----
LINKED_LIST GEOM MACRO TRACK SYSTEM FLUX EDIT ;
MODULE GEO: MAC: SNT: ASM: FLU: DELETE: END: ;
INTEGER s := 2 ;
PROCEDURE assertS ;
*
*--------
*GEOMETRY
*--------
GEOM := GEO: :: CAR3D 5 5 1
X- REFL X+ VOID
Y- REFL Y+ VOID
Z- VOID Z+ VOID
MIX 1 1 2 2 3
1 1 2 2 3
2 2 2 2 3
2 2 2 3 3
3 3 3 3 0
MESHX 0.0 40.0 80.0 120.0 160.0 200.0
MESHY 0.0 40.0 80.0 120.0 160.0 200.0
MESHZ 0.0 1000.0
SPLITX <<s>> <<s>> <<s>> <<s>> <<s>>
SPLITY <<s>> <<s>> <<s>> <<s>> <<s>>
;
*--------
*MACROLIB
*--------
MACRO := MAC: ::
EDIT 2 NGRO 1 NMIX 3 NIFI 1 ANIS 2
READ INPUT
MIX 1
TOTAL 0.025
NUSIGF 0.0155
CHI 1.0
SCAT 1 1 0.013
1 1 0.0
MIX 2
TOTAL 0.025
SCAT 1 1 0.024
1 1 0.006
MIX 3
TOTAL 0.075
SCAT 1 1 0.0
1 1 0.0
;
*--------
*TRACKING
*--------
*----
* HIGH-ORDER DIAMOND DIFFERENCE
*----
!
!FLAT ORDER
TRACK := SNT: GEOM ::
TITLE '3D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
EDIT 2 MAXR 5000 DIAM 0 SN 6 SCAT 2 QUAD 1 DSA 10 0 2 ;
SYSTEM := ASM: MACRO TRACK ::
EDIT 1 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM ::
EDIT 2 TYPE K EXTE 100 5.0E-8 ;
assertS FLUX :: 'K-EFFECTIVE' 1 0.9178295 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
!
!LINEAR ORDER
TRACK := SNT: GEOM ::
TITLE '3D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
EDIT 2 MAXR 5000 DIAM 1 SN 6 SCAT 2 QUAD 1 DSA 10 1 2 ;
SYSTEM := ASM: MACRO TRACK ::
EDIT 1 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM ::
EDIT 2 TYPE K EXTE 100 5.0E-8 ;
assertS FLUX :: 'K-EFFECTIVE' 1 0.9817208 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
!
!PARABOLIC ORDER
TRACK := SNT: GEOM ::
TITLE '3D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
EDIT 2 MAXR 5000 DIAM 2 SN 6 SCAT 2 QUAD 1 DSA 10 2 2 ;
SYSTEM := ASM: MACRO TRACK ::
EDIT 1 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM ::
EDIT 2 TYPE K EXTE 100 5.0E-8 ;
assertS FLUX :: 'K-EFFECTIVE' 1 0.9829153 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
*----
* DISCONTINUOUS GALERKIN
*----
!
!LINEAR ORDER
TRACK := SNT: GEOM ::
TITLE '3D ANISOTROPIC DG LINEAR BENCHMARK, 1 GROUP'
EDIT 2 MAXR 5000 SCHM 2 DIAM 1 SN 6 SCAT 2 QUAD 1 DSA 10 1 2 ;
SYSTEM := ASM: MACRO TRACK ::
EDIT 1 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM ::
EDIT 2 TYPE K EXTE 100 5.0E-8 ;
assertS FLUX :: 'K-EFFECTIVE' 1 0.9509456 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
!
!PARABOLIC ORDER
TRACK := SNT: GEOM ::
TITLE '3D ANISOTROPIC DG PARABOLIC BENCHMARK, 1 GROUP'
EDIT 2 MAXR 5000 SCHM 2 DIAM 2 SN 6 SCAT 2 QUAD 1 DSA 10 2 2 ;
SYSTEM := ASM: MACRO TRACK ::
EDIT 1 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM ::
EDIT 2 TYPE K EXTE 100 5.0E-8 ;
assertS FLUX :: 'K-EFFECTIVE' 1 0.9830175 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
!
!CUBIC ORDER
TRACK := SNT: GEOM ::
TITLE '3D ANISOTROPIC DG CUBIC BENCHMARK, 1 GROUP'
EDIT 2 MAXR 5000 SCHM 2 DIAM 3 SN 4 SCAT 2 QUAD 1 ;
SYSTEM := ASM: MACRO TRACK ::
EDIT 1 ARM ;
FLUX := FLU: MACRO TRACK SYSTEM ::
EDIT 2 TYPE K EXTE 100 5.0E-8 ;
assertS FLUX :: 'K-EFFECTIVE' 1 0.9827218 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
ECHO "test TCM49 completed" ;
END: ;
QUIT "LIST" .
|
