forked from idaholab/moose
/
sub_pin_depletion_id.i
119 lines (111 loc) · 2.52 KB
/
sub_pin_depletion_id.i
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
[Mesh]
[pin1_base]
type = PolygonConcentricCircleMeshGenerator
num_sides = 4
num_sectors_per_side = '4 4 4 4'
background_intervals = 2
polygon_size = 0.63
polygon_size_style = 'apothem'
ring_radii = '0.2 0.4 0.5'
ring_intervals = '2 2 1'
preserve_volumes = on
flat_side_up = true
sector_id_name = 'sector_id'
ring_id_name = 'ring_id'
[]
[pin1_a]
type = SubdomainExtraElementIDGenerator
input = pin1_base
subdomains = '1 2 3 4 5'
extra_element_id_names = 'material_id'
extra_element_ids = '1 1 1 8 9'
[]
[pin1]
type = RenameBoundaryGenerator
input = pin1_a
old_boundary = '10002 15002 10004 15004 10001 15001 10003 15003'
new_boundary = 'left left right right top top bottom bottom'
[]
[pin2_base]
type = PolygonConcentricCircleMeshGenerator
num_sides = 4
num_sectors_per_side = '4 4 4 4'
background_intervals = 2
polygon_size = 0.63
polygon_size_style = 'apothem'
ring_radii = '0.15 0.3 0.4'
ring_intervals = '2 3 1'
preserve_volumes = on
flat_side_up = true
[]
[pin2_a]
type = SubdomainExtraElementIDGenerator
input = pin2_base
subdomains = '1 2 3 4 5'
extra_element_id_names = 'material_id'
extra_element_ids = '2 2 2 8 9'
[]
[pin2]
type = RenameBoundaryGenerator
input = pin2_a
old_boundary = '10002 15002 10004 15004 10001 15001 10003 15003'
new_boundary = 'left left right right top top bottom bottom'
[]
[assembly]
type = CartesianIDPatternedMeshGenerator
inputs = 'pin1 pin2'
pattern = ' 1 0;
0 1'
assign_type = 'cell'
id_name = 'pin_id'
[]
[depletion_id]
type = DepletionIDGenerator
input = 'assembly'
id_name = 'pin_id sector_id ring_id'
material_id_name = 'material_id'
exclude_id_name = 'material_id ring_id'
exclude_id_value = '8 9; 0'
[]
[]
[Executioner]
type = Steady
[]
[Problem]
solve = false
[]
[AuxVariables]
[sector_id]
family = MONOMIAL
order = CONSTANT
[]
[ring_id]
family = MONOMIAL
order = CONSTANT
[]
[depletion_id]
family = MONOMIAL
order = CONSTANT
[]
[]
[AuxKernels]
[set_sector_id]
type = ExtraElementIDAux
variable = sector_id
extra_id_name = sector_id
[]
[set_ring_id]
type = ExtraElementIDAux
variable = ring_id
extra_id_name = ring_id
[]
[set_depletion_id]
type = ExtraElementIDAux
variable = depletion_id
extra_id_name = depletion_id
[]
[]
[Outputs]
exodus = true
execute_on = timestep_end
[]