模拟腔体内基板的辐射散热问题

liusong0138

腔体内有一基板,基板温度T1,腔体壁面的温度是T2,已知腔体壁面和基板的辐射率,要求考虑辐射和空气的热传导,模拟基板300s内的降温过程.
问题的关键是模型中必须对空气和玻璃基板同时建立模型,在模拟辐射时,我是用plae55单元和辐射求解器模拟辐射过程的,用了sfl,rdsf命令进行加载辐射率,虽然模拟出来了(命令流见下),但我总觉得结果有一些问题,主要问题是,用sfl,rdsf命令进行加载时,必须要设定一个space temp,在模型中,我将这个space temp设定为两个辐射面之间的初始温度15度.

但这个空间温度一般是对于不建立模型的区域来说的,也就是说两个辐射面之间必须是"space",才能设定空间温度.但在我的模型中,壁面和基板表面实际上有空气的,也就是说,两个辐射面之间不是"space",温度也是随着时间变化的,我尝试过改变space temp,发现这个温度对最后结果影响非常大,所以我觉得这种方法有问题.
我也尝试过使用link32生成热辐射矩阵的方法,但发现由于我的玻璃表面属于内表面,根本无法生成表面效应单元.
经过查阅,发现有一种方法,要用link31单元将一对辐射面相互连接起来,但具体方法我不清楚,请问有哪位高手遇到过此类问题,能否给予解答?
我是新手,主要做热流耦合方面的,有哪位同仁愿意可以加我qq446894158,欢迎一起讨论一起学习.下面附该模型的命令流:

h=0.003
l=1.3
hb=0.0735
lb=0.1
T1=15
T2=600
/prep7
!选择单元
et,1,plane55
!建立模型
rect,,lb,0,hb
rect,,lb,0-h,0
rect,,lb,0-h-hb,0-h
rect,lb,lb+l,0,hb
rect,lb,lb+l,0-h-hb,0-h
rect,lb+l,lb+l+lb,0,hb
rect,lb+l,lb+l+lb,0-h,0
rect,lb+l,lb+l+lb,0-h-hb,0-h
ALLSEL,ALL
aglue,all
cm,fluid,area
!流体区域建模完成
rect,lb,lb+l,0-h,0
allsel,all
cmsel,u,fluid
cm,glass,area
!固体区域建模完成
allsel,all
nummrg,all!合并关键点
numcmp,all! 压缩编号
!划分网格
lsel,,loc,x,,lb+l+lb
lsel,r,loc,y,1e-8,hb-1e-8
lesi,all,,,10,5.0
lsel,,loc,x,,lb+l+lb
lsel,r,loc,y,0-h-1e-8,0-h-hb+1e-8
lesi,all,,,10,5.0
lsel,s,loc,y,hb,0-h-hb
lsel,r,loc,x,1e-8,lb-1e-8
lesi,all,,,5
lsel,s,loc,x,lb,lb+l
lsel,r,loc,x,lb+2*1e-8,lb+l-2*1e-8
lesi,all,,,50
lsel,s,loc,x,lb+l,lb+l+lb
lsel,r,loc,x,lb+l+1e-8,lb+l+lb-1e-8
lesi,all,,,5
!定义材料
cmsel,s,fluid
mat,1
amesh,all
cmsel,s,glass
mat,2
amesh,all
!定义固体材料属性
mp,dens,2,2700
mp,kxx,2,0.75
mp,c,2,790
mptemp,1,260,400,420,460,480,600
mpdata,emis,2,1,1,1,0.95,0.85,0.77,0.8
!定义液体材料属性
mp,dens,1,1.29
mp,c,1,1.4
mp,kxx,1,0.03
!加载玻璃边界条件
cmsel,s,glass
ALLSEL,BELOW,AREA
d,all,temp,600
cmsel,s,glass
lsel,s,ext
lsel,r,loc,x,lb+2*1e-8,lb+l-2*1e-8
sfl,all,rdsf,-2,,1!加载辐射条件
!加载流体边界条件
lsel,s,loc,y,hb
nsll,s,1
d,all,temp,15
sfl,all,rdsf,0.85,,1!加载辐射条件
lsel,s,loc,y,0-h-hb
nsll,s,1
d,all,temp,15
sfl,all,rdsf,0.85,,1!加在辐射条件
allsel
!生成AUX12矩阵
/aux12
STEF,0.567E-07
TOFFST,273
spctemp,1,15!设定空间温度
radopt,,0.001
!求解
/solu
antype,trans
timint,off
time,0.01
deltim,0.01
allsel
solve
timint,on
time,300
nsubst,60
autots,on
ddele,all,temp
outres,,all
allsel
solve
finish

noisybaot

帮顶~

是不是要用到耦合场分析？？

liusong0138

因为空气是静止的,而且不考虑对流,所以用不到cfd单元.只有热场的作用
问题就是怎样考虑两个辐射面之间变化的空气温度的?如果规定了空间温度,计算辐射的时候应该就是按照指定的空间温度计算辐射热量的,可是在这里空气的温度是随着辐射传热变化的,感觉做起来有点棘手.