Lawine.cpp File Reference

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <sys/types.h>
#include <sys/stat.h>
#include "scr/Chute.h"

Classes
class	ChutePeriodic

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

Function Documentation

main()

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

{
    ChutePeriodic md;
    
    // Problem parameters
    md.setName("Lawine");
    md.fStatFile.setFileType(FileType::ONE_FILE);
    md.dataFile.setFileType(FileType::ONE_FILE);
    md.autoNumber();
    md.setTimeMax(100);
 
    // Particle properties
    md.setDensity(2400.0);
    md.setInflowParticleRadius(2.5e-3);//d=5mm
    md.speciesHandler.getObject(0)->setCollisionTimeAndRestitutionCoefficient(5e-4, 0.8);
    md.setTimeStep(1e-5); 
    md.setTimeMax(25); 
    md.speciesHandler.getObject(0)->setSlidingDissipation(md.get_dissipation());
    md.speciesHandler.getObject(0)->setSlidingStiffness(2./7.*md.speciesHandler.getObject(0)->getStiffness());
    md.speciesHandler.getObject(0)->setSlidingFrictionCoefficient(0.5);
    cout << "friction angle " << atan(md.getSlidingFrictionCoefficient())*180./pi << endl;
    md.setFixedParticleRadius(0);
    md.setFixedParticleRadius(md.getInflowParticleRadius());
    md.setRoughBottomType(MONOLAYER_ORDERED);
    
    // Chute properties
    md.setChuteAngle(0.0);
    md.setChuteWidth(450e-3);
    md.setZMax(150e-3);//initial height 150-200mm
    md.setChuteLength(1);
    md.setMaxFailed(1000);
    
    if (argc>1&&argv[1][0]!='-') {
        double ratio = pow(atof(argv[1]),1./3.);
        md.set_hopperLength_(md.getZMax()/mathsFunc::square(ratio));
        md.setChuteWidth(md.getChuteWidth()*ratio);
        md.setZMax(md.getZMax()*ratio);
        argc--; argv++;
    } else md.set_hopperLength_(md.getZMax());
        
    
    cout << "hopperLength_=" << md.get_hopperLength_() << endl; 
    cout << "ChuteWidth=" << md.getChuteWidth() << endl; 
    cout << "zmax=" << md.getZMax() << endl; 
    cout << "Maximum allowed speed of particles: " << md.getMaximumVelocity() << endl; 
    md.setSaveCount(.5e4);
    cout << "dt=" << md.getTimeStep() << endl;
    
    //This set to colouring based of size and small vectors
    md.setXBallsColourMode(7);
    md.setXBallsVectorScale(1);
        
    //solve
    md.readArguments(argc, argv);
    md.solve();
    
    md.write(std::cout,false);
    md.writeRestartFile();
    md.writeRestartFile();
    
}

References Eigen::bfloat16_impl::atan(), e(), md, MONOLAYER_ORDERED, ONE_FILE, constants::pi, Eigen::bfloat16_impl::pow(), and mathsFunc::square().