采用有限容积法还是LB方法求解NS方程哪？

cqulcm1980

如题！
现在需要求解一定区域里的NS方程，手中也有simple程序，可是个人觉得读这个程序也很痛苦，很多都不清楚为何如何了？
听说，LB求解Ns方程很简单，但是其原理很难了。
而求解NS方程只是为了获得速度场，只是一个中间步，求得速度场之后，是去求浓度场之类的了。而对流－扩散形式的浓度场已经编程完毕，是采用有限容积法了。现在只等着速度场的求解了。

onesupeng

有限差分有限元呢

cqulcm1980

楼上的，你好！
因为此前已经做了一些采用有限容积求浓度场的程序，所以还是打算采用有限容积法来求速度场，只是觉得simple较为麻烦了。

onesupeng

那何以又想到LB呢？

通流

我也不懂。为什么要琢磨LB呢？弄个商业软件不就行了吗？

cqulcm1980

楼上的几位，先谢谢关注此问题了。
我也想用fulent啊，但是因为涉及到求速度场之后，再求浓度场的时候，出现的方程比较复杂，因为涉及到多个场的耦合作用了。这里不赘述了噢。
之所以想用LB，因为看到LB求NS方程比较快，而且适用于复杂边界吧。
现在还在看，看哪一个更省力些啊。
不过simple，也觉得很麻烦。
Ps：羡慕西交的学生啊！

cqulcm1980

也就是说我只是用这个求速度场，然后再求其他了。
也考虑过采用fluent算速度场，但是因为求解域中涉及到液体区域和固体区域了，再求浓度场的时候，是采用整场求解了。而在fluent中，听说，液体固体交界面之处，是两个面了，觉得不容易和我现有的程序结合起来了。

通流

Ps：羡慕西交的学生啊！
------------------------------------

为什么？

cqulcm1980

因为每年陶老师都会亲自讲啊。
我们可是自学啊。
而且，各个觉得最难的就是NS方程的求解了！

通流

就那过时的SIMPLE方法？不知是福气还是不幸。

cqulcm1980

楼上的，对于我们最开始接触求解NS方程来说，可能去读一些别人的程序，消化完之后再自己编程好些吧。

现在目前所能找到的程序就是simple程序，还有peric的程序。

但是peric的书看了看，觉得东西很灵活，也并没有一些具体的离散格式，可能才看两天，但是对着他的程序，很多都看不懂了。

cqulcm1980

那你觉得自己编程的话，用什么好啊？我现在不管是哪种方法，只要能解决我的问题就好了，毕竟求出速度场只是我的一个初步的中间的结果了。

通流

我的感觉还是商业软件容易。可能你需要的是后处理。

我也不反对SIMPLE方法。我也用。

会SIMPLE方法的，也不是只有西交大有。很多年前，教我SIMPLE方法的是科大的老师。如果是二维问题，编程也不难。

写程序的人都知道，读别人的程序比自己写程序更累。

后面的事，只有靠你自己决定了。

cqulcm1980

我们这边simple基本上靠自学吧。

simple方法是pantanker提出的，只是陶老师的弟子做了个注释版的并挂在网上而已。

程序是读别人的不错，但是在自己一头雾水的情况下，从没接触过这样的情况下，也许看些别人的程序也好，当然不能是去用别人的了。
  
本想找个捷径，没想到还是回到原路了。

所以少上网，多看书是王道了。

网上继续推导公式了。

cqulcm1980

商业软件容易的地方在于本身是一个正确的框架，而自己编程一个小错误可能会导致结果错误了。
编了小半年的程序的感觉了。