From b6b47d10176fd16eb248adfdff393de4c52bacbf Mon Sep 17 00:00:00 2001 From: Hrvoje Jasak Date: Wed, 10 Jun 2015 11:43:55 +0100 Subject: [PATCH] Porting bug fix --- .../immersedBoundaryFvPatch/immersedBoundaryFvPatch.C | 6 +++--- .../immersedBoundaryFvPatchField.C | 5 +++-- 2 files changed, 6 insertions(+), 5 deletions(-) diff --git a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatch/immersedBoundaryFvPatch.C b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatch/immersedBoundaryFvPatch.C index df0d5aee1..1a19cc8c2 100644 --- a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatch/immersedBoundaryFvPatch.C +++ b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatch/immersedBoundaryFvPatch.C @@ -1127,7 +1127,7 @@ void Foam::immersedBoundaryFvPatch::makeIbCellCells() const // Note: the algorithm is originally written with inward-facing normals // and subsequently changed: IB surface normals point outwards // HJ, 21/May/2012 - const vectorField& ibn = ibNormals(); +// const vectorField& ibn = ibNormals(); forAll (cellCells, cellI) { @@ -1152,8 +1152,8 @@ void Foam::immersedBoundaryFvPatch::makeIbCellCells() const // Collect the cells within rM of the fitting cell if (r <= rM[cellI]) { - scalar angleLimit = - -Foam::cos(angleFactor_()*mathematicalConstant::pi/180); +// scalar angleLimit = +// -Foam::cos(angleFactor_()*mathematicalConstant::pi/180); vector dir = (C[curCell] - ibp[cellI]); dir /= mag(dir) + SMALL; diff --git a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchField/immersedBoundaryFvPatchField.C b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchField/immersedBoundaryFvPatchField.C index 432b4ccc5..33778c5dc 100644 --- a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchField/immersedBoundaryFvPatchField.C +++ b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchField/immersedBoundaryFvPatchField.C @@ -39,7 +39,8 @@ template const Foam::debug::optimisationSwitch immersedBoundaryFvPatchField::nBcIter_ ( - debug::optimisationSwitch("immersedBoundaryNBCIter", 5) + "immersedBoundaryNBCIter", + 5 ); @@ -251,7 +252,7 @@ immersedBoundaryFvPatchField::imposeDirichletCondition() const reduce(maxError, maxOp()); } - while (maxError > bcTolerance_() && counter < nBcIter_()); + while (maxError > bcTolerance_ && counter < nBcIter_); if (counter == nBcIter_() && debug) {