Conjugate natural convection from an array of discrete heat sources: part 2 — a numerical parametric study

Coupled conduction and natural convection transport within a discretely heated cavity have been investigated numerically. One vertical wall of the cavity is composed of discrete, isoflux heat sources mounted in a substrate of finite thermal conductivity. The opposite vertical wall and the horizontal walls are assumed to be isothermal and adiabatic, respectively. The governing steady-state partial differential equations for the fluid and solid region are solved simultaneously using a control volume formulation, coupled with an additive correction multigrid procedure that increases the convergence rate of the solution. The fluid Prandtl number and heater/fluid thermal conductivity ratio are fixed at 25 and 2350, respectively, corresponding to a dielectric fluid (FC-77) and heaters manufactured from silicon. With increasing modified Rayleigh number (10⁴ < Ra_Lz^* < 109), the cavity flow becomes more boundary layer-like along the vertical walls, and multiple fluid cells develop in the central region. Thermal spreading in the substrate increases with decreasing modified Rayleigh number and with increasing values of the substrate/fluid thermal conductivity ratio (10⁻¹ <- R_s ≤ 10³). For large R_s, the discrete heat sources lose their thermal identity, and the streamlines and isotherms resemble those associated with a differentially heated cavity. Thermal spreading in the substrate also has a significant effect on circulation in the cavity and on maximum surface temperatures.     

Numerical comparison of conjugate and non-conjugate natural convection for internally heated semi-circular pools

In this numerical investigation, a detailed comparison of the conjugate and non-conjugate natural convection within a semi-cylindrical cavity has been presented. The cavity is assumed to be filled with a fluid containing uniformly distributed internal heating sources. The bottom circular wall of the cavity is taken to be thick with finite conductive properties, while the top wall is considered to be isothermal. The Navier–Stokes and energy equations are solved numerically by using the SIMPLER algorithm. A Rayleigh number range from 3.2×10⁶ to 3.2×10¹¹ has been investigated and the effects of solid-to-fluid conductivity ratios of 1.0, 5.0 and 23.0 have been analysed. The present numerical results for a semi-circular cavity with entirely isothermal walls are compared with known results from the open literature. It was found that these results for the non-conjugate problem are in very good agreement. The present results for a conjugate cavity show a remarkable difference from the non-conjugate analysis. The average Nusselt number for the solid–fluid interface shows a decrease while the top wall average Nu number has increased. It has been concluded that these effects increase for a system with a low solid-to-fluid conductivity ratio. It is evident from the present conjugate results that the assumption of isothermal enclosing walls gives somewhat different results when the walls are thick and the solid-to-fluid conductivity ratio is small.     

The variation of heat transfer in a two‐sided lid‐driven differentially heated square cavity with nanofluids containing carbon nanotubes for physical properties of fluid dependent on temperature

The behavior of nanofluids containing cylindrical nanoparticles are investigated numerically inside a two‐sided lid‐driven differentially heated square cavity to gain insight into the convective recirculation and flow processes induced by a nanofluid. The physical properties of the base fluid such as viscosity, thermal conductivity and thermal expansion coefficient are, respectively, assumed to be temperature independent (taking the mean temperature of the left and right walls) and temperature dependent. A model is developed to analyze the behavior of nanofluids taking into account the nanoparticle volume fraction whereas the transport equations are solved numerically with finite volume approach using SIMPLEC algorithm. The left and right moving walls are maintained at different constant temperatures while the upper and bottom walls are thermally insulated. The directions of the moving walls were considered in a way that the force and natural convections aid each other. The governing parameter Richardson number was 0.1<Ri<50.0 but due to space constraints only the results for 0.1<Ri<10.0 from fluid flow are presented. It was found that the temperature dependency of physical properties at different Richardson numbers and nanoparticle volume fractions affects the fluid flow and heat transfer in the cavities. Finally, comparisons between the behaviors of the average Nusselt number at the left wall for two cases are presented. Copyright © 2011 John Wiley & Sons, Ltd.     

The effects of forced boundary conditions on flow within a cubic cavity using digital particle image thermometry and velocimetry (DPITV)

Buoyancy-driven convection within a cavity, whose sidewalls are heated and cooled, is a problem of great interest, because it has applications in heat transfer and mixing. Most studies to date have studied one of two cases: the steady-state case or the development of the transient flow as it approaches steady state. Our main concern was to study the response of the cavity to time-varying thermal boundary conditions. We therefore decided to observe the flow phenomena within a convection cavity under sinusoidal thermal forcing of the sidewalls. To map the flow properly, it is necessary to have simultaneous kinematic and thermal information. Therefore, the digital particle image thermometry and velocimetry (DPITV) is used to acquire data. Implementing this technique requires seeding the flow with encapsulated liquid crystal particles and illuminating a cross section of the flow with a sheet of white light. Extraction of the thermal and kinematic content is in two parts. For the first, the liquid crystals will reflect different colors of the visible spectrum, depending on the temperatures to which they are subjected. Therefore, calibrating their color reflection with temperature allows for the extraction of the thermal content. For the second part, the kinematic information is obtained through the use of a digital cross-correlation particle image velocimetry technique. With the use of DPITV, the flow within a convection cavity is mapped and studied under steady forcing and sinusoidally forced boundary conditions at the Brunt-Väisälä frequency. For the sinusoidally forced case, three cases are studied. In the first, the heating between the two walls is in phase. In the second, the heating between the two walls is 180° out of phase. In the third, the heating between the two walls is 90° out of phase. For steady forcing, the thermal plots show that the flow develops a linearly stratified profile within the center of the cell. At the sidewalls, however, owing to forcing, hot/cold thermal boundary layers develop at the left/right walls. These hot/cold thermal boundary layers then turn around the upper-left/lower-right corners and develop into intrusion layers that extend across the top and bottom walls. The vorticity and streamlines show that the bulk of the fluid motion is concentrated around the walls, whereas the fluid within the center of the cell remains stationary. For the sinusoidally forced cases, the thermal plots show the existence of many thermal “islands,” or pockets of fluid where the temperature is different with respect to its surroundings. The vorticity plots show that the center of the cell is mostly devoid of vorticity and that the vorticity is mainly confined to the sidewalls, with some vorticity at the top and bottom walls. For the 0° forcing, the streamlines show the development of two counterrotating rollers. For the 180° forcing, the streamlines show the development of only one roller. Finally, for the 90° forcing, the streamlines show the development of both a two-roller and a one-roller system, depending on the position within the forcing cycle.     

Non-homogeneous model for a side heated square cavity filled with a nanofluid

A side heated two dimensional square cavity filled with a nanofluid is here studied. The side heating condition is obtained by imposing two different uniform temperatures at the vertical boundary walls. The horizontal walls are assumed to be adiabatic and all boundaries are assumed to be impermeable to the base fluid and to the nanoparticles. In order to study the behavior of the nanofluid, a non-homogeneous model is taken into account. The thermophysical properties of the nanofluid are assumed to be functions of the average volume fraction of nanoparticles dispersed inside the cavity. The definitions of the nondimensional governing parameters (Rayleigh number, Prandtl number and Lewis number) are exactly the same as for the clear fluids. The distribution of the nanoparticles shows a particular sensitivity to the low Rayleigh numbers. The average Nusselt number at the vertical walls is sensitive to the average volume fraction of the nanoparticles dispersed inside the cavity and it is also sensitive to the definition of the thermophysical properties of the nanofluid. Highly viscous base fluids lead to a critical behavior of the model when the simulation is performed in pure conduction regime. The solution of the problem is obtained numerically by means of a Galerkin finite element method.     

Effect of an Inserted Porous Layer Located at a Wall of a Parallel Plate Channel on Forced Convection Heat Transfer

A theoretical study is performed on heat and fluid flow in partially porous medium filled parallel plate channel. A uniform symmetrical heat flux is imposed onto the boundaries of the channel partially filled with porous medium. The dimensional forms of the governing equations are solved numerically for different permeability and effective thermal conductivity ratios. Then, the governing equations are made dimensionless and solved analytically. The results of two approaches are compared and an excellent agreement is observed, indicating correctness of the both solutions. An overall Nusselt number is defined based on overall thermal conductivity and difference between the average temperature of walls and mean temperature to compare heat transfer in different channels with different porous layer thickness, Darcy number, and thermal conductivity ratio. Moreover, individual Nusselt numbers for upper and lower walls are also defined and obtained. The obtained results show that the maximum overall Nusselt number is achieved for thermal conductivity ratio of 1. At specific values of Darcy number and thermal conductivity ratio, individual Nusselt numbers approach to infinity since the value of wall temperatures approaches to mean temperature.     

Freezing of Water in a Differentially Heated Cubic Cavity

An experimental and numerical study has been made of transient natural convection of water freezing in a cube-shaped cavity. The effect of the heat transfer through the side walls is studied in two configurations: with the cavity surrounded by air and with the cavity immersed in an external water bath of constant temperature. The experimental data for the velocity and temperature fields are obtained using liquid crystal tracers. The transient development of the ice/water interface is measured. The collected data are used as an experimental benchmark and compared with numerical results obtained from a Finite-difference code with boundary fitted grid generation. The computational model has been adopted to simulate as closely as possible the physical experiment. Hence, fully variable fluid properties are implemented in the code, and, to improve modelling of the thermal boundary conditions, the energy equation is also solved inside the bounding walls. Although the general behaviour of the calculated ice front and its volume matches observations, several details of the flow structure do not. Observed discrepancies between experimental and numerical results indicate the necessity of verifying and improving the usual assumptions for modelling ice formation.     

Conjugate conduction-natural convection in an enclosure with time-periodic sidewall temperature and inclination

The unsteady conjugate conduction-natural convection in enclosure is of great theoretical significance and is widely encountered in engineering applications in the areas of fluid dynamics and heat transfer. However, there are relatively few efforts to investigate the unsteady flow physics and heat transfer characteristics in the inclined enclosure of finite thickness walls. In the present work, this problem is numerically investigated by a high accuracy multidomain temporal-spatial pseudospectral method. The enclosure is filled with Boussinesq fluid and is bounded by four finite thickness and conductive walls; one of the vertical sidewall is exposed to time-periodic temperature environment while the opposite sidewall holds constant temperature; the top and bottom walls are assumed to be adiabatic. Particular efforts are focused on the effects of three types of influential factors: the wall thermophysical properties, the time-periodic temperature patterns and the inclination, and the time-periodic flow patterns and heat transfer characteristics are presented. Numerical results reveal that within the present parameter range, the heat transfer rate increases almost linearly with the thermal conductivity ratio and thermal diffusivity ratio but decreases with the inclination angle. Moreover, the heat transfer could be enhanced or weakened by selecting different temperature pulsating period in the case of finite thickness wall, while it is always enhanced if the walls are zero thickness. The back heat transfer and heat transfer resonance phenomena are observed, and their relationships with the time-periodic flow patterns and temperature distributions are analyzed. The findings are helpful to the understandings of the fluid flow and heat transfer mechanisms in the related enclosure configurations, and may be of engineering use in thermal design improvement.     

Numerical and experimental study of heat transfer in a tall vertical closed cavity

In this work the numerical and experimental results of heat transfer in a vertical tall closed cavity are presented. The cavity has an aspect ratio of 20, one of the vertical walls receive a constant and uniform heat flux, while the opposite wall is kept at a constant temperature. The remaining walls are assumed adiabatic. The cavity is full of air. The computational fluid dynamics software Fluent 6.3 was used for the simulation and an experimental prototype was built to obtain the heat transfer coefficients. The air temperature and the fluid velocity values are higher when emissivity (ε) is 0.03 (almost pure natural convection). The experimental total heat transfer coefficient increases between 119.9 and 159.9 % when the emissivity of the walls changes from 0.03 to 0.95.     

Analytical approximation for solidification processes in PCM storage with internal fins: imposed heat flux

Uses of thermal energy storage systems have expanded notably in recent decades. In thermal energy systems, internal heat transfer enhancement techniques such as fins are often used because of the low thermal conductivity of the phase change materials (PCMs). In this paper, solidification of a PCM is studied in a rectangular storage with horizontal internal plate fins and an imposed constant heat flux on the vertical walls. A simplified analytical solution is presented and its results are compared to those for a numerical approach based on an enthalpy method. The fraction of solidified PCM in storage is calculated with the derived analytical model which determines how much of the storage is solidified after a certain time. The results show that the analytical model satisfactorily estimates the solid–liquid interface and the temperature distribution for the fin, which are useful in the design of PCM-based thermal energy storage or cooling systems.     

Effects of walls temperature variation on double diffusive natural convection of Al2O3–water nanofluid in an enclosure

The effect of wall temperature variations on double diffusive natural convection of Al₂O₃–water nanofluid in a differentially heated square enclosure with constant temperature hot and cold vertical walls is studied numerically. Transport mechanisms of nanoparticles including Brownian diffusion and thermophoresis that cause heterogeneity are considered in non-homogeneous model. The hot and cold wall temperatures are varied, but the temperature difference between them is always maintained 5 °C. The thermophysical properties such as thermal conductivity, viscosity and density and thermophoresis diffusion and Brownian motion coefficients are considered variable with temperature and volume fraction of nanoparticles. The governing equations are discretized using the control volume method. The results show that nanoparticle transport mechanisms affect buoyancy force and cause formation of small vortexes near the top and bottom walls of the cavity and reduce the heat transfer. By increasing the temperature of the walls the effect of transport mechanisms decreases and due to enhanced convection the heat transfer rate increases.     

Numerical study on heat/mass transfer by shear/thermal/solutal-convection in a cavity under counteracting buoyancies

This study provides physical insights into a lid-driven square cavity filled with a mixture of a solvent vapor and non-condensable gas, subjected to the vertically parallel thermal and solutal gradients. The top lid is maintained at constant speed while bottom lid and the other two walls are kept fixed. Zero heat and mass fluxes are imposed on the vertical side walls. The transport equations are solved numerically through a pressure-correction-based iterative algorithm (SIMPLE) with the QUICK scheme for convective terms. The diffusivities of heat and salt are assumed to be equal throughout this investigation. The essential details of flow, temperature and concentration fields are presented for the opposing buoyancy forces ratio (B < 0) with special attention being given for the values of parameters for which either the flow inside the cavity is operated by the mechanically induced convection; or the flow structure inside the cavity is akin to a single-diffusive thermal or solutal convection. The variations of average rates of heat and mass transfer are uniform with the Reynolds number, while the variations of these quantities against the solutal Richardson number (Ri _C ) and thermal Richardson number (Ri _T ) point to the existence of the local minimum/maximum. Finally, two linear relations between Ri _C and Ri _T at a constant sliding speed are proposed to identify the above points of high and low transport phenomena and justified with the exhibition of the flow structures inside the cavity.     

The effects of variable properties on MHD unsteady natural convection heat and mass transfer over a vertical wavy surface

A mathematical model will be analyzed in order to study the effects of variables viscosity and thermal conductivity on unsteady heat and mass transfer over a vertical wavy surface in the presence of magnetic field numerically by using a simple coordinate transformation to transform the complex wavy surface into a flat plate. The fluid viscosity is assumed to vary as a exponential function of temperature and thermal conductivity is assumed to vary linearly with temperature. An implicit marching Chebyshev collocation scheme is employed for the analysis. Numerical solutions are obtained for different values of variable viscosity, variable thermal conductivity and MHD variation parameter. Numerical results show that, variable viscosity, variable thermal conductivity and MHD variation parameter have significant influences on the velocity, temperature and concentration profiles as well as for the local skin friction, Nusselt number and Sherwood number.     

Natural convection of an alumina-water nanofluid inside an inclined wavy-walled cavity with a non-uniform heating using Tiwari and Das’ nanofluid model

The present study is devoted to numerical analysis of natural convective heat transfer and fluid flow of alumina-water nanofluid in an inclined wavy-walled cavity under the effect of non-uniform heating. A single-phase nanofluid model with experimental correlations for the nanofluid viscosity and thermal conductivity has been included in the mathematical model. The considered governing equations formulated in dimensionless stream function, vorticity, and temperature have been solved by the finite difference method. The cavity inclination angle and irregular walls (wavy and undulation numbers) are very good control parameters for the heat transfer and fluid flow. Nowadays, optimal parameters are necessary for the heat transfer enhancement in different practical applications. The effects of the involved parameters on the streamlines and isotherms as well as on the average Nusselt number and nanofluid flow rate have been analyzed. It has been found that the heat transfer rate and fluid flow rate are non-monotonic functions of the cavity inclination angle and undulation number.     

A Local Model for the Energy Transfer in a Saturated Static Mixture

In the present work the transient energy transfer in a nonsaturated porous medium is studied, using a mixture theory viewpoint. The porous matrix is assumed homogeneous, rigid and isotropic, while the fluid is a Newtonian incompressible one and both are assumed static. Since the homogeneous matrix is not saturated, gradients of concentration are present. The porous medium and the fluid (a liquid) will be regarded as continuous constituents of a mixture that will have also a third constituent, an inert gas, assumed with zero mass density and thermal conductivity. The problem is described by a set of two partial differential equations which represent the energy balances for the fluid and the solid constituents. Isovalues for these two constituents are plotted, considering representative time instants and selected values for the energy equations coefficients and for the saturation.     

Effects of variable fluid properties on natural convection of MHD fluid from a heated vertical wavy surface

In the present paper, the influence of temperature-dependent fluid properties, density, viscosity and thermal conductivity on MHD natural convection flow from a heated vertical wavy surface is studied. It is assumed that, the fluid density and the thermal conductivity vary as exponential and linear functions of temperature, respectively. However, the fluid viscosity is assumed to vary as a reciprocal of a linear function of temperature. The model analysis used here is more relevant to liquid flow. Using the appropriate variables, the wavy surface are transformed into a flat one. The transformed boundary layer equations are solved numerically, using implicit-Chebyshev pseudospectral method, for several sets of values of the physical parameters, namely, the temperature dependent fluid properties parameters, the magnetic parameter, the amplitude-wavelength ratio parameter, and the Prandtl number. The numerical values obtained for the velocity, temperature, shearing stress, and the Nusselt number are presented through graphs and tables for several sets of values of the parameters. The effects of the physical parameters on the flow and heat transfer characteristics are discussed. The results were compared with numerical solutions of previous works. The present results are found to be in good agreement.     

Unsteady two-fluid flow and heat transfer in a horizontal channel

The problem of unsteady oscillatory flow and heat transfer of two viscous immiscible fluids through a horizontal channel with isothermal permeable walls has been considered. The partial differential equations governing the flow and heat transfer are solved analytically using two-term harmonic and non-harmonic functions in both fluid regions of the channel. Effects of physical parameters such as viscosity ratio, conductivity ratio, Prandtl number and frequency parameter on the velocity and temperature fields are shown graphically. It is observed that the velocity and temperature decrease as the viscosity ratio increases, while they increase with increases in frequency parameter. The effect of increasing the thermal conductivity ratio also suppresses the temperature in both fluid regions. The effect of periodic frequency on the flow is depicted in tabular form. It is predicted that both the velocity and temperature profiles decrease as the periodic frequency increases.     

Level‐set based numerical simulation of a migrating and dissolving liquid drop in a cylindrical cavity

In the present paper the dissolution of a binary liquid drop having a miscibility gap and migrating due to thermo‐solutal capillary convection in a cylindrical cavity is studied numerically. The interest in studying this problem is twofold. From a side, in the absence of gravity, capillary migration is one of the main physical mechanisms to set into motion dispersed liquid phases and from the other side, phase equilibria of multi‐component liquid systems, ubiquitous in applications, often exhibit a miscibility gap. The drop capillary migration is due to an imposed temperature gradient between the cavity top and bottom walls. The drop dissolution is due to the fact that initial composition and volume values, and thermal boundary conditions are only compatible with a final single phase equilibrium state. In order to study the drop migration along the cavity and the coupling with dissolution, a previously developed planar two‐dimensional code is extended to treat axis‐symmetric geometries. The code is based on a finite volume formulation. A level‐set technique is used for describing the dynamics of the interface separating the different phases and for mollifying the interface discontinuities between them. The level‐set related tools of redistancing and off‐interface extension are used to enhance code resolution in the critical interface region. Migration speeds and volume variations are determined for different drop radii. Copyright © 2004 John Wiley & Sons, Ltd.     

Mixed convection in an inclined lid-driven cavity with non-uniform heating on both sidewalls

The present work reports a numerical simulation of mixed convection in an inclined square cavity. The vertical sidewalls are assumed to have a nonuniform temperature distribution. The finite volume method is used to solve dimensionless governing equations. Simulations are performed for different Richardson numbers, amplitude ratios, phase deviations, and cavity inclination angles. The results are presented graphically. The mean heat transfer significantly increases in the buoyancy-dominated mode on increasing cavity inclination angle if both walls have identical heating and cooling zones.     

Conjugate natural convection in a square porous cavity

 Steady-state conjugate natural convection in a square cavity filled with a porous medium is studied numerically in this paper. The enclosure consists of two horizontal conductive walls of finite thickness and two vertical walls at different uniform temperatures. The focus is on the role of solid-fluid conductivity ratio, k, on the flow and heat transfer characteristics and the average Nusselt number, , over the vertical hot and cold walls of the cavity for a limited set of particular parameters. It was shown that the interface temperature, θ_w, along the top of the solid wall decreases with the increase in the wall conductivity k. Also, the values of decreases with the increase of the values of the parameter k. Comparison with known results from the open literature when the wall thickness of the horizontal solid walls is neglected (non-conjugate problem) is excellent. Received on 4 April 2000