ANSYS中FLOTRAN的共轭传热分析

ANSYS中FLOTRAN的共轭传热分析
单根发热管敷设在排管中,排管埋设在土壤中,假设土壤四周温度恒定,分析排管内空气的流体温度场,以及土壤中温度场分布.
但是运行后显示错误,说温度出现负值,下面是源程序.
想知道是什么原因造成的.
/filname,paiguanliuti07
/title,paiguanliuti07
/prep7
TOFFST,273
建立模型
pcirc,5
pcirc,7
rectng,-20,20,-20,20
local,11,cylin,0,-3,0
cyl4,0,-3,0,0,1,360
allsel,all
aovlap,all
numcmp,all
定义单元类型
et,1,fluid141
定义材料属性
mp,kxx,2,400
mp,kxx,3,20
mp,kxx,4,30
进行暂态分析,要指定材料属性的密度和热容
给材料赋值
asel,s,1
aatt,2,1
asel,s,2
aatt,3,1
asel,s,4
aatt,4,1
划分网格
smrtsize,1
allsel,all
mshape,0,2D
mshkey,0
amesh,all 疑问,用流体单元后,空气部分未指定材料属性就可以划分网格
指定边界条件（暂时不考虑辐射能的传热作用）
csys,11
lsel,s,loc,x,1
nsll,s,1
d,all,vx,0,vy
csys,1
lsel,s,loc,x,5
nsll,s,1
d,all,vx,0,vy 压力载荷在排管内表面如何施加,空气面的边界温度需要指定吗
csys,0
lsel,s,11
sfl,all,conv,10,25
lsel,s,10,12,2
sfl,all,hflux,0
lsel,s,9
dl,all,temp,15
施加体积生热率
asel,s,1
bfa,all,hgen,150000
allsel,all
ACEL,0,9.81,0 设置重力加速度
FLDATA1,SOLU,TRAN,0
FLDA,SOLU,TEMP,T
FLDATA7,PROT,DENS,AIR-SI 空气属性国际制
FLDATA7,PROT,VISC,AIR-SI
FLDATA7,PROT,COND,AIR-SI
FLDATA8,NOMI,COND,-1
FLDATA7,PROT,SPHT,AIR-SI
FLDATA13,VARY,DENS,1 密度变化
FLDA,SOLU,TEMP,T 激活求解温度方程
FLDA,SOLU,FLOW,F 关闭流体方程
FLDA,RELX,TEMP,1.0 流体物性常值,温度松弛因子1,一个全局迭代步内即可求解温度方
程
不必解流体方程,两次全局迭代导致第二次结果有小的变化.
FLDA,ITER,EXEC,1 确定稳态分析的迭代和输出,迭代次数1
FLDA,METH,TEMP,3 求解病态传热问题选择预条件共轭残差法（PCR）
SAVE
/solu
SOLVE
FLData2,iter,exec,200
FLData2,iter,over,50
FLDATA5,OUTP,SUMF,50
FLDATA14,TEMP,NOMI,25 指定初始温度
FLDATA14,TEMP,BULK,T0 指定用于计算传热系数的环境温度
FLDATA15,PRES,REFE,1.0135E5 指定参考压力
FLData18,meth,pres,1 指定压力求解器
FLData18,meth,temp,4 指定温度求解器
FLDATA19,TDMA,PRES,100 指定压力用TDMA法求解时的迭代数
FLData21,conv,pres,1e-10 定义压力自由度的收敛准则
FLData21,conv,temp,1E-02 定义温度自由度的收敛准则
FLData25,Relx,Temp,0.8 设定温度自由度松弛系数
FLData25,Relx,Dens,0.5 设定流体密度属性松弛系数
FLData33,ADVM,TEMP,MSU
Flda,wadv,temp,on
/solu
solve

了然哦,徐麒雯,应该是的有｀｀