[anandsudhi]

Hi,

the improved version of limited interpolation Schemes for vectors like limitedLinearV, linearUpwindV are different from the the original versions (linearUpwind and limitedLinear).

If anyone can tell me the extra modification done to linearUpwindV to account for direction of the field. I understand there is a new variable 'maxCorr' in the code. But i am not able to understand its implementation.

Also does this modification make linearUpwindV bounded with out the need of specifying limited gradients?

Any comments will be very much appreciated.

Thank you

Anand

[anandsudhi]

forAll(faceFlux, facei)
{


vector maxCorr;

if (faceFlux[facei] > 0.0)
{
maxCorr =
(1.0 - w[facei])*(vf[nei[facei]] - vf[own[facei]]);

sfCorr[facei] =
(Cf[facei] - C[own[facei]]) & gradVf[own[facei]];
}
else
{
maxCorr =
w[facei]*(vf[own[facei]] - vf[nei[facei]]);

sfCorr[facei] =
(Cf[facei] - C[nei[facei]]) & gradVf[nei[facei]];
}

scalar sfCorrs = magSqr(sfCorr[facei]);
scalar maxCorrs = sfCorr[facei] & maxCorr;

if (sfCorrs > 0)
{
if (maxCorrs < 0)
{
sfCorr[facei] = vector::zero;
}
else if (sfCorrs > maxCorrs)
{
sfCorr[facei] *= maxCorrs/(sfCorrs + VSMALL);
}
}
else if (sfCorrs < 0)
{
if (maxCorrs > 0)
{
sfCorr[facei] = vector::zero;
}
else if (sfCorrs < maxCorrs)
{
sfCorr[facei] *= maxCorrs/(sfCorrs - VSMALL);
}
}
}