fluid_element_utilities.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Zorrilla
//                   Jordi Cotela
//
 
#if !defined(KRATOS_FLUID_ELEMENT_UTILITIES_H )
#define  KRATOS_FLUID_ELEMENT_UTILITIES_H
 
// External includes
 
// Project includes
#include "includes/define.h"
 
// Application includes
#include "custom_utilities/fluid_element_data.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
template< std::size_t TNumNodes >
class FluidElementUtilities {
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(FluidElementUtilities);
 
    using ShapeDerivatives2DType = typename FluidElementData<2,TNumNodes,false>::ShapeDerivativesType;
    using ShapeDerivatives3DType = typename FluidElementData<3,TNumNodes,false>::ShapeDerivativesType;
 
    constexpr static std::size_t VoigtVector2DSize = FluidElementData<2,TNumNodes,false>::StrainSize;
    constexpr static std::size_t VoigtVector3DSize = FluidElementData<3,TNumNodes,false>::StrainSize;
 
 
    // Deleted default constructor
    FluidElementUtilities() = delete;
 
    FluidElementUtilities(FluidElementUtilities const& rOther) = delete;
 
    ~FluidElementUtilities();
 
 
    FluidElementUtilities& operator=(FluidElementUtilities const& rOther) = delete;
 
 
    static void GetStrainMatrix(
        const ShapeDerivatives2DType& rDNDX,
        BoundedMatrix<double, VoigtVector2DSize, 3 * TNumNodes>& rStrainMatrix);
 
    static void GetStrainMatrix(
        const ShapeDerivatives3DType& rDNDX,
        BoundedMatrix<double, VoigtVector3DSize, 4 * TNumNodes>& rStrainMatrix);
 
    static void GetNewtonianConstitutiveMatrix(
        const double DynamicViscosity,
        BoundedMatrix<double, VoigtVector2DSize, VoigtVector2DSize>& rConstitutiveMatrix);
 
    static void GetNewtonianConstitutiveMatrix(
        const double DynamicViscosity,
        BoundedMatrix<double, VoigtVector3DSize, VoigtVector3DSize>& rConstitutiveMatrix);
 
    static void VoigtTransformForProduct(
        const array_1d<double,3>& rVector,
        BoundedMatrix<double, 2, VoigtVector2DSize>& rVoigtMatrix);
 
    static void VoigtTransformForProduct(
        const array_1d<double,3>& rVector,
        BoundedMatrix<double, 3, VoigtVector3DSize>& rVoigtMatrix);
 
    static void SetNormalProjectionMatrix(
        const array_1d<double, 3>& rUnitNormal,
        BoundedMatrix<double, 2, 2>& rNormalProjMatrix);
 
    static void SetNormalProjectionMatrix(
        const array_1d<double, 3>& rUnitNormal,
        BoundedMatrix<double, 3, 3>& rNormalProjMatrix);
 
    static void SetTangentialProjectionMatrix(
        const array_1d<double, 3>& rUnitNormal,
        BoundedMatrix<double, 2, 2>& rTangProjMatrix);
 
    static void SetTangentialProjectionMatrix(
        const array_1d<double, 3>& rUnitNormal,
        BoundedMatrix<double, 3, 3>& rTangProjMatrix);
 
    static void DenseSystemSolve(
        const BoundedMatrix<double,2,2> &rA,
        const array_1d<double,2> &rB,
        array_1d<double,2> &rX);
 
    static void DenseSystemSolve(
        const BoundedMatrix<double,3,3> &rA,
        const array_1d<double,3> &rB,
        array_1d<double,3> &rX);
 
 
};  // Class FluidElementUtilities
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_FLUID_ELEMENT_UTILITIES_H  defined