三维方腔内竖直平板自然对流的数值模拟

三维方腔内竖直平板自然对流的数值模拟杨沫，李少华，徐志明，杨善让（东北电力学院动力系吉林１３２０１２）陶文铨（西安交通大学动力系西安７１００４９）关键词封闭腔；自然对流；数值模拟１引言三维封闭腔内孤立物体的自然对流换热在工程上具有广泛应用背景。例如；...     

倾斜封闭立方腔内多块孤立平板自然对流换热的实验研究

本文研究了平板的相互位置、腔体倾角及Rα数对倾斜封闭立方腔内多块孤立平板自然对流换热的影响．实验查明：在一定Rα数下存在着一个使每块板换热强度相当的间距；平均换热特性对位置很不敏感；平板的换热强度比相同Rα数下的大空间自然对流更弱，且随倾角增加，Nu数随Rα数增加的增长幅度逐渐减慢．     

几何位置对封闭方腔内水平孤立平板自然对流换热的影响

本文以冰箱蒸发器等实际问题为背景，以空气为介质用数值模拟方法研究了封闭方腔内水平孤立平板层流自然对流时Ｒａ数及平板位置变化对Ｎｕｍ的影响，得出了可供工程应用参考的结论。同时用烟可视化的实验验证了计算结果的正确性。     

充满金属泡沫的方腔自然对流数值模拟

本文使用格子-Boltzmann方法模拟了充满金属泡沫的封闭方腔自然对流.通过构建适用于上下边界绝热,左右边界恒定温度的封闭方腔自然对流的双分布函数热模型,模拟了充满不同孔隙率和不同孔密度的金属泡沫,在不同Ra数下的封闭方腔自然对流,结果表明金属泡沫具有强化传热的作用,并随着Ra数和孔密度的增大,换热增强.     

自然对流的LB方法模拟及其非线性特性

采用LB方法和QUICK差分方法模拟了方腔内竖直平板自然对流和底部加热方腔内自然对流换热问题.实现了LB方法对具有孤立体的封闭空间内耦合的流动和传热问题的数值模拟,所得数值结果与QUICK差分方法的数值结果及已有的烟可视化实验结果一致;两种方法对底部加热方腔内自然对流问题的预测结果均出现了静态分岔和动态分岔;根据底部加热方腔内自然对流换热数值结果给出的极限环型速度相图和功率谱表明,两种方法得到的数值结果的非线性特性一致.     

封闭方腔热磁对流强化换热的实验研究

自然对流是由非保守体积力驱动的流体流动,这种体积力可以是重力、离心力以及电磁场力等;通过外加梯度磁场来影响流体的流动,产生强化换热效果的方法是一种新型的强化技术．本文利用永磁体产生0．3 g左右的加速度,采用马赫干涉仪对封闭方腔内空气的热磁对流的强化换热现象进行了实验研究,得到了温度场的干涉图像,并与理论模拟结果进行了比较,发现磁场对空气的自然对流换热有一定的强化作用．     

周期性边界条件下多孔介质方腔内自然对流换热数值研究

对充满多孔介质的倾斜方腔,在其上下壁面绝热、右侧壁面维持恒温To及左侧壁面温度基于To随时间按正弦规律变化的情况下,采用Brinkman扩展达西模型和SIMPLER算法对方腔内的自然对流与换热进行了非稳态数值模拟,方腔倾角α的范围为0°～90°,Pr数为1,Ra数为106.计算研究不同的振荡频率和方腔倾角对方腔内对流换热的影响.数值模拟结果表明:振荡频率f为60π,方腔倾角为43°时,方腔内的换热最强.     

方腔散射辐射与自然对流耦合传热的数值模拟

本文在自主软件GTEA基础上,针对封闭方腔内参与性介质的辐射与自然对流耦合换热进行模拟,计算结果与文献中经典算例结果吻合较好。并分析了普朗克数和介质散射反照率的变化对方腔内温度和速度的影响,结果表明:随着普朗克数的减小,增强了辐射换热的影响;考虑各向同性散射介质时,随着散射反照率的减小,中心处的温度有所增大,靠近绝热壁面处,速度U的变化更剧烈。     

局部热壁面多孔介质方腔内自然对流的数值研究

本文对上下壁面绝热、左侧壁面长度为b的嵌装加热器部分维持恒定温度T_h而剩余部分绝热,且右侧壁面维持恒定温度T_c的多孔介质方腔内的自然对流换热进行了数值研究.在热壁面无量纲长度B=0.5(B=b/L)的条件下,综合研究了左侧壁面受热部分中心距上壁面的无量纲长度D(D=d/L)、Da数、Ra数和孔隙率对腔体内自然对流换热的影响.数值计算结果表明,左侧壁面受热部分位置的不同对腔体内自然对流换热有很大的影响,D在0.6附近取值时,Na数最大.Da数、Ra数对腔体的自然对流换热影响较大,而孔隙率对换热的影响较小.     

非保守力与自然对流

为研究自然对流流动和换热的实质,对重力场中自然对流的驱动力进行了进一步分析,提出了自然对流产生的充分必要条件。以二维封闭方腔内空气自然对流为对象,分别对三种不同的自然对流形式进行了数值模拟。结果表明:只要流体所受体积力为非保守力,即可产生自然对流。这为自然对流的控制提供了一定的理论依据。     

NATURAL CONVECTION FROM DISCRETE HEAT SOURCES IN A VERTICALLY VENTED RECTANGULAR ENCLOSURE

Natural convection in a discretely heated, vertically vented enclosure has been investigated experimentally. A vertically vented enclosure is one whose top and bottom boundaries are partially open, allowing ambient fluid to be drawn in by buoyancy. Mach-Zehnder interferometry was used to visualize the temperature field within the enclosure and to determine the local and average heat transfer characteristics of the discrete heat sources. A smoke-generation technique was used to visualize the flow structure. The experimental parameters investigated include Grashof number, vent gap width, and heater locations for a dual heater configuration and a single-enclosure aspect ratio. The vent gap width was varied between a fully open condition (discretely healed parallel plates) and nearly closed. The Mach-Zehnder interferograms and local heat transfer results suggest a complex buoyancy-driven flow field that depends intimately on the rent gap width. For some geometric configurations, the heat transfer coefficient was found to be nearly uniform across the length of the heater, while for others traditional boundary-layer-type heat transfer characteristics were found. In general, the heat transfer coefficient increased with increased vent gap spacing. However, a maximum in average Nusselt number was observed for a dimensionless gap spacing of G/W – 0.67 for one of the two dual-heater placement configurations studied. The heat transfer characteristics were compared to those of a single isolated plate and unobstructed parallel-plate channels.     

Laminar free convection heat transfer between vertical isothermal plates

The paper represents results on numerical investigation of flow and heat transfer between two isothermal vertical plates under laminar natural convection. A system of complete Navier–Stokes equations is solved for a two-dimensional gas flow between the plates along with additional rectangular regions (connected to inlet and outlet sections), whose characteristic sizes are much greater than the spacing between the plates. The calculations were performed over very wide ranges of Rayleigh number Ra = 10 ÷ 10⁵ and a relative channel length AR = L/w = 1 ÷ 500. The influence of the input parameters on the gas-dynamic and thermal structure of thermogravitational convection, the local and mean heat transfer, and also the gas flow rate between the plates (convective draft. We determined sizes of the regions and regime parameters when the local heat flux on the walls tends to zero due to the gas temperature approach to the surface temperature. It is shown that the mean heat transfer decreases as the relative channel length AR grows, whereas the integral gas flow rate (convective draft) and Reynolds number in the channel Re = 2wUm/ν increase. The use of a modified Rayleigh number Ra* = Ra · (w/L) (Elenbaas number) leads to generalization of calculation data on mean heat transfer. These data are in good agreement with the correlations for heat transfer [1, 2] and gas flow rate [3]. The reasons of variation of the data in the range of low Rayleigh numbers are discussed in detail.     

水平板自然对流换热的非线性特性

采用SIMPLE算法,QUICK差分方案,对封闭方腔内水平板自然对流换热进行了数值模拟.数值结果显示,低Ra数时流动和换热处于稳态,当Rayleigh数超过某一临界值时,流动和换热就会发生非稳态振荡,此时流动和换热表现出非对称性.对不同Rayleigh数,流动和换热通过单周期分岔从稳态过渡到非稳态,并通过倍周期分岔过渡到混沌.在混沌区,仍然会出现周期性窗口,并且数值结果与初始条件有关.     

Heat transfer at a laminar free convection and separated flow past a rib in a vertical channel with isothermal walls

The results of numerical computations of a free laminar convection and heat transfer between two parallel isothermal plates in the presence of a single rib on the channel surface are presented. The investigations have been conducted for a channel with the aspect ratio AR = L/w = 10, where L is the channel height, and w is the distance between the plates. An infinitely thin adiabatic rib was located on one of the channel walls in the middle of its height. The relative rib height l/w was varied in the range 0÷0.8. The wall temperature was higher than the ambient temperature, and the Rayleigh number was varied in the range Ra = 10²÷10⁵. The main attention has been paid to the study of the influence of the rib height and the Rayleigh number on local and integral heat transfer and the Reynolds number in the channel (the convective thrust). A fundamental difference in the heat transfer over the channel height has been shown on the ribbed wall and on a smooth surface. The computational results have been compared with the case of a symmetric distribution of the ribs on the both walls with the integral height equal to a single rib.     

TURBULENT NATURAL CONVECTION IN AN ENCLOSURE WITH LOCALIZED HEATING FROM BELOW

Natural-convection heat transfer in a cylindrical enclosure, heated partially from below by a disk-shaped heating surface and cooled from the top and the side, was investigated experimentally and numerically. Heat transfer measurements are presented for the range of Rayleigh numbers from 10⁸ to 2 × 10¹⁰. The total acceleration normal to the heating surface was varied from I to 100 times the standard gravitational acceleration. The heat flux was varied from 0.19 MW / m1 to 1.5 MW / m². The test results were correlated by an equation of the form Nu_D = 0.2Ra^0.325 _D for Prandtl number 2 and aspect ratio I. The flow field was studied numerically using FLUENT code.     

Turbulent periodic flow and heat transfer in a rectangular channel with detached V-baffles

This paper presents a numerical analysis of turbulent periodic flow and heat transfer in a rectangular channel with detached V-baffles. The computations are based on the finite volume method with the SIMPLE algorithm for handling the pressure–velocity coupling and using the QUICK scheme for the convection terms. Air is used as the test fluid with the air flow rate in terms of Reynolds numbers ranging from 3000 to 20,000. The effects of different detached-clearance ratios (c/H, CR) of 0.0, 0.05, 0.1, 0.15, and 0.2, baffles-pitch to square channel-diameter ratio (pitch ratio (p/H), PR) is 1.0, baffles-height to square channel-diameter ratio (blockage ratio (b/H), BR) is 0.10, and attack angle (α) is 45? on heat transfer, friction factor and thermal enhancement factor are investigated numerically. It is found that a pair of counter-rotating vortices (P-vortex) caused by the baffles can induce impingement/attachment flows repeatedly on the rectangular channel walls leading to a greater increase in the heat transfer over the test channel. The maximum thermal performance and heat transfer are found to be about 1.5 and 3.3, respectively for CR = 0.05 and Re = 3000, while the highest pressure loss is about 21.5 in the case of CR = 0.2 and Re = 20,000.     

多孔介质传热传质中耦合扩散效应的研究

本文采用加罚函数有限元方法数值模拟在温度梯度和浓度梯度同时存在时双浮升力自然对流及其伴随的传热传质,着重探讨了多孔介质中传热传质的交叉耦合扩散Soret效应和Dufour效应的基本影响规律,展示了流场、温度场和浓度场随热附加扩散准则数和扩散附加热准则数的变化状况。     

超临界CO2-丙烷混合物管内传热特性数值研究

采用SST k-w湍流模型对超临界CO₂/丙烷混合工质水平管内的传热特性进行数值模拟研究。管径d=4 mm,加热段L₂=800 mm;混合工质浓度配比为100/0、95/5、90/10、85/15、80/20、75/25;质量流速为150～250 kg·m^?2·s^?1;热流密度为30～40 kW·m^?2,入口温度293 K,入口压力7.5～30 MPa。随着丙烷浓度的增加,CO₂/丙烷二元混合工质的临界压力降低,临界温度升高,丙烷浓度从5%增加到25%,换热系数峰值降低6.19%～31.45%,但增加丙烷浓度可提高拟临界温度后的换热效果。P=7.5～8.5 MPa,换热系数有明显峰值;P=20～30 MPa,换热系数变化规律无明显峰值,并随压力的升高而减小。混合工质的换热系数随质量流速的增大而增大。同一流体温度所对应的换热系数,随着热流密度的增加而减小。     

超临界CO2-丙烷混合物管内传热特性数值研究

采用SST k-w湍流模型对超临界CO₂/丙烷混合工质水平管内的传热特性进行数值模拟研究。管径d=4 mm,加热段L₂=800 mm;混合工质浓度配比为100/0、95/5、90/10、85/15、80/20、75/25;质量流速为150～250 kg·m^?2·s^?1;热流密度为30～40 kW·m^?2,入口温度293 K,入口压力7.5～30 MPa。随着丙烷浓度的增加,CO₂/丙烷二元混合工质的临界压力降低,临界温度升高,丙烷浓度从5%增加到25%,换热系数峰值降低6.19%～31.45%,但增加丙烷浓度可提高拟临界温度后的换热效果。P=7.5～8.5 MPa,换热系数有明显峰值;P=20～30 MPa,换热系数变化规律无明显峰值,并随压力的升高而减小。混合工质的换热系数随质量流速的增大而增大。同一流体温度所对应的换热系数,随着热流密度的增加而减小。