70 <<
"coefficient = " << k_
71 <<
" should be >= 0 and <= 1"
78 const scalar cdWeight,
79 const scalar faceFlux,
86 scalar phiP = Sf&PhiP;
87 scalar phiN = Sf&PhiN;
100 scalar phiCD = cdWeight*phiP + (1 - cdWeight)*phiN;
107 ((faceFlux - phiU)/
stabilise(phiCD - phiU, SMALL) + k_);
110 return max(
min(PLimiter, 1), 0);
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
Class with limiter function which returns the limiter for the Phi differencing scheme.
scalar limiter(const scalar cdWeight, const scalar faceFlux, const vector &PhiP, const vector &PhiN, const vector &Sf, const scalar &) const
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensionedScalar stabilise(const dimensionedScalar &x, const dimensionedScalar &y)
errorManipArg< error, int > exit(error &err, const int errNo=1)