MercuryDPM  Trunk
ThermalParticle Class Referencefinal

#include <ThermalParticle.h>

Inheritance diagram for ThermalParticle:

## Public Member Functions

ThermalParticle ()
Basic Particle constructor, creates a particle at (0,0,0) with radius, mass and inertia equal to 1. More...

ThermalParticle (const ThermalParticle &p)
Particle copy constructor, which accepts as input a reference to a Particle. It creates a copy of this Particle and all it's information. Usually it is better to use the copy() function for polymorfism. More...

~ThermalParticle () override
Particle destructor, needs to be implemented and checked if it removes tangential spring information. More...

ThermalParticlecopy () const override
Particle copy method. It calls to copy constructor of this Particle, useful for polymorfism. More...

void write (std::ostream &os) const override

std::string getName () const override

void read (std::istream &is) override

Mdouble getTemperature () const

void setTemperature (Mdouble temperature)

void addTemperature (Mdouble temperature)

void setTemperatureDependentDensity (const std::function< double(double)> &temperatureDependentDensity)

const std::function< double(double)> & getTemperatureDependentDensity () const

const std::function< double(double)> & getTimeDependentTemperature () const

void setTimeDependentTemperature (const std::function< double(double)> &timeDependentTemperature)

void actionsAfterTimeStep () override

bool isSphericalParticle () const override

Public Member Functions inherited from BaseParticle
BaseParticle ()
Basic Particle constructor, creates an Particle at (0,0,0) with radius, mass and inertia equal to 1. More...

BaseParticle (const BaseParticle &p)
Particle copy constructor, which accepts as input a reference to a Particle. It creates a copy of this Particle and all it's information. Usually it is better to use the copy() function for polymorfism. More...

BaseParticle (const ParticleSpecies *s)

~BaseParticle () override
Particle destructor, needs to be implemented and checked if it removes tangential spring information. More...

virtual Mdouble getVolume () const
Get Particle volume function, which required a reference to the Species vector. It returns the volume of the Particle. More...

void fixParticle ()
Fix Particle function. It fixes a Particle by setting its inverse mass and inertia and velocities to zero. More...

bool isFixed () const override
Is fixed Particle function. It returns whether a Particle is fixed or not, by checking its inverse Mass. More...

bool isMPIParticle () const
Indicates if this particle is a ghost in the MPI domain. More...

void setMPIParticle (bool flag)
Flags the mpi particle status. More...

bool isInMPIDomain ()
Indicates if the particle is in the communication zone of the mpi domain. More...

void setInMPIDomain (bool flag)
Flags the status of the particle if wether it is in the communication zone or not. More...

bool isInPeriodicDomain () const
Indicates if the particle is in the periodic boundary communication zone. More...

void setInPeriodicDomain (bool flag)
Flags the status of the particle wether it is in the periodic communication zone or not. More...

bool isPeriodicGhostParticle () const
Indicates if this particle is a ghost in the periodic boundary. More...

void setPeriodicGhostParticle (bool flag)
Flags the status of the particle to be a ghost in periodic boundary or not. More...

bool isMaserParticle () const
Indicates if this particle belongs to the maser boundary. More...

void setMaserParticle (bool flag)
Flags the status of the particle if it belongs to the maser boundary or not. More...

void setCommunicationComplexity (unsigned complexity)
Set the communication complexity of the particle. More...

unsigned getCommunicationComplexity ()
Obtains the communication complexity of the particle. More...

void setPeriodicComplexity (std::vector< int > complexity)
Set the periodic communication complexity of the particle. More...

void setPeriodicComplexity (int index, int value)
Set the periodic communication complexity of the particle. More...

const std::vector< int > & getPeriodicComplexity ()
Obtains the periodic communication complexity of the particle. More...

void setPreviousPeriodicComplexity (std::vector< int > complexity)
Set the previous periodic communication complexity of the paritcle. More...

const std::vector< int > & getPreviousPeriodicComplexity () const
Sets the previous periodic communication complexity of the particle. More...

int getPeriodicComplexity (int index)
Gets the periodic communication complexity of a certain boundary. More...

void unfix ()
Unfix Particle function, which required a reference to the Species vector. It unfixes a Particle by computing the Particles mass and inertia. More...

virtual void oldRead (std::istream &is)

virtual void setInfo (Mdouble info)
Sets some user-defined information about this object (by default, species ID). More...

virtual Mdouble getInfo () const
Returns some user-defined information about this object (by default, species ID). More...

void printHGrid (std::ostream &os) const
Adds particle's HGrid level and cell coordinates to an ostream. More...

unsigned int getHGridLevel () const
Returns particle's HGrid level. More...

BaseParticlegetHGridNextObject () const
Returns pointer to next object in particle's HGrid level & cell. More...

BaseParticlegetHGridPrevObject () const
Returns pointer to previous object in particle's HGrid level & cell. More...

int getHGridX () const
Returns particle's HGrid cell X-coordinate. More...

int getHGridY () const
Returns particle's HGrid cell Y-coordinate. More...

int getHGridZ () const
Returns particle's HGrid cell Z-coordinate. More...

MatrixSymmetric3D getInvInertia () const
Returns the inverse of the particle's inertia tensor. More...

Mdouble getInvMass () const override
Returns the inverse of the particle's mass. More...

Mdouble getCurvature (const Vec3D &labFixedCoordinates) const override

Mdouble getKineticEnergy () const
Calculates the particle's translational kinetic energy. More...

Mdouble getRotationalEnergy () const
Calculates the particle's rotational kinetic energy. More...

Mdouble getGravitationalEnergy () const
Calculates the particle's gravitational energy. More...

Mdouble getMass () const
Returns the particle's mass. More...

Vec3D getMomentum () const

MatrixSymmetric3D getInertia () const

Vec3D getAngularMomentum () const

BaseParticlegetPeriodicFromParticle () const
Returns the 'original' particle this one's a periodic copy of. More...

Mdouble getRadius () const
Returns the particle's radius. More...

Mdouble getMaxInteractionRadius () const
Returns the particle's interaction radius, which might be different from radius_ (e.g., when dealing with wet particles) More...

Mdouble getInteractionDistance (const BaseInteractable *i) const
Returns the interactionDistance_ of the mixed species of this particle and the particle or wall i. More...

Mdouble getSumOfInteractionRadii (const BaseParticle *particle) const
returns the sum of the radii plus the interactionDistance More...

Mdouble getWallInteractionRadius (const BaseWall *wall) const
returns the radius plus the interactionDistance More...

const Vec3DgetDisplacement () const
Returns the particle's displacement relative to the previous time step. More...

const Vec3DgetPreviousPosition () const
Returns the particle's position in the previous time step. More...

const Vec3D getDisplacement2 (Mdouble xmin, Mdouble xmax, Mdouble ymin, Mdouble ymax, Mdouble zmin, Mdouble zmax, Mdouble t) const

virtual void setInertia ()

void setInertia (MatrixSymmetric3D inertia)
Sets the particle's inertia_ (and adjusts invInertia_ accordingly) More...

void setInverseInertia (MatrixSymmetric3D inverseInertia)
Sets the particle's inertia_ (and adjusts invInertia_ accordingly) More...

void setInfiniteInertia ()
Sets the particle's inertia_ to 'infinite' (1e20) and its invInertia_ to 0. More...

void setPeriodicFromParticle (BaseParticle *p)
Assigns the pointer to the 'original' particle this one's a periodic copy of (used in periodic boundary condition implementations). More...

void setHGridX (const int x)
Sets the particle's HGrid cell X-coordinate. More...

void setHGridY (const int y)
Sets the particle's HGrid cell Y-coordinate. More...

void setHGridZ (const int z)
Sets the particle's HGrid cell Z-coordinate. More...

void setHGridLevel (const unsigned int level)
Sets the particle's HGrid level. More...

void setHGridNextObject (BaseParticle *p)
Sets the pointer to the next object in the particle's HGrid cell & level. More...

void setHGridPrevObject (BaseParticle *p)
Sets the pointer to the previous object in the particle's HGrid cell & level. More...

Sets the particle's radius_ (and adjusts the mass_ accordingly, based on the particle's species) More...

virtual Vec3D getAxes () const
Only ustilised in case of superquadric particles. Had to create a virtual function to allow function access in writeVTK function in the particle handler. More...

virtual Mdouble getExponentEps1 () const
Only ustilised in case of superquadric particles. Had to create a virtual function to allow function access in writeVTK function in the particle handler. More...

virtual Mdouble getExponentEps2 () const
Only ustilised in case of superquadric particles. Had to create a virtual function to allow function access in writeVTK function in the particle handler. More...

virtual void setAxes (const Vec3D &axes)
Only ustilised in case of superquadric particles. More...

virtual void setExponents (const Mdouble &eps1, const Mdouble &eps2)
Only ustilised in case of superquadric particles. More...

MERCURY_DEPRECATED void setMass (Mdouble mass)
Sets the particle's mass. More...

void setMassForP3Statistics (Mdouble mass)
Sets the particle's mass This function should not be used, but is necessary to extend the CG toolbox to non-spherical particles. More...

void setDisplacement (const Vec3D &disp)
Sets the particle's displacement (= difference between current position and that of the previous time step) More...

void setPreviousPosition (const Vec3D &pos)
Sets the particle's position in the previous time step. More...

void movePrevious (const Vec3D &posMove)
Adds a vector to the particle's previousPosition_. More...

void accelerate (const Vec3D &vel)
Increases the particle's velocity_ by the given vector. More...

void angularAccelerate (const Vec3D &angVel)
Increases the particle's angularVelocity_ by the given vector. More...

Adds a vector to the particle's displacement_. More...

void setHandler (ParticleHandler *handler)
Sets the pointer to the particle's ParticleHandler. More...

ParticleHandlergetHandler () const
Returns pointer to the particle's ParticleHandler. More...

BaseInteractiongetInteractionWith (BaseParticle *P, unsigned timeStamp, InteractionHandler *interactionHandler) override
Checks if particle is in interaction with given particle P, and if so, returns vector of pointer to the associated BaseInteraction object (else returns empty vector). More...

virtual bool isInContactWith (const BaseParticle *P) const
Get whether or not this particle is in contact with the given particle. More...

void integrateBeforeForceComputation (double time, double timeStep)
First step of Velocity Verlet integration. More...

void integrateAfterForceComputation (double time, double timeStep)
Second step of Velocity Verlet integration. More...

unsigned int getParticleDimensions () const
Returns the particle's dimensions (either 2 or 3). More...

MERCURY_DEPRECATED void setIndSpecies (unsigned int indSpecies) override

void setSpecies (const ParticleSpecies *species)

virtual unsigned getNumberOfFieldsVTK () const

virtual std::string getTypeVTK (unsigned i) const

virtual std::string getNameVTK (unsigned i) const

virtual std::vector< MdoublegetFieldVTK (unsigned i) const

const HGridCellgetHGridCell () const

virtual void computeMass (const ParticleSpecies &s)
Computes the particle's (inverse) mass and inertia. More...

Public Member Functions inherited from BaseInteractable
BaseInteractable ()
Default BaseInteractable constructor. More...

BaseInteractable (const BaseInteractable &p)
Copy constructor. More...

~BaseInteractable () override
Destructor, it simply destructs the BaseInteractable and all the objects it contains. More...

unsigned int getIndSpecies () const
Returns the index of the species associated with the interactable object. More...

const ParticleSpeciesgetSpecies () const
Returns a pointer to the species of this BaseInteractable. More...

void setSpecies (const ParticleSpecies *species)
Sets the species of this BaseInteractable. More...

const Vec3DgetForce () const
Returns the force on this BaseInteractable. More...

const Vec3DgetTorque () const
Returns the torque on this BaseInteractable. More...

void setForce (const Vec3D &force)
Sets the force on this BaseInteractable. More...

void setTorque (const Vec3D &torque)
Sets the torque on this BaseInteractable. More...

Adds an amount to the force on this BaseInteractable. More...

Adds an amount to the torque on this BaseInteractable. More...

void resetForceTorque (int numberOfOMPthreads)

void sumForceTorqueOMP ()

const Vec3DgetPosition () const
Returns the position of this BaseInteractable. More...

const QuaterniongetOrientation () const
Returns the orientation of this BaseInteractable. More...

void setPosition (const Vec3D &position)
Sets the position of this BaseInteractable. More...

void setOrientationViaNormal (Vec3D normal)
Sets the orientation of this BaseInteractable by defining the vector that results from the rotation of the (1,0,0) vector. More...

void setOrientationViaEuler (Vec3D eulerAngle)
Sets the orientation of this BaseInteractable by defining the euler angles. More...

void setOrientation (const Quaternion &orientation)
Sets the orientation of this BaseInteractable. More...

virtual void move (const Vec3D &move)
Moves this BaseInteractable by adding an amount to the position. More...

virtual void rotate (const Vec3D &angularVelocityDt)
Rotates this BaseInteractable. More...

const std::vector
< BaseInteraction * > &
getInteractions () const
Returns a list of interactions which belong to this interactable. More...

void addInteraction (BaseInteraction *I)
Adds an interaction to this BaseInteractable. More...

bool removeInteraction (BaseInteraction *I)
Removes an interaction from this BaseInteractable. More...

void copyInteractionsForPeriodicParticles (const BaseInteractable &p)
Copies interactions to this BaseInteractable whenever a periodic copy made. More...

void setVelocity (const Vec3D &velocity)
set the velocity of the BaseInteractable. More...

void setAngularVelocity (const Vec3D &angularVelocity)
set the angular velocity of the BaseInteractble. More...

void addVelocity (const Vec3D &velocity)
adds an increment to the velocity. More...

void addAngularVelocity (const Vec3D &angularVelocity)
add an increment to the angular velocity. More...

virtual const Vec3DgetVelocity () const
Returns the velocity of this interactable. More...

virtual const Vec3DgetAngularVelocity () const
Returns the angular velocity of this interactable. More...

void setPrescribedPosition (const std::function< Vec3D(double)> &prescribedPosition)
Allows the position of an infinite mass interactable to be prescribed. More...

void applyPrescribedPosition (double time)
Computes the position from the user defined prescribed position function. More...

void setPrescribedVelocity (const std::function< Vec3D(double)> &prescribedVelocity)
Allows the velocity of an infinite mass interactable to be prescribed. More...

void applyPrescribedVelocity (double time)
Computes the velocity from the user defined prescribed velocity function. More...

void setPrescribedOrientation (const std::function< Quaternion(double)> &prescribedOrientation)
Allows the orientation of the infinite mass interactbale to be prescribed. More...

void applyPrescribedOrientation (double time)
Computes the orientation from the user defined prescribed orientation function. More...

void setPrescribedAngularVelocity (const std::function< Vec3D(double)> &prescribedAngularVelocity)
Allows the angular velocity of the infinite mass interactable to be prescribed. More...

void applyPrescribedAngularVelocity (double time)
Computes the angular velocity from the user defined prescribed angular velocity. More...

virtual const Vec3D getVelocityAtContact (const Vec3D &contact) const
Returns the velocity at the contact point, use by many force laws. More...

void integrateBeforeForceComputation (double time, double timeStep)
This is part of integrate routine for objects with infinite mass. More...

void integrateAfterForceComputation (double time, double timeStep)
This is part of the integration routine for objects with infinite mass. More...

Public Member Functions inherited from BaseObject
BaseObject ()=default
Default constructor. More...

BaseObject (const BaseObject &p)=default
Copy constructor, copies all the objects BaseObject contains. More...

virtual ~BaseObject ()=default
virtual destructor More...

virtual void moveInHandler (unsigned int index)
Except that it is virtual, it does the same thing as setIndex() does. More...

void setIndex (unsigned int index)
Allows one to assign an index to an object in the handler/container. More...

void setId (unsigned long id)
Assigns a unique identifier to each object in the handler (container) which remains constant even after the object is deleted from the container/handler. More...

unsigned int getIndex () const
Returns the index of the object in the handler. More...

unsigned int getId () const
Returns the unique identifier of any particular object. More...

void setGroupId (unsigned groupId)

unsigned getGroupId () const

## Private Attributes

std::function< double(double
temperature)>
timeDependentTemperature_

Mdouble temperature_

## Additional Inherited Members

Protected Attributes inherited from BaseParticle

Mdouble invMass_

MatrixSymmetric3D invInertia_
Inverse Particle mass (for computation optimization) More...

## Detailed Description

Definition at line 35 of file ThermalParticle.h.

## Constructor & Destructor Documentation

 ThermalParticle::ThermalParticle ( )

Basic Particle constructor, creates a particle at (0,0,0) with radius, mass and inertia equal to 1.

default constructor, creates an Particle at (0,0,0) with radius, mass and heatinertia equal to 1

Definition at line 38 of file ThermalParticle.cc.

References temperature_.

Referenced by copy().

39 {
40  temperature_ = 0;
41  //temperatureDependentDensity_
42 }
 ThermalParticle::ThermalParticle ( const ThermalParticle & p )

Particle copy constructor, which accepts as input a reference to a Particle. It creates a copy of this Particle and all it's information. Usually it is better to use the copy() function for polymorfism.

Constructor that copies most of the properties of the given particle. Please note that not everything is copied, for example the position in the HGrid is not determined yet by the end of this constructor. It also does not copy the interactions and the pointer to the handler that handles this particle. Use with care.

Parameters
 [in,out] p Reference to the ThermalParticle this one should become a copy of.

Definition at line 52 of file ThermalParticle.cc.

References temperature_, and timeDependentTemperature_.

53  : BaseParticle(p)
54 {
57 }
std::function< double(double temperature)> timeDependentTemperature_
BaseParticle()
Basic Particle constructor, creates an Particle at (0,0,0) with radius, mass and inertia equal to 1...
Definition: BaseParticle.cc:33
 ThermalParticle::~ThermalParticle ( )
overridedefault

Particle destructor, needs to be implemented and checked if it removes tangential spring information.

Destructor. It asks the ParticleHandler to check if this was the smallest or largest particle and adjust itself accordingly.

## Member Function Documentation

 void ThermalParticle::actionsAfterTimeStep ( )
overridevirtual

Reimplemented from BaseParticle.

Definition at line 129 of file ThermalParticle.cc.

130 {
132  {
133  temperature_ = timeDependentTemperature_(getHandler()->getDPMBase()->getTime());
134  }
136  {
138  radius_ = getRadius() * cbrt(getMass() / (getVolume() * density));
139  }
140 }
virtual Mdouble getVolume() const
Get Particle volume function, which required a reference to the Species vector. It returns the volume...
double Mdouble
Definition: GeneralDefine.h:34
const std::function< double(double)> & getTemperatureDependentDensity() const
const ParticleSpecies * getSpecies() const
Returns a pointer to the species of this BaseInteractable.
ParticleHandler * getHandler() const
Returns pointer to the particle's ParticleHandler.
const std::function< double(double)> & getTemperatureDependentDensity() const
Mdouble getMass() const
Returns the particle's mass.
Definition: BaseParticle.h:322
Returns the particle's radius.
Definition: BaseParticle.h:345
std::function< double(double temperature)> timeDependentTemperature_
Definition: BaseParticle.h:652
 void ThermalParticle::addTemperature ( Mdouble temperature )

Definition at line 124 of file ThermalParticle.cc.

References temperature_.

125 {
126  temperature_ += temperature;
127 }
 ThermalParticle * ThermalParticle::copy ( ) const
overridevirtual

Particle copy method. It calls to copy constructor of this Particle, useful for polymorfism.

Copy method. Uses copy constructor to create a copy on the heap. Useful for polymorphism.

Returns
pointer to the particle's copy

Implements BaseParticle.

Definition at line 71 of file ThermalParticle.cc.

References ThermalParticle().

72 {
73  return new ThermalParticle(*this);
74 }
ThermalParticle()
Basic Particle constructor, creates a particle at (0,0,0) with radius, mass and inertia equal to 1...
 std::string ThermalParticle::getName ( ) const
overridevirtual

Returns the name of the object; in this case 'ThermalParticle'.

Returns
The object name.

Reimplemented from BaseParticle.

Definition at line 92 of file ThermalParticle.cc.

93 {
94  return "ThermalParticle";
95 }
 Mdouble ThermalParticle::getTemperature ( ) const

Definition at line 114 of file ThermalParticle.cc.

References temperature_.

115 {
116  return temperature_;
117 }
 const std::function& ThermalParticle::getTemperatureDependentDensity ( ) const

Referenced by actionsAfterTimeStep().

 const std::function< double(double)> & ThermalParticle::getTimeDependentTemperature ( ) const

Definition at line 142 of file ThermalParticle.cc.

References timeDependentTemperature_.

143 {
145 }
std::function< double(double temperature)> timeDependentTemperature_
 bool ThermalParticle::isSphericalParticle ( ) const
inlineoverridevirtual

Reimplemented from BaseParticle.

Definition at line 81 of file ThermalParticle.h.

81 {return true;}
 void ThermalParticle::read ( std::istream & is )
overridevirtual

Particle read function. Has an std::istream as argument, from which it extracts the radius_, invMass_ and invInertia_, respectively. From these the mass_ and inertia_ are deduced. An additional set of properties is read through the call to the parent's method BaseParticle::read().

Parameters
 [in,out] is input stream with particle properties.

Reimplemented from BaseParticle.

Definition at line 107 of file ThermalParticle.cc.

References BaseParticle::read(), and temperature_.

108 {
110  std::string dummy;
111  is >> dummy >> temperature_;
112 }
void read(std::istream &is) override
Particle read function, which accepts an std::istream as input.
 void ThermalParticle::setTemperature ( Mdouble temperature )

Definition at line 119 of file ThermalParticle.cc.

References temperature_.

120 {
121  temperature_ = temperature;
122 }
 void ThermalParticle::setTemperatureDependentDensity ( const std::function< double(double)> & temperatureDependentDensity )
 void ThermalParticle::setTimeDependentTemperature ( const std::function< double(double)> & timeDependentTemperature )

Definition at line 147 of file ThermalParticle.cc.

References INFO, logger, temperature_, and timeDependentTemperature_.

148 {
149  timeDependentTemperature_ = timeDependentTemperature;
150  temperature_ = timeDependentTemperature(0);
151  logger(INFO, "Setting initial temperature to %", temperature_);
152 }
Logger< MERCURY_LOGLEVEL > logger("MercuryKernel")
std::function< double(double temperature)> timeDependentTemperature_
 void ThermalParticle::write ( std::ostream & os ) const
overridevirtual

ThermalParticle print method, which accepts an os std::ostream as input. It prints human readable ThermalParticle information to the std::ostream.

Parameters
 [in,out] os stream to which the info is written

Reimplemented from BaseParticle.

Definition at line 82 of file ThermalParticle.cc.

References temperature_, and BaseParticle::write().

83 {
85  os << " temperature " << temperature_;
86 }
void write(std::ostream &os) const override
Particle print function, which accepts an std::ostream as input.

## Member Data Documentation

 Mdouble ThermalParticle::temperature_
private

Definition at line 90 of file ThermalParticle.h.

 std::function ThermalParticle::timeDependentTemperature_
private

Change this function to let the temperature be time-dependent.

Definition at line 88 of file ThermalParticle.h.

The documentation for this class was generated from the following files: