Foam::vector Foam::DampingModels::Relaxation<CloudType>::velocityCorrection ( typename CloudType::parcelType& p, const scalar deltaT ) const { const tetIndices tetIs(p.cell(), p.tetFace(), p.tetPt(), this->owner().mesh()); const scalar x = deltaT*oneByTimeScaleAverage_->interpolate(p.position(), tetIs); const vector u = uAverage_->interpolate(p.position(), tetIs); return (u - p.U())*x/(x + 2.0); }