/
composition-reaction.prm
113 lines (92 loc) · 3.17 KB
/
composition-reaction.prm
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
# A description of convection in a rectangular box where we prescribe
# temporally variable boundary conditions for the velocity along the top
# boundary. See the manual for more information.
############### Global parameters
set Dimension = 2
set Start time = 0
set End time = 20
set Use years in output instead of seconds = false
set Output directory = output-composition-reaction
############### Parameters describing the model
# Let us here choose again a box domain of size 2x1
# where we fix the temperature at the bottom and top,
# allow free slip along the bottom, left and right,
# and prescribe the velocity along the top using the
# `function' description.
subsection Geometry model
set Model name = box
subsection Box
set X extent = 2
set Y extent = 1
end
end
subsection Compositional fields
set Number of fields = 2
end
# We then set the temperature to one at the bottom and zero
# at the top:
subsection Boundary temperature model
set Fixed temperature boundary indicators = bottom, top
set List of model names = box
subsection Box
set Bottom temperature = 1
set Top temperature = 0
end
end
# The velocity along the top boundary models a spreading
# center that is moving left and right:
subsection Boundary velocity model
set Tangential velocity boundary indicators = left, right, bottom
set Prescribed velocity boundary indicators = top: function
subsection Function
set Variable names = x,z,t
set Function constants = pi=3.1415926
set Function expression = if(x>1+sin(0.5*pi*t), 1, -1); 0
end
end
# We then choose a vertical gravity model and describe the
# initial temperature with a vertical gradient. The default
# strength for gravity is one. The material model is the
# same as before.
subsection Gravity model
set Model name = vertical
end
subsection Initial temperature model
set Model name = function
subsection Function
set Variable names = x,z
set Function expression = (1-z)
end
end
subsection Initial composition model
set Model name = function
subsection Function
set Variable names = x,z
set Function expression = if(z<0.5,1,0);0
end
end
subsection Material model
set Model name = composition reaction
subsection Composition reaction model
set Thermal conductivity = 1e-6
set Thermal expansion coefficient = 1e-4
set Viscosity = 1
set Density differential for compositional field 1 = -5
set Density differential for compositional field 2 = 5
set Reaction depth = 0.2
end
end
# The final part of this input file describes how many times the
# mesh is refined and what to do with the solution once computed
subsection Mesh refinement
set Initial adaptive refinement = 0
set Initial global refinement = 5
set Time steps between mesh refinement = 0
end
subsection Postprocess
set List of postprocessors = visualization, temperature statistics, heat flux statistics
subsection Visualization
set Time between graphical output = 0.1
set List of output variables = density
end
end