再次求助simplec

大葱蘸酱 · 发表于 2004-6-4 09:21:10

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

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

x

计算参数 A，临近边界，粘性力那一项一开始为
a_Para[j]=a1+a2+a3+a4;
a4=Vis[j]*(Ve_u[j+1]-2*Ve_u[j]+Ve_u[j-1])/(dy[j]*dy[j]);
因为Ve_u[j+1]=1 Ve_u[j]=0，Ve_u[j-1]=0
所以 Ve_u[j+1]-2*Ve_u[j]+Ve_u[j-1])=1
粘度Vis[j]取1，dy[j]取0.01这一项就变成10000，时间步取为0.01
那么dtStep*a_Para[j]/Density[j]就变得很大很大
Ve_u_t[j]=Ve_u[j]+dtStep*a_Para[j]/Density[j]也变得很大
于是不收敛，怎么办呢

大葱蘸酱 · 发表于 2004-6-4 09:27:27

a1+a2+a3中只出现过一次Ve_u[j+1]，但是和v乘后变成零了，因此a4的作用突出
是不是密度不能太小

凤翼天翔 · 发表于 2004-8-30 11:18:59

老兄，你也太厉害了吧？曲高和寡啊

大葱蘸酱 · 发表于 2004-10-6 20:31:51

偶研究了半年，感觉simple大多数问题都知道了，有问题可以问我了，呵呵

liuhuafeifei · 发表于 2004-10-9 07:33:52

算算网格Pe书 density*V*d/vis
看看扩散项的占的比例就清楚了

zgp0816 · 发表于 2004-10-26 17:00:14

能给我推荐一本simplec的资料吗，我这里没有，还想学，不知道，有没有扶贫的啊，本人zgp0816@163.com,qq38798357
非常感谢

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

再次求助simplec

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。