fluid_element_data.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Jordi Cotela
//
 
#if !defined(KRATOS_FLUID_ELEMENT_DATA_H)
#define KRATOS_FLUID_ELEMENT_DATA_H
 
// External includes
 
// Project includes
#include "geometries/geometry.h"
#include "includes/define.h"
#include "includes/node.h"
#include "includes/element.h"
#include "includes/process_info.h"
#include "includes/constitutive_law.h"
 
namespace Kratos
{
 
 
template< size_t TDim, size_t TNumNodes, bool TElementIntegratesInTime >
class KRATOS_API(FLUID_DYNAMICS_APPLICATION) FluidElementData
{
public:
 
    using NodalScalarData = array_1d<double,TNumNodes>;
 
    using NodalVectorData = BoundedMatrix<double,TNumNodes,TDim>;
 
    using NodalTensorData = std::array<BoundedMatrix<double,TDim,TDim>,TNumNodes>;
 
    using ShapeFunctionsType = array_1d<double,TNumNodes>;
 
    using ShapeDerivativesType = BoundedMatrix<double,TNumNodes,TDim>;
 
    using MatrixRowType = MatrixRow< Matrix >;
 
    constexpr static unsigned int Dim = TDim;
 
    constexpr static unsigned int NumNodes = TNumNodes;
 
    constexpr static unsigned int StrainSize = (TDim-1)*3; // 3 in 2D, 6 in 3D
 
    constexpr static bool ElementManagesTimeIntegration = TElementIntegratesInTime;
 
 
    FluidElementData();
 
    virtual ~FluidElementData();
 
    FluidElementData& operator=(FluidElementData const& rOther) = delete;
 
    FluidElementData(FluidElementData const& rOther) = delete;
 
 
    virtual void Initialize(const Element& rElement, const ProcessInfo& rProcessInfo);
 
    static int Check(const Element& rElement, const ProcessInfo& rProcessInfo);
 
    virtual void UpdateGeometryValues(
        unsigned int IntegrationPointIndex,
        double NewWeight,
        const MatrixRowType& rN,
        const ShapeDerivativesType& rDN_DX);
 
 
    unsigned int IntegrationPointIndex;
 
    double Weight;
 
    ShapeFunctionsType N;
 
    ShapeDerivativesType DN_DX;
 
 
    Vector StrainRate;
 
 
    Vector ShearStress;
 
 
    Matrix C;
 
    ConstitutiveLaw::Parameters ConstitutiveLawValues;
 
    double EffectiveViscosity;
 
protected:
 
 
    //TODO: This needs to be removed
    void FillFromNodalData(
        NodalScalarData &rData,
        const Variable<double> &rVariable,
        const Geometry<Node> &rGeometry)
    {
        KRATOS_WARNING("FluidElementData") << "\'FillFromNodalData\' is deprecated. Use \'FillFromHistoricalNodalData\' instead." << std::endl;
        FillFromHistoricalNodalData(rData, rVariable, rGeometry);
    }
 
    //TODO: This needs to be removed
    void FillFromNodalData(
        NodalVectorData &rData,
        const Variable<array_1d<double, 3>> &rVariable,
        const Geometry<Node> &rGeometry)
    {
        KRATOS_WARNING("FluidElementData") << "\'FillFromNodalData\' is deprecated. Use \'FillFromHistoricalNodalData\' instead." << std::endl;
        FillFromHistoricalNodalData(rData, rVariable, rGeometry);
    }
 
    void FillFromHistoricalNodalData(
        NodalScalarData &rData,
        const Variable<double> &rVariable,
        const Geometry<Node> &rGeometry);
 
    void FillFromHistoricalNodalData(
        NodalVectorData &rData,
        const Variable<array_1d<double, 3>> &rVariable,
        const Geometry<Node> &rGeometry);
 
    void FillFromHistoricalNodalData(
        NodalTensorData& rData,
        const Variable<Matrix>& rVariable,
        const Geometry<Node>& rGeometry);
 
    void FillFromHistoricalNodalData(NodalScalarData& rData, const Variable<double>& rVariable, const Geometry<Node>& rGeometry, const unsigned int Step);
 
    void FillFromHistoricalNodalData(NodalVectorData& rData, const Variable<array_1d<double,3>>& rVariable, const Geometry<Node>& rGeometry, const unsigned int Step);
 
    void FillFromNonHistoricalNodalData(
        NodalScalarData& rData,
        const Variable<double>& rVariable,
        const Geometry<Node>& rGeometry);
 
    void FillFromNonHistoricalNodalData(
        NodalVectorData& rData,
        const Variable<array_1d<double,3>>& rVariable,
        const Geometry<Node>& rGeometry);
 
    void FillFromProcessInfo(double& rData, const Variable<double>& rVariable, const ProcessInfo& rProcessInfo);
 
    void FillFromProcessInfo(int& rData, const Variable<int>& rVariable, const ProcessInfo& rProcessInfo);
 
    void FillFromElementData(double& rData, const Variable<double>& rVariable, const Element& rElement);
 
    void FillFromElementData(Vector& rData, const Variable<Vector>& rVariable, const Element& rElement);
 
    void FillFromElementData(NodalScalarData& rData, const Variable<Vector>& rVariable, const Element& rElement);
 
    void FillFromProperties(double& rData, const Variable<double>& rVariable, const Properties& rProperties);
 
};
 
 
}
 
#endif // KRATOS_FLUID_ELEMENT_DATA_H