求助：flotran分析中的问题

liangge1069

小弟是第一次做流动分析，下面是命令流，solve之后总是出现问题，还请各位帮忙看看到底问题出在哪儿。
先谢谢了！


/COM,  Thermal  
/COM,  FLOTRAN CFD  
/PREP7
r=15 ! radius of spring
ro=1 !outer radius of the tube
ri=0.9
n=6 ! segment number in the tube
s=2.1 ! lead
a=1 !the increas of the radius by per circle
aa=0 ! the total increas of the radius
nc=2 !the number of the circle
pi=3.14159265
an=1 !the number of the area which you want to drag
!设置气体参数
nomidens=1.7533
nomivisc=18.63e-6
nomicond=0.143
nomispht=5.2
cof1dens=273
cof1visc=273
cof1cond=273
cof2dens=10*1.01e5
cof2visc=108.6
cof2cond=153

!设置铜管物性参数
cudens=8.92e3
mp,c,1,cp
mp,dens,2,cudens
MPTEMP,1,300,500,900,1000
MPDATA,c,2,1,4.185*92,4.185*93,4.185*106,4.185*115
MPDATA,kxx,2,1,4.185*95.6,4.185*93.1,4.185*87.4,4.185*80.3
CSYS,1
!draw the points on the centre line
*do,i,0,nc*n,1
  aa=aa+a/n
  x=r+aa
  y=i*60
  z=2.2*i/n
  k,i+1,x,y,z
*enddo
*do,j,1,nc*n,1
L,j,j+1
*enddo
!draw the cross section of the tube
LWPL,-1,1,0
cyl4,0,0,ro,,ri,360 ! draw section of tube
cyl4,0,0,ri,,,360 !draw section of fluid
aglue,all
numcmp,all
!draw the grid of the cross section(inlet or outlet)
ET,1,MESH200
ET,2,142
KEYOPT,1,1,6
KEYOPT,1,2,0
LESIZE,ALL,,,8,,,,,1  
MSHAPE,0,2D  
MSHKEY,0
aesize,1,1/30
AMESH,1
mshkey,1
amesh,2
!sweep the cross section and generate the volume mesh
*do,i,1,nc*n,1
mat,2
VDRAG,an,,,,,,i
mat,1
vdrag,an+1,,,,,,i
NUMMRG,all
numcmp,all
an=an+10
*enddo
aclear,1,2

/solu
fldata1,solu,temp,t
fldata1,solu,turb,t
fldata1,solu,comp,t
fldata2,iter,exec,200 !设置迭代次数
fldata2,iter,over,50
fldata7,prot,dens,gas!设置流体属性
fldata7,prot,visc,gas
fldata7,prot,cond,gas
fldata7,prot,spht,constant
fldata12,nomi,dens,nomidens
fldata12,nomi,visc,nomivisc
fldata12,nomi,cond,nomicond
fldata12,nomi,spht,nomispht
fldata12,cof1,dens,cof1dens
fldata12,cof1,visc,cof1visc
fldata12,cof1,cond,cof1cond
fldata12,cof2,dens,cof2dens
fldata12,cof2,visc,cof2visc
fldata12,cof2,cond,cof2cond
fldata12,cof3,dens,cof3dens
fldata12,cof3,visc,cof3visc
fldata12,cof3,cond,cof3cond
fldata13,vary,dens,1
fldata13,vary,visc,1
fldata13,vary,cond,1
fldata13,vary,spht,1

asel,s,,,3,113,10
asel,a,,,6,116,10
nsla,s,1
sf,all,hflux,10000
allsel
asel,s,,,122
nsla,s,1
d,all,temp,300
d,all,pres,1*1.01e5
d,all,vz,10
nsla,s,,,2
nsla,s,1
d,all,pres,1.01e5
allsel
solve