OrthotropicEmbeddedFracturePermeability.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 "MathLib/KelvinVector.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialPropertyLib
{
class Medium;
class Phase;
class Component;
 
template <int DisplacementDim>
class OrthotropicEmbeddedFracturePermeability final : public Property
{
private:
    Medium* _medium = nullptr;
    std::vector<double> const _a;
    std::vector<double> const _e0;
    Eigen::Matrix<double, 3, 3> const _n;
    ParameterLib::Parameter<double> const& _k;
    ParameterLib::Parameter<double> const& _phi_xy;
    ParameterLib::Parameter<double> const& _phi_yz;
    double const _jf;
 
public:
    OrthotropicEmbeddedFracturePermeability(
        std::string name,
        std::vector<double> const& mean_fracture_distances,
        std::vector<double> const& threshold_strains,
        Eigen::Matrix<double, 3, 3> const fracture_normals,
        ParameterLib::Parameter<double> const& intrinsic_permeability,
        ParameterLib::Parameter<double> const& fracture_rotation_xy,
        ParameterLib::Parameter<double> const& fracture_rotation_yz,
        double const jacobian_factor);
 
    using SymmetricTensor = Eigen::Matrix<
        double,
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim), 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