Skip to content

Reply To: add schillerNaumannDragForce

#7959
Rookie
Participant

I’ve solved the problem, and if I’m right, it looks like this:

SuperLatticeInterpPhysVelocity3D<T, DESCRIPTOR> getVel(sLattice, converter);
T dynVisc = converter.getPhysViscosity() * converter.getPhysDensity();
T physDensity = converter.getPhysDensity();

auto schillerNaumannDragForce = std::make_shared<SchillerNaumannDragForce3D<T, PARTICLE, DESCRIPTOR>>(getVel, dynVisc, physDensity);
supParticleSystem.addForce(schillerNaumannDragForce);