parabolic_profile_utilities.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Eduardo Soudah
//                   Ruben Zorilla
//
 
#if !defined(KRATOS_PARABOLIC_PROFILE_UTILITIES_H )
#define  KRATOS_PARABOLIC_PROFILE_UTILITIES_H
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "utilities/function_parser_utility.h"
 
// Application includes
 
namespace Kratos
{
 
 
 
class KRATOS_API(FLUID_DYNAMICS_BIOMEDICAL_APPLICATION) ParabolicProfileUtilities
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(ParabolicProfileUtilities);
 
    using NodeType = typename ModelPart::NodeType;
 
 
    // Deleted default constructor
    ParabolicProfileUtilities() = delete;
 
    ParabolicProfileUtilities(ParabolicProfileUtilities const& rOther) = delete;
 
    ~ParabolicProfileUtilities() = default;
 
 
    ParabolicProfileUtilities& operator=(ParabolicProfileUtilities const& rOther) = delete;
 
 
    static double CalculateInletArea(const ModelPart& rModelPart);
 
    static ModelPart& CreateAndFillInletAuxiliaryVolumeModelPart(ModelPart& rInletModelPart);
 
    static void CalculateWallParallelDistance(
        ModelPart& rWallModelPart,
        ModelPart& rFluidModelPart,
        const std::size_t WallDistanceLevels);
 
    static void ImposeParabolicInlet(
        ModelPart &rModelPart,
        const double MaxParabolaValue,
        const double MaxValueFactor = 1.0);
 
    static void ImposeParabolicInlet(
        ModelPart &rModelPart,
        const GenericFunctionUtility::Pointer rMaxParabolaValue,
        const double MaxValueFactor = 1.0);
 
    static void FreeParabolicInlet(ModelPart& rModelPart);
 
private:
 
    template<class TInputType>
    static void ImposeParabolicProfile(
        ModelPart &rModelPart,
        const TInputType& rMaxParabolaValue,
        const double MaxValueFactor);
 
    template<class TInputType>
    static double GetMaxParabolaValue(
        const double Time,
        const NodeType& rNode,
        TInputType& rMaxParabolaValue);
 
    static double CalculateBoundingBoxCharacteristicLength(const ModelPart& rModelPart);
 
};  // Class ParabolicProfileUtilities
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_PARABOLIC_PROFILE_UTILITIES_H  defined