From 7bd68a288dc4bdd0da03e98ac0764e4d9fb4f03a Mon Sep 17 00:00:00 2001
From: Ned Coltman <edward.coltman@iws.uni-stuttgart.de>
Date: Mon, 4 Oct 2021 09:54:37 +0200
Subject: [PATCH] [stokes_darcy][spatialparams] Runtime porosity parameter
---
appl/coupled/stokes_darcy/beaversjoseph/params.input | 1 +
appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh | 3 ++-
2 files changed, 3 insertions(+), 1 deletion(-)
diff --git a/appl/coupled/stokes_darcy/beaversjoseph/params.input b/appl/coupled/stokes_darcy/beaversjoseph/params.input
index 0219f08..1f0bfca 100644
--- a/appl/coupled/stokes_darcy/beaversjoseph/params.input
+++ b/appl/coupled/stokes_darcy/beaversjoseph/params.input
@@ -29,6 +29,7 @@ NumericDifference.BaseEpsilon = 1e-8
[Darcy]
Problem.Name = darcy
SpatialParams.Permeability = 1e-4
+SpatialParams.Porosity = 1.0
SpatialParams.AlphaBeaversJoseph = 1.0
[Vtk]
diff --git a/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh b/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh
index 18113c4..18bd803 100644
--- a/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh
+++ b/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh
@@ -56,6 +56,7 @@ public:
{
permeability_ = getParam<Scalar>("Darcy.SpatialParams.Permeability");
alphaBJ_ = getParam<Scalar>("Darcy.SpatialParams.AlphaBeaversJoseph");
+ porosity_ = getParam<Scalar>("Darcy.SpatialParams.Porosity");
}
/*!
@@ -72,7 +73,7 @@ public:
* \param globalPos The global position
*/
Scalar porosityAtPos(const GlobalPosition& globalPos) const
- { return 0.4; }
+ { return porosity_; }
/*! \brief Defines the Beavers-Joseph coefficient in [-].
*
--
GitLab