EmbeddedFracturePermeability.h

Go to the documentation of this file.

 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
#include "MathLib/KelvinVector.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialPropertyLib
{
template <int DisplacementDim>
class EmbeddedFracturePermeability final : public Property
{
private:
    Medium* _medium = nullptr;
    Eigen::Matrix<double, 3, 1> const _n;
    bool const _n_const;
    double const _k;
    double const _b0;
    double const _a;
    double const _e0;
    ParameterLib::Parameter<double> const& _phi_xy;
    ParameterLib::Parameter<double> const& _phi_yz;
    double const _jf;
 
public:
    EmbeddedFracturePermeability(
        std::string name,
        Eigen::Matrix<double, 3, 1> const fracture_normal,
        bool const fracture_normal_is_constant,
        double const intrinsic_permeability,
        double const initial_aperture,
        double const mean_fracture_distance,
        double const threshold_strain,
        ParameterLib::Parameter<double> const& fracture_rotation_xy,
        ParameterLib::Parameter<double> const& fracture_rotation_yz,
        double const jacobian_factor);
 
    static int const KelvinVectorSize =
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
    using SymmetricTensor = Eigen::Matrix<double, KelvinVectorSize, 1>;
 
    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