/
composition_active_with_melt.prm
119 lines (91 loc) · 2.97 KB
/
composition_active_with_melt.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
114
115
116
117
118
119
#########################################################
# This is a variation of the platelike-boundary.prm file.
# It is identical to the composition-active test, except
# for the second compositional field, which is called
# porosity and set to zero, and Include melt transport is
# set to true. Result can be compared to the test
# composition-active-without-melt.
set Dimension = 2
set Start time = 0
set End time = 1
set Use years in output instead of seconds = false
subsection Discretization
set Stokes velocity polynomial degree = 2
set Composition polynomial degree = 1
end
subsection Geometry model
set Model name = box
subsection Box
set X extent = 2
set Y extent = 1
end
end
subsection Model settings
set Fixed temperature boundary indicators = 2, 3
set Zero velocity boundary indicators =
set Tangential velocity boundary indicators = 0, 1, 2
set Prescribed velocity boundary indicators = 3: function
set Include melt transport = true
end
subsection Melt settings
set Melt transport threshold = 0.0
set Use discontinuous compaction pressure = false
end
subsection Heating model
set List of model names = shear heating
end
subsection Boundary temperature model
set List of model names = box
subsection Box
set Bottom temperature = 1
set Top temperature = 0
end
end
subsection Boundary velocity model
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
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
# Compared to the passive material model, we here make
# the density composition dependent by letting it depend
# linearly on the value of the first compositional field.
subsection Material model
set Model name = simple with melt
end
subsection Mesh refinement
set Initial adaptive refinement = 0
set Initial global refinement = 3
set Time steps between mesh refinement = 0
end
subsection Postprocess
set List of postprocessors = temperature statistics, composition statistics
end
# This is the new part: We declare that there will
# be two compositional fields that will be
# advected along. Their initial conditions are given by
# a function that is one for the lowermost 0.2 height
# units of the domain and zero otherwise in the first case,
# and one in the top most 0.2 height units in the latter.
subsection Compositional fields
set Number of fields = 2
set Names of fields = C_1, porosity
end
subsection Initial composition model
set Model name = function
subsection Function
set Variable names = x,y
set Function expression = if(y<0.2, 1, 0) ; 0.0
end
end