[rxgrch]

Hello,

I want to implement a high-order difference sheme which is 4 or 5 order, but it seemed that the difference order is just up to 3 order, the difference template in OpenFOAM now is just two cells, but the high order（4or5） will need more than 3 neighbor cells, so who can tell me how to implement the high order difference with the smallest modify with the current C++ code?

[mkraposhin]

Hello, i tried to implement high-order schemes in OpenFOAM (WENO + k-exact) one year ago, but i found, that i don't know how to select proper stencil for scheme. If You will help with mathematics, then i can help you with implementation

[rxgrch]

Originally Posted by mkraposhin

Hello, i tried to implement high-order schemes in OpenFOAM (WENO + k-exact) one year ago, but i found, that i don't know how to select proper stencil for scheme. If You will help with mathematics, then i can help you with implementation

Thank you for your reply, and I don't quite get what do you mean " i don't know how to select proper stencil for scheme"? the scheme that we want to implement is WCNS scheme. So can you give me some advice or with a example code of high-order sheme? Thanks very much!

[rxgrch]

Originally Posted by mkraposhin

Hello, i tried to implement high-order schemes in OpenFOAM (WENO + k-exact) one year ago, but i found, that i don't know how to select proper stencil for scheme. If You will help with mathematics, then i can help you with implementation

what you said had been solved in a paper "High Order Schemes on Three-Dimensional General Polyhedral Meshes—Application to the Level Set Method", it maybe help you.