ChuteWithVerticalHopper.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.
 
#include "Chute.h"
#include "Statistics.h"
 
class ChuteWithVerticalHopper : public Chute {
public:
    double inflowLength;
    void set_inflowLength(double new_){inflowLength=new_;}
 
    double dragCoefficient;
    void set_dragCoefficient(double new_){dragCoefficient=new_;}
 
    Vec3D targetVelocity;
    void set_targetVelocity(Vec3D new_){targetVelocity=new_;}
 
private:
    void setupInitialConditions() override {}
    
    void computeExternalForces(int CI) {
        MD::computeExternalForces(CI);
 
        if (Particles[CI].Position.X<inflowLength) {
            Particles[CI].Force += dragCoefficient*(targetVelocity-Particles[CI].Velocity);
        }
    }
    
    void create_inflow_particle()
    {
        //the following formula yields polydispersed particle radii:
        P0.Radius = random(MinInflowParticleRadius,MaxInflowParticleRadius);
        P0.computeMass(Species);
        P0.Velocity.set_zero();
        P0.Position = Vec3D(
            random(getXMin(),inflowLength/2.0),
            random(getYMin(),getYMax()),
            random(getZMin()+1.5*P0.Radius,getZMax()));
    }
    
};