songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

fvMatrix: Further improvement to relaxation procedure

Browse Source
This commit is contained in:
Henry 2012-05-30 15:23:33 +01:00
parent ea6a8d4b6b
commit d5b8f549b9
1 changed files with 10 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
src/finiteVolume/fvMatrices/fvMatrix/fvMatrix.C
View File

@ -570,24 +570,11 @@ void Foam::fvMatrix<Type>::relax(const scalar alpha)
}
else
{
// For non-coupled boundaries subtract the diagonal
// contribution off-diagonal sum which avoids having to remove
// it from the diagonal later.
// Also add the source contribution from the relaxation
// For non-coupled boundaries add the maximum magnitude diagonal
// contribution to ensure stability
forAll(pa, face)
{
// Type iCoeff0 = iCoeffs[face];
// iCoeffs[face] = cmptMag(iCoeffs[face]);
// sumOff[pa[face]] -= cmptMin(iCoeffs[face]);
// iCoeffs[face] /= alpha;
D[pa[face]] +=
cmptMag(cmptMin(iCoeffs[face]))
- cmptMin(iCoeffs[face]);
sumOff[pa[face]] +=
cmptMag(cmptMin(iCoeffs[face]))
- cmptMin(iCoeffs[face]);
// S[pa[face]] +=
// cmptMultiply(iCoeffs[face] - iCoeff0, psi_[pa[face]]);
D[pa[face]] += cmptMax(cmptMag(iCoeffs[face]));
}
}
}
@ -653,6 +640,13 @@ void Foam::fvMatrix<Type>::relax(const scalar alpha)
D[pa[face]] -= component(iCoeffs[face], 0);
}
}
else
{
forAll(pa, face)
{
D[pa[face]] -= cmptMin(iCoeffs[face]);
}
}
}
}