embedded_transonic_perturbation_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 Nuñez, Eloisa Baez Jones, Inigo Lopez and Riccardo Rossi
//
 
#if !defined(KRATOS_EMBEDDED_TRANSONIC_PERTURBATION_POTENTIAL_FLOW_ELEMENT_H)
#define KRATOS_EMBEDDED_TRANSONIC_PERTURBATION_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 "transonic_perturbation_potential_flow_element.h"
#include "custom_utilities/potential_flow_utilities.h"
 
namespace Kratos
{
 
 
 
 
 
template <int TDim, int TNumNodes>
class EmbeddedTransonicPerturbationPotentialFlowElement : public TransonicPerturbationPotentialFlowElement<TDim,TNumNodes>
{
public:
 
    typedef TransonicPerturbationPotentialFlowElement<TDim,TNumNodes> BaseType;
 
    typedef typename BaseType::SizeType SizeType;
    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;
    typedef typename BaseType::GeometriesArrayType GeometriesArrayType;
    typedef PotentialFlowUtilities::ElementalData<TNumNodes, TDim> ElementalData;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(EmbeddedTransonicPerturbationPotentialFlowElement);
 
 
    // Constructors.
 
 
    explicit EmbeddedTransonicPerturbationPotentialFlowElement(IndexType NewId = 0)
    {
    }
 
    EmbeddedTransonicPerturbationPotentialFlowElement(IndexType NewId,
                                        const NodesArrayType& ThisNodes)
        : TransonicPerturbationPotentialFlowElement<TDim,TNumNodes>(NewId, ThisNodes)
    {
    }
 
    EmbeddedTransonicPerturbationPotentialFlowElement(IndexType NewId,
                                        typename GeometryType::Pointer pGeometry)
        : TransonicPerturbationPotentialFlowElement<TDim,TNumNodes>(NewId, pGeometry)
    {
    }
 
    EmbeddedTransonicPerturbationPotentialFlowElement(IndexType NewId,
                                     typename GeometryType::Pointer pGeometry,
                                     typename PropertiesType::Pointer pProperties)
        : TransonicPerturbationPotentialFlowElement<TDim,TNumNodes>(NewId, pGeometry, pProperties)
    {
    }
 
    EmbeddedTransonicPerturbationPotentialFlowElement(EmbeddedTransonicPerturbationPotentialFlowElement const& rOther) = delete;
 
    EmbeddedTransonicPerturbationPotentialFlowElement(EmbeddedTransonicPerturbationPotentialFlowElement&& rOther) = delete;
 
    ~EmbeddedTransonicPerturbationPotentialFlowElement() override
    {
    }
 
 
    EmbeddedTransonicPerturbationPotentialFlowElement& operator=(EmbeddedTransonicPerturbationPotentialFlowElement const& rOther) = delete;
 
    EmbeddedTransonicPerturbationPotentialFlowElement& operator=(EmbeddedTransonicPerturbationPotentialFlowElement&& 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 Initialize(const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(VectorType& rRightHandSideVector,
                                const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLeftHandSide(MatrixType& rLeftHandSideMatrix,
                               const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculateOnIntegrationPoints(const Variable<double>& rVariable,
                                     std::vector<double>& rValues,
                                     const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(const Variable<int>& rVariable,
                                     std::vector<int>& rValues,
                                     const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(const Variable<array_1d<double, 3>>& rVariable,
                                     std::vector<array_1d<double, 3>>& rValues,
                                     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:
 
 
    void CalculateLeftHandSideKuttaWakeElement(MatrixType& rLeftHandSideMatrix,
                                          const ProcessInfo& rCurrentProcessInfo);
 
    void CalculateRightHandSideKuttaWakeElement(VectorType& rRightHandSideVector,
                                          const ProcessInfo& rCurrentProcessInfo);
 
    void CalculateLeftHandSideContribution(BoundedMatrix<double, TNumNodes, TNumNodes>& rLhs_total,
                                         const ProcessInfo& rCurrentProcessInfo,
                                         const array_1d<double, TDim>& rVelocity,
                                         const ElementalData& rData) override;
 
    void CalculateRightHandSideContribution(BoundedVector<double, TNumNodes>& rRhs_total,
                                        const double rDensity,
                                        const array_1d<double, TDim>& rVelocity) override;
 
 
    void AssembleSupersonicLeftHandSide(MatrixType& rLeftHandSideMatrix,
                                        const double densityDerivativeWRTVelocity,
                                        const double densityDerivativeWRTUpwindVelocity,
                                        const array_1d<double, TDim> velocity,
                                        const array_1d<double, TDim> upwindVelocity,
                                        const ProcessInfo& rCurrentProcessInfo) override;
 
    ModifiedShapeFunctions::Pointer pGetModifiedShapeFunctions(Vector& rDistances);
 
    void FindUpwindElement(const ProcessInfo& rCurrentProcessInfo) override;
 
private:
 
    void SelectUpwindElement(std::vector<IndexType>& rUpwindElementNodesIds,
                            GlobalPointersVector<Element>& rUpwindElementCandidates,
                            const ProcessInfo& rCurrentProcessInfo);
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
 
 
 
}; // Class EmbeddedTransonicPerturbationPotentialFlowElement
 
 
 
 
 
} // namespace Kratos.
 
#endif // KRATOS_TRANSONIC_PERTURBATION_POTENTIAL_FLOW_ELEMENT_H  defined