NewtonsCradleSelfTest.cpp File Reference

#include "Mercury3D.h"
#include "StatisticsVector.h"
#include "Particles/BaseParticle.h"
#include "Walls/InfiniteWall.h"
#include <cmath>
#include <iostream>
#include <iomanip>
#include <Species/LinearViscoelasticSpecies.h>

Go to the source code of this file.

Classes
class	NewtonsCradleSelfTest

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

Function Documentation

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

Definition at line 76 of file NewtonsCradleSelfTest.cpp.

References BaseHandler< T >::copyAndAddObject(), Files::getName(), DPMBase::getTimeMax(), NewtonsCradleSelfTest::N, constants::pi, StatisticsVector< T >::set_h(), StatisticsVector< T >::setCGShape(), StatisticsVector< T >::setCGTimeMin(), StatisticsVector< T >::setCGWidth(), DPMBase::setGravity(), File::setName(), Files::setName(), StatisticsVector< T >::setNY(), DPMBase::setParticleDimensions(), Files::setSaveCount(), StatisticsVector< T >::setSuperExact(), DPMBase::setSystemDimensions(), DPMBase::setTimeMax(), StatisticsVector< T >::setTimeMaxStat(), DPMBase::setTimeStep(), StatisticsVector< T >::setXMaxStat(), StatisticsVector< T >::setXMinStat(), StatisticsVector< T >::setYMaxStat(), StatisticsVector< T >::setYMinStat(), StatisticsVector< T >::setZMaxStat(), StatisticsVector< T >::setZMinStat(), DPMBase::solve(), DPMBase::speciesHandler, Files::statFile, and StatisticsVector< T >::statistics_from_fstat_and_data().

{
    NewtonsCradleSelfTest problem;
    problem.setName("NewtonsCradleSelfTest");
    auto species=problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    problem.N=5;    //set the number of particles
    problem.setSystemDimensions(3);
    problem.setParticleDimensions(3);
    species->setDensity(6./constants::pi/8*(2*problem.N+.5)*1.05*1.05);
    problem.setGravity(Vec3D(0.,0.,-1.));
    species->setCollisionTimeAndRestitutionCoefficient(.01,.1,1.);
    problem.setTimeStep(.0001);
    problem.setTimeMax(8.0);
    problem.setSaveCount(1000);
    std::cout << "Relax the packing" << std::endl;
    problem.solve();
    
    std::cout << "Get statistics" << std::endl; 
    std::string str[3];
    str[0]="Gaussian";
    str[1]="HeavisideSphere";
    str[2]="Lucy";
    
    int i=0;
    double w=0.46/3.;
    StatisticsVector<Z> stats1("NewtonsCradleSelfTest");
    stats1.setZMinStat(-.5);
    stats1.setZMaxStat(2*problem.N);
    stats1.setXMinStat(0.5);
    stats1.setXMaxStat(1.55);
    stats1.setYMinStat(0.5);
    stats1.setYMaxStat(1.55);
    std::cout << stats1.getName() + "_" + str[i] + "T_XYZ.stat" << std::endl;
    stats1.statFile.setName(stats1.getName() + "_" + str[i] + "T_XYZ.stat");
    stats1.set_h(0.0005);
    stats1.setCGWidth(w);
    stats1.setSuperExact(false);
    stats1.setCGShape(str[i].c_str());
    stats1.setCGTimeMin(problem.getTimeMax()*1.00000999999);
    stats1.setTimeMaxStat(1e20);
    stats1.statistics_from_fstat_and_data();
    
    for (int i=0; i<3; i++) {
        double w=0.125;
        if (i==1) w*=2;
        if (i==2) w*=2*sqrt(3);
        
        StatisticsVector<Z> stats0("NewtonsCradleSelfTest");
        stats0.setZMinStat(-0.5);
        stats0.statFile.setName(stats0.getName() + "_" + str[i] + "T_XYZ.stat");
        stats0.set_h(0.02);
        stats0.setCGWidth(w);
        stats0.setSuperExact(true);
        stats0.setCGShape(str[i].c_str());
        stats0.setCGTimeMin(problem.getTimeMax()*.999999);
        stats0.setTimeMaxStat(1e20);
        stats0.statistics_from_fstat_and_data();

        StatisticsVector<XYZ> stats2("NewtonsCradleSelfTest");
        stats2.setZMinStat(-0.5);
        stats2.statFile.setName(stats2.getName() + "_" + str[i] + "T_XYZ.stat");
        stats2.set_h(0.02);
        stats2.setNY(1);
        stats2.setCGWidth(w);
        stats2.setSuperExact(true);
        stats2.setCGShape(str[i].c_str());
        stats2.setCGTimeMin(problem.getTimeMax()*.999999);
        stats2.setTimeMaxStat(1e20);
        stats2.statistics_from_fstat_and_data();
        
        /* Current the 2D statistics is not implement this tets will be added when it is.
        StatisticsVector<XZ> stats1("NewtonsCradleSelfTest");
        stats1.setZMinStat(-0.5);
            stats1.statFile.setName(str[i] + "T_XZ.stat");
        stats1.set_h(0.02);
        stats1.setCGWidth(w);
        stats1.setSuperExact(true);
        stats1.setCGShape(str[i].c_str());
        stats1.setCGTimeMin(problem.getTimeMax()*.999999);
        stats1.setTimeMaxStat(1e20);
        stats1.statistics_from_fstat_and_data();
        */
        // should give you Density 5
    }

    return(0);
}