PermeabilityOrthotropicPowerLaw.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"
 
namespace ParameterLib
{
struct CoordinateSystem;
}
 
namespace MaterialPropertyLib
{
class Medium;
class Phase;
class Component;

template <int DisplacementDim>
class PermeabilityOrthotropicPowerLaw final : public Property
{
private:
    std::array<double, DisplacementDim> const k_;
    std::array<double, DisplacementDim> const lambda_;
    ParameterLib::CoordinateSystem const* const local_coordinate_system_;
 
public:
    PermeabilityOrthotropicPowerLaw(
        std::string name,
        std::array<double, DisplacementDim>
            intrinsic_permeabilities,
        std::array<double, DisplacementDim>
            exponents,
        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;
};
 
}  // namespace MaterialPropertyLib