41 #ifndef limitedCubicV_H
42 #define limitedCubicV_H
55 template<
class LimiterFunc>
72 <<
"coefficient = " << k_
73 <<
" should be >= 0 and <= 1"
78 twoByk_ = 2.0/
max(k_, SMALL);
83 const scalar cdWeight,
84 const scalar faceFlux,
85 const typename LimiterFunc::phiType& phiP,
86 const typename LimiterFunc::phiType& phiN,
87 const typename LimiterFunc::gradPhiType& gradcP,
88 const typename LimiterFunc::gradPhiType& gradcN,
92 scalar twor = twoByk_*LimiterFunc::r
94 faceFlux, phiP, phiN, gradcP, gradcN, d
97 vector fV = cdWeight*phiP + (1.0 - cdWeight)*phiN;
99 scalar fVphiP = fV & phiP;
100 scalar fVphiN = fV & phiN;
115 cdWeight*(fVphiP - 0.25*(fV & (d & gradcN)))
116 + (1 - cdWeight)*(fVphiN + 0.25*(fV & (d & gradcP)));
118 scalar fVphiCD = cdWeight*fVphiP + (1 - cdWeight)*fVphiN;
121 scalar cubicLimiter =
122 (fVphif - fVphiU)/
stabilise(fVphiCD - fVphiU, SMALL);
125 return max(
min(
min(twor, cubicLimiter), 2), 0);