InsertionBoundary.h

Go to the documentation of this file.

// This file is part of the MercuryDPM project (https://www.mercurydpm.org).
// Copyright (c), The MercuryDPM Developers Team. All rights reserved.
// License: BSD 3-Clause License; see the LICENSE file in the root directory.
 
#ifndef BOUNDARIES_INSERTIONBOUNDARY_H
#define BOUNDARIES_INSERTIONBOUNDARY_H
 
#include "BaseBoundary.h"
#include <Math/PSD.h>
 
class MD;
 
class RNG;
 
class PSD;
 
class InsertionBoundary : public BaseBoundary
{
public:
 
//    /*!
// * \brief Defines a custom particle size distribution; distribution_ will always be used, unless particleSizeDistributionVector_ is non-empty
// */
//    enum class Distribution
//    {
//        Uniform,
//        Normal_1_5
//        // TODO add LogNormal distribution
//        // LogNormal
//    };
 
 
    InsertionBoundary();
    
    InsertionBoundary(const InsertionBoundary& other);
    
    ~InsertionBoundary() override;
 
 
//    /*!
//     * \brief Sets the properties of the InsertionBoundary for mutliple different particle types
//     */
//    virtual void set(std::vector<BaseParticle*> particleToCopy, unsigned int maxFailed, Vec3D velMin, Vec3D velMax,
//                    double radMin, double radMax)=0;
//
//    /*!
//     * \brief Sets the properties of the InsertionBoundary for a single particle type
//     */
//    virtual void set(BaseParticle* particleToCopy, unsigned int maxFailed, Vec3D velMin, Vec3D velMax, double radMin,
//                  double radMax)=0;
    
    virtual BaseParticle* generateParticle(RNG& random);
 
    virtual void placeParticle(BaseParticle* p, RNG& random) = 0;
 
    virtual void shiftBoundary(Vec3D shift);
 
    virtual void rotateBoundary(Vec3D angle);
 
    void checkBoundaryBeforeTimeStep(DPMBase* md) override;
 
    void insertParticles(DPMBase* md);
 
    unsigned int getNumberOfParticlesInserted() const;
    
    Mdouble getMassOfParticlesInserted() const;
 
   Mdouble getVolumeOfParticlesInserted() const;
 
    void reset();
 
    void activate();
 
    void deactivate();
 
    bool isActivated();
 
 
    
    unsigned int getMaxFailed() const;
    
    void setParticleToCopy(std::vector<BaseParticle*> particleToCopy);
 
    void setParticleToCopy(BaseParticle* particleToCopy);
    
    std::vector<BaseParticle*> getParticleToCopy();
    
    void read(std::istream& is) override;
    
    void write(std::ostream& os) const override;
 
    Mdouble getVolumeFlowRate() const;
 
    void setVolumeFlowRate(Mdouble volumeFlowRate);
 
    Mdouble getInitialVolume() const;
 
    void setInitialVolume(Mdouble initialVolume);
 
    void setPSD(const PSD psd);
 
    void setPSD(std::vector<PSD> psd, std::vector<Mdouble> probability);
 
    std::vector<PSD> getPSD();
 
    void
    setVariableVolumeFlowRate(const std::vector<Mdouble>& variableCumulativeVolumeFlowRate, Mdouble samplingInterval);
 
    bool insertParticle(Mdouble time);
 
    bool getCheckParticleForInteraction() const;
 
    void setCheckParticleForInteraction(bool checkParticleForInteraction);
 
    void setManualInsertion(bool manualInsertion);
 
//    /*!
//     * \brief write Distribution class to file.
//     */
//    friend std::ostream& operator<<(std::ostream& os, InsertionBoundary::Distribution type);
//
//    /*!
//     * \brief read Distribution class from file.
//     */
//    friend std::istream& operator>>(std::istream& is, InsertionBoundary::Distribution& type);
 
protected:
 
    std::vector<BaseParticle*> particleToCopy_;
    
    unsigned int maxFailed_;
    
    unsigned int numberOfParticlesInserted_;
 
    Mdouble massInserted_;
 
    Mdouble volumeInserted_;
 
    bool isActivated_;
 
    Mdouble volumeFlowRate_;
 
    Mdouble initialVolume_;
 
    std::vector<Mdouble> variableCumulativeVolumeFlowRate_;
 
    Mdouble samplingInterval_;
 
    bool checkParticleForInteraction_;
 
    std::vector<PSD> particleSizeDistributionVector_;
 
    Vec3D velMin_, velMax_;
 
    bool isManuallyInserting_;
 
    std::vector<Mdouble> probability_;
    
    int chosenSpecies_;
 
    BaseParticle* p0= nullptr;
};
 
#endif