add schillerNaumannDragForce

[Rookie]

Dear OpenLB developers,

If I add schillerNaumannDragForce like this, an error will be reported saying that the number of parameters is wrong, how can I correct it

SuperLatticeInterpPhysVelocity3D<T, DESCRIPTOR> getVel(sLattice, converter);
SuperLatticePhysViscosity3D<T, DESCRIPTOR> dynVisc(sLattice, converter);
SuperLatticeDensity3D<T, DESCRIPTOR> physDensity(sLattice);

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

Regards,
Rookie

[Rookie]

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);

[Adrian]

Looks ok and matches the constructor

SchillerNaumannDragForce3D(SuperLatticeInterpPhysVelocity3D<T, DESCRIPTOR>& getVel, T dynVisc, T physDensity);

[Author]

Posts