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

[stokes_darcy][spatialparams] Runtime porosity parameter

parent e2f73d26
No related branches found
No related tags found
1 merge request!14Navier slip condition as an option for interface condition
...@@ -29,6 +29,7 @@ NumericDifference.BaseEpsilon = 1e-8 ...@@ -29,6 +29,7 @@ NumericDifference.BaseEpsilon = 1e-8
[Darcy] [Darcy]
Problem.Name = darcy Problem.Name = darcy
SpatialParams.Permeability = 1e-4 SpatialParams.Permeability = 1e-4
SpatialParams.Porosity = 1.0
SpatialParams.AlphaBeaversJoseph = 1.0 SpatialParams.AlphaBeaversJoseph = 1.0
[Vtk] [Vtk]
......
...@@ -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 [-].
* *
......
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