dynamic_vms.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_DYNAMIC_VMS_H_INCLUDED )
#define  KRATOS_DYNAMIC_VMS_H_INCLUDED
 
// System includes
#include <string>
#include <iostream>
 
 
// External includes
 
 
// Project includes
#include "containers/array_1d.h"
#include "includes/define.h"
#include "includes/element.h"
#include "includes/serializer.h"
#include "geometries/geometry.h"
#include "utilities/math_utils.h"
 
// Application includes
#include "fluid_dynamics_application_variables.h"
 
namespace Kratos
{
 
 
 
 
 
 
 
 
template< unsigned int TDim >
class DynamicVMS : public Element
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(DynamicVMS);
 
    typedef Kratos::Vector ShapeFunctionsType;
 
    typedef Kratos::Matrix ShapeDerivativesType;
 
 
    //Constructors.
 
 
    DynamicVMS(IndexType NewId, const NodesArrayType& ThisNodes);
 
 
    DynamicVMS(IndexType NewId, GeometryType::Pointer pGeometry);
 
 
    DynamicVMS(IndexType NewId, GeometryType::Pointer pGeometry, const GeometryData::IntegrationMethod ThisIntegrationMethod);
 
 
    DynamicVMS(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    DynamicVMS(IndexType NewId, GeometryType::Pointer pGeometry, Properties::Pointer pProperties, const GeometryData::IntegrationMethod ThisIntegrationMethod);
 
    ~DynamicVMS() override;
 
 
 
 
 
    Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes,
                                    PropertiesType::Pointer pProperties) const override;
 
 
    Element::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        PropertiesType::Pointer pProperties) const override;
 
    void Initialize(const ProcessInfo &rCurrentProcessInfo) override;
 
    void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void InitializeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) 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 EquationIdVector(EquationIdVectorType& rResult,
                          const ProcessInfo& rCurrentProcessInfo) const override;
 
 
    void GetDofList(DofsVectorType& rElementalDofList,
                    const ProcessInfo& rCurrentProcessInfo) const override;
 
 
    void GetFirstDerivativesVector(Vector& rValues, int Step = 0) const override;
 
 
    void GetSecondDerivativesVector(Vector& rValues, int Step = 0) const override;
 
 
//    virtual void Calculate(const Variable<double>& rVariable,
//                           double& rOutput,
//                           const ProcessInfo& rCurrentProcessInfo);
 
    void Calculate(const Variable< array_1d<double,3> >& rVariable,
                           array_1d<double,3>& rOutput,
                           const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculateOnIntegrationPoints(const Variable<array_1d<double, 3 > >& rVariable,
            std::vector<array_1d<double, 3 > >& rOutput,
            const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(const Variable<double>& rVariable,
            std::vector<double>& rValues,
            const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(const Variable<array_1d<double, 6 > >& rVariable,
            std::vector<array_1d<double, 6 > >& rValues,
            const ProcessInfo& rCurrentProcessInfo) override
    {}
 
    void CalculateOnIntegrationPoints(const Variable<Vector>& rVariable,
            std::vector<Vector>& rValues,
            const ProcessInfo& rCurrentProcessInfo) override
    {}
 
    void CalculateOnIntegrationPoints(const Variable<Matrix>& rVariable,
            std::vector<Matrix>& rValues,
            const ProcessInfo& rCurrentProcessInfo) override
    {}
 
    void SetValuesOnIntegrationPoints(const Variable<double> &rVariable, const std::vector<double> &rValues, const ProcessInfo &rCurrentProcessInfo) override;
 
 
 
    GeometryData::IntegrationMethod GetIntegrationMethod() const override;
 
 
 
 
    std::string Info() const override;
 
    void PrintInfo(std::ostream& rOStream) const override;
 
    void PrintData(std::ostream &rOStream) const override;
 
 
 
 
protected:
 
 
 
 
 
 
 
    virtual void CalculateGeometryData();
 
 
    virtual void UpdateSubscale(const ProcessInfo& rCurrentProcessInfo);
 
    virtual void LinearUpdateSubscale(const ProcessInfo& rCurrentProcessInfo);
 
    virtual void CalculateASGSVelocityContribution(MatrixType& rDampingMatrix,
                                                   VectorType& rRightHandSideVector,
                                                   const ProcessInfo& rCurrentProcessInfo);
 
    virtual void CalculateOSSVelocityContribution(MatrixType& rDampingMatrix,
                                                  VectorType& rRightHandSideVector,
                                                  const ProcessInfo& rCurrentProcessInfo);
 
    virtual void LumpedMassMatrix(MatrixType& rMassMatrix);
 
    virtual void ConsistentMassMatrix(MatrixType& rMassMatrix);
 
 
    template< typename TValueType >
    void EvaluateInPoint(TValueType& rValue,
                         const Kratos::Variable< TValueType >& rVariable,
                         const ShapeFunctionsType& rN)
    {
        GeometryType& rGeom = this->GetGeometry();
        const unsigned int NumNodes = rGeom.PointsNumber();
 
        rValue = rN[0] * rGeom[0].FastGetSolutionStepValue(rVariable);
        for (unsigned int i = 1; i < NumNodes; i++)
            rValue += rN[i] * rGeom[i].FastGetSolutionStepValue(rVariable);
    }
 
    template< typename TValueType >
    void EvaluateInPoint(TValueType& rValue,
                         const Kratos::Variable< TValueType >& rVariable,
                         const ShapeFunctionsType& rN,
                         const unsigned int Step)
    {
        const GeometryType& rGeom = this->GetGeometry();
        const unsigned int NumNodes = rGeom.PointsNumber();
 
        rValue = rN[0] * rGeom[0].FastGetSolutionStepValue(rVariable,Step);
        for (unsigned int i = 1; i < NumNodes; i++)
            rValue += rN[i] * rGeom[i].FastGetSolutionStepValue(rVariable,Step);
    }
 
    virtual void EvaluateConvVelocity(array_1d<double,3>& rConvection,
                                      const ShapeFunctionsType& rN);
 
    virtual void EvaluateConvVelocity(array_1d<double,3>& rConvection,
                                      const array_1d<double,3>& rSubscaleVel,
                                      const ShapeFunctionsType& rN);
 
    virtual void EvaluateViscosity(double& rViscosity,
                                   const ShapeFunctionsType& rN);
 
    virtual void ConvectionOperator(Vector& rResult,
                                    const array_1d<double,3>& rConvVel);
 
 
    virtual double TauOne(const double Density,
                          const double Viscosity,
                          const double ConvVel);
 
 
    virtual double TauTwo(const double Density,
                          const double Viscosity,
                          const double ConvVel);
 
 
    virtual double TauTime(const double Density,
                           const double Viscosity,
                           const double ConvVel,
                           const double Dt);
 
    virtual double InvTauTime(const double Density,
                              const double Viscosity,
                              const double ConvVel,
                              const double Dt);
 
 
    virtual void ASGSMomentumResidual(array_1d<double,3>& rResult,
                                      const double Density,
                                      const array_1d<double,3>& rConvVel,
                                      const ShapeFunctionsType& rN);
 
    virtual void OSSMomentumResidual(array_1d<double,3>& rResult,
                                     const double Density,
                                     const array_1d<double,3>& rConvVel,
                                     const ShapeFunctionsType& rN);
 
    virtual void MassResidual(double& rResult);
 
    virtual void AddViscousTerm(MatrixType& rDampingMatrix,
                                const double Weight,
                                const ShapeDerivativesType& rDN_DX);
 
    virtual void DenseSystemSolve(const Matrix& rA,
                                  const Vector& rB,
                                  Vector& rX);
 
    void EvaluateVorticity(array_1d<double,3>& rVorticity,
                           const ShapeDerivativesType& rDN_DX);
 
 
 
 
 
 
 
 
 
private:
 
    static const double mSubscaleTol;
 
    static const double mSubscaleRHSTol;
 
 
    GeometryData::IntegrationMethod mIntegrationMethod;
 
    ShapeDerivativesType mDN_DX;
 
    double mDetJ;
 
    double mElemSize;
 
    std::vector< array_1d<double,3> > mSubscaleVel;
 
    std::vector< array_1d<double,3> > mOldSubscaleVel;
 
    std::vector< unsigned int > mIterCount;
 
//    std::vector< array_1d<double,3> > mExtraTerm;
 
 
    friend class Serializer;
 
    DynamicVMS();
 
    void save(Serializer& rSerializer) const override;
 
 
    void load(Serializer& rSerializer) override;
 
 
 
 
 
 
 
 
 
 
 
    DynamicVMS & operator=(DynamicVMS const& rOther);
 
    DynamicVMS(DynamicVMS const& rOther);
 
 
}; // Class DynamicVMS
 
 
 
 
 
 
template< unsigned int TDim >
inline std::istream& operator >>(std::istream& rIStream,
                                 DynamicVMS<TDim>& rThis)
{
    return rIStream;
}
 
template< unsigned int TDim >
inline std::ostream& operator <<(std::ostream& rOStream,
                                 const DynamicVMS<TDim>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
 
} // namespace Kratos.
 
#endif // KRATOS_DYNAMIC_VMS__H_INCLUDED  defined