particle_dynamics_dissipation.h

Go to the documentation of this file.

/* -------------------------------------------------------------------------*
 *                              SPHinXsys                                   *
 * --------------------------------------------------------------------------*
 * SPHinXsys (pronunciation: s'finksis) is an acronym from Smoothed Particle    *
 * Hydrodynamics for industrial compleX systems. It provides C++ APIs for   *
 * physical accurate simulation and aims to model coupled industrial dynamic *
 * systems including fluid, solid, multi-body dynamics and beyond with SPH  *
 * (smoothed particle hydrodynamics), a meshless computational method using *
 * particle discretization.                                                 *
 *                                                                          *
 * SPHinXsys is partially funded by German Research Foundation              *
 * (Deutsche Forschungsgemeinschaft) DFG HU1527/6-1, HU1527/10-1                *
 * and HU1527/12-1.                                                         *
 *                                                                           *
 * Portions copyright (c) 2017-2020 Technical University of Munich and      *
 * the authors' affiliations.                                               *
 *                                                                           *
 * Licensed under the Apache License, Version 2.0 (the "License"); you may   *
 * not use this file except in compliance with the License. You may obtain a *
 * copy of the License at http://www.apache.org/licenses/LICENSE-2.0.        *
 *                                                                           *
 * --------------------------------------------------------------------------*/
#ifndef PARTICLE_DYNAMICS_DISSIPATION_H
#define PARTICLE_DYNAMICS_DISSIPATION_H

#include "all_particle_dynamics.h"

namespace SPH
{
    typedef DataDelegateInner<SPHBody, BaseParticles, BaseMaterial> DissipationDataInner;
    typedef DataDelegateContact<SPHBody, BaseParticles, BaseMaterial,
                                SPHBody, BaseParticles, BaseMaterial, DataDelegateEmptyBase>
        DissipationDataContact;
    typedef DataDelegateContact<SPHBody, BaseParticles, BaseMaterial,
                                SolidBody, SolidParticles, Solid, DataDelegateEmptyBase>
        DissipationDataWithWall;

    template <typename VariableType>
    struct ErrorAndParameters
    {
        VariableType error_;
        Real a_, c_;
        ErrorAndParameters(Real zero = 0.0) : error_(zero), a_(zero), c_(zero){};
    };

    template <typename VariableType>
    class DampingBySplittingInner : public InteractionDynamicsSplitting, public DissipationDataInner
    {
    protected:
    public:
        DampingBySplittingInner(BaseBodyRelationInner &inner_relation, const std::string &variable_name, Real eta);
        virtual ~DampingBySplittingInner(){};
        void resetDampingCoefficient(Real reset_ratio) { eta_ *= reset_ratio; };

    protected:
        Real eta_; 
        StdLargeVec<Real> &Vol_, &mass_;
        StdLargeVec<VariableType> &variable_;

        virtual ErrorAndParameters<VariableType> computeErrorAndParameters(size_t index_i, Real dt = 0.0);
        virtual void updateStates(size_t index_i, Real dt, const ErrorAndParameters<VariableType> &error_and_parameters);
        virtual void Interaction(size_t index_i, Real dt = 0.0) override;
    };

    template <typename VariableType>
    class DampingBySplittingComplex : public DampingBySplittingInner<VariableType>, public DissipationDataContact
    {
    public:
        DampingBySplittingComplex(ComplexBodyRelation &complex_relation, const std::string &variable_name, Real eta);
        virtual ~DampingBySplittingComplex(){};

    protected:
        virtual ErrorAndParameters<VariableType> computeErrorAndParameters(size_t index_i, Real dt = 0.0) override;
        virtual void updateStates(size_t index_i, Real dt, const ErrorAndParameters<VariableType> &error_and_parameters) override;

    private:
        StdVec<StdLargeVec<Real> *> contact_Vol_, contact_mass_;
        StdVec<StdLargeVec<VariableType> *> contact_variable_;
    };

    template <typename VariableType,
              template <typename BaseVariableType>
              class BaseDampingBySplittingType>
    class DampingBySplittingWithWall : public BaseDampingBySplittingType<VariableType>, public DissipationDataWithWall
    {
    public:
        DampingBySplittingWithWall(ComplexBodyRelation &complex_wall_relation, const std::string &variable_name, Real eta);
        virtual ~DampingBySplittingWithWall(){};

    protected:
        virtual ErrorAndParameters<VariableType> computeErrorAndParameters(size_t index_i, Real dt = 0.0) override;

    private:
        StdVec<StdLargeVec<Real> *> wall_Vol_;
        StdVec<StdLargeVec<VariableType> *> wall_variable_;
    };

    template <typename VariableType>
    class DampingPairwiseInner : public InteractionDynamicsSplitting, public DissipationDataInner
    {
    public:
        DampingPairwiseInner(BaseBodyRelationInner &inner_relation, const std::string &variable_name, Real eta);
        virtual ~DampingPairwiseInner(){};
        void resetDampingCoefficient(Real reset_ratio) { eta_ *= reset_ratio; };

    protected:
        StdLargeVec<Real> &Vol_, &mass_;
        StdLargeVec<VariableType> &variable_;
        Real eta_; 
        virtual void Interaction(size_t index_i, Real dt = 0.0) override;
    };

    template <typename VariableType>
    class DampingPairwiseComplex : public DampingPairwiseInner<VariableType>, public DissipationDataContact
    {
    public:
        DampingPairwiseComplex(BaseBodyRelationInner &inner_relation,
                               BaseBodyRelationContact &contact_relation, const std::string &variable_name, Real eta);
        DampingPairwiseComplex(ComplexBodyRelation &complex_relation, const std::string &variable_name, Real eta);
        virtual ~DampingPairwiseComplex(){};

    protected:
        virtual void Interaction(size_t index_i, Real dt = 0.0) override;

    private:
        StdVec<StdLargeVec<Real> *> contact_Vol_, contact_mass_;
        StdVec<StdLargeVec<VariableType> *> contact_variable_;
    };

    template <typename VariableType,
              template <typename BaseVariableType> class BaseDampingPairwiseType>
    class DampingPairwiseWithWall : public BaseDampingPairwiseType<VariableType>,
                                    public DissipationDataWithWall
    {
    public:
        DampingPairwiseWithWall(BaseBodyRelationInner &inner_relation,
                                BaseBodyRelationContact &contact_relation, const std::string &variable_name, Real eta);
        DampingPairwiseWithWall(ComplexBodyRelation &complex_wall_relation, const std::string &variable_name, Real eta);
        virtual ~DampingPairwiseWithWall(){};

    protected:
        virtual void Interaction(size_t index_i, Real dt = 0.0) override;

    private:
        StdVec<StdLargeVec<Real> *> wall_Vol_;
        StdVec<StdLargeVec<VariableType> *> wall_variable_;
    };

    template <typename VariableType>
    class DampingPairwiseFromWall : public InteractionDynamicsSplitting,
                                  public DataDelegateContact<SPHBody, BaseParticles, BaseMaterial,
                                                             SolidBody, SolidParticles, Solid>
    {
    public:
        DampingPairwiseFromWall(BaseBodyRelationContact &contact_relation, const std::string &variable_name, Real eta);
        virtual ~DampingPairwiseFromWall(){};

    protected:
        virtual void Interaction(size_t index_i, Real dt = 0.0) override;

    private:
        Real eta_; 
        StdLargeVec<Real> &Vol_, &mass_;
        StdLargeVec<VariableType> &variable_;
        StdVec<StdLargeVec<Real> *> wall_Vol_;
        StdVec<StdLargeVec<VariableType> *> wall_variable_;
    };

    template <class DampingAlgorithmType>
    class DampingWithRandomChoice : public DampingAlgorithmType
    {
    protected:
        Real random_ratio_;
        bool RandomChoice();

    public:
        template <typename... ConstructorArgs>
        DampingWithRandomChoice(Real random_ratio, ConstructorArgs &&...args);
        virtual ~DampingWithRandomChoice(){};

        virtual void exec(Real dt = 0.0) override;
        virtual void parallel_exec(Real dt = 0.0) override;
    };
}
#endif // PARTICLE_DYNAMICS_DISSIPATION_H