Finite element simulation of magnetoconvection inside a sinusoidal corrugated enclosure with discrete isoflux heating from below

A numerical investigation of the steady magnetoconvection in a sinusoidal corrugated enclosure has been performed. In this analysis, two vertical sinusoidal corrugated walls are maintained at a constant low temperature whereas a constant heat flux source whose length is varied from 20 to 80% of the reference length of the enclosure is discretely embedded at the bottom wall. The Penalty finite element method has been used to solve the governing Navier–Stokes and energy conservation equation of the fluid medium in the enclosure in order to investigate the effect of discrete heat source sizes on heat transfer for different values of Grashof number and Hartmann number. The values of the governing parameters are the Grashof number Gr (10³ to 10⁶), Hartmann number Ha (0 to 100) and Prandtl number Pr (0.71). The present numerical approach is found to be consistent and the solution is obtained in terms of stream functions and isotherm contours.     

Unsteady natural convection within a differentially heated enclosure of sinusoidal corrugated side walls

Numerically investigation of natural convection within a differentially heated modified square enclosure with sinusoidally corrugated side walls has been performed for different values of Rayleigh number. The fluid inside the enclosure considered is air and is quiescent, initially. The top and bottom surfaces are flat and considered as adiabatic. Results reveal three main stages: an initial stage, a transitory or oscillatory stage and a steady stage for the development of natural convection flow inside the corrugated cavity. The numerical scheme is based on the finite element method adapted to triangular non-uniform mesh element by a non-linear parametric solution algorithm. Investigation has been performed for the Rayleigh number, Ra ranging from 10⁵ to 10⁸ with variation of corrugation amplitude and frequency. Constant physical properties for the fluid medium have been assumed except for the density where Boussinesq’s approximation has been considered. Results have been presented in terms of the isotherms, streamlines, temperature plots, average Nusselt numbers, traveling waves and thermal boundary layer thickness plots, temperature and velocity profiles. The effects of sudden differential heating and its consequent transient behavior on fluid flow and heat transfer characteristics have been observed for the range of governing parameters. The present results show that the transient phenomena are greatly influenced by the variation of the Rayleigh number with corrugation amplitude and frequency.     

Numerical Analysis of the Effects of Selected Geometrical Parameters and Fluid Properties on MHD Natural Convection Flow in an Inclined Elliptic Porous Enclosure with Localized Heating

Magnetohydrodynamic (MHD) natural convection flow and associated heat convection in an oriented elliptic enclosure has been investigated with numerical simulations. A magnetic field was applied to the cylindrical wall of the configuration, the top and bottom walls of the enclosure were circumferentially cooled and heated, respectively, while the extreme ends along the cross‐section of the elliptic duct were considered adiabatic. The full governing equations in terms of continuity, momentum, and energy transport were transformed into nondimensional form and solved numerically using finite difference method adopting Gauss–Seidel iteration technique. The selected geometrical parameters and flow properties considered for the study were eccentricity (0, 0.2, 0.4, 0.6, and 0.8), angle of inclination (0°, 30°, 60°, and 90°), Hartmann number (0, 25, and 50), Grashof number (10⁴, 10⁵, and 10⁶), and Darcy number (10^?3, 10^?4, and 10^?5). The Prandtl number was held constant at 0.7. Numerical results were presented by velocity distributions as well as heat transfer characteristics in terms of local and average Nusselt numbers (i.e., rate of heat transfer). The optimum heat transfer rate was attained at e value of 0.8. Also, the heat transfer rate increased significantly between the angles of inclination 58° and 90°. In addition, Hartmann number increased with decreased heat transfer rate and flow circulation. A strong flow circulation (in terms of velocity distribution) was observed with increased Grashof and Darcy numbers. The combination of the geometric and fluid properties therefore can be used to regulate the circulation and heat transfer characteristics of the flow in the enclosure.     

Effects of corrugation frequency and aspect ratio on natural convection within an enclosure having sinusoidal corrugation over a heated top surface

Natural convection of a two-dimensional laminar steady-state incompressible fluid flow in a modified rectangular enclosure with sinusoidal corrugated top surface has been investigated numerically. The present study has been carried out for different corrugation frequencies on the top surface as well as aspect ratios of the enclosure in order to observe the change in hydrodynamic and thermal behavior with constant corrugation amplitude. A constant flux heat source is flush mounted on the top sinusoidal wall, modeling a wavy sheet shaded room exposed to sunlight. The flat bottom surface is considered as adiabatic, while the both vertical side walls are maintained at the constant ambient temperature. The fluid considered inside the enclosure is air having Prandtl number of 0.71. The numerical scheme is based on the finite element method adapted to triangular non-uniform mesh element by a non-linear parametric solution algorithm. The results in terms of isotherms, streamlines and average Nusselt numbers are obtained for the Rayleigh number ranging from 10³ to 10⁶ with constant physical properties for the fluid medium considered. It is found that the convective phenomena are greatly influenced by the presence of the corrugation and variation of aspect ratios.     

Effect of Prandtl Number on Free Convection from Two Cylinders in a Square Enclosure

This study explores the effect of Prandtl number on the laminar natural convection heat transfer to Newtonian fluids in a square enclosure consisting of one hot circular cylinder and one cold circular cylinder. The walls of the square enclosure are maintained isothermal and at the same temperature as the cold cylinder and the fluid medium. The governing partial differential equations have been solved numerically over the following ranges of conditions: Grashof number, 10 to 10⁵; Prandtl number, 0.7 to 100 (or the range of Rayleigh numbers as 7 to 10⁷); and relative positioning of the cylinders, ?0.25 to 0.25. However, the ratio of the radius of the cylinder to the side of the enclosure is held fixed at 0.2. Extensive results on the streamline and isotherm contours, the local Nusselt number distribution, and the average Nusselt number are discussed to delineate the influence of Grashof and Prandtl numbers on them for a given location with respect to the horizontal center line. The surface-averaged Nusselt number shows a positive dependence on Grashof and Prandtl numbers for a fixed location of the two cylinders. The heat transfer results have been correlated as a function of the Rayleigh number and geometric parameters, thereby enabling its prediction in a new application.     

Free Convection in Confined Power-Law Fluids From Two Side-by-Side Cylinders in a Square Enclosure

Laminar free convection heat transfer in power-law fluids from two side-by-side cylinders (one hot and one cold) confined in a square duct has been studied numerically in the two-dimensional flow regime. For a fixed value of the ratio of cylinder radius to size of the enclosure, the effect of geometrical placement of the cylinders is studied on the resulting velocity and temperature fields in the laminar free convection regime by considering six asymmetric locations of the two cylinders. In particular, extensive results reported herein span the range of conditions of Grashof number, 10 to 10⁵; Prandtl number, 0.7 to 100, thereby yielding the range of the Rayleigh number as 7 to 10⁷; power-law index, 0.3 to 1.8; and the relative positions (dimensionless) of the cylinders with respect to the centerline, –0.25 to 0.25. The heat transfer characteristics are analyzed in terms of the local Nusselt number along the surfaces of the two cylinders and the enclosure walls. Overall, the average Nusselt number shows a positive dependence on both the Grashof number and the Prandtl number irrespective of the values of power-law index and relative positioning of the cylinders. Also, all else being equal, shear-thinning fluid behavior promotes heat transfer with reference to that in Newtonian fluids. When the two cylinders are situated close to the bottom wall, the rate of heat transfer is augmented with reference to that for the symmetric positioning of the cylinders along the horizontal mid-plane of the enclosure. Conversely, heat transfer deteriorates as the cylinders are located above the centerline of the enclosure. The present numerical results have been consolidated via the use of a modified Rayleigh number, thereby enabling the estimation of the average Nusselt number in a new application.     

Effect of eccentricity and radius ratio on fluid flow and heat transfer inside an eccentric semicircular enclosure

latroductionNatural convechon heat transfer has gainedconsiderable attention because of itS' numerousaPPlications in the areas of energy conservation, coolingo f electrical and electronic components, design of solarcoil~, bed exchangers, and many others. Heattransfer inside annular space, air-filled cavity or annularsector has wide aPPlication in many engineeringProblems. In our earlier work["n, we have shown theeffeCt of eccentricity on heat transfer and flow field forradius ratio R'=2.0 f…     

Conjugate natural convection in a square enclosure: effect of conduction in one of the vertical walls

Steady, laminar, natural convection flow in a square enclosure has been analysed numerically. One vertical wall of the enclosure is thick, with a finite thermal conductivity, while the other three walls are taken to be of zero thickness. The problem is conjugate and the main focus of the study is on examining the effect of conduction in the wall on the natural convection flow in the enclosure. Three separate models to account for the wall conduction are investigated : (i) the complete conjugate case in which conduction in the thick vertical wall is assumed to be fully two-dimensional; (ii) a one-dimensional model in which conduction in the wall is assumed to be in the horizontal direction only; and (iii) a lumped parameter approach which assumes the solid-fluid interface temperature to be uniform. A Boussinesq fluid with Prandtl number of 0.7 (air) and Grashof numbers ranging from 10³ to 10⁷ are considered. For Grashof number > 10⁵, the temperature distribution in the wall shows significant two-dimensional effects and the solid-fluid interface temperature is found to be quite non-uniform. This non-uniformity tends to make the flow pattern in the enclosure asymmetric. In the parametric range investigated, all three models predict nearly the same value for the overall heat transfer.     

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.     

Numerical simulation of natural convection heat transfer from a heated square cylinder in a square cavity filled with micropolar fluid

A numerical investigation has been performed to visualize the magnetohydrodynamic natural convective heat transfer from a heated square cylinder situated within a square enclosure subjected to nonuniform temperature distributions on the left wall. The flow inside the enclosure is unsteady, incompressible, and laminar and the working fluid is micropolar fluid with constant Prandtl number (Pr = 7). The governing equations of the flow problem are the conservation of mass, energy, and linear momentum, as well as the angular momentum equations. Governing equations formulated in dimensionless velocity and pressure form has been solved by Marker and Cell method with second-order accuracy finite difference scheme. Comprehensive verification of the utilized numerical method and mathematical model has shown a good agreement with numerical data of other authors. The results are discussed in terms of the distribution of streamlines and isotherms and surface-averaged Nusselt number, for combinations of Rayleigh number, Ra (10³–10⁶), Vortex viscosity parameter, K (0–5), and Ha parameter (0–50). It has been shown that an increase in the vortex viscosity parameter leads to attenuation of the convective flow and heat transfer inside the cavity.     

Magnetohydrodynamic Natural Convection Flow in an Enclosure with a Finite Length Heater Using the Differential Quadrature (DQ) Method

Analysis of heat and fluid flow transport due to natural convection and magnetohydrodynamic (MHD) flows in a square enclosure with a finite length heater has been performed using the differential quadrature (DQ) technique. The heater with constant heat flux is located on the bottom wall of the enclosure and isothermal boundary conditions are applied to the right vertical wall while the remaining walls are adiabatic. The effects of heater length (0.2 ≤ ? ≤ 0.8), heater location (0.1 ≤ c/L ≤ 0.9), and direction of magnetic force (0° ≤ φ ≤ 90°) for different values of Grashof (10³ ≤ Gr ≤ 10⁶) and Hartmann numbers (0 ≤ Ha ≤ 100) on the heat and fluid flow in the enclosure are studied. According to the results obtained, heat transfer reduces when increasing the Hartmann number. The rate of reduction is higher for high values of Grashof number. The heat transfer rate for the heater closer to the cold wall is considerably higher than the heaters far from the right wall.     

Lattice Boltzmann Simulation of Natural Convection in a Differentially Heated Square Enclosure Containing Heat Generating Low‐Prandtl Number Fluid

Lattice Boltzmann simulations were conducted for the free convective flow of a low‐Prandtl number (Pr = 0.0321) fluid with internal heat generation in a square enclosure having adiabatic top and bottom walls and isothermal side walls. The problem of free convection with volumetric heat source has represented itself in connection with advanced engineering applications, such as water‐cooled lithium–lead breeder blankets for nuclear fusion reactors and liquid metal sources of spallation neutrons for subcritical fission systems. A single relaxation time (SRT) thermal lattice Boltzmann method (LBM) was employed. While applying SRT, a D2Q9 model was used to simulate the flow field and temperature field. Results have been obtained for various Rayleigh numbers characterizing internal and external heating from 10³ to 10⁶. Flow and temperature fields in terms of stream function and isotherms in the enclosure were predicted for these cases. The temperature of the fluid in the enclosure was found higher than the heated wall temperature at high values of internal Rayleigh numbers. The internal heat generation affected the rate of heat transfer significantly as two convection loops are observed in the enclosure. The average Nusselt number at the heated and cold wall was determined for all the cases.     

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.     

NATURAL-CONVECTIVE HEAT TRANSFER IN A SQUARE CAVITY WITH TIME-VARYING SIDE-WALL TEMPERATURE

ABSTRACT

The natural-convective heat transfer in an inclined square enclosure is studied numerically. The top and bottom horizontal walls are adiabatic, and the right side wall is maintained at a constant temperature T ₀. The temperature of the opposing vertical wall varies by sine law with time about a mean value T ₀. The system of Navier–Stokes Equations in Boussinesq approximation is solved numerically by the control-volume method with SIMPLER algorithm. The enclosure is filled with air (Pr = 1) and results are obtained in the range of inclination angle 0° ≤ α ≤ 90° for two values of Grashof number (2 × 10⁵ and 3 × 10⁵). It can be noted that there is a nonzero time-averaged heat flux through the enclosure at α ≠ 0°. The dependencies of time-averaged heat flux on oscillation frequency and inclination angle are depicted. It is found that the maximal heat transfer corresponds to the values of inclination angle α = 54^○ and dimensionless frequency f = 20π for both Grashof numbers studied (2 × 10⁵ and 3 × 10⁵).     

HEAT TRANSFER FROM THE COLD WALL OF A SQUARE CAVITY TO THE HOT ONE BY OSCILLATORY NATURAL CONVECTION

ABSTRACT

The periodic natural convection in an inclined square enclosure is studied numerically. One of the cavity sides is maintained at a constant temperature, and the temperature of the opposite wall is varied by the sine law. The mean value of this temperature in a period was lower than the given constant temperature. The other two walls of the enclosure were adiabatic. The system of Navier–Stokes Equations in Boussinesq approximation is solved numerically by the control-volume method with the SIMPLER algorithm. Solutions are obtained for Grashof number equal to 5 × 10⁵ and Prandtl number equal to unity in a wide range of oscillation frequency variation. The heat transfer dependency on oscillation frequency is studied for different values of inclination angle. It is found that there is a possibility of heat transferring from the colder wall to the hotter one in the whole range of problem parameters. The possible causes of such a phenomenon are analyzed.     

Conjugate natural convection in an enclosure with a heat source of constant heat transfer rate

A numerical study of two-dimensional transient natural convection in a rectangular enclosure having finite thickness heat-conducting walls with a heat source of constant heat transfer rate located on the inner side of the left wall in conditions of convection–radiation heat exchange with an environment on one of the external boundaries has been performed. Mathematical simulation has been carried out in terms of the dimensionless variables such as stream function – vorticity – temperature. Stream function, vorticity and energy equations have been solved by finite difference numerical method. The relevant governing parameters were: the Grashof number from 10⁶ to 10⁸, the Prandtl number, Pr = 0.7 and the conductivity ratio. Detailed results including streamlines and temperature profiles have been obtained.     

Lattice Boltzmann simulation of natural convection in tall enclosures using water/SiO2 nanofluid

Natural convection in enclosures using water/SiO₂ nanofluid is simulated with Lattice Boltzmann method (LBM). This investigation compared with other numerical methods and found to be in excellent agreement. This study has been carried out for the pertinent parameters in the following ranges: the Rayleigh number of base fluid, Ra = 10³-10⁵, the volumetric fraction of nanoparticles between 0 and 4% and aspect ratio (A) of the enclosure between 0.5 and 2. The thermal conductivity of nanofluids is obtained on basis of experimental data. The comparisons show that the average Nusselt number increases with volume fraction for the whole range of Rayleigh numbers and aspect ratios. Also the effect of nanoparticles on heat transfer augments as the enclosure aspect ratio increases.     

Numerical study of mixed convective heat transfer in a lid‐driven enclosure with rotating cylinders

A numerical simulation is performed to characterize the mixed convective transport in a three‐dimensional square lid‐driven enclosure with two rotating cylinders. The top wall is moving in the positive x‐direction, and the bottom wall is at a higher fixed temperature compared with all other isothermal walls. Both cylinders are rotating in its own plane about their centroidal axis. On the basis of rotation of both cylinders in clockwise or counter‐clockwise directions, four rotational models are studied. Various controlling parameters considered in the present study are Grashof number (10 ³ < Gr < 10 ⁵), rotating speed of the cylinder (5 < ω < 50), and the Reynolds number based on top wall movement is fixed to 100. The effect of cylinder rotation on the heat transfer of bottom wall is reported with the help of streamlines, contour plots of z‐component of vorticity, averaged and local Nusselt number, ratios of secondary flow and drag coefficient. It is observed that the heat transfer at the bottom wall is substantially dependent on the rotational model and rotational speed of the cylinder.     

Transient natural convective heat transfer of a low-Prandtl-number fluid from a heated horizontal circular cylinder to its coaxial triangular enclosure

Transient natural convective heat transfer of liquid gallium, which has a Prandtl number of 0.023 at 310 K, from a heated horizontal circular cylinder to its coaxial triangular enclosure is studied numerically by employing the control volume method. Two orientations of the triangular cylinder are investigated and the Grashof number is varied from 10⁴ to 10⁷. Development of natural convection is presented by means of the evolutions of the average Nusselt number over the outer triangular wall. Temporal stages during the course of development are identified and demonstrated through representative snapshots of streamlines and isotherms. The time-averaged Nusselt number is scaled with Grashof number for both conduction- and convection-dominated regimes. It is found that by placing horizontally the top side of the triangular cylinder, the convective flow becomes more stable and the overall heat transfer is enhanced. In addition, pitchfork bifurcation is explored quantitatively and its onset times are predicted as well.     

Effects of inclination angle on conduction—natural convection in divided enclosures filled with different fluids

Laminar natural convection in inclined enclosures filled with different fluids was studied by a numerical method. The enclosure was divided by a solid impermeable divider. One side of partition of enclosure was filled with air and the other side had water. The enclosure was heated from one vertical wall and cooled from the other while horizontal walls were adiabatic. The governing equations which were written in stream function–vorticity form were solved using a finite difference technique. Results were presented by streamlines, isotherms, mean and local Nusselt numbers for different thermal conductivity ratios of solid impermeable material (plywood or concrete), inclination angle (0° ≤ ? ≤ 360°) and Grashof numbers (10³ ≤ Gr ≤ 10⁶). The code was validated by earlier studies, which are available in the literature on conjugate natural convection heat transfer. Analytical solutions were obtained for low Grashof numbers. Obtained results showed that both heat transfer and flow strength strongly depended on thermal conductivity ratio of the solid material of partition, inclination angle and Grashof numbers. The heat transfer was lower in the air side of the enclosure than that of the water side.