SquareDrumTest.cpp File Reference

#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;
}

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