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.     

Convective heat transfer and temperature stratification in a sphere completely filled with a liquid,with a given heat flux

A number of articles have been devoted to the theoretical and experimental investigation of natural convection in spherical vessels completely filled with a liquid [1–6]. Analytical solutions are known, obtained by the expansion of the sought function in series in powers of the Rayleigh number (see, for example, [1]), valid for very small values of this number. A numerical solution of the nonlinear Boussinesq equations can be used to obtain solutions with larger Rayleigh numbers, but the existing data for spherical regions [2, 3] embrace a relatively narrow range of Rayleigh numbers. The experimental data with a given heat flux, published in [4–6], were obtained with relatively large Rayleigh numbers (Ra*=10⁹?10¹¹) and Prandtl numbers (P= 3?1500). Data on the characteristics of convection in spherical vessels are still not very numerous and, in a number of cases, contradictory. This relates, in particular, to the boundaries of unsteady-state conditions. The present article, continuing [7–9], expounds a method and gives the results of a calculation of convection in a sphere with a thinwalled shell, in a range of Rayleigh and Fourier numbers embracing the principal conditions of unsteady-state laminar convection with a given heat flux.     

Entropy generation due to natural convection in a partially heated cavity by local RBF-DQ method

The Local Radial Basis Function-Differential Quadrature (RBF-DQ) method is applied to two-dimensional incompressible Navier-Stokes equations in primitive form. Numerical results of heatlines and entropy generation due to heat transfer and fluid friction have been obtained for laminar natural convection. The variations of the entropy generation for different Rayleigh numbers are also investigated. Comparison between the present results and previous works demonstrated excellent agreements which verify the accuracy and flexibility of the method in simulating the fluid mechanics and heat transfer problems.     

Comparison of low Mach number models for natural convection problems

 We investigate in this paper two numerical methods for solving low Mach number compressible flows and their application to single-phase natural convection flow problems. The first method is based on an asymptotic model of the Navier–Stokes equations valid for small Mach numbers, whereas the second is based on the full compressible Navier–Stokes equations with particular care given to the discretization at low Mach numbers. These models are more general than the Boussinesq incompressible flow model, in the sense that they are valid even for cases in which the fluid is subjected to large temperature differences, that is when the compressibility of the fluid manifests itself through low Mach number effects. Numerical solutions are computed for a series of test problems with fixed Rayleigh number and increasing temperature differences, as well as for varying Rayleigh number for a given temperature difference. Numerical difficulties associated with low Mach number effects are discussed, as well as the accuracy of the approximations. Received on 17 January 2000     

Robust semidirect finite difference methods for solving the Navier–Stokes and energy equations

Semidirect solution techniques can be an effective alternative to the more conventional iterative approaches used in many finite difference methods. This paper summarizes several semidirect techniques which generally have not been applied to the Navier–Stokes and energy equations in finite difference form. The methods presented use both successive substitution and Jacobian-based updates as well as two variations of Broyden's full matrix update. A hybrid method is also presented, as is a norm-reducing search technique that can be used to enhance the convergence characteristics of any semidirect approach. These methods have been compared with the well known iterative methods SIMPLE and SIMPLER. The comparison was performed on the natural convection and driven cavity problems. The semidirect methods proved to be reliably convergent without the need for a priori specification of variable under-relaxation factors, which was necessary with the iterative methods. Natural convection and driven cavity solutions have been readily obtained with the proposed methods for Rayleigh and Reynolds numbers up to 10⁹ and 10⁶ respectively. Of the semidirect techniques, the hybrid approach was the most robust. From an arbitrary zero initial guess this method was able to obtain a solution to the natural convection problem for Rayleigh numbers three orders of magnitude larger than was possible with the Newton-Raphson update. The computational effort required by the semidirect methods is comparable to that required by the iterative methods; however, the memory requirements can be significantly greater.     

Stability of Thermosolutal Natural Convection in Superposed Fluid and Porous Layers

This article deals with the onset of thermosolutal natural convection in horizontal superposed fluid and porous layers. A linear stability analysis is performed using the one-domain approach. As in the thermal convection case, the results show a bimodal nature of the marginal stability curves where each mode corresponds to a different convective instability. At small wave numbers, the convective flow occurs in the whole cavity (“porous mode”) while perturbations of large wave numbers lead to a convective flow mainly confined in the fluid layer (“fluid mode”). Furthermore, it is shown that the onset of thermosolutal natural convection is characterized by a multi-cellular flow in the fluid region for negative thermal Rayleigh numbers. For positive thermal Rayleigh numbers, the convective flow takes place both in the fluid and porous regions. The influence of the depth ratio and thermal diffusivity ratio is also investigated for a wide range of the thermal Rayleigh numbers.     

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.     

Entropy generation and natural convection in square cavities with wavy walls

The aim of the present work is to study the entropy generation in the natural convection process in square cavities with hot wavy walls through numerical simulations for different undulations and Rayleigh numbers, while keeping the Prandtl number constant. The results show that the hot wall geometry affects notably the heat transfer rate in the cavity. It has been found in the present numerical study that the mean Nusselt number in the case of heat transfer in a cavity with wavy walls is lower, as compared to heat transfer in a cavity without undulations. Based on the obtained dimensionless velocity and temperature values, the distributions of the local entropy generation due to heat transfer and fluid friction, the local Bejan number, and the local entropy generation are determined and plotted for different undulations and Rayleigh numbers. The study is performed for Rayleigh numbers 10³ < Ra < 10⁵, irreversibility coefficients 10^?4 < φ < 10^?2, and Prandtl numbers Pr = 0.71. The total entropy generation is found to increase with increasing undulation number.     

Parameter-free numerical method for modeling thermal convection in square cavities in a wide range of Rayleigh numbers

Some numerical results for the two- and three-dimensional de Vahl Davis benchmark are presented. This benchmark describes thermal convection in a square (cubic) cavity with vertical heated walls in a wide range of Rayleigh numbers (10⁴ to 10¹⁴), which covers both laminar and highly turbulent f lows. Turbulent f lows are usually described using a turbulence model with parameters that depend on the Rayleigh number and require adjustment. An alternative is Direct Numerical Simulation (DNS) methods, but they demand extremely large computational grids. Recently, there has been an increasing interest in DNS methods with an incomplete resolution, which, in some cases, are able to provide acceptable results without resolving Kolmogorov scales. On the basis of this approach, the so-called parameter-free computational techniques have been developed. These methods cover a wide range of Rayleigh numbers and allow computing various integral properties of heat transport on relatively coarse computational grids. In this paper, a new numerical method based on the CABARET scheme is proposed for solving the Navier–Stokes equations in the Boussinesq approximation. This technique does not involve a turbulence model or any tuning parameters and has a second-order approximation scheme in time and space on uniform and nonuniform grids with a minimal computational stencil. Testing the technique on the de Vahl Davis benchmark and a sequence of refined grids shows that the method yields integral heat f luxes with a high degree of accuracy for both laminar and highly turbulent f lows. For Rayleigh numbers up to 10¹⁴, a several percent accuracy is achieved on an extremely coarse grid consisting of 20 × 20 cells refined toward the boundary. No definite or comprehensive explanation of this computational phenomenon has been given. Cautious optimism is expressed regarding the perspectives of using the new method for thermal convection computations at low Prandtl numbers typical of liquid metals.     

Natural convection heat transfer along vertical rectangular ducts

Experimental investigations have been reported on steady state natural convection from the outer surface of vertical rectangular and square ducts in air. Seven ducts have been used; three of them have a rectangular cross section and the rest have square cross section. The ducts are heated using internal constant heat flux heating elements. The temperatures along the vertical surface and the peripheral directions of the duct wall are measured. Axial (perimeter averaged) heat transfer coefficients along the side of each duct are obtained for laminar and transition to turbulent regimes of natural convection heat transfer. Axial (perimeter averaged) Nusselt numbers are evaluated and correlated using the modified Rayleigh numbers for laminar and transition regime using the vertical axial distance as a characteristic length. Critical values of the modified Rayleigh numbers are obtained for transition to turbulent. Furthermore, total overall averaged Nusselt numbers are correlated with the modified Rayleigh numbers and the area ratio for the laminar regimes. The local axial (perimeter averaged) heat transfer coefficients are observed to decrease in the laminar region and increase in the transition region. Laminar regimes are obtained at the lower half of the ducts and its chance to appear decreases as the heat flux increases.     

Numerical convergence studies of the mixed finite element method for natural convection flow in a fluid-saturated porous medium

This paper describes a numerical approximation scheme for the natural convection (NC) flow in a fluid-saturated porous medium. Our formulation of the problem is based on the mixed finite element method (FEM). Using the so-called consistent splitting scheme, a second-order accuracy in time and in space is verified by the numerical calculation. The resulting flow patterns and heat transfer for different Rayleigh numbers, Darcy numbers and porosities are numerically studied by the proposed scheme.     

Convective stability analysis of a micropolar fluid layer by variational method

This paper studies Rayleigh-Bénard convection of micropolar fluid layer heated from below with realistic boundary conditions. A specific approach for stability analysis of a convective problem based on variational principle is applied to characterize the Rayleigh number for quite general nature of bounding surfaces. The analysis consists of replacing the set of field equations by a variational principle and the expressions for Rayleigh number are then obtained by using trial function satisfying the essential boundary conditions. Further, the values of the Rayleigh number for particular cases of large and small values of the microrotation coefficient have been obtained. The effects of wave number and micropolar parameter on the Rayleigh numbers for onset of stationary instability for each possible combination of the bounding surfaces are discussed and illustrated graphically. The present analysis establishes that the nature of bounding surfaces combination and microrotation have significant effect on the onset of convection.     

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⁷.     

Effect of adiabatic wall on the natural convection heat transfer from a wavy surface created by attached horizontal cylinders

The effect of a vertical adiabatic wall on the natural convection heat transfer from vertical array of attached cylinders, which can be considered as wavy surface, was investigated experimentally and numerically. The experiments were carried out using Mach-Zehnder interferometer and the commercial FLUENT code was used for numerical study. This paper focuses on the effect of wall-wavy surface spacing and Rayleigh number variation on the local and average free convection heat transfer coefficients from the each cylinder and the wavy surface. Rayleigh number ranges from 2400 to 10,000 and from 300,000 to 1,250,000 based on cylinder diameter and wavy surface height respectively. The local and average Nusselt numbers were determined for the different Rayleigh numbers, and the ratio of wall- wavy surface spacing to cylinder diameter 0.75, 1, 1.5, 2, 3, 4, 5, and ∞. Results are indicated with a single correlation which gives the average Nusselt number as a function of the ratio of the wall-wavy surface spacing to cylinder diameter and the Rayleigh numbers. There is an optimum distance between the wall and wavy surface in which the Nusselt number attain its maximum value. This optimum distance depends on the Rayleigh number.     

Laminar natural convection in inclined enclosures bounded by a solid wall

Laminar natural convection has been studied in enclosures bounded by a solid wall with its outer boundary at constant temperature while the opposing side has a constant heat flux. Two-dimensional equations of conservation of mass, momentum and energy, with the Boussinesq approximation are solved using a finite difference method. The numerical procedure adopted is based on the SIMPLER algorithm. Various parameters were: Rayleigh number (from 10³ to 10⁶), dimensionless conductivity of bounding wall (from 1 to 10) and dimensionless wall width (from 0.15 to 0.5), aspect ratio (from 0.5 to 1) and the inclination angle (from 30^° to 180^°). The results are reduced in terms of the normalized Nusselt number as a function of the Rayleigh number, and other dimensionless parameters. The isotherms and streamlines are produced for various Rayleigh numbers and geometrical conditions. It is found that the heat transfer is an increasing function of the Rayleigh number, wall to fluid conductivity ratio, enclosure aspect ratio and a decreasing function of the wall thickness. It passes from a maximum for the inclination angle of about 80^°.     

Interferometric study of buoyancy-driven convection in a differentially heated circular fluid layer

The present study is concerned with buoyancy-driven convection experiments in a circular horizontal differentially heated layer of air. The radius-to-height ratio of 14, and Rayleigh numbers of 5,861 and 12,124 have been considered. A Mach–Zehnder interferometer has been used to visualize the convection patterns in the fluid layer. The fluid layer has been imaged at view angles of 0, 45 and 90°. Results obtained show that fringe patterns appropriate for a cavity square in plan are seen in the fluid layer during the early stages of the experiments. After the passage of the initial transients, steady fringes have been observed in the fluid layer for a Rayleigh number of 5,861. At Ra=12,124, a dominant pattern was detectable combined with mild unsteadiness. The steady thermal field at Ra=5,861 displayed symmetry with respect to the viewing angle. A stronger three dimensionality was seen at the higher Rayleigh number. The average steady state Nusselt numbers were found to vary with view angle from 1.91 to 2.04 at Ra=5,861 and 2.28 to 2.43 at Ra = 12,124. The cavity-averaged Nusselt numbers are in good agreement with the available correlations.     

非定常流函数涡量方程的一种数值解法的研究

对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解．以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数（Ra）条件下的数值研究．结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟．     

Thermal stratification studies in a side heated water pool for advanced heavy water reactor applications

Strong heat source at the isolation condenser wall of an Advanced Heavy Water Reactor, results in natural convection in gravity driven water pool, which leads to a thermally stratified pool. Governing equations simulating fluid flow and heat distribution are solved numerically by a general purpose Computational Fluid Dynamics solver developed at Indian Institute of Technology, Kanpur. Incompressible finite volume method with non-staggered grid arrangement is used in this exercise. This algorithm is fully implicit and semi-coupled. Turbulent natural convection in a boundary layer for high Rayleigh numbers is analyzed by the Lam–Bremhorst k − ε turbulence model. Analysis of unsteady laminar natural convection in a side-heated water cavity is also done for different values of Rayleigh number. Results show a warm fluid layer floating on the top of gradually colder layer (along the vertical direction) that indicates the presence of thermal stratification phenomenon. This fact necessitates additional safety features in such a system so that the detrimental effect such as stratification is minimized.     

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.     

LES-based filter-matrix lattice Boltzmann model for simulating turbulent natural convection in a square cavity

In this paper, a novel thermal filter-matrix lattice Boltzmann model based on large eddy simulation (LES) is proposed for simulating turbulent natural convection. In this study, the Vreman subgrid-scale eddy-viscosity model is introduced into the present framework of LES to accurately predict the flow in near-wall region. Two dimensional numerical simulations of natural convection in a square cavity were performed at high Rayleigh number varying from 10⁷ to 10¹⁰ with a fixed Prandtl number of Pr = 0.71. The influences of the higher-order terms upon the present results at high Rayleigh numbers are examined, taking Ra = 10⁷ and 10⁸ as the example, revealing that the proper minimization of the higher-order terms can improve numerical accuracy of present model for high Rayleigh convective flow. For the turbulent convective flow, the time-averaged quantities in the median lines are presented and compared against those available results from previous studies. The general structure of turbulent boundary layers is well predicted. All numerical results exhibit good agreement with the benchmark solutions available in the previous literatures.