辐射下加热腔内自然对流传热特性研究.pdf
《辐射下加热腔内自然对流传热特性研究.pdf》由会员分享,可在线阅读,更多相关《辐射下加热腔内自然对流传热特性研究.pdf(11页珍藏版)》请在文库网上搜索。
1、Modeling and Simulation 建模与仿真建模与仿真,2023,12(5),4605-4615 Published Online September 2023 in Hans.https:/www.hanspub.org/journal/mos https:/doi.org/10.12677/mos.2023.125420 文章引用文章引用:葛忠义,周瑞睿.辐射下加热腔内自然对流传热特性研究J.建模与仿真,2023,12(5):4605-4615.DOI:10.12677/mos.2023.125420 辐射下加热腔内自然对流传热特性研究辐射下加热腔内自然对流传热特性研究 葛忠
2、义葛忠义,周瑞睿周瑞睿*上海理工大学能源与动力工程学院,上海 收稿日期:2023年7月28日;录用日期:2023年9月11日;发布日期:2023年9月18日 摘摘 要要 辐射吸收引发的自然对流现象广泛存在于自然界以及工业领域,为了探究辐射引发的自然对流传热特性,辐射吸收引发的自然对流现象广泛存在于自然界以及工业领域,为了探究辐射引发的自然对流传热特性,采用理论分析与数值模拟相结合,对由下方入射辐射加热二维方腔内半透明流体所引发的自然对流进行采用理论分析与数值模拟相结合,对由下方入射辐射加热二维方腔内半透明流体所引发的自然对流进行研究。结果表明:同一瑞利数下,流体失稳时间随光学厚度研究。结果表明
3、:同一瑞利数下,流体失稳时间随光学厚度 的增加而减小,当的增加而减小,当20=时,流体失稳时时,流体失稳时间最短;通过尺度分析推导出当光学厚度间最短;通过尺度分析推导出当光学厚度0时,热边界层厚度与瑞利数及光学厚度之间满足标度率:时,热边界层厚度与瑞利数及光学厚度之间满足标度率:()IRaPr0.20.20.211+,并将数值模拟结果与之对比,并将数值模拟结果与之对比,IRa数的标度指数基本吻合;同时发现,不同数的标度指数基本吻合;同时发现,不同光学厚度下,光学厚度下,Nu数随着数随着IRa数的增加而逐渐增大,随着光学厚度的增加,数的增加而逐渐增大,随着光学厚度的增加,INuRa 中的标度指数
4、中的标度指数 反反而减小;最后采用数值模拟方法,得到了而减小;最后采用数值模拟方法,得到了Nu数与光学厚度及数与光学厚度及IRa数之间的标度关系:数之间的标度关系:INuRa0.9180.302 ,并利用尺度分析方法推导出光学厚度并利用尺度分析方法推导出光学厚度0时,时,Nu数与光学厚度数与光学厚度 满足:满足:Nu 14 。关键词关键词 辐射加热,光学厚度,自然对流,热边界层,传热标度率辐射加热,光学厚度,自然对流,热边界层,传热标度率 Characterization of Natural Convection Heat Transfer in a Heated Cavity under
5、Radiation Zhongyi Ge,Ruirui Zhou*School of Energy and Power Engineering,The University of Shanghai for Science and Technology,Shanghai Received:Jul.28th,2023;accepted:Sep.11th,2023;published:Sep.18th,2023 Abstract The phenomenon of natural convection induced by radiation absorption exists widely in
6、nature as well as in industry.In order to investigate the heat transfer characteristics of radiation-induced *通讯作者。葛忠义,周瑞睿 DOI:10.12677/mos.2023.125420 4606 建模与仿真 natural convection,a combination of theoretical analyses and numerical simulations is used to investigate the natural convection induced
7、by incident radiation from below heating a two-di-mensional square cavity with a semi-transparent fluid.The results show that at a constant Ray-leigh number,the fluid destabilization time decreases with the increase of the optical thickness,when 20=,with a large fluctuation of the temperature at the
8、 lower wall surface;through the scale analysis,it is deduced that when the optical thickness 0,the scale rate is satisfied be-tween the temperature boundary layer thickness and the Rayleigh number and the optical thick-ness:()IRaPr0.20.20.211+,and the numerical simulation results are compared with i
9、t,and the scale exponent of the Rayleigh number IRa is basically coincident.It is also found that the Nu number gradually increases with increasing IRa number for different optical thicknesses,and the scaling index in INuRa decreases instead with increasing optical thickness.By using the numerical s
10、imulation method,the scaling relationship between the Nu number and the optical thickness and IRa number is obtained:INuRa0.9180.302 ,and the scaling analysis method is used to deduce that the Nu number and the optical thickness satisfy:Nu14 when the optical thickness 0.Keywords Radiation Heating,Op
11、tical Thickness,Natural Convection,Thermal Boundary Layer,Heat Transfer Scale Rate Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 自然对流是自然界和各类工程问题中的常见物理现象
12、,与人们的日常生活息息相关。其不仅存在于冰川湖泊的融化1、大气及海洋环流2 3、地幔对流板块运动4等自然现象中;还广泛应用于核反应堆冷却5、芯片散热6、磁流体传热和新型储能材料7 8等工业问题中。而 Rayleigh-Bnard(RB)对流系统9,则是从众多复杂的自然对流现象以及工程问题中简化出来的研究对流问题的经典模型。由于 RB系统中存在多种热量传递方式,因此它也是传热领域经典的课题,研究 RB 系统对实际情况下各类对流过程中的传热机理有重要的意义。由于经典的RB对流模型其上下边界温度恒定,这与一些实际系统的热边界条件并不相符。如太阳辐射在地面上反射对大气进行加热现象10、水冷壁受辐射热升
13、温11等由辐射内热源驱动的热对流问题。在辐射驱动的热对流问题中,众多学者对其开展了模拟、理论和实验等方面研究。在理论分析方面,Lei等12对三角封闭腔进行标度分析,根据瑞利数Ra、底部三角倾斜角tan和临界瑞利数cRa对流动状态进行了描述。Mao等13对小斜坡水库吸收辐射引发的自然对流进行了尺度分析研究,获得瑞利数相对于水平位置的两个临界函数,并根据底部坡度和最大水深对流体流动区域进行划分,建立位置相关标度对不同区域的流量进行说明。Hattori等14通过线性理论与直接稳定性分析,对辐射驱动的自然对流进行了探究,并结合水平通流和横向约束,研究了两种不同形式的流动不稳定性,并对每种不稳定性模式形
14、成的本征问题进行数值求解,获得了本征函数与失稳条件。在模拟方面:Miquel等15研究了努塞尔数NuOpen AccessOpen Access葛忠义,周瑞睿 DOI:10.12677/mos.2023.125420 4607 建模与仿真 与普朗特数Pr之间的依赖关系,对于无应力下表面,传热标度率满足1 2NuPr,而对于无滑移下表面,可能呈现两种传热标度率,1 2NuPr与1 6NuPr。Hattori等16采用理论分析与模拟验证相结合的方法,对上方直接入射辐射所形成的稳定热分层以及由底部边界热通量所产生的羽流进行了相关研究,并将标度关系与数值模拟进行了对比。在实验方面:Bouillaut等
15、17搭建了辐射对流实验装置,通过改变光学厚度,对腔体内流体流动机理进行研究,结果表明:随着光学厚度的变化,TNuRa的指数可在1/3至1/2范围内变动,该标度指数范围与RB系统所得到的范围相一致,标度指数由1/3转变为1/2的临界点为36TRaPr与12NuPr。但仍有部分学者认为该标度区间并不适用于含有内热源的物理问题。Klein等18通过ODT模型对文献17进行更大范围的数值模拟,发现其标度指数上界限为0.55,并非是Bouillaut等所得到的标度指数。以上学者均对辐射所引发的自然对流进行了研究。但在其研究过程中,未能充分考虑光学厚度对传热特性的影响。因此,本文将采用理论分析与数值模拟相
16、结合的方式,探讨不同光学厚度、不同瑞利数IRa下辐射引发的自然对流传热特性。接下来第 4.1 节主要探讨光学厚度对流体失稳时间的影响,采用模拟方法对失稳时间进行探究。第 4.2 节简要分析光学厚度对温度边界层厚度的影响。第 4.3 节讨论努塞尔数、光学厚度以及瑞利数三者之间的标度指数关系。最后就辐射对自然对流的影响进行总结。2.问题描述及控制方程问题描述及控制方程 2.1.物理模型物理模型 考虑如图 1 所示几何模型。方腔左、右边界为周期边界,下表面为透明、无滑移绝热表面,上表面为透明、无应力绝热表面,内部为半透明流体。辐射从下方垂直均匀射入方腔内,流体吸收辐射后温度升高,形成密度差,近而在浮
17、力的作用下发生自然对流。Figure 1.Geometry model 图图 1.几何模型 2.2.数学模型数学模型 对于二维方腔区域,辐射从下方垂直均匀射入,辐射能随着流体吸收的深度逐渐衰弱,满足比尔朗伯定律:()0exp,0QIyy=(1)其中,为流体的辐射吸收系数,0I为入射辐射强度。下方垂直辐射加热效果等效于一个内热源的作用,葛忠义,周瑞睿 DOI:10.12677/mos.2023.125420 4608 建模与仿真 其引发的自然流动现象可由 Navier-Stokes 方程组和能量方程描述。其中,浮升力的影响采用 Boussinesq近似描述:0uvxy+=(2)22221uuvp
18、uuuvtxyxxy+=+(3)()22221vvvpvvuvg TTtxyyxy+=+(4)22022eypITTTTTuvatxycxy+=+(5)其中,为流体密度,为流体运动粘度,为热膨胀系数,pc为流体比热容,T为流体平均温度,p为流体压力,a为流体热扩散系数。对应边界条件为:000,xx Lxx Lxx LuuvvTT=(6)0,0,0Tuvyy=(7)0,0,uvTyHxyy=(8)上述系统中,流动在一定时间后可进入准稳态阶段,而系统平均温度由于壁面绝热随时间线性增加:()000000e1 ed dLHypHpy xtHIcITTTtcLH+=+=(9)其中:0T为系统初始温度。采
19、用以下标尺对控制方程进行无量纲化:()*2200,ppTTc apacxytavHuHxytvupHHaaI HHH g I=,(10)以及无量纲控制参数,瑞利数IRa、普朗特数Pr以及光学厚度 420=,Ipg HRa,PHc aIra=(11)代入原方程,整理后略去变量上标星号,可得无量纲控制方程组为:0uvxy+=(12)2222IuuupuuuvRa PrPrtxyxxy+=+(13)2222IIvvvpvvuvRa PrPrRa Prtxyyxy+=+(14)葛忠义,周瑞睿 DOI:10.12677/mos.2023.125420 4609 建模与仿真 ()2222e1 eyuvtx
20、yxy+=+(15)记LH=,则相应无量纲边界条件如下:000,xxxxxxuuvv=(16)0,0,0uvyy=(17)0,0,0,1uvyxy=(18)易知对于任意时刻均有100d d0 x y=。3.数值求解方法与无关性验证数值求解方法与无关性验证 3.1.数值求解方法数值求解方法 采用有限体积法对控制方程(12)(15)进行离散,基于SIMPLE算法处理压力与速度耦合,扩散项采用中心差分格式离散,对流项采用QUICK格式离散,时间项应用二阶隐式迎风的离散方式。采用的无量纲控制参数见表1。Table 1.Dimensionless parameter 表表 1.无量纲参数 Pr Ra L
21、/H 7 107/108/109/1010 4 3.2.网格无关性与时间独立性验证网格无关性与时间独立性验证 分别采用 4 种网格尺寸 400 130、640 185、880 246、1120 310 进行网格无关性验证,选取瑞利数1010IRa=,通过时间空间平均的努塞尔数 Nu 进行网格无关性与时间独立性检验。详细信息如表 2 所示。1Nu=(19)其中,102yy=为下壁面与中心位置的时均温差,表示空间平均温度,表示时间平均 温差。x方向网格均匀布置,y方向网格在边界层处进行加密,确保边界层内部至少存在15个网格。流体光学厚度取为无穷大,即此时流体不透明,辐射加热效果等效于下壁面恒热流加
22、热。Table 2.Grid information 表表 2.网格信息 网格 第一层网格尺寸 节点 Nx Ny Nu 相对误差%1 9.4 104 400130 3.23%2 7.4 104 640185 0.98%3 5.8 104 880246 1.03%4 4.7 104 1120310 葛忠义,周瑞睿 DOI:10.12677/mos.2023.125420 4610 建模与仿真 (a)(b)Figure 2.(a)Comparison plot of Nu number for different grid sizes(b)Comparison plot of Nu number
23、for different time steps 图图 2.(a)不同网格尺寸下的 Nu 数对比图(b)不同时间步长下的 Nu 数对比图 图 2(a)为四组网格系统在准稳态阶段 Nu 数随时间的发展。其中,通过计算可以得到,网格 1、网格2 以及网格 3 与网格 4 的 Nu 误差相差分别为 3.23%、0.98%、1.03%,因此为了保证精度以及节省计算成本,本文所有工况选取网格 2 进行模拟计算。图 2(b)为三组时间步长在准稳态阶段 Nu 数随时间的发展。时间步长为 5.75 106,与时间步长为 1.44 106相比,Nu 误差相差 5.83%,时间步长为 2.87 106与时间步长为
24、 1.44 106相比,Nu 误差相差 0.43%,本文选取时间步长为 2.87 106进行计算可以更好的节约时间成本。3.3.程序有效性验证程序有效性验证 文献17采用实验的方式,测得多组工况下时间空间平均努塞尔数 Nu 与瑞利数TRa的标度关系,本文对照文献17中标度关系对程序进行有效验证。选取其中20 Pr7=、时的工况,文献中努塞尔数 Nu与瑞利数TRa的标度指数关系近似满足0.5TNuRa,图 3 给出了本文努塞尔数 Nu 与瑞利数TRa的标度指数与文献中的标度指数结果对比图。对比结果表明,当前研究所用到的程序准确可靠。Figure 3.Scaling relationship be
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 辐射 加热 自然 对流 传热 特性 研究