Skip to content
Snippets Groups Projects
Commit 6a8957d4 authored by Timo Koch's avatar Timo Koch
Browse files

[spatialparams] Get Forchheimer coeff with getParam and default

parent 983d0f9e
No related branches found
No related tags found
1 merge request!617[WIP] Next
...@@ -361,9 +361,6 @@ private: ...@@ -361,9 +361,6 @@ private:
params["TimeLoop.MaxTimeStepSize"] = "1e300"; params["TimeLoop.MaxTimeStepSize"] = "1e300";
params["TimeLoop.MaxTimeStepDivisions"] = "10"; params["TimeLoop.MaxTimeStepDivisions"] = "10";
// parameters in the spatial params group
params["SpatialParams.ForchCoeff"] = "0.55";
// parameters in the vtk group // parameters in the vtk group
params["Vtk.AddVelocity"] = "false"; params["Vtk.AddVelocity"] = "false";
params["Vtk.AddProcessRank"] = "true"; params["Vtk.AddProcessRank"] = "true";
......
...@@ -49,6 +49,7 @@ class ImplicitSpatialParamsOneP ...@@ -49,6 +49,7 @@ class ImplicitSpatialParamsOneP
using GridView = typename GET_PROP_TYPE(TypeTag, GridView); using GridView = typename GET_PROP_TYPE(TypeTag, GridView);
using Implementation = typename GET_PROP_TYPE(TypeTag, SpatialParams); using Implementation = typename GET_PROP_TYPE(TypeTag, SpatialParams);
using SubControlVolume = typename GET_PROP_TYPE(TypeTag, SubControlVolume); 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 ElementSolutionVector = typename GET_PROP_TYPE(TypeTag, ElementSolutionVector);
using Element = typename GridView::template Codim<0>::Entity; using Element = typename GridView::template Codim<0>::Entity;
...@@ -287,14 +288,12 @@ public: ...@@ -287,14 +288,12 @@ public:
* porous medium. Taking it as a constant is only a first approximation * porous medium. Taking it as a constant is only a first approximation
* (Nield, Bejan, Convection in porous media, 2006, p. 10 \cite nield2006 ) * (Nield, Bejan, Convection in porous media, 2006, p. 10 \cite nield2006 )
* *
* \param element The current finite element * \param scv The sub-control volume face where the
* \param fvGeometry The current finite volume geometry of the element * intrinsic velocity ought to be calculated.
* \param scvIdx The index 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; return forchCoeff;
} }
......
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