西交大陶文铨老师simple教学程序的讨论

cnlimin

MatLab的循环速度很慢的，用这个编写的SIMPLE会是什么执行速度？不敢想。

cnlimin

因为要进行大量的迭代计算，我觉得你最好还是用FORTRAN编吧。

zhyb176

通过COMMON F（N，N，NF），使得诸变量共享；
可在SUBROUTINE SUPPLY中，给出
EQUIVALENCE (F(1,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3),PC(1,1))
     &,(F(1,1,11),P(1,1)),(F(1,1,12),RHO(1,1)),(F(1,1,13),GAM(1,1))
使得，LPRINT（11）=.TRUE.时，压力得到输出。
可能有些是我改变了源程序，还请大家给与反馈。

xuyu

原始程序中的p,rho,gam都没有进行等价，也被打印了，请问怎么解释呢

zlf

cnlimin，你好，你说的NHT是不是数值传热学？那这个simple程序就是数值传热学里讲的压力-速度耦合的计算程序，是吗？不过对于编程我还是连菜鸟都不是的，别笑话我，希望多多的向大家学习！

通流

simple是过时的东西。用它来教学生CFD是不合适的。希望大家不要浪费你们的宝贵时间来研究这东西。
有功夫，可以学习Jameson的方法，或者Bob Ni的论文。用时间推进，虽然收敛的慢一些，但结果好的多。
我说的比较直，望大家谅解。

shirazbj

我从这下的：
http://nht.xjtu.edu.cn/nht/downpage.asp

其中下载28说要用power station4.0编译，开始没有。但以前有个g95,是从www.g95.org下的，用它在dos终端里一编译，没问题。能运行。

但编译下载30有个错，它是把TITLE声明成双精度类型了，改称字符型就编译过去了。

shirazbj

有三种坐标。我只感兴趣直角坐标，解2d速度和压力场。要是只针对这个就好了。

我现在读到setup1了，y方向的设置比较乱。

还有后面的插值系数也没搞清呢。是从p网格插到u网格么？还是反过来？还是都需要。插值不能取两侧值的平均值么？怎么还要涉及x向控制体的宽度什么的？

[ 本帖最后由 shirazbj 于 2009-12-3 12:13 编辑 ]