[Chrisi1984]

Hello all,

I would like to be able to stop lagrangian particles in a specified cell zone.

In the lagrangian folder "lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel" I can find the calculation of the particle velocity "this->U_" in the file "ReactingMultiphaseParcel.C".

How to implement there anything that I gain "this->U_=0" for all particles those are in a specified cell zone?

Is this the right environment to assign the velocity zero to all particles in the specified cell zone (if I will use parcels of the type ReactingMultiphaseParcel)?

Kind regards,
Chrisi

[Chrisi1984]

Hello,

it seems that nobody has an idea how one can stop parcels in a specified cell zone.

Due to the case that in a porous media the flow is orientated in one direction. It would also help me if I would be able to introduce an if clause that stops all parcels those are in a cell where the gas flow velocity is nearly zero in a given direction.

Does anybody of you has an idea how to do so?

Kind regards

Chrisi

[Chrisi1984]

Hello all,

related to the same topic I want to ask how should the code look like to make a query whether a cell belongs to a specified porous zone or not.

Thanks in advance and have a nice weekend!

[Chrisi1984]

Hi,

regarding my probleme it would help me to know, how I can make a query wheter the cell where the parcel is belongs to a specified cellZone or not.
I tried to do it with a loop over all cells. But that is too expansive.

forAll(mesh.cellZones(), iZone) //loop over all cellZones
{
if ((mesh.cellZones()[iZone].name() == "porous"))

{

forAll(mesh.cellZones()[iZone], iCell) //loop over all cells in Zone
{

if (iCell=cellI)

{
Unew = vector(0, 0, 0);
}


else
{
Unew = Ures.value();
}



}

}

else
{
Unew = Ures.value();
}

}

Any ideas?

Kind regards
Chrisi