Natural convection of a viscoplastic fluid in vertical circular channels

In contrast to a Newtonian fluid, a viscoplastic fluid can be in a state of mechanical equilibrium when heat is supplied from the side. Therefore, natural convection in a viscoplastic fluid heated from the side occurs only when the determining parameters exceed certain threshold values. The threshold conditions for the onset of convection for a flat vertical layer have been investigated several times [1–4]. The present paper is an investigation of the conditions of occurrence of plane-parallel natural convection of a viscoplastic fluid in regions with cylindrical symmetry: in a vertical annular layer and in a vertical circular tube.     

An adaptive remeshing technique for laminar natural convection problems

 An adaptive remeshing procedure based on interpolation error has been developed for solving laminar natural convection problems. A simple relation between the Rayleigh number and the equilibration constant is introduced to increase the efficiency of the method. Unstructured meshes have been regenerated based on an initial solution and used in the calculations. This avoids the expensive mesh sensitivity study. Two typical natural convection problems are solved to demonstrate the present technique. Excellent performance has been observed at moderate and high Rayleigh numbers although present method gives good results over the whole range of Rayleigh numbers considered. Received on 27 October 2000 / Published online: 29 November 2001     

Higher order stability effects in a natural convection boundary layer over a vertical heated wall

 The stability of a laminar boundary layer flow under natural convection on a vertical isothermally heated wall is studied analytically. The analysis is performed by using two different two-dimensional linear models: (1) The non-parallel flow model in which the steady mean flow as well as the disturbance amplitude functions can change in the streamwise direction; (2) The parallel flow model in which the effects of the mean flow and disturbance changes in the streamwise direction are neglected. The linear non-parallel stability analysis is based on the so-called parabolised stability equations (PSEs) which have been successfully applied to the stability analysis of forced convection boundary layers. In this study the PSE equations are applied to natural convection boundary layers in order to show the difference between parallel and non-parallel stability analysis. A second part of this study deals with the effects of variable properties, which are always present in natural convection flows. They are analysed by an extended version of the Orr–Sommerfeld equation (EOSE). Received on 31 May 2000     

A transient finite element analysis of natural convection around a horizontal hot cylinder

This paper presents an advanced method for a 2-dimensional analysis of transient natural convection by finite element method. The present method, based on stream function—vorticity formulation, could get rid of numerical errors and constraint of perpendicular mesh subdivision, since we excluded a finite difference approximation of vorticity on no-slip boundaries. A considerable effect of upwind weighting function was examined. The method was successfully applied to a problem of natural convection around a horizontal hot cylinder.     

Numerical investigation of transient natural convection heat transfer of freezing water in a square cavity

In this study, a transient heat transfer process of freezing water inside a two-dimensional square cavity has been investigated numerically. Water was used as a phase-change medium, and the numerical model has been created with control volume approach by using C++ programming language. To be able to accelerate the numerical calculations, CUT (Consistent-Update-Technique) algorithm has been implemented in the numerical code. Span-wise variations of the vertical component of the velocity have been represented in comparison with the experimental measurements from the literature at various vertical positions to examine the accuracy of the numerical scheme. The influence of natural convection has been considered by comparing the conduction and convection dominated solidification under same boundary conditions. Comparative results have been obtained regarding time-wise variations of the cold wall temperature and the dimensionless effectiveness. Moreover, the streamlines and isotherms have been represented to understand the differences between the conduction and convection driven phase change processes.Results indicate that natural convection becomes remarkable and has different forms at the initial periods of the phase change process. Increasing the effect of natural convection in the cavity increases the cooling rate of water. Near the density inversion temperature of water (4°C), temperature variations fluctuate and counter currents observed in the domain.     

Numerical simulation of natural convection of latent heat phase-change-material microcapsulate slurry packed in a horizontal rectangular enclosure heated from below and cooled from above

A two-dimensional numerical simulation of natural convection in a rectangular enclosure heated from below and cooled from above has been conducted with non-Newtonian phase-change-material (PCM) microcapsulate slurry with latent heat capacities. The formulation of the mathematical model in dimensionless co-ordinates and discretization of the governing equations have been done using the finite volume method. Both natural convection and heat transfer characteristics are discussed about natural convection with PCM microcapsulate slurry, which exhibits the pseudoplastic non-Newtonian fluid behavior and a peak value in the specific heat capacity with latent heat. The viscosity of the present PCM microcapsulate slurry is assumed to follow the Ostwald-de Waele power law fluid model with the power-law index n and the consistency coefficient K. The effects of phase-change material, the mass concentration, and the aspect ratio Ar on the natural convection heat transfer are described, respectively. By comparing with the results of microcapsule slurry without phase change, the enhancement in heat transfer is found in microcapsule slurry with phase change during the phase change temperature range. Numerical simulations are performed in the following parametric ranges: the width–height aspect ratio of the enclosure Ar from 2 to 20, the mass concentrations C _m of the slurry from 10 to 40%, power law index n of the slurry from 0.89 to 1.0 and Rayleigh numbers Ra ranges from 10³ to 10⁷.     

Hydrodynamics in weak force fields onset of steady thermocapillary convection in a spherical fluid layer under zero-g conditions

In the study of cellular convection in an infinite plane fluid layer with a free surface, both the Archimedes and thermocapillary forces [1–3] have been cited as reasons for the onset of convection. This has also been confirmed experimentally [4], When mass forces are absent or negligibly small it is natural to pose the question of the onset of pure thermocapillary convection or convection caused only by the surface tension gradients (see [2–3]). In the present paper, this problem is examined for a spherical fluid layer under zero-g conditions.     

The Effect of Heavy Oil Viscosity Reduction by Solvent Dissolution on Natural Convection in the Boundary Layer of VAPEX

We have studied the effect of viscosity on natural convection in the boundary layer of the vapor extraction (VAPEX) process. VAPEX is a heavy oil recovery method that uses solvents to reduce oil viscosity, and is a potential process in reservoirs where thermal recovery methods cannot be applied. Natural convection may happen in VAPEX if the solvents that are used to decrease oil viscosity increase the density of the oil. This can especially occur with $\text{ CO }_{2}$ CO 2 -based solvents. Reduction of the oil viscosity due to solvent dissolution can have a large impact on the onset of convection by decreasing the critical Rayleigh number. When the viscosity reduction is significant, the critical Rayleigh can decrease up to two orders of magnitude. The transverse Peclet number is also a crucial parameter in determining the critical Rayleigh and onset of convection. Our analysis shows that the longitudinal Peclet does not have a significant impact on the natural convention in VAPEX. When oil viscosity reduction is included in the analysis of boundary layer instability in VAPEX, natural convection may occur in high-permeable reservoirs (where Rayleigh number is high) leading to a greater oil production rate compared with current models where the effect of boundary layer instability has been ignored.     

Étude expérimentale et numérique de la convection thermomagnétique dans un fluide non électroconducteur

We present an experimental and numerical study of convection driven by a central heat source in a vertical and cylindrical cavity with a higth aspect ration, filled by air and submitted to the action of non homogeneous intence magnetic field. We show how the magnetization force applying on oxygene is able to intensify, reduce, stop or even reverse the natural convection fluid flow.     

A Multiple-Relaxation-Time Lattice-Boltzmann Analysis for Double-Diffusive Natural Convection in a Cavity with Heating and Diffusing Plate Inside Filled with a Porous Medium

Transport in Porous Media - In the present study, a multiple-relaxation-time lattice-Boltzmann method is considered to investigate double-diffusive natural convection in a cavity with heating and...     

Magnetic and Gravitational Convection of Air in a Porous Cubic Enclosure with a Coil Inclined Around the Y Axis

The natural convection heat transfer of air in a porous media can be controlled by gradient magnetic field. Thermomagnetic convection of air in a porous cubic enclosure with an electric coil inclined around the $Y$ axis was numerically investigated. The Biot–Savart law was used to calculate the magnetic field. The governing equations in primitive variables were discretized by the finite-volume method and solved by the SIMPLE algorithm. The flow and temperature fields for the air natural convection were presented and the mean Nusselt number on the hot wall was calculated and compared. The results show that both the magnetic force and coil inclination have significant effect on the flow field and heat transfer in a porous cubic enclosure, the natural convection heat transfer of air can be enhanced or controlled by applying gradient magnetic field.     

Natural convection heat transfer on surfaces of copper micro-wires

The natural convection heat transfer characteristics and mechanism for copper micro-wires in water and air were investigated experimentally and numerically. The wires with diameters of 39.9, 65.8 and 119.1 μm were placed horizontally in water inside of a sealed tube and in air of a large room, respectively. Using Joule heating, the heat transfer coefficients and Nusselt numbers of natural convection for micro-wires in ultra pure water and air were obtained. A three dimensional incompressible numerical model was used to investigate the natural convection, and the prediction with this model was in reasonable accordance with the experimental results. With the decrease of micro-wire diameter, the heat transfer coefficient of natural convection on the surface of micro-wire becomes larger, while the Nu number of natural convection decreases in water and air. Besides, the change rate of Nu number in water decreases apparently with the increase of heat flux and the decrease of wire diameter, which is larger than that in air. The thickness of boundary layer on the wall of micro-wire becomes thinner with the decrease of diameter in both water and air, but the ratio of boundary layer thickness in water to the diameter increases. However, there is almost no change of this ratio for natural convection in air. As a result, the proportion of conduction in total heat transfer of natural convection in water increases, while the convective heat transfer decreases. The velocity distribution, temperature field and the boundary layer in the natural convection were compared with those of tube with conventional dimension. It was found that the boundary layer around the micro-wire is an oval-shaped film on the surface, which was different from that around the conventional tube. This apparently reduces the convection strength in the natural convection, thus the heat transfer presents a conduction characteristic.     

Experimental study of free convection in an elliptical annular enclosure in blunt and slender orientations

In this paper, free convection heat transfer between two elliptical cylinders having different elliptical ratios in blunt and slender situations was studied experimentally. Three pairs of elliptical cylinders having the same radius ratio of 2.0, the same surface area and different elliptical ratios of 0.662, 0.866 and 0.968 were cut using the computed numerically controlled wire-cut machining. The tests were carried out by keeping a constant heat flux on the inner cylinder while cooling the outer one to be isothermal. The effects of vertical eccentricity, lateral eccentricity, angle of attack of the inner cylinder on natural convection for both blunt and slender situations of each pair were investigated. Empirical correlation was deduced within an acceptable uncertainty for the experimental results. Compatible and satisfactory to the conscience agreement was found in the comparison among the results of present and previous works. In the vision of the comparison, it was found that; the vertical eccentricity can enhance free convection by about 15% than the concentric case, the horizontal eccentricity can enhance natural convection by about 10% than concentric case and slender situation offers about 40% enhancement in free convection than the blunt situation for the same elliptical ratio.     

Experimental Study on Oscillatory Natural Convection in a Hele-Shaw Cell due to Unstably Heated Side

The oscillatory motion of natural convection in a porous medium has been investigated experimentally using a Hele-Shaw cell technique. The cell has been heated on the lower half and cooled on the upper half along the same vertical sidewall. Flows have been visualized using the pH indicator method. Photographs of natural convection patterns as well as average Nusselt number data have been presented for different Rayleigh numbers. Oscillatory motion of natural convection has been observed for large enough Rayleigh numbers and the critical Rayleigh number has been estimated to be between 120 and 450. Scaling analysis has been conducted to understand the heat transfer and the oscillating mechanism. According to the scaling analysis, it has been found that the average Nusselt number is proportional to the square root of the Rayleigh number, and that the oscillation frequency is proportional to the Rayleigh number. Obtained experimental data support the scaling analysis.     

Experimental and numerical study on unsteady natural convection heat transfer in helically coiled tube heat exchangers

Both of experimental and numerical investigations were performed to understand unsteady natural convection from outer surface of helical coils. Four helical coils with two different curvature ratios were used. Each coil was mounted in the shell both vertically and horizontally. The cold water was entered the coil and the hot water in the shell was cooling by unsteady natural convection. A CFD code was developed to simulate natural convection heat transfer. Equations of tube and shell are solved simultaneously. Statistical analyses have been done on data points of temperature and natural convection Nusselt number. It was revealed that shell-side fluid temperature and the Nusselt number of the outer surface of coils are functions of in-tube fluid mass flow rate, specific heat of fluids and geometrical parameters including length, inner diameter of the tube and the volume of the shell, and time.     

Linear Stability Analysis of Double-Diffusive Convection in Porous Media,with Application to Geological Storage of CO<Subscript>2</Subscript>

Onset of double-diffusive buoyancy-driven flow resulted from vertical temperature and concentration gradients in a horizontal layer of a saturated and homogenous porous medium is investigated using amplification factor theory. After injection of CO₂ into a deep saline aquifer, the density of the brine saturated with CO₂ increases slightly. This increase in density induces natural convection. The effect of geothermal gradient is also considered in this work as a second incentive for convection and the double-diffusion convection was studied. Linear stability analysis is used to predict the inception of instabilities and initial wavelength of the convective instabilities. The analysis presented is applied to acid gas injection (as an analogue for CO₂ storage) into saline aquifers in the Alberta basin. It is found that the geothermal gradient does not have significant effect on the onset of convection for these aquifers. It is shown that the geothermal effects on the onset of natural convection are negligible as compared to the solutal effects induced by dissolution and diffusion of CO₂ in deep saline aquifers. Therefore, the linear stability analysis and the long-term numerical simulation of CO₂ sequestration into such saline aquifers may be assumed to be isothermal in terms of natural convection occurrence.     

ISPH modeling of natural convection heat transfer with an analytical kernel renormalization factor

The objective of this study is to extend the attention of the incompressible smoothed particle hydrodynamics method (ISPH) in the heat transfer field. The ISPH method for the natural convection heat transfer under the Boussinesq approximation in various environments: pure-fluid, nanofluid, and non-Darcy porous medium is introduced. We adopted the improved analytical method for calculating the kernel renormalization factor and its gradient based on a quintic kernel function for the wall boundary treatment in the ISPH method. The proposed method requires no dummy particle layer to meet the impermeability condition and makes the heat flux over the wall boundary easy to implement. We performed four different numerical simulations of natural convection in cavities with increasing complexity in modeling and implementation: the natural convection in a square cavity with constant differentially heated wall temperature, natural convection with the heat flux from the bottom wall for a wide range of Rayleigh numbers, natural convection in a non-Darcy porous cavity fully filled with nanofluid in different flow regimes, and natural convection in a partially layered porous cavity. The results showed excellent agreement with results from literatures and the in-house P1–P1 finite element method code.     

Modeling,numerical simulation and experimental verification of the unsteady cooling of a solid body in quiescent ambient air

The scope of the present article is two-fold. Firstly, to conduct an experiment to provide the temperature-time history of the cooling of a hot ball bearing in quiescent ambient air. Secondly, to predict the temporal variation of the bearing under the hypothesis of natural convection, radiation or natural convection coexists with radiation for a non-vanishing total hemispherical emissivity of the surface of the bearing. Numerical solutions of the three governing nonlinear lumped heat equations were carried out with a Runge-Kutta-Fehlberg (RKF45) algorithm accounting for automatic step size control. The experimental data was obtained with chrome steel ball bearings of diameter 0.953 cm (7/16 in) heated in an electric oven to a pre-set temperature. The heated bearing was exposed later to ambient air at atmospheric temperature and pressure.     

Heat transfer from extended surfaces subject to variable heat transfer coefficient

The present article investigates the effect of locally variable heat transfer coefficient on the performance of extended surfaces (fins) subject to natural convection. Fins of different profiles have been investigated. The fin profiles presently considered are namely; straight and pin fin with rectangular (constant diameter), convex parabolic, triangular (conical) and concave parabolic profiles and radial fins with constant profile with different radius ratios. The local heat transfer coefficient was considered as function of the local temperature and has been obtained using the available correlations of natural convection for each pertinent extended surface considered. The performance of the fin has been expressed in terms of the fin efficiency. Comparisons between the present results for all fins considered and the results obtained for the corresponding fins subject to constant heat transfer coefficient along the fin are presented. Comparisons, i.e. showed an excellent agreement with the experimental results available in the literature. Results show that there is a considerable deviation between the fin efficiency calculated based on constant heat transfer coefficient and that calculated based on variable heat transfer coefficient and this deviation increases with the dimensionless parameter m.     

侧加热腔内的自然对流

开展侧加热腔内自然对流的研究具有重大的环境及工业应用背景. 总结侧加热腔内水平温差驱动的自然对流的最新研究进展, 并概述相应的流动性质、动力机制和传热特性以及对不同无量纲控制参数的依赖也有重要的科学价值. 已取得的研究结果显示突然侧加热的腔内自然对流的发展可包括初始阶段、过渡阶段和定常或准定常阶段. 不同发展阶段的流动依赖于瑞利数、普朗特数及腔体的高宽比, 且定常或准定常阶段的流态可以是定常层流流动、非定常周期性流动或者湍流流动. 此外, 回顾了对流流动失稳机制的研究成果以及湍流自然对流方面的新进展. 最后, 展望了侧加热腔内的自然对流研究的前景.