fluid_characteristic_numbers_utilities.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Zorrilla
//
//
 
#ifndef KRATOS_FLUID_CHARACTERISTIC_NUMBERS_UTILITIES_H
#define KRATOS_FLUID_CHARACTERISTIC_NUMBERS_UTILITIES_H
 
// System includes
#include <utility>
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "includes/kratos_parameters.h"
 
// Application includes
 
 
namespace Kratos
{
 
 
 
class KRATOS_API(FLUID_DYNAMICS_APPLICATION) FluidCharacteristicNumbersUtilities
{
public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION(FluidCharacteristicNumbersUtilities);
 
    typedef std::function<double(const Geometry<Node>&)> ElementSizeFunctionType;
 
 
    FluidCharacteristicNumbersUtilities() = delete;
 
 
    static void CalculateLocalCFL(ModelPart& rModelPart);
 
    static double CalculateElementCFL(
        const Element &rElement,
        const ElementSizeFunctionType& rElementSizeCalculator,
        const double Dt);
 
 
    static double CalculateElementCFLWithSoundVelocity(
        const Element &rElement,
        const ElementSizeFunctionType& rElementSizeCalculator,
        const double Dt);
 
 
    template<bool ConsiderArtificialMagnitudes>
    static double CalculateElementPrandtlNumber(const Element& rElement);
 
    template<bool ConsiderArtificialMagnitudes, bool DensityIsNodal>
    static std::pair<double,double> CalculateElementPecletNumbers(
        const Element& rElement,
        const ElementSizeFunctionType& rElementSizeCalculator);
 
    template<bool ConsiderArtificialMagnitudes, bool DensityIsNodal>
    static double CalculateElementViscousPecletNumber(
        const Element& rElement,
        const ElementSizeFunctionType& rElementSizeCalculator);
 
    template<bool ConsiderArtificialMagnitudes, bool DensityIsNodal>
    static double CalculateElementThermalPecletNumber(
        const Element& rElement,
        const ElementSizeFunctionType& rElementSizeCalculator);
 
    template<bool ConsiderArtificialMagnitudes, bool DensityIsNodal>
    static std::pair<double,double> CalculateElementFourierNumbers(
        const Element& rElement,
        const ElementSizeFunctionType& rElementSizeCalculator,
        const double Dt);
 
    template<bool ConsiderArtificialMagnitudes, bool DensityIsNodal>
    static double CalculateElementViscousFourierNumber(
        const Element& rElement,
        const ElementSizeFunctionType& rElementSizeCalculator,
        const double Dt);
 
    template<bool ConsiderArtificialMagnitudes, bool DensityIsNodal>
    static double CalculateElementThermalFourierNumber(
        const Element& rElement,
        const ElementSizeFunctionType& rElementSizeCalculator,
        const double Dt);
 
    static double CalculateElementMachNumber(const Element& rElement);
 
    static ElementSizeFunctionType GetMinimumElementSizeFunction(const Geometry<Node>& rGeometry);
 
    static ElementSizeFunctionType GetAverageElementSizeFunction(const Geometry<Node>& rGeometry);
 
 
private:
 
 
    template<bool IsNodal>
    static double AuxiliaryGetDensity(const Element& rElement);
 
    template<bool AddArtificialValues>
    static std::pair<double,double> GetDiffusivityValues(const Element& rElement);
 
    template<bool AddArtificialValues>
    static double GetDynamicViscosityValue(const Element& rElement);
 
    template<bool AddArtificialValues>
    static double GetConductivityValue(const Element& rElement);
 
};
 
 
 
} // namespace Kratos.
 
 
#endif  /* KRATOS_FLUID_CHARACTERISTIC_NUMBERS_UTILITIES_H */