conservative_element.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_CONSERVATIVE_ELEMENT_H_INCLUDED
#define KRATOS_CONSERVATIVE_ELEMENT_H_INCLUDED
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "wave_element.h"
 
namespace Kratos
{
 
 
 
 
 
 
template<std::size_t TNumNodes>
class ConservativeElement : public WaveElement<TNumNodes>
{
public:
 
    typedef std::size_t IndexType;
 
    typedef Node NodeType;
 
    typedef Geometry<NodeType> GeometryType;
 
    typedef WaveElement<TNumNodes> WaveElementType;
 
    typedef typename WaveElementType::NodesArrayType NodesArrayType;
 
    typedef typename WaveElementType::PropertiesType PropertiesType;
 
    typedef typename WaveElementType::ElementData ElementData;
 
    typedef typename WaveElementType::LocalMatrixType LocalMatrixType;
 
    typedef typename WaveElementType::LocalVectorType LocalVectorType;
 
    typedef typename WaveElementType::ShapeFunctionsGradientsType ShapeFunctionsGradientsType;
 
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(ConservativeElement);
 
 
    ConservativeElement() : WaveElementType(){}
 
    ConservativeElement(IndexType NewId, const NodesArrayType& ThisNodes) : WaveElementType(NewId, ThisNodes){}
 
    ConservativeElement(IndexType NewId, GeometryType::Pointer pGeometry) : WaveElementType(NewId, pGeometry){}
 
    ConservativeElement(IndexType NewId, GeometryType::Pointer pGeometry, typename PropertiesType::Pointer pProperties) : WaveElementType(NewId, pGeometry, pProperties){}
 
    ~ ConservativeElement() override {};
 
 
    Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, typename PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<ConservativeElement<TNumNodes>>(NewId, this->GetGeometry().Create(ThisNodes), pProperties);
    }
 
    Element::Pointer Create(IndexType NewId, GeometryType::Pointer pGeom, typename PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<ConservativeElement<TNumNodes>>(NewId, pGeom, pProperties);
    }
 
    Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override
    {
        Element::Pointer p_new_elem = Create(NewId, this->GetGeometry().Create(ThisNodes), this->pGetProperties());
        p_new_elem->SetData(this->GetData());
        p_new_elem->Set(Flags(*this));
        return p_new_elem;
    }
 
 
    const Parameters GetSpecifications() const override;
 
 
    std::string Info() const override
    {
        return "ConservativeElement";
    }
 
 
protected:
 
    static constexpr IndexType mLocalSize = WaveElementType::mLocalSize;
 
 
    const Variable<double>& GetUnknownComponent(int Index) const override;
 
    LocalVectorType GetUnknownVector(const ElementData& rData) const override;
 
    void GetNodalData(ElementData& rData, const GeometryType& rGeometry, int Step = 0) override;
 
    void UpdateGaussPointData(ElementData& rData, const array_1d<double,TNumNodes>& rN) override;
 
    void CalculateArtificialViscosity(
        BoundedMatrix<double,3,3>& rViscosity,
        BoundedMatrix<double,2,2>& rDiffusion,
        const ElementData& rData,
        const array_1d<double,TNumNodes>& rN,
        const BoundedMatrix<double,TNumNodes,2>& rDN_DX) override;
 
    void CalculateArtificialDamping(
        BoundedMatrix<double,3,3>& rFriction,
        const ElementData& rData) override;
 
    double StabilizationParameter(const ElementData& rData) const override;
 
    void CalculateGradient(array_1d<double,2>& rGradient, const GeometryType& rGeometry);
 
    void CalculateEdgeUnitNormal(array_1d<double,2>& rNormal, const GeometryType& rGeometry);
 
 
private:
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Element);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Element);
    }
 
 
 
 
 
 
}; // Class ConservativeElement
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_CONSERVATIVE_ELEMENT_H_INCLUDED  defined