流场-温度场分析的问题

yuguoxingqin

最近做一个流场温度场的仿真，边界条件是：空气内边界施加固定温度80度，外边界固定30度，内外边界都施加无滑移边界条件，另外就是重力加速度。最后结果不是对称的，哪位高手帮忙分析分析，谢谢！QQ:327574871

FINISH
/CLEAR
/FILNAME,GISSIMU
/TITLE,SINGLE BUS GIS
/UNITS,SI
/PREP7
RW0=0.045  !导体外半径
RN1=0.246  !外壳内半径
TOFFST,273
ET,1,FLUID141
/PNUM,LINE,1  !生成几何模型
/PNUM,AREA,1
PCIRC,RW0,RN1,0,180
PCIRC,RW0,RN1,180,360
AGLUE,ALL
NUMCMP,ALL
LSEL,S,,,2,4,2
LESIZE,ALL,,,40
ESIZE,,40
MSHAPE,0,2D
MSHKEY,1
!MAT,1
!TYPE,1
ASEL,S,,,1,2
AMESH,ALL
ALLSEL,ALL
LSEL,S,,,1,5,4
LSEL,A,,,3,6,3
NSLL,S,1
D,ALL,,,,,,VX,VY  !无滑移边界条件
ALLSEL,ALL
LSEL,S,,,3
NSLL,S,1
D,ALL,TEMP,80   !导体温度
LSEL,S,,,6
NSLL,S,1
D,ALL,TEMP,80   !导体温度
ALLSEL,ALL
LSEL,S,,,1
NSLL,S,1
D,ALL,TEMP,30   !外壳温度
LSEL,S,,,5
NSLL,S,1
D,ALL,TEMP,30   !外壳温度
FLDATA2,ITER,EXEC,200  !总迭代次数
FLDATA2,ITER,OVER,50  !文件覆盖频率
FLDATA5,OUTP,SUMF,50  !输出总结的频率
FLDATA12,PROP,DENS,4
FLDATA13,VARY,DENS,1  !密度可变
FLDATA12,PROP,VISC,4
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,4
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,4
FLDATA13,VARY,SPHT,0
FLDATA7,PROT,DENS,AIR-SI !密度
FLDATA8,NOMI,DENS,-1
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,AIR-SI !粘性
FLDATA8,NOMI,VISC,-1
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA7,PROT,COND,AIR-SI !导热系数
FLDATA8,NOMI,COND,-1
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,AIR-SI !比热
FLDATA8,NOMI,SPHT,-1
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
ACEL,0,9.8,0   !Y方向重力加速度
ALLSEL,ALL
FINISH
/SOLU
ALLSEL,ALL
FLDATA1,SOLU,TEMP,1  !设置求解温度场
FLDATA1,SOLU,FLOW,1  !设置求解压力场
FLDATA18,METH,PRES,3  !指定求解器
FLDATA18,METH,TEMP,6
ALLSEL,ALL
SOLVE
FINISH
/POST1
SET,LAST
PLNSOL,TEMP
!PLNSOL,V,X
!PLNSOL,V,Y
!PLVECT,V,,,,VECT,ELEM
FINISH