axisymmetric_eulerian_convection_diffusion.h

Go to the documentation of this file.

// KRATOS ___ ___  _  ___   __   ___ ___ ___ ___
//       / __/ _ \| \| \ \ / /__|   \_ _| __| __|
//      | (_| (_) | .` |\ V /___| |) | || _|| _|
//       \___\___/|_|\_| \_/    |___/___|_| |_|  APPLICATION
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Zorrilla
//
 
#pragma once
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "geometries/geometry_data.h"
#include "includes/define.h"
#include "includes/serializer.h"
 
// Application includes
#include "eulerian_conv_diff.h"
 
 
 
 
 
 
 
namespace Kratos
{
 
template< unsigned int TDim, unsigned int TNumNodes>
class KRATOS_API(CONVECTION_DIFFUSION_APPLICATION) AxisymmetricEulerianConvectionDiffusionElement
    : public EulerianConvectionDiffusionElement<TDim, TNumNodes>
{
public:
 
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(AxisymmetricEulerianConvectionDiffusionElement);
 
    using BaseType = EulerianConvectionDiffusionElement<TDim, TNumNodes>;
 
    using GeometryType = typename BaseType::GeometryType;
 
    using PropertiesType = typename BaseType::PropertiesType;
 
    using IndexType = typename BaseType::IndexType;
 
    using SizeType = typename BaseType::SizeType;
 
    using VectorType = typename BaseType::VectorType;
 
    using MatrixType = typename BaseType::MatrixType;
 
    using NodesArrayType = typename BaseType::NodesArrayType;
 
    using ShapeFunctionsGradientsType = typename GeometryType::ShapeFunctionsGradientsType;
 
 
    AxisymmetricEulerianConvectionDiffusionElement() : BaseType()
    {
    }
 
    AxisymmetricEulerianConvectionDiffusionElement(
        IndexType NewId,
        typename GeometryType::Pointer pGeometry)
        : BaseType(NewId, pGeometry)
    {}
 
    AxisymmetricEulerianConvectionDiffusionElement(
        IndexType NewId,
        typename GeometryType::Pointer pGeometry,
        typename PropertiesType::Pointer pProperties)
        : BaseType(NewId, pGeometry, pProperties)
    {}
 
    virtual ~AxisymmetricEulerianConvectionDiffusionElement(){};
 
 
 
 
    Element::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        typename PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<AxisymmetricEulerianConvectionDiffusionElement<TDim, TNumNodes>>(NewId, this->GetGeometry().Create(ThisNodes), pProperties);
    }
 
    Element::Pointer Create(
        IndexType NewId,
        typename GeometryType::Pointer pGeom,
        typename PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<AxisymmetricEulerianConvectionDiffusionElement<TDim, TNumNodes>>(NewId, pGeom, pProperties);
    }
 
    void CalculateLocalSystem(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
 
    int Check(const ProcessInfo &rCurrentProcessInfo) const override;
 
 
    std::string Info() const override
    {
        return "AxisymmetricConvectionDiffusion #";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << this->Info() << this->Id();
    }
 
 
 
 
protected:
 
 
 
 
 
 
 
    double CalculateTau(
        const typename BaseType::ElementVariables &rVariables,
        const double NormVelocity,
        const double ElementSize);
 
 
 
 
 
 
 
private:
 
    static constexpr GeometryData::IntegrationMethod mIntegrationMethod = GeometryData::IntegrationMethod::GI_GAUSS_2;
 
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, BaseType);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, BaseType);
    }
 
 
 
 
    void CalculateGaussPointData(
        const array_1d<double, TNumNodes> &rN,
        const BoundedMatrix<double, TNumNodes, TDim> &rDNDX,
        typename BaseType::ElementVariables &rVariables,
        double &rRadius,
        array_1d<double, TDim> &rVelocity,
        array_1d<double, TNumNodes> &rConvectionOperator,
        BoundedMatrix<double, TDim, TDim> &rVelocityGradient) const;
 
 
 
 
 
 
 
}; // Class AxisymmetricConvectionDiffusion
 
 
 
} // namespace Kratos.