47 "linearUpwindCorrection(" + vf.
name() +
')',
71 gradVf = gradScheme_().grad(vf);
75 if (faceFlux[facei] > 0)
77 label own = owner[facei];
78 sfCorr[facei] = (Cf[facei] - C[own]) & gradVf[own];
82 label nei = neighbour[facei];
83 sfCorr[facei] = (Cf[facei] - C[nei]) & gradVf[nei];
89 GeometricBoundaryField& bSfCorr = sfCorr.boundaryField();
105 gradVf.boundaryField()[
patchi].patchNeighbourField();
123 label own = pOwner[facei];
125 if (pFaceFlux[facei] > 0)
127 pSfCorr[facei] = (pCf[facei] - C[own]) & gradVf[own];
132 (pCf[facei] - pd[facei] - C[own]) & pGradVfNei[facei];