A Study on Natural Convection in a Cold Square Enclosure With Two Vertical Eccentric Square Heat Sources Using the IB–LBM Scheme

Author(s):  
S. M. Dash ◽  
S. Sahoo

In this article, the natural convection process in a two-dimensional cold square enclosure is numerically investigated in the presence of two inline square heat sources. Two different heat source boundary conditions are analyzed, namely, case 1 (when one heat source is hot) and case 2 (when two heat sources are hot), using the in-house developed flexible forcing immersed boundary–thermal lattice Boltzmann model. The isotherms, streamlines, local, and surface-averaged Nusselt number distributions are analyzed at ten different vertical eccentric locations of the heat sources for Rayleigh number between 103 and 106. Distinct flow regimes including primary, secondary, tertiary, quaternary, and Rayleigh–Benard cells are observed when the mode of heat transfer is changed from conduction to convection and heat sources eccentricity is varied. For Rayleigh number up to 104, the heat transfer from the enclosure is symmetric for the upward and downward eccentricity of the heat sources. At Rayleigh number greater than 104, the heat transfer from the enclosure is better for downward eccentricity cases that attain a maximum when the heat sources are near the bottom enclosure wall. Moreover, the heat transfer rate from the enclosure in case 2 is nearly twice that of case 1 at all Rayleigh numbers and eccentric locations. The correlations for heat transfer are developed by relating Nusselt number, Rayleigh number, and eccentricity of the heat sources.

2020 ◽  
Author(s):  
Siva Subrahmanyam Mendu ◽  
P.K. Das

Abstract The present paper reports the numerical investigations for steady-state natural convection in power-law fluids inside a square enclosure embedded with bottom discrete heaters. The Lattice Boltzmann Method (LBM) is employed to model the flow and heat transfer phenomenon at different combinations of power-law index, Rayleigh number, and heat source length for a constant Prandtl number. The buoyancy force is incorporated in the collision term of the LBM through Boussinesq approximation. Simulation outcomes are furnished using streamlines and, temperature contours, velocity profiles and variation of heat transfer on the non-adiabatic walls to probe natural convection phenomena. The results indicate that the temperature and the flow fields in the enclosure are symmetric about the vertical centerline. The detailed physical interpretations have been provided for the reported results. Further, the increase in the power-law index means a rise in viscosity and a decrease in thermal energy transport for other constant parameters. The outcomes also specify that the intensity of circulation and heat transfer enhances with the increase of Rayleigh number and size of the localized heater. Finally, though, a rise in the size of the confined heat source enhances the rate of total thermal transport, it does not change the trend of fluid flow and local heat transfer rate.


2020 ◽  
Vol 847 ◽  
pp. 114-119
Author(s):  
Barbie Leena Barhoi ◽  
Ramesh Chandra Borah ◽  
Sandeep Singh

The present study relates to numerical investigation of natural convection heat transfer in a nanofluid filled square enclosure. One side of the enclosure is maintained at high temperature and the other side at a low temperature; while the top and bottom sides are adiabatic. The commercial CFD software ANSYS-FLUENT© was used to solve this numerical problem with the governing differential equations discretized by a control volume approach. nanofluids of Cu-water, Al2O3-water and TiO2-water have been simulated for a range of Rayleigh numbers and volume fractions. The results were obtained in the form of streamlines and isotherms. Interpretations of the results are done based on heat transfer rates, volume fraction, Rayleigh number and Nusselt number. It is to be noted that addition of nanoparticles enhances the heat transfer rate. It is also observed that the Nusselt number is highly affected by volume fraction and Rayleigh number.


2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Chithra Devaraj ◽  
Eswaramurthi Muthuswamy ◽  
Sundararaj Kandasamy

Natural convection heat transfer in a two-dimensional square enclosure at various angles of inclination is investigated numerically using a finite volume based computational procedure. The heat transfer is from a constant temperature heat source of finite length centred at one of the walls to the cold wall on the opposite side while the remaining walls are insulated. The effect of area ratio of the heat sourceAfrom 0.2 to 1.0, Rayleigh number Ra from 103to 107, and angle of inclination of the enclosureθvarying from 0° to 360° on the flow field and heat transfer characteristics are investigated. Streamline and isothermal line patterns are found to be similar at low Rayleigh numbers whereas at high Rayleigh numbers the differences are significant due to the influence of the parameters considered. Average Nusselt number decreases drastically as the position of the heat source is moved above the horizontal centre line of the enclosure. Correlation of the average Nusselt number which depends on the parameters of interest is obtained in the general formCRamAn. The correlation coefficients are determined by multiple regression analysis for the entire range of Rayleigh numbers analysed and the values found by correlation equations are in good agreement with the numerical results.


Author(s):  
Lubhani Mishra ◽  
R. P. Chhabra

Abstract Laminar natural convection in Bingham plastic fluids has been investigated from two differentially heated cylinders arranged either one above the other or along the diagonal of the square enclosure. The coupled momentum and energy equations have been solved to elucidate the effect of Rayleigh number (104–106), Prandtl number (10–100), Bingham number (0.01 to Bnmax), and the gap between the two cylinders in terms of the geometric parameters (0 to −0.25 for vertical alignment and 0.15 to 0.35 for diagonal alignment) on the detailed structure of the flow field and the overall heat transfer characteristics of the system. New extensive results are visualized in terms of streamlines, isotherm contours, and variation of the local Nusselt number along various surfaces. Additional insights are developed by examining the shear-rate contours and the yield surfaces delineating the fluid-like and solid-like regions in the flow domain. At high values of the Bingham number, the average Nusselt number reaches its asymptotic value corresponding to the conduction limit. The increasing Rayleigh number promotes fluid-like behavior which promotes heat transfer. The augmentation in heat transfer depends on the volume of fluid participating in the buoyancy-induced flow. For the vertical arrangement, the average Nusselt number (for the heated cylinder) decreases a little as these are moved slightly away from the center of the enclosure, followed by an increase as the two cylinders approach one of the sidewalls; this is so even in the conduction limit. In contrast, when the two cylinders are arranged along the diagonal, the Nusselt number progressively decreases as the gap between the two cylinders increases. Finally, predictive correlations have been developed for the average Nusselt number and the limiting Bingham number thereby enabling their estimation in a new application.


2013 ◽  
Vol 135 (12) ◽  
Author(s):  
M. A. Hassan ◽  
Manabendra Pathak ◽  
Mohd. Kaleem Khan

Viscoplastic fluids are special kind of non-Newtonian materials which deform or flow only when applied stresses are more than a critical value known as yield stress. In this work, a numerical investigation of natural convection in a square enclosure filled with viscoplastic fluids has been reported. The enclosure has been partially heated from the bottom wall by a heating source and symmetrically cooled from both the side walls. The rheology of the viscoplastic fluids has been modeled with Bingham fluid model. A scaling analysis has been presented to establish the gross dependence of heat transfer on different values of operating parameters of the problem. The effects of yield stress of the fluid on heat and fluid transport inside the enclosure have been investigated for different values of temperature difference, across the hot and cold surfaces and also for different fluids. The effects of different lengths of heated zone on the flow phenomena and heat transfer characteristics have been investigated for three different values of the heated zones. All the important results have been expressed in terms of Bingham number (Bn), Rayleigh number (Ra), and Prandtl number (Pr). It has been observed that with the increase in Bingham number, the buoyancy induced fluid circulation and convection effect decreases inside the enclosure. For each Rayleigh number, there correspond a critical Bingham number for which the heat transfer inside the enclosure takes place solely by conduction mode. This critical value increases with the increase in Rayleigh number. For fixed value of Bingham number, i.e., fixed value of yield stress, the effects of Rayleigh number and heated length on heat transfer have been observed similar to the case of natural convection in Newtonian fluid. It has been also observed that at high Bingham number the effect of increase in Rayleigh number on average Nusselt number is lesser compared to the effect of increasing Rayleigh number at low Bingham number. Using the present numerical results, a correlation of average Nusselt number as a function of other nondimensional numbers has been established.


1991 ◽  
Vol 113 (3) ◽  
pp. 268-274 ◽  
Author(s):  
G. Refai Ahmed ◽  
M. M. Yovanovich

A numerical study is carried out to investigate the influence of discrete heat sources on natural convection heat transfer in a square enclosure filled with air. The enclosure has two vertical boundaries of height H; one of them is cooled at Tc and the other has discrete heat sources [isoflux (q = c) or isothermal (Th = c)]. The enclosure has two horizontal adiabatic boundaries of length L. Results are reported for 0 ≤ Ra ≤ 106, Pr = 0.72, A = 1, aspect ratio ε, the relative size of the heat source to the total height, lies in the range 0.25 ≤ ε ≤ 1 and the discrete heat sources are located at the top or the bottom of the enclosure. Verification of numerical results is obtained at Ra = 0 (conduction limit) with analytical conduction solutions. In addition, a comparison with experimental and numerical data is made which also shows good agreement. The relationships between both Nu, ΔNu (change of thermal conductance) and Ra based on scale length (the size of the heat source S divided by the aspect ratio A) are also investigated here. A relationship Nu and Ra, based on scale length obtained from analytical solutions is correlated as Nu = Nu(Ra, ε). In addition, extrapolation correlations of Nu over the very high range of Rayleigh numbers (Ra ≥ 108) are developed.


2007 ◽  
Vol 129 (11) ◽  
pp. 1592-1599 ◽  
Author(s):  
C. D. Sankhavara ◽  
H. J. Shukla

Numerical investigation is carried out for natural convection in square enclosures consisting of partitions protruding from the end walls with different thermal boundary conditions at the end walls and partitions. The vertical walls were maintained isothermal at different temperatures. The Rayleigh number varies from 104 to 106 and the Prandtl number is 0.71. The thickness of the partition is fixed and is equal to one-tenth of the width of the enclosure. Their nondimensional length (l∕H) varies from 0 (a nonpartitioned enclosure) to 0.5 (two separate enclosures). A good agreement was found between the results in the present study and those published previously. The partitions were found to significantly influence the convective heat transfer. The average Nusselt number is less in the presence of partitions, and it decreases with increasing partition length (l∕H) from 0 to 0.5.


1970 ◽  
Vol 39 (1) ◽  
pp. 1-7 ◽  
Author(s):  
Sumon Saha ◽  
Noman Hasan ◽  
Chowdhury Md Feroz

A numerical study has been carried out for laminar natural convection heat transfer within a two-dimensional modified square enclosure having a triangular roof. The vertical sidewalls are differentially heated considering a constant flux heat source strip is flush mounted with the left wall. The opposite wall is considered isothermal having a temperature of the surrounding fluid. The rest of the walls are adiabatic. Air is considered as the fluid inside the enclosure. The solution has been carried out on the basis of finite element analysis by a non-linear parametric solver to examine the heat transfer and fluid flow characteristics. Different heights of the triangular roof have been considered for the present analysis. Fluid flow fields and isotherm patterns and the average Nusselt number are presented for the Rayleigh numbers ranging from 103 to 106 in order to show the effects of these governing parameters. The average Nusselt number computed for the case of isoflux heating is also compared with the case of isothermal heating as available in the literature. The outcome of the present investigation shows that the convective phenomenon is greatly influenced by the inclined roof height. Keywords: Natural convection, triangular roof, Rayleigh number, isoflux heating. Doi:10.3329/jme.v39i1.1826 Journal of Mechanical Engineering, vol. ME39, No. 1, June 2008 1-7


2013 ◽  
Vol 136 (2) ◽  
Author(s):  
Yanwei Hu ◽  
Yurong He ◽  
Shufu Wang ◽  
Qizhi Wang ◽  
H. Inaki Schlaberg

An experimental and numerical investigation on natural convection heat transfer of TiO2–water nanofluids in a square enclosure was carried out for the present work. TiO2–water nanofluids with different nanoparticle mass fractions were prepared for the experiment and physical properties of the nanofluids including thermal conductivity and viscosity were measured. Results show that both thermal conductivity and viscosity increase when increasing the mass fraction of TiO2 nanoparticles. In addition, the thermal conductivity of nanofluids increases, while the viscosity of nanofluids decreases with increasing the temperature. Nusselt numbers under different Rayleigh numbers were obtained from experimental data. Experimental results show that natural convection heat transfer of nanofluids is no better than water and even worse when the Rayleigh number is low. Numerical studies are carried out by a Lattice Boltzmann model (LBM) coupling the density and the temperature distribution functions to simulate the convection heat transfer in the enclosure. The experimental and numerical results are compared with each other finding a good match in this investigation, and the results indicate that natural convection heat transfer of TiO2–water nanofluids is more sensitive to viscosity than to thermal conductivity.


Author(s):  
Mustapha Faraji ◽  
El Mehdi Berra

Abstract This paper reported the mathematical modeling and numerical simulation of natural convection flow of Cu/water nanofluid in a square enclosure using the lattice Boltzmann method (LBM). The cavity is heated from below by heat source and cooled by the top wall. The vertical walls are adiabatic. After validating the numerical code against the numerical and experimental data, simulations were performed for different Rayleigh numbers (104–0.5 × 107), nanoparticles volume fractions (0–8%), and cavity inclination angle (0 deg–90 deg). The effects of the studied parameters on the streamlines, on isotherms distributions within the enclosure, and on the local and average Nusselt numbers are investigated. It was found that heat transfer and fluid flow structure depend closely on the nanoparticle concentration. Results show differences in stream separation between a base fluid and the nanofluid. Also, adding small nanoparticles fractions, less than 6%, to the base fluid enhances the heat transfer for higher Rayleigh numbers and cavity inclination angle less than 30 deg. It is concluded that the optimal dilute suspension of copper nanoparticles can be applied as a passive way to enhance heat transfer in natural convection engineering applications.


Sign in / Sign up

Export Citation Format

Share Document