granular_jet_thomas.cpp File Reference

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

Classes
class	GranularJet

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

Function Documentation

main()

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

{
    GranularJet problem;
 
    // Particle properties
    problem.set_silbert();  
    // Problem parameters
    problem.setName("granular_jet");
 
    //Corrections
    problem.setRoughBottomType(MONOLAYER_DISORDERED); //to save base particles
    problem.setAlignBase(false);
    
    // Chute properties
    problem.setChuteAngle(20);
    problem.setChuteLength(700); //700=40cm
    problem.setChuteWidth(400); //400=24cm
    problem.setMaxFailed(6);
    problem.makeChutePeriodic();
    problem.setHopperDimension(2);
    problem.setHopperLift(250); //250=15cm
    double ExitHeight = 25.0, //25=15mm
           ExitLength = 3.0 * ExitHeight, 
           hopperAngle_ = 45.0, 
           hopperLength_ = 4.0 * ExitLength;
    problem.set_Hopper(ExitLength,ExitHeight,hopperAngle_,hopperLength_);
    problem.setZMax(200);
    
    //solve
    problem.setTimeStep(); 
    problem.setSaveCount(5000);
    problem.dataFile.setFileType(FileType::ONE_FILE);
    problem.fStatFile.setFileType(FileType::NO_FILE);
    problem.restartFile.setFileType(FileType::ONE_FILE);
    problem.setXBallsAdditionalArguments("-v0 ");
    cout << "Maximum allowed speed of particles: " << problem.getMaximumVelocity() 
    << endl; // speed allowed before particles move through each other!
    cout << "dt=" << problem.getTimeStep() << endl;
    
    problem.auto_number();
 
    //solve
    problem.readArguments(argc, argv);
    problem.solve();
}

References MONOLAYER_DISORDERED, NO_FILE, ONE_FILE, and problem.