[ckpark]

Dear Foamers,


I'm trying to simulate a flow of Re~38k in around a triangular cylinder in a narrow channel using k-omega SST.I am using the PIMPLE solver in OpenFoam for this simulation. I do have some experience with the program but I am constantly facing an issue in my simulation results.


At the early phase of the simulation when the flow has yet to stabilise, the high amount of TKE behind the cylinder is pretty well presented.
However, when the flow starts to form smaller vortices and way more rapid separations are observed in the velocity field (as expected), the TKE field shows a weird behaviour.


The flow right behind the cylinder has a tiny amount of TKE (almost none) but as it proceeds downstream, the magnitude of TKE grows.


I have tried to change many variables in the OpenFoam settings such as boundary/initial conditions, various numerical methods in fvSchemes and different solver/relaxation factor/pimple settings in fvSolution but the results are either it diverges and crash or show the exact same behaviour.


Could someone advise me on a solution or different approach I should take to tackle this problem?
I have attached the screenshots of velocity and TKE field together with the 0th time variables, fvSchemes, fvSolution and controlDict



Thank you