GasPressureDependentPermeability.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
#include "MaterialLib/MPL/VariableType.h"
 
namespace ParameterLib
{
struct CoordinateSystem;
template <typename T>
struct Parameter;
}  // namespace ParameterLib
 
namespace MaterialPropertyLib
{
template <int DisplacementDim>
class GasPressureDependentPermeability final : public Property
{
public:
    GasPressureDependentPermeability(
        std::string name, ParameterLib::Parameter<double> const& k0,
        double const a1, double const a2, double const pressure_threshold,
        double const minimum_permeability, double const maximum_permeability,
        ParameterLib::CoordinateSystem const* const local_coordinate_system);
 
    void checkScale() const override;
 
    PropertyDataType value(VariableArray const& variable_array,
                           ParameterLib::SpatialPosition const& pos,
                           double const t, double const dt) const override;
    PropertyDataType dValue(VariableArray const& variable_array,
                            Variable const variable,
                            ParameterLib::SpatialPosition const& pos,
                            double const t, double const dt) const override;
 
private:
    ParameterLib::Parameter<double> const& k0_;
    double const a1_;
    double const a2_;
    double const pressure_threshold_;
    double const minimum_permeability_;
    double const maximum_permeability_;
    ParameterLib::CoordinateSystem const* const local_coordinate_system_;
};
 
extern template class GasPressureDependentPermeability<2>;
extern template class GasPressureDependentPermeability<3>;
 
}  // namespace MaterialPropertyLib