boussinesq_condition.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Miguel Maso Sotomayor
//
 
#ifndef KRATOS_BOUSSINESQ_CONDITION_H_INCLUDED
#define KRATOS_BOUSSINESQ_CONDITION_H_INCLUDED
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "primitive_condition.h"
 
namespace Kratos
{
 
 
 
 
template<std::size_t TNumNodes>
class BoussinesqCondition : public PrimitiveCondition<TNumNodes>
{
public:
 
    typedef std::size_t IndexType;
 
    typedef Node NodeType;
 
    typedef Geometry<NodeType> GeometryType;
 
    typedef PrimitiveCondition<TNumNodes> BaseType;
 
    typedef typename BaseType::NodesArrayType NodesArrayType;
 
    typedef typename BaseType::PropertiesType PropertiesType;
 
    typedef typename BaseType::EquationIdVectorType EquationIdVectorType;
 
    typedef typename BaseType::DofsVectorType DofsVectorType;
 
    typedef typename BaseType::ConditionData ConditionData;
 
    typedef typename BaseType::LocalVectorType LocalVectorType;
 
    typedef typename BaseType::LocalMatrixType LocalMatrixType;
 
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(BoussinesqCondition);
 
 
    BoussinesqCondition() : BaseType(){}
 
    BoussinesqCondition(IndexType NewId, const NodesArrayType& ThisNodes) : BaseType(NewId, ThisNodes){}
 
    BoussinesqCondition(IndexType NewId, GeometryType::Pointer pGeometry) : BaseType(NewId, pGeometry){}
 
    BoussinesqCondition(IndexType NewId, GeometryType::Pointer pGeometry, typename PropertiesType::Pointer pProperties) : BaseType(NewId, pGeometry, pProperties){}
 
    ~ BoussinesqCondition() override {};
 
 
    Condition::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, typename PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<BoussinesqCondition<TNumNodes>>(NewId, this->GetGeometry().Create(ThisNodes), pProperties);
    }
 
    Condition::Pointer Create(IndexType NewId, GeometryType::Pointer pGeom, typename PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<BoussinesqCondition<TNumNodes>>(NewId, pGeom, pProperties);
    }
 
    Condition::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override
    {
        Condition::Pointer p_new_cond = Create(NewId, this->GetGeometry().Create(ThisNodes), this->pGetProperties());
        p_new_cond->SetData(this->GetData());
        p_new_cond->Set(Flags(*this));
        return p_new_cond;
    }
 
    void InitializeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    const Parameters GetSpecifications() const override;
 
 
    std::string Info() const override
    {
        return "BoussinesqCondition";
    }
 
 
 
 
protected:
 
    static constexpr IndexType mLocalSize = BaseType::mLocalSize;
 
 
    void AddDispersionProjection(
        LocalVectorType& rDispersionH,
        LocalVectorType& rDispersionU,
        const GeometryType& rParentGeometry,
        const ConditionData& rData,
        const array_1d<double,TNumNodes>& rN,
        const Matrix& rDN_DX,
        const double Weight = 1.0);
 
    void CalculateShapeFunctionDerivatives(
        Matrix& rDN_DX,
        const GeometryType& rParentGeometry,
        const Point& rPoint);
 
 
private:
 
 
 
 
 
    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 BoussinesqCondition
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_BOUSSINESQ_CONDITION_H_INCLUDED  defined