Computes the normal forces due to linear plastic visco elastic interaction.
79 NormalForceInteraction::computeNormalForce();
80 Mdouble radius = 2.0 * NormalForceInteraction::getEffectiveRadius();
81 Mdouble contactArea =
constants::pi * radius * std::max(0.0,NormalForceInteraction::getOverlap());
82 const SpeciesType* species =
static_cast<const SpeciesType*
>(NormalForceInteraction::getBaseSpecies()->getNormalForce());
86 if (pParticle && iParticle)
89 Mdouble heatTransfer = species->getThermalConductivity()
90 * (pParticle->getTemperature() - iParticle->getTemperature())
91 * contactArea / NormalForceInteraction::getDistance();
93 Mdouble mdT = heatTransfer / species->getHeatCapacity()
94 * NormalForceInteraction::getHandler()->getDPMBase()->getTimeStep();
96 iParticle->addTemperature(mdT * iParticle->getInvMass());
double Mdouble
Definition: GeneralDefine.h:34
HeatFluidCoupledSpecies< typename NormalForceInteraction::SpeciesType > SpeciesType
Definition: HeatFluidCoupledInteraction.h:40
Class of particles that store both temperature and liquid volume, which is adapted for the CFD-DEM st...
Definition: HeatFluidCoupledParticle.h:46
void addTemperature(Mdouble temperature)
Definition: ThermalParticle.h:112
const Mdouble pi
Definition: ExtendedMath.h:45