Implementation of new convective discretization scheme
 

Hello everyone,

I would like to implement a new discretization scheme for the convective term. For each cell of my domain, I have the corresponding coefficients for the north, south, west and east cells gathered on a file (.txt, .csv or in OpenFOAM format as a new internalField).

I would like to construct the matrix of coefficients using these custom convective coefficients rather than using the tradicional Upwind or QUICK coefficients.

Anyone could give me some guidelines or some examples to implement this modification? Thank you in advance.

For theory: https://link.springer.com/book/10.10...-3-319-16874-6

For implementation/coding: https://github.com/UnnamedMoose/Basi...OFtutorial15.H