CreateGasPressureDependentPermeability.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateGasPressureDependentPermeability.h"
 
#include <string>
 
#include "BaseLib/ConfigTree.h"
#include "BaseLib/Error.h"
#include "GasPressureDependentPermeability.h"
#include "Parameter.h"
#include "ParameterLib/CoordinateSystem.h"
#include "ParameterLib/Parameter.h"
#include "ParameterLib/Utils.h"
 
namespace MaterialPropertyLib
{
std::unique_ptr<Property> createGasPressureDependentPermeability(
    int const geometry_dimension,
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    ParameterLib::CoordinateSystem const* const local_coordinate_system)
{
    if ((geometry_dimension != 2) && (geometry_dimension != 3))
    {
        OGS_FATAL(
            "The GasPressureDependentPermeability is implemented only for 2D "
            "or 3D problems");
    }

    config.checkConfigParameter("type", "GasPressureDependentPermeability");
 
    // Second access for storage.
    auto property_name = config.peekConfigParameter<std::string>("name");
 
    DBUG("Create GasPressureDependentPermeability property {:s}.",
         property_name);
 
    std::string const& parameter_name =
        config.getConfigParameter<std::string>("initial_permeability");
    auto const& parameter_k0 = ParameterLib::findParameter<double>(
        parameter_name, parameters, 0, nullptr);
 
    auto const a1 =
        config.getConfigParameter<double>("a1");
    auto const a2 =
        config.getConfigParameter<double>("a2");
    auto const pressure_threshold =
        config.getConfigParameter<double>("pressure_threshold");
    auto const minimum_permeability =
        config.getConfigParameter<double>("minimum_permeability");
    auto const maximum_permeability =
        config.getConfigParameter<double>("maximum_permeability");
 
    if (minimum_permeability > maximum_permeability)
    {
        OGS_FATAL(
            "The value of minimum_permeability of {:e} is larger that the "
            "value of maximum_permeability of {:e} in "
            "GasPressureDependentPermeability",
            minimum_permeability, maximum_permeability);
    }
 
    if (geometry_dimension == 2)
    {
        return std::make_unique<GasPressureDependentPermeability<2>>(
            std::move(property_name), parameter_k0, a1, a2, pressure_threshold,
            minimum_permeability, maximum_permeability,
            local_coordinate_system);
    }
 
    return std::make_unique<GasPressureDependentPermeability<3>>(
        std::move(property_name), parameter_k0, a1, a2, pressure_threshold,
        minimum_permeability, maximum_permeability, local_coordinate_system);
}
}  // namespace MaterialPropertyLib