Numerical simulation of natural convective heat transfer and fluid flow around a heated cylinder inside an enclosure

 A collocated, non-orthogonal grid based finite volume technique has been applied for investigating the two dimensional natural convective flow and heat transfer around a heated cylinder kept in a square enclosure. The effects of different enclosure wall thermal boundary conditions, fluid Prandtl number and the ratio between enclosure and cylinder dimensions (aspect ratio) upon the flow and thermal features, have been systematically studied. It is observed that the patterns of recirculatory flow and thermal stratification in the fluid are significantly modified, if any of these parameters is varied. The overall heat transfer rates are also affected due to the changes in the flow and temperature patterns. The study presents useful observations regarding the variation of local Nusselt number along each wall, for the different cases considered. Received on 2 August 2000 / Published online: 29 November 2001     

Transport through baffles in bottom heated top cooled enclosures: parametric studies

This paper presents parametric studies on the heat transfer and fluid exchange through single-hole baffles located at the median height in bottom heated top cooled enclosures. Results indicate that when the baffle area-opening ratio is smaller than 2%, the heat transfer in the enclosure is dominated by the transport through the baffle opening. Even with such small baffle openings, increasing the enclosure aspect ratio still enhances the transport across the baffle. The characteristic length scale of flow in the enclosure is a combination of baffle opening diameter and the chamber height. The Nusselt number that characterize the heat transfer through the baffle-hole is linearly correlated with the Rayleigh number based on baffle opening diameter and the temperature difference between the bulk temperatures in the two chambers, while no effects of Prandtl numbers are observed. The mechanism of transport across the baffle opening varies from conduction dominated, combined conduction and convection, and convection dominated regimes as Rayleigh number increases.     

Numerical study of steady magneto-convection around an adiabatic body inside a square enclosure in low Prandtl numbers

In this paper, the effects of Prandtl number on the steady magneto-convection around a centrally located adiabatic body inside a square enclosure are numerically investigated. Two-dimensional nonlinear governing equations are discretized using the control volume method and hybrid scheme. The equations are solved using SIMPLER algorithm. The results are displayed in the form of streamlines and isotherms when the Rayleigh number varies between 10³ and 10⁶, the Hartmann number changes between 0 and 100 and the Prandtl number ranges between 0.005 and 0.1. The ratio of the buoyancy force to the Lorentz force (Ra/Ha ²) is introduced as an index to compare the contribution of natural convection and magnetic field strength on heat transfer. The results obtained from numerical modeling show that the Prandtl number has not considerable effect on heat transfer at low Rayleigh numbers. The effect of magnetic field strength on convection is increased by increasing Prandtl number. The effect of Prandtl number on the average Nusselt number in the presence of a magnetic field is less than the case without a magnetic field.     

MHD natural convection flow and heat transfer in a laterally heated partitioned enclosure

This study looks at MHD natural convection flow and heat transfer in a laterally heated enclosure with an off-centred partition. Governing equations in the form of vorticity–stream function formulation are solved using the polynomial differential quadrature (PDQ) method. Numerical results are obtained for various values of the partition location, Rayleigh, Prandtl and Hartmann numbers. The results indicate that magnetic field significantly suppresses flow, and thus heat transfer, especially for high Rayleigh number values. The results also show that the x-directional magnetic field is more effective in damping convection than the y-directional magnetic field, and the average heat transfer rate decreases with an increase in the distance of the partition from the hot wall. The average heat transfer rate decreases up to 80% if the partition is placed at the midpoint and an x-directional magnetic field is applied. The results also show that flow and heat transfer have little dependence on the Prandtl number.     

Natural convection in cavities containing internal sources

 Numerical predictions are reported for buoyancy-induced circulations in a 2D closed cavity with internal heat sources. Two cases are considered: (A) two vertical plates with uniform heat generation, forming a short vertical channel within the enclosure; and, (B) a rectangular heating block with uniform wall temperature, placed central in the enclosure. Air, with a Prandtl number 0.71, is considered as the working fluid. The vertical enclosure walls are isothermal, while the horizontal enclosure walls are adiabatic. Results are presented for two values of the Grashof number, one below the stability limit for laminar flow, and one well above it. In those latter cases, the long-term behaviour of the numerical solution is time-dependent, i.e. no steady-state can be reached. Heat transfer results are compared with predictions from standard correlations for isolated surfaces. Received on 17 January 2000     

Heat transfer over an unsteady stretching surface

Similarity solution of the laminar boundary layer equations corresponding to an unsteady stretching surface have been studied. The governing time-dependent boundary layer are transformed to ordinary differential equations containg Prandtl number and unsteadiness parameter. The effect of various govern-ing parameters such as Prandtl number and unsteadiness param-eter which determine the velocity and temperature profiles and heat transfer coefficient are studied.     

Simultaneous heat and mass transfer under unsteady mixed convection along a vertical slender cylinder embedded in a porous medium

Unsteady laminar mixed convection flow (combined free and forced convection flow) along a vertical slender cylinder embedded in a porous medium under the combined buoyancy effect of thermal and species diffusion has been studied. The effect of the permeability of the medium as well as the magnetic field has been included in the analysis. The partial differential equations with three independent variables governing the flow have been solved numerically using a implicit finite difference scheme in combination with the quasilinearization technique. Computations have been carried out for accelerating, decelerating and oscillatory free stream velocity distributions. The effects of the permeability of the medium, buoyancy forces, transverse curvature and magnetic field on skin friction, heat transfer and mass transfer have been studied. It is found that the effect of free stream velocity distribution is more pronounced on the skin friction than on the heat and mass transfer. The permeability and magnetic parameters increase the skin friction, but reduce the heat and mass transfer. The skin friction, heat transfer and mass transfer are enhanced due to the buoyancy forces and curvature parameter. The heat transfer is strongly dependent on the viscous dissipation parameter and the Prandtl number, and the mass transfer on the Schmidt number.     

Natural convection in an inclined rectangular porous cavity with uniform heat flux from the side

The problem of natural convection in an inclined rectangular porous layer enclosure is studied numerically. The enclosure is heated from one side and cooled from the other by a constant heat flux while the two other walls are insulated. The effect of aspect ratio, inclination angle and Rayleigh number on heat transfer is studied. It is found that the enclosure orientation has a considerable effect on the heat transfer. The negative orientation sharply inhibits the convection and consequently the heat transfer and a positive orientation maximizes the energy transfer. The maximum temperature within the porous medium can be considerably higher than that induced by pure conduction when the cavity is negatively oriented. The peak of the average Nusselt number depends on the Rayleigh number and the aspect ratio. The heat transfer between the two thermally active boundaries is sensitive to the effect of aspect ratio. For an enclosure at high or low aspect ratio, the convection is considerably decreased and the heat transfer depends mainly on conduction.     

Heat transfer characteristics of internally heated liquid pools at high Rayleigh numbers

Detailed numerical analysis is presented for buoyancy driven flow of a Newtonian fluid contained in a square enclosure for high Rayleigh (Ra) numbers. Natural convection is due to internal heating sources, which are assumed to be uniformly distributed within the enclosure. All walls of the cavity are maintained at constant temperature. Flow and heat transfer characteristics are investigated for a Ra number range of 10⁷ to 10¹² while Prandtl (Pr) number is taken to be 7.0. Governing equations (in primitive variables) are discretised using control volume technique based on staggered grid formulation. These equations are solved using SIMPLER algorithm of Patankar. Flow and heat transfer characteristics, streamlines, isotherms and average wall Nusselt (Nu) number, are presented for whole range of Ra number considered. Finally, present results for average wall Nu numbers are compared with experimental observations obtained from open literature. It is concluded that both results are in very good agreement, which confirmed the accuracy of the scaling used for present investigation. Received on 15 November 1999     

Thin film flow and heat transfer on an unsteady stretching sheet with internal heating

This research studied the influence of internal heat generation on flow and heat transfer in a thin liquid film on an unsteady stretching sheet. The velocity and temperature fields were solved using the Homotopy Analysis Method (HAM), taking a general surface temperature into consideration. The analytical series solution are presented and the numerical results obtained are tabulated. The effects of unsteadiness parameter, Prandtl number and temperature-dependent parameter in this study are discussed and presented graphically via the velocity and temperature profiles.     

Unsteady flow with heat and mass transfer due to impingement of slot jet on an inclined plate

The unsteady laminar incompressible boundary layer flow due to a two-dimensional slot jet on a flat plate at an angle of attack has been studied. The unsteadiness in the flow field is due to the free stream velocity distribution or wall temperature (concentration) which varies with time. The governing partial differential equations in primitive variables have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The effect of the variation of the free stream velocity distribution with time is found to be more pronounced on the skin friction than on the heat or mass transfer. The Prandtl number and the variation of the wall temperature with time strongly affect the heat transfer. Similarly, the Schmidt number and the variation of the concentration at the wall with time strongly affect the mass transfer. Beyond a certain critical value of the viscous dissipation parameter, the plate gets heated instead of being cooled.     

Laminar natural convection in a fully partitioned enclosure containing fluid with nonlinear thermophysical properties

The effects of a heat conducting partition on the laminar natural convection heat transfer and fluid flow were obtained by comparing the numerical and experimental results for a cubic enclosure without and with a partition. The two opposite vertical walls of the enclosure were isothermal at different temperatures. The working fluid was glycerol. The complete vertical partition, made of Plexiglass, was positioned in the middle of the enclosure. The visualizations of the velocity and temperature fields were obtained by using respectively, Plexiglass and liquid crystal particles as tracers. A middle plane perpendicular to the partition was numerically modeled. The steady two-dimensional model accounted for the variable thermophysical properties of the fluid. The finite volume method based on the finite difference approach was applied. The convective terms were approximated using a deferred correction central difference scheme. The velocity and temperature fields and the distribution of the local and average Nusselt numbers were found as a function of the Rayleigh (38 000 <Ra <369 000) and Prandtl (2700 < Pr < 7000) numbers.     

Numerical analysis on flow dynamics and heat transfer of falling liquid films with interfacial waves

Flow dynamics and heat transfer of falling liquid films with interfacial waves flowing on a vertical plate have been studied with originally proposed numerical simulation method. To discretize basic equations a staggered grid fixed on a physical space is employed. A small amplitude disturbance generated at inflow boundary develops to a solitary wave which consists of a large amplitude roll wave and small amplitude capillary waves. Instantaneous streamwise velocity profiles at the wave crest and trough are very different from a laminar flow. A circulation flow occurs in the roll wave and it affects temperature distributions, especially the strong effect is observed for high Prandtl number liquids. The interfacial wave enhances the heat transfer by two kinds of effects which are a film thinning effect and a convection effect. The dominating effect depends on the Prandtl number. Received on 23 December 1998     

Steady mixed convection flow of Maxwell fluid over an exponentially stretching vertical surface with magnetic field and viscous dissipation

The steady mixed convection flow and heat transfer from an exponentially stretching vertical surface in a quiescent Maxwell fluid in the presence of magnetic field, viscous dissipation and Joule heating have been studied. The stretching velocity, surface temperature and magnetic field are assumed to have specific exponential function forms for the existence of the local similarity solution. The coupled nonlinear ordinary differential equations governing the local similarity flow and heat transfer have been solved numerically by Chebyshev finite difference method. The influence of the buoyancy parameter, viscous dissipation, relaxation parameter of Maxwell fluid, magnetic field and Prandtl number on the flow and heat transfer has been considered in detail. The Nusselt number increases significantly with the Prandtl number, but the skin friction coefficient decreases. The Nusselt number slightly decreases with increasing viscous dissipation parameter, but the skin friction coefficient slightly increases. Maxwell fluid reduces both skin friction coefficient and Nusselt number, whereas buoyancy force enhances them.     

Convection of a heat-generating fluid in a rotating horizontal cylinder

Thermal convection of a fluid in a horizontal cylinder rotating about its own axis with uniformly volume-distributed internal heat sources is experimentally investigated. The enclosure boundary temperature was kept constant. The threshold of the excitation of convective flows and their structure are studied as functions of the heat-release intensity and the rotation velocity. The experiments are performed with water and water-glycerin solutions. It is shown that rapidly rotating fluid is in a stable quasiequilibrium state, namely, the temperature distribution is axisymmetric and has a maximum at the center of the enclosure. It is found that with decrease in the rotation velocity a convective flow arises thresholdwise, in the form of vortex cells periodically arranged along the axis. The thermal convection in the rotating enclosure is shown to be determined by the effects of two different mechanisms. One of these is due to the centrifugal force of inertia and plays the stabilizing role, while the other, thermovibrational mechanism is connected with nonisothermal fluid oscillations under the action of gravity in the enclosure-fitted reference frame and is responsible for the occurrence of mean thermal convection. The boundaries of the convection generation are plotted in the plane of the governing dimensionless parameters and the heat transfer in the supercritical region is studied.     

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation(DNS) of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed.The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers(P r = 0.71,1.5,and 3.0) and a shear Reynolds number(Reτ = 180).Some typical thermal statistics,including normalized mean temperature and their fluctuations,turbulent heat fluxes,Nusselt number and so on,are analyzed.The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number.Two reasons can explain this.First,the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases,and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field.Second,the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number,and the thermal feedback of the particles to turbulence becomes weak.     

Statistical analysis of instantaneous turbulent heat transfer in circular pipe flows

Turbulent heat transfer in circular pipe flow with constant heat flux on the wall is investigated numerically via Large Eddy Simulations for frictional Reynolds number Re _τ  = 180 and for Prandtl numbers in the range 0.1 ≤ Pr ≤ 1.0. In our simulations we employ a second-order finite difference scheme, combined with a projection method for the pressure, on a collocated grid in cylindrical coordinates. The predicted statistical properties of the velocity and temperature fields show good agreement with available data from direct numerical simulations. Further, we study the local thermal flow structures for different Prandtl numbers. As expected, our simulations predict that by reducing the Prandtl number, the range of variations in the local heat transfer and the Nusselt number decrease. Moreover, the thermal flow structures smear in the flow and become larger in size with less sharpness, especially in the vicinity of the wall. In order to characterize the local instantaneous heat transfer, probability density functions (PDFs) for the instantaneous Nusselt number are derived for different Prandtl number. Also, it is shown that these PDFs are actually scaled by the square root of the Prandtl number, so that a single PDF can be employed for all Prandtl numbers. The curve fits of the PDFs are presented in two forms of log-normal and skewed Gaussian distributions.     

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.     

Heat transfer in turbulent pipe flow based on a new mixing length model

A new model for the heat transfer in turbulent pipe flow is presented based on a modified form of the mixing length theory developed by Cebeci [1] for boundary layer flow problems. The model predicts the velocity and temperature distributions and the Nusselt number for fluids with low, medium and high Prandtl numbers (Pr=.02 to 15) and fits the available experimental data very accurately for values of Reynolds number exceeding 10⁴. Expressions for the eddy conductivity and for the turbulent Prandtl number are presented and shown to be dependent upon the Reynolds number, the Prandtl number, and the distance from the tube wall.     

Transient free convection flow of an incompressible viscous dissipative fluid

A finite-difference solution of transient free convection flow of a viscous dissipative fluid past an infinite vertical plate, on taking into account viscous dissipative heat is presented. Velocity profiles, temperature profiles are shown for different values of Pr, the Prandtl number and E, the Eckert number. The numerical values of the skin-friction and the rate of heat transfer are entered in a Table. It is observed that greater viscous dissipative heat causes a rise in the velocity, temperature and the skin-friction and a fall in the rate of heat transfer. An increase in Pr leads to a fall in the velocity, temperature and the skin-friction but the rate of heat transfer increases with increasing Pr.