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

Classes
class	RotatingDrum

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

Function Documentation

◆ main()

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

 {
     RotatingDrum problem;
     /* Make sure to add an (unique) name here */
     problem.setName("SquareDrumTest");
  
     //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 = 10.0;
  
     //vibration amplitude and frequency
     double fVib = 0.0;
     double aVib = 0.0;
  
     //*******************Frictional Properties*****************************************
     //sliding friction for species 1, 2 and wall
     double muSWall = 0.6;
     double muS1 = 0.16;
     double muS2 = 0.16;
  
     //rolling friction for species 1, 2 and wall
     double muRWall = 0.06;
     double muR1 = 0.016;
     double muR2 = 0.016;
  
     //torsion friction for species 1, 2 and wall
     double muTWall = 0.0;
     double muT1 = 0.000;
     double muT2 = 0.000;
  
     //the density ratio of particles
     double rhoRatio = 5000.0/2500.0;
     //the size ratio of particles
     double dRatio = 1;
  
     // computes rpm based on froude number and drumRad
     //the fraction of species 1 particles
     double specFrac = 0.5;
  
     double drumRad = 0.150/2; //diameter is 0.0150 m so drumRad is halve of the width
     double drumLength = 0.024; // depth of the drum is 0.024 m
  
     problem.setTimeMax(5.0);
     problem.setTimeStep(1.0/(2000.0 * 50.0));
     double froudeNumber = 0.22; //was 0.22
     double rotRate = 20; //pow(froudeNumber*9.81/drumRad,0.5);
     logger(INFO, "rotRate = %", rotRate);
  
     //*******************Setting Up Filename******************************************
  
     problem.setDrumRadius(drumRad);// in meters
  
     problem.setXMin(-(pow(2*pow(drumRad,2.0),0.5)));
     problem.setYMin(-drumLength/2);
     problem.setZMin(-(pow(2*pow(drumRad,2.0),0.5)));
  
     problem.setXMax((pow(2*pow(drumRad,2.0),0.5)));
     problem.setYMax(drumLength/2);// in meters
     problem.setZMax((pow(2*pow(drumRad,2.0),0.5)));
  
     problem.setGravity(Vec3D(0.,0.,-9.81));
     problem.setCOR(0.97,0.97,0.97);//drumWall, species1, species2 was all 0.97
     problem.setSizeAndDensityRatio(dRatio,rhoRatio);//size ratio and density ratio
     problem.setFractionalPolydispersity(0.00);
     problem.setDrumFillFraction(0.4);// At 0.5 the drum is 3/4 filled.
     problem.setSpeciesVolumeFraction(specFrac);//Species1 volume fraction
  
     problem.setSlidingFriction(muSWall,muS1,muS2); //wall, species1, species2
     problem.setRollingFriction(muRWall,muR1,muR2); //wall, species1, species2
     problem.setTorsionFriction(muTWall,muT1,muT2); //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.wallHandler.setWriteVTK(true);
     problem.setParticlesWriteVTK(true);
     problem.solve();
     return 0;
 }

Classes

Functions

Function Documentation

◆ main()