圆管自然对流计算和模拟.docx
- 文档编号:6038810
- 上传时间:2023-05-09
- 格式:DOCX
- 页数:22
- 大小:554.36KB
圆管自然对流计算和模拟.docx
《圆管自然对流计算和模拟.docx》由会员分享,可在线阅读,更多相关《圆管自然对流计算和模拟.docx(22页珍藏版)》请在冰点文库上搜索。
圆管自然对流计算和模拟
水平管和竖直管自然对流计算汇总
1.计算工况表
计算结果
温度工况
100°C
150C
200C
250C
300C
传热系数
w/(m2
h
水平管
7.958
9.115
10.045
10.803
11.527
K)
竖直管
4.715
5.369
5.899
6.335
6.754
换热量
*
水平管
75.962
141.388
215.734
296.472
385.128
W
竖直管
45.008
83.390
126.703
173.860
225.649
最大速度
Umax
水平管
0.476
0.537
0.585
0.697
0.736
m/s
竖直管
0.840
1.050
1.180
1.290
1.390
2.
变化曲线图
worisor200r250r300=c
確温(戈)
换热量(W)
oooooooo
IMIII.IMIII.
ood」50“cMood25CTC309C
餐VC)
以大速皮(m/s)
o
bo
1—^s
-j-JIHH「sod」50dMood250dwood
NraBa(n)
圆管自然对流的计算和数值模拟
已知条件如图1所示:
将一圆管分别水平放置和垂直放置在大空间中进行自
然对流换热,圆管外径D=38mm,长度L=1000mm,空气温度=20C,恒壁
温条件Tw=100,150,200,250,300C,求解自然对流换热系数和换热量以及对流换热时的空气最大速度。
i
A\
U.J
一、数值计算
1.自然对流换热系数和换热量的计算
1)圆管水平放置计算
以壁温Tw-100C为例,计算过程如下:
特征长度:
D=0.038m;
定性温度tm=tw匕2二1002060C;
查空气物性:
■=0.029W:
mK;>=20.110-6m2s;Pr=0.696
空气的体积膨胀系数:
=1tm-2731=1;60•2731=1:
333K'
格拉晓夫数Gr:
g:
\tw7:
D39.81/333100-200.03835
Gr2=62一3.210
v2(20.M10)
大空间自然对流的实验关联式为:
丿Ln
Nu=CGrPr(1-1)
根据计算的格拉晓夫数Gr选择合适的常数C和n(表1):
表1式(1-1)中的常数C和n
加热表面形
流动情况示
流态
系数C和指数n
Gr数适用范围
状与位置
意图
C
n
层流
0.48
1/4
"48
10~5.76X10
横圆管
[&■
过渡流
0.0445
0.37
5.76X108~4.65X109
湍流
0.10
1/3
>4.65X109
由式(1-1)和表1可得:
1/4
Nu=0.48GrPr
51/4
=0.483.21050.696=10.428
:
』=hAtw-t二严7.9583.140.0381100-20i=75.962W
按照以上相同的步骤,在给定恒壁温100,150,200,250,300C的情况下,可以
计算出相应的自然对流的换热系数和换热量,计算结果列于表2中:
表2水平管计算工况表
温度工况
150C
200C
250C
300C
计算参数
100C
hWm2K
7.958
9.115
10.045
10.803
11.527
W
75.962
141.388
215.734
296.472
385.128
2)圆管垂直放置计算
以壁温Tw=100C为例,计算过程如下:
特征长度:
H=1m
定性温度垢二twt:
:
2二1002060C;
查空气物性:
’=0.029WmK;、.=20.110-6m2.s;Pr=0.696
空气的体积膨胀系数:
=1tm-273=160•273=1;333K‘
格拉晓夫数Gr:
厲‘八严L3=9.81/333100:
20「5.83109
v2(20.1"0冷2
对于竖圆柱按照竖壁同用一个关联式必须满足:
d、35
HGrH/4
经验算,并不满足情况,应该按照文献【杨世铭•细长竖圆柱外及竖圆管内
自然对流传热】中的关联式进行计算。
表3竖圆柱自然对流关联式
加热表面
流动情况
Nu/RCE.52严冷
、一1
1/4D
0.1 H 竖圆管 由表1可得: 1/4四8=9.59 1 先计算Ra1/4—=(GrPr)/4D=(5.83心09汉0.696) HH 1/4 0.696 N」0.59十0.52〔Ra 0.59+0.52疋(9.595.83汉109 =162.594 0.029162.594 -4.715Wm2 故水平圆管换热量: : yhAtw七[=4.7153.140.0381100-20=45.008W 按照以上相同的步骤,在给定恒壁温100,150,200,250,300C的情况下,可以 计算出相应的自然对流的换热系数和换热量,计算结果列于表2中: W 45.00883.390 126.703173.860225.649 2. 9A('CJ 超. W0- 2M- 1恥 10Q- 50< —i11■1'■■p— iMrisocao&rss&tsoot-・■CtI 水平管(H)和竖直管(V)自然对流换热系数和换热量的对比图形 3.计算结果分析 由图2和图3可知: 1)水平放置的圆管自然对流的换热系数和换热量都明显高于竖直放置的圆管; 2)随着温度的增加,两者换热系数和换热量都逐渐呈线性增长; 3)水平圆管自然对流换热系数相对增加较多。 二、数值模拟 1.水平圆管的数值模拟 1)物理模型 如图4所示,本文采用的物理模型为大空间自然对流,外边界设置为压力出口边界,与大气相通,内边界为高温管道壁面,圆管直径按照实际尺寸设计。 用 ICEM-CED建立的模型长为380mm,宽为380mm,圆管直径38mm,位于中心位 20C空气 2)网格划分 本次模拟的网格为结构化网格,ICEM网格划分需要对物理模型进行分块处理(block),块的划分采用O-block,O—block易于对内边界做网格加密处理,块的划分和网格的生成如图5和图6所示。 图5 4)网格质量和网格无关性验证 经网格无关性验证后,网格质量符合要求,网格划分合理。 5)计算结果与分析 自然对流是由于空气温度差引起的密度差,从而产生浮升力推动空气运动的现象,实质属于可压缩流动。 在Fluent中气体模型采用Boussinesq可以得到比较好的模拟结果。 Boussinesq近似是将动量方程中密度定义为时间的函数,而能量方程中的密度视为常量。 在Fluent中设置好参数和边界条件后,计算结果如下: ①壁温100C模拟结果 1IJT^-DI 1两£1 3z? ? ni 2怪01 j 3.7*e-l>l 鱼I JU卷・D1 2.|le-DI £>J»ci£D 图9速度云图 图10旋涡 图14旋涡 最大速度可在云图中直接读出: 0.476m/s ②壁温150C模拟结果 图13速度云图 最大速度可在云图中直接读出: 0.537m/s ②壁温200C模拟结果 图17速度云图 图18旋涡 最大速度可在速度云图中直接读出: 0.585m/s ②壁温250C模拟结果 2J01-S-E2 UJ 4,-JU-U.i 1: =8MJPJ 七.J*4U 图19温度云图 图20压力云图 r.ZZt-OI 4.9Er«^1l £加皿 图22旋涡 图21速度云图 最大速度可在云图中直接读出: 0.697m/s ②壁温300C模拟结果 €.176102 7恥個 -2.5J=-O2 -424e4JlZ 3AU-0: 站4 图25速度云图图26旋涡 最大速度可在云图中直接读出: 0.736m/s 速度随着温度变化的汇总表: 表5水平管最大速度计算工况表 尿度工况 计算参数. 100°C s_ 150C 200C 250C 300C Umax,m/s 0.476 0.537 0.585 0.697 0.736 结论分析: 1)自然对流换热强弱取决于高温壁面温度与周围流体温度差的大小,温差越大,换热发展越迅速,流动越强烈; 2)随着壁面温度的增加,最大空气流速也在随之增加; 3)在温差的驱动下形成上升流,并在压差作用下上升流两侧形成漩涡。 2.竖直圆管的数值模拟 1)物理模型 如图4所示,本文采用的物理模型为大空间自然对流,由于物理模型左右对 称,故只需模拟其中的一侧即可,同样外边界设置为压力出口边界,与大气相通,内边界为高温管道壁面,圆管直径按照实际尺寸设计。 用ICEM-CED建立的模型 长为2000mm,宽为570mm,圆管直径38mm,位于中心位置。 绝热 绝热 2) 本次模拟的网格为结构化网格,ICEM网格划分需要对物理模型进行分块处 理(block),内边界的网格加密处理,块的划分和网格的生成如图28和图29所 14 17 图28 示。 网格加密 图29 6)网格质量和网格无关性验证 经网格无关性验证后,网格质量符合要求,网格划分合理。 7)计算结果与分析 自然对流是由于空气温度差引起的密度差,从而产生浮升力推动空气运动的 现象,实质属于可压缩流动。 在Fluent中气体模型采用Boussinesq可以得到比 较好的模拟结果。 Boussinesq近似是将动量方程中密度定义为时间的函数,而能量方程中的密度视为常量。 在Fluent中设置好参数和边界条件后,计算结果如下: ①壁温100C模拟结果 4.1kM 1.73eD2 541e+4)£ 图30温度云图 e.eae-os ■3.0^-02 ■6.44^-02 7.83e-tt2 图31压力云图 a.3i^-ci 0.65e-01 493e^1 3.32fl-d1 1,6Ge-D1 00(k+CD 图32速度云图 aaa-%a.-hTn.,i■■-L■• 、0'-le»■V-«a,x-fcnv«、jXf,a、、h・、f-W! 078^-01 -■■JB*■&1**1! ■-b-b"11■! 1'b-b1iWH1-■■! ■''.,V |—牛科}匕(汗,一1、2・7 -补i-i'i1! ■-^-^*iv■■nviIB■-^■ 、-'、、、、、-、 h气Tdl•补F*1讥1UPL*jFITU-iI BB! —.ai-n.-M1/i.M.in! S、、叭・、[BypiH 川uwi八w■si.mm■ii订[1111*7 -、“审・乍*初刘•卄-r、i]刖、“气[和鸟 沙n、wzw,钠啊飞xvn]\门勺1>、 ・・・2"T2电UJIu11、U、I・■"■%・.'骤 ■•十占炸临、沾i沙代、「山他]\门、■*■* '»■・・■1'・tt・l5•••-,»■LIUV11«V-R.R-"""a*JUJPW^i"■JL 图33旋涡 最大速度可在云图中直接读出: 0.831m/s ②壁温150C模拟结果 3.45e443J 2.92*4)2 图34温度云图 1.OSe-kOD S.44e^1 Q.33H31 472e~D1 62^02 1.74e-D2 -L76e-D2 -5.27e-02 -8a77e-Da 图35压力云图 2.11&-D1 O.DOt+OD 图36速度云图 最大速度可在云图中直接读出: ②壁温200C模拟结果 —— ■■■■IBMrS■■■■■ 1.05m/s 图37旋涡 图38温度云图 -O.^-DZ 图40速度云图 图39压力云图 图41旋涡 最大速度可在速度云图中直接读出: 1.18m/s ②壁温250C模拟结果 -6 5*4-01 图43压力云图 7-72 图45旋涡 1LW •.w....•. 二: 薦站;: 『: ;「;,; 图42温度云图 图44速度云图 最大速度可在云图中直接读出: 1.29m/s②壁温300C模拟结果 图46温度云图 725C-D2 -S.flDte-ftZ U^h-UI -l.7^=-OI 图47压力云图 •.42&+0C : Q站辽: ・"g ・‘亠・•"■fah■■ 1%U'hfciriV'fc“I" 3-1 图48速度云图图49旋涡 最大速度可在云图中直接读出: 1.39m/s 速度随着温度变化的汇总表: 表6竖直管最大速度计算工况表 温度工况 计算参数. 100°C 150C 200C 250C 300C Umax,m/s 0.84 1.05 1.18 1.29 1.39 结论分析: 1)自然对流换热强弱取决于咼温壁面温度与周围流体温度差的大小,温差越大,换热发展越迅速,流动越强烈; 2)随着壁面温度的增加,最大空气流速也在随之增加; 3)在温差的驱动下形成上升流,并在压差作用下上升流两侧形成漩涡。 3.水平圆管和竖直圆管自然对流的最大速度对比 1)现将模拟的最大速度汇总,如表格7: 表7最大速度对比表格 温度工况 计算参数 100C 150C 200C 250C 300C 水平管 Umax,m/s 0.476 0.537 0.585 0.697 0.736 垂直管 Umax,m/s 0.84 1.05 1.18 1.29 1.39 2)水平圆管和竖直圆管自然对流的最大速度曲线图,如图50: ■水平管 —竖直菅 ].4T 1.2 1.0 Q-4r|1I1I1IvI1 100T! 150r2D0C25CTC30CTC 内壁消M o.a 0.6 图50 结论: ①竖直管的自然对流最大速度明显高于水平管自然对流的最大的速度; ②随着温度的增加,两者的最大速度都逐渐呈线性增加,且增加的幅 度越来越小; 1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 圆管 自然 对流 计算 模拟
![提示](https://static.bingdoc.com/images/bang_tan.gif)