diff --git a/dumux/common/parameters.hh b/dumux/common/parameters.hh index 88e7cf5cdc975dad1472b08f0fed896da6bd3f3b..5f952529e7fa972a941d8dd6f966c02e9bd4e954 100644 --- a/dumux/common/parameters.hh +++ b/dumux/common/parameters.hh @@ -361,9 +361,6 @@ private: params["TimeLoop.MaxTimeStepSize"] = "1e300"; params["TimeLoop.MaxTimeStepDivisions"] = "10"; - // parameters in the spatial params group - params["SpatialParams.ForchCoeff"] = "0.55"; - // parameters in the vtk group params["Vtk.AddVelocity"] = "false"; params["Vtk.AddProcessRank"] = "true"; diff --git a/dumux/material/spatialparams/implicit1p.hh b/dumux/material/spatialparams/implicit1p.hh index 5d41b887b61bb0eb3925e24a0c772eec1984305c..bcc06400ac346c5b40ae6a094dc37506c866005e 100644 --- a/dumux/material/spatialparams/implicit1p.hh +++ b/dumux/material/spatialparams/implicit1p.hh @@ -49,6 +49,7 @@ class ImplicitSpatialParamsOneP using GridView = typename GET_PROP_TYPE(TypeTag, GridView); using Implementation = typename GET_PROP_TYPE(TypeTag, SpatialParams); using SubControlVolume = typename GET_PROP_TYPE(TypeTag, SubControlVolume); + using SubControlVolumeFace = typename GET_PROP_TYPE(TypeTag, SubControlVolumeFace); using ElementSolutionVector = typename GET_PROP_TYPE(TypeTag, ElementSolutionVector); using Element = typename GridView::template Codim<0>::Entity; @@ -287,14 +288,12 @@ public: * porous medium. Taking it as a constant is only a first approximation * (Nield, Bejan, Convection in porous media, 2006, p. 10 \cite nield2006 ) * - * \param element The current finite element - * \param fvGeometry The current finite volume geometry of the element - * \param scvIdx The index sub-control volume face where the - * intrinsic velocity ought to be calculated. + * \param scv The sub-control volume face where the + * intrinsic velocity ought to be calculated. */ - Scalar forchCoeff(const SubControlVolume &scv) const + Scalar forchCoeff(const SubControlVolumeFace &scvf) const { - static Scalar forchCoeff = getParamFromGroup<Scalar>(GET_PROP_VALUE(TypeTag, ModelParameterGroup), "SpatialParams.ForchCoeff"); + static Scalar forchCoeff = getParamFromGroup<Scalar>(GET_PROP_VALUE(TypeTag, ModelParameterGroup), "SpatialParams.ForchCoeff", 0.55); return forchCoeff; }