From 52800b72c5e47ca05f1986aa80acd6e99b42053a Mon Sep 17 00:00:00 2001
From: Timo Koch <timokoch@uio.no>
Date: Sun, 6 Oct 2024 16:01:38 +0200
Subject: [PATCH] [hyperelastic] Add dynamic default spatial parameters adding
 density

---
 .../geomechanics/hyperelastic/spatialparams.hh | 18 ++++++++++++++++++
 1 file changed, 18 insertions(+)

diff --git a/dumux/geomechanics/hyperelastic/spatialparams.hh b/dumux/geomechanics/hyperelastic/spatialparams.hh
index 8b5cf2cac1..6f0859eee6 100644
--- a/dumux/geomechanics/hyperelastic/spatialparams.hh
+++ b/dumux/geomechanics/hyperelastic/spatialparams.hh
@@ -51,6 +51,24 @@ private:
     Scalar E_, nu_, mu_, K_, lambda_;
 };
 
+template<class GridGeometry, class Scalar>
+class DefaultDynamicHyperelasticSpatialParams
+: public DefaultHyperelasticSpatialParams<GridGeometry, Scalar>
+{
+    using ParentType = DefaultHyperelasticSpatialParams<GridGeometry, Scalar>;
+public:
+    DefaultDynamicHyperelasticSpatialParams(std::shared_ptr<const GridGeometry> gridGeometry)
+    : ParentType(gridGeometry)
+    , rho_(getParam<Scalar>("SpatialParams.SolidDensity"))
+    {}
+
+    Scalar solidDensity() const
+    { return rho_; }
+
+private:
+    Scalar rho_;
+};
+
 } // end namespace Dumux
 
 #endif
-- 
GitLab