[jgs55]

Hi everyone,

I want to implement a custom surface interpolation scheme for the convective term, so that:

div(phi,T)      Gauss myCustomScheme;

For each face (both internal or boundary) I have an array of 6 coefficients representing the weights () of my custom scheme. With them, I obtain the value of T at the face as a linear combination of 6 neighbouring cell center values:

, being: the transport property, the face to get the interpolated value, the corresponding weight and each neighboring cell.

These 6 cells correspond to the E, W, NE, NW, SE, and SW neighbouring cells of the face . Besides, these coefficients are time dependent, so I will have new coefficients every timestep.

I am using the code of the upwing scheme as template for my implementation, however, the upwing scheme only uses previous or following cell (E or W) instead of the 6 neighbours. I also tried the QUICK scheme as template, but I did not get any good result.

Could anyone help me on how i should modify a existing scheme to fullfil my requirements? Thank you very much in advance.