From 549d78a924bfd5eca589e587404592e100555f5e Mon Sep 17 00:00:00 2001
From: Maziar Veyskarami <maziar.veyskarami@iws.uni-stuttgart.de>
Date: Wed, 9 Feb 2022 11:45:41 +0100
Subject: [PATCH] [examples][porenetwork_upscaling] modify spatial params to
include momentum coefficient needed in non-creeping flow calculation
---
examples/porenetwork_upscaling/spatialparams.hh | 15 ++++++++++++++-
1 file changed, 14 insertions(+), 1 deletion(-)
diff --git a/examples/porenetwork_upscaling/spatialparams.hh b/examples/porenetwork_upscaling/spatialparams.hh
index a0267a274a..d4c6c7c256 100644
--- a/examples/porenetwork_upscaling/spatialparams.hh
+++ b/examples/porenetwork_upscaling/spatialparams.hh
@@ -102,11 +102,24 @@ public:
return M_PI * r * r;
}
+ // dimensionless kinetic-energy coeffiecient which for non-creeping flow
+ template<class ElementSolutionVector>
+ Scalar kineticEnergyCoefficient(const Element& element,
+ const SubControlVolume& scv,
+ const ElementSolutionVector& elemSol) const
+ { return 1.0; }
+
+ // dimensionless momentum coeffiecient which for non-creeping flow
+ template<class ElementSolutionVector>
+ Scalar momentumCoefficient(const Element& element,
+ const SubControlVolume& scv,
+ const ElementSolutionVector& elemSol) const
+ { return 1.0; }
+
private:
std::vector<Scalar> poreShapeFactor_;
};
-
} // namespace Dumux::PoreNetwork
#endif
--
GitLab