Commit 731ebc91 authored by Andrew Heather's avatar Andrew Heather
Browse files

ENH: Particle interaction lists - do not include partially open wall faces. See #957

parent 0f4c6572
......@@ -324,10 +324,15 @@ void Foam::InteractionLists<ParticleType>::buildInteractionLists()
if (isA<wallPolyPatch>(patch))
{
localWallFaces.append
(
identity(patch.size(), patch.start())
);
const scalarField areaFraction(patch.areaFraction());
forAll(areaFraction, facei)
{
if (areaFraction[facei] > 0.5)
{
localWallFaces.append(facei + patch.start());
}
}
}
}
......@@ -335,10 +340,12 @@ void Foam::InteractionLists<ParticleType>::buildInteractionLists()
forAll(wallFaceBbs, i)
{
wallFaceBbs[i] = treeBoundBox
(
mesh_.faces()[localWallFaces[i]].points(mesh_.points())
);
wallFaceBbs[i] =
treeBoundBox
(
mesh_.points(),
mesh_.faces()[localWallFaces[i]]
);
}
// IAndT: index and transform
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment