fluid_constitutive_law.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Jordi Cotela
//
 
#if !defined (KRATOS_FLUID_CONSTITUTIVE_LAW)
#define  KRATOS_FLUID_CONSTITUTIVE_LAW
 
// System includes
 
// External includes
 
// Project includes
#include "includes/constitutive_law.h"
 
 
namespace Kratos
{
 
 
 
class KRATOS_API(FLUID_DYNAMICS_APPLICATION) FluidConstitutiveLaw : public ConstitutiveLaw
{
public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION(FluidConstitutiveLaw);
 
    typedef ConstitutiveLaw         BaseType;
    typedef std::size_t             SizeType;
 
 
    FluidConstitutiveLaw();
 
    FluidConstitutiveLaw (const FluidConstitutiveLaw& rOther);
 
    ~FluidConstitutiveLaw() override;
 
 
    ConstitutiveLaw::Pointer Clone() const override;
 
 
    void CalculateMaterialResponseCauchy(ConstitutiveLaw::Parameters& rValues) override;
 
 
    int Check(
        const Properties& rMaterialProperties,
        const GeometryType& rElementGeometry,
        const ProcessInfo& rCurrentProcessInfo) const override;
 
 
    int& CalculateValue(ConstitutiveLaw::Parameters& rParameters, const Variable<int>& rThisVariable, int& rValue) override;
 
 
    double& CalculateValue(ConstitutiveLaw::Parameters& rParameters, const Variable<double>& rThisVariable, double& rValue) override;
 
    Vector& CalculateValue(ConstitutiveLaw::Parameters& rParameters, const Variable<Vector>& rThisVariable, Vector& rValue) override;
 
    Matrix& CalculateValue(ConstitutiveLaw::Parameters& rParameters, const Variable<Matrix>& rThisVariable, Matrix& rValue) override;
 
    array_1d<double, 3 > & CalculateValue(ConstitutiveLaw::Parameters& rParameters, const Variable<array_1d<double, 3 > >& rThisVariable,array_1d<double, 3 > & rValue) override;
 
    array_1d<double, 6 > & CalculateValue(ConstitutiveLaw::Parameters& rParameters, const Variable<array_1d<double, 6 > >& rThisVariable, array_1d<double, 6 > & rValue) override;
 
 
 
    SizeType WorkingSpaceDimension() override;
 
 
    SizeType GetStrainSize() const override;
 
 
    std::string Info() const override;
 
    void PrintInfo(std::ostream& rOStream) const override;
 
    void PrintData(std::ostream& rOStream) const override;
 
 
protected:
 
 
 
    void NewtonianConstitutiveMatrix2D(double EffectiveViscosity, Matrix& rC);
 
 
    void NewtonianConstitutiveMatrix3D(double EffectiveViscosity, Matrix& rC);
 
 
    virtual double GetEffectiveViscosity(ConstitutiveLaw::Parameters& rParameters) const;
 
    virtual double GetValueFromTable(
        const Variable<double> &rVariableInput,
        const Variable<double> &rVariableOutput,
        ConstitutiveLaw::Parameters &rParameters) const;
 
 
    private:
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
}; // Class FluidConstitutiveLaw
 
}  // namespace Kratos.
 
#endif // KRATOS_FLUID_CONSTITUTIVE_LAW  defined