wave_condition.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Miguel Maso Sotomayor
//
 
#ifndef KRATOS_WAVE_CONDITION_H_INCLUDED
#define KRATOS_WAVE_CONDITION_H_INCLUDED
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "includes/condition.h"
#include "includes/serializer.h"
 
namespace Kratos
{
 
 
 
 
template<std::size_t TNumNodes>
class WaveCondition : public Condition
{
public:
 
    typedef std::size_t IndexType;
 
    typedef Node NodeType;
 
    typedef array_1d<double, 3*TNumNodes> LocalVectorType;
 
    typedef BoundedMatrix<double, 3*TNumNodes, 3*TNumNodes> LocalMatrixType;
 
    typedef GeometryType::ShapeFunctionsGradientsType ShapeFunctionsGradientsType;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(WaveCondition);
 
 
 
    WaveCondition() : Condition(){}
 
    WaveCondition(IndexType NewId, const NodesArrayType& ThisNodes) : Condition(NewId, ThisNodes){}
 
    WaveCondition(IndexType NewId, GeometryType::Pointer pGeometry) : Condition(NewId, pGeometry){}
 
    WaveCondition(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties) : Condition(NewId, pGeometry, pProperties){}
 
    ~ WaveCondition() override {};
 
 
    Condition::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<WaveCondition<TNumNodes>>(NewId, this->GetGeometry().Create(ThisNodes), pProperties);
    }
 
    Condition::Pointer Create(IndexType NewId, GeometryType::Pointer pGeom, PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<WaveCondition<TNumNodes>>(NewId, pGeom, pProperties);
    }
 
    Condition::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override
    {
        Condition::Pointer p_new_elem = Create(NewId, this->GetGeometry().Create(ThisNodes), this->pGetProperties());
        p_new_elem->SetData(this->GetData());
        p_new_elem->Set(Flags(*this));
        return p_new_elem;
    }
 
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
    void EquationIdVector(EquationIdVectorType& rResult, const ProcessInfo& rCurrentProcessInfo) const override;
 
    void GetDofList(DofsVectorType& rConditionalDofList, const ProcessInfo& rCurrentProcessInfo) const override;
 
    void GetValuesVector(Vector& rValues, int Step = 0) const override;
 
    void GetFirstDerivativesVector(Vector& rValues, int Step = 0) const override;
 
    void GetSecondDerivativesVector(Vector& rValues, int Step = 0) const override;
 
    void Calculate(
        const Variable<array_1d<double,3>>& rVariable,
        array_1d<double,3>& rOutput,
        const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLocalSystem(MatrixType& rLeftHandSideMatrix, VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateMassMatrix(MatrixType& rMassMatrix, const ProcessInfo& rCurrentProcessInfo) override;
 
 
    const Parameters GetSpecifications() const override;
 
 
    std::string Info() const override
    {
        return "WaveCondition";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << Info() << " : " << Id();
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        rOStream << GetGeometry();
    }
 
 
 
 
protected:
 
    static constexpr IndexType mLocalSize = 3 * TNumNodes;
 
 
    struct ConditionData
    {
        bool integrate_by_parts;
        double stab_factor;
        double relative_dry_height;
        double length;
        double gravity;
 
        double depth;
        double height;
        array_1d<double,3> velocity;
 
        double v_neumann;
        double h_dirichlet;
        array_1d<double,3> flux;
 
        array_1d<double,TNumNodes> nodal_f;
        array_1d<double,TNumNodes> nodal_h;
        array_1d<double,TNumNodes> nodal_z;
        array_1d<array_1d<double,3>,TNumNodes> nodal_v;
        array_1d<array_1d<double,3>,TNumNodes> nodal_q;
 
        array_1d<double,3> normal;
    };
 
 
    virtual const Variable<double>& GetUnknownComponent(int Index) const;
 
    virtual LocalVectorType GetUnknownVector(ConditionData& rData);
 
    virtual void InitializeData(
        ConditionData& rData,
        const ProcessInfo& rProcessInfo);
 
    virtual void CalculateGaussPointData(
        ConditionData& rData,
        const IndexType PointIndex,
        const array_1d<double,TNumNodes>& rN);
 
    static void CalculateGeometryData(
        const GeometryType& rGeometry,
        Vector &rGaussWeights,
        Matrix &rNContainer);
 
    static const array_1d<double,3> VectorProduct(
        const array_1d<array_1d<double,3>,TNumNodes>& rV,
        const array_1d<double,TNumNodes>& rN);
 
 
    void AddFluxTerms(
        LocalVectorType& rVector,
        const ConditionData& rData,
        const array_1d<double,TNumNodes>& rN,
        const double Weight = 1.0);
 
    void AddMassTerms(
        LocalMatrixType& rMatrix,
        const ConditionData& rData,
        const array_1d<double,TNumNodes>& rN,
        const double Weight);
 
 
private:
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Condition);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Condition);
    }
 
 
 
 
}; // Class WaveCondition
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_WAVE_CONDITION_H_INCLUDED  defined