实战回报：FDM简单方腔流程序

onesupeng

看到这个版上很多新人想编程无从下手，看别人程序有时候很难读，我空的话做些最简单程序，意在帮助那些想编程的人。我提供的程序将会非常简单，但我的目的并不是代替大家解决自己的问题，而是从最简单的例子出发，引导新手进行发展自己的程序，使新手有个概念。

对于大多数教学程序，我将用我的风格自己写或者把别人的程序经过修改，增加可读性，并加以验证。每个程序都将留有一些的改进空间，让新手读过之后，可以动手进行一定的操作和改进，以达到锻炼的目的。

这次提供的是投影法有限差分算不可压N-S方程方腔流的简单程序。对流项采用二阶迎风格式，粘性项采用二阶中心差分，时间推进采用A-B格式，压力采用G-S超松弛点迭代方法。

新手可考虑改进以下几个地方：
1、压力边界条件高阶格式；
2、三阶，甚至四阶迎风格式
3、时间推进可改为3阶TVD或者一些隐式推进方法
4、改算槽道流动，贴加重力，设置滑移边界等等
5、有一定基础之后改进压力迭代方式等
XXXXXXXXXX

[ 本帖最后由 onesupeng 于 2010-3-8 16:32 编辑 ]

通流

程序写的倒是清楚。可是没有解释，想要搞懂也不是那么简单的。反正我是没有看懂。似乎是预测-校正的时间离散。既然是不可压的，那么一部分的空间离散应该是隐式得（或者全场耦合的），这是不是就是pressure那个子程序做的事？

我的建议，你最好再写个关于理论方面的文件作为补充。

onesupeng

是在没有时间写详细的东西，过两天给参考文献

shirazbj

好，学习。

我觉得教学程序可以针对一两个概念，程序未必优化，算的未必高效。不必通用。要是能手算几个点才好呢。

在网上搜了一下投影法, 这么说：the projection method scheme to solve navier-stokes equation using finite difference is for un-steady state.

这篇介绍了一点：
http://www.stanford.edu/class/me469b/handouts/incompressible.pdf

[ 本帖最后由 shirazbj 于 2010-3-17 19:35 编辑 ]

mfc08

不是交错网格，有没有交错的