embedded_compressible_potential_flow_element.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:        BSD License
//                  Kratos default license: kratos/license.txt
//
//  Main authors:    Marc Nunez, based on Inigo Lopez and Riccardo Rossi work
//
 
#if !defined(KRATOS_EMBEDDED_COMPRESSIBLE_POTENTIAL_FLOW_ELEMENT_H)
#define KRATOS_EMBEDDED_COMPRESSIBLE_POTENTIAL_FLOW_ELEMENT_H
 
// System includes
 
// External includes
 
// Project includes
#include "modified_shape_functions/triangle_2d_3_modified_shape_functions.h"
#include "modified_shape_functions/tetrahedra_3d_4_modified_shape_functions.h"
#include "compressible_potential_flow_element.h"
 
namespace Kratos
{
template <int Dim, int NumNodes>
class EmbeddedCompressiblePotentialFlowElement : public CompressiblePotentialFlowElement<Dim,NumNodes>
{
public:
 
    typedef CompressiblePotentialFlowElement<Dim,NumNodes> BaseType;
 
    typedef typename BaseType::IndexType IndexType;
    typedef typename BaseType::GeometryType GeometryType;
    typedef typename BaseType::PropertiesType PropertiesType;
    typedef typename BaseType::NodesArrayType NodesArrayType;
    typedef typename BaseType::VectorType VectorType;
    typedef typename BaseType::MatrixType MatrixType;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(EmbeddedCompressiblePotentialFlowElement);
 
 
    // Constructors.
 
 
    explicit EmbeddedCompressiblePotentialFlowElement(IndexType NewId = 0){}
 
    EmbeddedCompressiblePotentialFlowElement(IndexType NewId, const NodesArrayType& ThisNodes)
        : CompressiblePotentialFlowElement<Dim,NumNodes>(NewId, ThisNodes){}
 
    EmbeddedCompressiblePotentialFlowElement(IndexType NewId, typename GeometryType::Pointer pGeometry)
        : CompressiblePotentialFlowElement<Dim,NumNodes>(NewId, pGeometry){}
 
    EmbeddedCompressiblePotentialFlowElement(IndexType NewId,
                                       typename GeometryType::Pointer pGeometry,
                                       typename PropertiesType::Pointer pProperties)
        : CompressiblePotentialFlowElement<Dim,NumNodes>(NewId, pGeometry, pProperties){}
 
    EmbeddedCompressiblePotentialFlowElement(EmbeddedCompressiblePotentialFlowElement const& rOther) = delete;
 
    EmbeddedCompressiblePotentialFlowElement(EmbeddedCompressiblePotentialFlowElement&& rOther) = delete;
 
    ~EmbeddedCompressiblePotentialFlowElement() override{}
 
 
    EmbeddedCompressiblePotentialFlowElement& operator=(EmbeddedCompressiblePotentialFlowElement const& rOther) = delete;
 
    EmbeddedCompressiblePotentialFlowElement& operator=(EmbeddedCompressiblePotentialFlowElement&& rOther) = delete;
 
    Element::Pointer Create(IndexType NewId,
                            NodesArrayType const& ThisNodes,
                            typename PropertiesType::Pointer pProperties) const override;
 
    Element::Pointer Create(IndexType NewId,
                            typename GeometryType::Pointer pGeom,
                            typename PropertiesType::Pointer pProperties) const override;
 
    Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
 
    void CalculateLocalSystem(MatrixType& rLeftHandSideMatrix,
                              VectorType& rRightHandSideVector,
                              const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(VectorType& rRightHandSideVector,
                                const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLeftHandSide(MatrixType& rLeftHandSideMatrix,
                               const ProcessInfo& rCurrentProcessInfo) override;
 
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
    std::string Info() const override;
 
    void PrintInfo(std::ostream& rOStream) const override;
 
    void PrintData(std::ostream& rOStream) const override;
 
 
protected:
 
    ModifiedShapeFunctions::Pointer pGetModifiedShapeFunctions(Vector& rDistances);
 
private:
 
    void CalculateEmbeddedLocalSystem(MatrixType& rLeftHandSideMatrix,
                              VectorType& rRightHandSideVector,
                              const ProcessInfo& rCurrentProcessInfo);
 
    void CalculateKuttaWakeLocalSystem(MatrixType& rLeftHandSideMatrix,
                              VectorType& rRightHandSideVector,
                              const ProcessInfo& rCurrentProcessInfo);
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
}; // Class EmbeddedCompressiblePotentialFlowElement
 
} // namespace Kratos.
 
#endif // KRATOS_EMBEDDED_INCOMPRESSIBLE_POTENTIAL_FLOW_ELEMENT_H  defined