50 "linearUpwindV::correction(" + vf.name() +
')',
51 mesh.time().timeName(),
81 > tgradVf = gradScheme_().grad(vf, gradSchemeName_);
88 >& gradVf = tgradVf();
94 if (faceFlux[facei] > 0.0)
97 (1.0 - w[facei])*(vf[nei[facei]] - vf[own[facei]]);
100 (Cf[facei] -
C[own[facei]]) & gradVf[own[facei]];
105 w[facei]*(vf[own[facei]] - vf[nei[facei]]);
108 (Cf[facei] -
C[nei[facei]]) & gradVf[nei[facei]];
111 scalar sfCorrs =
magSqr(sfCorr[facei]);
112 scalar maxCorrs = sfCorr[facei] & maxCorr;
118 sfCorr[facei] =
Zero;
120 else if (sfCorrs > maxCorrs)
122 sfCorr[facei] *= maxCorrs/(sfCorrs + VSMALL);
129 Boundary& bSfCorr = sfCorr.boundaryFieldRef();
138 mesh.boundary()[patchi].faceCells();
147 gradVf.boundaryField()[patchi].patchNeighbourField()
160 label own = pOwner[facei];
164 if (pFaceFlux[facei] > 0)
166 pSfCorr[facei] = (pCf[facei] -
C[own]) & gradVf[own];
168 maxCorr = (1.0 - pW[facei])*(pVfNei[facei] - vf[own]);
173 (pCf[facei] - pd[facei] -
C[own]) & pGradVfNei[facei];
175 maxCorr = pW[facei]*(vf[own] - pVfNei[facei]);
178 scalar pSfCorrs =
magSqr(pSfCorr[facei]);
179 scalar maxCorrs = pSfCorr[facei] & maxCorr;
185 pSfCorr[facei] =
Zero;
187 else if (pSfCorrs > maxCorrs)
189 pSfCorr[facei] *= maxCorrs/(pSfCorrs + VSMALL);