[cojua8]

Hello,

I'm working on LPT with basicKinematicCollidingCloud and I'm trying to program a new Magnus force.

for this, I created new files based on existing force, and changing the return value to:

 
value.Su() =  Cl*0.75*mass*(p.omega()^(p.U()-p.Uc()));

The problem is, when compiling I get the following error:

error: ‘const parcelType {aka const class Foam::KinematicParcel<Foam::particle>}’ has no member named ‘omega’
     vector velang= p.omega();

After looking in the files why I get this, I found it is because the omega member belongs to CollidingCloud and not KinematicCloud.

basicKinematicCollidingCloud has this form:

namespace Foam
{
    typedef CollidingCloud
    <
        KinematicCloud
        <
            Cloud
            <
                basicKinematicCollidingParcel
            >
        >
    > basicKinematicCollidingCloud;
}

So I think (I don't have a big understanding on this) that there must be a way to use p.omega() in a particle force.

Thank you for your help

EDIT:

I found the solution and it is quite easy:

Create a copy of makeParcelForces.H and add the force with the omega term.

Then, in file makeBasicKinematicCollidingParcelSubmodels.C change

#include "makeParcelForcesPM.H"

with the new file name.

This way, when calling wmake libso there will be no error.

The cause of the error was that all the parcel types call the same files for building themselves, and some of them don't have the properties required.