75 gradVf = gradScheme_().grad(vf);
81 if (faceFlux[facei] > 0.0)
85 *(vf[nei[facei]] - vf[own[facei]]);
88 (Cf[facei] - C[own[facei]]) & gradVf[own[facei]];
93 w[facei]*(vf[own[facei]] - vf[nei[facei]]);
96 (Cf[facei] - C[nei[facei]]) & gradVf[nei[facei]];
99 scalar sfCorrs =
magSqr(sfCorr[facei]);
100 scalar maxCorrs = sfCorr[facei] & maxCorr;
106 sfCorr[facei] = vector::zero;
108 else if (sfCorrs > maxCorrs)
110 sfCorr[facei] *= maxCorrs/(sfCorrs + VSMALL);
113 else if (sfCorrs < 0)
117 sfCorr[facei] = vector::zero;
119 else if (sfCorrs < maxCorrs)
121 sfCorr[facei] *= maxCorrs/(sfCorrs - VSMALL);