epsilon_element_data.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_K_EPSILON_HIGH_RE_ELEMENT_DATA_EPSILON_ELEMENT_DATA_H_INCLUDED)
#define KRATOS_K_EPSILON_HIGH_RE_ELEMENT_DATA_EPSILON_ELEMENT_DATA_H_INCLUDED
 
// System includes
 
// Project includes
#include "containers/variable.h"
#include "geometries/geometry_data.h"
#include "includes/node.h"
#include "includes/process_info.h"
#include "includes/ublas_interface.h"
 
// Application includes
#include "custom_elements/convection_diffusion_reaction_element_data.h"
 
namespace Kratos
{
 
namespace KEpsilonElementData
{
template <unsigned int TDim>
class EpsilonElementData : public ConvectionDiffusionReactionElementData<TDim>
{
public:
 
    using BaseType = ConvectionDiffusionReactionElementData<TDim>;
 
    using NodeType = Node;
 
    using GeometryType = typename BaseType::GeometryType;
 
 
    EpsilonElementData(
        const GeometryType& rGeometry,
        const Properties& rProperties,
        const ProcessInfo& rProcessInfo)
        : BaseType(rGeometry, rProperties, rProcessInfo)
    {
    }
 
    ~EpsilonElementData() override = default;
 
 
    static const Variable<double>& GetScalarVariable();
 
    static void Check(
        const Element& rElement,
        const ProcessInfo& rCurrentProcessInfo);
 
    static const std::string GetName()
    {
        return "KEpsilonEpsilonElementData";
    }
 
 
    void CalculateConstants(
        const ProcessInfo& rCurrentProcessInfo);
 
    void CalculateGaussPointData(
        const Vector& rShapeFunctions,
        const Matrix& rShapeFunctionDerivatives,
        const int Step = 0);
 
 
protected:
 
    using BaseType::mEffectiveVelocity;
    using BaseType::mEffectiveKinematicViscosity;
    using BaseType::mReactionTerm;
    using BaseType::mSourceTerm;
 
    BoundedMatrix<double, TDim, TDim> mVelocityGradient;
 
    double mC1;
    double mC2;
    double mCmu;
    double mGamma;
    double mTurbulentKineticEnergy;
    double mTurbulentKinematicViscosity;
    double mKinematicViscosity;
    double mVelocityDivergence;
    double mInvEpsilonSigma;
    double mDensity;
 
};
 
 
} // namespace KEpsilonElementData
 
} // namespace Kratos
 
#endif