Merge pull request #23 from wenqing/new_BM

Added one for modelling the excavation within H(Richards)M coupled processes
ufz · Dec 21, 2018 · a1cbb1a · a1cbb1a
2 parents bfb7323 + 524ea07
commit a1cbb1a
Show file tree

Hide file tree

Showing 12 changed files with 6,731 additions and 0 deletions.
diff --git a/HM/excavation/BoreholeExcavation/borehole_excav.bc b/HM/excavation/BoreholeExcavation/borehole_excav.bc
@@ -0,0 +1,177 @@
+GINA - Boundary Condition
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  SURFACE upper
+ $DIS_TYPE
+  CONSTANT   5e6
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  SURFACE right
+ $DIS_TYPE
+  CONSTANT   5e6
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  SURFACE back
+ $DIS_TYPE
+  CONSTANT   5e6
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  SURFACE bottom
+ $DIS_TYPE
+  CONSTANT   5e6
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  SURFACE left
+ $DIS_TYPE
+  CONSTANT   5e6
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  SURFACE front
+ $DIS_TYPE
+  CONSTANT   5e6
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $DIS_TYPE
+  CONSTANT   1
+ $TIM_TYPE
+  CURVE 4
+ $EXCAVATION
+  1 1
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  DISPLACEMENT_X1
+ $GEO_TYPE
+  SURFACE left
+ $DIS_TYPE
+  CONSTANT   0
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  DISPLACEMENT_X1
+ $GEO_TYPE
+  SURFACE right
+ $DIS_TYPE
+  CONSTANT   0
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  DISPLACEMENT_Y1
+ $GEO_TYPE
+  SURFACE back
+ $DIS_TYPE
+  CONSTANT   0
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  DISPLACEMENT_Y1
+ $GEO_TYPE
+  SURFACE front
+ $DIS_TYPE
+  CONSTANT   0
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  DISPLACEMENT_Z1
+ $GEO_TYPE
+  SURFACE bottom
+ $DIS_TYPE
+  CONSTANT   0
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  DISPLACEMENT_Z1
+ $GEO_TYPE
+  SURFACE upper
+ $DIS_TYPE
+  CONSTANT   0
+#STOP
+
+
+
+
+
+
+
+
+
+
+
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  POLYLINE LEFT
+ $DIS_TYPE
+  CONSTANT   9.e5
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  POLYLINE RIGHT
+ $DIS_TYPE
+  CONSTANT   9.e5
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  POLYLINE TOP
+ $DIS_TYPE
+  CONSTANT   9.e5
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  POLYLINE BOTTOM
+ $DIS_TYPE
+  CONSTANT   9.e5
+#BOUNDARY_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+
+ $DIS_TYPE
+  CONSTANT   -2.9596E+07
diff --git a/HM/excavation/BoreholeExcavation/borehole_excav.gli b/HM/excavation/BoreholeExcavation/borehole_excav.gli
@@ -0,0 +1,157 @@
+#POINTS
+0 -.5 0 -0.5 $NAME POINT0
+1 .5 0 -.50 $NAME POINT1
+2 .5 0 .5 $NAME POINT2
+3 -.5 0 .5 $NAME POINT3
+4 -0.5 1.5 -.50 $NAME POINT4
+5 .5 1.5 -.50 $NAME POINT5
+6 .5 1.5 .5 $NAME POINT6
+7 -.5 1.5 .5 $NAME POINT7
+8 0 0 .0429 $NAME POINT8
+9 0 10 .0429 $NAME POINT9
+10 0 0.5 0 $NAME POINT10
+11 0.5 0.5 0 $NAME POINT11
+12 0 0.5 0.5 $NAME POINT12
+#POLYLINE
+ $NAME
+  front
+ $POINTS
+  0
+  1
+  2
+  3
+  0
+ $MAT_GROUP
+  0
+#POLYLINE
+ $NAME
+  left
+ $POINTS
+  0
+  3
+  7
+  4
+  0
+ $MAT_GROUP
+  0
+#POLYLINE
+ $NAME
+  upper
+ $POINTS
+  3
+  2
+  6
+  7
+  3
+ $MAT_GROUP
+  0
+#POLYLINE
+ $NAME
+  bottom
+ $POINTS
+  0
+  1
+  5
+  4
+  0
+ $MAT_GROUP
+  0
+#POLYLINE
+ $NAME
+  back
+ $POINTS
+  4
+  5
+  6
+  7
+  4
+ $MAT_GROUP
+  0
+#POLYLINE
+ $NAME
+  right
+ $POINTS
+  2
+  1
+  5
+  6
+  2
+ $MAT_GROUP
+  0
+#POLYLINE
+ $NAME
+  horizon
+ $POINTS
+  10
+  11
+#POLYLINE
+ $NAME
+  vertikal
+ $POINTS
+  10
+  12
+#SURFACE
+ $NAME
+  left
+ $EPSILON
+  0.001
+ $POLYLINES
+  left
+ $MAT_GROUP
+  0
+ $TYPE
+  0
+#SURFACE
+ $NAME
+  upper
+ $EPSILON
+  0.001
+ $POLYLINES
+  upper
+ $MAT_GROUP
+  0
+ $TYPE
+  0
+#SURFACE
+ $NAME
+  bottom
+ $EPSILON
+  0.001
+ $POLYLINES
+  bottom
+ $MAT_GROUP
+  0
+ $TYPE
+  0
+#SURFACE
+ $NAME
+  back
+ $EPSILON
+  0.001
+ $POLYLINES
+  back
+ $MAT_GROUP
+  0
+ $TYPE
+  0
+#SURFACE
+ $NAME
+  right
+ $EPSILON
+  0.001
+ $POLYLINES
+  right
+ $MAT_GROUP
+  0
+ $TYPE
+  0
+#SURFACE
+ $NAME
+  front
+ $EPSILON
+  0.001
+ $POLYLINES
+  front
+ $TYPE
+  0
+#STOP
diff --git a/HM/excavation/BoreholeExcavation/borehole_excav.ic b/HM/excavation/BoreholeExcavation/borehole_excav.ic
@@ -0,0 +1,41 @@
+GINA - Initial Condition
+#INITIAL_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  STRESS_XX
+ $GEO_TYPE
+  SUB_DOMAIN
+    2
+    0   -9e6
+    1   -9e6
+#INITIAL_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  STRESS_YY
+ $GEO_TYPE
+  SUB_DOMAIN
+    2
+    0   -9e6
+    1   -9e6
+#INITIAL_CONDITION
+ $PCS_TYPE
+  DEFORMATION
+ $PRIMARY_VARIABLE
+  STRESS_ZZ
+ $GEO_TYPE
+  SUB_DOMAIN
+    2
+    0   -9e6
+    1   -9e6
+#INITIAL_CONDITION
+ $PCS_TYPE
+  RICHARDS_FLOW
+ $PRIMARY_VARIABLE
+  PRESSURE1
+ $GEO_TYPE
+  DOMAIN 
+ $DIS_TYPE
+  CONSTANT   5e6
+#STOP
diff --git a/HM/excavation/BoreholeExcavation/borehole_excav.mfp b/HM/excavation/BoreholeExcavation/borehole_excav.mfp
@@ -0,0 +1,11 @@
+#FLUID_PROPERTIES
+ $FLUID_TYPE
+  LIQUID
+ $PCS_TYPE
+  PRESSURE1
+ $DENSITY
+  1   1e3
+ $VISCOSITY
+  1 1.000000e-003     
+ $NON_GRAVITY
+#STOP