time_integrated_qsvms_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_TIME_INTEGRATED_QSVMS_DATA_H)
#define KRATOS_TIME_INTEGRATED_QSVMS_DATA_H
 
#include "fluid_dynamics_application_variables.h"
#include "custom_utilities/qsvms_data.h"
 
namespace Kratos {
 
 
 
template< size_t TDim, size_t TNumNodes >
class TimeIntegratedQSVMSData : public QSVMSData<TDim,TNumNodes,true>
{
public:
 
 
using NodalScalarData = typename FluidElementData<TDim,TNumNodes, false>::NodalScalarData;
using NodalVectorData = typename FluidElementData<TDim,TNumNodes, false>::NodalVectorData;
 
 
NodalVectorData Velocity_OldStep1;
NodalVectorData Velocity_OldStep2;
 
double bdf0;
double bdf1;
double bdf2;
 
 
void Initialize(const Element& rElement, const ProcessInfo& rProcessInfo) override
{
    QSVMSData<TDim,TNumNodes,true>::Initialize(rElement,rProcessInfo);
    const Geometry< Node >& r_geometry = rElement.GetGeometry();
    this->FillFromHistoricalNodalData(Velocity_OldStep1,VELOCITY,r_geometry,1);
    this->FillFromHistoricalNodalData(Velocity_OldStep2,VELOCITY,r_geometry,2);
 
    const Vector& BDFVector = rProcessInfo[BDF_COEFFICIENTS];
    bdf0 = BDFVector[0];
    bdf1 = BDFVector[1];
    bdf2 = BDFVector[2];
}
 
 
static int Check(const Element& rElement, const ProcessInfo& rProcessInfo)
{
    int out = QSVMSData<TDim,TNumNodes,true>::Check(rElement,rProcessInfo);
 
    return out;
}
 
 
};
 
 
 
}
 
#endif