【讨论】一个涉及流固耦合,流体对流,固体传导,固体间的辐射的典型问题

fossiler

大家好,我现在做一个温度场的分析,里边涉及到了流畅,还有辐射,的问题,辐射能量比较大,所以不能简化,全部因素都可以考虑进去了,但是算出来的结果有问题,我感觉是辐射那里出的问题,因为我加载辐射算出来的温度与不加辐射算出来的相差零点几度,当然我感觉辐射是起做用了,但是可能是没有加正确还是什么原因呀,还是ansys的问题,请大家多多指教呀,命令流如下

fossiler

模型图

fossiler

铜导体一方面向铜外壳辐射热量,同时通过中间的流体空气向外客传递能量,
外壳外表面想空间辐射能量

fossiler

mx_h=0.016
mx_wj=0.53/2
wk_h=0.009
wk_wj=0.98/2  
mx_jr=324  
wk_jr=245  
   
/PREP7
et,1,FLUID141
   
mp,kxx,2,400
mp,c,2,400
mp,dens,2,8890
   
PCIRC,MX_WJ-MX_H,0,0,360,  
PCIRC,MX_WJ,0,0,360,
PCIRC,WK_WJ-WK_H,0,0,360,  
PCIRC,WK_WJ,0,0,360,  
asel,s,,,1,4,1  
AOVLAP,all
lsel,s,,,1,16,1
lesize,all,,,40,,,,,1
asel,s,,,5,7,2
AATT,       2, ,   1,       0,  
asel,s,,,1,6,5
AATT,       1, ,   1,       0,   
!ADELE,1, , ,1   
  asel,s,,,1
  asel,a,,,5,7,1
!asel,s,,,5,7,1
mshape,0,2d
mshkey,0
amesh,all
ALLSEL,ALL
ASUM, 0  
*GET,area5,AREA,5,AREA  
*GET,area7,AREA,7,AREA  
allsel,all
!定义空间节点
!*GET,node1,NODE,,NUM,MAX, , , ,  
!n,node1+1,0,1,0  
/SOL
toffst,273
!加焦耳热
!asel,s,,,5
!nsla,s,1
!BF,all,HGEN,mx_jr/area5
!asel,s,,,7
!nsla,s,1
!BF,all,HGEN,wk_jr/area7
BFA,5,HGEN,mx_jr/area5
BFA,7,HGEN,wk_jr/area7
allsel,all
lsel,s,,,1,4,1
dl,all,1,VX,0,1   !定义无滑移边界条件
dl,all,1,VY,0,1
allsel,all
lsel,s,,,5,12,1
dl,all,6,VX,0,1   !定义无滑移边界条件
dl,all,6,VY,0,1
  
!加对流换热系数
lsel,s,,,13,16,1
nsll,s,1
SF,all,CONV,4,40
allsel,all
!加内部辐射
  
asel,s,,,5
esla,s,1
lsel,s,,,5,8,1
nsll,s,1   
sf, all, rdsf, 0.85,2
allsel,all  
asel,s,,,7
esla,s,1
lsel,s,,,9,12,1
nsll,s,1   
sf, all, rdsf, 0.85,2
allsel,all
  
!加外部辐射
asel,s,,,7
esla,s,1  
lsel,s,,,13,16,1
nsll,s,1
sf, all, rdsf, 1, 1
ALLSEL,ALL
  /AUX12  
STEF,1.12e-11   
RADOPT,0.5,0.01,0,  
SPCTEMP,1,40
!STEF,1.19e-11      
RADOPT,0.5,0.01,0,  
   
!SPCNOD,1,node1+1
HEMIOPT,10      
vfcalc,flo,vf      
/SOL  
!D,node1+1, ,40, , , ,TEMP, , , , ,  
allsel,all  
   
!流体性质  
FLDATA1,SOLU,TRAN,0  
FLDATA1,SOLU,FLOW,1  
FLDATA1,SOLU,TEMP,1  
FLDATA1,SOLU,TURB,1
FLDATA1,SOLU,COMP,0  
FLDATA1,SOLU,VOF,0   
FLDATA1,SOLU,SFTS,0  
FLDATA1,SOLU,IVSH,0  
FLDATA1,SOLU,SWRL,0  
FLDATA1,SOLU,SPEC,0  
FLDATA1,SOLU,ALE,0   
FLDATA1,SOLU,RDSF,1
/COM,,Steady State Analysis,0   
FLDATA2,ITER,EXEC,300,   
FLDATA2,ITER,OVER,30,   
FLDATA2,ITER,APPE,30,   
FLDATA3,TERM,VX,0.01,   
FLDATA3,TERM,VY,0.01,   
FLDATA3,TERM,VZ,0.01,   
FLDATA3,TERM,PRES,1e-008,   
FLDATA3,TERM,TEMP,1e-008,   
FLDATA3,TERM,ENKE,0.01,  
FLDATA3,TERM,ENDS,0.01,  
FLDATA5,OUTP,SUMF,10,     
  
   
FLDATA12,PROP,DENS,4
FLDATA13,VARY,DENS,1
FLDATA12,PROP,VISC,4
FLDATA13,VARY,VISC,0
FLDATA12,PROP,COND,4
FLDATA13,VARY,COND,1
FLDATA12,PROP,SPHT,4
FLDATA13,VARY,SPHT,0
  
FLDATA7,PROT,DENS,AIR-SI
FLDATA8,NOMI,DENS,-1
FLDATA9,COF1,DENS,0  
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,AIR-SI
FLDATA8,NOMI,VISC,-1
FLDATA9,COF1,VISC,0  
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS   
FLDATA7,PROT,COND,AIR-SI
FLDATA8,NOMI,COND,-1
FLDATA9,COF1,COND,0  
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,AIR-SI
FLDATA8,NOMI,SPHT,-1
FLDATA9,COF1,SPHT,0  
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
   
FLDATA15,PRES,REFE,101350,   
FLDATA16,BULK,BETA,1e+015,   
FLDATA17,GAMM,COMP,1.4,  
FLDATA14,TEMP,NOMI,40,   
   
FLDATA14,TEMP,BULK,40,  
   
ACEL,0,9.8,0,   
  FLDATA18,METH,PRES,1
  FLDATA19,TDMA,PRES,100,  
!FLDATA18,METH,temp,1
!FLDATA19,TDMA,temp,100,   
  FLDATA18,METH,TEMP,2
  FLDATA22,MAXI,TEMP,1000,
  FLDATA20,SRCH,TEMP,2,   
  FLDATA21,CONV,TEMP,1e-012,   
  FLDATA23,DELT,TEMP,1e-010,  
solve
  /post1
  set,last  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  
  

fossiler

大家讨论一下把
我看好象没有这种流固耦合中辐射的例子
ansys的流畅动力学资料里边只是提到了,
"""因为辐射能穿透流体区域并照射到固体上,所以可以在流固交界面和模型的外边界上施加面载荷,但没有说具体怎么做,
是不是我辐射载荷加的有问题呀
请高手门不吝赐教.
  
感谢感谢

fossiler

请大家多讨论讨论呀，
辐射与对流可以同时加吗

ainite

好例子，谢谢楼主分享，大家多参与讨论定有收获。
  
对于楼主的问题，我印象中似乎在不同的地方是可以分别加对流和辐射，但同一个边界上只能加其中之一。

fossiler

那位高手指点一下呀

fossiler

问题解决了，

tripstones

老大，问题怎么解决的，麻烦说一下
让我们这伙虾虾学习学习
你这样做是不是不用耦合空气和铜的相邻节点？
我现在也在做这方面的东西
辐射也加不上去
急死了

xielihui

问题怎么解决的？

flchflchflch

请问一下，辐射边界条件怎么加载的，以前我做的仅仅是对流传热，加上对流系数就可以了，现在用辐射传热条件但不知道怎么加载。假如是两维的，听说还要什么表面效应单元，请问怎么实现，如果有这方面的事例或教程请给个，由于本人刚接触这个，不懂命令\

xushfu008

问题怎么解决的？

crazy_leo

大家看好这个贴子的时间,如果相差太多,就不要抱太多希望了,如果二维的,且需要做对空气的辐射可以选用surf151,自己看看一个相关的例子就很容易解决的,如果是三维的就选相应的单元,感觉"热分析教程与实例解析"上就有这样的例子.

lgyd

做面-面辐射的可以采用radiosity求解器，可以不用表面效应单元。关于辐射的内容建议专门看一下相关教程，自己瞎摸索是很浪费时间的。

yuguoxingqin

fossiler 怎么解决的麻烦说下，我做的问题与你的相似，也碰到类似问题。QQ：327574871

yuguoxingqin

楼主的辐射是怎么加的啊？急！！