-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathoriginal.prm
174 lines (144 loc) · 5.37 KB
/
original.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
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
set Dimension = 2
set Output directory = outputs/MR_RG_VK_case9_AGG
set Use years in output instead of seconds = true
subsection Discretization
set Stokes velocity polynomial degree = 2
set Use locally conservative discretization = true
end
subsection Formulation
set Formulation = Boussinesq approximation
end
# Load the prescribed velocity and pressure library
# set Additional shared libraries = ./libprescribed_velocity_pressure.so
set Additional shared libraries = /home/staff/agrima/aspect/release_210316/prescribed_velocity_pressure_plugin/libprescribed_velocity_pressure.so
subsection Geometry model
set Model name = box
subsection Box
set X extent = 6.6e5
set Y extent = 6.0e5
set X repetitions = 66
set Y repetitions = 60
end
end
set Pressure normalization = no
set Surface pressure = 0
set Use conduction timestep = true
set End time = 2.4e8
set Resume computation = false
subsection Termination criteria
set End step = 1
set Termination criteria = end time
set Checkpoint on termination = true
end
subsection Checkpointing
set Steps between checkpoint = 200
end
subsection Initial temperature model
set Model name = function
subsection Function
set Variable names = x,y,t
set Function constants = H=6.0e5, Hl=5.0e4, TS=273.0, T0=1573.0, kappa=0.7272e-6, t50=1.5778463e15
set Function expression = if(x<=(H-y),TS+(T0-TS)*(1.0-erfc( (H-y)/(2.0*sqrt(kappa*t50))) ), \
if(x>(H-y) & y>(H-Hl), TS+(T0-TS)*(H-y)/Hl, \
T0 ) )
end
end
# The parameters below this comment were created by the update script
# as replacement for the old 'Model settings' subsection. They can be
# safely merged with any existing subsections with the same name.
subsection Boundary velocity model
set Prescribed velocity boundary indicators = left: function, top:function
end
subsection Boundary traction model
set Prescribed traction boundary indicators = right: zero traction, bottom:zero traction
end
subsection Boundary temperature model
set Fixed temperature boundary indicators = top, left, right
end
subsection Boundary temperature model
set List of model names = initial temperature
end
subsection Gravity model
set Model name = vertical
subsection Vertical
set Magnitude = 0.0
end
end
subsection Material model
set Model name = simple
subsection Simple model
set Viscosity = 1e21
set Thermal conductivity = 3.0
set Reference density = 3300.0
set Reference specific heat = 1250.0
set Thermal expansion coefficient = 0.0
end
end
subsection Mesh refinement
set Initial global refinement = 1
set Initial adaptive refinement = 0
set Time steps between mesh refinement = 0
set Run postprocessors on initial refinement = true
set Refinement fraction = 1.0
set Coarsening fraction = 0.0
set Strategy = minimum refinement function
subsection Minimum refinement function
set Function expression = 7.0
end
set Additional refinement times = 2.0e8, 2.2e8
end
subsection Boundary velocity model
subsection Function
set Variable names = x,y,t
set Function constants = H=6.0e5, v=0.035355339, eps=0.01, width=50000
set Function expression = if( x<=(H-y)+eps, v, if(x<=(H-y)+width*sqrt(2), v*(0.5 * (1 + cos(pi*(x-(H-y))/(width*sqrt(2))))), 0.0)); \
if( x<=(H-y)+eps,-v, if(x<=(H-y)+width*sqrt(2),-v*(0.5 * (1 + cos(pi*(x-(H-y))/(width*sqrt(2))))), 0.0))
end
end
# Turn prescribed velocities and pressure on
set Prescribe internal velocities = true
set Prescribe internal pressure = true
subsection Prescribed fields
set Use x velocity indicator for pressure indicator = true
set Use x velocity indicator for pressure indicator = true
subsection Velocity indicator function
set Variable names = x,y,t
# Return where to prescribe u_x; u_y; u_z
# (last one only used if dimension = 3)
# 1 if velocity should be prescribed, 0 otherwise
set Function constants = pi=3.14159, H=6.0e5, W=6.6e5, eps=0.01
set Function expression = if( x<=(H-y)+eps | ( y>=5.5e5-eps ), 1, 0); \
if( x<=(H-y)+eps | ( y>=5.5e5-eps ), 1, 0)
end
subsection Velocity function
set Variable names = x,y,t
# Return u_x; u_y; u_z (u_z only used if in 3d)
set Function constants = H=6.0e5, v=0.035355339, eps=0.01, width=50000
set Function expression = if( x<=(H-y)+eps, v, if(x<=(H-y)+width*sqrt(2), v*(0.5 * (1 + cos(pi*(x-(H-y))/(width*sqrt(2))))), 0.0)); \
if( x<=(H-y)+eps,-v, if(x<=(H-y)+width*sqrt(2),-v*(0.5 * (1 + cos(pi*(x-(H-y))/(width*sqrt(2))))), 0.0))
end
subsection Pressure indicator function
set Function constants = pi=3.14159, H=6.0e5, W=6.6e5, eps=0.01
set Function expression = if( x<=(H-y)+eps | ( y>=5.5e5-eps ), 1, 0)
end
subsection Pressure function
set Function expression = 0.0
end
end
subsection Postprocess
set List of postprocessors = visualization, temperature statistics, velocity statistics
subsection Visualization
set List of output variables = volumetric strain rate
# set Output format = vtu
set Output format = hdf5
set Interpolate output = false
set Time between graphical output = 0
end
end
subsection Solver parameters
subsection Stokes solver parameters
set Maximum number of expensive Stokes solver steps = 10000
set Linear solver tolerance = 1.0e-7
set Use direct solver for Stokes system = false
end
end