-
Notifications
You must be signed in to change notification settings - Fork 1
/
THeat.i
118 lines (107 loc) · 2.01 KB
/
THeat.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
!include Parameters.i
[Mesh]
type = FileMesh
file = ../mesh/vac_oval_coil_solid_target_coarse.e
second_order = true
[]
[Variables]
[T]
family = LAGRANGE
order = FIRST
initial_condition = ${room_temperature}
[]
[]
[AuxVariables]
[P]
family = MONOMIAL
order = CONSTANT
[]
[]
[Kernels]
[HeatConduction]
type = HeatConduction
variable = T
[]
[TimeDerivative]
type = HeatConductionTimeDerivative
variable = T
[]
[HeatSource]
type = CoupledForce
variable = T
v = P
[]
[]
[Materials]
[copper]
type = GenericConstantMaterial
prop_names = 'thermal_conductivity specific_heat density'
prop_values = '${copper_tconductivity} ${copper_capacity} ${copper_density}'
block = 'coil target'
[]
[vacuum]
type = GenericConstantMaterial
prop_names = 'thermal_conductivity specific_heat density'
prop_values = '${vacuum_tconductivity} ${vacuum_capacity} ${vacuum_density}'
block = vacuum_region
[]
[]
[BCs]
[plane]
type = DirichletBC
variable = T
boundary = 'coil_in coil_out terminal_plane'
value = ${room_temperature}
[]
[]
[Postprocessors]
[P(total){W}]
type = ElementIntegralVariablePostprocessor
variable = P
block = target
[]
[P(Max){W.m-3}]
type = ElementExtremeValue
variable = P
block = target
[]
[T(average){K}]
type = ElementAverageValue
variable = T
block = target
[]
[T(Max){K}]
type = NodalExtremeValue
variable = T
block = target
[]
[]
[Executioner]
type = Transient
solve_type = LINEAR
petsc_options_iname = -pc_type
petsc_options_value = hypre
start_time = 0.0
end_time = ${end_t}
dt = ${delta_t}
[]
[Outputs]
exodus = true
csv = true
[]
[MultiApps]
[AForm]
type = TransientMultiApp
input_files = AForm.i
execute_on = timestep_begin
clone_parent_mesh = true
[]
[]
[Transfers]
[pull_power]
type = MultiAppCopyTransfer
from_multi_app = AForm
source_variable = P
variable = P
[]
[]