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.     

Numerical simulation of steady-state and transient natural convection in an isothermal open cubic cavity

In this work the numeric results, of the steady-state and transient heat transfer by natural convection in a horizontal isothermal open cubic cavity are presented. The most important assumptions in the mathematical formulation are two, the flow is laminar and the Boussinesq approximation is valid. The conservation equations in primitive variables are solved using the finite volume method and the SIMPLEC algorithm. The advective terms are approximated by the SMART scheme and the diffusive terms are approximated using the central differencing scheme. The results are obtained for a Rayleigh number range from 10⁴ to 10⁷.The numerical model predicted flow instabilities and Nusselt number oscillations for high Rayleigh numbers.     

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.     

Finite volume multigrid prediction of laminar natural convection: Bench-mark solutions

A finite volume multigrid procedure for the prediction of laminar natural convection flows is presented, enabling efficient and accurate calculations on very fine grids. The method is fully conservative and uses second-order central differencing for convection and diffusion fluxes. The calculations start on a coarse (typically 10 × 10 control volumes) grid and proceed to finer grids until the desired accuracy or maximum affordable storage is reached. The computing times increase thereby linearly with the number of control volumes. Solutions are presented for the flow in a closed cavity with side walls at different temperatures and insulated top and bottom walls. Rayleigh numbers of 10⁴, 10⁵ and 10⁶ are considered. Grids as fine as 640 × 640 control volumes are used and the results are believed to be accurate to within 0–01%. Second-order monotonic convergence to grid-independent values is observed for all predicted quantities.     

Time-dependent natural convection in a square cavity: Application of a new finite volume method

A new finite volume (FV) approach with adaptive upwind convection is used to predict the two-dimensional unsteady flow in a square cavity. The fluid is air and natural convection is induced by differentially heated vertical walls. The formulation is made in terms of the vorticity and the integral velocity (induction) law. Biquadratic interpolation formulae are used to approximate the temperature and vorticity fields over the finite volumes, to which the conservation laws are applied in integral form. Image vorticity is used to enforce the zero-penetration condition at the cavity walls. Unsteady predictions are carried sufficiently forward in time to reach a steady state. Results are presented for a Prandtl number (Pr) of 0-71 and Rayleigh numbers equal to 10³, 10⁴ and 10⁵. Both 11 × 11 and 21 × 21 meshes are used. The steady state predictions are compared with published results obtained using a finite difference (FD) scheme for the same values of Pr and Ra and the same meshes, as well as a numerical bench-mark solution. For the most part the FV predictions are closer to the bench-mark solution than are the FD predictions.     

Two numerical modelings of free convection heat transfer using nanofluids inside a square enclosure

This work describes the numerical simulation of natural convection heat transfer of Cu–water nanofluids in a square enclosure for Rayleigh numbers varying from 10³ up to 10⁵. Two different numerical approaches were used: the finite volume method and the finite element method. The nanofluids were assumed to be single-phase fluids with modified thermal properties obtained from experimental results and theoretical models. The results showed that the Nusselt number for nanofluids was basically the same as that obtained for the base fluid. Therefore, the enhancement observed in the heat transfer coefficient was significant due to the augmentation in the thermal conductivity.     

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.     

Natural convection in a rectangular cavity with wall temperature decreasing at a uniform rate

The transient natural convection in a fluid contained in a rectangular enclosure, the wall of which is maintained at a uniform temperature which changes at a steady rate, is approached by a numerical method. Numerical solutions are obtained forPr=0.73, 7.3 and 73 and a range of Rayleigh numbersRa=10² ～ 10⁸. At relatively low Rayleigh numbers the flow is characterized by the development of double cells with flow up the center and down the sidewalk However it was found that an increase of the Rayleigh number leads to the development of strong secondary circulation on the axis of symmetry of the cavity near the top wall. Thus, as the Rayleigh number is increased the secondary cells grow in size. The effects of the secondary cells on the temperature field and heat transfer coefficients are discussed. Most results are obtained for the case of a square cavity (E=2) but the influence of the aspect ratio of the cavity is also studied forE=1 and 4.     

A least-squares finite element method for time-dependent incompressible flows with thermal convection

The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l²-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 10⁶, lid-driven cavity flow at Reynolds numbers up to 10⁴ and flow over a square obstacle at Reynolds number 200, are presented to validate the method.     

Numerical investigation of transient laminar natural convection of air in a tall cavity

Transient laminar natural convection of air in a tall cavity has been studied numerically. The Navier-Stokes and Energy equations were solved by the accurate projection method (PmIII), in which the derived Poisson equation for pressure potential was solved by the approximate factorization one method (AF1). The aspect ratio of the tall cavity is 16, and the Prandtl number of air filled in the tall cavity is 0.71. To obtain the numerical results of heat transfer by natural convection of air in the tall cavity, the second order schemes for the space and time discretizations were utilized. The availability of the numerical algorithm was also assessed by considering the natural convection of air in a square cavity which is differentially heated from side walls. It was found that the overall Nusselt numbers for the Rayleigh numbers covering the range from 1000 to 100000 reveal a good agreement with measured data. When Ra takes the value in the range from 100000 to 600000, the overall Nusselt number have a relative deviation less than 18% from the experimental data. For the suddenly heating mode, the multicellular flow pattern occurs when Rayleigh number belongs to the range of Ra from 7000 to 20000. or greater than 115000. At the critical number of cats' eye instability, the cell distance is just twice of the cavity width. This is rather similar to the observed result for Bénard problem. When Ra is over 115000, a further increase of heat flux across the tall cavity causes serious cell-breaking. There are 8 cells when Ra = 600000.     

The numerical solution of the unsteady natural convection flow in a square cavity at high Rayleigh number using SADI method

The unsteady natural convection flow in a square cavity at high Rayleigh number Ra=10 ⁷ and 2×10 ⁷ has been computed using cubic spline integration. The required solutions to the two dimensional Navier-Stokes and energy equations have been obtained using two alternate numerical formulations on non-uniform grids. The main features of the transient flow have been briefly discussed. The results obtained by using the present method are in good agreement with the theoretical predictions ^[1,2].The steady state results have been compared with accurate solutions presented recently for Ra=10 ⁷.     

Natural convection and conduction heat transfer in open shallow cavities with bounding walls

A numerical study has been carried out on inclined open shallow cavities, which are formed by a wall and horizontal fins. Constant heat flux is applied on the surface of the wall inside the cavity while its other surface was kept isothermal. The wall and the fins are conductive. Conjugate heat transfer by natural convection and conduction is studied by numerically solving equations of mass, momentum and energy. Streamlines and isotherms are produced, heat and mass transfer is calculated. A parametric study is carried out using following parameters: Rayleigh number from 10⁶ to 10¹², conductivity ratio from 1 to 60, open cavity aspect ratio from 1 to 0.125, dimensionless end wall thickness from 0.05 to 0.20, horizontal walls from 0.01 to 0.15 and inclination of the end wall from 90° to 45°. It is found that the volume flow rate and Nusselt number are a decreasing function of the cavity aspect ratio, horizontal fin thickness and conductivity ratio. They are an increasing function of end wall thickness and inclination angle, except in the latter case optima exist at high Rayleigh numbers.     

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

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

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.     

Numerical study of convection of a liquid heated from below

The equilibrium of a liquid heated from below is stable only for small values of the vertical temperature gradient. With increase of the temperature gradient a critical equilibrium situation occurs, as a result of which convection develops. If the liquid fills a closed cavity, then there is a discrete sequence of critical temperature gradients (Rayleigh numbers) for which the equilibrium loses stability with respect to small characteristic disturbances. This sequence of critical gradients and motions may be found from the solution of the linear problem of equilibrium stability relative to small disturbances. If the temperature gradient exceeds the lower critical value, then (for steady-state heating conditions) there is established in the liquid a steady convective motion of a definite amplitude which depends on the magnitude of the temperature gradient. Naturally, the amplitude of the steady convective motion cannot be determined from linear stability theory; to find this amplitude we must solve the problem of convection with heating from below in the nonlinear formulation. A nonlinear study of the steady motion of a liquid in a closed cavity with heating from below was made in [1]. In that study it was shown that for Rayleigh numbers R which are less than the lower critical value R_c steady-state motions of the liquid are not possible. With R>R_c a steady convection arises, whose amplitude near the threshold is small and proportional to (R–R_c)^1/2 (the so-called soft instability)-this is in complete agreement with the results of the phenom-enological theory of Landau [2, 3].Primarily, various versions of the method of expansion in powers of the amplitude [4–8] have been used, and, consequently, the results obtained in those studies are valid only for values of R which are close to R_c, i. e., near the convection threshold.It is apparent that the study of developed convective motion far from the threshold can be carried out only numerically, with the use of digital computers. In [9, 10] the numerical methods have been successfully used for the study of developed convection in an infinite plane horizontal liquid layer.The present paper undertakes the numerical study of plane convective motions of a liquid in a closed cavity of square section. The complete nonlinear system of convection equations is solved by the method of finite differences on a digital computer for various values of the Rayleigh number, the maximal value exceeding by a factor of 40 the minimal critical value R_c. The numerical solution permits following the development of the steady motion which arises with R>R_c in the course of increase of the Rayleigh number and permits study of the oscillatory motions which occur at some value of the parameter R. The heat transfer through the cavity is studied. The corresponding linear problem on equilibrium stability is solved approximately by the Galerkin method.     

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^°.     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

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