Funnel.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 FUNNEL_H
#define FUNNEL_H
#include "Chute.h"
 
 
class Funnel : public Chute {
    public:
    
    Funnel() {constructor();}
    
    //Funnel(MD& other) : MD(other), Mercury3D(other), Chute() {constructor();}
    //Funnel(HGRID_base& other) :  MD(other), Mercury3D(other), Chute() {constructor();}
    //Funnel(Mercury3D& other) : MD(other), Mercury3D(other), Chute() {constructor();}
 
    void write(std::ostream& os, bool writeAllParticles = false) const override;
    
    void read(std::istream& is, ReadOptions opt = ReadOptions::ReadAll) override;
    
    void constructor();
    
    void cleanChute();
    
    void set_funnz(double funnz_){funnz=funnz_;}
    
    double get_funnz() const {return funnz;}
    
    void set_funfr(double funfr_){
        if (funfr_ < 0.0 || funfr_ > 1.0)
        {
            logger(WARN, "Filling Ratio is below 0 or above 1, default value of 0.33 is used!");
        }
        else
        {
            funfr = funfr_;
        }
    }
    
    double get_funfr() const {return funfr;}
    
    double get_funr(){return funr;}
    
    void set_funO(double x, double y){funO[0]=x; funO[1]=y;}
    void set_funO(double* x){funO[0]=x[0]; funO[1]=x[1];}
    
    const double* get_funO() const {return funO;}
    
    double get_funOx() const {return funO[0];}
    
    double get_funOy() const {return funO[1];}
    
    double get_fundiag() const {return fundiag;}
    
    double get_funrmax() const {return funrmax;}
    
    void set_funa(double funa_){funa=funa_*constants::pi/180;}
    
    double get_funa() const {return funa;}
        
    double get_funH() const {return funH;}
    
    void set_funHf(double funHf_){funHf=funHf_;}
    
    double get_funHf() const {return funHf;}
    
    void set_funD(double funD_){funD=funD_;}
    
    double get_funD() const {return funD;}
    
    void setName_();
 
protected:
    SphericalParticle inflowParticle_;
    
    void set_funH(double funH_){funH=funH_;}
 
    void set_funrmax(double funrmax_){funrmax=funrmax_;}
 
    void set_fundiag(double fundiag_){fundiag=fundiag_;}
 
    void set_funr(double funr_){funr=funr_;}
 
    virtual void setupInitialConditions() override;
    
    virtual void create_inflow_particle();
    
    virtual void create_funnel();
    
    virtual void update_funnel();
    
    virtual void check_funnel();
    
    virtual void create_walls();
    
    virtual bool readNextArgument(int& i, int argc, char *argv[]) override;
 
    double funr; // Funnel radius.
    double funO[2]; // Origin of the funnel is the (x,y) location of the center of the top of the funnel.
    double funa; // Angle of the funnel
    double funH; // Heigth of the funnel
    double funHf; // Falling heigth
    double funD; // Funnel diameter at the downside of the funnel.
    double funnz; //Number of particles along the heigth of the funnel
    double funfr; // Filling ratio of the funnel
    double fundiag; //The diagonal of the filling region
    double funrmax; //The maximum range for r
    
};
#endif