HGridUpdateUnitTest.cpp File Reference
#include <iostream>
#include <Species/LinearViscoelasticSpecies.h>
#include "Chute.h"

Classes

class  DPM
 In this file a cubic packing of 5^3 particles in a tri-axial box is created and allowed to settle under small gravity. After that Z statistics are calculated. More...
 

Functions

int main ()
 

Function Documentation

◆ main()

int main ( )
50 {
51 
52  // Problem parameters
53  DPM problem;
54  problem.setFileType(FileType::NO_FILE);
55  problem.setHGridCellOverSizeRatio(2);
56  problem.setHGridUpdateEachTimeStep(false);
57  problem.setName("HGridUpdateUnitTest"); // data output file name
58  problem.setTimeStep(1e-4);
59  problem.setTimeMax(10);
60  problem.setDomain({0,0,0},{1,1,1});
61 
62  // Particle species
63  LinearViscoelasticSpecies species; // initialise species
64  species.setDensity(2000); // particle density
65  species.setCollisionTimeAndRestitutionCoefficient( 5e-3, 0.8, 1);
66  auto s = problem.speciesHandler.copyAndAddObject(species); // assign species to problem species handler
67 
68  SphericalParticle particle;
69  particle.setSpecies(s);
70  particle.setRadius(0.5);
71  particle.setVelocity({1,0,0});
72  auto p = problem.particleHandler.copyAndAddObject(particle);
73 
74  //logger(INFO,"Cell %",)
75 
76 
77  //solve
78  problem.solve();
79  problem.hGridInfo();
80  return 0;
81 }
Array< double, 1, 3 > e(1./3., 0.5, 2.)
@ NO_FILE
file will not be created/read
float * p
Definition: Tutorial_Map_using.cpp:9
void setVelocity(const Vec3D &velocity)
set the velocity of the BaseInteractable.
Definition: BaseInteractable.cc:328
virtual void setRadius(Mdouble radius)
Sets the particle's radius_ (and adjusts the mass_ accordingly, based on the particle's species)
Definition: BaseParticle.cc:548
virtual void setSpecies(const ParticleSpecies *species)
Definition: BaseParticle.cc:798
In this file a cubic packing of 5^3 particles in a tri-axial box is created and allowed to settle und...
Definition: HGridUpdateUnitTest.cpp:9
void setCollisionTimeAndRestitutionCoefficient(Mdouble tc, Mdouble eps, BaseParticle *p)
Sets k, disp such that it matches a given tc and eps for a collision of two copies of particle p.
Definition: LinearViscoelasticNormalSpecies.cc:191
void setDensity(Mdouble density)
Definition: ParticleSpecies.cc:88
A spherical particle is the most simple particle used in MercuryDPM.
Definition: SphericalParticle.h:16
RealScalar s
Definition: level1_cplx_impl.h:130
Constructor for SteadyAxisymAdvectionDiffusion problem
Definition: steady_axisym_advection_diffusion.cc:213

References e(), NO_FILE, p, problem, s, LinearViscoelasticNormalSpecies::setCollisionTimeAndRestitutionCoefficient(), ParticleSpecies::setDensity(), BaseParticle::setRadius(), BaseParticle::setSpecies(), and BaseInteractable::setVelocity().