Heat flow visualization for natural convection in rhombic enclosures due to isothermal and non-isothermal heating at the bottom wall

Analysis has been carried out for the energy distribution and thermal mixing in steady laminar natural convective flow through the rhombic enclosures with various inclination angles, φ for various industrial applications. Simulations are carried out for various regimes of Prandtl (Pr) and Rayleigh (Ra) numbers. Dimensionless streamfunctions and heatfunctions are used to visualize the flow and energy distribution, respectively. Multiple flow circulations are observed at Pr = 0.015 and 0.7 for all φs at Ra = 10⁵. On the other hand, two asymmetric flow circulation cells are found to occupy the entire cavity for φ = 75° at higher Pr (Pr = 7.2 and 1000) and Ra (Ra = 10⁵). Heatlines are found to be parallel circular arcs connecting the cold and hot walls for the conduction dominant heat transfer at Ra = 10³. The enhanced convective heat transfer is explained with dense heatlines and convective loop of heatlines at Ra = 10⁵. Heatlines clearly demonstrate that the left wall receives heat from the bottom wall as heatlines directly connect both the walls whereas the convective heat circulation cells play lead role to distribute the heat along the right wall, especially for smaller φs. On the other hand, the heat flow is evenly distributed to both side walls at higher φs via convection as well as direct conductive transport. Significant convective heat transfer from the bottom hot wall to the left cold wall occurs for φ = 30° cavity whereas the heat transfer to the right cold wall is maximum for φ = 75° irrespective of Pr. Average Nusselt number studies also show that φ = 30° cavity gives maximum heat transfer rate from the bottom to left wall irrespective of Pr in isothermal heating case. On the other hand, enhanced thermal mixing occurs at φ = 75° for both isothermal and non-isothermal heating strategies except at Pr = 0.015 in isothermal heating case.     

Heatlines based natural convection analysis in tilted isosceles triangular enclosures with linearly heated inclined walls: effect of various orientations

Natural convection in isosceles triangular enclosures with various configurations (case 1 — inverted, case 2 — straight and case 3 — tilted) is studied via heatline analysis for linear heating of inclined walls. Detailed analysis and comparison for various base angles (φ = 45°, 60°) of triangular enclosures have been carried out for a range of fluids (Pr = 0.015 − 1000) within Ra = 10³ − 10⁵ using Galerkin finite element method. The heat flow distributions indicate conduction dominant heat transfer at low Ra (Ra = 10³) for case 1 and case 2 whereas in case 3, convective heat flow is observed due to high buoyancy force. As Ra increases, enhanced thermal mixing is observed at the core of the cavity. Wall to wall heat transfer occurs at walls AB and AC due to linear heating boundary condition in all the cases. Although the distributions of fluid flow and heat flow are qualitatively similar for φ = 45° and 60°, the intensity of fluid flow and heat flow decreases as φ increases. Strength of fluid flow and heat flow circulation cells is found to be higher in case 3 for identical parameters. Results show that upper side wall (AC) for case 3 exhibits higher heat transfer rates whereas heat transfer rates for walls AB and AC are the same for case 1 and case 2. Also Nu_AB is higher for case 2 followed by case 1 and case 3 at the middle portion of wall AB. Thus to achieve high heat transfer from fluid to wall at the central region, case 2 and case 3 configurations may be recommended at high Ra (Ra = 10⁵) and Pr, irrespective of φ.     

Analysis of Bejan’s heatlines on visualization of heat flow and thermal mixing in tilted square cavities

This article analyzes the detailed heat transfer phenomena during natural convection within tilted square cavities with isothermally cooled walls (BC and DA) and hot wall AB is parallel to the insulated wall CD. A penalty finite element analysis with bi-quadratic elements has been used to investigate the results in terms of streamlines, isotherms and heatlines. The present numerical procedure is performed over a wide range of parameters (10³ ? Ra ? 10⁵,0.015 ? Pr ? 1000,0° ? φ ? 90°). Secondary circulations cells are observed near corner regions of cavity for all φ’s at Pr = 0.015 with Ra = 10⁵. Two asymmetric flow circulation cells are found to occupy the entire cavity for φ = 15° at Pr = 0.7 and Pr = 1000 with Ra = 10⁵. Heatlines indicate that the cavity with inclination angle φ = 15° corresponds to large convective heat transfer from the wall AB to wall DA whereas the heat transfer to wall BC is maximum for φ = 75°. Heat transfer rates along the walls are obtained in terms of local and average Nusselt numbers and they are explained based on gradients of heatfunctions. Average Nusselt number distributions show that heat transfer rate along wall DA is larger for lower inclination angle (φ = 15°) whereas maximum heat transfer rate along wall BC occur for higher inclination angle (φ = 75°).     

Effect of a magnetic field on natural convection in an open cavity subjugated to water/alumina nanofluid using Lattice Boltzmann method

In this paper, the effect of a magnetic field on natural convection in an open enclosure which subjugated to water/alumina nanofluid using Lattice Boltzmann method has been investigated. The cavity is filled with water and nanoparticles of Al₂O₃ at the presence of a magnetic field. Calculations were performed for Rayleigh numbers (Ra = 10⁴–10⁶), volume fractions of nanoparticles (φ = 0,0.02,0.04 and .0.06) and Hartmann number (0 ≤ Ha ≤ 90) with interval 30 while the magnetic field is considered horizontally. Results show that the heat transfer decreases by the increment of Hartmann number for various Rayleigh numbers and volume fractions. The magnetic field augments the effect of nanoparticles at Rayleigh number of Ra = 10⁶ regularly. Just as the most effect of nanoparticles for Ra = 10⁴ is observed at Ha = 30, so the most influence of nanoparticles occurs at Ha = 60 for Ra = 10⁵.     

Investigation of natural circulation in cavities with uniform heat generation for different Prandtl number fluids

Natural convection in enclosures with uniform heat generation and isothermal side walls is studied here. For the rectangular enclosure, two-dimensional conservation equations are solved using SIMPLE algorithm. Parametric studies are conducted to examine the effects of orientation of the cavity, fluid properties (Pr number), and aspect ratio for Rayleigh numbers up to 10⁶. For a horizontal square cavity, the flow becomes periodically oscillating at Ra = 5 × 10⁴ and chaotic at Ra = 8 × 10⁵. With a slight increase in the inclination angle, the oscillations die and for inclination angles greater than 15⁰, the flow attain a steady state over a range of Ra. It is found that for tall cavities (aspect ratio > 1), the steady-state solution is obtained for all values of Ra considered here. However, for wide cavities (aspect ratio < 1), an oscillatory flow regime is observed. The maximum temperature within the cavity is calculated for the range of Ra, aspect ratio and Pr number. Correlations for the maximum cavity temperature is presented here. The values of critical Rayleigh number at which the convection sets in the rectangular cavity are also studied and two distinct criteria are determined to evaluate the critical Rayleigh number. Further, a three-dimensional simulation is performed for a cubic cavity. It is found that the steady state solutions are obtained for all Rayleigh number, except at Ra = 10⁶. This is in contrast to the predictions for a two-dimensional square cavity, which has an oscillatory zone from Ra = 5 × 10⁴ onwards.     

Natural convection effects on heat and mass transfer in a curvilinear triangular cavity

Finite element method is used in this study to analyze the effects of buoyancy ratio and Lewis number on heat and mass transfer in a triangular cavity with zig-zag shaped bottom wall. Buoyancy ratio is defined as the ratio of Grashof number of solutal and thermal. Inclined walls of the cavity have lower temperature and concentration according to zig-zag shaped bottom wall. Enclosed space consists mostly of an absorber plate and two inclined glass covers that form a cavity. Both high temperature and high concentrations are applied to bottom corrugated wall. Computations were done for different values of buoyancy ratio (?10 ? Br ? 10), Lewis number (0.1 ? Le ? 20) and thermal Rayleigh number (10⁴ ? Ra_T ? 10⁶). Streamlines, isotherms, iso-concentration, average Nusselt and Sherwood numbers are obtained. It is found that average Nusselt and Sherwood numbers increase by 89.18% and 101.91% respectively as Br increases from ?10 to 20 at Ra_T = 10⁶. Also, average Nusselt decreases by 16.22% and Sherwood numbers increases by 144.84% as Le increases from 0.1 to 20 at this Rayleigh number.     

Turbulent Rayleigh–Bénard convection of water in cubical cavities: A numerical and experimental study

Experimental measurements and numerical simulations of natural convection in a cubical cavity heated from below and cooled from above are reported at turbulent Rayleigh numbers using water as a convective fluid (Pr = 6.0). Direct numerical simulations were carried out considering the Boussinesq approximation with a second-order finite volume code (10⁷ ⩽ Ra ⩽ 10⁸). The particle image velocimetry technique was used to measure the velocity field at Ra = 10⁷, Ra = 7 × 10⁷ and Ra = 10⁸ and there was general agreement between the predicted time averaged local velocities and those experimentally measured if the heat conduction through the sidewalls was considered in the simulations.     

Numerical investigation of several physical and geometric parameters in the natural convection into trapezoidal cavities

Natural convection in trapezoidal cavities, especially those with two internal baffles in conjunction with an insulated floor, inclined top surface, and isothermal left-heated and isothermal right-cooled vertical walls, has been investigated numerically using the Element based Finite Volume Method (EbFVM). In numerical simulations, the effect of three inclination angles of the upper surface as well as the effect of the Rayleigh number (Ra), the Prandtl number (Pr), and the baffle’s height (H_b) on the stream functions, temperature profiles, and local and average Nusselt numbers has been investigated. A parametric study was performed for a wide range of Ra numbers (10³ ? Ra ? 10⁶) H_b heights (H_b = H¹/3, 2H¹/3, and H¹), Pr numbers (Pr = 0.7, 10 and 130), and top angle (θ) ranges from 10 to 20. A correlation for the average Nusselt number in terms of Pr and Ra numbers, and the inclination of the upper surface of the cavity is proposed for each baffle height investigated.     

The temperature statistics of a surfactant-covered air/water interface during mixed convection heat transfer and evaporation

Surface temperature fields were measured of an air/water interface where heat was transferred from the water to the air under mixed convection conditions. The interfacial temperature field was measured using an infrared (IR) camera for mean wind speeds ranging from 0 to 4.0 m/s, in 1.0 m/s increments. Statistics of these surface temperature fields, specifically, the root mean square (rms) and the skewness were obtained. Plots of the rms versus the heat flux showed linear behavior for low wind speeds (U = 0–3 m/s), and the skewness was also found to increase with heat flux for U = 0–3 m/s, although these data exhibited significant scatter. The scaled root mean square temperature was revealed to be governed by the ratio Ra^1/3/(Re^14/5Pr^1/3) where Ra is the Rayleigh number, Re¹ the Reynolds number based on water side friction velocity and Pr is the Prandtl number.     

Ignition properties of n-butane and iso-butane in a rapid compression machine

Autoignition delay times of n-butane and iso-butane have been measured in a Rapid Compression Machine in the temperature range 660–1010 K, at pressures varying from 14 to 36 bar and at equivalence ratios φ = 1.0 and φ = 0.5. Both butane isomers exhibit a negative-temperature-coefficient (NTC) region and, at low temperatures, two-stage ignition. At temperatures below ～900 K, the delay times for iso-butane are longer than those for the normal isomer, while above this temperature both butanes give essentially the same results. At temperatures above ～720 K the delay times of the lean mixtures are twice those for stoichiometric compositions; at T < 720 K, the equivalence ratio is seen to have little influence on the ignition behavior. Increasing the pressure from 15 bar to 30 bar decreases the amplitude of the NTC region, and reduces the ignition delay time for both isomers by roughly a factor of 3. In the region in which two-stage ignition is observed, 680–825 K, the duration of the first ignition stage decreases sharply in the range 680–770 K, but is essentially flat above 770 K. Good quantitative agreement is found between the measurements and calculations for n-butane using a comprehensive model for butane ignition, including both delay times in the two-stage region, with substantial differences being observed for iso-butane, particularly in the NTC region.     

Numerical study of turbulent double-diffusive natural convection in a square cavity by LES-based lattice Boltzmann model

Turbulent double-diffusive natural convection in a square cavity represents numerous important problems in practice as well as in fundamental. However up to date the study on it is quite sparse and most previous studies just focus on laminar regime. To the best knowledge of the present authors, only several k–? models were developed to investigate turbulent double-diffusive convection and there is no attempt to use Large Eddy Simulation (LES). In order to deepen our knowledge on turbulent double-diffusive convection in a square cavity, we propose a novel LES-based lattice Boltzmann (LB) model to simulate such turbulent convectional flow. Previous LES-based LB models can be recovered from the present model. We find that the symmetry of the fluid circulation becomes broken since the Rayleigh number Ra = 10⁸, although the asymmetry is more clear when Ra ? 10¹⁰. More important, in the present study we find the power-law relationship among the Nusselt (Nu), the ratio of buoyancy forces (N) and the Rayleigh number (Ra) still exists in turbulent regime. The formula among them can be concluded as Nu = a × (Ra × ∣1 ? N∣)^b + c. The values of parameters a, b and c are given in this work.     

Flow visualization of natural convection in a tall,air-filled vertical cavity

Natural convection of air in a tall vertical cavity was studied using a smoke patterns and interferometry. Experiments covered Rayleigh numbers of 4850 < Ra < 54,800 and aspect ratio A ≈ 40. Secondary cells were noted at Ra as low as 6228. The flow was stable at Ra < 10⁴. As Ra exceeded 10⁴ the flow became irregular, the core flow became increasingly unsteady and 3-D motion became evident. Interferometry showed that most of the temperature drop exists in boundary layers near the walls. The core is well mixed and of relatively uniform temperature with little or no vertical stratification.     

Steady natural convection flows in a square cavity with linearly heated side wall(s)

The present numerical study deals with natural convection flow in a closed square cavity when the bottom wall is uniformly heated and vertical wall(s) are linearly heated whereas the top wall is well insulated. Non-linear coupled PDEs governing the flow have been solved by penalty finite element method with bi-quadratic rectangular elements. Numerical results are obtained for various values of Rayleigh number (Ra) (10³ ⩽ Ra ⩽ 10⁵) and Prandtl number (Pr) (0.7 ⩽ Pr ⩽ 10). Results are presented in the form of streamlines, isotherm contours, local Nusselt number and the average Nusselt as a function of Rayleigh number.     

Natural convection in a square cavity filled with a porous medium: Effects of various thermal boundary conditions

Natural convection flows in a square cavity filled with a porous matrix has been studied numerically using penalty finite element method for uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls. Darcy–Forchheimer model is used to simulate the momentum transfer in the porous medium. The numerical procedure is adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁶, Darcy number Da, 10⁻⁵ ⩽ Da ⩽ 10⁻³, and Prandtl number Pr, 0.71 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous thermal boundary conditions. Numerical results are presented in terms of stream functions, temperature profiles and Nusselt numbers. Non-uniform heating of the bottom wall produces greater heat transfer rate at the center of the bottom wall than uniform heating case for all Rayleigh numbers but average Nusselt number shows overall lower heat transfer rate for non-uniform heating case. It has been found that the heat transfer is primarily due to conduction for Da ⩽ 10⁻⁵ irrespective of Ra and Pr. The conductive heat transfer regime as a function of Ra has also been reported for Da ⩾ 10⁻⁴. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes the power law correlations between average Nusselt number and Rayleigh numbers are presented.     

The effect of surface thermal radiation on entropy generation in an open cavity with natural convection

In this work, a numerical study on entropy generation in a square open cavity with natural convection and surface thermal radiation is presented. The overall continuity, momentum, and energy equations for the gas phase in the open cavity were solved numerically by means of the finite-volume method. Temperature-dependent fluid properties were considered. During the calculations, the values of the Rayleigh number (Ra) were set in the range of 10⁴–10⁶. The temperature difference between the hot wall and the bulk fluid (ΔT) was varied between 50 and 400 K, and was represented by a dimensionless temperature difference (φ) for the purpose of generalization of the present study. The results of this investigation indicate that surface thermal radiation increases the overall entropy generation rate between 33.52% and 560.87%, and thus cannot be neglected in the analysis of this type of system.     

Wide-optical bandgap with improved conductivity p-μc-Si:Ox:H films prepared by Cat-CVD

Boron-doped hydrogenated microcrystalline silicon oxide (p-μc-Si:O_x:H) films have been deposited using catalytic chemical vapor deposition (Cat-CVD). The single-coiled tungsten catalyst temperature (T_fil) was varied from 1850 to 2100 °C and films were deposited on glass substrates at the temperatures (T_sub) of 100–300 °C. Different catalyst-to-substrate distances of 3–5 cm and deposition pressures from 0.1 to 0.6 Torr were considered.Optical and electrical characterizations have been made for the deposited samples. The sample transmittance measurement shows an optical-bandgap (Eg_opt) variation from 1.74 to 2.10 eV as a function of the catalyst and substrate temperatures. One of the best window materials was obtained at T_sub=100 °C and T_fil=2050 °C, with Eg_opt=2.10 eV, dark conductivity of 3.0×10^?3 S cm^?1 and 0.3 nm s^?1 deposition rate.     

Unsteady conjugate natural convection in a square enclosure filled with a porous medium

Mathematical simulation of unsteady natural convection modes in a square cavity filled with a porous medium having finite thickness heat-conducting walls with local heat source in conditions of heterogeneous heat exchange with an environment at one of the external boundaries has been carried out. Numerical analysis was based on Darcy–Forchheimer model in dimensionless variables such as a stream function, a vorticity vector and a temperature. The special attention was given to analysis of Rayleigh number effect Ra = 10⁴, 10⁵, 10⁶, of Darcy number effect Da = 10^?5, 10^?4, 10^?3, ∞, of the transient factor effect 0 < τ < 1000 and of the heat conductivity ratio k_2,1 = 3.7 × 10^?2, 5.7 × 10^?4, 6.8 × 10^?5 on the velocity and temperature fields. The influence scales of the defining parameters on the average Nusselt number have been detected.     

Laminar natural convection of non-Newtonian power-law fluids between concentric circular cylinders

The two-dimensional steady-state natural convection of power-law fluids is studied numerically between two concentric horizontal cylinders with different constant temperatures. The governing equations are discretized using finite volume technique based on second order upwind and are solved using the SIMPLE algorithm. The effects of Rayleigh number (10³ ≤ Ra ≤ 10⁵) and Prandtl number (10 ≤ Pr ≤ 10³) on the dimensionless velocity and temperature are investigated for both pseudoplastic and dilatant fluids. Also the mean Nusselt number for various values of governing parameters is obtained and discussed. The results indicate that with increasing the power-law index from 0.6 to 1.4, the mean Nusselt number decreases. In the best case among the range of parameters considered here the heat transfer rate for pseudo-plastic fluid (n = 0.6) is 170% higher than the Newtonian one and for dilatant fluid (n = 1.4) the heat transfer rate is 43% lower than the Newtonian fluid. So the pseudoplastic and dilatant fluids are more efficient than Newtonian fluids for cooling and insulating purposes, respectively. It is shown that as the Rayleigh number increases the cooling effect of pseudoplastic fluid and the insulating effect of dilatant fluid become more pronounced.     

Treatment of turbulent heat fluxes with the elliptic-blending second-moment closure for turbulent natural convection flows

In this paper a comparative study on the treatment of the turbulent heat fluxes with the elliptic-blending second-moment closure for natural convection flows is presented. Three different cases for treating the turbulent heat fluxes are considered. Those are the generalized gradient diffusion hypothesis (GGDH), the algebraic flux model (AFM) and the differential flux model (DFM). These models are implemented in a computer code especially designed for an evaluation of turbulent models. Calculations are performed for turbulent natural convection flows in an 1:5 rectangular cavity (Ra = 4.3 × 10¹⁰) and in a square cavity with conducting top and bottom walls (Ra = 1.58 × 10⁹). The calculated results are compared with the available experimental data. The results show that the GGDH, AFM and DFM models produce sufficiently accurate solutions for the turbulent natural convection in an 1:5 rectangular cavity where the strength of the thermal stratification is weak in a central region of the cavity. However, the GGDH model produces very erroneous solutions for the turbulent natural convection in a square cavity with conducting walls where the Rayleigh number is relatively small and the thermally stratified region is dominant. The AFM and DFM produce very accurate solutions for both cases without invoking any numerical problems.     

Entropy generation in natural convection in a symmetrically and uniformly heated vertical channel

In this study numerical predictions of local and global entropy generation rates in natural convection in air in a vertical channel symmetrically heated at uniform heat flux are reported. Results of entropy generation analysis are obtained by solving the entropy generation equation based on the velocity and temperature data. The analyzed regime is two-dimensional, laminar and steady state. The numerical procedure expands an existing computer code on natural convection in vertical channels. Results in terms of fields and profiles of local entropy generation, for various Rayleigh number, Ra, and aspect ratio values, L/b, are given. The distributions of local values show different behaviours for the different Ra values. A correlation between global entropy generation rates, Rayleigh number and aspect ratio is proposed in the ranges 10³ ⩽ Ra ⩽ 10⁶ and 5 ⩽ L/b ⩽ 20.