[讨论]一个涉及流场温度场耦合，辐射，对流换热的例子

yuguoxingqin

请大家帮忙看看问题出现在哪里？我怀疑是边界条件有问题，可是没有检查出来。
FINISH
/CLEAR
/FILNAME,GISSIMU
/TITLE,SINGLE BUS GIS
/UNITS,SI
/PREP7
RN0=0.53/2-0.016 !导体内半径
RW0=0.53/2  !导体外半径
RN1=0.98/2-0.009 !外壳内半径
RW1=0.98/2  !外壳外半径
HGN=317.47
HGW=139.49
TOFFST,273
ET,1,FLUID141
MP,KXX,2,500  !母线的导热系数
MP,C,2,400
MP,DENS,2,2700
MP,KXX,3,500  !外壳的导热系数
MP,C,3,400
MP,DENS,3,2700
/PNUM,LINE,1  !生成几何模型
/PNUM,AREA,1
PCIRC,0,RN0
PCIRC,0,RW0
PCIRC,0,RN1
PCIRC,0,RW1  
AOVLAP,ALL
NUMCMP,AREA
NUMCMP,LINES
NUMCMP,NODES
ASEL,S,,,2
AATT,2,,1,0
ASEL,S,,,4
AATT,3,,1,0
ASEL,S,,,1,3,2
AATT,1,,1,0
LSEL,S,,,1,16  !划分母线导体单元网格
LESIZE,ALL,,,40  
MSHAPE,0,2D
MSHKEY,0
ALLSEL,ALL
AMESH,ALL
ASUM,FINE
*GET,S2,AREA,2,AREA
*GET,S4,AREA,4,AREA   
FINISH
/SOLU
!--------------辐射----------------
ALLSEL,ALL
LSEL,S,,,13,16
SFL,ALL,RAD,1,,27.3
ALLSEL,ALL   !流场内辐射
LSEL,S,,,5,8
SFL,ALL,RDSF,0.85,,-1
LSEL,S,,,9,12
SFL,ALL,RDSF,0.85,,-1
ALLSEL,ALL
FINISH
!------------------------------------
/AUX12
STEF,5.67E-8
RADOPT,0.5,0.01
HEMIOPT,10
ALLSEL,ALL
FINISH
!-------------------------------
/SOLU
ASEL,S,,,2   !施加焦耳热
BFA,ALL,HGEN,HGN/S2
ASEL,S,,,4   !施加焦耳热
BFA,ALL,HGEN,HGW/S4
ALLSEL,ALL
LSEL,S,,,13,16   !外壳对流换热
NSLL,S,1
SF,ALL,CONV,3.59,27.3
ALLSEL,ALL
LSEL,S,,,1,12   !无滑移边界条件   
NSLL,S,1
D,ALL,,,,,,VX,VY
ALLSEL,ALL
FLDATA1,SOLU,TEMP,1  !设置求解温度场
FLDATA1,SOLU,TURB,1  !设置求解湍流
FLDATA2,ITER,EXEC,1000  !总迭代次数
FLDATA2,ITER,OVER,50  !文件覆盖频率
FLDATA3,TERM,VX,0.01   
FLDATA3,TERM,VY,0.01   
FLDATA3,TERM,VZ,0.01   
FLDATA3,TERM,PRES,1e-008   
FLDATA3,TERM,TEMP,1e-008   
FLDATA3,TERM,ENKE,0.01  
FLDATA3,TERM,ENDS,0.01
FLDATA5,OUTP,SUMF,50  !输出总结的频率

FLDATA12,PROP,DENS,4
FLDATA13,VARY,DENS,1  !密度可变
FLDATA12,PROP,VISC,4
FLDATA13,VARY,VISC,1
FLDATA12,PROP,COND,4
FLDATA13,VARY,COND,1
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方向重力加速度
FLDATA15,PRES,REFE,101350*3  
FLDATA16,BULK,BETA,1E15   
FLDATA17,GAMM,COMP,1.4  
FLDATA14,TEMP,NOMI,27.3   
FLDATA14,TEMP,BULK,27.3
FLDATA18,METH,PRES,3  !指定求解器
FLDATA18,METH,TEMP,6
ALLSEL,ALL
SOLVE
FINISH
/POST1
SET,LAST
PLNSOL,TEMP
FINISH

yuguoxingqin

怎么没有人顶？这个问题困扰很久了，内、外辐射这个加法对不对，求真相！