define_2d_wake_process.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Inigo Lopez
//
 
#if !defined(KRATOS_DEFINE_2D_WAKE_PROCESS_H_INCLUDED )
#define  KRATOS_DEFINE_2D_WAKE_PROCESS_H_INCLUDED
 
// Project includes
#include "includes/model_part.h"
#include "processes/process.h"
#include "compressible_potential_flow_application_variables.h"
 
namespace Kratos
{
 
 
class KRATOS_API(COMPRESSIBLE_POTENTIAL_FLOW_APPLICATION) Define2DWakeProcess
    : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(Define2DWakeProcess);
 
    typedef Node NodeType;
 
 
    Define2DWakeProcess(ModelPart& rBodyModelPart,const double Tolerance);
 
    Define2DWakeProcess(Define2DWakeProcess const& rOther) = delete;
 
    ~Define2DWakeProcess() override = default;
 
 
    Define2DWakeProcess& operator=(Define2DWakeProcess const& rOther) = delete;
 
 
    void ExecuteInitialize() override;
 
 
    std::string Info() const override
    {
        return "Define2DWakeProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "Define2DWakeProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
private:
 
    // The airfoil model part conatining the trailing edge
    ModelPart& mrBodyModelPart;
    // Tolerance to avoid nodes laying exactly on the wake
    const double mTolerance;
    NodeType* mpTrailingEdgeNode;
    BoundedVector<double, 3> mWakeDirection;
    BoundedVector<double, 3> mWakeNormal;
    // Vector to store the trailing edge elements ids
    std::vector<std::size_t> mTrailingEdgeElementsOrderedIds;
    void InitializeTrailingEdgeSubModelpart() const;
 
    void InitializeWakeSubModelpart() const;
 
    void SetWakeDirectionAndNormal();
 
    void SaveTrailingEdgeNode();
 
    void MarkWakeElements();
 
    void CheckIfTrailingEdgeElement(Element& rElement);
 
    bool CheckIfPotentiallyWakeElement(const Element& rElement) const;
 
    const BoundedVector<double, 3> ComputeNodalDistancesToWake(const Element& rElement) const;
 
    void AddTrailingEdgeAndWakeElements(std::vector<std::size_t>& rWakeElementsOrderedIds);
 
    void MarkKuttaElements() const;
 
    void MarkWakeTrailingEdgeElement() const;
 
    bool CheckIfTrailingEdgeElementIsCutByWake(const Element& rElement) const;
 
    const BoundedVector<double, 3> ComputeDistanceFromTrailingEdgeToPoint(const Point& rInputPoint) const;
 
}; // Class Define2DWakeProcess
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  Define2DWakeProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const Define2DWakeProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
}  // namespace Kratos.
 
#endif // KRATOS_DEFINE_2D_WAKE_PROCESS_H_INCLUDED  defined