exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Porting bug fix

Browse source
This commit is contained in:
Hrvoje Jasak 2015-06-10 11:43:55 +01:00
parent a8bc92b12a
commit b6b47d1017
2 changed files with 6 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatch/immersedBoundaryFvPatch.C
View file

@ -1127,7 +1127,7 @@ void Foam::immersedBoundaryFvPatch::makeIbCellCells() const
// Note: the algorithm is originally written with inward-facing normals // Note: the algorithm is originally written with inward-facing normals
// and subsequently changed: IB surface normals point outwards // and subsequently changed: IB surface normals point outwards
// HJ, 21/May/2012 // HJ, 21/May/2012
const vectorField& ibn = ibNormals(); // const vectorField& ibn = ibNormals();
forAll (cellCells, cellI) forAll (cellCells, cellI)
{ {
@ -1152,8 +1152,8 @@ void Foam::immersedBoundaryFvPatch::makeIbCellCells() const
// Collect the cells within rM of the fitting cell // Collect the cells within rM of the fitting cell
if (r <= rM[cellI]) if (r <= rM[cellI])
{ {
scalar angleLimit = // scalar angleLimit =
-Foam::cos(angleFactor_()*mathematicalConstant::pi/180); // -Foam::cos(angleFactor_()*mathematicalConstant::pi/180);
vector dir = (C[curCell] - ibp[cellI]); vector dir = (C[curCell] - ibp[cellI]);
dir /= mag(dir) + SMALL; dir /= mag(dir) + SMALL;

5
src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchField/immersedBoundaryFvPatchField.C
View file

@ -39,7 +39,8 @@ template<class Type>
const Foam::debug::optimisationSwitch const Foam::debug::optimisationSwitch
immersedBoundaryFvPatchField<Type>::nBcIter_ immersedBoundaryFvPatchField<Type>::nBcIter_
( (
debug::optimisationSwitch("immersedBoundaryNBCIter", 5) "immersedBoundaryNBCIter",
5
); );
@ -251,7 +252,7 @@ immersedBoundaryFvPatchField<Type>::imposeDirichletCondition() const
reduce(maxError, maxOp<scalar>()); reduce(maxError, maxOp<scalar>());
} }
while (maxError > bcTolerance_() && counter < nBcIter_()); while (maxError > bcTolerance_ && counter < nBcIter_);
if (counter == nBcIter_() && debug) if (counter == nBcIter_() && debug)
{ {