convection_diffusion_reaction_element.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Suneth Warnakulasuriya
//
 
#if !defined(KRATOS_CONVECTION_DIFFUSION_REACTION_ELEMENT_H_INCLUDED)
#define KRATOS_CONVECTION_DIFFUSION_REACTION_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/element.h"
 
// Application includes
 
namespace Kratos
{
 
template <
    unsigned int TDim,
    unsigned int TNumNodes,
    class TConvectionDiffusionReactionData>
class ConvectionDiffusionReactionElement : public Element
{
public:
 
    using BaseType = Element;
 
    using NodeType = Node;
 
    using GeometryType = Geometry<NodeType>;
 
    using NodesArrayType = Geometry<NodeType>::PointsArrayType;
 
    using VectorType = Vector;
 
    using MatrixType = Matrix;
 
    using IndexType = std::size_t;
 
    using EquationIdVectorType = std::vector<IndexType>;
 
    using DofsVectorType = std::vector<Dof<double>::Pointer>;
 
    using ShapeFunctionDerivativesArrayType = GeometryType::ShapeFunctionsGradientsType;
 
    using PropertiesType = typename BaseType::PropertiesType;
 
    using ConvectionDiffusionReactionDataType = TConvectionDiffusionReactionData;
 
 
    using CurrentElementType =
        ConvectionDiffusionReactionElement<TDim, TNumNodes, TConvectionDiffusionReactionData>;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(ConvectionDiffusionReactionElement);
 
 
    explicit ConvectionDiffusionReactionElement(
        IndexType NewId = 0)
    : Element(NewId)
    {
    }
 
    ConvectionDiffusionReactionElement(
        IndexType NewId,
        const NodesArrayType& ThisNodes)
    : Element(NewId, ThisNodes)
    {
    }
 
    ConvectionDiffusionReactionElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry)
    : Element(NewId, pGeometry)
    {
    }
 
    ConvectionDiffusionReactionElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry,
        typename PropertiesType::Pointer pProperties)
    : Element(NewId, pGeometry, pProperties)
    {
    }
 
    ConvectionDiffusionReactionElement(
        ConvectionDiffusionReactionElement const& rOther)
    : Element(rOther)
    {
    }
 
    ~ConvectionDiffusionReactionElement() override = default;
 
 
    Element::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        typename PropertiesType::Pointer pProperties) const override
    {
        KRATOS_TRY
        return Kratos::make_intrusive<CurrentElementType>(
            NewId, this->GetGeometry().Create(ThisNodes), pProperties);
        KRATOS_CATCH("");
    }
 
    Element::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        typename PropertiesType::Pointer pProperties) const override
    {
        KRATOS_TRY
        return Kratos::make_intrusive<CurrentElementType>(NewId, pGeom, pProperties);
        KRATOS_CATCH("");
    }
 
    Element::Pointer Clone(
        IndexType NewId,
        NodesArrayType const& ThisNodes) const override
    {
        KRATOS_TRY
        return Kratos::make_intrusive<CurrentElementType>(
            NewId, this->GetGeometry().Create(ThisNodes), this->pGetProperties());
        KRATOS_CATCH("");
    }
 
    void EquationIdVector(
        EquationIdVectorType& rResult,
        const ProcessInfo& CurrentProcessInfo) const override;
 
    void GetDofList(
        DofsVectorType& rElementalDofList,
        const ProcessInfo& CurrentProcessInfo) const override;
 
    void GetValuesVector(
        Vector& rValues,
        int Step = 0) const override;
 
    void GetFirstDerivativesVector(
        Vector& rValues,
        int Step = 0) const override;
 
    void GetSecondDerivativesVector(
        Vector& rValues,
        int Step = 0) const override;
 
    void CalculateLocalSystem(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLocalVelocityContribution(
        MatrixType& rDampingMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateMassMatrix(
        MatrixType& rMassMatrix,
        const ProcessInfo& rCurrentProcessInfo) override;
    void CalculateDampingMatrix(
        MatrixType& rDampingMatrix,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
    GeometryData::IntegrationMethod GetIntegrationMethod() const override;
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "ConvectionDiffusionReactionElement #" << Id();
        return buffer.str();
    }
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "CDR" << TConvectionDiffusionReactionData::GetName();
    }
 
 
protected:
 
    void GetValuesArray(
        BoundedVector<double, TNumNodes>& rValues,
        const int Step = 0) const;
 
    virtual void CalculateGeometryData(
        Vector& rGaussWeights,
        Matrix& rNContainer,
        ShapeFunctionDerivativesArrayType& rDN_DX) const;
 
    void AddLumpedMassMatrix(
        Matrix& rMassMatrix,
        const double Mass) const;
 
    void AddDampingMatrixGaussPointContributions(
        Matrix& rDampingMatrix,
        const double ReactionTerm,
        const double EffectiveKinematicViscosity,
        const Vector& rVelocityConvectiveTerms,
        const double GaussWeight,
        const Vector& rGaussShapeFunctions,
        const Matrix& rGaussdNa_dNb) const;
 
 
private:
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_TRY
 
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Element);
 
        KRATOS_CATCH("");
    }
    void load(Serializer& rSerializer) override
    {
        KRATOS_TRY
 
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Element);
 
        KRATOS_CATCH("");
    }
 
 
}; // Class ConvectionDiffusionReactionElement
 
 
template <unsigned int TDim, unsigned int TNumNodes, class TConvectionDiffusionReactionData>
inline std::istream& operator>>(
    std::istream& rIStream,
    ConvectionDiffusionReactionElement<TDim, TNumNodes, TConvectionDiffusionReactionData>& rThis);
 
template <unsigned int TDim, unsigned int TNumNodes, class TConvectionDiffusionReactionData>
inline std::ostream& operator<<(
    std::ostream& rOStream,
    const ConvectionDiffusionReactionElement<TDim, TNumNodes, TConvectionDiffusionReactionData>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << " : " << std::endl;
    rThis.PrintData(rOStream);
    return rOStream;
}
 
 
} // namespace Kratos.
 
#endif // KRATOS_CONVECTION_DIFFUSION_REACTION_ELEMENT_H_INCLUDED defined