Nonstationary laminar heat convection in a closed domain for a given heat flux

Characteristic modes of the time development of nonstationary heat convection in a closed planar domain upon a sudden supply of heat from the lateral surface are considered for Rayleigh numbers 10³–10⁷. Estimates of the boundaries of the beginning of the influence of convection on the temperature field and the buildup of a quasistationary convection mode in the range of Rayleigh and Fourier numbers are given. Characteristics of the circulation flow, the singularities of the temperature-field configuration and of the heat transfer from the wall to the fluid, are investigated. The mechanism for the origination and disappearance of vertical temperature differences, caused by convection, and the dependence of the vertical temperature differences on the Rayleigh and Fourier numbers, on the thermal mode of the boundary, and the domain geometry, are considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 109–117, July–August, 1970.The author is grateful to T. D. Pirumov and T. V. Volokitin for assistance in performing the computations.     

Transient natural convection between concentric and vertically eccentric spheres with mixed boundary conditions

Transient analysis has been investigated numerically to determine heat transfer by natural convection between concentric and vertically eccentric spheres with constant heat flux on the inner wall and a specified isothermal temperature on the outer wall. The governing equations, in terms of vorticity, stream function and temperature are expressed in a spherical polar coordinate system. The alternating direction implicit method and the successive over-relaxation techniques are applied to solve the finite difference form of governing equations. A physical model is introduced which accounts for the effects of fluid buoyancy as well as eccentricity of the outer sphere. Transient solutions of the entire flow field are obtained for a range of modified Rayleigh number (10³<Ra?<5×10⁵), for a Prandtl number of 0.7 and a radius ratio of 2.0, with the outer sphere near the top and bottom of the inner sphere (ε=±0.625). Results of the parametric study conducted further reveal that the heat and flow fields are primarily dependent on the modified Rayleigh number and the eccentricity of the spherical annulus. The results of average Nusselt numbers are also compared with the results obtained for flow between two isothermal spheres.     

Turbulent penetrative convection with an internal heat source

An infinite fluid with a vertical cubic temperature profile in the absence of fluid motion is considered as a model for penetrative convection in which a central unstably stratified fluid layer is bounded above and below by stably stratified layers. Turbulence statistics from direct and large eddy numerical simulations for the mean temperature gradient, the velocity and temperature variances and the heat flux are presented for Rayleigh numbers R up to four orders of magnitude above critical. By means of a simplified second-moment closure, analytical scaling laws for the statistics are determined. For high Rayleigh numbers, the mean temperature gradient approaches zero in a central well-mixed layer, a reduced positive (stable) value in upper and lower partially mixed layers, and an unmixed value far above and below. The temperature variance is a factor of R^1/3 larger in the partially mixed layers compared to the well-mixed layer; the velocity variance and heat flux scales the same in both layers. Approximation of the three layers by a two layer model yields an estimate for the height of the mixed layer: the height decreases slowly with increasing Rayleigh number and at the highest Rayleigh number simulated is approximately 30% longer than the unstable layer in the absence of fluid motion.     

Natural convection in an inclined enclosure containing internal energy sources and cooled from below

Natural convection in an inclined enclosure from below and containing internally heated fluid has been investigated using a finite difference calculation procedure. Results have been obtained for Rayleigh number values up to 10⁶ and for inclination angles of 30 and 60°. For internal Rayleigh numbers that are much larger than the external Rayleigh number, the flow rises in the interior and moves down both the hot and cold walls. On the other hand, if the external Rayleigh number has a larger magnitude, the flow moves upwards along the hot surface and downwards along the cold surface. For the latter situation, the inner core is multicellular in nature at large external Rayleigh numbers. The average heat flux ratio along the cold surface (convective heat flux/corresponding conduction heat flux) increases with increasing external Rayleigh number and decreasing internal ratio is non-monotonic in nature. The heat flux ratio along both surfaces is observed to be strongly dependent on the inclination angle at high external Rayleigh numbers. A maximum in the local heat flux along the cold surface is obtained in the vicinity of x/L = 1 where hot fluid, either from the interior or directly from the opposite hot wall, meets the surface. Along the hot wall, a maximum in the heat flux ra flo     

Two numerical modelings of free convection heat transfer using nanofluids inside a square enclosure

This work describes the numerical simulation of natural convection heat transfer of Cu–water nanofluids in a square enclosure for Rayleigh numbers varying from 10³ up to 10⁵. Two different numerical approaches were used: the finite volume method and the finite element method. The nanofluids were assumed to be single-phase fluids with modified thermal properties obtained from experimental results and theoretical models. The results showed that the Nusselt number for nanofluids was basically the same as that obtained for the base fluid. Therefore, the enhancement observed in the heat transfer coefficient was significant due to the augmentation in the thermal conductivity.     

Influence of boundary conditions on the convective stability of a rotating horizontal layer of liquid

A theoretical study is made of the critical curves for the onset of convection in a plane horizontal layer of liquid rotating with constant angular velocity for different conditions on the boundary of the layer. It is shown that, in contrast to Chandresekhar's curves [1] obtained under the condition of constancy of the temperature on the boundaries, the curves for a constant heat flux lie significantly lower, so that convection occurs earlier for all Taylor numbers. At large Taylor numbers all the stability curves, as in [1], tend to the asymptotes R_C Ta^2/3, where Ta is the Taylor number and R_C is the critical Rayleigh number. A similar investigation for a nonrotating liquid was made in [2].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 124–129, May–June, 1988.We thank G. S. Golitsyn for the proposed subject and constant interest in the work.     

Natural convection in a rectangular cavity with wall temperature decreasing at a uniform rate

The transient natural convection in a fluid contained in a rectangular enclosure, the wall of which is maintained at a uniform temperature which changes at a steady rate, is approached by a numerical method. Numerical solutions are obtained forPr=0.73, 7.3 and 73 and a range of Rayleigh numbersRa=10² ～ 10⁸. At relatively low Rayleigh numbers the flow is characterized by the development of double cells with flow up the center and down the sidewalk However it was found that an increase of the Rayleigh number leads to the development of strong secondary circulation on the axis of symmetry of the cavity near the top wall. Thus, as the Rayleigh number is increased the secondary cells grow in size. The effects of the secondary cells on the temperature field and heat transfer coefficients are discussed. Most results are obtained for the case of a square cavity (E=2) but the influence of the aspect ratio of the cavity is also studied forE=1 and 4.     

Experimental research on heat transfer of natural convection in vertical rectangular channels with large aspect ratio

This work presents the experimental research on the steady laminar natural convection heat transfer of air in three vertical thin rectangular channels with different gap clearance. The much higher ratio of width to gap clearance (60–24) and the ratio of length to gap clearance (800–320) make the rectangular channels similar with the coolant flow passage in plate type fuel reactors. The vertical rectangular channels were composed of two stainless steal plates and were heated by electrical heating rods. The wall temperatures were detected with the K-type thermocouples which were inserted into the blind holes drilled in the steal plates. Also the air temperatures at the inlet and outlet of the channel were detected. The wall heat fluxes added to the air flow were calculated by the Fourier heat conduction law. The heat transfer characteristics were analyzed, and the average Nusselt numbers in all the three channels could be well correlated with the Rayleigh number or the modified Rayleigh number in a uniform correlation. Furthermore, the maximum wall temperatures were investigated, which is a key parameter for the fuel’s integrity during some accidents. It was found that even the wall heat flux was up to 1500 W/m², the maximum wall temperature was lower than 350 °C. All this work is valuable for the plate type reactor’s design and safety analysis.     

Natural convection and conduction heat transfer in open shallow cavities with bounding walls

A numerical study has been carried out on inclined open shallow cavities, which are formed by a wall and horizontal fins. Constant heat flux is applied on the surface of the wall inside the cavity while its other surface was kept isothermal. The wall and the fins are conductive. Conjugate heat transfer by natural convection and conduction is studied by numerically solving equations of mass, momentum and energy. Streamlines and isotherms are produced, heat and mass transfer is calculated. A parametric study is carried out using following parameters: Rayleigh number from 10⁶ to 10¹², conductivity ratio from 1 to 60, open cavity aspect ratio from 1 to 0.125, dimensionless end wall thickness from 0.05 to 0.20, horizontal walls from 0.01 to 0.15 and inclination of the end wall from 90° to 45°. It is found that the volume flow rate and Nusselt number are a decreasing function of the cavity aspect ratio, horizontal fin thickness and conductivity ratio. They are an increasing function of end wall thickness and inclination angle, except in the latter case optima exist at high Rayleigh numbers.     

Natural convection heat transfer in a uniformly heated horizontal pipe

Natural convection heat transfers inside horizontal pipes were measured. The Rayleigh numbers were varied from 6.8 × 10⁸ to 1.5 × 10¹², while the Prandtl number was fixed at 2,094. Based on the analogy concept, a copper sulfate electroplating system was adopted to measure mass transfer rates in place of heat transfer rates. Test results using single-piece electrodes were in good agreement with the work of Sarac and Korkut. The angle-dependent mass transfer rates, measured using piecewise electrodes, were compared with the results of studies on natural convection in concentric annuli, and showed similar trends. The experiments were expanded to the turbulent region, and a transition criterion was proposed. Angle-dependent natural convection heat transfer correlations for the laminar and turbulent regions were derived.     

Convection and heat transfer in a spherical vessel,partly filled with liquid,under low-gravity conditions

In [1] the problem of natural and thermocapillary convection in a spherical vessel containing a bubble under low-gravity conditions, i.e., at low Bond numbers (Bo 1), was examined in one of the limiting cases — where the bubble is located in the center of the vessel. The results of [1] and experimental data, however, indicate that when heat is supplied from outside over a long period, the most probable location of the bubble under low-gravity conditions is at the vessel wall. In this paper, which is a continuation of [1], convection and heat transfer in the latter case are investigated. Possible locations of the bubble at the top and bottom of the vessel relative to the resultant of the weak mass forces are discussed. It is shown that natural and thermocapillary convection contribute to the increase in the mean free-surface temperature, which determines the increase in pressure in the closed vessel for a prescribed heat flux. The rates of increase of this temperature are compared in the cases considered here and in [1–4], where there is a fuller bibliography relating to convective heat and mass transfer under low-gravity conditions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 154–159, March–April, 1976.     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

Thermal effects of the micropolar fluid parameters on the laminar free convection in spherical annuli with mixed boundary conditions

The effects of micro-rotation and vortex viscosity in micropolar fluids have been investigated numerically to determine heat transfer by natural convection between concentric and vertically eccentric spheres with specified mixed boundary conditions. Calculations were carried out systematically for several different eccentricities and a range of modified Rayleigh numbers to determine the average Nusslet numbers which are affected by the micropolar parameters (F) of the flow and temperature fields. The skin friction stress on the walls has also been studied and discussed. The governing equations, in terms of vorticity, stream function, temperature and angular momentum are expressed in a spherical polar coordinate system. Results were obtained for steady heat-transfer in spherical annuli at a Prandtl number of 0.7, with the modified Rayleigh number ranging from 10³ to 5 × 10⁵, for a radius ratio of 2.0 and eccentricities varying from −0.625 to +0.625. Comparisons are attempted between the Newtonian fluid and micropolar fluid.     

Mean square velocities of bénard convection and heat transfer

The approximate formula K a^–2R(N–1), where a is a constant near 9 and R and N are the Rayleigh and Nusselt numbers, was proposed in [1] for the dimensionless kinetic energy K of convection in a horizontal layer of liquid. It is shown in the present paper that this expression is exact in linear and weakly nonlinear convection theory when the velocity and temperature fields are represented analytically [2–4]. The valuea is found to be 8.76 when the upper and lower boundaries of the layer are solid walls. The results are given of numerical calculations of the kinetic energy of the convection and the heat transfer in a wide range of Rayleigh numbers (up to 44 000) and Prandtl numbers (0.025 P 15). Analysis of the results shows that a is in fact a weak function of both R and P. If this is also the case at large R, it indicates a certain breaking of scaling of the mean convection characteristics at sufficiently large values of the Rayleigh number. It also indicates why laboratory experiments give values of n in the dependence N Rⁿ which are generally slightly less than the theoretical value n = 1/3.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 11–18, September–October, 1981.We should like to thank N. F. Vel'tishchev for providing first data of the numerical experiments of [13–15].     

Effects of surface radiation on natural convection in a rayleigh-benard square enclosure: steady and unsteady conditions

Coupled laminar natural convection with radiation in air-filled square enclosure heated from below and cooled from above is studied numerically for a wide variety of radiative boundary conditions at the sidewalls. A numerical model based on the finite difference method was used for the solution of mass, momentum and energy equations. The surface-to-surface method was used to calculate the radiative heat transfer. Simulations were performed for two values of the emissivities of the active and insulated walls (ɛ₁=0.05 or 0.85, ɛ₂=0.05 or 0.85) and Rayleigh numbers ranging from 10³ to 2.3×10⁶ . The influence of those parameters on the flow and temperature patterns and heat transfer rates are analyzed and discussed for different steady-state solutions. The existing ranges of these solutions are reported for the four different cases considered. It is founded that, for a fixed Ra, the global heat transfer across the enclosure depends only on the magnitude of the emissivity of the active walls. The oscillatory behavior, characterizing the unsteady-state solutions during the transitions from bicellular flows to the unicellular flow are observed and discussed.     

Transient behavior of heat removal from a cylindrical heat storage vessel packed with spherical porous particles

The transient heat transfer behavior in the case of heat removal from a cylindrical heat storage vessel packed with spherical particles was investigated experimentally for various factors (flow rate, diameter of spherical particles packed, temperature difference between flowing cold air and spherical particles accumulating heat, and physical properties of spherical particles). The experiments were covered in ranges of Reynolds number based on the mean diameter of spherical particles packed Re_d = 10.3–2200, porosity?=0.310 to 0.475, ratio of spherical particle diameter to cylinder diameterd/D = 0.0075–0.177 and ratio of length of the cylinder to cylinder diameterL/D=2.5–10. It was found that especially the flow rate and the dimension of spherical particles played an important role in estimating the transient local heat transfer characteristics near the wall of the cylindrical vessel in the present heat storage system. As flow rate and diameter of spherical particles were increased under a given diameter of the cylinder heat storage vessel, the mean heat transfer coefficient between the flow cold air and the hot spherical particles increased and the time period to finish removing heat from the vessel reduced. In addition, the useful experimental correlation equations of mean heat transfer coefficient between both phases and the time period to finish removing heat from the vessel were derived with the functional relationship of Nusselt numberNu _d=f [modified Prandtl numberPr ^* (d/D), Re_d) and Fourier numberFo = f(d/D, L/D, Pr^*, Re_d).     

Effect of the Coriolis force on the onset of convection in a layer with a fixed heat flux on the boundaries

The effect of the Coriolis force on the onset of convection in a plane horizontal layer of viscous fluid with a fixed heat flux on the rigid lower and free upper boundaries is investigated. Expressions for the critical Rayleigh numbers and wave number are obtained analytically in the rapid rotation limit.Perm'. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–46, May–June, 1994.     

R410A condensation inside a commercial brazed plate heat exchanger

This paper presents the heat transfer coefficients and the pressure drop measured during HFC-410A condensation inside a commercial brazed plate heat exchanger: the effects of saturation temperature, refrigerant mass flux and vapour super-heating are investigated. The heat transfer coefficients show weak sensitivity to saturation temperature and great sensitivity to refrigerant mass flux and vapour super-heating. At low refrigerant mass flux (<20 kg/m² s) the saturated vapour condensation heat transfer coefficients are not dependent on mass flux and are well predicted by Nusselt [W. Nusselt, Die oberflachenkondensation des wasserdampfes, Energy 60 (1916) 541–546, 569–575] analysis for vertical surface: the condensation process is gravity controlled. For higher refrigerant mass flux (>20 kg/m²s) the saturated vapour condensation heat transfer coefficients depend on mass flux and are well predicted by Akers et al. [W.W. Akers, H.A. Deans, O.K. Crosser, Condensing heat transfer within horizontal tubes, Chem. Eng. Prog. Symp. Series 55 (1959) 171–176] equation: forced convection condensation occurs. In the forced convection condensation region the heat transfer coefficients show a 30% increase for a doubling of the refrigerant mass flux. The condensation heat transfer coefficients of super-heated vapour are 8–10% higher than those of saturated vapour and are well predicted by Webb [R.L. Webb, Convective condensation of superheated vapor, ASME J. Heat Transfer 120 (1998) 418–421] model. A simple linear equation based on the kinetic energy per unit volume of the refrigerant flow is proposed for the computation of the frictional pressure drop.