embedded_laplacian_element.h

Go to the documentation of this file.

// KRATOS ___ ___  _  ___   __   ___ ___ ___ ___
//       / __/ _ \| \| \ \ / /__|   \_ _| __| __|
//      | (_| (_) | .` |\ V /___| |) | || _|| _|
//       \___\___/|_|\_| \_/    |___/___|_| |_|  APPLICATION
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Franziska Wahl
//
 
#ifndef KRATOS_EMBEDDED_LAPLACIAN_ELEMENT_H_INCLUDED
#define  KRATOS_EMBEDDED_LAPLACIAN_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/element.h"
#include "modified_shape_functions/modified_shape_functions.h"
 
// Application includes
#include "laplacian_element.h"
 
namespace Kratos
{
 
 
 
 
 
 
namespace EmbeddedLaplacianInternals {
 
    template <size_t TDim, size_t TNumNodes>
    ModifiedShapeFunctions::Pointer GetContinuousShapeFunctionCalculator(
        const Element &rElement,
        const Vector &rNodalDistances);
 
    // Data container class
    template<std::size_t TDim>
    class EmbeddedElementData
    {
    public:
 
        static constexpr std::size_t NumNodes = TDim + 1;
 
        typedef GeometryData::ShapeFunctionsGradientsType ShapeFunctionsGradientsType;
        typedef std::vector<array_1d<double,3>> InterfaceNormalsType;
        typedef array_1d<double,NumNodes> NodalScalarData;
 
 
        double PenaltyCoefficient;
 
        NodalScalarData NodalDistances;
 
        Matrix PositiveSideN;
        ShapeFunctionsGradientsType PositiveSideDNDX;
        Vector PositiveSideWeights;
 
        Matrix PositiveInterfaceN;
        ShapeFunctionsGradientsType PositiveInterfaceDNDX;
        Vector PositiveInterfaceWeights;
        InterfaceNormalsType PositiveInterfaceUnitNormals;
 
        size_t NumPositiveNodes;
        size_t NumNegativeNodes;
 
 
        void Initialize(const Element& rElement);
 
        bool IsSplit();
 
    };
 
} //namespace EmbeddedLaplacianInternals
 
 
template<std::size_t TDim>
class EmbeddedLaplacianElement : public LaplacianElement
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(EmbeddedLaplacianElement);
 
    typedef LaplacianElement BaseType;
    typedef GeometryData::ShapeFunctionsGradientsType ShapeFunctionsGradientsType;
 
    static constexpr std::size_t NumNodes = TDim + 1;
 
    using NodalScalarData = array_1d<double,NumNodes>;
    using EmbeddedElementData = EmbeddedLaplacianInternals::EmbeddedElementData < TDim >;
 
 
    EmbeddedLaplacianElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry);
 
    EmbeddedLaplacianElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry,
        PropertiesType::Pointer pProperties);
 
    virtual ~EmbeddedLaplacianElement();
 
 
 
 
    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;
 
    void CalculateLocalSystem(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLeftHandSide(
        MatrixType& rLeftHandSideMatrix,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
 
 
 
 
 
 
 
 
 
protected:
 
 
 
 
 
    void InitializeGeometryData(
        EmbeddedElementData& rData);
 
    void NormalizeInterfaceNormals(
        typename EmbeddedElementData::InterfaceNormalsType& rNormals,
        double Tolerance) const;
 
    void AddPositiveElementSide(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo,
        const EmbeddedElementData& rData);
 
    void AddPositiveInterfaceTerms(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo,
        const EmbeddedElementData& rData);
 
    void AddNitscheBoundaryTerms(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo,
        const EmbeddedElementData& rData);
 
 
 
 
 
 
    // Protected default constructor necessary for serialization
    EmbeddedLaplacianElement() : LaplacianElement()
    {
    }
 
 
private:
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, LaplacianElement);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, LaplacianElement);
    }
 
 
 
 
 
 
 
 
 
 
    //EmbeddedLaplacianElement& operator=(const EmbeddedLaplacianElement& rOther);
 
    //EmbeddedLaplacianElement(const EmbeddedLaplacianElement& rOther);
 
 
}; // Class EmbeddedLaplacianElement
 
 
 
 
 
 
/*  inline std::istream& operator >> (std::istream& rIStream,
                    EmbeddedLaplacianElement& rThis);
*/
/*  inline std::ostream& operator << (std::ostream& rOStream,
                    const EmbeddedLaplacianElement& rThis)
    {
      rThis.PrintInfo(rOStream);
      rOStream << std::endl;
      rThis.PrintData(rOStream);
 
      return rOStream;
    }*/
 
}  // namespace Kratos.
 
#endif // KRATOS_EMBEDDED_LAPLACIAN_ELEMENT_H_INCLUDED  defined