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

xuyu

    我最近一直在学习这个程序，相信也有很多人也在通过研究它来进行最基本的CFD学习，在源代码学习过程中，遇到很多理论和编程技术的问题，前面也看过不少关于该程序的讨论，但是太分散，能不能集中讨论一下呢？
    我在这里先提个小问题，由于是初学者，希望大家多多指点：
    数组F(I,J,NF)代表变量U,V等，程序中NF=11设置为压力P，但是在程序的计算过程中P(I,J)并没有传输给F(I,J,11)，可是在结尾的打印输出时还是输出了F(I,J,11)，不知道为什么，请各位执教！
    希望和SIMPLE算法的研究者交流，谢谢！

CFDNHT

怎么没传送？...
前面的说明语句就有说明等价语句。
如果你不改nf＝10，那么就是默认nf＝11为rho

xuyu

等价语句是：
EQUIVALENCE (F(1,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3),PC(1,1))
这里只是把U,V和校正压力P';传给了F,但是压力P没有传过去阿
“如果你不改nf＝10，那么就是默认nf＝11为rho”这句我没看懂

zlf

这个程序哪里有啊？是算什么的呀！在编程方面，我是菜鸟都算不上的，希望各位前辈们不吝赐教！

zla2003

能把整段代码哪出来看看的吗
我最近用simple做了个三维的计算的呀

cnlimin

[这个贴子最后由cnlimin在 2005/03/15 02:49pm 第 1 次编辑]

小弟不才，用F95格式改写了该程序，有兴趣可以到交叉学科论坛看看，还请各位同行多多提意见。西交大网上的程序，在visualfortran6.5版本运行有错误发生——打印错误。必须改一改。

cnlimin

天哪，SIMPLE是NHT里的主要内容啊，4楼的，难道你们没学吗？

xuyu

cnlimin ，你改写的就是西交大的程序吗？你改写得是不是主要就在定义部分，希望能和你多多交流

cnlimin

[这个贴子最后由cnlimin在 2005/03/15 08:36pm 第 1 次编辑]

xuyu，我只是用新的语法把西交大给的程序写了一遍，至于定义部分我觉得它的一个地方不大合适，在VF6.5下运行会报数组溢出的错误，而在POWERSTATION4.0下运行则不会报错。我也希望能和你多多交流。

xuyu

你对我上面提到的数据传递的问题是怎么看得？请赐教

cnlimin

我觉得，在声明里要加上相应的等价声明。通用程序只给出了前3个，后面的都要用户具体指定。不知道，你用的FORTRAN是哪个版本，如果用fortran6.5不改教学程序是打印不出压力的，如果压力的NF〉10。

cnlimin

我对这程序能够在POWERSTATION4中顺利打印出压力，感到很奇怪。

xuyu

我就是非常奇怪阿，用了4.0就是能打印出来，但是我没有看到P(i,j)是怎么传输到F(i,j,11)的，是不是有什么特殊语句阿，我只看到了U,V和Pc是等价过去的

dongshanyan

我也在编制SIMPLE算法，不过是用MATLAB在编制。我觉得西交大的FORTRAN方法太复杂(对我而言)，有C语言编制的SIPLE版本吗？
另外，我有一些小问题想请教大家，第一源项的中Sp、Sc的给定有何标准，以及源项中使用何种差分格式，另外关于湍流的壁面函数如何设置。请各位高手指教。如果有人正在做或是做过，请大家相互讨论或者帮助我一下吧。我的Email:dongaiyan@163.net. QQ：10181694。

xuyu

Sp Sc的选择标准计算传热书上有一些说明，关于源项的差分格式你指的是不是当源项比较复杂有偏微分项时的处理方法