vanalbada limiter
 

Hi Foamers.

Does anybody know why the vanAlbada limiter (src/finiteVolume / interpolation / surfaceInterpolation / limitedSchemes / vanAlbada/vanAlbada.H)

isimplemented as :

r*(r + 1)/(sqr(r) + 1)

Where the van albada limiter is actually,
r*(r + 1)/(r**2 + 1) (see for example wikipedia or Berger and Aftosmis,
nalysis of Slope Limiters on Irregular Grids).

Another troubling aspect is that the OF implementation does not
fulfill Sweby conditions for limiters.

The limiter is greater than 2 when r is greater then ~ 1.7.

Hi,

sqr(r) == r**2

(square root is sqrt in C)

oh my mistake I thought it is square root...