利用FLUID 141单元进行热分析，运行时报错

yadan99

请教：利用FLUID 141单元进行热分析，运行时报错
*** ERROR ***                           CP =     318.844   TIME= 23:16:49
Negative or zero temperatures may be avoided by mesh refinement         
especially if Petrov-Galerkin formulation is used for the energy        
equation (FLDA,ENRG,SUPG,T).

命令流如下：

问题描述：利用喷枪火焰加热钢板，求解流场和温度场

!!计算流场至稳定状态
!前处理过程
/prep7 et,1,141  !定义单元类型为FLUID 141
keyopt,1,3,2  !利用关于x轴称的二维分析
!!创建关键点
k,1,0,0
k,2,1.2,0
k,3,9.7,0
k,4,10,0
k,5,0,0.25
k,6,1.2,0.25
k,7,9.7,0.25
k,8,10,0.25
k,9,1.2,10
k,10,9.7,10
k,11,10,10
!!创建线
l,1,2
l,2,3
l,3,4
l,5,6
l,6,7
l,7,8
l,1,5
l,2,6
l,3,7
l,4,8
l,6,9
l,7,10
l,8,11
l,9,10
l,10,11
  !!创建面
a,1,2,6,5
a,2,3,7,6
a,6,7,10,9
a,3,4,8,7
a,7,8,11,10
  !!对线段进行划分
lsel,s,,,7,10
lesize,all,,,6,0.75
lsel,s,,,1,4,3
lesize,all,,,15,-2
lsel,s,,,2,5,3
lsel,a,,,14
lesize,all,,,80,-3
lsel,s,,,3,6,3
lsel,a,,,15
lesize,all,,,5,-1.5
lsel,s,,,11,13
lesize,all,,,80,4  
!!对流体区域划分网格
asel,s,,,1,3
mat,1  !材料号为1
amesh,all
!!对固体区域划分网格
asel,s,,,4,5
mat,2  !材料号为2
amesh,all alls
!!定义边界条件及初始条件
!!对火焰喷口上边界、中心对称线和金属与流体的交界处定义y方向速度为零
lsel,s,,,1,2
lsel,a,,,4,9,5
lsel,a,,,12
nsll,,1
d,all,vy,0
!!对火焰喷口上边界和金属与流体的交界处定义x方向速度为零
lsel,s,,,4,9,5
lsel,a,,,12
nsll,,1
d,all,vx,0
!!定义气体进口速度为3500cm/s，进口处在x方向有速度，y方向没有速度
lsel,s,,,7
nsll,,1
d,all,vx,3500
d,all,vy,0
!!定义计算空气的边界处压力为零
lsel,s,,,11,14,3
nsll,,1
d,all,Pres,0  
alls  !选取全部节点
!!定义固体材料的属性
mp,dens,2,7.8 !密度
mp,kxx,2,0.03 !热传导系数
mp,c,2,0.5 !比热容
!!定义流体属性，选取流体的属性为空气，单位为cm/g/s制
fldata,prot,dens,air-cm  !密度
fldata,prot,visc,air-cm  !粘性
fldata,prot,cond,air-cm  !传热系数
fldata,prot,spht,air-cm  !比热
fldata,vary,dens,t  !设置密度可变

!!设置流动环境参数
fldata,pres,refe,1.0325e+6  !设置参考压力为一个大气压
fldata,temp,ttot,350  !设置总温为350
toffst,273  !设置偏移温度  
!!设置求解选项和CFD求解器
fldata,solu,turb,t  !设置湍流
fldata,solu,comp,t  !设置压缩流
fldata,meth,pres,3  !设置压力半直接求解器
fldata,conv,pres,1e-12  !设置收敛准则
fldata,relX,dens,1  !设置密度松弛系数为1
fldata,stab,mome,1  !设置动量惯性为1
fldata,stab,pres,1e-4  !设置压力惯IE-4
fldata,stab,visc,1  !设置人工粘性为1
fldata,iter,exec,20  !设置总体迭代次数为20
save  !保存
!求解
/SOLU
SOLVE
!!!逐步减小人工粘性直到0来求解稳定流场
fldata,stab,visc,0.1
solve
fldata,stab,visc,0.01
solve
fldata,stab,visc,0.001
solve
fldata,stab,visc,0.0001
solve
fldata,stab,visc,0.00001
solve
fldata,stab,visc,0
fldata,iter,exec,140
solve  
!!!增大动量惯性和压力惯性继续求解
fldata,stab,mome,10
fldata,stab,Pres,0.01
solve
fldata,stab,mome,1e+15
fldata,stab,Pres,1e+15
fldata,iter,exec,300!设置总体迭代步数为300
solve
FINIsH !!至此,可以得到一个稳定流动的流场  
!!在喷口处施加上温度
/PREP7
NPLOT  !显示节点
FLST,2,84,1,ORDE,6  !选取在喷口处的节点
!!!界定选取节点的范围
FITEM,2,122
FITEM,2,-133
FITEM,2,207
FITEM,2,-218
FITEM,2,318
FITEM,2,-377
/GO  !!!在选取的节点上施加温度条件
D,P51X,TEMP,3000  !定义温度为3000K
ALLSEL,ALL  !选取全部节点
!!!求初始温度解
FLDATA1,SOLU,TEMP,1
FLDATA1,SOLU,FLOW,0
FLDATA25,RELX,TEMP,1.0,
FLDATA2,ITER,EXEC,1,
/SOLU
SOLVE !!求解温度
finish

注意：之前运行结果自认为正确，但之后的命令流运行结果不正确，求教！！

/PREP7  !!设置求解选项
FLDATA1,SOLU,TRAN,1  !设置为瞬态求解
FLDATA1,SOLU,FLOW,1  n!求解流场
FLDATA1,SOLU,TEMP,1  !求解温度
FLDATA1,SOLU,TURB,1  !湍流项
FLDATA1,SOLU,COMP,1  !可压缩
/COM,,Transient Analysis,1  !!设置输出控制
FLDATA4,TIME,STEP,-1  !设置为时间控制
FLDATA4,TIME,ISTEP,0.001,  !分步长为0.001秒
FLDATA4,TIME,TEND,10,  !计算的时间长为10秒
FLDATA4,TIME,OVER,0.1,  !设置文件覆盖的间隔0.1秒
FLdATA4,TIME,APPE,0.1,  !设置附加频率间隔0.1秒
FLdATA4,TIME,SUMF,0.1,  !设置总结频率间隔0.1秒
FINISH
/SOLU
SOLVE