Modeling of Unsteady Natural Convection for Double-Pipe in a Partially Cooled Enclosure

In this article, a model is developed for unsteady natural convection heat transfer and fluid flow in a partially cooled enclosure with a hollow cylinder through it. The right vertical wall of the enclosure is cooled partially. The location of the partial cooling is set up in three different configurations; namely, bottom (P ₁), middle (P ₂), and top (P ₃). A hollow cylinder is located at the middle of the enclosure to simulate a double-pipe heat exchanger. Three values of Grashof number are applied in this work, i.e., 10⁴, 10⁵ and 10⁶, and three lengths of the cooler, i.e., 0.2, 0.4 and 0.6. Finite element method was utilized to solve the unsteady dimensionless conservation equations of mass, momentum and energy. It is found that the length and location of cooler does not have a significant effect on the natural convection for the case of the low Grashof number. The maximum heat transfer rate is reached when the cooler is located at the middle of the vertical wall.     

Natural Convection in a Square Enclosure with Two Horizontal Cylinders

This study investigates natural convection in a cooled square enclosure with two inner heated circular cylinders with the same diameter. The centers of two equidiameter cylinders are placed at those of the lower and upper half of the enclosure, respectively. The immersed boundary method (IBM) to model the inner circular cylinders based on the finite volume method is used to study a two-dimensional natural convection for different Rayleigh numbers varying in the range of 10³ ≤ Ra ≤ 10⁵. The effect of the radius of inner circular cylinders in an enclosure on heat transfer and fluid flow at different Rayleigh numbers has been examined. As the Rayleigh number increases, the horizontal symmetry is broken and the asymmetry occurred from the smaller radius. As the radius decreases, the dependence of the convection on the Rayleigh number is considerable. The dependence of the Nusselt number on the radius and the Rayleigh number is presented.     

Numerical Investigation of the Entropy Generation Due to Natural Convection in a Partially Heated Square Cavity Filled With Nanofluids

A comprehensive numerical investigation has been carried out on the heat transfer performance and entropy generation within a rectangular cavity containing nanofluid. The cavity consists of two heat sources located on the bottom and a side wall. The effects of influential parameters including type and concentration of nanoparticles, radius of corner, width and thickness of heaters, heater distance from corners and aspect ratio of the enclosure were studied. The results showed that the Nusselt number enhanced by increasing the aspect ratio of the cavity, the distance of heaters from the corners, and concentration of nanoparticle and applying Cu as nanoparticle while it reduced by increasing the radius of the corner and the width and thickness of the heat sources. The entropy generation was found to be profoundly minimized by lowering the Rayleigh number. In addition, the entropy generation was attenuated by increasing the Eckert number, corner radius, the distance from the corner and concentration of nanoparticles and using Al₂O₃ as nanoparticle. On the other hand, increasing the aspect ratio of the cavity, width and thickness of the heaters augmented the entropy generation. Interestingly, the entropy generation of the system was lowered by just increasing the distance of one heater from the corner, whereas increasing the thickness and width of one heater resulted in larger entropy generation. This study provides valuable insight into the change in the amount of heat transfer and entropy by altering the geometry as well as fluid properties.     

Natural Convection in an Enclosure with a Partially Active Magnetic Field

Chaotic natural convection flow of a molten gallium in a square enclosure with the upper and lower surfaces being insulated was studied by two-dimensional numerical simulation. Constant temperatures are imposed along the left and right walls of the enclosure with a volumetrically heated enclosure. In addition, a nonuniform partially active magnetic field is applied in a vertical direction. The flux lines spread out into a fringing field so the effective cross-sectional area of the gap is larger than that of the pole face. A chaotic regime is considered under steady state boundary condition. This study was done for an internal Rayleigh number of 10⁷, external Rayleigh number of 10⁵, and Prandtl number of 0.024. The study covers various magnet pole effect widths of 1/4, 1/2, and 3/4 from enclosure width and the magnetic field strength ranges 0.0 ≤ B _o  ≤ 10 Tesla. The transport equations for continuity, momentum, and energy are solved. The numerical results are reported for the effect of the partially active magnetic field on the velocity vectors, counters of temperature, streamline, and heat transfer coefficient. The numerical study shows that a magnetic field is damping chaotic oscillation behavior and decreases the amplitude of oscillation. Also, at a certain magnetic field strength the chaotic flow tend to becomes periodic flow at certain amplitude and frequency, and at high magnetic field strength the flow in the square enclosure flow tends to become steady laminar flow with stable average Nusselt number values; so, the random oscillation behavior disappeared. The effect of a nonuniform magnetic field tends to push the fluid to flow away from magnetic field region.     

Natural Convection Heat Transfer Enhancement in a Square Cavity Periodically Cooled from Above

The present article reports numerical results of natural convection within an air filled square cavity with its horizontal walls submitted to different heating models. The temperature of the bottom horizontal surface (hot temperature) is maintained constant, while that of the opposite surface (cold temperature) is varied sinusoidally with time. The remaining vertical walls are considered adiabatic. The parameters governing the problem are the amplitude (0 ≤ a ≤ 0.8) and the period (τ ≥ 0.001) of the variable temperature, the Rayleigh number (10³ ≤ Ra ≤ 7 × 10⁶), and the Prandtl number (Pr = 0.71). In constant cooling conditions (a = 0), up to three different solutions (monocellular flow MF, bicellular vertical flow BVF, and bicellular horizontal flow BHF) are obtained. Their existence ranges are delineated and, in the limits of the existence range of each solution, the transitions observed are identified and described. In the variable cooling conditions, the effect of the amplitude and the period of the exciting temperature on fluid flow and heat transfer is examined in the case of the MF, and BHF for specific values of Ra. Results are presented in terms of Ψ _max (t), Ψ _min (t), Nu(t) and streamlines, heatlines, and isotherms during the evolutions of selected flow cycles. In comparison with the constant heating conditions, it is found that the variable cooling temperature could lead to a drastic change in the flow structure and the corresponding heat transfer, especially at specific low periods of the cold variable temperature. This leads to a resonance phenomenon characterized by an important increase in heat transfer by about 46.1% compared to the case of a constant cold temperature boundary condition.     

Lattice Boltzmann Method Based Simulation of Natural Convection in Partially Open Enclosure

A thermal lattice Boltzmann method‐based analysis was performed to numerically investigate the heat transfer by natural convection from an enclosure with a large vertical side opening. The height of the opening was less than the enclosure height and the vertical wall opposite to the opening was maintained at constant temperature. A parametric study was carried out for different values of Rayleigh number (Ra) ranging from 10³ to 10⁵ with air as the working fluid for three opening sizes and three opening locations. The Prandtl number was fixed at 0.71 and the enclosure aspect ratio was also fixed at 2 in all calculations. With Boussinesq approximation, the temperature distribution and stream functions in the enclosure were predicted. The profile of the normal velocity component at the opening location was determined. The opening size affects the stratification and recirculation pattern within the enclosure. The average Nusselt number at the heated wall was determined for all cases. © 2013 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.21110     

Natural Convection in a Square Enclosure with a Square Heat Source at Different Horizontal and Diagonal Eccentricities

The present article aims to extend a previous numerical study on the natural convection process in a square enclosure with a vertical eccentric square heat source (cylinder) [1]. Here, we investigated buoyancy-induced convective flow and heat transfer for horizontal and diagonal eccentric displacement in a square cylinder. Numerical studies are performed for Rayleigh numbers in the range 10³–10⁶, using our recently developed flexible forcing IB–thermal lattice Boltzmann method scheme [1 S. M. Dash, T. S. Lee, and H. Huang Natural Convection from an Eccentric Square Cylinder Using a Novel Flexible Forcing Ib-Lbm Method, Numer. Heat Transfer, Part A: Appl., vol. 65, pp. 531–555, 2013.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]]. Detailed analysis of isotherms, streamlines, and Nusselt number distribution as functions of Rayleigh number and eccentricity is provided in this paper.     

Entropy Generation Minimization Versus Thermal Mixing Due to Natural Convection in Differentially and Discretely Heated Square Cavities

Uniform temperature distribution is a key parameter in many thermal processing applications. A considerable amount of additional energy is used to enhance the fluid mixing in order to maintain the temperature uniformity, but that affects the overall efficiency of the process. In this article, an alternate approach is proposed for maintaining uniform temperature via various distributed/discrete heating strategies while maintaining the minimal entropy generation. The system of laminar natural convection in differentially and discretely heated square cavities filled with various materials (molten metals, air, aqueous solutions, oils) is considered, and finite element simulations are performed for a range of Rayleigh numbers (Ra = 10³–10⁵). Entropy generation is evaluated using finite-element basis sets for the first time in this work, and the derivatives at particular nodes are estimated based on the functions within adjacent elements. Analysis of entropy generation in each case is carried out and a detailed investigation of entropy production due to local heat transfer and fluid friction irreversibilities is presented. It is found that a high thermal mixing may not be the optimal situation for achieving uniform temperature distribution based on entropy production. A greater degree of temperature uniformity with moderate thermal mixing may correspond to minimum entropy generation with distributed heating. Further, based on entropy generation minimization approach, it has been thermodynamically established that the distributed heating methodology with multiple heat sources may be the energy efficient strategy for attaining adequate uniform temperature distribution with minimum entropy generation.     

含热源方腔的耦合自然对流换热的三维数值模拟

对具有内热源方腔的稳态层流耦合自然对流换热进行了三维的数值模拟，采用的模拟代码基于连续介质计算力学的开源库OpenFoam，解决了自然对流换热与固体传热的耦合问题。对外壁面为常温、方腔内充满含体积热源流体的自然对流计算结果表明，温度场、速度场与非耦合的工况有很大差异。     

Entropy Generation due to Combined Natural Convection and Thermal Radiation within a Rectangular Enclosure

ABSTRACT

The present work investigates entropy production due to coupled natural convection/radiation heat transfer phenomenon in an inclined rectangular enclosure, isothermally heated from the bottom side and isothermally cooled from the other sides. The discrete-ordinate method is used in modeling the radiative transport equation while the statistical narrow band correlated-k model is adopted to deduce the radiative properties of the medium. The influence of pertinent parameters such as aspect ratio, inclination angle and walls emissivities on entropy generation is studied. It is found that the volumetric entropy generation is reduced when increasing the inclination angle of the enclosure. Moreover, it is shown that the minimum entropy production due to radiation heat transfer in participating media occurs at aspect ratio equal to unity.     

Effect of Fluid Yield Stress on Natural Convection From Horizontal Cylinders in a Square Enclosure

In this study, laminar natural convection heat transfer to Bingham plastic fluids from two differentially heated isothermal cylinders confined in a square enclosure (with isothermal walls) has been investigated numerically. The governing partial differential equations have been solved over the ranges of the dimensionless parameters, namely, Rayleigh number, 10² to 10⁶, Prandtl number, 10 to 100, and Bingham number, 0.01 to 100, for seven locations of inner cylinders as ±0.25, ±0.2, ±0.1 and 0. These values correspond to the range of Grashof number varying from 10 to 10⁵. The detailed flow and temperature fields are visualized in terms of the streamlines and isotherm contours. Further insights are developed by examining the iso-shear rate contours and the yield surfaces delineating the fluid-like and solid-like regions. The corresponding heat transfer results are analyzed in terms of the distribution of the local Nusselt number along the cylinder surface together with its surface averaged value as functions of the Rayleigh number, Prandtl number, Bingham number, and positions of the cylinders. It is found that the average Nusselt number increases with the increasing values of the Rayleigh number and decreases with the increasing Bingham number. For sufficiently large values of the Bingham number, the average Nusselt number reaches its asymptotic value wherein heat transfer takes place solely by conduction. Based on the present numerical results, simple correlations for the prediction of the average Nusselt number and the limiting Bingham number have been developed. Also, a dimensionless criterion denoting the cessation of convection regime is outlined for this configuration.     

Natural Convection from a Heated Elliptic Cylinder with a Different Axis Ratio in a Square Enclosure

A detailed study about the free convection over a heated elliptic cylinder, placed at the center of a square cavity having cooled walls, is performed. Simulations are carried out for three Rayleigh numbers (10⁴, 10⁵, and 10⁶) and two cavity aspect ratios (CR = 2.5 and 5.0) for different axis ratio (AR). The effect of AR on fluid flow and heat transfer characteristics for varying Rayleigh number and cavity aspect ratio are analyzed. The influence of AR is phenomenal at higher Ra and lower CR. At higher Ra, thermal plumes are observed above the cylinder for different ARs. Bicellular vortices are formed at low Ra by changing CR. The surface-averaged Nusselt number (Nu _avg ) increases with increasing AR and Ra. The value of Nu _avg increases with decreasing CR, and a correlation for Nu _avg in terms of AR is obtained for each CR.     

Natural Convection in Power-Law Fluids in a Square Enclosure from Two Differentially Heated Horizontal Cylinders

Laminar natural convection has been numerically investigated from two differentially heated horizontal cylinders in a square enclosure filled with power-law fluids. Two basic configurations, namely, vertical- and diagonal-alignment of the cylinders at various locations have been considered. The coupled continuity, momentum and energy equations have been solved numerically to elucidate the effect of the Grashof number (10²–10⁴), Prandtl number (0.7–100) and power-law index (0.2–2) for a range of symmetric and asymmetric locations of the cylinders. The velocity and temperature fields are visualized in terms of streamlines, isothermal contours and plots of the local and average Nusselt number for different positions of the cylinders. The occurrence of the power-law index in the definitions of the Grashof and Prandtl numbers accentuates the interplay between the viscous, inertial and buoyancy forces thereby leading to nonlinearity in the observed trends. The presence of the dead zones coupled with the dominance of conduction under certain conditions strongly influences the overall heat transfer. All else being equal, it is possible to improve heat transfer for asymmetric positioning of the cylinders, especially at high values of the Prandtl number and Grashof number in shear-thinning fluids. A predictive correlation has been developed thereby enabling the estimation of the heat transfer coefficient in a new application in terms of the geometric and kinematic parameters.     

Laminar Natural Convection of Bingham Fluids in a Square Enclosure with Vertical Walls Subjected to Constant Heat Flux

In this study, two-dimensional steady-state simulations of laminar natural convection in square enclosures with vertical sidewalls subjected to constant heat flux have been carried out, where the enclosures are considered to be completely filled with a yield-stress fluid obeying the Bingham model. Yield stress effects on heat and momentum transport are investigated for nominal values of Rayleigh number (Ra) in the range 10³–10⁶ and a Prandtl number (Pr) range of 0.1–100. It is found that the mean Nusselt number Nu increases with increasing values of Rayleigh number for both Newtonian and Bingham fluids. However, Nu values obtained for Bingham fluids are smaller than that obtained in the case of Newtonian fluids with the same nominal value of Rayleigh number Ra due to weakening of convective transport. The mean Nusselt number Nu in the case of Bingham fluids is found to decrease with increasing Bingham number, and for large values of Bingham number Bn, the value settles to unity (Nu = 1.0) as heat transfer takes place principally due to thermal conduction. The Nu values for the vertical walls subjected to constant heat flux are smaller than the corresponding values in the same configuration with constant vertical wall temperatures (for identical values of nominal Rayleigh, Prandtl, and Bingham numbers). However, the value of Bingham number at which Nu approaches to unity remains the same for both constant wall temperature and constant wall heat flux configurations. It is demonstrated that for small values of Bingham number Nu increases with increasing Prandtl number, but the opposite behavior occurs for large values of Bingham number. New correlations are proposed for the mean Nusselt number Nu for both Newtonian and Bingham fluids for square enclosures with vertical walls subjected to constant heat flux, which are shown to satisfactorily capture the correct qualitative and quantitative behavior of Nu in response to changes in Ra, Pr, and Bn.     

Numerical Investigation on Natural Convection Inside the Side Ventilated Square Enclosure with Vertical Mid-Partition

A finite volume-based computational study of steady laminar natural convection inside the square enclosure with cold partition wall centrally placed on top and bottom is presented. The fluid considered is air with Prandtl number 0.71. Except the partition walls, all other walls were assumed as hot. The heights of inlet and outlet ports are constantly fixed as 20% of height of the enclosure. The height of the opening in the partition walls were 10%, 20%, and 30% of height of enclosure. The buoyancy-driven heat transfer mechanism inside the domain is influenced by the percentage of opening in the partition, Rayleigh number, and geometrical position of inlet and outlet. The cold partition walls attract the flow that subsequently influences the thermal modifications around the partition. The fluid flow and heat transfer were investigated for 81 cases with different positions of inlet and outlet and varying height of openings for Ra = 10³, 10⁴, and 10⁵. The cross flow between the vertical walls through the opening in the central partition wall was considered in this study. The formation of vortices and their sizes depends on the configuration of inlet and outlet ports. The higher temperature gradient occurs near the inlet and outlet port of vertical walls. Local Nusselt number is maximum just below the inlet for all cases. Due to the dominance of buoyancy forces, heat transfer rate increases when Rayleigh number increases for all the cases. The hydrodynamic block effect by the openings on the partition wall has significant effect on the velocity profile than on the heat transfer.     

Magnetohydrodynamic Natural Convection in a Square Enclosure with Four Circular Cylinders Positioned at Different Rectangular Locations

We deploy a finite volume numerical computation to investigate the two-dimensional hydromagnetic natural convection in a cooled square enclosure in the presence of four inner heated circular cylinders with identical shape. The inner circular cylinders are placed in a rectangular array with equal distance away from each other within the enclosure and moving along the diagonals of the enclosure. All the walls of the enclosure are kept isothermal with temperatures less than that of the cylinders. A uniform magnetic field is applied along the horizontal direction normal to the vertical wall. All solid walls are assumed electrically insulated. Simulations are performed for a range of the controlling parameters such as the Rayleigh number 10³ to 10⁶, Hartmann number 0 to 50, and the dimensionless horizontal and vertical distance from the center of a cylinder to center of another cylinder 0.3 to 0.7. The study specifically aims to understand the effects of the location of the cylinders in the enclosure on the magnetoconvective transport, when they moved along the diagonals of the enclosure. It is observed that the unsteady behavior of the flow and thermal fields at relatively larger Rayleigh numbers and for some cylinder position are suppressed by imposition of the magnetic field. The heat transfer strongly depends on the position of the cylinders and the strength of the magnetic field. Hence, by controlling the position of the objects and the magnetic field strength, a significant control on the hydrodynamic and thermal transport can be achieved.     

Entropy Generation Due to Natural Convection in a Partially Open Cavity with a Thin Heat Source Subjected to a Nanofluid

This article presents a numerical study of natural convection cooling of a heat source mounted inside the cavity, with special attention being paid to entropy generation. The right vertical wall is partially open and is subjected to copper–water nanofluid at a constant low temperature and pressure, while the other boundaries are assumed to be adiabatic. The governing equations have been solved using the finite volume approach, using SIMPLE algorithm on the collocated arrangement. The study has been carried out for a Rayleigh number in the range 10³ < Ra < 10⁶, and for solid volume fraction 0 <? <0.05. In order to investigate the effect of the heat source and open boundary location, six different configurations are considered. The effects of Rayleigh numbers, heat source and open boundary locations on the streamlines, isotherms, local entropy generation, Nusselt number, and total entropy generation are investigated. The results indicate that when open boundary is located up, the fluid flow augments and hence the heat transfer and Nusselt number increase and total entropy generation decreases.     

Natural Convection in a Square Enclosure by a Finite-Element,Penalty Function Method Using Primitive Fluid Variables

Abstract

Numerical solutions for the problem of laminar natural convection in a square enclosure using the penalty function, finite-element method are presented. Solutions are obtained for values of the Rayleigh number up to 10⁷ using primitive fluid variables, and the efficacy of the method is demonstrated through a qualitative and quantitative evaluation of the results. The simplicity and general applicability of the method are shown especially in the context of extensions to three-dimensional geometries and irregular computational grids.     

Laminar Free Convection in a Square Enclosure Driven by the Lorentz Force

A numerical study is presented of laminar free convection flow driven by magnetic forces. An external magnetic field with one spatially varying component is applied to an electrically conducting fluid in a square enclosure. This magnetically-driven flow is controlled by the intensity and the wave number of the applied magnetic forcing. In addition, when the enclosure is heated laterally in a non-zero gravity environment, the resulting buoyant forces may contribute or resist the magnetically-driven fluid motion. The present results show that a strong magnetic field can even reverse the buoyant flow. The circulation intensity of the flow and the heat transfer from the sidewalls is increased with increasing magnetic field or with decreasing magnetic Reynolds number. The wave number of the magnetic forcing is also an important parameter that determines the vortex patterns and, consequently, the convection heat transfer.     

Numerical Investigation of Unsteady Natural Convection in an Inclined Square Enclosure With Heat-Generating Porous Medium

The unsteady laminar natural convection in an inclined square enclosure with heat-generating porous medium whose heat varies by a cosine function is investigated by a thermal equilibrium model and the Brinkman–Darcy–Forchheimer model numerically, with the four cooled walls of closure as isothermal. The numerical code based on the finite-volume method has been validated by reference data before it was adopted. Influence of dimensionless frequency and inclination angle on heat transfer characteristics in a square enclosure, such as flow distribution, isotherm, averaged Nusselt number on each wall, and time-averaged Nusselt number, are discussed, with specified value for Rayleigh number = 10⁸, Darcy number = 10^?4, Prandtl number = 7, porosity = 0.4, and specific heat ratio = 1. It is found that when the internal heat source varies by cosine, the Nusselt numbers of the four walls oscillate with the same frequency as the internal heat source; however, phase difference occurs. Moreover, frequency has little impact on time-averaged Nusselt number of the four walls, which is different from the phenomenon discovered in natural convection with suitable periodic varying wall temperature boundary condition. Moreover, inclination angle plays an important role in the heat transfer characteristics of the walls studied.