侧加热腔内的自然对流

开展侧加热腔内自然对流的研究具有重大的环境及工业应用背景. 总结侧加热腔内水平温差驱动的自然对流的最新研究进展, 并概述相应的流动性质、动力机制和传热特性以及对不同无量纲控制参数的依赖也有重要的科学价值. 已取得的研究结果显示突然侧加热的腔内自然对流的发展可包括初始阶段、过渡阶段和定常或准定常阶段. 不同发展阶段的流动依赖于瑞利数、普朗特数及腔体的高宽比, 且定常或准定常阶段的流态可以是定常层流流动、非定常周期性流动或者湍流流动. 此外, 回顾了对流流动失稳机制的研究成果以及湍流自然对流方面的新进展. 最后, 展望了侧加热腔内的自然对流研究的前景.      

Quasi-periodicity and chaos in a differentially heated cavity

Convective flows of a small Prandtl number fluid contained in a two-dimensional vertical cavity subject to a lateral thermal gradient are studied numerically. The chosen geometry and the values of the material parameters are relevant to semiconductor crystal growth experiments in the horizontal configuration of the Bridgman method. For increasing Rayleigh numbers we find a transition from a steady flow to periodic solutions through a supercritical Hopf bifurcation that maintains the centro-symmetry of the basic circulation. For a Rayleigh number of about ten times that of the Hopf bifurcation, the periodic solution loses stability in a subcritical Neimark–Sacker bifurcation, which gives rise to a branch of quasiperiodic states. In this branch, several intervals of frequency locking have been identified. Inside the resonance horns the stable limit cycles lose and gain stability via some typical scenarios in the bifurcation of periodic solutions. After a complicated bifurcation diagram of the stable limit cycle of the 1:10 resonance horn, a soft transition to chaos is obtained. PACS 44.25.+f, 47.20.Ky, 47.52.+j     

Experimental observations of the thermal flow around a square obstruction on a vertical wall in a differentially heated cavity

The transient thermal boundary layer flow around a square obstruction placed at the middle of the hot wall in a differentially heated cavity is visualized using a shadowgraph technique. The results show that the thermal boundary layer flow, which is blocked by the obstruction, firstly forms an intrusion head under the obstruction (the lower intrusion head). Subsequently, the lower intrusion head bypasses the obstruction and reattaches to the down-stream boundary. During the reattachment process, a more complicated flow is induced, and eventually both the lower intrusion head and the thermal boundary layer destabilize. After the lower intrusion head is convected away, the thermal boundary layer flow re-stabilizes. At the quasi-steady state, the thermal boundary layer forms a double-layer structure, which is split into two sections by the obstruction. It is demonstrated that both the transient processes and the quasi-steady state flow structures of the thermal boundary layer are significantly altered by the obstruction in comparison with the case without the obstruction.     

Transient natural convection of low-Prandtl-number fluids in a differentially heated cavity

The transient convective motion in a two-dimensional square cavity driven by a temperature gradient is analysed. The cavity is filled with a low-Prandtl-number fluid and the vertical walls are maintained at constant but different temperatures, while the horizontal boundaries are adiabatic. A control volume approach with a staggered grid is employed to formulate the finite difference equations. Numerically accurate solutions are obtained for Prandtl numbers of 0·001, 0·005 and 0·01 and for Grashof numbers up to 1 × 10⁷. It was found that the flow field exhibits periodic oscillation at the critical Grashof numbers, which are dependent on the Prandtl number. As the Prandtl number is decreased, the critical Grashof number and the frequency of oscillation decrease. Prior to the oscillatory flow, steady state solutions with an oscillatory transient period were predicted. In addition to the main circulation, four weak circulations were predicted at the corners of the cavity.     

Natural convection in partially cooled tilted cavities

Two-dimensional numerical simulations of laminar natural convection in a partially cooled, differentially heated inclined cavities are performed. One of the cavity walls is entirely heated to a uniformly high temperature (heat source) while the opposite wall is partially cooled to a lower temperature (heat sink). The remaining walls are adiabatic. The tilt angle of the cavity is varied from 0° (heated from left) to −90° (heated from top). The fast false implicit transient scheme (FITS) algorithm, developed earlier by the same authors, is modified to solve the derived variables vorticity-streamfunction formulation. The effects of aspect ratio (AR), sink–source ratio and tilt angle on the average Nusselt number are examined through a parametric study; solutions are obtained for two Grashof numbers, 10⁵ and 10⁷. Flow patterns and isotherms are used to investigate the heat transfer and fluid flow mechanisms inside the cavity. © 1998 John Wiley & Sons, Ltd.     

An efficient Legendre-Galerkin spectral method for the natural convection in two-dimensional cavities

An efficient numerical method is developed for solving the natural convection in two-dimensional cavities. The numerical scheme is proposed by using second-order projection scheme in time direction and Legendre-spectral in spatial variable of the incompressible flow. Finally, a series of numerical examples are presented to demonstrate the efficiency of our algorithm. The numerical strategies developed in this article could be readily applied to study other incompressible fluid problems.     

Experiments on the cooling by natural convection of an array of vertical heated plates with constant heat flux

This paper presents experimental results that aim to document the phenomenon of cooling by natural convection of an array of vertical plates with uniform and equal heat fluxes. The working fluid is air. The effect of several factors on the plate temperature distribution was determined. These factors are the spacing of the plates, the existence of a floorlike flow obstruction near the entrance of the vertical channel, and the existence of a ceilinglike flow obstruction near the exit of the vertical channel. In several cases, these factors had a paramount effect on the plate temperature distribution. Examined also was the impact on the effectiveness of the natural convection cooling of a second row of plates positioned under the row of plates under investigation in an aligned or a staggered fashion.     

Unsteady flow and heat transfer adjacent to the sidewall wall of a differentially heated cavity with a conducting and an adiabatic fin

The unsteady natural convection flow adjacent to the finned sidewall of a differentially heated cavity is numerically investigated through comparisons between the cases with a conducting fin and an adiabatic fin. The results show that the flow and temperature structures in the transition to a periodic flow induced by a conducting fin are considerably different from those by an adiabatic fin. Based on the present numerical results, the temporal development and spatial structures of the flow adjacent to the finned sidewall are described, and instabilities are characterized. It is found that the conducting fin improves the transient convective flows in the cavity and enhances heat transfer across the cavity (by up to 52% in comparison with the case without a fin).     

Investigation of low-intensity convection regimes in a rectangular cavity with a heat flux on the boundary

The characteristics of heat transfer during natural thermogravitational fluid convection of low intensity in a rectangular cavity heated from below (cooled from above) are investigated. Local convection effects are studied. The dependence of local superheating (supercooling) on the Grashof number and the cavity side ratio is found for single-, two-and three-vortex steady motions. The limits of the convection regimes are estimated.     

On the reversed flow and oscillating wake in an asymmetrically heated channel

The detailed processes of flow reversal in a buoyancy-induced flow through a one-side-heated vertical channel of finite height were simulated numerically. It is of interest to note that the wake above the heated plate is oscillatory at high Rayleigh number and there exists a minimum in the transient variation of the average Nusselt number. Additionally, the predicted steady average Nusselt number and induced flow rate are correlated by empirical equations.     

Nonlinear dynamics between two differentially heated vertical plates in the presence of stratification

We consider the numerical simulation of the flow between infinite, differentially heated vertical plates with positive stratification. We use a two-dimensional Boussinesq approximation, with periodic boundary conditions in the vertical direction. The relative stratification parameter ${\gamma=(\frac{1}{4}Ra S)^{1/4}}$ , where Ra is the Rayleigh number and S the adimensional stratification, is kept constant and equal to 8. The Prandtl number is 0.71. We derive a complex Ginzburg-Landau equation from the equations of motion. Coefficients are computed analytically, but we find that the domain of validity of these coefficients is small and rely on the numerical simulation to adjust the coefficients over a wider range of Rayleigh numbers. We show that the Ginzburg-Landau equation is able to accurately predict the characteristics of the periodic solution at moderate Rayleigh numbers. Above the primary bifurcation at Ra = 1.63 × 10⁵, the Ginzburg-Landau model is found to be Benjamin-Feir unstable and to be characterized by modulated traveling waves and phase-defect chaos, which is supported by evidence from the DNS. As the Rayleigh number is increased beyond Ra = 2.7 × 10⁵, nonlinearities become strong and the flow is characterized by cnoidal waves.     

Effects of anisotropy on the free convection from a vertical plate embedded in a porous medium

The effects of anisotropy on the steady laminar boundary-layer free convection over a vertical impermeable surface are analysed by using the method of integral relations. If the permeability in the direction orthogonal to the plate is greater than the permeability along the plate, then there is an increase in the temperature field.     

Numerical study of flow and thermal behaviour of lid‐driven flows in cavities of small aspect ratios

Numerical study has been performed to investigate the effects of cavity shape on flow and heat transfer characteristics of the lid‐driven cavity flows. Dependence of flow and thermal behaviour on the aspect ratio of the cavities is also evaluated. Three types of the cross‐sectional shape, namely, circular, triangular, and rectangular, and four aspect ratios, 0.133, 0.207, 0.288, and 0.5, are taken into account to construct twelve possible combinations; however, attention is focused on the small‐aspect‐ratio situations. Value of the Reynolds number considered in this study is varied between 100 and 1800. For the cases considered in this study a major clockwise vortex driven by the moving lid prevailing in the cavity is always observed. When the Reynolds number is fixed, the rectangular cavity produces strongest lid‐driven flow, and the triangular cavity weakest. For the cases at small aspect ratio and low Reynolds number, the streamlines appear symmetric fore‐and‐aft with respect to the central line at x/L = 0.5. Data for the local and average Nusselt numbers are also provided. For rectangular cavities, it is observed that case 1/5R produces the highest average Nusselt number at any Reynolds number. Among the twelve possible geometric cases considered herein, the highest and lowest average Nusselt numbers are found with cases 1/6T and 1/2C, respectively. Copyright © 2006 John Wiley & Sons, Ltd.     

Natural convection in an inclined quadrantal cavity heated and cooled on adjacent walls

Natural convective flow and heat transfer in an inclined quadrantal cavity is studied experimentally and numerically. The particle tracing method is used to visualize the fluid motion in the enclosure. Numerical solutions are obtained via a commercial CFD package, Fluent. The working fluid is distilled water. The effects of the inclination angle, ? and the Rayleigh number, Ra on fluid flow and heat transfer are investigated for the range of angle of inclination between 0° ? ? ? 360°, and Ra from 10⁵ to 10⁷. It is disclosed that heat transfer changes dramatically according to the inclination angle which affects convection currents inside, i.e. flow physics inside. A fairly good agreement is observed between the experimental and numerical results.     

Transient natural convection in a rectangular enclosure with one heated side wall

This paper describes a numerical and theoretical study of the transient natural convection heating of a two-dimensional rectangular enclosure filled with fluid. The heating is applied suddenly along one of the side walls, while the remaining three walls are maintained insulated. It is shown that the process has two distinct phases, an early period dominated by conduction and a late period dominated by convection. The scaling laws for the heat transfer rate and the effectiveness (energy storage fraction) are determined based on scale analysis. These theoretical results are confirmed by numerical experiments conducted in the domain Ra = 10³−10⁶, Pr = 7, A = 1, where Ra is the Rayleigh number based on height and initial temperature difference, Pr is the Prandtl number, and A is the height/length ratio of the enclosure. Correlations for heat transfer rate and effectiveness are constructed by comparing the theoretical scaling laws with the numerical results.     

An experimental study of the effect of wall temperature nonuniformity on natural convection in an enclosure heated from the side

An experimental study of natural convection in a parallelepipedal enclosure induced by a single vertical wall is described. The upper half of this wall was warm and the lower half cold. The other enclosure walls were insulated. The temperature and flow measurements were performed in the high Rayleigh number regime (10¹⁰<Ra<5×10¹⁰) by using water as the working fluid. The Rayleigh number was based on the enclosure height and the temperature difference between the warm and the cold part of the driving wall. The flow field featured two flat cells, one filled with warm fluid along the top horizontal wall, and the other filled with cold fluid along the bottom horizontal wall. Each of these cells was surrounded by an additional cell as tall as half the enclosure height. The above flow structure prohibited extensive thermal contact between warm and cold fluid, thus limiting the role of convection on the heat transfer process in the cavity. The findings of this study differ significantly from the findings of previous studies based on the ‘classical’ enclosure model possessing two isothermal vertical walls, the one warm and the other cold, and support the view that the use of ‘more realistic‘ temperature boundary conditions in enclosure natural convection needs careful examination.     

Non-linear dynamics and pattern formation in a vertical fluid layer heated from the side

We study both experimentally and numerically the convective flow in a tall vertical slot with differently heated walls. The flow is investigated for the fluid with the Prandtl number Pr=26, which is large enough to ensure the traveling waves as primary instability and small enough to prevent boundary layer convection. The flow evolution is determined on the base of the visual observations, power spectra and amplitude analysis. In the numerical simulations of two- and three-dimensional flows, we accept an assumption of an infinite fluid layer. The satisfactory agreement with experiment is observed, and the sequence of convection states is discovered. It starts with a plane-parallel flow as primary solution, which becomes unstable to two counter-propagating waves. It is followed by a tertiary three-dimensional flow in the form of wavy traveling waves. As the Grashof number is increased even further, a chaotically oscillating cellular pattern consisting of the pieces of broken waves arises. The formation of a structure in the form of the vertical rolls chaotically modulated along axes concludes this complicated picture.     

Numerical study of natural convection flow in a vertical slot

The stability of convective motion, generated by a lateral temperature difference across a vertical slot, is studied numerically over a range ofGr=5000 to 1.5 × 10⁵,Pr=0.01 to 10, andA=8,16 and 20. Various cellular flow structures and temperature patterns are illustrated. Several branches of solutions characterized by different numbers of the cells in the flow patterns as well as by both steady and unsteady multicellular patterns are found for low-Prandtl-number fluid in the vertical slot. Meanwhile, the behaviors of the temperature variation and heat transfer are also discussed. The project supported by the National Natural Science Foundation of China (59776011) and by the Returnee from Abroad Funding of Academia Sinica.     

Natural convective heat transfer in water enclosed between pairs of differentially heated vertical plates

This paper presents the results of an experimental study of buoyancy-driven convective heat transfer between three parallel vertical plates, symmetrically spaced with water as the intervening medium. The centre plate was electrically heated, while the other side plates were water-cooled forming two successive parallel vertical channels of dimensions 20 cm × 3.5 cm × 35 cm (length W, gap L, height H) each. Top, bottom and sides of the channels were open to water in the chamber which is the novel aspect of this study. Plate surface temperature and bath temperature at different levels of height from the bottom of channel were measured by K-type thermocouples. Experimental data have been correlated as under: