Skip to content
Snippets Groups Projects
Commit 2f193857 authored by Stefanie Kiemle's avatar Stefanie Kiemle
Browse files

spatial params fixed

parent 29134e9e
No related branches found
No related tags found
1 merge request!116spatial params fixed
...@@ -51,7 +51,7 @@ class ExerciseFluidsystemSpatialParams ...@@ -51,7 +51,7 @@ class ExerciseFluidsystemSpatialParams
static constexpr int dim = GridView::dimension; static constexpr int dim = GridView::dimension;
static constexpr int dimWorld = GridView::dimensionworld; static constexpr int dimWorld = GridView::dimensionworld;
using Element = typename GridView::template Codim<0>::Entity; using Element = typename GridView::template Codim<0>::Entity;
using GlobalPosition = typename Element::Geometry::GlobalCoordinate; using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
using PcKrSwCurve = FluidMatrix::BrooksCoreyDefault<Scalar>; using PcKrSwCurve = FluidMatrix::BrooksCoreyDefault<Scalar>;
...@@ -69,7 +69,7 @@ public: ...@@ -69,7 +69,7 @@ public:
, K_(0) , K_(0)
, KLens_(0) , KLens_(0)
, pcKrSwCurve_("SpatialParams") , pcKrSwCurve_("SpatialParams")
, lensPcKrSwCurve_("Lens.SpatialParams") , lensPcKrSwCurve_("SpatialParams.Lens")
{ {
//set main diagonal entries of the permeability tensor to a value //set main diagonal entries of the permeability tensor to a value
//setting to one value means: isotropic, homogeneous //setting to one value means: isotropic, homogeneous
...@@ -111,10 +111,10 @@ public: ...@@ -111,10 +111,10 @@ public:
* \param globalPos The global coordinates for the given location * \param globalPos The global coordinates for the given location
*/ */
auto fluidMatrixInteractionAtPos(const GlobalPosition& globalPos) const auto fluidMatrixInteractionAtPos(const GlobalPosition& globalPos) const
{ {
if (isInLens(globalPos)) if (isInLens(globalPos))
return makeFluidMatrixInteraction(pcKrSwCurve_); return makeFluidMatrixInteraction(pcKrSwCurve_);
return makeFluidMatrixInteraction(lensPcKrSwCurve_); return makeFluidMatrixInteraction(lensPcKrSwCurve_);
} }
/*! /*!
...@@ -150,7 +150,7 @@ private: ...@@ -150,7 +150,7 @@ private:
Dune::FieldMatrix<Scalar, dim, dim> K_; Dune::FieldMatrix<Scalar, dim, dim> K_;
Dune::FieldMatrix<Scalar, dim, dim> KLens_; Dune::FieldMatrix<Scalar, dim, dim> KLens_;
const PcKrSwCurve pcKrSwCurve_; const PcKrSwCurve pcKrSwCurve_;
const PcKrSwCurve lensPcKrSwCurve_; const PcKrSwCurve lensPcKrSwCurve_;
}; };
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment