xxj2005 发表于 2007-6-5 11:20:50

大家帮忙看看这个命令流!谢谢

一个很简单的模型,做流场和温度场分析。结果流场计算结果与想象中的不太一样,怎么回流那么不明显呢?是不是哪个地方设置错了?请大家多多指点,谢谢!!
/BATCH                        !
/CLEAR
/TITLE,
/PREP7
EMUNIT,MKS                  !定义国际单位制,米
*AFUN,DEG
et,1,142
CYLIND,0.037,0.008,0,-0.055,0,30,   !生成溶液
CYLIND,0.008,0,0,-0.055,0,30,
NUMMRG,KP, , , ,LOW
VGLUE,ALL
NUMCMP,VOLUME
NUMCMP,AREA
/REPLOT
mptemp,1,0,1490,1526,1560,   !定义的热性能
mpdata,kxx,1,1,29.4,29.4,23.26,23.26
mpdata,c,1,1,0,670,842,842
VSEL,S,VOLUM, ,2
ESIZE,,10
VMESH,ALL
VSEL,S,VOLUM, ,1
LSEL,s,LINE,,1
LESIZE,ALL,,,30
VMESH,ALL
ALLSEL
asel,s, , ,8         !施加入口边界条件
da,all,Vz,-0.2,1
da,all,Vy,0,1
da,all,VX,0,1
da,ALL,TEMP,1550
asel,s, , ,3      !壁面的边界条件
da,all,Vz,0,1
da,all,Vy,0,1
da,all,VX,0,1
asel,s, , ,2    !液面
da,all,Vz,0,1
csys,1
asel,s,,,5,6,
asel,a,,,9,10,
da,all,Vy,0,1
csys,0
asel,s, , ,1,7,6       !出口
Da,all,PRES,0,1
asel,s, , ,3
sfa,all,,conv,5000000,1280
ALLSEL,ALL
finish
/SOLU
FLDATA1,SOLU,TURB,1             !湍流分析
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,1
fldata1,solu,rdsf,0            
FLDATA2,ITER,EXEC,200         !总体迭代次数
FLDATA3,TERM,VX,0.001         
FLDATA3,TERM,VY,0.001         
FLDATA3,TERM,PRES,0.001      
FLDATA3,TERM,ENKE,0.001      
FLDATA3,TERM,ENDS,0.001         
FLDATA3,TERM,teMP,0.001         
FLDATA7,PROT,DENS,CONSTANT
FLDATA8,NOMI,DENS,7100
FLDATA7,PROT,VISC,CONSTANT
FLDATA8,NOMI,VISC,0.0055
FLDATA7,PROT,COND,TABLE
FLDATA7,PROT,SPHT,TABLE
FLDATA13,VARY,COND,1
FLDATA13,VARY,SPHT,1
FLDATA12,PROP,DENS,0
FLDATA12,PROP,VISC,0
FLDATA12,PROP,COND,0
FLDATA12,PROP,SPHT,0
FLDATA24,TURB,MODL,1            
FLDATA24,TURB,CMU,0.09
FLDATA24,TURB,C1,1.44
FLDATA24,TURB,C2,1.92
FLDATA24,TURB,SCTK,1
FLDATA24,TURB,SCTD,1.25
FLDATA24,TURB,WALL,eqlb
FLDATA24,TURB,cks,0.5
fldata24,turb,inin,0.02
fldata24,turb,insf,0.02
FLDATA24,TURB,RATI,2
FLDATA25,RELX,VX,0.3         !设置流体性质松弛系数
FLDATA25,RELX,VY,0.3
FLDATA25,RELX,PRES,0.5
FLDATA25,RELX,ENKE,0.3
FLDATA25,RELX,ENDS,0.3
FLDATA25,RELX,TEMP,0.5
FLDATA26,STAB,TURB,0.01      !设置有助于求解稳定的参数   
FLDATA26,STAB,VISC,100         !指定人工粘度系数

FLDATA34,MIR,MOME,0.2
FLDATA34,MIR,TURB,0.2
FLDATA18,METH,PRES,3          !选择每个自由度的解算方法
FLDATA18,METH,vx,1
FLDATA18,METH,vy,1
FLDATA18,METH,enke,1
FLDATA18,METH,ends,1
FLDATA18,METH,TEMP,4
flda,capp,pres,t             !设置自由度变量的限值以防止出现不合理的结果
flda,capp,pmax,1.0e+8
flda,capp,pmin,-1.0e+8

flda,quad,prsd,2          !控制压力扩散项的面积积分阶次
flda,quad,prss,2          !控制压力源项的面积积分阶次

fldata33,advm,mome,supg   
fldata33,advm,turb,supg
fldata33,advm,pres,colg
flda,turb,rati,1.0e3                              
flda,relx,evis,0.1         
FLDATA22,MAXI,PRES,100         
FLDATA20,SRCH,PRES,12            
FLDATA21,CONV,PRES,1e-02         
FLDATA23,DELT,PRES,1e-02,      
FLDATA22,MAXI,TEMP,100      
FLDATA20,SRCH,TEMP,12         
FLDATA21,CONV,TEMP,1e-02         !定义温度方程求解收敛准则
FLDATA23,DELT,TEMP,1e-02,
SOLVE
FINISH

/POST1
SET,LAST
/VSCALE,1,1,0   
PLVECT,V,,,,VECT,ELEM,ON,0
PLNSOL,TEMP,,0

xxj2005 发表于 2007-6-5 20:25:02

怎么都没有人看得啊?自己顶一下吧!
高手帮帮忙啊!
页: [1]
查看完整版本: 大家帮忙看看这个命令流!谢谢