WENO-OF / WENOEXT

Weighted essentially non-oscillatory library for the framework of OpenFOAM
GNU General Public License v3.0
119 stars 64 forks source link

Does WENOEXT work on 1d advection test case? #53

Open wangsen992 opened 1 year ago

wangsen992 commented 1 year ago

Hi there,

I am trying to test the package by trying it on a 1D test case running with scalarTransportFoam. When using the package directly by replacing the divScheme for div(phi, T) to "Gauss WENOUpwindFit 4 0", or with any other order number, the result is extremely diffusively, behaving like a basic upwind scheme.

After checking the constant/WENOBase* folder, particularly the file "StencilIDs", I realize that for most elements in the file, it is empty, image

whereas for your tutorials, it looks like this. image

Therefore, I am wondering if there is some limit/check in the implementation that prevents the scheme from working on 1D cases, therefore it simply behaves as an upwind scheme (as in no explicit correction is applied)?

Thank you very much and looking forward to hearing your replies!

Best wishes, Sen