[IsAv]

Hello, this must be something very simple but I am lost reading the api description, the question is, is the LiftForce included by default in the particleForces for mppiccloud? I am using of12. I checked the tutorials and all of them include at the most drag and gravity. Thanks a lot.

[IsAv]

Ok there is no default lift force, I followed some other posts, and arrived to these dictionaries modifications in cloudProperties:

    
interpolationSchemes
    {
        rho.air         cell;
        U.air           cellPointWallModified;
        mu.air          cell;
        alpha.air       cell;
        curlUcDt        cell; //<- this
   }

and

    particleForces
    {
        ErgunWenYuDrag
        {
            alphac alpha.air;
        }
         gravity;
       //this ->

         SaffmanMeiLiftForce
        {
            U    U.air;
        }
    }

to add the SaffmanMeiLiftForce. This runs as longer as there are no particles in the domain, as soon as particles enter, it crushes with the error:

--> FOAM FATAL ERROR: 
Carrier phase curlUc interpolation object not set

I thought I was setting the interpolation object by giving the interpolation scheme of curlUcDt with the code above but it seems not. Does anybody know how I should set the object?
Thanks a lot

[dscian]

You probably need to define curlUc instead of curlUcDt in interpolationSchemes.

[IsAv]

Originally Posted by dscian

You probably need to define curlUc instead of curlUcDt in interpolationSchemes.

I have tried that in several ways I still get the same fatal error

[IsAv]

Originally Posted by IsAv

Ok there is no default lift force

Ok I am going back may be, there is a LiftForce that calculates (from LiftForce.C):

     forceSuSp value(Zero, 0.0); 
     const vector curlUc =
         curlUcInterp().interpolate
          (
              p.coordinates(),
              p.currentTetIndices(td.mesh)
          );
 scalar Cl = this->Cl(p, td, curlUc, Re, muc);
value.Su() = mass/p.rho()*td.rhoc()*Cl*((td.Uc() - p.U())^curlUc)

to which the TomoyamaLift and SaffmanMeiLift contribute, I am prone to think that LiftForce is there by default since it is not offered as an option by ParticleForces but still not sure. I don't what Cl is so I can not confirm, it is defined with value zero in LiftForce.C:

 Foam::scalar Foam::LiftForce<CloudType>::LiftForce::Cl
  (
      const typename CloudType::parcelType& p,
      const typename CloudType::parcelType::trackingData& td,
      const vector& curlUc,
      const scalar Re,
      const scalar muc
  ) const
  {
      // dummy
      return 0.0;
  }

So the LiftForce is zero? How could I confirm that the lift is there?

[IsAv]

For the record, it works with collidingCloud, that's enough for me.