CreateRelPermGeneralizedPowerNonwettingPhase.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "BaseLib/ConfigTree.h"
#include "RelPermGeneralizedPowerNonwettingPhase.h"
 
namespace MaterialPropertyLib
{
std::unique_ptr<RelPermGeneralizedPowerNonwettingPhase>
createRelPermGeneralizedPowerNonwettingPhase(BaseLib::ConfigTree const& config)
{
    config.checkConfigParameter(
        "type", "RelativePermeabilityGeneralizedPowerNonwettingPhase");
 
    // Second access for storage.
    auto property_name = config.peekConfigParameter<std::string>("name");
 
    DBUG("Create RelPermGeneralizedPowerNonwettingPhase medium property {:s}.",
         property_name);
 
    auto const residual_liquid_saturation =
        config.getConfigParameter<double>("residual_liquid_saturation");
    auto const residual_gas_saturation =
        config.getConfigParameter<double>("residual_gas_saturation");
    auto const min_relative_permeability =
        config.getConfigParameter<double>("min_relative_permeability");
    double const a =
        config.getConfigParameter<double>("enhancement_factor", 1000);
    double const lambda =
        config.getConfigParameter<double>("exponent", 0.8);
 
    if (min_relative_permeability < 0)
    {
        OGS_FATAL("Minimal relative permeability must be non-negative.");
    }
 
    return std::make_unique<RelPermGeneralizedPowerNonwettingPhase>(
        std::move(property_name), residual_liquid_saturation,
        residual_gas_saturation, min_relative_permeability, a, lambda);
}
}  // namespace MaterialPropertyLib