MonodispersedDrum.cpp File Reference

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cmath>
#include "Mercury3D.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"
#include "Walls/AxisymmetricIntersectionOfWalls.h"
#include "Walls/InfiniteWall.h"
#include "Boundaries/PeriodicBoundary.h"
#include <chrono>

Classes
class	DrumRot

Functions
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

{
 
    // Start measuring elapsed time
    //std::chrono::time_point<std::chrono::system_clock> startClock, endClock;
    //startClock = std::chrono::system_clock::now();
 
    
    DrumRot problem;
    problem.setName("MonodisperesedDrum");
 
    //setting locally-available variables to define key
    //parameters such that they can be automatically included in
    //the name of the file
    //*******************Excitation Properties*****************************************
    //Drum rotation rate (rpm)
    double rotRateBasal = 1.0;
 
    //*******************Frictional Properties*****************************************
    //sliding friction for species 1, 2 and wall
    double muSWall = 0.61;
    double muS1 = 0.19;
 
    //rolling friction for species 1, 2 and wall
    double muRWall = 0.06;
    double muR1 = 0.01;
 
    //torsion friction for species 1, 2 and wall
    double muTWall = 0.0;
    double muT1 = 0.000;
 
    double dimDrumRad = 15.0;
    double drumRad = dimDrumRad*0.0005; // 45 * smallest particle radius -> dimDrumRad particles fit in drum diameter
    //the dimensionless drum length (L/d_l)
    double dimDrumLength = 4.0;
    double drumLength = dimDrumLength*0.001; // dimDrumLength * smallest particle diameter
 
    // computes rpm based on froude number and drumRad
    problem.setTimeMax(3.0);
    problem.setTimeStep(1.0/(4000.0 * 50.0));
    double froudeNumber = 0.22;
    double rotRate = pow(froudeNumber*9.81/drumRad,0.5);//*(60.0/(2.0*3.1415926535));
    logger(INFO, "rotRate = %", rotRate);
    //Set the number of domains for parallel decomposition
    //problem.setNumberOfDomains({2,1,1});
 
    //*******************Setting Up Filename******************************************
    double dRatio;
    double rhoRatio;
 
    
 
    problem.setDrumRadius(drumRad);// in meters
 
    problem.setXMin(-problem.getDrumRadius());
    problem.setYMin(-0.5*drumLength);
    problem.setZMin(-problem.getDrumRadius());
 
    problem.setXMax(problem.getDrumRadius());
    problem.setYMax(0.5*drumLength);// in meters
    problem.setZMax(problem.getDrumRadius());
 
    problem.setGravity(Vec3D(0.,0.,-9.81));
    problem.setCOR(0.831,0.831);//drumWall, species1, species2
    problem.setFractionalPolydispersity(0.0);//10% dispersity
    problem.setDrumFillFraction(0.2);// At 0.5 the drum is 3/4 filled.
    //redundant
    //problem.setFrictionCoeff(.6,.19);//(particle-wall-friction,part-part-friction)
    problem.setSlidingFriction(muSWall,muS1); //wall, species1, species2
    problem.setRollingFriction(muRWall,muR1); //wall, species1, species2
    problem.setTorsionFriction(muTWall,muT1); //wall, species1, species2
 
    problem.setRevolutionSpeed(rotRate);//rpm
 
    problem.setSaveCount(2000);
    problem.readArguments(argc,argv);
 
    problem.dataFile.setFileType(FileType::ONE_FILE);
    problem.restartFile.setFileType(FileType::ONE_FILE);
    problem.fStatFile.setFileType(FileType::ONE_FILE);
    problem.eneFile.setFileType(FileType::ONE_FILE);
 
    problem.solve();
 
    return 0;
}

References INFO, logger, ONE_FILE, Eigen::ArrayBase< Derived >::pow(), and problem.