BUG: snappyHexMesh: put internal faces only in layerFaces

This commit is contained in:
mattijs
2014-02-19 15:00:51 +00:00
parent c8446568cb
commit 1eb9ab1238

View File

@ -2463,8 +2463,8 @@ void Foam::autoLayerDriver::getLayerCellsFaces
if (layer.size()) if (layer.size())
{ {
// Layer contains both original boundary face and new boundary // Layer contains both original boundary face and new boundary
// face so is nLayers+1 // face so is nLayers+1. Leave out old internal face.
forAll(layer, i) for (label i = 1; i < layer.size(); i++)
{ {
faceRealThickness[layer[i]] = realThickness; faceRealThickness[layer[i]] = realThickness;
} }
@ -2608,7 +2608,7 @@ bool Foam::autoLayerDriver::writeLayerData
} }
{ {
label nAdded = 0; label nAdded = 0;
forAll(faceRealThickness, faceI) for (label faceI = 0; faceI < mesh.nInternalFaces(); faceI++)
{ {
if (faceRealThickness[faceI] > 0) if (faceRealThickness[faceI] > 0)
{ {
@ -2617,7 +2617,7 @@ bool Foam::autoLayerDriver::writeLayerData
} }
faceSet layerFacesSet(mesh, "layerFaces", nAdded); faceSet layerFacesSet(mesh, "layerFaces", nAdded);
forAll(faceRealThickness, faceI) for (label faceI = 0; faceI < mesh.nInternalFaces(); faceI++)
{ {
if (faceRealThickness[faceI] > 0) if (faceRealThickness[faceI] > 0)
{ {
@ -3405,7 +3405,7 @@ void Foam::autoLayerDriver::addLayers
addedCellSet.write(); addedCellSet.write();
faceSet layerFacesSet(newMesh, "layerFaces", newMesh.nFaces()/100); faceSet layerFacesSet(newMesh, "layerFaces", newMesh.nFaces()/100);
forAll(faceRealThickness, faceI) for (label faceI = 0; faceI < newMesh.nInternalFaces(); faceI++)
{ {
if (faceRealThickness[faceI] > 0) if (faceRealThickness[faceI] > 0)
{ {