[MatL] Fracture: ConstantPermeability model.

MaterialLib/FractureModels/Permeability/ConstantPermeability.cpp

+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2019, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#include "ConstantPermeability.h"
+
+namespace MaterialLib::Fracture::Permeability
+{
+ConstantPermeability::ConstantPermeability(double const permeability)
+    : _permeability(permeability)
+{
+}
+
+double ConstantPermeability::permeability(
+    PermeabilityState const* const /*state*/,
+    double const /*aperture0*/,
+    double const /*aperture_m*/) const
+{
+    return _permeability;
+}
+
+double ConstantPermeability::dpermeability_daperture(
+    PermeabilityState const* const /*state*/,
+    double const /*aperture0*/,
+    double const /*aperture_m*/) const
+{
+    return 0.;
+}
+}  // namespace MaterialLib::Fracture::Permeability

MaterialLib/FractureModels/Permeability/ConstantPermeability.h

+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2019, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#pragma once
+
+#include "Permeability.h"
+
+namespace MaterialLib::Fracture::Permeability
+{
+/// A constant permeability model.
+class ConstantPermeability final : public Permeability
+{
+public:
+    ConstantPermeability(double const permeability);
+
+private:
+    double permeability(PermeabilityState const* const /*state*/,
+                        double const /*aperture0*/,
+                        double const /*aperture_m*/) const override;
+
+    double dpermeability_daperture(PermeabilityState const* const /*state*/,
+                                   double const /*aperture0*/,
+                                   double const /*aperture_m*/) const override;
+
+private:
+    double const _permeability;
+};
+}  // namespace MaterialLib::Fracture::Permeability

MaterialLib/FractureModels/Permeability/CreateConstantPermeability.cpp

+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2019, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#include "CreateConstantPermeability.h"
+#include "BaseLib/ConfigTree.h"
+#include "BaseLib/Error.h"
+
+#include "ConstantPermeability.h"
+
+namespace MaterialLib::Fracture::Permeability
+{
+std::unique_ptr<Permeability> createConstantPermeability(
+    BaseLib::ConfigTree const& config)
+{
+    //! \ogs_file_param{material__fracture_properties__permeability_model__type}
+    config.checkConfigParameter("type", "ConstantPermeability");
+
+    auto const permeability =
+        //! \ogs_file_param{material__fracture_properties__permeability_model__ConstantPermeability__value}
+        config.getConfigParameter<double>("value");
+
+    if (permeability < 0)
+    {
+        OGS_FATAL(
+            "The permeability parameter must be non-negative. Given value %g.",
+            permeability);
+    }
+
+    return std::make_unique<ConstantPermeability>(permeability);
+}
+}  // namespace MaterialLib::Fracture::Permeability

MaterialLib/FractureModels/Permeability/CreateConstantPermeability.h

+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2019, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#pragma once
+
+#include <memory>
+
+namespace BaseLib
+{
+class ConfigTree;
+}
+namespace MaterialLib::Fracture::Permeability
+{
+class Permeability;
+}
+
+namespace MaterialLib::Fracture::Permeability
+{
+std::unique_ptr<Permeability> createConstantPermeability(
+    BaseLib::ConfigTree const& config);
+}  // namespace MaterialLib::Fracture::Permeability