zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

fvMatrix: Improved the relaxation of the matrix in the case that the diagonal is initially negative

Browse Source
This commit is contained in:
Henry
2012-02-17 12:45:10 +00:00
parent e7b828d69c
commit f71567a71d
1 changed files with 7 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -596,7 +596,7 @@ void Foam::fvMatrix<Type>::relax(const scalar alpha)
scalar sumNon = 0.0; scalar sumNon = 0.0;
forAll(D, celli) forAll(D, celli)
{ {
scalar d = (sumOff[celli] - D[celli])/D[celli]; scalar d = (sumOff[celli] - D[celli])/mag(D[celli]);
if (d > 0) if (d > 0)
{ {
@ -621,8 +621,11 @@ void Foam::fvMatrix<Type>::relax(const scalar alpha)
// Ensure the matrix is diagonally dominant... // Ensure the matrix is diagonally dominant...
// (assumes that the central coefficient is positive) // Assumes that the central coefficient is positive and ensures it is
max(D, D, sumOff); forAll(D, celli)
{
D[celli] = max(mag(D[celli]), sumOff[celli]);
}
// ... then relax // ... then relax
D /= alpha; D /= alpha;