-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathparameters_biot.prm
executable file
·153 lines (135 loc) · 9.07 KB
/
parameters_biot.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
# Listing of Parameters
# ---------------------
#################################################################################
# Discretization
set degree = 1
set refinements = 3
set grid_flag = 0
set time_step = 0.1
set num_time_steps = 3
#################################################################################
# Physical parameters
set alpha = 1.0
set Storativity = 1.0
subsection permeability 2D
set Function constants =
# Order: 00,01,10,11
set Function expression = 1;0;0;1
set Variable names = x,y,t
end
subsection permeability 3D
set Function constants =
# Order: 00,01,10,11
set Function expression = 1; 0; 0; 0; 1; 0; 0; 0; 1
set Variable names = x,y,z,t
end
subsection lambda 2D
set Function constants =
set Function expression = 100.0
set Variable names = x,y,t
end
subsection lambda 3D
set Function constants =
set Function expression = 123.0
set Variable names = x,y,z,t
end
subsection mu 2D
set Function constants =
set Function expression = 100.0
set Variable names = x,y,t
end
subsection mu 3D
set Function constants =
set Function expression = 79.3
set Variable names = x,y,z,t
end
#################################################################################
# Right hand sides
# Darcy RHS
subsection Darcy RHS 2D
set Function constants =
set Function expression = exp(t)*(cos(y*pi)*sin(x*pi) + 10) + exp(t)*(2*x + 3*x^2*y^4 + cos(y - 1)*cos((x - 1)*(y - 1))*(y - 1)) - exp(t)*(3*(x - 1)^4*(y - 1)^2 - 2*y + x*sin(x*y)*sin(x) + 2) + 2*pi^2*exp(t)*cos(y*pi)*sin(x*pi)
set Variable names = x,y,t
end
subsection Darcy RHS 3D
set Function constants =
set Function expression = x+y+exp(t)*sin(x+y)*2.0+1.0*exp(t)*sin(x+y)
set Variable names = x,y,z,t
end
# Elasticity RHS
subsection Elasticity RHS 2D
set Function constants =
set Function expression = 100*exp(t)*(cos(y - 1)*sin((x - 1)*(y - 1)) - 12*x^3*y^2 + cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)^2 + 2*sin(y - 1)*cos((x - 1)*(y - 1))*(x - 1)) - 300*exp(t)*(6*x*y^4 - cos(y - 1)*sin((x - 1)*(y - 1))*(y - 1)^2 + 2) + 200*exp(t)*(sin(x*y)*sin(x) + 12*(x - 1)^3*(y - 1)^2 + x*sin(x*y)*cos(x) + x*y*cos(x*y)*sin(x)) + pi*exp(t)*cos(x*pi)*cos(y*pi); 300*exp(t)*(3*(2*y - 2)*(x - 1)^4 + x^2*cos(x*y)*sin(x) - 2) - 200*exp(t)*(cos(y - 1)*cos((x - 1)*(y - 1)) + 12*x^2*y^3 - sin(y - 1)*cos((x - 1)*(y - 1))*(y - 1) - cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)*(y - 1)) + 100*exp(t)*(12*(x - 1)^2*(y - 1)^3 + cos(x*y)*sin(x) + 2*y*sin(x*y)*cos(x) + y^2*cos(x*y)*sin(x)) - pi*exp(t)*sin(x*pi)*sin(y*pi)
set Variable names = x,y,t
end
subsection Elasticity RHS 3D
set Function constants =
set Function expression = 123.0*t+79.3*t*2.0-exp(t)*cos(x+y); 123.0*t+79.3*t*2.0-exp(t)*cos(x+y)
set Variable names = x,y,z,t
end
#################################################################################
# Boundary conditions
# Darcy BC
subsection Darcy BC 2D
set Function constants =
set Function expression = exp(t)*(cos(y*pi)*sin(x*pi) + 10)
set Variable names = x,y,t
end
subsection Darcy BC 3D
set Function constants =
set Function expression = exp(t)*sin(x+y)
set Variable names = x,y,z,t
end
# Elasticity BC
subsection Elasticity BC 2D
set Function constants =
set Function expression = exp(t)*(x^3*y^4 + cos(y - 1)*sin((x - 1)*(y - 1)) + x^2); exp(t)*((y - 1)^2 - (x - 1)^4*(y - 1)^3 + cos(x*y)*sin(x))
set Variable names = x,y,t
end
subsection Elasticity BC 3D
set Function constants =
set Function expression = t*0.5*x^2; t*0.5*y^2
set Variable names = x,y,z,t
end
#################################################################################
# Initial conditions (only pressure and stress values are used)
subsection Initial conditions 2D
set Function constants =
# Order: velocity, pressure, stress, displacement, rotation
set Function expression = 0; 0; cos(y*pi)*sin(x*pi) + 10; 600*x + 200*y + 900*x^2*y^4 - cos(y*pi)*sin(x*pi) - 300*(x - 1)^4*(y - 1)^2 - 100*x*sin(x*y)*sin(x) + 300*cos(y - 1)*cos((x - 1)*(y - 1))*(y - 1) - 210; 400*x^3*y^3 - 100*sin(y - 1)*sin((x - 1)*(y - 1)) - 400*(x - 1)^3*(y - 1)^3 + 100*cos(x*y)*cos(x) - 100*y*sin(x*y)*sin(x) + 100*cos(y - 1)*cos((x - 1)*(y - 1))*(x - 1); 400*x^3*y^3 - 100*sin(y - 1)*sin((x - 1)*(y - 1)) - 400*(x - 1)^3*(y - 1)^3 + 100*cos(x*y)*cos(x) - 100*y*sin(x*y)*sin(x) + 100*cos(y - 1)*cos((x - 1)*(y - 1))*(x - 1); 200*x + 600*y + 300*x^2*y^4 - cos(y*pi)*sin(x*pi) - 900*(x - 1)^4*(y - 1)^2 - 300*x*sin(x*y)*sin(x) + 100*cos(y - 1)*cos((x - 1)*(y - 1))*(y - 1) - 610; 0; 0; 0
set Variable names = x,y,t
end
subsection Initial conditions 3D
set Function constants =
# Order: velocity, pressure, stress, displacement, rotation
set Function expression = 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0
set Variable names = x,y,z,t
end
#################################################################################
# Exact solution, values
subsection Exact solution 2D
set Function constants =
# Order: velocity, pressure, stress, displacement, rotation
set Function expression = -pi*exp(t)*cos(x*pi)*cos(y*pi); pi*exp(t)*sin(x*pi)*sin(y*pi); exp(t)*(cos(y*pi)*sin(x*pi) + 10); 300*exp(t)*(2*x + 3*x^2*y^4 + cos(y - 1)*cos((x - 1)*(y - 1))*(y - 1)) - exp(t)*(cos(y*pi)*sin(x*pi) + 10) - 100*exp(t)*(3*(x - 1)^4*(y - 1)^2 - 2*y + x*sin(x*y)*sin(x) + 2); 100*exp(t)*(4*x^3*y^3 - sin(y - 1)*sin((x - 1)*(y - 1)) + cos(y - 1)*cos((x - 1)*(y - 1))*(x - 1)) - 100*exp(t)*(4*(x - 1)^3*(y - 1)^3 - cos(x*y)*cos(x) + y*sin(x*y)*sin(x)); 100*exp(t)*(4*x^3*y^3 - sin(y - 1)*sin((x - 1)*(y - 1)) + cos(y - 1)*cos((x - 1)*(y - 1))*(x - 1)) - 100*exp(t)*(4*(x - 1)^3*(y - 1)^3 - cos(x*y)*cos(x) + y*sin(x*y)*sin(x)); 100*exp(t)*(2*x + 3*x^2*y^4 + cos(y - 1)*cos((x - 1)*(y - 1))*(y - 1)) - exp(t)*(cos(y*pi)*sin(x*pi) + 10) - 300*exp(t)*(3*(x - 1)^4*(y - 1)^2 - 2*y + x*sin(x*y)*sin(x) + 2); exp(t)*(x^3*y^4 + cos(y - 1)*sin((x - 1)*(y - 1)) + x^2); exp(t)*((y - 1)^2 - (x - 1)^4*(y - 1)^3 + cos(x*y)*sin(x)); (exp(t)*(4*x^3*y^3 - sin(y - 1)*sin((x - 1)*(y - 1)) + cos(y - 1)*cos((x - 1)*(y - 1))*(x - 1)))/2 + (exp(t)*(4*(x - 1)^3*(y - 1)^3 - cos(x*y)*cos(x) + y*sin(x*y)*sin(x)))/2
set Variable names = x,y,t
end
subsection Exact solution 3D
set Function constants =
# Order: velocity, pressure, stress, displacement, rotation
set Function expression = 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0
set Variable names = x,y,z,t
end
#################################################################################
# Exact solution, gradient values (only velocity and stress values are used)
subsection Exact gradient 2D
set Function constants =
# Order: velocity, pressure, stress, displacement, rotation
set Function expression = pi^2*exp(t)*cos(y*pi)*sin(x*pi); pi^2*exp(t)*cos(x*pi)*sin(y*pi); pi^2*exp(t)*cos(x*pi)*sin(y*pi); pi^2*exp(t)*cos(y*pi)*sin(x*pi); 0; 0; 300*exp(t)*(6*x*y^4 - cos(y - 1)*sin((x - 1)*(y - 1))*(y - 1)^2 + 2) - 100*exp(t)*(sin(x*y)*sin(x) + 12*(x - 1)^3*(y - 1)^2 + x*sin(x*y)*cos(x) + x*y*cos(x*y)*sin(x)) - pi*exp(t)*cos(x*pi)*cos(y*pi); 300*exp(t)*(cos(y - 1)*cos((x - 1)*(y - 1)) + 12*x^2*y^3 - sin(y - 1)*cos((x - 1)*(y - 1))*(y - 1) - cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)*(y - 1)) - 100*exp(t)*(3*(2*y - 2)*(x - 1)^4 + x^2*cos(x*y)*sin(x) - 2) + pi*exp(t)*sin(x*pi)*sin(y*pi); 100*exp(t)*(cos(y - 1)*cos((x - 1)*(y - 1)) + 12*x^2*y^3 - sin(y - 1)*cos((x - 1)*(y - 1))*(y - 1) - cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)*(y - 1)) - 100*exp(t)*(12*(x - 1)^2*(y - 1)^3 + cos(x*y)*sin(x) + 2*y*sin(x*y)*cos(x) + y^2*cos(x*y)*sin(x)); - 100*exp(t)*(cos(y - 1)*sin((x - 1)*(y - 1)) - 12*x^3*y^2 + cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)^2 + 2*sin(y - 1)*cos((x - 1)*(y - 1))*(x - 1)) - 100*exp(t)*(sin(x*y)*sin(x) + 12*(x - 1)^3*(y - 1)^2 + x*sin(x*y)*cos(x) + x*y*cos(x*y)*sin(x)); 100*exp(t)*(cos(y - 1)*cos((x - 1)*(y - 1)) + 12*x^2*y^3 - sin(y - 1)*cos((x - 1)*(y - 1))*(y - 1) - cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)*(y - 1)) - 100*exp(t)*(12*(x - 1)^2*(y - 1)^3 + cos(x*y)*sin(x) + 2*y*sin(x*y)*cos(x) + y^2*cos(x*y)*sin(x)); - 100*exp(t)*(cos(y - 1)*sin((x - 1)*(y - 1)) - 12*x^3*y^2 + cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)^2 + 2*sin(y - 1)*cos((x - 1)*(y - 1))*(x - 1)) - 100*exp(t)*(sin(x*y)*sin(x) + 12*(x - 1)^3*(y - 1)^2 + x*sin(x*y)*cos(x) + x*y*cos(x*y)*sin(x)); 100*exp(t)*(6*x*y^4 - cos(y - 1)*sin((x - 1)*(y - 1))*(y - 1)^2 + 2) - 300*exp(t)*(sin(x*y)*sin(x) + 12*(x - 1)^3*(y - 1)^2 + x*sin(x*y)*cos(x) + x*y*cos(x*y)*sin(x)) - pi*exp(t)*cos(x*pi)*cos(y*pi); 100*exp(t)*(cos(y - 1)*cos((x - 1)*(y - 1)) + 12*x^2*y^3 - sin(y - 1)*cos((x - 1)*(y - 1))*(y - 1) - cos(y - 1)*sin((x - 1)*(y - 1))*(x - 1)*(y - 1)) - 300*exp(t)*(3*(2*y - 2)*(x - 1)^4 + x^2*cos(x*y)*sin(x) - 2) + pi*exp(t)*sin(x*pi)*sin(y*pi); 0; 0; 0; 0; 0; 0
set Variable names = x,y,t
end
subsection Exact gradient 3D
set Function constants =
# Order: velocity, pressure, stress, displacement, rotation
set Function expression = 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0
set Variable names = x,y,z,t
end