Skip to content

add schillerNaumannDragForce

Viewing 3 posts - 1 through 3 (of 3 total)
  • Author
    Posts
  • #7956
    Rookie
    Participant

    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

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

    #7965
    Adrian
    Keymaster

    Looks ok and matches the constructor

    
    SchillerNaumannDragForce3D(SuperLatticeInterpPhysVelocity3D<T, DESCRIPTOR>& getVel, T dynVisc, T physDensity);
    
Viewing 3 posts - 1 through 3 (of 3 total)
  • You must be logged in to reply to this topic.