blob: 259aaa5130cd90ce11cdea4745afa5597fbde033 (
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
|
*----
* TEST CASE TCM35
* NSE 2D CARTESIAN PROBLEM
* 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 ;
*
GEOM := GEO: :: CAR2D 5 5
X- REFL X+ VOID
Y- REFL Y+ 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
SPLITX <<s>> <<s>> <<s>> <<s>> <<s>>
SPLITY <<s>> <<s>> <<s>> <<s>> <<s>>
;
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
;
!Level-Symmetric quadrature
TRACK := SNT: GEOM ::
TITLE '2D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
EDIT 2 MAXR 2500 DIAM 1 SN 20 SCAT 2
LIVO 3 3 DSA 0 1 1 MAXI 100 EPSI 1.E-5 GMRES 10
QUAD 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.9923077 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
!Legendre-Chebyshev quadrature
TRACK := SNT: GEOM ::
TITLE '2D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
EDIT 2 MAXR 2500 DIAM 1 SN 44 SCAT 2
LIVO 3 3 DSA 0 1 1 MAXI 100 EPSI 1.E-5 GMRES 10
QUAD 5 ;
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.9940005 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
!Quadruple range quadrature
TRACK := SNT: GEOM ::
TITLE '2D ANISOTROPIC NSE BENCHMARK, 1 GROUP'
EDIT 2 MAXR 2500 DIAM 1 SN 74 SCAT 2
LIVO 3 3 DSA 0 1 1 MAXI 100 EPSI 1.E-5 GMRES 10
QUAD 6 ;
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.9898672 ;
TRACK SYSTEM FLUX := DELETE: TRACK SYSTEM FLUX ;
ECHO "test TCM35 completed" ;
END: ;
QUIT "LIST" .
|