[aylalisa]

Hi Foamers,

I have a cell label of a cell that belongs to the inner mesh and is located next to a boundary patch.

Can I directly access the face label ID of the adjacent face of that cell which belongs to the boundary mesh, without looping through all faces of that boundary patch.

Aylalisa

[raptorishungry]

Hello aylalisa,

Did you find out how to do this? I am also interested

adithya

[Tobi]

It should not be possible as there is no function for that (as far as I know).
The best solution is to check the face ID before the solver starts and save the ID.
Therefore, you do not have to search for the face ID again and again.

[raptorishungry]

Hello Tobi,

The Lagrangian solver (eg. icoUncoupledKinematicParcelFoam) writes the positions of the particles in /<time>/lagrangian/kinematicCloud/positions.

Aside from writing the x, y and z coordinates, it also writes the cellID of the cell the particle is located at that time step according to the following code fragment in particleIO.C ($FOAM_SRC/lagrangian/basic/particle/particleIO.C)

void Foam::particle::writePosition(Ostream& os) const
{
    if (os.format() == IOstream::ASCII)
    {
        os  << position_ << token::SPACE << celli_;
    }
    else
    {
        os.write(reinterpret_cast<const char*>(&position_), sizeofPosition_);
    }

    // Check state of Ostream
    os.check(FUNCTION_NAME);
}

Now, I want to write the face ID instead of cell ID. The facei_ has the value of faceID as per the following code from $FOAM_SRC/lagrangian/basic/particle/particle.H

        //- Face index if the particle is on a face otherwise -1
        label facei_;

So I changed the lines

os  << position_ << token::SPACE << celli_;

to

os  << position_ << token::SPACE << facei_;

This works in serial (and in every individual processor), but it does not work in parallel. In the sense that, while reconstructing the results from processors, I get -1 everywhere (although its not showing -1 for some particles in the individual processors).

Do you have any idea what is going on? or Do you suggest any work-around to get the solver to write the faceID instead of cellID along with positions?

Thanks.

-adithya

[Tobi]

I am sorry. I do not have any background for parallel-computing in FOAM (only a small one - not worth to talk about).

[Behradskn]

Hi Tobi,

First of all, thank you for all you have done for the OpenFOAM community. I am desperately looking for how to get access to surface values of the face ids produced by faceSet in every time step. Do you have any idea about how can I solve this problem?

Thank you so much!
Behrad