- 积分
- 0
- 注册时间
- 2004-9-24
- 仿真币
-
- 最后登录
- 1970-1-1
|
我在做铸件充型过程的流场和温度场的耦合模拟
我想先做一个比较简单的例子比较一下网格和时间步长对结果的影响,现在我做了
一个非常简单的例子,如下图,尺寸为20cm×60cm。在初始时刻,假定空腔上面
1/3的体积为金属流体,速度为零,然后整个流体自由下降撞击到下壁面。
因为对命令流不是特别熟悉,所以采用了命令流和GUI混合采用的方式,我做的过
程是如下:
1、在prefrences里面选择flotran cfd。
2、preprocessor——element type——add fluid141
3、preprocessor——modeling——create areas——rectangle 0,0.2,0,0.6
,生成了0.2m×0.6m的面积
4、preprocessor——meshing——mesh tool——line set——pick all,再输入
0.01,把整个面积划分为0.01m×0.01m的网格
5、preprocessor——loads——difine loads——apply——thermal——temp,然
后选中上面1/3面积的节点,输入温度值1623K.
6、preprocessor——loads——difine loads——apply——cfd——velocity——
on lines,然后选中底边和左右边共3条线,vx和vy都输入0,固壁边界条件,比较
简单。
7、输入一下命令流:
nsel,s,loc,x,,0.2
nsel,r,loc,y,0.4,0.6
esln,s,1
ice,all,vfrc,1.0
alls
ic,all,vx
ic,all,vy
grav=9.81
/solu
acel,0,grav,0
flda,nomi,dens,7800 !铁合金
flda,nomi,visc,0.025
flda,nomi,cond,22.26
flda,nomi,spht,726.6
flda,solu,flow,t
flda,solu,temp,t
flda,solu,vof,t
flda,solu,tran,t
flda,time,glob,5
flda,time,numb,30
flda,time,step,0.01
flda,time,appe,0.1
flda,relx,vx,0.2
flda,relx,vy,0.2
flda,relx,pres,0.5
flda,advm,mome,supg
flda,mir,mome,1.0
solve
然后进入general postproc,read the next step,就是0.1秒的时候的温度分布
,结果见图,查询节点的值也还比较合理,然后再读下一个step,即0.2秒时的温
度分布,就发现温度的分布显然已经不合理了,到了0.3秒这种趋势就更明显了,
液体的前端和气体接触的部分温度应该是最低的,但是却变成最高了,结果同样见
图。然后我想可能是因为速度增大时间步长过大引起的,但是反复的细化时间步长
再从0.1秒时重启也不能解决问题,总是会在某个时间部长之后液体最前端温度就
变成最高的了。图里面是把时间细化为0.005秒的结果,再喜欢为0.001秒我也做过
,还是会最终变成类似的情况。真诚请各位虾虾指教,不胜感激!
还有就是我试着在一开始就把时间步长设定小一些,比如设为0.001,但是就会出
现液体一直留不下去的情况,看流体的速度,发现居然是向上的,实在是莫不着头
脑,希望虾虾们一并指教。
谢谢谢谢!!!!!!!!!!!!!!! |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|