辐射下加热腔内自然对流传热特性研究.pdf

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 辐射下加热腔内自然对流传热特性研究辐射下加热腔内自然对流传热特性研究 葛忠义葛忠义，周瑞睿周瑞睿*上海理工大学能源与动力工程学院，上海 收稿日期：2023年7月28日；录用日期：2023年9月11日；发布日期：2023年9月18日 摘摘 要要 辐射吸收引发的自然对流现象广泛存在于自然界以及工业领域，为了探究辐射引发的自然对流传热特性，辐射吸收引发的自然对流现象广泛存在于自然界以及工业领域，为了探究辐射引发的自然对流传热特性，采用理论分析与数值模拟相结合，对由下方入射辐射加热二维方腔内半透明流体所引发的自然对流进行采用理论分析与数值模拟相结合，对由下方入射辐射加热二维方腔内半透明流体所引发的自然对流进行研究。结果表明：同一瑞利数下，流体失稳时间随光学厚度研究。结果表明：同一瑞利数下，流体失稳时间随光学厚度 的增加而减小，当的增加而减小，当20=时，流体失稳时时，流体失稳时间最短；通过尺度分析推导出当光学厚度间最短；通过尺度分析推导出当光学厚度0时，热边界层厚度与瑞利数及光学厚度之间满足标度率：时，热边界层厚度与瑞利数及光学厚度之间满足标度率：()IRaPr0.20.20.211+，并将数值模拟结果与之对比，并将数值模拟结果与之对比，IRa数的标度指数基本吻合；同时发现，不同数的标度指数基本吻合；同时发现，不同光学厚度下，光学厚度下，Nu数随着数随着IRa数的增加而逐渐增大，随着光学厚度的增加，数的增加而逐渐增大，随着光学厚度的增加，INuRa 中的标度指数中的标度指数 反反而减小；最后采用数值模拟方法，得到了而减小；最后采用数值模拟方法，得到了Nu数与光学厚度及数与光学厚度及IRa数之间的标度关系：数之间的标度关系：INuRa0.9180.302 ，并利用尺度分析方法推导出光学厚度并利用尺度分析方法推导出光学厚度0时，时，Nu数与光学厚度数与光学厚度 满足：满足：Nu 14 。关键词关键词 辐射加热，光学厚度，自然对流，热边界层，传热标度率辐射加热，光学厚度，自然对流，热边界层，传热标度率 Characterization of Natural Convection Heat Transfer in a Heated Cavity under 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 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 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 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 it,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 simulation 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,Optical 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.引言引言 自然对流是自然界和各类工程问题中的常见物理现象，与人们的日常生活息息相关。其不仅存在于冰川湖泊的融化1、大气及海洋环流2 3、地幔对流板块运动4等自然现象中；还广泛应用于核反应堆冷却5、芯片散热6、磁流体传热和新型储能材料7 8等工业问题中。而 Rayleigh-Bnard(RB)对流系统9，则是从众多复杂的自然对流现象以及工程问题中简化出来的研究对流问题的经典模型。由于 RB系统中存在多种热量传递方式，因此它也是传热领域经典的课题，研究 RB 系统对实际情况下各类对流过程中的传热机理有重要的意义。由于经典的RB对流模型其上下边界温度恒定，这与一些实际系统的热边界条件并不相符。如太阳辐射在地面上反射对大气进行加热现象10、水冷壁受辐射热升温11等由辐射内热源驱动的热对流问题。在辐射驱动的热对流问题中，众多学者对其开展了模拟、理论和实验等方面研究。在理论分析方面，Lei等12对三角封闭腔进行标度分析，根据瑞利数Ra、底部三角倾斜角tan和临界瑞利数cRa对流动状态进行了描述。Mao等13对小斜坡水库吸收辐射引发的自然对流进行了尺度分析研究，获得瑞利数相对于水平位置的两个临界函数，并根据底部坡度和最大水深对流体流动区域进行划分，建立位置相关标度对不同区域的流量进行说明。Hattori等14通过线性理论与直接稳定性分析，对辐射驱动的自然对流进行了探究，并结合水平通流和横向约束，研究了两种不同形式的流动不稳定性，并对每种不稳定性模式形成的本征问题进行数值求解，获得了本征函数与失稳条件。在模拟方面：Miquel等15研究了努塞尔数NuOpen AccessOpen Access葛忠义，周瑞睿 DOI:10.12677/mos.2023.125420 4607 建模与仿真 与普朗特数Pr之间的依赖关系，对于无应力下表面，传热标度率满足1 2NuPr，而对于无滑移下表面，可能呈现两种传热标度率，1 2NuPr与1 6NuPr。Hattori等16采用理论分析与模拟验证相结合的方法，对上方直接入射辐射所形成的稳定热分层以及由底部边界热通量所产生的羽流进行了相关研究，并将标度关系与数值模拟进行了对比。在实验方面：Bouillaut等17搭建了辐射对流实验装置，通过改变光学厚度，对腔体内流体流动机理进行研究，结果表明：随着光学厚度的变化，TNuRa的指数可在1/3至1/2范围内变动，该标度指数范围与RB系统所得到的范围相一致，标度指数由1/3转变为1/2的临界点为36TRaPr与12NuPr。但仍有部分学者认为该标度区间并不适用于含有内热源的物理问题。Klein等18通过ODT模型对文献17进行更大范围的数值模拟，发现其标度指数上界限为0.55，并非是Bouillaut等所得到的标度指数。以上学者均对辐射所引发的自然对流进行了研究。但在其研究过程中，未能充分考虑光学厚度对传热特性的影响。因此，本文将采用理论分析与数值模拟相结合的方式，探讨不同光学厚度、不同瑞利数IRa下辐射引发的自然对流传热特性。接下来第 4.1 节主要探讨光学厚度对流体失稳时间的影响，采用模拟方法对失稳时间进行探究。第 4.2 节简要分析光学厚度对温度边界层厚度的影响。第 4.3 节讨论努塞尔数、光学厚度以及瑞利数三者之间的标度指数关系。最后就辐射对自然对流的影响进行总结。2.问题描述及控制方程问题描述及控制方程 2.1.物理模型物理模型 考虑如图 1 所示几何模型。方腔左、右边界为周期边界，下表面为透明、无滑移绝热表面，上表面为透明、无应力绝热表面，内部为半透明流体。辐射从下方垂直均匀射入方腔内，流体吸收辐射后温度升高，形成密度差，近而在浮力的作用下发生自然对流。Figure 1.Geometry model 图图 1.几何模型 2.2.数学模型数学模型 对于二维方腔区域，辐射从下方垂直均匀射入，辐射能随着流体吸收的深度逐渐衰弱，满足比尔朗伯定律：()0exp,0QIyy=(1)其中，为流体的辐射吸收系数，0I为入射辐射强度。下方垂直辐射加热效果等效于一个内热源的作用，葛忠义，周瑞睿 DOI:10.12677/mos.2023.125420 4608 建模与仿真 其引发的自然流动现象可由 Navier-Stokes 方程组和能量方程描述。其中，浮升力的影响采用 Boussinesq近似描述：0uvxy+=(2)22221uuvpuuuvtxyxxy+=+(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为系统初始温度。采用以下标尺对控制方程进行无量纲化：()*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 eyuvtxyxy+=+(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/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个网格。流体光学厚度取为无穷大，即此时流体不透明，辐射加热效果等效于下壁面恒热流加热。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 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与时间步长为 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 between Nusselt number and Rayleigh number(=20)图图 3.(=20)时努塞尔数 Nu 与瑞利数 RaT的标度关系 葛忠义，周瑞睿 DOI:10.12677/mos.2023.125420 4611 建模与仿真 4.计算结果与分析计算结果与分析 4.1.光学厚度对流动失稳时间影响光学厚度对流动失稳时间影响 流体流动可分为，初始阶段、过渡阶段、准稳态阶段三个部分。在初始阶段，流体速度较低，此时能量主要通过导热方式传递。由于腔体下方辐射的作用，下壁面温度随着时间逐渐升高，当温度边界层温差达到一定值时，下表面将产生向上的热羽流，此时流体热边界层被破坏，流体将从初始阶段转为过渡阶段。流体失稳需要满足热边界层温差超过一定范围，即当地瑞利数LRa超过临界值TcRa。3LIbTcRaRaRa=(20)其中，b为热边界层温差，为热边界层厚度。对于小光学厚度引发的对流，失稳临界瑞利数影响较小，并且由公式(15)可得：t (21)由公式(20)与(21)得：32TcLIbRaRatRa=(22)(a)(b)(c)(d)Figure 4.Change of critical time under different Rayleigh number 图图 4.不同瑞利数下临界时间变化 葛忠义，周瑞睿 DOI:10.12677/mos.2023.125420 4612 建模与仿真 图 4 给出了不同瑞利数、光学厚度下，腔体下表面流体平均温度随时间的变化情况。可以看出，在同一瑞利数下，随着光学厚度的增加，临界失稳时间逐渐减小。其原因在于，当光学厚度较大时，大部分辐射能被边界层吸收，导致边界层温度升高较快，并以热羽流形式带走边界层能量，导致流体达到失稳状态。而当光学厚度较小时，大部分辐射能将透过边界层对腔体内部流体进行加热，此时边界层内部温升较慢，流体达到失稳时间较长。同时，该结论可由公式(22)得到验证。表 3 给出了不同工况下的具体失稳时间。可以看出，瑞利数710IRa=时的失稳时间相比于其他工况更长。这是由于低瑞利数下流体的粘性力较强，流动强度较弱，因此流动失稳时间较长。Table 3.Change of critical instability time under different Rayleigh number and optical thickness 表表 3.不同瑞利数、光学厚度下失稳临界时间变化 瑞利数 RaI 光学厚度 失稳时间 tcr 107 0.5 0.0582 1 0.0349 10 0.0125 20 0.0107 108 0.5 0.0201 1 0.0132 10 0.005 20 0.0038 109 0.5 0.0027 1 0.0018 10 0.0006 20 0.0004 1010 0.5 0.0027 1 0.0018 10 0.0007 20 0.0005 4.2.光学厚度对热边界层厚度的影响光学厚度对热边界层厚度的影响 早期Chill等19人通过实验的方式得到681010TRa，4Pr=范围内热流边界层的厚度，得出结论，温度边界层内部的温度基本满足理论的线性变化。因此在提取边界层方法中，通常采用斜率法来取得温度边界层的厚度，温度场取时间空间平均，在下壁面处可以得到近似直线的温度分布。同时，在腔体中心2yH=处与下壁面温度边界层之间，温度分布同样可以近似的看成线性变化，两条切线的交点到下壁面的距离可被看作温度边界层厚度。图5给出了不同光学厚度下，瑞利数IRa与热边界层厚度之间的关系，可以看出随着瑞利数IRa的增加，边界层厚度逐渐减小。结论式(23)与分析推导结果式(27)可以验证这一结论。与此同时，发现不同光学厚度下，瑞利数IRa的标度指数几乎一致。为了验证其标度率的正确性，下文将采用理论分析方法对其进行验证。为方便模拟与理论分析结果对比，模拟结果将取平均标度指数进行标度探究，如公式(23)所示：0.216IRa (23)葛忠义，周瑞睿 DOI:10.12677/mos.2023.125420 4613 建模与仿真 Figure 5.Relation between Rayleigh number and thickness scal-ing rate of thermal Boundary layer 图图 5.瑞利数 RaI与热边界层厚度 标度率关系 流体在准稳态阶段，由动量方程(14)，热流体脱离热边界层时惯性力与粘性力及浮力之间达到平衡：22IbvvPrRa Pr+(24)其中：b为热边界层温差。又由能量方程(15)，对流项、扩散项以及源项三者之间达到平衡：()2e1 ebbv (25)若0成立，则式(25)可进一步简化：2bbv (26)由式(24)、(26)可得：()0.20.20.211IRaPr+(27)由模拟方法得到的结论式(23)与尺度分析推导的结果式(27)相比可得，标度指数量级基本相当，进一步说明数值模拟结果的准确性，同时也展现了瑞利数IRa对流体流动有重要影响。4.3.光学厚度对标度指数的影响光学厚度对标度指数的影响 在湍流流动中，热量如何运输是研究热对流的重点关注问题，即努塞尔数与瑞利数之间的变化规律。对于有辐射驱动的对流系统，热量运输也与光学厚度有紧密的关系。因此，本节将主要探讨光学厚度对努塞尔数的影响，试图找到小光学厚度下INuRa的表达式。图 6(a)中给出了由模拟得到不同光学厚度下，Nu 数随IRa数的变化关系，随着IRa数的增加，Nu 数逐渐增大。随着光学厚度的增加，INuRa中的标度指数反而减小。在0.520范围，有0.3020.341。从图 6(b)可看出，以20=时结果为基准做归一化，可以得到关系式：0.9180.9180.30220INuNuRa (28)葛忠义，周瑞睿 DOI:10.12677/mos.2023.125420 4614 建模与仿真 (a)(b)Figure 6.(a)Relationship between Rayleigh number and Nusselt number under different optical thicknesses(b)Relationship between Nusselt number and optical thickness 图图 6.(a)不同光学厚度下瑞利数 RaI与努塞尔数 Nu 的变化关系(b)努塞尔数 Nu 与光学厚度 之间关系 实际上根据式(15)，可对努塞尔数 Nu 与光学厚度之间的关系Nu做进一步理论推导，由式(15)可得：()2224e1 evHH (29)当0时，式(29)进一步可以得到：4 (30)根据式(19)，有：14Nu (31)由理论推导结果式(31)与模拟结果式(28)对比可以看出，光学厚度的指数几乎相当，也进一步说明模拟结果的准确性。5.结果与讨论结果与讨论 本文对下方垂直入射辐射引发的自然对流问题进行了相关研究，主要探讨不同瑞利数IRa(7101010IRa)与不同光学厚度(0.520)对流体传热特性的影响，得到如下结论：1)流体在初始阶段，恒定瑞利数下，流体失稳时间随着光学厚度的增加而减小，当20=时，下壁面温度波动最大，失稳时间最小。当瑞利数710IRa=时，流体失稳时间相比与其他工况下时间更长。2)通过尺度分析方法，推导出了光学厚度0时温度边界层厚度与瑞利数IRa之间的标度关系，两者之间满足：()0.20.20.211IRaPr+。同时将模拟结果(0.216IRa)与尺度分析结果进行对比，结果 吻合较好，验证了数值模拟的准确性。3)通过数值模拟方法得到了光学厚度0时，努塞尔数 Nu 与瑞利数IRa以及光学厚度三者之间标度关系：0.9180.302INuRa。并采用尺度分析方法推导出努塞尔数Nu与光学厚度之间满足：14Nu。两种结果中光学厚度的标度指数几乎相当。同时，由模拟得到不同光学厚度下，随着 Nu 数的增加，IRa葛忠义，周瑞睿 DOI:10.12677/mos.2023.125420 4615 建模与仿真 数逐渐增大。随着光学厚度的增加，INuRa中的标度指数反而减小。基金项目基金项目 国家自然科学基金青年基金(52206103)、上海市青年科技英才扬帆计划(21YF1430400)。参考文献参考文献 1 Polashenski,C.,Perovich,D.and Courville,Z.(2012)The Mechanisms of Sea Ice Melt Pond Formation and Evolution.Journal of Geophysical Research Oceans,117,C01001.https:/doi.org/10.1029/2011JC007231 2 Mei,S.J.and Yuan,C.(2022)Urban Buoyancy-Driven Air Flow and Modelling Method:A Critical Review.Building and Environment,210,Article ID:108708.https:/doi.org/10.1016/j.buildenv.2021.108708 3 Riehl,H.,Yeh,T.C.H.,Malkus,J.S.and la Seur,N.E.(2010)The North-East Trade of the Pacific Ocean.Quarterly Journal of the Royal Meteorological Society,77,598-626.https:/doi.org/10.1002/qj.49707733405 4 Henshall,T.(2019)Mantle Convection in Super-Earths.Nature Reviews Physics,1,181.https:/doi.org/10.1038/s42254-019-0035-5 5 Sato,T.,Maeda,K.,Nagataki,S.,et al.(2021)High-Entropy Ejecta Plumes in Cassiopeia A from Neutrino-Driven Convection.Nature,592,537-540.https:/doi.org/10.1038/s41586-021-03391-9 6 Myyhren,J.A.and Holmberg,S.(2011)Improving the Thermal Performance of Ventilation RadiatorsThe Role of Internal Convection Fins.International Journal of Thermal Sciences,50,115-123.https:/doi.org/10.1016/j.ijthermalsci.2010.10.011 7 Li,M.G.,Zheng,C.,Zhao,Q.,Chen,X.and Wu,W.T.(2021)Anisotropic Heat Transfer of Ferro-Nanofluid in Partially Heated Rectangular Enclosures under Magnetic Field.Case Studies in Thermal Engineering,26,Article ID:101145.https:/doi.org/10.1016/j.csite.2021.101145 8 Zhang,S.,Feng,D.L.,Shi,L.,et al.(2021)A Review of Phase Change Heat Transfer in Shape-Stabilized Phase Change Materials(ss-PCMs)Based on Porous Supports for Thermal Energy Storage.Renewable and Sustainable Energy Reviews,135,Article ID:110127.https:/doi.org/10.1016/j.rser.2020.110127 9 周全,夏克青.Rayleigh-Bnard 湍流热对流研究的进展、现状及展望J.力学进展,2012,42(3):231-251.10 Ingersoll,P.(1996)Natural Convection as a Heat Engine:A Theory for CAPE.Journal of the Atmospheric Sciences,53,572-585.11 Mies,J.J.,Hammaker,R.G.,Madding,R.P.and Sunderland,J.E.(2023)Analysis of Thermal Radiation in Coal-Fired Furnaces.Proceedings of SPIE,3056,20-32.12 Lei,C.and Patterson,J.C.(2002)Unsteady Natural Convection in a Triangular Enclosure Induced by Absorption of Radiation.Journal of Fluid Mechanics,460,181-209.https:/doi.org/10.1017/S0022112002008091 13 Mao,Y.,Lei,C.and Patterson,J.C.(2009)Unsteady Natural Convection in a Triangular Enclosure Induced by Ab-sorption of Radiation a Revisit by Improved Scaling Analysis.Journal of Fluid Mechanics,622,75-102.https:/doi.org/10.1017/S0022112008005077 14 Hattori,T.,Patterson,J.C.and Lei,C.(2015)On the Stability of Internally Heated Natural Convection Due to the Absorp-tion of Radiation in a Laterally Confined Fluid Layer with a Horizontal Throughflow.International Journal of Heat and Mass Transfer,81,846-861.https:/doi.org/10.1016/j.ijheatmasstransfer.2014.11.010 15 Miquel,B.,Bouillaut,V.,Aumatre,S.and Gallet,B.(2020)On the Role of the Prandtl Number in Convection Driven by Heat Sources and Sinks.Journal of Fluid Mechanics,900,R1.https:/doi.org/10.1017/jfm.2020.485 16 Hattori,T.,Patterson,J.C.and Lei,C.(2015)Scaling and Direct Stability Analyses of Natural Convection Induced by Absorption of Solar Radiation in a Parallelepiped Cavity.International Journal of Thermal,88,19-32.https:/doi.org/10.1016/j.ijthermalsci.2014.09.005 17 Bouillaut,V.,Lepot,S.,Aumatre,S.and Gallet,B.(2019)Transition to the Ultimate Regime in a Radiatively Driven Convection Experiment.Journal of Fluid Mechanics,861,R5.https:/doi.org/10.1017/jfm.2018.972 18 Marten,K.(2021)Transition to the Ultimate Regime in a Stochastic Model for Thermal Convection with Internal Sources.19 Chilla,F.,Ciliberto,S.,Innocenti,C.and Pampaloni,E.(1993)Boundary Layer and Scaling Properties Inturbulent Thermal Convection.Il Nuovo Cimento D,15,1229-1249.https:/doi.org/10.1007/BF02451729