spalart_allmaras.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Jordi Cotela, Riccardo Rossi
//
 
#if !defined(KRATOS_SPALART_ALLMARAS_H_INCLUDED )
#define  KRATOS_SPALART_ALLMARAS_H_INCLUDED
 
// System includes
#include <string>
#include <iostream>
 
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "includes/serializer.h"
#include "includes/element.h"
#include "geometries/geometry_data.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
 
class SpalartAllmaras : public Element
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(SpalartAllmaras);
 
    typedef Kratos::Vector ShapeFunctionsType;
 
    typedef Kratos::Matrix ShapeDerivativesType;
 
 
    SpalartAllmaras(IndexType NewId = 0) :
        Element(NewId),
        mIntegrationMethod(GeometryData::IntegrationMethod::GI_GAUSS_1)
    {}
 
    SpalartAllmaras(IndexType NewId, GeometryType::Pointer pGeometry) :
        Element(NewId, pGeometry),
        mIntegrationMethod(pGeometry->GetDefaultIntegrationMethod())
    {}
 
    SpalartAllmaras(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties) :
        Element(NewId, pGeometry, pProperties),
        mIntegrationMethod(pGeometry->GetDefaultIntegrationMethod())
    {}
 
    SpalartAllmaras(IndexType NewId, GeometryType::Pointer pGeometry, const Element::IntegrationMethod& ThisIntegrationMethod) :
        Element(NewId, pGeometry),
        mIntegrationMethod(ThisIntegrationMethod)
    {}
 
    SpalartAllmaras(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties, const Element::IntegrationMethod& ThisIntegrationMethod) :
        Element(NewId,pGeometry,pProperties),
        mIntegrationMethod(ThisIntegrationMethod)
    {}
 
    ~SpalartAllmaras() override
    {}
 
 
 
 
    Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties)  const override;
 
    Element::Pointer Create(IndexType NewId, GeometryType::Pointer pGeom, PropertiesType::Pointer pProperties)  const override;
 
 
    int Check(const ProcessInfo &rCurrentProcessInfo) const override;
 
    void Initialize(const ProcessInfo &rCurrentProcessInfo) override;
 
    void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
//    /// Compute projection of convective term for stabilization
//    virtual void InitializeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo);
 
 
    void CalculateLocalSystem(MatrixType& rLeftHandSideMatrix, VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
//    {
//        KRATOS_THROW_ERROR(std::logic_error, "SplartAllmaras::CalculateRightHandSide method not implemented", "");
//    }
 
    void GetDofList(DofsVectorType& rElementalDofList, const ProcessInfo& rCurrentProcessInfo) const override;
 
    void EquationIdVector(Element::EquationIdVectorType& rResult, const ProcessInfo& rCurrentProcessInfo) const override;
 
    void GetValuesVector(Vector& rValues, int Step = 0) const override;
 
 
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "SpalartAllmaras" << this->GetGeometry().WorkingSpaceDimension() << "D #" << Id();
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "SpalartAllmaras" << this->GetGeometry().WorkingSpaceDimension() << "D #" << Id() << std::endl;
        rOStream << "Number of Nodes: " << this->GetGeometry().PointsNumber() << std::endl;
        rOStream << "Integration method: " << static_cast<int>(this->mIntegrationMethod);
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        this->PrintInfo(rOStream);
        rOStream << "Geometry Data: " << std::endl;
        this->GetGeometry().PrintData(rOStream);
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
    void InitializeElementData();
 
    void AddMassTerm(MatrixType &rMassMatrix,
                     const ShapeFunctionsType& N,
                     const double Weight);
 
    void AddConvection(MatrixType& rLHS,
                       const ShapeFunctionsType& N,
                       const Vector& UGradN,
                       const double Weight);
 
    void AddModelTerms(MatrixType& rLHS,
                       const double MolecularViscosity,
                       const double LastEddyViscosity,
                       const array_1d<double,3>& rLastEddyViscosityGradient,
                       const double Distance,
                       const array_1d<double,3>& rVelocity,
                       const ShapeFunctionsType& N,
                       const ShapeDerivativesType& DN_DX,
                       const double Weight);
 
    void AddStabilization(MatrixType& rLHS,
                          VectorType& rRHS,
                          const double Tau,
                          const ShapeFunctionsType& N,
                          const ShapeDerivativesType& DN_DX,
                          const Vector& UGradN,
                          const double Weight);
 
    void EvaluateConvection(Vector& rResult,
                            const array_1d<double,3>& rConvVel,
                            const ShapeDerivativesType& DN_DX);
 
    template< class TVariableType >
    void EvaluateInPoint(TVariableType& rResult,
                         const Kratos::Variable<TVariableType>& Var,
                         const ShapeFunctionsType& rShapeFunc)
    {
        rResult = rShapeFunc[0] * this->GetGeometry()[0].FastGetSolutionStepValue(Var);
 
        for(SizeType i = 1; i < rShapeFunc.size(); i++)
        {
            rResult += rShapeFunc[i] * this->GetGeometry()[i].FastGetSolutionStepValue(Var);
        }
    }
 
    double MeasureOfVorticity(const ShapeDerivativesType& DN_DX);
 
    void VelocityGradientNorms(double& rNormS,
                               double& rNormOmega,
                               const ShapeDerivativesType& DN_DX);
 
    double CalculateTau(double ElementSize, const ProcessInfo& rCurrentProcessInfo);
 
    double ElementSize();
 
 
 
 
 
 
 
 
 
private:
 
 
 
    Element::IntegrationMethod mIntegrationMethod;
 
    std::vector< ShapeDerivativesType > mDN_DX;
 
 
    double mDetJ;
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_TRY;
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Element );
        //enum IntegrationMethod {GI_GAUSS_1, GI_GAUSS_2, GI_GAUSS_3, GI_GAUSS_4, GI_GAUSS_5,NumberOfIntegrationMethods};
        unsigned int IntMethod = 0;
        switch(mIntegrationMethod)
        {
        case GeometryData::IntegrationMethod::GI_GAUSS_1:
            IntMethod = 1;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_2:
            IntMethod = 2;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_3:
            IntMethod = 3;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_4:
            IntMethod = 4;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_5:
            IntMethod = 5;
            break;
        default:
            KRATOS_ERROR << "Unknown integration method encountered on serializer save for SpalartAllmaras element: " << static_cast<int>(mIntegrationMethod) << std::endl;
            break;
        }
        rSerializer.save("IntMethod",IntMethod);
        rSerializer.save("mDN_DX",mDN_DX);
        rSerializer.save("mDetJ",mDetJ);
        KRATOS_CATCH("");
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_TRY;
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer,Element);
        unsigned int IntMethod = 0;
        rSerializer.load("IntMethod",IntMethod);
        switch(mIntegrationMethod)
        {
        case GeometryData::IntegrationMethod::GI_GAUSS_1:
            mIntegrationMethod = GeometryData::IntegrationMethod::GI_GAUSS_1;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_2:
            mIntegrationMethod = GeometryData::IntegrationMethod::GI_GAUSS_2;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_3:
            mIntegrationMethod = GeometryData::IntegrationMethod::GI_GAUSS_3;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_4:
            mIntegrationMethod = GeometryData::IntegrationMethod::GI_GAUSS_4;
            break;
        case GeometryData::IntegrationMethod::GI_GAUSS_5:
            mIntegrationMethod = GeometryData::IntegrationMethod::GI_GAUSS_5;
            break;
        default:
            KRATOS_ERROR << "Unknown integration method encountered on serializer load for SpalartAllmaras element: " << static_cast<int>(IntMethod) << std::endl;
            break;
        }
        rSerializer.load("mDN_DX",mDN_DX);
        rSerializer.load("mDetJ",mDetJ);
        KRATOS_CATCH("");
    }
 
 
 
 
 
 
 
 
 
 
    SpalartAllmaras& operator=(SpalartAllmaras const& rOther);
 
    SpalartAllmaras(SpalartAllmaras const& rOther);
 
 
 
}; // Class SpalartAllmaras
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  SpalartAllmaras& rThis)
{
    return rIStream;
}
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const SpalartAllmaras& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_SPALART_ALLMARAS_H_INCLUDED  defined