[alainislas]

Hi OpenFOAM friends,


I'm struggling with making my runTime function to work under parallel running. I know there are the returnReduce and gatherList functions, but don't know how to make it work appropriately. Under serial run is OK.

Basically I have an injection of coalParcels trraveling through my domain. Whenever they fall inside a specified cellZone I want the function to remove the coalParcels.
I'm wondering how to loop over the processors to get all the coalParcels and all the cells from my cellZone ?

Any help / ideas are kindly appreciated

 codeExecute
#{

        word cellZoneName = "fluid_porousCells";

         const label cellZoneID = mesh().cellZones().findZoneID(cellZoneName);
        const labelList& cellLabels = mesh().cellZones()[cellZoneID];

        const reactingMultiphaseCloud& a = obr_.lookupObject<reactingMultiphaseCloud>("coalCloud1");
        basicReactingMultiphaseCloud * q = (basicReactingMultiphaseCloud*) &a;
        Info << returnReduce(q->size(), sumOp<scalar>()) << endl;

        typedef typename Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::KinematicParcel<Foam::particle>>>> parcelType;

        forAllIter(basicReactingMultiphaseCloud, *q, pIter)
         {
            if ( mesh().cellZones()[cellZoneID].whichCell(pIter().cell()) != -1 )
            {
                parcelType& p = pIter();
                q->deleteParticle(p);
            }
        }
#};

[alainislas]

https://gifyu.com/image/S87H4