-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathgeomech.in
158 lines (147 loc) · 2.09 KB
/
geomech.in
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
SIMULATION
SIMULATION_TYPE GEOMECHANICS_SUBSURFACE
PROCESS_MODELS
SUBSURFACE_FLOW
MODE RICHARDS
/
GEOMECHANICS_SUBSURFACE geom
/
END
#---
SUBSURFACE
#--- Spatial discretization
GRID
TYPE STRUCTURED
NXYZ 15 15 3
DXYZ
1524
1524
10 10 10
/
END
#--- Layers property
MATERIAL_PROPERTY layer1
ID 1
POROSITY 0.15
CHARACTERISTIC_CURVES satu
PERMEABILITY
PERM_HORIZONTAL 9.07e-11
VERTICAL_ANISOTROPY_RATIO 1
/
END
MATERIAL_PROPERTY layer2
ID 2
POROSITY 0.15
CHARACTERISTIC_CURVES satu
PERMEABILITY
PERM_HORIZONTAL 9.07e-16
VERTICAL_ANISOTROPY_RATIO 1
/
END
MATERIAL_PROPERTY layer3
ID 3
POROSITY 0.15
CHARACTERISTIC_CURVES satu
PERMEABILITY
PERM_HORIZONTAL 9.07e-12
VERTICAL_ANISOTROPY_RATIO 1
/
END
#--- Characteristic curves
CHARACTERISTIC_CURVES satu
SATURATION_FUNCTION VAN_GENUCHTEN
M 1.d-30
ALPHA 1.9401d-4
LIQUID_RESIDUAL_SATURATION 0.115
/
PERMEABILITY_FUNCTION MUALEM_VG_LIQ
M 1.d-30
LIQUID_RESIDUAL_SATURATION 0.115
/
END
#--- Regions
REGION ly1
COORDINATES
0 0 0
22860 22860 10
/
FACE TOP
END
REGION ly2
COORDINATES
0 0 10
22860 22860 20
/
END
REGION ly3
COORDINATES
0 0 20
22860 22860 30
/
END
REGION all
COORDINATES
0 0 0
22860 22860 30
/
END
#--- Strata
STRATA
MATERIAL layer1
REGION ly1
END
STRATA
MATERIAL layer2
REGION ly2
END
STRATA
MATERIAL layer3
REGION ly3
END
#--- Flow condition
FLOW_CONDITION rch
TYPE
RATE VOLUMETRIC_RATE NEIGHBOR_PERM
/
RATE 1 m^3/day
# RATE 0.d0 1.d-5 0.d0 kg/s kg/s MW
END
FLOW_CONDITION chd
TYPE
LIQUID_PRESSURE DIRICHLET
/
LIQUID_PRESSURE 0.
END
#--- Boundary conditions
#BOUNDARY_CONDITION
# FLOW_CONDITION rch
# REGION ly1
#END
SOURCE_SINK
FLOW_CONDITION rch
REGION ly1
END
#--- INITIAL CONDITION
INITIAL_CONDITION
FLOW_CONDITION chd
REGION all
END
#--- Time
TIME
FINAL_TIME 24 h
INITIAL_TIMESTEP_SIZE 0.01 h
MAXIMUM_TIMESTEP_SIZE 5.d-2 y
END
#--- Output
OUTPUT
TIMES y 0.01 0.1 1.0
FORMAT TECPLOT POINT
END
END_SUBSURFACE
#---Geomechanics
GEOMECHANICS
GEOMECHANICS_GRID
TYPE unstructured mesh_for_pflotran.umi
GRAVITY 0.0 0.0 0.0
END
END_GEOMECHANICS