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

Tight tolerances on mesh motion error

Browse source
This commit is contained in:
Hrvoje Jasak 2019-03-27 18:12:22 +00:00
parent 62228e7189
commit 6c2263a2bf
1 changed files with 4 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -135,7 +135,7 @@ void Foam::immersedBoundaryFvPatch::updatePhi
scalarField& phiIn = phi.internalField(); scalarField& phiIn = phi.internalField();
const labelList& deadFaces = ibPolyPatch_.deadFaces(); const labelList& deadFaces = ibPolyPatch_.deadFaces();
forAll(deadFaces, dfI) forAll (deadFaces, dfI)
{ {
const label faceI = deadFaces[dfI]; const label faceI = deadFaces[dfI];
if (mesh.isInternalFace(faceI)) if (mesh.isInternalFace(faceI))
@ -145,7 +145,7 @@ void Foam::immersedBoundaryFvPatch::updatePhi
} }
const labelList& cutFaces = ibPolyPatch_.ibFaces(); const labelList& cutFaces = ibPolyPatch_.ibFaces();
forAll(cutFaces, cfI) forAll (cutFaces, cfI)
{ {
const label faceI = cutFaces[cfI]; const label faceI = cutFaces[cfI];
if (mesh.isInternalFace(faceI)) if (mesh.isInternalFace(faceI))
@ -225,7 +225,8 @@ void Foam::immersedBoundaryFvPatch::updatePhi
// HJ, 22/Dec/2017 // HJ, 22/Dec/2017
forAll (magDivPhi, cellI) forAll (magDivPhi, cellI)
{ {
if (magDivPhi[cellI] > SMALL) // if (magDivPhi[cellI] > SMALL)
if (magDivPhi[cellI] > 1e-40)
{ {
// Attempt to correct via old volume // Attempt to correct via old volume
scalar corrOldVol = newVols[cellI] - divPhi[cellI]*deltaT; scalar corrOldVol = newVols[cellI] - divPhi[cellI]*deltaT;