[srinath]

Hello

I was reading Hrv's thesis and have a couple of doubts on the interpolation routine in OpenFoam

For the NVD schemes, there is an expression on page 109, which gives
\phi c \tilda = 1 - (gradphi)f.d/(2*grad phi)c.d

And a nice algorithm to give the value of phi on the face, using blended/gamma differencing.

But how do we calculate (grad phi)c?
Suppose i had rho at cell centers and wanted to obtain a limited value of rho on the surface, what would i do?
Should i use a least squared approach to find (grad phi)c?

Also i can't seem to find a formulation of TVD using slope limiters on unstructured grids.

Could someone please help me on this?

Thanks
Srinath