d_vms_dem_coupled.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Joaquin Gonzalez-Usua
//
 
#ifndef KRATOS_D_VMS_DEM_COUPLED_H
#define KRATOS_D_VMS_DEM_COUPLED_H
 
#include "includes/define.h"
#include "includes/element.h"
#include "includes/serializer.h"
#include "geometries/geometry.h"
 
#include "includes/cfd_variables.h"
#include "custom_elements/qs_vms_dem_coupled.h"
#include "custom_elements/d_vms.h"
#include "fluid_dynamics_application_variables.h"
 
namespace Kratos
{
 
 
 
 
 
 
 
template< class TElementData >
class DVMSDEMCoupled : public DVMS<TElementData>
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(DVMSDEMCoupled);
 
    typedef Node NodeType;
 
    typedef Geometry<NodeType> GeometryType;
 
    typedef Geometry<NodeType>::PointsArrayType NodesArrayType;
 
    typedef Vector VectorType;
 
    typedef Matrix MatrixType;
 
    typedef std::size_t IndexType;
 
    typedef std::size_t SizeType;
 
    typedef std::vector<std::size_t> EquationIdVectorType;
 
    typedef std::vector< Dof<double>::Pointer > DofsVectorType;
 
    typedef PointerVectorSet<Dof<double>, IndexedObject> DofsArrayType;
 
    typedef Kratos::Vector ShapeFunctionsType;
 
    typedef Kratos::Matrix ShapeFunctionDerivativesType;
 
    typedef GeometryType::ShapeFunctionsGradientsType ShapeFunctionDerivativesArrayType;
 
    typedef GeometryType::ShapeFunctionsSecondDerivativesType ShapeFunctionsSecondDerivativesType;
 
    constexpr static unsigned int Dim = DVMS<TElementData>::Dim;
    constexpr static unsigned int NumNodes = DVMS<TElementData>::NumNodes;
    constexpr static unsigned int BlockSize = DVMS<TElementData>::BlockSize;
    constexpr static unsigned int LocalSize = DVMS<TElementData>::LocalSize;
    constexpr static unsigned int StrainSize = DVMS<TElementData>::StrainSize;
 
 
    //Constructors.
 
 
    DVMSDEMCoupled(IndexType NewId = 0);
 
 
    DVMSDEMCoupled(IndexType NewId, const NodesArrayType& ThisNodes);
 
 
    DVMSDEMCoupled(IndexType NewId, GeometryType::Pointer pGeometry);
 
 
    DVMSDEMCoupled(IndexType NewId, GeometryType::Pointer pGeometry, Properties::Pointer pProperties);
 
    ~DVMSDEMCoupled();
 
 
 
 
 
 
    Element::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        Properties::Pointer pProperties) const override;
 
 
    Element::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        Properties::Pointer pProperties) const override;
 
    void Initialize(const ProcessInfo& rCurrentProcessInfo) override;
 
    GeometryData::IntegrationMethod GetIntegrationMethod() const 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>& rOutput,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(
        Variable<Matrix> const& rVariable,
        std::vector<Matrix>& rValues,
        ProcessInfo const& rCurrentProcessInfo) override;
 
    void FinalizeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
    void InitializeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) override;
 
    void FinalizeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
    std::string Info() const override;
 
 
    void PrintInfo(std::ostream& rOStream) const override;
 
 
 
 
 
protected:
 
 
    int mInterpolationOrder = 1;
    DenseVector <BoundedMatrix<double,Dim,Dim>> mViscousResistanceTensor;
    // Velocity subscale history, stored at integration points
    DenseVector< array_1d<double,Dim> > mPredictedSubscaleVelocity;
    DenseVector< array_1d<double,Dim> > mOldSubscaleVelocity;
    DenseVector< array_1d<double,Dim> > mPreviousVelocity;
 
 
 
 
    // Protected interface of FluidElement ////////////////////////////////////
 
    void AlgebraicMomentumResidual(
        const TElementData& rData,
        const array_1d<double,3> &rConvectionVelocity,
        array_1d<double,3>& rResidual) const override;
 
    void MomentumProjTerm(
        const TElementData& rData,
        const array_1d<double,3>& rConvectionVelocity,
        array_1d<double,3> &rMomentumRHS) const override;
 
    void AddVelocitySystem(
        TElementData& rData,
        MatrixType& rLocalLHS,
        VectorType& rLocalRHS) override;
 
    void AddReactionStabilization(
        TElementData& rData,
        BoundedMatrix<double,NumNodes*(Dim+1),NumNodes*(Dim+1)>& rLHS,
        VectorType& rLocalRHS);
 
    void AddMassStabilization(
        TElementData& rData,
        MatrixType& rMassMatrix) override;
 
    void CalculateProjections(const ProcessInfo &rCurrentProcessInfo) override;
 
    void CalculateStabilizationParameters(
        const TElementData& rData,
        const array_1d<double,3> &Velocity,
        BoundedMatrix<double,Dim,Dim> &TauOne,
        double &TauTwo) const;
 
    void UpdateIntegrationPointDataSecondDerivatives(
        TElementData& rData,
        unsigned int IntegrationPointIndex,
        double Weight,
        const typename TElementData::MatrixRowType& rN,
        const typename TElementData::ShapeDerivativesType& rDN_DX,
        const typename TElementData::ShapeFunctionsSecondDerivativesType& rDDN_DDX) const;
 
    void SubscaleVelocity(
        const TElementData& rData,
        array_1d<double,3>& rVelocitySubscale) const override;
 
    void SubscalePressure(
        const TElementData& rData,
        double& rPressureSubscale) const override;
 
    void CalculateMassMatrix(MatrixType& rMassMatrix,
                            const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLocalVelocityContribution(MatrixType& rDampMatrix,
                                            VectorType& rRightHandSideVector,
                                            const ProcessInfo& rCurrentProcessInfo) override;
 
    array_1d<double,3> FullConvectiveVelocity(
        const TElementData& rData) const override;
 
    void UpdateSubscaleVelocityPrediction(
        const TElementData& rData) override;
 
    void UpdateSubscaleVelocity(
        const TElementData& rData);
 
    void CalculateResistanceTensor(
        const TElementData& rData);
 
    void AddMassLHS(
        TElementData& rData,
        MatrixType& rMassMatrix) override;
 
    void MassProjTerm(
        const TElementData& rData,
        double& rMassRHS) const override;
 
 
 
 
 
 
 
private:
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
 
 
 
 
 
 
 
    DVMSDEMCoupled& operator=(DVMSDEMCoupled const& rOther);
 
    DVMSDEMCoupled(DVMSDEMCoupled const& rOther);
 
 
 
}; // Class DVMSDEMCoupled
 
 
 
 
 
 
template< class TElementData >
inline std::istream& operator >>(std::istream& rIStream,
                                 DVMSDEMCoupled<TElementData>& rThis)
{
    return rIStream;
}
 
template< class TElementData >
inline std::ostream& operator <<(std::ostream& rOStream,
                                 const DVMSDEMCoupled<TElementData>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
} // namespace Kratos.
 
#endif // KRATOS_D_VMS_DEM_COUPLED_H