HorizontalMixerWalls.cpp File Reference

#include "HorizontalMixer.h"

Classes
class	HorizontalMixerWalls

Functions
int	main ()

Function Documentation

main()

int main

(

)

{
    Mdouble revolutionsPerSecond = 0.25;
    Mdouble rotationSpeed = 0.25*constants::pi*revolutionsPerSecond;
    Mdouble particleRadius = 0.2;
    Mdouble timeMin = 0.0;
    Mdouble fillHeight = 0;
    HorizontalMixerWalls mixer(particleRadius, rotationSpeed, timeMin, fillHeight);
 
    //set name
    mixer.setName("HorizontalMixerWalls");
 
    //remove old files (note this is a bit dangerous)
    //mixer.removeOldFiles();
 
    //set species
    auto s = mixer.speciesHandler.copyAndAddObject(LinearViscoelasticSlidingFrictionSpecies());
    s->setDensity(2000);
    Mdouble mass = s->getMassFromRadius(particleRadius);
    s->setCollisionTimeAndRestitutionCoefficient(0.01, 0.5, mass);
    s->setSlidingFrictionCoefficient(0.5);
    s->setSlidingStiffness(2.0 / 7.0 * s->getStiffness());
    s->setSlidingDissipation(2.0 / 7.0 * s->getDissipation());
    
    //set timestep
    mixer.setTimeStep(0.2 * s->getCollisionTime(mass));
    //save every 2 collision times for smooth viewable output
    mixer.setSaveCount((unsigned) (20.0 * s->getCollisionTime(mass) / mixer.getTimeStep()));
    logger(INFO, "Savecount: %", mixer.dataFile.getSaveCount());
    
    mixer.setTimeMax(100);
    
    mixer.solve();
    return 0;
}

References BaseHandler< T >::copyAndAddObject(), DPMBase::dataFile, File::getSaveCount(), DPMBase::getTimeStep(), INFO, logger, constants::pi, DPMBase::setName(), DPMBase::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::solve(), and DPMBase::speciesHandler.