Computes the normal forces due to linear plastic visco elastic interaction.
59 NormalForceInteraction::computeNormalForce();
60 Mdouble radius = 2.0 * NormalForceInteraction::getEffectiveRadius();
62 const SpeciesType* species =
static_cast<const SpeciesType*
>(NormalForceInteraction::getBaseSpecies()->getNormalForce());
66 if (pParticle && iParticle)
69 Mdouble heatTransfer = species->getThermalConductivity()
70 * (pParticle->getTemperature() - iParticle->getTemperature())
71 * contactArea / NormalForceInteraction::getDistance();
73 Mdouble mdT = heatTransfer / species->getHeatCapacity()
74 * NormalForceInteraction::getHandler()->getDPMBase()->getTimeStep();
75 pParticle->addTemperature(-mdT * pParticle->getInvMass());
76 iParticle->addTemperature(mdT * iParticle->getInvMass());
HeatFluidCoupledSpecies< typename NormalForceInteraction::SpeciesType > SpeciesType
Definition: HeatFluidCoupledInteraction.h:19
Class of particles that store both temperature and liquid volume, which is adapted for the CFD-DEM st...
Definition: HeatFluidCoupledParticle.h:25
#define max(a, b)
Definition: datatypes.h:23
const Mdouble pi
Definition: ExtendedMath.h:23