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

Least square vectors next to fixed gradient boundaries work better with short distance

Browse source
This commit is contained in:
Hrvoje Jasak 2010-09-30 23:01:55 +01:00
parent db4a63ed07
commit f675358ffd
1 changed files with 17 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

18
src/finiteVolume/finiteVolume/gradSchemes/leastSquaresGrad/leastSquaresVectors.C
View file

@ -135,8 +135,24 @@ void Foam::leastSquaresVectors::makeLeastSquaresVectors() const
const unallocLabelList& faceCells = p.patch().faceCells();
// Build the d-vectors
// Original version: closest distance to boundary
vectorField pd =
mesh_.Sf().boundaryField()[patchi]
/(
mesh_.magSf().boundaryField()[patchi]
*mesh_.deltaCoeffs().boundaryField()[patchi]
);
if (!mesh_.orthogonal())
{
pd -= mesh_.correctionVectors().boundaryField()[patchi]
/mesh_.deltaCoeffs().boundaryField()[patchi];
}
// Better version of d-vectors: Zeljko Tukovic, 25/Apr/2010
vectorField pd = p.delta();
// Experimental: review fixed gradient condition. HJ, 30/Sep/2010
// vectorField pd = p.delta();
if (p.coupled())
{