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

ENH: limitedSnGrad: added debug writing

Browse Source
This commit is contained in:
mattijs
2019-03-25 15:42:29 +00:00
committed by Andrew Heather
parent 33894bf893
commit 1561a46927
1 changed files with 35 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

35
src/finiteVolume/finiteVolume/snGradSchemes/limitedSnGrad/limitedSnGrad.C
View File

@ -30,6 +30,7 @@ License
#include "volFields.H" #include "volFields.H"
#include "surfaceFields.H" #include "surfaceFields.H"
#include "localMax.H" #include "localMax.H"
#include "fvcCellReduce.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -82,6 +83,40 @@ limitedSnGrad<Type>::correction
<< "limiter min: " << min(limiter.primitiveField()) << "limiter min: " << min(limiter.primitiveField())
<< " max: "<< max(limiter.primitiveField()) << " max: "<< max(limiter.primitiveField())
<< " avg: " << average(limiter.primitiveField()) << endl; << " avg: " << average(limiter.primitiveField()) << endl;
if (fv::debug & 2)
{
static scalar oldTime = -1;
static label subIter = 0;
if (vf.mesh().time().value() != oldTime)
{
oldTime = vf.mesh().time().value();
subIter = 0;
}
else
{
++subIter;
}
word fieldName("limiter_" + Foam::name(subIter));
GeometricField<scalar, fvPatchField, volMesh> volLimiter
(
IOobject
(
fieldName,
vf.mesh().time().timeName(),
vf.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE,
false
),
fvc::cellReduce(limiter, minEqOp<scalar>(), scalar(1.0))
);
Info<< "Writing limiter field to " << volLimiter.objectPath()
<< endl;
volLimiter.write();
}
} }
return limiter*corr; return limiter*corr;