Merge pull request #106 from pmgbergen/fractured_mpsa

Fractured mpsa

Author: rbe051

examples/papers/arXiv-1712_00460/low_pressure_stimulation/Simultaion.ipynb

@@ -0,0 +1,8 @@
+1200, 2200, 2000, 1500, 1500,  0, 100, 80, 25
+1800, 1800, 2150,  700,  700,  30, 170, 70, 13
+500, 2500, 1850,  850,  850,  0, 130, 70, 12
+2100, 2100, 2000, 1000, 1000,  0,  10, 80, 16
+2000, 2400, 1500,  900,  800,  5, 110, 80, 15
+1500, 2500, 1800,  500,  800,  0,  20, 90, 12
+2000, 3000, 2500,  700,  500, 30,  50, 80, 10
+

examples/papers/arXiv-1712_00460/low_pressure_stimulation/fig/aperture_change.gif

@@ -48,9 +48,10 @@ def __init__(self, deltaT):
             def matrix_rhs(self, g, data):
                 lhs, rhs = mass_matrix.MassMatrix.matrix_rhs(self, g, data)
                 return lhs * data['compressibility'], rhs * data['compressibility']
-        single_dim_discr = TimeDisc(self.time_step())
-        multi_dim_discr = Coupler(single_dim_discr)
-        return multi_dim_discr
+        time_discretization = TimeDisc(self.time_step())
+        if self.is_GridBucket:
+            time_discretization = Coupler(time_discretization)
+        return time_discretization
 
     def pressure(self, pressure_name='pressure'):
         if self.is_GridBucket:
@@ -62,6 +63,13 @@ def discharge(self, d_name='discharge', p_name='pressure'):
         self.pressure(p_name)
         fvutils.compute_discharges(self.grid(),d_name=d_name, p_name=p_name)
 
+    def pressure(self, pressure_name='pressure'):
+        self.pressure_name = pressure_name
+        if self.is_GridBucket:
+            self.split(self.pressure_name)
+        else:
+            self._data[self.pressure_name] = self._solver.p
+
 
 class SlightlyCompressibleDataAssigner(ParabolicDataAssigner):
     '''