请问有没有人算过非牛顿流体 (Non-Newtonian)？

ethan · 发表于 2006-5-19 04:42:09

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

x

请问有没有人算过非牛顿流体 (Non-Newtonian)？我用VISMOL.F，粘性(viscosity)随SECINV变化，可是好像不论什么流速， SECINV都是零。测了粘性公式，没错。请指教，谢谢。

飞翔鸟 · 发表于 2006-5-19 10:16:39

你如何检查SECINV的值？
这个值好象只能在粘性子程序中调用，是不能输出的吧。
你检查你的粘性，应该是变化的。

ethan · 发表于 2006-5-19 23:16:13

根据粘性公式，在高速流，粘性值会趋进某一值，在低速流，粘性值会趋进另某一值。
我假设SECINV=0时，粘性=某常数，SECINV 不等于0时，粘性随SECINV变化。
可是，不论什么流速，得出的粘性都是在SECINV=0时的常数。
有办法检查SECINV的值吗?谢谢．

飞翔鸟 · 发表于 2006-5-20 11:12:37

把你的子程序传上来我看一下吧

ethan · 发表于 2006-5-25 00:25:58

就这简单的.
   IF (SECINV .EQ. 0) THEN
      VISM=0.004
   ELSE
VISM=(SQRT(0.004*(1-EXP(-100.*ABS(SECINV)))/ABS(SECINV))+SQRT(0.004))**2.
   ENDIF
谢谢!

飞翔鸟 · 发表于 2006-5-26 18:08:58

你换一个判别条件再试一下，判别条件中不要用这个变量

ethan · 发表于 2006-6-3 00:42:29

我换成了
　if (u,v,w = 0) then vism = 0.004
else ....
还是不对。。　：（

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