Go to the documentation of this file.
62 template<
class LimiterFunc>
75 const scalar cdWeight,
76 const scalar faceFlux,
77 const typename LimiterFunc::phiType& phiP,
78 const typename LimiterFunc::phiType& phiN,
79 const typename LimiterFunc::gradPhiType& gradcP,
80 const typename LimiterFunc::gradPhiType& gradcN,
84 scalar phiCD = cdWeight*phiP + (1 - cdWeight)*phiN;
91 phif = 0.5*(phiCD + phiP + (1 - cdWeight)*(d & gradcP));
96 phif = 0.5*(phiCD + phiN - cdWeight*(d & gradcN));
100 scalar QLimiter = (phif - phiU)/
stabilise(phiCD - phiU, SMALL);
103 return max(
min(QLimiter, 2), 0);
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
tmp< DimensionedField< scalar, GeoMesh > > stabilise(const DimensionedField< scalar, GeoMesh > &dsf, const dimensioned< scalar > &ds)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Class with limiter function which returns the limiter for the quadratic-upwind differencing scheme.
scalar limiter(const scalar cdWeight, const scalar faceFlux, const typename LimiterFunc::phiType &phiP, const typename LimiterFunc::phiType &phiN, const typename LimiterFunc::gradPhiType &gradcP, const typename LimiterFunc::gradPhiType &gradcN, const vector &d) const