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

Improvement for non-orthogonality > 90 deg

Browse source
This commit is contained in:
Hrvoje Jasak 2011-03-13 23:17:19 +00:00
parent c764bdbdde
commit e290850d40
4 changed files with 5 additions and 17 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
src/finiteVolume/fvMesh/fvPatches/basic/coupled/coupledFvPatch.C
View file

@ -56,7 +56,9 @@ void coupledFvPatch::makeCorrVecs(vectorField& cv) const
vectorField patchDeltas = delta(); vectorField patchDeltas = delta();
vectorField n = nf(); 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);
} }

2
src/finiteVolume/fvMesh/fvPatches/constraint/overlapGgi/overlapGgiFvPatch.C
View file

@ -86,7 +86,7 @@ void Foam::overlapGgiFvPatch::makeCorrVecs(vectorField& cv) const
{ {
// Non-orthogonality correction on a ggi interface // Non-orthogonality correction on a ggi interface
// MB, 7/April/2009 // MB, 7/April/2009
// Calculate correction vectors on coupled patches // Calculate correction vectors on coupled patches
const scalarField& patchDeltaCoeffs = deltaCoeffs(); const scalarField& patchDeltaCoeffs = deltaCoeffs();

7
src/finiteVolume/fvMesh/fvPatches/derived/cohesive/cohesiveFvPatch.C
View file

@ -42,13 +42,6 @@ namespace Foam
void Foam::cohesiveFvPatch::makeCorrVecs(vectorField& cv) const 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; cv = vector::zero;
} }

9
src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C
View file

@ -115,7 +115,7 @@ tmp<vectorField> fvPatch::Cn() const
const unallocLabelList& faceCells = this->faceCells(); const unallocLabelList& faceCells = this->faceCells();
// get reference to global cell centres // Get reference to global cell centres
const vectorField& gcc = boundaryMesh().mesh().cellCentres(); const vectorField& gcc = boundaryMesh().mesh().cellCentres();
forAll (faceCells, faceI) forAll (faceCells, faceI)
@ -165,13 +165,6 @@ void fvPatch::makeDeltaCoeffs(scalarField& dc) const
void fvPatch::makeCorrVecs(vectorField& cv) 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; cv = vector::zero;
} }