BUG: snappyHexMesh: incorrect feature angle. Fixes #699.

This commit is contained in:
mattijs
2018-01-08 11:50:33 +00:00
parent 34a71c19dc
commit f868b11a31
2 changed files with 7 additions and 4 deletions

View File

@ -668,7 +668,7 @@ void Foam::snappyLayerDriver::handleFeatureAngle
( (
const indirectPrimitivePatch& pp, const indirectPrimitivePatch& pp,
const labelList& meshEdges, const labelList& meshEdges,
const scalar minCos, const scalar minAngle,
pointField& patchDisp, pointField& patchDisp,
labelList& patchNLayers, labelList& patchNLayers,
List<extrudeMode>& extrudeStatus List<extrudeMode>& extrudeStatus
@ -676,7 +676,10 @@ void Foam::snappyLayerDriver::handleFeatureAngle
{ {
const fvMesh& mesh = meshRefiner_.mesh(); const fvMesh& mesh = meshRefiner_.mesh();
Info<< nl << "Handling feature edges ..." << endl; const scalar minCos = Foam::cos(degToRad(minAngle));
Info<< nl << "Handling feature edges (angle < " << minAngle
<< ") ..." << endl;
if (minCos < 1-SMALL) if (minCos < 1-SMALL)
{ {
@ -3647,7 +3650,7 @@ void Foam::snappyLayerDriver::addLayers
( (
pp, pp,
meshEdges, meshEdges,
degToRad(layerParams.featureAngle()), layerParams.featureAngle(),
patchDisp, patchDisp,
patchNLayers, patchNLayers,

View File

@ -218,7 +218,7 @@ private:
( (
const indirectPrimitivePatch& pp, const indirectPrimitivePatch& pp,
const labelList& meshEdges, const labelList& meshEdges,
const scalar minCos, const scalar minAngle,
pointField& patchDisp, pointField& patchDisp,
labelList& patchNLayers, labelList& patchNLayers,
List<extrudeMode>& extrudeStatus List<extrudeMode>& extrudeStatus