63 const scalar faceFlux,
71 scalar gradf = phiN - phiP;
84 if (
mag(gradf) >= 1000*
mag(gradcf))
86 return 1 - 0.5*1000*
sign(gradcf)*
sign(gradf);
90 return 1 - 0.5*gradf/gradcf;
97 const scalar faceFlux,
105 scalar gradf = phiN - phiP;
118 if (
mag(gradcf) >= 1000*
mag(gradf))
120 return 2*1000*
sign(gradcf)*
sign(gradf) - 1;
124 return 2*(gradcf/gradf) - 1;
scalar phict(const scalar faceFlux, const scalar phiP, const scalar phiN, const vector &gradcP, const vector &gradcN, const vector &d) const
scalar r(const scalar faceFlux, const scalar phiP, const scalar phiN, const vector &gradcP, const vector &gradcN, const vector &d) const
dimensionedScalar sign(const dimensionedScalar &ds)
dimensioned< scalar > mag(const dimensioned< Type > &)