van_genuchten_law.h

Go to the documentation of this file.

// KRATOS___
//     //   ) )
//    //         ___      ___
//   //  ____  //___) ) //   ) )
//  //    / / //       //   / /
// ((____/ / ((____   ((___/ /  MECHANICS
//
//  License:         geo_mechanics_application/license.txt
//
//  Main authors:    Vahid Galavi
//
 
#pragma once
 
// System includes
#include <string>
#include <iostream>
#include "includes/define.h"
 
// External includes
 
// Project includes
#include "includes/serializer.h"
#include "custom_retention/retention_law.h"
 
// Application includes
#include "geo_mechanics_application_variables.h"
 
namespace Kratos
{
class KRATOS_API(GEO_MECHANICS_APPLICATION) VanGenuchtenLaw
    : public RetentionLaw
{
public:
    using BaseType = RetentionLaw;
 
    using GeometryType = Geometry<Node>;
 
    using SizeType = std::size_t;
 
    KRATOS_CLASS_POINTER_DEFINITION( VanGenuchtenLaw );
 
    RetentionLaw::Pointer Clone() const override;
 
    void InitializeMaterial(const Properties& rMaterialProperties,
                            const GeometryType& rElementGeometry,
                            const Vector& rShapeFunctionsValues) override;
 
    void Initialize(Parameters &rParameters) override;
 
    void InitializeSolutionStep(Parameters &rParameters) override;
 
    double CalculateSaturation(Parameters &rParameters) override;
 
    double CalculateEffectiveSaturation(Parameters &rParameters) override;
 
    double CalculateDerivativeOfSaturation(Parameters &rParameters) override;
 
    double CalculateRelativePermeability(Parameters &rParameters) override;
 
    double CalculateBishopCoefficient(Parameters &rParameters) override;
 
    void Finalize(Parameters &rParameters) override;
 
    void FinalizeSolutionStep(Parameters &rParameters) override;
 
    double& CalculateValue(RetentionLaw::Parameters& rParameterValues,
                           const Variable<double>& rThisVariable,
                           double& rValue) override;
 
 
    int Check(const Properties& rMaterialProperties,
              const ProcessInfo& rCurrentProcessInfo) override;
 
private:
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, RetentionLaw )
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, RetentionLaw)
    }
 
}; // Class VanGenuchtenLaw
}  // namespace Kratos.