自编二维的矩形管道流动的问题探讨

xiaojun218

问题是一个二维的矩形管道流动，在x方向长度为2，y方向宽度为0.2，划分的网格分别为400和60，所划分的网格均为标准的四方形，其中左侧为入口，右侧为出口，入口速度为4米每秒。这个问题是用来检验我所编译的一个SST K-omega双方程模型的计算程序的一个算例。
    整个问题的假设是：通过二维泊肃叶流动的解析解给定整个内部的流场流速并保持不变，内部压力也不变（为大气压）
    目的：通过上面的假设来得到我期望的k和omega分布

    详细的内容都已经写在了后面的附件内，由于都是些公式，无法直接黏贴到上面来，见谅
    我的第一个问题是我这样的假设合理吗？至少到目前为止我还是没有得到满意的结果，最终总是出错了。当中的原因自己也不是很清楚。不过还是先提一下自己的一些想法，在调试过程中出错最多的便是对流项，而且入口顶部靠近入口和壁面的网格的omega值同时受到来流和壁面的影响，和后面的临近壁面的网格的omega值相差很大，在采用了附件中列举的对流项处理方法后容易造成对流项较大从而引起出错。当然，这只是我自己的看法，还要请教高手来最终的判断。

在附件中，包括空间和时间的离散、对流项扩散项的处理都已经列举出来了，希望高手能帮忙看一下，并提出修改意见，谢谢了。

maok

刚接触CFD几个月。一看楼主的附件，才知道什么叫做高手。楼主博士？

shirazbj

你的假设是整个内部的流场流速不变，把k和e输运过来. 我觉得做一步还行,时间一步进可能就不对了.k和e变了，速度也应该变呀。