Pw_normal_flux_condition.hpp

Go to the documentation of this file.

// KRATOS___
//     //   ) )
//    //         ___      ___
//   //  ____  //___) ) //   ) )
//  //    / / //       //   / /
// ((____/ / ((____   ((___/ /  MECHANICS
//
//  License:         geo_mechanics_application/license.txt
//
//
//  Main authors:    Ignasi de Pouplana,
//                   Vahid Galavi,
//                   Aron Noordam
//
 
 
#if !defined(KRATOS_GEO_PW_NORMAL_FLUX_CONDITION_H_INCLUDED )
#define  KRATOS_GEO_PW_NORMAL_FLUX_CONDITION_H_INCLUDED
 
// Project includes
#include "includes/serializer.h"
 
// Application includes
#include "custom_conditions/Pw_condition.hpp"
#include "custom_utilities/element_utilities.hpp"
#include "custom_utilities/condition_utilities.hpp"
#include "geo_mechanics_application_variables.h"
 
namespace Kratos
{
 
template< unsigned int TDim, unsigned int TNumNodes >
class KRATOS_API(GEO_MECHANICS_APPLICATION) PwNormalFluxCondition : public PwCondition<TDim,TNumNodes>
{
 
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( PwNormalFluxCondition );
    
    using IndexType = std::size_t;
    using PropertiesType = Properties;
    using NodeType = Node;
    using GeometryType = Geometry<NodeType>;
    using NodesArrayType = GeometryType::PointsArrayType;
    using VectorType = Vector;
    using MatrixType = Matrix;
    
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    // Default constructor
    PwNormalFluxCondition() : PwCondition<TDim,TNumNodes>() {}
    
    // Constructor 1
    PwNormalFluxCondition( IndexType NewId, GeometryType::Pointer pGeometry ) : PwCondition<TDim,TNumNodes>(NewId, pGeometry) {}
    
    // Constructor 2
    PwNormalFluxCondition( IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties ) : PwCondition<TDim,TNumNodes>(NewId, pGeometry, pProperties) {}
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    Condition::Pointer Create(IndexType NewId,NodesArrayType const& ThisNodes,PropertiesType::Pointer pProperties ) const override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
protected:
 
    struct NormalFluxVariables
    {
        double NormalFlux;
        double IntegrationCoefficient;
        array_1d<double,TNumNodes> Np;
        array_1d<double,TNumNodes> PVector;
    };
    
    // Member Variables
    
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
                                    
    void CalculateRHS(VectorType& rRightHandSideVector,
                      const ProcessInfo& CurrentProcessInfo) override;
    
    void CalculateAndAddRHS(VectorType& rRightHandSideVector, NormalFluxVariables& rVariables);
 
    virtual double CalculateIntegrationCoefficient(const Matrix& Jacobian, const double& Weight);
    
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
private:
    
    // Member Variables
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    // Serialization
    
    friend class Serializer;
    
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, Condition )
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, Condition )
    }
    
}; // class PwNormalFluxCondition.
 
} // namespace Kratos.
 
#endif // KRATOS_GEO_PW_NORMAL_FLUX_CONDITION_H_INCLUDED defined