flux_condition.h

Go to the documentation of this file.

// KRATOS ___ ___  _  ___   __   ___ ___ ___ ___
//       / __/ _ \| \| \ \ / /__|   \_ _| __| __|
//      | (_| (_) | .` |\ V /___| |) | || _|| _|
//       \___\___/|_|\_| \_/    |___/___|_| |_|  APPLICATION
//
//  License: BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Jordi Cotela
//
 
#if !defined(KRATOS_FLUX_CONDITION_H_INCLUDED )
#define  KRATOS_FLUX_CONDITION_H_INCLUDED
 
// System includes
#include <string>
#include <iostream>
 
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "includes/node.h"
#include "includes/process_info.h"
#include "includes/properties.h"
#include "includes/condition.h"
#include "geometries/geometry.h"
#include "includes/variables.h"
 
namespace Kratos
{
 
 
namespace FluxConditionInternals
{
 
template< unsigned int TNodeNumber >
class IntegrationData
{
public:
 
    IntegrationData(
        Geometry< Node >& rGeometry,
        const Variable<double>& rFluxVar
        ):
        mGaussPoint(0),
        mNodalFluxes(TNodeNumber,0.0)
    {
        NumGauss = rGeometry.IntegrationPointsNumber(GeometryData::IntegrationMethod::GI_GAUSS_2);
        Vector DetJ = ZeroVector(NumGauss);
        rGeometry.DeterminantOfJacobian(DetJ,GeometryData::IntegrationMethod::GI_GAUSS_2);
 
        mShapeFunctionValues.resize(NumGauss,TNodeNumber);
        mIntegrationWeights.resize(NumGauss);
 
        noalias(mShapeFunctionValues) = rGeometry.ShapeFunctionsValues(GeometryData::IntegrationMethod::GI_GAUSS_2);
 
        const auto& IntegrationPoints = rGeometry.IntegrationPoints(GeometryData::IntegrationMethod::GI_GAUSS_2);
 
        for (unsigned int g = 0; g < NumGauss; g++)
        {
            mIntegrationWeights[g] = DetJ[g]*IntegrationPoints[g].Weight();
        }
 
        for (unsigned int i = 0; i < TNodeNumber; i++)
        {
            mNodalFluxes[i] = rGeometry[i].FastGetSolutionStepValue(rFluxVar);
        }
    }
 
    void SetCurrentGaussPoint(unsigned int g)
    {
        mGaussPoint = g;
    }
 
    double N(unsigned int i) const
    {
        return mShapeFunctionValues(mGaussPoint,i);
    }
 
    double NodalFlux(unsigned int i) const
    {
        return mNodalFluxes[i];
    }
 
    double GaussPointFlux() const
    {
        double flux = mNodalFluxes[0]*mShapeFunctionValues(mGaussPoint,0);
        for (unsigned int i = 1; i < TNodeNumber; i++)
        {
            flux += mNodalFluxes[i]*mShapeFunctionValues(mGaussPoint,i);
        }
        return flux;
    }
 
    double IntegrationWeight() const
    {
        return mIntegrationWeights[mGaussPoint];
    }
 
    unsigned int NumGauss;
 
private:
 
    unsigned int mGaussPoint;
 
    array_1d<double,TNodeNumber> mNodalFluxes;
 
    Matrix mShapeFunctionValues;
 
    Vector mIntegrationWeights;
};
 
 
}
 
 
 
template< unsigned int TNodeNumber >
class FluxCondition: public Condition
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(FluxCondition);
 
    typedef Condition::MatrixType MatrixType;
    typedef Condition::VectorType VectorType;
 
 
    FluxCondition(
        IndexType NewId,
        Geometry< Node >::Pointer pGeometry);
 
    FluxCondition(
        IndexType NewId,
        Geometry< Node >::Pointer pGeometry,
        Properties::Pointer pProperties);
 
    ~FluxCondition() override;
 
 
    Condition::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        Properties::Pointer pProperties) const override;
 
    Condition::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        Properties::Pointer pProperties) const override;
 
    void CalculateLocalSystem(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void EquationIdVector(
        EquationIdVectorType& rResult,
        const ProcessInfo& rCurrentProcessInfo) const override;
 
    void GetDofList(
        DofsVectorType& ConditionalDofList,
        const ProcessInfo& CurrentProcessInfo) const override;
 
    GeometryData::IntegrationMethod GetIntegrationMethod() const override;
 
    void CalculateOnIntegrationPoints(
        const Variable<array_1d<double, 3 > >& rVariable,
        std::vector<array_1d<double, 3 > >& rValues,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(
        const Variable<double>& rVariable,
        std::vector<double>& rValues,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(
        const Variable<array_1d<double, 6 > >& rVariable,
        std::vector<array_1d<double, 6 > >& rValues,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(
        const Variable<Vector>& rVariable,
        std::vector<Vector>& rValues,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateOnIntegrationPoints(
        const Variable<Matrix>& rVariable,
        std::vector<Matrix>& rValues,
        const ProcessInfo& rCurrentProcessInfo) override;
 
 
    std::string Info() const override;
 
    void PrintInfo(std::ostream& rOStream) const override;
 
    void PrintData(std::ostream& rOStream) const override;
 
 
protected:
 
 
    // Default constructor necessary for serialization
    FluxCondition();
 
 
    void AddIntegrationPointRHSContribution(
        VectorType& rRightHandSideVector,
        const FluxConditionInternals::IntegrationData<TNodeNumber>& rData);
 
    void CalculateNormal(array_1d<double,3>& An );
 
 
private:
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
    FluxCondition& operator=(FluxCondition const& rOther);
 
    FluxCondition(FluxCondition const& rOther);
 
 
}; // Class FluxCondition
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_FLUX_CONDITION_H_INCLUDED  defined