57 template<
class LimiterFunc>
70 const scalar cdWeight,
71 const scalar faceFlux,
72 const typename LimiterFunc::phiType& phiP,
73 const typename LimiterFunc::phiType& phiN,
74 const typename LimiterFunc::gradPhiType& gradcP,
75 const typename LimiterFunc::gradPhiType& gradcN,
79 scalar phiCD = cdWeight*phiP + (1 - cdWeight)*phiN;
86 phif = 0.5*(phiCD + phiP + (1 - cdWeight)*(d & gradcP));
91 phif = 0.5*(phiCD + phiN - cdWeight*(d & gradcN));
95 scalar QLimiter = (phif -
phiU)/
stabilise(phiCD - phiU, SMALL);
98 return max(
min(QLimiter, 2), 0);