Skip to content
Snippets Groups Projects
Commit c3b2882a authored by Ned Coltman's avatar Ned Coltman
Browse files

[spatialparams] Add porosity as input parameter

parent 2c2489f6
No related branches found
No related tags found
1 merge request!18Use Box for fancy interface conditions
...@@ -27,6 +27,7 @@ NumericDifference.BaseEpsilon = 1e-8 ...@@ -27,6 +27,7 @@ NumericDifference.BaseEpsilon = 1e-8
Problem.Name = darcy Problem.Name = darcy
SpatialParams.Permeability = 1e-4 SpatialParams.Permeability = 1e-4
SpatialParams.AlphaBeaversJoseph = 1.0 SpatialParams.AlphaBeaversJoseph = 1.0
SpatialParams.Porosity = 0.4
[Vtk] [Vtk]
WriteFaceData = false WriteFaceData = false
......
...@@ -56,6 +56,7 @@ public: ...@@ -56,6 +56,7 @@ public:
{ {
permeability_ = getParam<Scalar>("Darcy.SpatialParams.Permeability"); permeability_ = getParam<Scalar>("Darcy.SpatialParams.Permeability");
alphaBJ_ = getParam<Scalar>("Darcy.SpatialParams.AlphaBeaversJoseph"); alphaBJ_ = getParam<Scalar>("Darcy.SpatialParams.AlphaBeaversJoseph");
porosity_ = getParam<Scalar>("Darcy.SpatialParams.Porosity");
} }
/*! /*!
...@@ -72,7 +73,7 @@ public: ...@@ -72,7 +73,7 @@ public:
* \param globalPos The global position * \param globalPos The global position
*/ */
Scalar porosityAtPos(const GlobalPosition& globalPos) const Scalar porosityAtPos(const GlobalPosition& globalPos) const
{ return 0.4; } { return porosity_; }
/*! \brief Defines the Beavers-Joseph coefficient in [-]. /*! \brief Defines the Beavers-Joseph coefficient in [-].
* *
...@@ -85,6 +86,7 @@ public: ...@@ -85,6 +86,7 @@ public:
private: private:
Scalar permeability_; Scalar permeability_;
Scalar alphaBJ_; Scalar alphaBJ_;
Scalar porosity_;
}; };
} // end namespace Dumux } // end namespace Dumux
......
...@@ -27,6 +27,7 @@ NumericDifference.BaseEpsilon = 1e-8 ...@@ -27,6 +27,7 @@ NumericDifference.BaseEpsilon = 1e-8
Problem.Name = darcy_filter Problem.Name = darcy_filter
SpatialParams.Permeability = 1e-2 SpatialParams.Permeability = 1e-2
SpatialParams.AlphaBeaversJoseph = 1.0 SpatialParams.AlphaBeaversJoseph = 1.0
SpatialParams.Porosity = 0.4
[Vtk] [Vtk]
WriteFaceData = false WriteFaceData = false
......
...@@ -56,6 +56,7 @@ public: ...@@ -56,6 +56,7 @@ public:
{ {
permeability_ = getParam<Scalar>("Darcy.SpatialParams.Permeability"); permeability_ = getParam<Scalar>("Darcy.SpatialParams.Permeability");
alphaBJ_ = getParam<Scalar>("Darcy.SpatialParams.AlphaBeaversJoseph"); alphaBJ_ = getParam<Scalar>("Darcy.SpatialParams.AlphaBeaversJoseph");
porosity_ = getParam<Scalar>("Darcy.SpatialParams.Porosity");
} }
/*! /*!
...@@ -72,7 +73,7 @@ public: ...@@ -72,7 +73,7 @@ public:
* \param globalPos The global position * \param globalPos The global position
*/ */
Scalar porosityAtPos(const GlobalPosition& globalPos) const Scalar porosityAtPos(const GlobalPosition& globalPos) const
{ return 0.4; } { return porosity_; }
/*! \brief Defines the Beavers-Joseph coefficient in [-]. /*! \brief Defines the Beavers-Joseph coefficient in [-].
* *
...@@ -85,6 +86,7 @@ public: ...@@ -85,6 +86,7 @@ public:
private: private:
Scalar permeability_; Scalar permeability_;
Scalar alphaBJ_; Scalar alphaBJ_;
Scalar porosity_;
}; };
} // end namespace Dumux } // end namespace Dumux
......
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