Domino.cpp File Reference

#include "Mercury3D.h"
#include "Walls/InfiniteWall.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"
#include "Particles/ClumpParticle.h"
#include "../ClumpHeaders/ClumpInput.h"
#include "../ClumpHeaders/Mercury3DClump.h"
#include <stdlib.h>
#include "Boundaries/PeriodicBoundary.h"

Classes
struct	dominoes
class	ChangingTOIParticle

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

Variables
dominoes	D

Function Documentation

main()

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

{
 
    int NumParams = 1;
 
    // Get parameters passed through the command line
    /*
    std::vector <Mdouble> params(0);
    for (int i = 0; i<NumParams; i++) {
        std::string a;
        Mdouble param = stod(helpers::readFromCommandLine(argc, argv, "-p"+std::to_string(i), a));
        params.push_back(param);
        std::cout<<params[i]<<std::endl;
    }
*/
    ChangingTOIParticle problem;
 
    // D.S_dom = 0.2 * (2 * D.k_peb * D.S_peb) + 0.8 * (2 * D.k_peb * D.S_peb) * (params[0]/100);
 
    // Domino species
    auto species0 = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
 
    // Wall species
    auto species1 = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
 
    // Mixed
    auto species01 = problem.speciesHandler.getMixedObject(species0, species1);
 
 
 
    species0->setDensity(1.0); // sets the species type-0 density
    //species->setConstantRestitution(0);
    species0->setSlidingFrictionCoefficient(0.0);
    species0->setSlidingStiffness(5e5);
    species0->setRollingFrictionCoefficient(0.0);
    species0->setRollingStiffness(5e5);
    species0->setDissipation(1.0);
    species0->setStiffness(1e6);
    const Mdouble collisionTime = species0->getCollisionTime(D.mass);
 
    species01->setSlidingFrictionCoefficient(0.2);
    species01->setSlidingStiffness(5e5);
    species01->setRollingFrictionCoefficient(0.1);
    species01->setDissipation(1.0);
    species01->setRollingStiffness(5e5);
    species01->setStiffness(1e6);
 
 
 
    problem.setClumpDamping(5);
    problem.setTimeStep(collisionTime / 50.0);
 
    // Quick demonstration
    problem.setSaveCount(500);
    problem.setTimeMax(10);
 
 
    // For time averaging featuring equipartition - takes a while
    // problem.setSaveCount(50000);
    // problem.setTimeMax(100000);
 
    problem.removeOldFiles();
    problem.solve();
 
    if ((D.started)&&(D.finished)) D.velocity = D.S_peb * (D.N_fin - D.N_ini)  / (D.T2-D.T1);
 
    std::cout<<"DOMINO VELOCITY: "<<D.velocity<<std::endl;
    std::cout<<"N_ini: "<<D.N_ini<<std::endl;
    std::cout<<"N_fin: "<<D.N_fin<<std::endl;
    std::cout<<"T1: "<<D.T1<<std::endl;
    std::cout<<"T2: "<<D.T2<<std::endl;
 
    // Add velocity to log
    std::ofstream funct("velocity.txt", std::ios_base::app | std::ios_base::out);
    funct << D.velocity <<"\n";
    funct.close();
 
 
 
    return 0;
}

References D, dominoes::finished, dominoes::mass, dominoes::N_fin, dominoes::N_ini, out(), problem, dominoes::S_peb, dominoes::started, dominoes::T1, dominoes::T2, and dominoes::velocity.

Variable Documentation

D

dominoes D

Referenced by ChangingTOIParticle::actionsAfterTimeStep(), DomainHandler::addObject(), aliasing_with_resize(), bdcsvd_check_convergence(), ChangingTOIParticle::ChangingTOIParticle(), oomph::RefineableQElement< 3 >::check_integrity(), oomph::OcTreeForest::construct_up_right_equivalents(), PSD::convertPdfPhiToPdfMeter(), RachelsAdvectionDiffusion::diff_function(), Global_Physical_Variables::diff_function(), oomph::OcTree::edge_neighbour_is_face_neighbour(), EIGEN_DECLARE_TEST(), UnsteadyHeatProblem< ELEMENT >::enforce_time_periodic_boundary_conditions(), oomph::OcTree::faces_of_common_edge(), oomph::AdvectionDiffusionReactionEquations< NREAGENT, DIM >::fill_in_generic_residual_contribution_adv_diff_react(), oomph::RefineableAdvectionDiffusionReactionEquations< NREAGENT, DIM >::fill_in_generic_residual_contribution_adv_diff_react(), oomph::GeneralisedAdvectionDiffusionEquations< DIM >::fill_in_generic_residual_contribution_cons_adv_diff(), oomph::RefineableGeneralisedAdvectionDiffusionEquations< DIM >::fill_in_generic_residual_contribution_cons_adv_diff(), oomph::RefineableGeneralisedAxisymAdvectionDiffusionEquations::fill_in_generic_residual_contribution_cons_axisym_adv_diff(), oomph::OcTreeForest::find_neighbours(), generalized_eigensolver_real(), oomph::GeneralisedAdvectionDiffusionEquations< DIM >::get_diff_cons_adv_diff(), oomph::get_diff_cons_axisym_adv_diff(), RachelsAdvectionDiffusion::get_exact_u(), oomph::RefineableQElement< 3 >::get_face_bcs(), oomph::RefineableSolidQElement< 3 >::get_face_solid_bcs(), oomph::GeneralisedAdvectionDiffusionEquations< DIM >::get_total_flux(), oomph::get_total_flux(), IntersectionOfWalls::getDistanceAndNormal(), oomph::OcTree::gteq_face_neighbour(), oomph::RefineableQElement< 3 >::interpolated_zeta_on_face(), oomph::WarpedCubeDomain::macro_element_boundary(), oomph::QuarterPipeDomain::macro_element_boundary(), oomph::QuarterTubeDomain::macro_element_boundary(), oomph::TubeDomain::macro_element_boundary(), oomph::EighthSphereDomain::macro_element_boundary(), main(), oomph::PRefineableQElement< 3, INITIAL_NNODE_1D >::node_created_by_neighbour(), oomph::RefineableQElement< 3 >::node_created_by_neighbour(), oomph::PRefineableQElement< 3, INITIAL_NNODE_1D >::node_created_by_son_of_neighbour(), oomph::RefineableQElement< 3 >::oc_hang_helper(), oomph::PRefineableQElement< 3, INITIAL_NNODE_1D >::oc_hang_helper(), oomph::GeneralisedAdvectionDiffusionEquations< 2 >::output_fct(), oomph::RefineableQSpectralElement< 3 >::rebuild_from_sons(), oomph::OcTree::rotate(), Eigen::internal::compute_inverse_size4< Architecture::Target, float, MatrixType, ResultType >::run(), oomph::RefineableQElement< 3 >::setup_hang_for_value(), oomph::RefineableQElement< 3 >::setup_hanging_nodes(), oomph::OcTree::setup_static_data(), AxisymmetricHopper::setupInitialConditions(), ChangingTOIParticle::setupInitialConditions(), test1(), test2(), test_dynamic_bool(), test_lazy_single(), and Eigen::CholmodSimplicialLDLT< MatrixType_, UpLo_ >::vectorD().