对初学CFD朋友们的一个建议。

通流

看到好多人都在学CFD。这是好事。但把SIMPLE方法当成入门的东西，可能就不是好事了。老的东西不一定就是基础的东西。

我觉得，对初学者，还是先学时间推进容易一些。尤其是非定常的时间推进。那东西是直接解欧拉方程，出来的结果，能和流体力学对的上。尤其是，在没收敛的时候，也能用基本的流体力学来理解你看到的现象。过去，计算机太慢，现在不一样了。至少七八年前，Jameson就敢在作讲座时，当场在他的1公斤重的 Laptop上优化3元跨音机翼。

Jameson的格式就是一个很好的入门格式。隆格库塔法几乎使空间离散与时间离散相对独立。初学者可以试试不同阶的时间离散，及不同的空间离散。甚至可以试试有限元法。用数学理论来看格式的精度，稳定性，耗散，频散等都容易做。

这儿，搞计算的高手不少。希望提提你们的高见。

shirazbj

参考着学.

通流

（1）其实，做CFD的人都知道，读别人的程序，比自己写程序难。
（2）可以发个帖，看看有人有个现成的，解二维欧拉方程的时间推进程序。
（3）网上看看，应该能找到这样的程序。
（4）我有时间，看看能不能找到过去的程序。那东西还是在VAX上运行的。还有人知道VAX机器吗？

shirazbj

读别人的程序确实难。我在读那个一维管流的simple程序,一大堆问题.看原理有时觉得明白了,但一看程序又不知道哪跟哪了.

还是喜欢直接解N-S方程的.

通流

SIMPLE是比Jameson方法复杂。我一直说，老的东西不一定简单，也不一定基础。找Jameson的1981年的文章吧。里面讲的很详细。对于不可压流，可以用假如拟可压项来解决。格式是很简单的。难一点的是进出口边界条件。要理解的话，需要解线性化的欧拉方程，求特征值，及特征向量来解决。这可能会需要的时间来做，但做好后，整个程序会非常皮实。

onesupeng

很头头是道嘛

但我以为首先搞最简单的一维热方程波动方程-〉二维-〉其他如流动。而且每做一步，数值计算的验证都要做好掉

否则就像前面别版老表们讨论的，用Fluent做得多牛鼻，结果算个圆柱没有算对，那就没有意义了，，，，，，

做这个切忌别不会爬就想飞，没有用，经点算例的对比验证是必须的

通流

从一维开始做的建议不出啊。能不能给出验证用的一个或几个算例？

不过，我还是要强调，CFD麻烦的地方还是在边界。多花点时间在这方面是必须的。处理的办法最好是从数学和物理都能讲得通的办法（其实，数学与物理是一致的）。试凑的办法是万万不可取的。

其实解不可压流比解可压流要困难不少，所以，没有必要，甚至不应该，从不可压做起。原因是不可压流方程的刚性比较大，刚性大的方程搞数值解当然就难些。

shirazbj

只是对流动感兴趣。那个一维管流simple程序挺有意思。

找到了一个日本人的2d euler程序，说要用dvf5.0编译，我用g95编译没问题。300X100的网格，推了1000步，画出图来象个有中央射流的流动。隆格库塔法是4步？看着里面有4个步骤。上下壁是SLIP CONDITION，脑袋里只有non-slip的概念。进口给了压力，温度，密度，U1,U2有值,V1是0，还给了E1,E2.出口条件也给了. Roe's Average,Flux Difference Splitting是矢通分裂?这程序是2000年的. 是不是老了点?设完网格后有Jacobian calculation,对矩形网格还要算这个?不是贴体网格才变换么呢?到是感兴趣用贴体网格算个台阶的例子.

通流

难道没有信心自己写一个吗？

shirazbj

自己写程序要推方程。

看别人的程序时很难看清别人是怎么推的。即使有介绍，也是简化了的。还有介绍时是原理性的，比如那个simple程序,介绍时用的p,w,e表示节点。编程时要用i，i+1。如果这个i被交错网格共用，有时就搞不清在那个网格上了。如果原理时一个量叫ae,编程时叫做aim,读的时候就更不知道哪跟哪了。

感觉推方程一步一步清楚了，程序也就差不多了。

花了100流币，下了一个fvm的pdf,啃两口，看看对理解fvm离散两维动量方程有没有帮助。

通流

支持你的努力。如果有可能的话，把你的推导过程，及程序都贴上来。我想这东西对很多人有用。

niuaiming

总是感觉用CFD 模拟出来的东西不太准确

tennet

自立更生，丰衣足食！！！！

quaintchewster

哪位大侠写本CFD编程的傻瓜书啊

liuchen

傻瓜书就是直接给你一个调通的程序，但是你也就别指望能增加点啥功能了