general_bounding.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 and HU1527/12-4.                                          *
*                                                                           *
* Portions copyright (c) 2017-2022 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 GENERAL_BOUNDING_H
#define GENERAL_BOUNDING_H

#include "general_dynamics.h"

namespace SPH
{
    class BoundingInAxisDirection : public ParticleDynamics<void>, public GeneralDataDelegateSimple
    {
    protected:
        const int axis_;                 
        BoundingBox body_domain_bounds_; 
        StdLargeVec<Vecd> &pos_;
        BaseCellLinkedList *cell_linked_list_;
        Real cut_off_radius_max_; 
    public:
        BoundingInAxisDirection(RealBody &real_body, int axis_direction);
        virtual ~BoundingInAxisDirection(){};
    };

    class PeriodicConditionInAxisDirection
    {
    protected:
        Vecd periodic_translation_;
        StdVec<CellLists> bound_cells_;
        void setPeriodicTranslation(BoundingBox &body_domain_bounds, int axis_direction);

        class PeriodicBounding : public BoundingInAxisDirection
        {
        protected:
            Vecd &periodic_translation_;
            StdVec<CellLists> &bound_cells_;

            virtual void checkLowerBound(size_t index_i, Real dt = 0.0);
            virtual void checkUpperBound(size_t index_i, Real dt = 0.0);

        public:
            PeriodicBounding(Vecd &periodic_translation,
                             StdVec<CellLists> &bound_cells, RealBody &real_body, int axis_direction)
                : BoundingInAxisDirection(real_body, axis_direction), periodic_translation_(periodic_translation),
                  bound_cells_(bound_cells){};
            virtual ~PeriodicBounding(){};

            virtual void exec(Real dt = 0.0) override;
            virtual void parallel_exec(Real dt = 0.0) override;
        };

        class PeriodicCondition : public BoundingInAxisDirection
        {
        protected:
            Vecd &periodic_translation_;
            StdVec<CellLists> &bound_cells_;

            virtual void checkLowerBound(ListData &list_data, Real dt = 0.0) = 0;
            virtual void checkUpperBound(ListData &list_data, Real dt = 0.0) = 0;

        public:
            PeriodicCondition(Vecd &periodic_translation,
                              StdVec<CellLists> &bound_cells, RealBody &real_body, int axis_direction)
                : BoundingInAxisDirection(real_body, axis_direction), periodic_translation_(periodic_translation),
                  bound_cells_(bound_cells){};
            virtual ~PeriodicCondition(){};

            virtual void exec(Real dt = 0.0) override;
            virtual void parallel_exec(Real dt = 0.0) override { exec(); };
        };

    public:
        PeriodicConditionInAxisDirection(RealBody &real_body, int axis_direction);
        virtual ~PeriodicConditionInAxisDirection(){};
    };

    class PeriodicConditionInAxisDirectionUsingCellLinkedList : public PeriodicConditionInAxisDirection
    {
    protected:
        class PeriodicCellLinkedList : public PeriodicCondition
        {
        protected:
            virtual void checkLowerBound(ListData &list_data, Real dt = 0.0) override;
            virtual void checkUpperBound(ListData &list_data, Real dt = 0.0) override;

        public:
            PeriodicCellLinkedList(Vecd &periodic_translation,
                                   StdVec<CellLists> &bound_cells, RealBody &real_body, int axis_direction)
                : PeriodicCondition(periodic_translation, bound_cells, real_body, axis_direction){};
            virtual ~PeriodicCellLinkedList(){};
        };

    public:
        PeriodicConditionInAxisDirectionUsingCellLinkedList(RealBody &real_body, int axis_direction)
            : PeriodicConditionInAxisDirection(real_body, axis_direction),
              bounding_(this->periodic_translation_, this->bound_cells_, real_body, axis_direction),
              update_cell_linked_list_(this->periodic_translation_, this->bound_cells_, real_body, axis_direction){};
        virtual ~PeriodicConditionInAxisDirectionUsingCellLinkedList(){};

        PeriodicBounding bounding_;
        PeriodicCellLinkedList update_cell_linked_list_;
    };

    class OpenBoundaryConditionInAxisDirection
    {
    protected:
        StdVec<CellLists> bound_cells_;

        class ParticleTypeTransfer : public BoundingInAxisDirection
        {
        protected:
            StdVec<CellLists> &bound_cells_;
            ParticleFunctor checking_bound_;
            virtual void checkLowerBound(size_t index_i, Real dt = 0.0);
            virtual void checkUpperBound(size_t index_i, Real dt = 0.0);

        public:
            ParticleTypeTransfer(StdVec<CellLists> &bound_cells, RealBody &real_body, int axis_direction, bool positive)
                : BoundingInAxisDirection(real_body, axis_direction),
                  bound_cells_(bound_cells)
            {
                checking_bound_ = positive ? std::bind(&OpenBoundaryConditionInAxisDirection::ParticleTypeTransfer::checkUpperBound, this, _1, _2)
                                           : std::bind(&OpenBoundaryConditionInAxisDirection::ParticleTypeTransfer::checkLowerBound, this, _1, _2);
            };
            virtual ~ParticleTypeTransfer(){};

            virtual void exec(Real dt = 0.0) override;
            virtual void parallel_exec(Real dt = 0.0) override { exec(); };
        };

    public:
        OpenBoundaryConditionInAxisDirection(RealBody &real_body, int axis_direction, bool positive);
        virtual ~OpenBoundaryConditionInAxisDirection(){};

        ParticleTypeTransfer particle_type_transfer;
    };

    class PeriodicConditionInAxisDirectionUsingGhostParticles : public PeriodicConditionInAxisDirection
    {
    protected:
        StdVec<IndexVector> ghost_particles_;

        class CreatPeriodicGhostParticles : public PeriodicBounding
        {
        protected:
            StdVec<IndexVector> &ghost_particles_;
            virtual void setupDynamics(Real dt = 0.0) override;
            virtual void checkLowerBound(size_t index_i, Real dt = 0.0) override;
            virtual void checkUpperBound(size_t index_i, Real dt = 0.0) override;

        public:
            CreatPeriodicGhostParticles(Vecd &periodic_translation, StdVec<CellLists> &bound_cells,
                                        StdVec<IndexVector> &ghost_particles, RealBody &real_body, int axis_direction)
                : PeriodicBounding(periodic_translation, bound_cells, real_body, axis_direction),
                  ghost_particles_(ghost_particles){};
            virtual ~CreatPeriodicGhostParticles(){};

            virtual void parallel_exec(Real dt = 0.0) override { exec(); };
        };

        class UpdatePeriodicGhostParticles : public PeriodicBounding
        {
        protected:
            StdVec<IndexVector> &ghost_particles_;
            void checkLowerBound(size_t index_i, Real dt = 0.0) override;
            void checkUpperBound(size_t index_i, Real dt = 0.0) override;

        public:
            UpdatePeriodicGhostParticles(Vecd &periodic_translation, StdVec<CellLists> &bound_cells,
                                         StdVec<IndexVector> &ghost_particles, RealBody &real_body, int axis_direction)
                : PeriodicBounding(periodic_translation, bound_cells, real_body, axis_direction),
                  ghost_particles_(ghost_particles){};
            virtual ~UpdatePeriodicGhostParticles(){};

            virtual void exec(Real dt = 0.0) override;
            virtual void parallel_exec(Real dt = 0.0) override;
        };

    public:
        PeriodicConditionInAxisDirectionUsingGhostParticles(RealBody &real_body, int axis_direction)
            : PeriodicConditionInAxisDirection(real_body, axis_direction),
              bounding_(this->periodic_translation_, this->bound_cells_, real_body, axis_direction),
              ghost_creation_(this->periodic_translation_, this->bound_cells_, this->ghost_particles_, real_body, axis_direction),
              ghost_update_(this->periodic_translation_, this->bound_cells_, this->ghost_particles_, real_body, axis_direction)
        {
            ghost_particles_.resize(2);
        };

        virtual ~PeriodicConditionInAxisDirectionUsingGhostParticles(){};

        PeriodicBounding bounding_;
        CreatPeriodicGhostParticles ghost_creation_;
        UpdatePeriodicGhostParticles ghost_update_;
    };

    class MirrorBoundaryConditionInAxisDirection : public BoundingInAxisDirection
    {
    protected:
        CellLists bound_cells_;
        IndexVector ghost_particles_;

        class MirrorBounding : public BoundingInAxisDirection
        {
        protected:
            CellLists &bound_cells_;
            virtual void checkLowerBound(size_t index_i, Real dt = 0.0);
            virtual void checkUpperBound(size_t index_i, Real dt = 0.0);
            ParticleFunctor checking_bound_;

            StdLargeVec<Vecd> &vel_;
            void mirrorInAxisDirection(size_t particle_index_i, Vecd body_bound, int axis_direction);

        public:
            MirrorBounding(CellLists &bound_cells, RealBody &real_body, int axis_direction, bool positive);
            virtual ~MirrorBounding(){};
            virtual void exec(Real dt = 0.0) override;
            virtual void parallel_exec(Real dt = 0.0) override;
        };

        class CreatingGhostParticles : public MirrorBounding
        {
        protected:
            IndexVector &ghost_particles_;
            virtual void setupDynamics(Real dt = 0.0) override { ghost_particles_.clear(); };
            virtual void checkLowerBound(size_t index_i, Real dt = 0.0) override;
            virtual void checkUpperBound(size_t index_i, Real dt = 0.0) override;

        public:
            CreatingGhostParticles(IndexVector &ghost_particles, CellLists &bound_cells,
                                   RealBody &real_body, int axis_direction, bool positive);
            virtual ~CreatingGhostParticles(){};
            virtual void parallel_exec(Real dt = 0.0) override { exec(); };
        };

        class UpdatingGhostStates : public MirrorBounding
        {
        protected:
            IndexVector &ghost_particles_;
            void checkLowerBound(size_t index_i, Real dt = 0.0) override;
            void checkUpperBound(size_t index_i, Real dt = 0.0) override;
            ParticleFunctor checking_bound_update_;

        public:
            UpdatingGhostStates(IndexVector &ghost_particles, CellLists &bound_cells,
                                RealBody &real_body, int axis_direction, bool positive);
            virtual ~UpdatingGhostStates(){};

            virtual void exec(Real dt = 0.0) override;
            virtual void parallel_exec(Real dt = 0.0) override;
        };

    public:
        MirrorBounding bounding_;
        CreatingGhostParticles creating_ghost_particles_;
        UpdatingGhostStates updating_ghost_states_;

        MirrorBoundaryConditionInAxisDirection(RealBody &real_body, int axis_direction, bool positive);
        virtual ~MirrorBoundaryConditionInAxisDirection(){};

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