[cojua8]

Hello,

I am trying to create a new Magnus force for particles.

for this purpose I need to access the angular velocity in CollidingParcelI.H:

template<class ParcelType>
inline Foam::vector Foam::CollidingParcel<ParcelType>::omega() const
{
    return angularMomentum_/this->momentOfInertia();
}

And the code of the force is this:

template<class CloudType>
Foam::forceSuSp Foam::MagnusForce<CloudType>::calcCoupled
(
    const typename CloudType::parcelType& p,
    const scalar dt,
    const scalar mass,
    const scalar Re,
    const scalar muc
    const vector omega
) const
{
    forceSuSp value(Zero, 0.0);

    vector curlUc =
        curlUcInterp().interpolate(p.coordinates(), p.currentTetIndices());

    scalar Cl = this->Cl(p, curlUc, Re, muc);

    value.Su() =  Cl*0.75*mass*(omega^(p.U()-p.Uc()));
    
    Info << "Calculate Lift force" << endl;

    return value;
}

the compiler compiles the new force, but when running the solver there ir no "Calculate Lift Force" output, which means the solver doesn't get into the force code.

When I delete all the "omega's" the code works correctly (but obviously it doesnt work as needed), so I guess the omega call is not okay