自己编写的方腔流动程序（两维，非定常，交错网格，用时间分裂法计算）

yfjma · 发表于 2004-5-9 13:46:31

您需要登录才可以下载或查看，没有账号？注册

x

参看附件。

hoodlum1999 · 发表于 2004-5-9 15:14:10

正在下载中,,,,帮你顶
看了程序再说,不懂的地方请指教

ustcsunl · 发表于 2004-5-9 17:34:44

让我意外的是，居然是用c写的程序，500行，辛苦辛苦。俺的程序（算法
可能与你不同）也差不多如此长，可是Orlandi写了个同样功能的程序，居然
用了几千行，佩服佩服。

simply · 发表于 2004-5-11 10:08:15

求教楼主关于方腔边界条件如何施加
经典的方腔流动以流函数Psi为基本变量，
则边界条件为 Psi＝0 & D(Psi)/D(N)=0 or 1
若用u,p 做基本变量呢，p的边界条件如何加？

yfjma · 发表于 2004-5-11 10:58:52

我是用原始变量求解的，对于流函数的边界条件我没有研究过。

simply · 发表于 2004-5-11 15:45:31

若用原始变量，在等稳情况下，应该是速度U和压力P
u的条件倒好加，p的边界条件如何加呢？

yfjma · 发表于 2004-5-12 10:52:37

p的采用Nueman条件，即p(i,1)=p(i,2);p(i,jdm) =p(i,jdm-1)
p(1,j)=p(2,j);p(idm,j)=p(idm-1,j)
其中，idm,jdm表示x和y方向p的个数。

yuarmy · 发表于 2004-5-13 10:55:23

老兄，你给出的参考文献太老了，我也找不到阿。请问你知道还有那些文献介绍的方法和你的程序采用的方法相近的吗？谢谢。

simply · 发表于 2004-5-13 11:40:28

感谢楼主，步骤非常详细
另外在Zienkiewicz的“有限元”5th，也提到了方腔流，直接令左下角的节点压力为零，看来对于稳态问题，中间过程是否施加Nuemann条件影响不大

yfjma · 发表于 2004-5-13 22:09:27

你到google中输入fractional step可以找到一些东西。

leaningtower · 发表于 2004-5-25 23:15:09

你用的是什么方法？

a9002 · 发表于 2005-4-20 11:51:51

看起来好像是涡量流函数法，不过好慢啊。

xiaobin · 发表于 2005-5-15 10:20:51

请问有没有用三阶紧致格式算的啊
帮帮忙啊

CFDNHT · 发表于 2005-5-15 21:32:08

***** 版主模式 *****
该贴子是管理员从<a href=forums.cgi?forum=2>CFD基础理论</a>转移过来的！

gotowuhan · 发表于 2006-2-15 17:21:14

不知道这位仁兄的程序能算到多大雷诺数的流动

		自动登录	找回密码
密码			注册