Natural convection cooling of a localised heat source at the bottom of a nanofluid-filled enclosure

This article presents a numerical study of natural convection cooling of a heat source embedded on the bottom wall of an enclosure filled with nanofluids. The top and vertical walls of the enclosure are maintained at a relatively low temperature. The transport equations for a Newtonian fluid are solved numerically with a finite volume approach using the SIMPLE algorithm. The influence of pertinent parameters such as Rayleigh number, location and geometry of the heat source, the type of nanofluid and solid volume fraction of nanoparticles on the cooling performance is studied. The results indicate that adding nanoparticles into pure water improves its cooling performance especially at low Rayleigh numbers. The type of nanoparticles and the length and location of the heat source proved to significantly affect the heat source maximum temperature.     

Non-linear approximations for natural convection in a horizontal annulus

Two different problems are proposed as approximations of the usual system modelling natural convection under the Oberbeck-Boussinesq assumptions. The error is evaluated by means of the norm of its gradient in the Hilbert space. The average Nusselt number is also estimated.     

Experimental investigation of turbulent natural convection flow in a channel

This paper presents the results of velocity measurements of natural convection in symmetrically heated vertical channel using the particle image velocimetry (PIV) system. Velocity measurements were conducted at three different sections on the horizontal plane to validate the flow two-dimensionality and at three different heights in the vertical plane to establish vertical mean velocity profiles. The results indicate a considerable influence of the Rayleigh number and aspect ratio on the mean velocity profile. The results also indicate significant diffusion rates of the vertical mean velocity component and normal Reynolds stresses towards the center of the channel.     

Natural convection flow under a magnetic field in an inclined square enclosure differentialy heated on adjacent walls

Steady, laminar, natural-convection flow in the presence of a magnetic field in an inclined square enclosure differentially heated along the bottom and left vertical walls while the other walls are kept isothermal was considered. The governing equations were solved numerically for the stream function, vorticity and temperature ratio using the differential quadrature method for various Grashof and Hartmann numbers, inclination angle of the enclosure and direction of the magnetic field. The orientation of the enclosure changes the temperature gradient inside and has a significant effect on the flow pattern. Magnetic field suppresses the convective flow and its direction also influences the flow pattern, causing the appearance of inner loops and multiple eddies. The surface heat flux along the bottom wall is slightly increased by clockwise inclination and reduced by half by the counterclockwise inclination. The surface heat flux along the upper portion of the left side wall is reversed by the rise of warmer fluids due to the convection currents for no inclination and clockwise inclination of the enclosure.     

Numerical simulation of laminar natural convection in a laterally heated vertical cylindrical enclosure: application to crystal growth

Laminar natural convection has been studied in a laterally heated vertical cylindrical enclosure with a free insulated surface and a centrally located constant temperature wall at the top. These conditions are a simplification of the conditions existing in a Czochralski crystal pulling system. The laminar, axisymmetric flow of a Newtonian, constant physical properties fluid under Boussinesq’s approximation has been considered. Governing equations in primitive variable form are solved numerically by control volume method. SIMPLE algorithm due to Patankar has been used for the numerical simulation. The effects of the constant wall heat flux boundary condition at the side wall have been investigated whereas the bottom wall is considered to be insulated. Streamlines and isotherms are presented for various Rayleigh numbers and Prandtl numbers. Heat flux vectors through the melt are plotted for selected cases. The axial velocity and temperature variations at different horizontal sections of the crucible have been presented graphically to explain the transport processes inside the crucible. It has been observed that in case of low Pr and high Ra, flow separation occurs at the vertical wall of the crucible which leads to an oscillatory flow as Ra increases. The investigation has been extended to the oscillatory regime of flow in the zone of supercritical Rayleigh numbers and some unsteady results are also presented. Finally a heat transfer correlation has been developed for steady-state case.     

Numerical simulation of two-dimensional convection: Role of boundary conditions

The problem of two-dimensional, periodic in the horizontal coordinate, convection of an incompressible fluid heated from below between two horizontal planes is considered. The problem is solved in two formulations: with (stress-)free and hard (no-slip) boundary conditions on the horizontal planes. It is shown that at small supercriticalities the two-dimensional convection calculation leads to more correct results with hard than with free boundary conditions. It is established that the difference between the free and hard conditions is most strongly manifested in the pulsations of the vertical velocity component, whereas the dependence of the Nusselt number and the pulsations of the horizontal velocity component on the boundary conditions is more weakly expressed.     

Laminar natural convection in a partially divided rectangular cavity at high Rayleigh number

Finite element predictions of two-dimensional laminar natural convection in a partially divided rectangular cavity at high Rayleigh number are presented. The walls are differentially heated, the horizontal surfaces are insulated and the cavity contains a partial vertical divider which is centrally located and whose height is varied. Detailed results are presented for an aluminium half-divider in water for Rayleigh number up to 10¹¹ and compared directly with recent experiments in a cavity of aspect ratio 1/2. The predicted flow and heat transfer are in good agreement with the measurements and confirm the existence of a high Rayleigh number regime with characteristic behaviour that differs significantly from that found at lower Rayleigh number. In addition, the effects of the divider height, the divider conductivity, the fluid Prandtl number and the cavity aspect ratio are studied. The results show that a direct simulation of the complex flow and heat transfer that occurs in partially divided cavities is possible for realistic physical conditions.     

Turbulent natural convection scaling in a vertical channel

Using direct numerical simulation (DNS) data, this study appraises existing scaling laws in literature for turbulent natural convection of air in a differentially heated vertical channel. The present data is validated using past DNS studies, and covers a range of Rayleigh number, Ra between 5.4 × 10⁵ and 2.0 × 10⁷. We then appraise and compare the various scaling laws proposed by Versteegh and Nieuwstadt, 1999, Hölling and Herwig, 2005, Shiri and George, 2008, George and Capp, 1979 with the profiles of the mean temperature defect, mean streamwise velocity, normal velocity fluctuations, temperature fluctuations and Reynolds shear stress. Based on the arguments of an inner (near-wall) and outer (channel centre) region, the data is found to support a minus one-third power law for the mean temperature in an overlap region. Using the inner and outer temperature profiles, an implicit heat transfer equation is obtained and we show that a correction term is non-negligible for the present Ra range when compared with explicit equations found in literature. In addition, we determined that the mean streamwise velocity and normal velocity fluctuations collapse in the inner region when using the outer velocity scale. We also find that the temperature fluctuations scale in inner coordinates, in contrast to the outer scaling behaviour reported in the past. Lastly, we show evidence of an incipient proportional relationship between friction velocity, u_τ, and the outer velocity scale, u_o, with increasing Ra.     

Conjugate conduction-natural convection in an enclosure with time-periodic sidewall temperature and inclination

The unsteady conjugate conduction-natural convection in enclosure is of great theoretical significance and is widely encountered in engineering applications in the areas of fluid dynamics and heat transfer. However, there are relatively few efforts to investigate the unsteady flow physics and heat transfer characteristics in the inclined enclosure of finite thickness walls. In the present work, this problem is numerically investigated by a high accuracy multidomain temporal-spatial pseudospectral method. The enclosure is filled with Boussinesq fluid and is bounded by four finite thickness and conductive walls; one of the vertical sidewall is exposed to time-periodic temperature environment while the opposite sidewall holds constant temperature; the top and bottom walls are assumed to be adiabatic. Particular efforts are focused on the effects of three types of influential factors: the wall thermophysical properties, the time-periodic temperature patterns and the inclination, and the time-periodic flow patterns and heat transfer characteristics are presented. Numerical results reveal that within the present parameter range, the heat transfer rate increases almost linearly with the thermal conductivity ratio and thermal diffusivity ratio but decreases with the inclination angle. Moreover, the heat transfer could be enhanced or weakened by selecting different temperature pulsating period in the case of finite thickness wall, while it is always enhanced if the walls are zero thickness. The back heat transfer and heat transfer resonance phenomena are observed, and their relationships with the time-periodic flow patterns and temperature distributions are analyzed. The findings are helpful to the understandings of the fluid flow and heat transfer mechanisms in the related enclosure configurations, and may be of engineering use in thermal design improvement.     

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.     

Free convection heat transfer from an isothermal wavy surface in a porous enclosure

The coupled streamfuction–temperature equations governing the Darcian flow and convection process in a fluid-saturated porous enclosure with an isothermal sinusoidal bottom sun face, has been numerically analyzed using a finite element method (FEM). No restrictions have been imposed on the geometrical non-linearity arising from the parameters like wave amplitude (a), number of waves per unit length (N), wave phase (Φ), aspect ratio (A) and also on the flow driving parameter Rayleigh number (Ra). The numerical simulations for varying values of Ra bring about interesting flow features, like the transformation of a unicellular flow to a multicellular flow. Both with increasing amplitude and increasing number of waves per unit length, owing to the shift in the separation and reattachment points, a row–column pattern of multicellular flow transforms to a simple row of multicellular flow. A cycle of n celluar and n+1 cellular flows, with the flow in adjacent cells in the opposite direction, periodically manifest with phase varying between 0 and 360°. The global heat transfer into the system has been found to decrease with increasing amplitude and increasing number of waves per unit length. Only marginal changes in the global heat flux are observed, either with increasing Ra or varying Φ. Effectively, sinusoidal bottom surface undulations of the isothermal wall of a porous enclosure reduces the heat transfer into the system. © 1998 John Wiley & Sons, Ltd.     

Double diffusive natural convection in a partially heated enclosure with Soret and Dufour effects

The effect of double-diffusive natural convection of water in a partially heated enclosure with Soret and Dufour coefficients around the density maximum is studied numerically. The right vertical wall has constant temperature θ_c, while left vertical wall is partially heated θ_h, with θ_h > θ_c. The concentration in right wall is maintained higher than left wall (C_c < C_h) for case I, and concentration is lower in right wall than left wall (C_h > C_c) for case II. The remaining left vertical wall and the two horizontal walls are considered adiabatic. Water is considered as the working fluid. The governing equations are solved by control volume method using SIMPLE algorithm with QUICK scheme. The effect of the various parameters (thermal Rayleigh number, center of the heating location, density inversion parameter, Buoyancy ratio number, Schmidt number, and Soret and Dufour coefficients) on the flow pattern and heat and mass transfer has been depicted. Comprehensive Nusselt and Sherwood numbers data are presented as functions of the governing parameters mentioned above.     

Natural convection from a vertical cylinder in a thermally stratified porous medium

The problem of non-Darcy natural convection adjacent to a vertical cylinder embedded in a thermally stratified porous medium has been analyzed. Nonsimilarity solutions are obtained for the case that the ambient temperature increases linearly with height of the cylinder. A generalized flow model was used in the present study to include the effects of the macroscopic viscous term and the microscopic inertial force. Also, the thermal dispersion effect is considered in the energy equation. Thus, the main aim of this work is to examine the effects of thermal stratification and non-Darcy flow phenomena on the free convection flow and heat transfer characteristics. It was found that the present problem depends on six parameters, namely, the local thermal stratification parameter ξ, the boundary effect parameter Bp, the modified Grashof number Gr^*, wall temperature exponent m, the curvature parameter ω, and the modified Rayleigh number based on pore diameter Ra _d . The impacts of these governing parameters on the local heat transfer parameter are discussed in great detail. Also, representative velocity and temperature profiles are presented at selected values of the thermal stratification parameter. In general, the local heat transfer parameter is increased with increasing the values of m, ω, and Ra _d ; while it is decreased with increasing the values of ξ, Bp, and Gr^*. Received on 19 May 1998     

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.     

Numerical estimation of mixed convection heat transfer in vertical helically coiled tube heat exchangers

A numerical investigation of the mixed convection heat transfer from vertical helically coiled tubes in a cylindrical shell at various Reynolds and Rayleigh numbers, various coil‐to‐tube diameter ratios and non‐dimensional coil pitches was carried out. The particular difference in this study compared with other similar studies is the boundary conditions for the helical coil. Most studies focus on constant wall temperature or constant heat flux, whereas in this study it was a fluid‐to‐fluid heat exchanger. The purpose of this article is to assess the influence of the tube diameter, coil pitch and shell‐side mass flow rate on shell‐side heat transfer coefficient of the heat exchanger. Different characteristic lengths were used in the Nusselt number calculations to determine which length best fits the data and finally it has been shown that the normalized length of the shell‐side of the heat exchanger reasonably demonstrates the desired relation. Copyright © 2010 John Wiley & Sons, Ltd.     

Experimental investigation of natural convection in a vertical rib-roughened channel with asymmetric heating

An experimental study in an open-ended vertical channel is carried out in order to describe the fluid dynamics and heat transfer of transient free convection inside a vertical rib-roughened channel asymmetrically heated at various uniform heat fluxes (650, 700, and 780 W/m²) corresponding to various modified Rayleigh numbers (3.65 × 10⁶, 3.93 × 10⁶ and 4.4 × 10⁶). Two ribs are symmetrically located on each wall. The investigations focused more specifically on the influence of the ribs positions inside the channel and the modified Rayleigh number (Ra^*) both in steady-state regime and during the transitional phase occurring just after the start of the heating on the flow structure and the heat transfer performance. The results showed the appearance of large-scale flow instabilities which will develop and propagate until the development of the pocket-like vortex (reversed flow). Also, the formation and breakup of recirculation eddies, vortex banishment, besides that a separation and shifting of the boundary layer from one wall to another are identified. The best position of the ribs for heat extraction depends on the magnitude of the Rayleigh number. In that case, the top position is the optimal position for the small and the moderate modified Rayleigh numbers.     

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.     

Natural convection in a vertical strip immersed in a porous medium

In this work, the conjugated heat transfer characteristics of a thin vertical strip of finite length, placed in a porous medium has been studied using numerical and asymptotic techniques. The nondimensional temperature distribution in the strip and the reduced Nusselt number at the top of the strip are obtained as a function of the thermal penetration parameter s, which measures the thermal region where the temperature of the strip decays to the ambient temperature of the surrounding fluid. The numerical values of this nondimensional parameter permits to classify the different physical regimes, showing different solutions: a thermally long behaviour, an intermediate transition and a short strip limit.