46 typename std::deque<PARTICLETYPE<T> >::iterator p,
int pInt,
52 T fluidVel[3] = {0., 0., 0.};
53 _getVel(fluidVel, &p->getPos()[0], p->getCuboid());
57 T particleRe = 2. * p->getRad() *
abs(fluidVelAbs - partVelAbs) * _physDensity / _dynVisc;
58 T coeffSN = 1. + 0.15 *
util::pow(particleRe, 0.687);
60 for (
int i = 0; i < 3; i++) {
62 force[i] = -1. * 6 *
M_PI * p->getRad() * (p->getVel()[i] - fluidVel[i]) * _dynVisc * coeffSN;
63 p->getForce()[i] += force[i] ;