为什么在FLOTRAN CFD中辐射加了没作用(附命令流)

tripstones

我参照版友freefly的程序，做了一个简单的模型，可是辐射加和没加，基体的温度相差不到1度，大家帮我看看，是哪里出了问题？



/prep7
grav =9.81                !重力加速度
Dens_H2=0.0899                !氢气密度
Vis_H2=0.000022                !氢气动力粘性系数
R_sub=0.04                !衬底半径
h_sub=0.02                    !衬底厚度
n=14                        !热丝根数
r_fila=0.0003                !热丝半径
sp_btw=0.006                !热丝间距
d_fila_sub=-0.005        !热丝－衬底距离
r_st=0.006                !进气口直径
v_air=2.36                !氢气流速
wid_out=0.02                !出气口直径
fila_size=0.0001        !热丝单元大小
h1=0.04                        !进气口位置
h2=0.0006
h3=-h_sub-d_fila_sub        !衬底底边位置
h4=-0.06                !出气口位置
kpnum1=0
kpnum2=0
kpnum3=0
kpnum4=0
kpnum5=0
kpnum6=0
lnum1=0
lnum2=0
anum1=0
anum2=0


!定义单元类型及单元属性设置
et,1,141
!定义材料导热系数
mp,kxx,1,0.2             !氢气
mp,kxx,2,109                  !基体：钼
mp,kxx,3,60                 !热丝：钽丝
!建立流体模型
!热丝的关键点
k,,,0
k,,2*r_fila
k,,2*r_fila,h2
k,,,h2
k,,-2*r_fila,h2
k,,-2*r_fila
k,,-2*r_fila,-h2
k,,,-h2
k,,2*r_fila,-h2
circlr,1,r_fila,,,360,8
!热丝周围的面
*do,i,1,7,1
l,1+i,2+i
*enddo
l,9,2
*do,i,1,8,1
l,1+i,9+i
*enddo
*do,i,1,7,1
al,i,16+i,8+i,17+i
*enddo
al,8,24,16,17

rectng,-2*r_fila,-4*r_fila,2*r_fila,0
rectng,-2*r_fila,-4*r_fila,0,-2*r_fila
rectng,2*r_fila,4*r_fila,2*r_fila,0
rectng,2*r_fila,4*r_fila,0,-2*r_fila
allsel
NUMMRG,ALL
NUMCMP,ALL

!生成热丝上面的区域
kpnum1=kp(-4*r_fila,h2,0)
kpnum2=kp(4*r_fila,h2,0)
*get,kpnum3,kp,0,num,max
*get,lnum1,line,0,num,max
k,,-4*r_fila,8*r_fila,
k,,-0.5*r_fila,8*r_fila,
k,,0.5*r_fila,8*r_fila,
k,,4*r_fila,8*r_fila,

l,kpnum1,kpnum3+1
l,kpnum3+1,kpnum3+2
l,kpnum3+2,kpnum3+3
l,kpnum3+3,kpnum3+4
l,kpnum3+4,kpnum2

lsel,s,,,lnum1+1
lsel,a,loc,y,8*r_fila
lsel,a,,,lnum1+5
lsel,a,loc,y,h2
al,all

!生成热丝下面的区域
kpnum4=kp(-4*r_fila,-h2,0)
kpnum5=kp(4*r_fila,-h2,0)
*get,kpnum6,kp,0,num,max

k,,-4*r_fila,-8*r_fila,
k,,-0.5*r_fila,-8*r_fila,
k,,0.5*r_fila,-8*r_fila,
k,,4*r_fila,-8*r_fila,

l,kpnum4,kpnum6+1
l,kpnum6+1,kpnum6+2
l,kpnum6+2,kpnum6+3
l,kpnum6+3,kpnum6+4
l,kpnum6+4,kpnum5

lsel,s,,,lnum1+6
lsel,a,loc,y,-8*r_fila
lsel,a,,,lnum1+10
lsel,a,loc,y,-h2
al,all

rectng,2*r_fila,-2*r_fila,-4*r_fila,-6*r_fila
*get,anum1,area,0,num,max
asba,anum1-1,anum1

ALLSEL,ALL
CM,FLUID_WHOLE,AREA

allsel
lesize,all,fila_size
amesh,all

!流体建模完毕
!
!建立固体模型

lsel,s,,,1,8,,1
al,all
ALLSEL,ALL
CMSEL,U,FLUID_WHOLE
CM,filaments,AREA
lsla,s
lesize,all,fila_size
AATT,3
amesh,all

lsel,s,,,45,48,,1
al,all
*get,anum2,area,0,num,max
asel,s,,,anum2,,,1
CM,substrate,AREA,
allsel
lesize,all,fila_size
AATT,2
amesh,all

!固体建模完毕
!
ALLSEL,ALL
FINI

/SOLU
!
!施加热丝温度
cmsel,s,filaments
!da,all,temp,2700
lsla,s
nsll,s
!BFA,14,HGEN,1e8
CM,fila_node,NODE
D,fila_node,temp,2700

!施加边界条件
lsel,s,loc,x,-4*r_fila
lsel,a,loc,x,4*r_fila
lsel,a,,,lnum1+2,lnum1+4,2
lsel,a,,,lnum1+7,lnum1+9,2
nsll,S,1
CM,bj_node,NODE
D,bj_node,vx,0
D,bj_node,vy,0
D,bj_node,temp,100

!施加进气口条件
lsel,s,,,lnum1+3
nsll,S,1
CM,in_node,NODE
D,in_node,vx,0
D,in_node,vy,-2.36
D,in_node,temp,100
sfl,lnum1+3,vfrc,1,,1
!施加出气口条件
lsel,s,,,lnum1+8,,1
nsll,S,1
CM,out_node,NODE
D,out_node,pres,0
D,out_node,temp,100

ALLSEL,ALL
!施加辐射边界条件
lsel,s,,,1,4,1,1
nsll,s
SFl,ALL,rad,0.65,,2700
lsel,s,,,46,48,2,1
nsll,s
SFl,ALL,rad,0.82,,822,
cmsel,s,filaments        !热丝周围线段
!lsla,s
lsel,s,,,5,8,1,1
nsll,s
SFl,ALL,rdsf,0.65,,1,

cmsel,s,substrate     !基体上表面线段
lsel,s,loc,y,-4*r_fila
nsll,s
SFl,ALL,RDSF,0.82,,1,

!施加基体底面温度
ALLSEL,ALL
cmsel,s,substrate
lsel,s,loc,y,-6*r_fila
nsll,s,1
d,all,temp,822

!边界条件施加完毕

ALLSEL,ALL

!耦合结点
!cpintf,all,1e-8
FINI

/SOL
!设置参考温度为摄氏温度
TOFFST,273
!热分析
FLDATA1,SOLU,TEMP,1
!
flda37,algr,segr,simplen
flda25,relx,vx,0.95
flda25,relx,vy,0.95
flda25,relx,pres,0.95
flda25,relx,temp,0.95
!
!设置氢气属性

fldata8,nomi,dens,Dens_H2   
fldata8,nomi,visc,Vis_H2
fldata8,nomi,cond,0.2  
fldata8,nomi,spht,1.4  

!设置环境温度
!
FLDATA14,TEMP,NOMI,100,   
FLDATA14,TEMP,TTOT,100,   
FLDATA14,TEMP,BULK,100,   

!设置辐射参数

radopt,0.5,1.0e-5,0,10000,,0.9
v2dopt,0,0,0,400
STEF,0.567E-07   
SPCTEMP,1,100  
  
!重力加速度
ACEL,0,9.81,0,
solve
FINI
!
/post1
allsel
cmsel,s,substrate!filaments!
nsla,s,1
SET,LAST
etable,,erase,1
etable,tt,temp
pletab,tt,avg