-
Notifications
You must be signed in to change notification settings - Fork 1
/
jaguar2d.i
99 lines (89 loc) · 1.7 KB
/
jaguar2d.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
[Mesh]
type = GeneratedMesh
dim = 2
nx = 625
ny = 125
xmax = 0.005 # Length of test block
ymax = 0.001 # Width of test block
[]
[Variables]
[./temperature]
initial_condition = 296 # Initial temperature
[../]
[]
[Functions]
# [./location_x]
# type = ParsedFunction
# value = '0.0 + t'
# # vars = 'alpha'
# # vals = '0.05'
# [../]
# [./location_y]
# type = ParsedFunction
# value = '0.25'
# # vars = 'alpha'
# # vals = '0.05'
# [../]
# [./location_z]
# type = ParsedFunction
# value = '0.0'
# # vars = 'alpha'
# # vals = '0.05'
# [../]
[]
[Kernels]
[./heat_conduction]
type = HeatConduction
variable = temperature
[../]
[./heat_conduction_time_derivative]
type = HeatConductionTimeDerivative
variable = temperature
[../]
[]
[DiracKernels]
[./point_heat_source]
type = MovingDirac2d
variable = temperature
value = 3560
point = '0 0 0'
[../]
[]
[BCs]
# [./inlet_temperature]
# type = DirichletBC
# variable = temperature
# boundary = top
# value = 600 # (K)
# function = heatsource_func
# [../]
[./outlet_temperature]
type = DirichletBC
variable = temperature
boundary = bottom
value = 296 # (K)
[../]
[]
[Materials]
[./Ti64]
type = GenericConstantMaterial
prop_names = 'thermal_conductivity specific_heat density'
prop_values = '6.7 526 4430' # W/m*K, J/kg-K, kg/m^3 @ 296K
[../]
[]
[Problem]
type = FEProblem
coord_type = RZ
rz_coord_axis = X
[]
[Executioner]
type = Transient
num_steps = 100
end_time = 0.005
solve_type = 'PJFNK'
petsc_options_iname = '-pc_type -pc_hypre_type'
petsc_options_value = 'hypre boomeramg'
[]
[Outputs]
exodus = true
[]