createPermeabilityModel.cpp

Go to the documentation of this file.

 
#include "createPermeabilityModel.h"
 
#include <cassert>
 
#include "BaseLib/ConfigTree.h"
#include "BaseLib/Error.h"
#include "MaterialLib/PorousMedium/Permeability/DupuitPermeability.h"
#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
#include "ParameterLib/ConstantParameter.h"
#include "ParameterLib/SpatialPosition.h"
#include "ParameterLib/Utils.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
std::unique_ptr<Permeability> createPermeabilityModel(
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters)
{
    auto const type = config.getConfigParameter<std::string>("type");
 
    if (type == "Constant")
    {
        auto const& permeability_parameter = ParameterLib::findParameter<
            double>(
            config,
            "permeability_tensor_entries", parameters, 0);
 
        int dimension = static_cast<int>(
            std::sqrt(permeability_parameter.getNumberOfGlobalComponents()));
        if (permeability_parameter.getNumberOfGlobalComponents() !=
            dimension * dimension)
        {
            OGS_FATAL(
                "The given parameter has {:d} components, but the permeability "
                "tensor is defined for a {:d} dimensional problem.",
                permeability_parameter.getNumberOfGlobalComponents(),
                dimension);
        }
 
        return std::make_unique<Permeability>(permeability_parameter,
                                              dimension);
    }
 
    if (type == "Dupuit")
    {
        auto const& permeability_parameter = ParameterLib::findParameter<
            double>(
            config,
            "permeability_tensor_entries", parameters, 0);
 
        int dimension = static_cast<int>(
            std::sqrt(permeability_parameter.getNumberOfGlobalComponents()));
        if (permeability_parameter.getNumberOfGlobalComponents() !=
            dimension * dimension)
        {
            OGS_FATAL(
                "The given parameter has {:d} components, but the permeability "
                "tensor is defined for a {:d} dimensional problem.",
                permeability_parameter.getNumberOfGlobalComponents(),
                dimension);
        }
 
        return std::make_unique<DupuitPermeability>(permeability_parameter,
                                                    dimension);
    }
    OGS_FATAL("The permeability type '{:s}' is unavailable.\n",
              "The available types are \n\tConstant.",
              type.data());
}
 
}  // namespace PorousMedium
}  // namespace MaterialLib