Improvement for non-orthogonality > 90 deg

src/finiteVolume/fvMesh/fvPatches/basic/coupled/coupledFvPatch.C

@@ -56,7 +56,9 @@ void coupledFvPatch::makeCorrVecs(vectorField& cv) const
     vectorField patchDeltas = delta();
     vectorField n = nf();
 
-    cv = n - patchDeltas*patchDeltaCoeffs;
+    // If non-orthogonality is over 90 deg, kill correction vector
+    // HJ, 27/Feb/2011
+    cv = pos(n & patchDeltas)*(n - patchDeltas*patchDeltaCoeffs);
 }
 
 

src/finiteVolume/fvMesh/fvPatches/derived/cohesive/cohesiveFvPatch.C

@@ -42,13 +42,6 @@ namespace Foam
 
 void Foam::cohesiveFvPatch::makeCorrVecs(vectorField& cv) const
 {
-    // Correction vector
-//     vectorField d = delta();
-//     vectorField n = nf();
-//     cv = n - d/(n & d);
-
-    // Old version: no non-orthogonal correction on the wall
-    // Testing.  HJ, 9/Apr/2009
     cv = vector::zero;
 }
 

src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C

@@ -115,7 +115,7 @@ tmp<vectorField> fvPatch::Cn() const
 
     const unallocLabelList& faceCells = this->faceCells();
 
-    // get reference to global cell centres
+    // Get reference to global cell centres
     const vectorField& gcc = boundaryMesh().mesh().cellCentres();
 
     forAll (faceCells, faceI)
@@ -165,13 +165,6 @@ void fvPatch::makeDeltaCoeffs(scalarField& dc) const
 
 void fvPatch::makeCorrVecs(vectorField& cv) const
 {
-    // Correction vector
-//     vectorField d = delta();
-//     vectorField n = nf();
-//     cv = n - d/(n & d);
-
-    // Old version: no non-orthogonal correction on the wall
-    // Testing.  HJ, 9/Apr/2009
     cv = vector::zero;
 }