Added: virtual method for computing the total permeability tensor

Author: kmokstad

PoroFracture.C

@@ -12,6 +12,7 @@
 //==============================================================================
 
 #include "PoroFracture.h"
+#include "PoroMaterial.h"
 #include "FractureElasticityVoigt.h"
 #include "FiniteElement.h"
 #include "Tensor.h"
@@ -207,7 +208,11 @@ double PoroFracture::formCrackedPermeabilityTensor (SymmTensor& Kcrack,
   if (d < 0.0)
     return d; // Invalid crack phase field
   else if (d < d_min)
-    return 0.0; // The crack does not affect the permeability tensor
+  {
+    // The crack does not affect the permeability tensor at this point
+    Kcrack.zero();
+    return 0.0;
+  }
 
   double d2 = gradD.length2();
   if (d2 <= 0.0) return -1.0; // Should not happen
@@ -236,3 +241,22 @@ double PoroFracture::formCrackedPermeabilityTensor (SymmTensor& Kcrack,
   Kcrack *= pow(d,eps)*Kc*w*w*F.det();
   return w;
 }
+
+
+bool PoroFracture::formPermeabilityTensor (SymmTensor& K,
+                                           const Vectors& eV,
+                                           const FiniteElement& fe,
+                                           const Vec3& X) const
+{
+  if (this->formCrackedPermeabilityTensor(K,eV,fe,X) < 0.0)
+    return false;
+
+  const PoroMaterial* pmat = dynamic_cast<const PoroMaterial*>(material);
+  if (!pmat) return false;
+
+  Vec3 permeability = pmat->getPermeability(X);
+  for (size_t i = 1; i < K.dim(); i++)
+    K(i,i) += permeability(i);
+
+  return true;
+}

PoroFracture.h

@@ -91,6 +91,16 @@ class PoroFracture : public PoroElasticity
                                       const FiniteElement& fe,
                                       const Vec3& X) const;
 
+  //! \brief Evaluates the permeability tensor at a quadrature point.
+  //! \param[out] K The permeability tensor
+  //! \param[in] eV Element solution vectors
+  //! \param[in] fe Finite element data of current integration point
+  //! \param[in] X Cartesian coordinates of current integration point
+  virtual bool formPermeabilityTensor(SymmTensor& K,
+                                      const Vectors& eV,
+                                      const FiniteElement& fe,
+                                      const Vec3& X) const;
+
   //! \brief Computes the permeability tensor of the broken material.
   //! \param[out] Kcrack Permeability tensor of the broken material
   //! \param[in] eV Element solution vectors