对经典例题的修改出现了问题，求高手指点，万分感谢

cuizhikun00

这是一位仁兄的原题：
求解圆管中流体以及管壁的温度分布，考虑流体与管内壁的对流，管外壁与
周围环境对流。已知条件如下：
管内流体的物理属性以及入口条件
导热系数          0.6W/(m*c)
密度              1000kg/m3
比热              4183J/(kg*c)
入口温度          20c
入口流速          3m/s
壁管的物理属性和几何属性
导热系数          100W/(m*c)
管外径            5/16inch
壁厚              0.083inch
长度              12inch
管外对流条件
环境温度          100c
对流系数          10000w/(m2*c)
管内对流条件      
对流系数          15000w/(m2*c)
  
/prep7
!转换数据为SI
pi=acos(-1)
ri=((5/16)-0.083)*0.0254 =0.0058293 !管内半径
r0=(5/16)*0.0254=0.0079375           !管外半径
l=12*0.0254 =0.3048              !管长度
nthick=1                   !径向单元个数
size=(r0-ri)/nthick =0.0021082       !径向单元长度
nlength=nint(l/size) =145      !长度方向单元个数
  
ktube=100                  !管导热系数
kwater=0.6                 !流体导热系数
denswatr=1000              !流体密度
cwater=4183                !流1
  
vinput=3                   !入口流速
carea=pi*(ri**2)=0.           !流通面积
massflow=vinput*carea*denswatr  =0.3204 !质量流率
  
tinlet=20                  !流体入口温度
tbulk=100                  !管外环境温度
hi=15000                   !管内壁与流体的对流系数
h0=10000                   !管外壁对流系数
!定义三类单元
et,1,116,1,1               !KEYOPT(1)＝1自由度为温度
!KEYOPT(2)＝1       使用SURF151模拟管内对流
et,2,151,,,1,1,1           !KEYOPT(3)＝1，轴对称
                     !KEYOPT(4)＝1，无中间节点
                     !KEYOPT(5)＝1，有额外节点，在FLUID116上
keyopt,2,8,2              !KEYOPT(8)＝2，考虑对流
et,3,55,,,1                !KEYOPT(3)＝1，轴对称
r,1,2*ri,carea,1           !FLUID116的实常数
  
mp,kxx,1,kwater
mp,c,1,cwater
mp,dens,1,denswatr
mp,kxx,3,ktube
!建模
k,1,0,0,0
k,2,0,-l,0                !注意不是-1，而是-L
k,3,0,-(l+size),0         !注意不是-(1+size)，而是-(L+size)
l,1,2
l,2,3                      !线1和2模拟管内流体流动
cm,waterlin,lines          !创建组元，用于生成表面单元
rectng,ri,r0,0,-l          !管壁,-L
type,1
real,1
mat,1
esize,,nlength
lmesh,1
esize,,1
lmesh,2
lsel,s,loc,x,(ri+r0)/2     !选择管两端
lesize,all,,,nthick        !厚度方向单元个数
esize,,nlength             !单元大小
type,3
mat,3
allsel,all
amesh,1
  
lsel,s,loc,x,ri            !选择管内壁
cm,tubelin,lines
type,2
lfsurf,'tubelin','waterlin' !使用宏LFSURF创建表面效应单元 （我修改时出现问题的地方）
  
/solu
dk,1,temp,tinlet           !流体入口温度
lsel,s,loc,x,r0
sfl,all,conv,h0,,tbulk     !管外壁对流边界
esel,s,ename,,surf151
sfe,all,,conv,0,hi         !管内壁对流系数
esel,s,ename,,fluid116
sfe,all,,hflux,,massflow   !管内流体质量流率，以热流密度方式加载
allsel,all
solve
finish
/post1
nsel,s,loc,x,0             !选择FLUID116
prnsol,temp                !列出流体温度分布
nsel,s,loc,x,r0            !选择管外壁节点
prnsol,temp                !列出管外壁节点温度
##################################################################################################################################

我准备在原题上修改，主要有在以下几处进行修改：
（1）        原题是水在管道里流，我把它改为水在长方体里的圆孔里流（模型见附图，取1/2模型，图中直线表示水）

（2）        改后模型的传热情况是：长方体上表面有对流换热（与空气换热），水与圆孔表面换热，长方体其他面都为绝热面（包括对称中心面）。
（3）流体的物理属性以及入口条件 ：
导热系数          0.6W/(m*c)   
密度              1000kg/m3
比热              4183J/(kg*c)
入口温度         40
入口流速         0.5
长方体的物理属性和几何属性 ：
导热系数         0.52
比热             1840J/Kg℃
密度               
长方体长       10m
长方体 高         0.5m（圆孔轴线位置在0.25m处）
长方体 宽         0.5m（模型的一半为0.25m）
长方体上方环境温度         25
长方体上方对流系数          1000w/(m2*c)
圆孔半径           0.06m
圆孔里对流系数         87408.6w/(m2*c)
求解：在保持水进口温度恒定，流速恒定的情况下，求水沿管长方向温度变换，以及周围长方体每一截面的温度变化？（相关的参数在我下面的命令流里）

我已经修改到下面一步(见附件)了，后面遇到技术问题已经进行不下去，望大侠鼎立相助！
或者大侠们能帮忙写一下命令流，谢谢！

cuizhikun00

高手请不要躲啊，诚请帮忙．．．．继续等．．．

cuizhikun00

斑竹帮帮忙啊，小弟在此谢过，顺祝大家元旦快乐，圣诞快乐

cuizhikun00

07年到了，先祝大家学业有成，心想事成．可惜我０６年的遗憾还没弥补，顶一下，大家帮忙啊

cuizhikun00

各位大哥大姐，行行好啊．我急啊，赶紧帮个忙吧，哪怕给个意见

superyellor

你好，我也在做这方面的研究，能不能跟你交流下

tanly17

楼主能不能具体阐述下，对流的我做个一点