71 int locIC = _sLattice.getLoadBalancer().loc(globicFull[0]);
74 int latticeRoundedPosP[3] = { globicFull[1], globicFull[2], globicFull[3] };
77 T omega = _sLattice.getBlock(locIC).get (
78 latticeRoundedPosP[0],
79 latticeRoundedPosP[1],
80 latticeRoundedPosP[2] ).template getField<descriptors::OMEGA>();
83 T nu = this->_converter.getPhysViscosity (
84 (1./omega - 0.5) / descriptors::invCs2<T,Lattice>() );
86 T ReP = 2. * p->getRad() * magU / nu;
87 return ReP > this->_RePmin ? ReP : this->_RePmin;
virtual T operator()(Particle< T > *p, T magU, int globicFull[]) override
Returns the particle Reynolds number. globicFull = { globic, latticeRoundedP[0, .....
NewtonianParticleReynoldsNumber(UnitConverter< T, Lattice > &converter)
Constructor.
Abstract class for particle Reynolds number computation within drag model.
ParticleReynoldsNumberBase(UnitConverter< T, Lattice > &converter)
Constructor.
virtual T operator()(Particle< T > *p, T magU, int globicFull[]) override
Returns the particle Reynolds number. globicFull = { globic, latticeRoundedP[0, .....
PowerLawParticleReynoldsNumber(UnitConverter< T, Lattice > &converter, SuperLattice< T, Lattice > &sLattice)
Constructor.