请问Simple法模拟中出现此问题可能是什么造成的？

lyzhang98

我模拟上表面由驱使速度的流动，长相对宽长的多，粘度很小。结果：
1，前一段，只有上面第一个内节点的速度为正（边界条件），下面的速度均是负的。改变竖直方向的网格数到原来的3倍，还是这样。我想可能不再网格的划分上
2，我使用的是NS方程组。粘度很小，当把粘度变大（10三次方）时，上面的情况不存在了。这说明模型还有一定的正确性

jackliucheng

你好象没有提及收敛条件,如△t应该取多少.

lyzhang98

收敛判据：方程残差相对值小于某一值（10的-3次方）。这里模拟的是定常流

maples1

Reynolds number changes a lot. There is no fault with your model.

lyzhang98

谢谢maplesl.但是，我的结果只有上面第一个内节点的速度为正（边界条件），下面的速度均是负的。改变竖直方向的网格数到原来的3倍，还是上面第一个节点速度为正的，请帮忙分析一下，谢谢！！

maples1

Can you tell us how much your RE is?

lyzhang98

of course.我研究的Re计算方法：utop*H/v最后得：6400是金属融液

maples1

Can you change the Re number to 100 or something like that?
If the problem is still there, there could be something wrong with your code.
Otherwise, it may physically true.

lyzhang98

我将粘度变大，可以克服这种情况。但是，不好意思，我想问一下，为什么网格变细3倍，还是只有一个点的速度是正的呢？按道理讲，第一各节点是正的，当我加倍网格后，应该前一个和它两个节点是正的呀。我知道我可能是在想当然，但还是真心希望能帮我解释一下，不胜感激·

maples1

How many grid points do you have in y-direction?
网格变细3倍 does not tell this information.
BTW, the laminar assumption may not hold for Re = 6400.

lyzhang98

我还想问，网格变细3倍的话，模型没错，那么原来正的地方和这个点的上面都应该使正的呀！这种想法错在哪里呢？
2。6400不算很大呀，是金属液体，刚刚呈现湍流，有前人研究属于刚刚呈现湍流的状态，因为加上k-e方程，运行太慢，所以我去掉了。竖直方向：高0.1米，划分了20格网格。谢谢你！！

maples1

If I were u, I will try 40, 60 or 80 grids.
Good luck.

lyzhang98

我试了，但是没有用。我想主要原因在于，我的计算区域形状效应。因为它的长宽比为：300：1。当我把长缩小3倍时，就不会有这种情况。那么这种计算结果可靠吗？因为实际测量数据没有办法去得那么多，我怎么才能判断计算结果的好坏呢？

zliu

你不妨先用一些经典的低Re数工况,如Re=100,1000检验你的程序,并做一些定量的比较,以排除程序的问题.然后在此基础上,逐步提高Re数,观察有关规律.
从你说的网格情况看,高宽比显然太大了,通常要限制在1:1.3以内.

lyzhang98

不是的，是我计算的区域300：1,网格我取得相等。
谢谢您的建议
正在调试....