stabilisation for zero-area faces

This commit is contained in:
mattijs
2009-08-10 11:30:52 +01:00
parent 5ac6840d06
commit c635b44ae5

View File

@ -626,7 +626,8 @@ bool Foam::primitiveMesh::checkFaceSkewness
vector d = cellCtrs[nei[faceI]] - cellCtrs[own[faceI]];
// Skewness vector
vector sv = Cpf - ((fAreas[faceI] & Cpf)/(fAreas[faceI] & d))*d;
vector sv =
Cpf - ((fAreas[faceI] & Cpf)/((fAreas[faceI] & d) + SMALL))*d;
vector svHat = sv/(mag(sv) + VSMALL);
// Normalisation distance calculated as the approximate distance