zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: wallBoundedStreamline: handle patches with zero faces. See #502.

Browse Source
This commit is contained in:
mattijs
2017-06-21 14:05:26 +01:00
parent cbbf08dc0d
commit c462d48196
2 changed files with 34 additions and 30 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/functionObjects/field/streamLine/streamLineBase.C
View File

@ -767,7 +767,7 @@ bool Foam::functionObjects::streamLineBase::write()
// Convert scalar values
if (allScalars_.size() > 0)
if (allScalars_.size() > 0 && tracks.size() > 0)
{
List<List<scalarField>> scalarValues(allScalars_.size());
@ -811,7 +811,7 @@ bool Foam::functionObjects::streamLineBase::write()
// Convert vector values
if (allVectors_.size() > 0)
if (allVectors_.size() > 0 && tracks.size() > 0)
{
List<List<vectorField>> vectorValues(allVectors_.size());

60
src/functionObjects/field/wallBoundedStreamLine/wallBoundedStreamLine.C
View File

@ -139,39 +139,43 @@ void Foam::functionObjects::wallBoundedStreamLine::track()
forAll(seedPoints, i)
{
const point& seedPt = seedPoints[i];
label celli = seedPoints.cells()[i];
tetIndices ids(findNearestTet(isWallPatch, seedPt, celli));
if (ids.face() != -1 && isWallPatch[ids.face()])
if (celli != -1)
{
//Pout<< "Seeding particle :" << nl
// << " seedPt:" << seedPt << nl
// << " face :" << ids.face() << nl
// << " at :" << mesh_.faceCentres()[ids.face()] << nl
// << " cell :" << mesh_.cellCentres()[ids.cell()] << nl
// << endl;
const point& seedPt = seedPoints[i];
tetIndices ids(findNearestTet(isWallPatch, seedPt, celli));
particles.addParticle
(
new wallBoundedStreamLineParticle
if (ids.face() != -1 && isWallPatch[ids.face()])
{
//Pout<< "Seeding particle :" << nl
// << " seedPt:" << seedPt << nl
// << " face :" << ids.face() << nl
// << " at :" << mesh_.faceCentres()[ids.face()]
// << nl
// << " cell :" << mesh_.cellCentres()[ids.cell()]
// << nl << endl;
particles.addParticle
(
mesh_,
ids.faceTri(mesh_).centre(),
ids.cell(),
ids.face(), // tetFace
ids.tetPt(),
-1, // not on a mesh edge
-1, // not on a diagonal edge
lifeTime_ // lifetime
)
);
}
else
{
Pout<< type() << " : ignoring seed " << seedPt
<< " since not in wall cell." << endl;
new wallBoundedStreamLineParticle
(
mesh_,
ids.faceTri(mesh_).centre(),
ids.cell(),
ids.face(), // tetFace
ids.tetPt(),
-1, // not on a mesh edge
-1, // not on a diagonal edge
lifeTime_ // lifetime
)
);
}
else
{
Pout<< type() << " : ignoring seed " << seedPt
<< " since not in wall cell." << endl;
}
}
}
}