〖求助〗圆柱绕流的一个问题？

yuarmy

我用c语言编写了一个计算稳态圆柱绕流问题的程序，离散方法是由于计算的不可压缩流体所以，采用simple算法对压力和速度进行解耦，对流项是采用中心差分和迎风差分结合的延迟修正（只有一阶精度），就是对peric的程序进行了一些修改。没有采用贴体网格，只是用结构化网格进行处理，边界处理，就是将圆柱的速度条件直接移动到临近网格点形成的。但是，完成程序后，发现只能计算雷诺数在10以下的情况，当雷诺数再高的时候，计算就发散了。并且在雷诺数为10时，可以看到，流体绕圆柱流过，但是并没有见到，圆柱后面的涡流，请问各位大侠，下一步我应该做哪些工作来修改我的程序，使其能计算更高雷诺数的流体，并且计算出圆柱后面的涡流来。谢谢。