65 const scalar faceFlux,
73 scalar gradf = phiN - phiP;
86 if (
mag(gradf) >= 1000*
mag(gradcf))
88 return 1 - 0.5*1000*
sign(gradcf)*
sign(gradf);
92 return 1 - 0.5*gradf/gradcf;
99 const scalar faceFlux,
107 scalar gradf = phiN - phiP;
120 if (
mag(gradcf) >= 1000*
mag(gradf))
122 return 2*1000*
sign(gradcf)*
sign(gradf) - 1;
126 return 2*(gradcf/gradf) - 1;
dimensionedScalar sign(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
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
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...