问题描述:
内部铜块发热,引起封闭腔内部空气流动,希望得到流固的温度场分布
几何模型由3部分构成:A3为copper,A4为steel,A5为air。
边界条件:
L2,L4施加对流边界条件:50w/(m*m*k),bulk temp=20;
A3(铜块)生热率是100w;
L5,6,7,8,9,10,11,12施加边界无滑移条件:vx=0,vy=0;
我用141单元直接求解的,但是得不到结果,不知哪里出了问题,还是这种方法根本不可行,请大家帮忙指点,不胜感激!
/BATCH
/PREP7
/units,si
toff,273
!define element
ET,1,FLUID141
! define material
MP,KXX,2,386 !define copper
MP,KXX,3,20 !define steel
!creat model
RECTNG,0,1,0,1,
RECTNG,0.1,0.9,0.1,0.9,
RECTNG,0.3,0.4,0.3,0.4,
AOVLAP,1,2,3,
nummrg,all
!mesh contorl
lsel,s,,,1,4,1
lesize,all,,,50,1
lsel,s,,,5,8,1
lesize,all,,,40,1
lsel,s,,,9,12,1
lesize,all,,,5,1
alls
!!Boundary Conditions - Solid Model
lsel,s,,,5,12,l
dl,all,,vx,0
dl,all,,vy,0
lsel,s,line,,2,4,2
sfl,all,conv,50,20
!nsll,s,1
!sf,all,conv,50,293
bfa,3,hgen,100
!asel,s,area,3
!nsla,s,1
!bf,all,hgen,100
alls
asel,s,area,,3
aatt,2,1,1
amesh,all
asel,s,area,,4
aatt,3,1,1
amesh,all
asel,s,area,,5
aatt,1,1,1
amesh,all
finish
alls
/SOLU
FLDATA1,SOLU,TEMP,1 ! Set solution option to thermal
FLDATA2,ITER,EXEC,200 ! Set FLOTRAN execution controls
FLDATA2,ITER,OVER,50
FLDATA5,OUTP,SUMF,50
!*
FLDATA13,VARY,DENS,1 ! Set fluid properties
FLDATA7,PROT,DENS,AIR-SI
FLDATA7,PROT,VISC,AIR-SI
FLDATA7,PROT,COND,AIR-SI
!FLDATA8,NOMI,COND,-1
FLDATA7,PROT,SPHT,AIR-SI
!*
ACEL,0,9.81,0 ! Set flow environment gravity
FLDATA18,METH,PRES,1 ! Set FLOTRAN solver controls
FLDATA19,TDMA,PRES,100
solve
fini |