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 clamp(PLimiter, zero_one{});
errorManipArg< error, int > exit(error &err, const int errNo=1)
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
scalar limiter(const scalar cdWeight, const scalar faceFlux, const vector &PhiP, const vector &PhiN, const vector &Sf, const scalar &) const
dimensionedScalar stabilise(const dimensionedScalar &x, const dimensionedScalar &y)
#define FatalIOErrorInFunction(ios)
Report an error message using Foam::FatalIOError.
Class with limiter function which returns the limiter for the Phi differencing scheme.
dimensionSet clamp(const dimensionSet &a, const dimensionSet &range)
IOerror FatalIOError
Error stream (stdout output on all processes), with additional 'FOAM FATAL IO ERROR' header text and ...