Turbulent natural convection of sodium in a cylindrical enclosure with multiple internal heat sources: A conjugate heat transfer study

Turbulent natural convection of liquid metal in a cylindrical enclosure with locally distributed heat source has been investigated for Boussinesq number in the range of 4.5 × 0¹¹– 6.25 × 10¹². The enclosure considered is an ideal model of the lower plenum of a fast reactor with tray(s) holding damaged core debris, which continuously generate heat. The focus of the study has been to assess the heat dissipation capacity of single and multiple trays in respecting the specific temperature limits on the tray(s). Heat conduction in the metallic trays and in the impervious core debris mixture and turbulent natural convection of the liquid sodium are solved as a conjugate heat transfer problem. The equations that govern the various heat transfer processes in 2-D axi-symmetric cylindrical polar coordinate system have been solved by the finite volume method. Turbulence has been modeled by the k–ε turbulence model, without the use of wall functions. The predictions of the numerical model have been validated against benchmark data reported in open literature. Also, experiments have been conducted in an ideal water model towards validation of the computational model. For typical enclosure dimensions representing a 500 MWe fast reactor, it is seen that the critical parameters are heat dissipation area of the source (area of trays) and the thickness of the heat source (debris thickness). With increase in the number of trays, the heat transfer area increases while the debris thickness reduces. Both these effects lead to reduction in the tray as well as source temperatures. The heat dissipation capacity exhibits a non-linear relationship with the number of plates.     

Effects of thermal boundary conditions on natural convection flows within a square cavity

A numerical study to investigate the steady laminar natural convection flow in a square cavity with uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls has been performed. A penalty finite element method with bi-quadratic rectangular elements has been used to solve the governing mass, momentum and energy equations. The numerical procedure adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁵ and Prandtl number Pr, 0.7 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous Dirichlet boundary conditions. Non-uniform heating of the bottom wall produces greater heat transfer rates at the center of the bottom wall than the uniform heating case for all Rayleigh numbers; however, average Nusselt numbers show overall lower heat transfer rates for the non-uniform heating case. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes, power law correlations between average Nusselt number and Rayleigh numbers are presented.     

Numerical modeling of natural convection in an open cavity with two vertical thin heat sources subjected to a nanofluid

This paper presents a numerical study of natural convection cooling of two heat sources vertically attached to horizontal walls of a cavity. The right opening boundary is subjected to the copper–water nanofluid at 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 the Rayleigh number in the range 10⁴ ≤ Ra ≤ 10⁷, and for solid volume fraction 0 ≤ φ ≤ 0.05. In order to investigate the effect of heat source location, three different placement configurations of heat sources are considered. The effects of both Rayleigh numbers and heat source locations on the streamlines, isotherms, Nusselt number are investigated. The results indicate that the flow field and temperature distributions inside the cavity are strongly dependent on the Rayleigh numbers and the position of the heat sources. The results also indicate that the Nusselt number is an increasing function of the Rayleigh number, the distance between two heat sources, and distance from the wall. In addition it is observed that the average Nusselt number increases linearly with the increase in the solid volume fraction of nanoparticles.     

Geometric multigrid technique for solving heat convection-diffusion and phase change problems

Natural heat convection, mixed heat convection and heat transfer by conduction and convection with solidification of a ternary alloy are described by the finite volume method using a geometric multigrid approach. The objective of this paper is to analyze the effects of the multigrid technique on the accuracy and efficiency in describing convective heat transfer in closed and open cavities with and without liquid-solid phase changes of Newtonian and shear-thinning non-Newtonian fluids. It is found that the multigrid scheme reduces the computation time in natural convection in a square cavity from two times for Rayleigh number Ra = 10⁵ up to seven times for Ra = 10³, between 50% and 2.7 times for mixed convection with an inner solid when the Richardson number decreases from Ri = 10 to 0.1 and the Reynolds number Re = 100, and 20% for natural convection/heat conduction in solidification of a ternary aluminum alloy (Ra = 10⁴) with a shear-thinning rheology and a power index equal to 0.5.     

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.     

Prandtl number dependence of laminar natural convection heat transfer in a horizontal cylindrical enclosure with an inner coaxial triangular cylinder

A numerical study of laminar convection heat transfer from a horizontal triangular cylinder to its concentric cylindrical enclosure is performed to investigate the Prandtl number effect on flow and heat transfer characteristics. The Prandtl number over several orders of magnitude (10^?2 < Pr < 10³) as well as different aspect ratios (AR = 1.2 and 2.0) and different Rayleigh numbers (Ra = 10³, 10⁴, 10⁵, and 10⁶) are considered. The finite volume approach is used to solve the governing equations, in which buoyancy is modeled via the Boussinesq approximation. The computed flow patterns and temperature fields are shown by means of streamlines and isotherms, respectively, and the local and average heat transfer coefficients are also presented. It is found that the flow and heat transfer characteristics for a low Prandtl number fluid (Pr = 0.03) are unique and they are almost independent of Prandtl number when Pr ? 0.7. The entire spectrum of Prandtl number investigated can be divided into three sections based on the variations of average heat transfer coefficients. In each section, correlating equations of the average Nusselt number to the Rayleigh number are proposed with the maximum deviation less than 3%.     

Steady forced convection heat transfer from a heated circular cylinder to power-law fluids

Forced convection heat transfer to incompressible power-law fluids from a heated circular cylinder in the steady cross-flow regime has been investigated numerically by solving the momentum and thermal energy equations using a finite volume method and the QUICK scheme on a non-uniform Cartesian grid. The dependence of the average Nusselt number on the Reynolds number (5 ⩽ Re ⩽ 40), power-law index (0.6 ⩽ n ⩽ 2) and Prandtl number (1 ⩽ Pr ⩽ 1000) has been studied in detail. The numerical results are used to develop simple correlations as functions of the pertinent dimensionless variables. In addition to the average Nusselt number, the effects of Re, Pr and n on the local Nusselt number distribution have also been studied to provide further physical insights. The role of the two types of thermal boundary conditions, namely, constant temperature and uniform heat flux on the surface of the cylinder has also been presented.     

Investigation of natural convection induced outer side heat transfer rate of coiled-tube heat exchangers

Natural convection induced heat transfer has been studied over the outer surface of helically coiled-tube heat exchangers. Several different geometrical configurations (curvature ratio δ ε [0.035, 0.082]) and a wide range of flow parameters (60 <= T_tank <= 90, T_in = 19 and 60 <= T_in <= 90, T_tank = 20, 4000 <= Re <= 45000) have been examined to broaden the validity of the results gained from this research. A fluid-to-fluid boundary condition has been applied in the numerical calculations to create the most realistic flow configurations. Validity of the numerical calculations has been tested by experiments available in the open literature. Calculated results of the inner side heat transfer rate have also been compared to existing empirical formulas and experimental results to test the validity of the numerical computation in an independent way from the outer side validation of common helical tube heat exchangers. Water has been chosen to the working fluid inside and outside of the coiled tube (3 < Pr < 7). Outer side heat transfer rate along the helical tube axis has been investigated to get information about the performance of the heat transport process at different location of the helical tube. It was found that the outer side heat transfer rate is slightly dependent on the inner flow rate of any helical tube in case of increasing temperature differences between the tank working fluid temperature and the coil inlet temperature. A stable thermal boundary layer has been found along the axial direction of the tube.In addition to this the qualitative behavior of the peripherally averaged Nusselt number versus the axial location along the helical tube function is strongly dependent on the direction of the heat flow (from the tube to the storage tank and the reversed direction). Inner side heat transfer rate of helical coils have also been investigated in case of fluid-to-fluid boundary conditions and the calculation results have been compared with different prediction formulas published in the last couples of decades.     

Investigation of mixed convection heat transfer in a horizontal channel with discrete heat sources at the top and at the bottom

Mixed convection heat transfer from arrays of discrete heat sources inside a horizontal channel has been investigated experimentally. Each of the lower and upper surfaces of the channel was equipped with 8 × 4 flush mounted heat sources subjected to uniform heat flux. Sidewalls, lower and upper walls are insulated and adiabatic. The experimental parametric study was made for aspect ratios of AR = 2, 4 and 10, at various Reynolds and Grashof numbers. From the experimental measurements, row-average surface temperature and Nusselt number distributions of the discrete heat sources were obtained and effects of Reynolds and Grashof numbers on these numbers were investigated. From these results, the buoyancy affected secondary flow and the onset of instability have been discussed. Results show that top and bottom heater surface temperatures increase with increasing Grashof number. The top heater average-surface temperatures for AR = 2 are greater than those of bottom ones. For high values of Grashof numbers where natural convection is the dominant heat transfer regime (Gr¹/Re² ≫ 1), temperatures of top heaters can have much greater values. The variation of the row-average Nusselt numbers for the aspect ratio of AR = 4, show that with the increase in the buoyancy affected secondary flow and the onset of instability, values of Nusselt number level off and even rise as a result of heat transfer enhancement especially for low Reynolds numbers.     

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.     

Natural convective flow in an inclined lid-driven enclosure with a heated thin plate in the middle

Natural convection heat transfer from a heated thin plate located in the middle of a lid-driven inclined square enclosure has been analyzed numerically. Left and right of the cavity are adiabatic, the two horizontal walls have constant temperature lower than the plate’s temperature. The study is formulated in terms of the vorticity-stream function procedure and numerical solution was performed using a fully higher-order compact (FHOC) finite difference scheme on the 9-point 2D stencil. Air was chosen as a working fluid (Pr = 0.71). Two cases are considered depending on the position of heated thin plate (Case I, horizontal position; Case II, vertical position). Governing parameters, which are effective on flow field and temperature distribution, are Rayleigh number values (Ra) ranging from 10³ to 10⁵ and inclination angles γ (0° ? γ < 360°). The fluid flow, heat transfer and heat transport characteristics were illustrated by streamlines, isotherms and Nusselt number (Nu). It is found that fluid flow and temperature fields strongly depend on Rayleigh numbers and inclination angles. Further, for the vertical located position of thin plate heat transfer becomes more enhanced with lower γ at various Rayleigh numbers.     

Lattice Boltzmann simulation of natural convection heat transfer in an elliptical-triangular annulus

A numerical study for steady-state, laminar natural convection in a horizontal annulus between a heated triangular inner cylinder and cold elliptical outer cylinder was investigated using lattice Boltzmann method. Both inner and outer surfaces are maintained at the constant temperature and air is the working fluid. Study is carried out for Rayleigh numbers ranging from 1.0 × 10³ to 5.0 × 10⁵. The effects of different aspect ratios and elliptical cylinder orientation were studied at different Rayleigh numbers. The local and average Nusselt numbers and percentage of increment heat transfer rate were presented. The average Nusselt number was correlated. The results show that by decreasing the value of aspect ratio and/or increasing the Rayleigh number, the Nusselt number increases. Also the heat transfer rate increases when the ellipse positioned vertically.     

Natural convection in a rectangular enclosure containing an oval-shaped heat source and filled with Fe3O4/water nanofluid

This work concerns with the study of natural convection heat transfer in rectangular cavities with an inside oval-shaped heat source filled with Fe₃O₄/water nanofluid. The finite element method is employed to solve the governing equations for this problem. Average Nusselt numbers are presented for a wide range of Rayleigh number (10³ ≤ Ra ≤ 10⁵), volume fraction of nanoparticles (0 ≤ ϕ ≤ 14%), and four different size and shapes of the heat source. Depending on concentration of the nanoparticle, geometry of the heat source, and the value of Rayleigh number different behaviors are monitored for average Nusselt numbers. Configuration of the heat source dictates a significant change on the behavior of the average Nusselt number, while addition of the nanoparticles has a negative effect on the magnitude of Nusselt number for this problem.     

Experimental investigation of mixed convection heat transfer in a rectangular channel with discrete heat sources at the top and at the bottom

Mixed convection heat transfer in a top and bottom heated rectangular channel with discrete heat sources has been investigated experimentally for air. The lower and upper surfaces of the channel were equipped with 8 × 4 flush-mounted heat sources subjected to uniform heat flux. Sidewalls, the lower and upper walls were insulated and adiabatic. The experimental study was made for an aspect ratio of AR = 6, Reynolds numbers 955 ≤ Re_Dh ≤ 2220 and modified Grashof numbers Gr* = 1.7 × 10⁷ to 6.7 × 10⁷. From experimental measurements, surface temperature and Nusselt number distributions of the discrete heat sources were obtained for different Grashof numbers. Furthermore, Nusselt number distributions were calculated for different Reynolds numbers. Results show that surface temperatures increase with increasing Grashof number. The row-averaged Nusselt numbers first decrease with the row number and then, due to the increase in the buoyancy affected secondary flow and the onset of instability, they show an increase towards the exit as a result of heat transfer enhancement.     

Assessment of dissipative particle dynamics to simulate combined convection heat transfer: Effect of compressibility

The current study explored the capability of a discrete particle method known as dissipative particle dynamics with energy conservation (eDPD) to simulate combined convection heat transfer in a vertical lid driven cavity. The study investigated two cases of aiding and opposing buoyancy mechanisms in the lid driven cavity. The eDPD results were compared against the finite volume solutions for the range of Richardson number, 10^− 2 ≤ Ri ≤ 10². The method showed good comparison for the range of Richardson number 10^− 2 ≤ Ri ≤ 10¹. However, the eDPD method showed deviation from the FV solutions for a high value of Richardson number, Ri = 10², and this deviation is attributed to the compressibility of eDPD system experienced at such high value of Richardson number. Parametric study on the influence of the Richardson number (Ri) on the eDPD compressibility was conducted and presented via temperature isotherms, streamlines, velocity contours, velocity vectors, temperature and velocity profiles.     

Turbulent large-scale structures in natural convection vertical channel flow

We analyzed a database of a direct numerical simulation of natural convection in a vertical channel. The flow is driven by a constant temperature difference imposed at the walls (Ra = 5.4 × 10⁵, Pr = 0.7). The averaged flow and turbulent statistics are in good agreement with previous direct numerical simulations reported in the literature. Contrary to forced convection flows, the fluctuations of the heat transfer rate are uncorrelated with the fluctuations of the wall shear stress, which exhibit a symmetric probability density function. At the low Rayleigh number considered, the large-scale structures, which consist mainly in two counter-rotating vortices, with sizes comparable to the separation of the walls, are responsible for the extreme fluctuations of the wall heat transfer rate. The occurrence and the averaged topology of these structures have been determined using a conditional sampling technique.     

Numerical simulation of natural convection heat transfer in a porous cavity heated from below using a non-Darcian and thermal non-equilibrium model

The present paper investigates the numerical simulation of steady laminar incompressible natural convection heat transfer in an enclosed cavity that is filled with a fluid-saturated porous medium. The bottom wall is subjected to a relatively higher temperature than the top wall while the vertical walls are considered to be insulated. The flow field is modeled upon incorporating different non-Darcian effects, such as the convective term, Brinkman effect and Forchhiemer quadratic inertial effect. Moreover the two-equation model is used to separately account for the local fluid and solid temperatures. The numerical solution is obtained through the application of the finite volume method. The appraisals of the sought objectives are performed upon identifying key dimensionless groups of parameters. These dimensionless groups along with their operating domains are: Rayleigh number 1 ⩽ Ra ⩽ 400, Darcy number 10⁻⁴ ⩽ Da ⩽ 10⁻³, effective fluid-to-solid thermal conductivity ratio 0.1 ⩽ κ ⩽ 1.0, and the modified Biot number 1 ⩽ χ ⩽ 100. The non-Darcian effects are first examined over a broad range of Rayleigh number. Next, the implications of the group of parameters on the flow circulation intensity, local thermal non-equilibrium (LTNE) and average Nusselt number are highlighted and pertinent observations are documented.     

Numerical simulation of turbulent natural convection in a rectangular enclosure having finite thickness walls

Turbulent natural convection in a rectangular enclosure having finite thickness heat-conducting walls at local heating at the bottom of the cavity provided that convective-radiative heat exchange with an environment on one of the external borders has been numerically studied. Mathematical simulation has been carried out in terms of the dimensionless Reynolds averaged Navier–Stokes (RANS) equations in stream function–vorticity formulations. The formulation comprises the standard two equation k–ε turbulence model with wall functions, along with the Boussinesq approximation, for the flow and heat transfer. The special attention was paid to the effects of the Grashof number 10⁸ ? Gr < 10¹⁰, the transient factor 0 < τ < 1000 and the thermal conductivity ratio k_2,1 = 5.7 × 10^?4, 6.8 × 10^?5 both on local and on integral problem parameters. Detailed results including stream lines, temperature profiles and correlation for the average Nusselt number in terms of Grashof number have been obtained.     

Experimental investigation of the turbulent heavy liquid metal heat transfer in the thermal entry region of a vertical annulus with constant heat flux on the inner surface

The thermally developing turbulent lead bismuth (Pb⁴⁵Bi⁵⁵) flow within a vertical annulus with the inner surface uniformly heated is experimentally investigated. The ratio between internal and external diameter is 0.136. The investigated Reynolds numbers vary between 1.48 × 10⁴ and 2.37 × 10⁵ and the applied heat flux from 135 kW/m² to 905 kW/m², covering both the forced and mixed convection regime. For forced convection, the local heat transfer convective coefficient is evaluated and compared with literature correlations for low Prandtl number fluids. Additionally, a new empirical relation for the local Nusselt number in the thermal entry region is proposed. The applicability of criteria addressing the forced to mixed convection transition in medium to high Prandtl number fluids is reviewed and extended to account also for low Prandtl number fluids such as heavy liquid metals. Here, the transition from turbulent-diffusion to molecular-diffusion dominated mixed convection is clearly shown.     

Effects of moving lid direction on MHD mixed convection in a linearly heated cavity

Effects of moving lid-direction on MHD mixed convection in a cavity with the bottom wall being linearly heated are analyzed using a numerical technique. Vertical walls of the enclosure are adiabatic and the sliding wall at the top has constant temperature. The lid moves in the negative and positive x-direction. Finite volume method has been used to solve the governing equations. Results are presented for different values of Hartmann number (0 ? Ha ? 30), Reynolds number (100 ? Re ? 1000) and Grashof number (10⁴ ? Gr ? 10⁶). It is found that direction of lid is more effective on heat transfer and fluid flow in the case of mixed convection than it is the case in forced convection. Heat transfer is also decreased with increasing of magnetic field for all studied parameters.