[moritzhoefert]

Hello,

I would like to write a post processing tool for Euler-Lagrange-Simulations that reads particle properties and stores for each cell the particle mass weighted velocity of all particles located in the respective cell.

I created a vector field:

    
volVectorField effU
    (   IOobject
        (   "effU",
            mesh.time().timeName(),
            U.mesh(),
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        mesh,
        dimensionedVector("effU", dimensionSet(0,1,-1,0,0,0,0), vector(0,0,0))
    );

I was also able to create a loop through something that is related to the particles:

    
forAll(parcels, pI) 
    {   
        Info << pI << endl;
    }

Within that loop I would like to do someting like this ("pseudocode")

forAll(parcels, pI)
{  Foam::label cellI = parcels[pI].cell();
   effU[cellI] += parcels[pI].Volume()*parcels[pI].rho()*parcels[parcelI].U();
   w[cellI] += parcels[pI].Volume()*parcels[pI].rho();
}
effU = effU/w;

How can I access particle velocities, diameters and IDs of the cells the particles are located in within the above loop or rather how can I make the above idea work in OF?

I am sorry that I ask such a basic question but I have never programmed on that level in OpenFOAM. Thank you for any help!

All the best,
Moritz

Remark:
The template solver for my post processing tool is porousExplicitSourceReactingParcelFoam from OF 2.0.x.

[schf]

hello,
I know its a really old thread, but I have exactly the same problem..

So here again: How is it possible access the particle data in the main function?

int main(int argc, char *argv[])
{
    ...
    solidParticleCloud particles(mesh);
    ...
    forAll(particles, pIter)
    {
     Info << "pIter " << pIter << endl;
     Info << "Particle velocity " << particles.U()[pIter] << endl; //pseudo code
    }
   ...
}

Is it just a syntax problem or do I have a fundamental misunderstanding about accessing the values in the Lagrangian field?

Because the field U is defined in solidParticle.H, but when compiling the code above, I get the error, that

»class Foam::solidParticleCloud« has no element named »U«

Thanks for your help in advance!

[schf]

Maybe this helps somebody:

const Cloud<solidParticle>& a = U.mesh().lookupObject<Cloud<solidParticle>>("defaultCloud");

solidParticleCloud * q = (solidParticleCloud*) &a;
forAllIter(solidParticleCloud, *q, pIter)
{
  Info << "    " << pIter().position() << endl;
  Info << "    " << pIter().d() << endl;
}

It works, but isn't perfect. So if you have a better solution I would be happy for your help!