A high-accuracy temporal–spatial pseudospectral method for time-periodic unsteady fluid flow and heat transfer problems

A temporal–spatial pseudospectral (TSP) method is proposed for the high-accuracy solutions of time-periodic unsteady fluid flow and heat transfer problems. In this method, both the spatial and temporal derivative terms in the governing equations are computed by pseudospectral method. The spatial derivatives are computed through Chebyshev and Lagrange polynomials while the time derivatives are computed by Fourier series. The TSP method is capable of directly finding out the periodic state solutions without the necessity to resolve the initial transient state solutions, hence holds high computational efficiency and high numerical accuracy properties for the time-periodic problems. This method is validated by three 2D benchmark problems: the time-periodic incompressible flow with exact solutions; the natural convection in enclosure with time-periodic temperature on one sidewall, and on both sidewalls. The TSP results fit well the exact solutions or the benchmark solutions and the TSP accuracy is much higher than the time marching spatial pseudospectral accuracy. Some time-dependent fluid flow and heat transfer characteristic parameters are analysed. The proposed TSP method could be further extended to more complex time-periodic unsteady fluid flow and heat transfer problems where high-accuracy results are required.     

Conjugate Natural Convection in a Porous Enclosure with Non-Uniform Heat Generation

Effects of a conductive wall on natural convection in a square porous enclosure having internal heating at a rate proportional to a power of temperature difference is studied numerically in this article. The horizontal heating is considered, where the vertical walls heated isothermally at different temperatures while the horizontal walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and finite difference method is applied to solve the dimensionless governing equations. The governing parameters considered are the Rayleigh number (0 ???Ra ???1000), the internal heating and the local exponent parameters (0 ????? ???5), (1 ????? ???3), the wall to porous thermal conductivity ratio (0.44 ???Kr ???9.9) and the ratio of wall thickness to its width (0.02 ???D ???0.5). The results are presented to show the effect of these parameters on the fluid flow and heat transfer characteristics. It is found a strong internal heating can generate significant maximum fluid temperature more than the conductive solid wall. Increasing value thermal conductivity ratio and/or decreasing the thickness of solid wall can increase the maximum fluid temperature. It is also found that at very low Rayleigh number, the heat transfer across the porous enclosure remain stable for any values of the thermal conductivity ratio.     

Numerical simulation of natural convective heat transfer and fluid flow around a heated cylinder inside an enclosure

 A collocated, non-orthogonal grid based finite volume technique has been applied for investigating the two dimensional natural convective flow and heat transfer around a heated cylinder kept in a square enclosure. The effects of different enclosure wall thermal boundary conditions, fluid Prandtl number and the ratio between enclosure and cylinder dimensions (aspect ratio) upon the flow and thermal features, have been systematically studied. It is observed that the patterns of recirculatory flow and thermal stratification in the fluid are significantly modified, if any of these parameters is varied. The overall heat transfer rates are also affected due to the changes in the flow and temperature patterns. The study presents useful observations regarding the variation of local Nusselt number along each wall, for the different cases considered. Received on 2 August 2000 / Published online: 29 November 2001     

Effect of variable thermal conductivity models on the combined convection heat transfer in a square enclosure filled with a water–alumina nanofluid

In this numerical study, the effects of variable thermal conductivity models on the combined convection heat transfer in a two-dimensional lid-driven square enclosure are investigated. The fluid in the square enclosure is a water-based nanofluid containing alumina nanoparticles. The top and bottom horizontal walls are insulated, while the vertical walls are kept at different constant temperatures. Five different thermal conductivity models are used to evaluate the effects of various parameters, such as the nanofluid bulk temperature, nanoparticle size, nanoparticle volume fraction, Brownian motion, interfacial layer thickness, etc. The governing stream–vorticity equations are solved by using a second-order central finite difference scheme coupled with the conservation of mass and energy. It is found that higher heat transfer is predicted when the effects of the nanoparticle size and bulk temperature of the nanofluid are taken into account.     

Natural convection in inclined partitioned enclosures

The problem of steady, laminar, natural convective flow of a viscous fluid in an inclined enclosure with partitions is considered. Transverse gradient of temperature is applied on the two opposing regular walls of the inclined enclosure while the other walls are maintained adiabatic. The problem is formulated in terms of the vorticity-stream function procedure. A numerical solution based on the finite volume method is obtained. Representative results illustrating the effects of the enclosure inclination angle and the degree of irregularity on the contour maps of the streamlines and temperature are reported and discussed. In addition, results for the average Nusselt number at the heated wall of the enclosure and the difference of extreme stream-function values are presented and discussed for various Rayleigh numbers, inclination angles and dimensionless partition heights.     

Natural convection flow under a magnetic field in an inclined square enclosure differentialy heated on adjacent walls

Steady, laminar, natural-convection flow in the presence of a magnetic field in an inclined square enclosure differentially heated along the bottom and left vertical walls while the other walls are kept isothermal was considered. The governing equations were solved numerically for the stream function, vorticity and temperature ratio using the differential quadrature method for various Grashof and Hartmann numbers, inclination angle of the enclosure and direction of the magnetic field. The orientation of the enclosure changes the temperature gradient inside and has a significant effect on the flow pattern. Magnetic field suppresses the convective flow and its direction also influences the flow pattern, causing the appearance of inner loops and multiple eddies. The surface heat flux along the bottom wall is slightly increased by clockwise inclination and reduced by half by the counterclockwise inclination. The surface heat flux along the upper portion of the left side wall is reversed by the rise of warmer fluids due to the convection currents for no inclination and clockwise inclination of the enclosure.     

Natural convection flow of water near its density maximum in a rectangular enclosure having isothermal walls with heat generation

We consider unsteady laminar natural convection flow of water subject to density inversion in a rectangular cavity formed by isothermal vertical walls with internal heat generation. The top and bottom horizontal walls are considered to be adiabatic, whereas the temperature of the left vertical wall is assumed to be greater than that of the right vertical wall. The equations are non-dimensionalized and are solved numerically by an upwind finite difference method together with a successive over-relaxation (SOR) technique. The effects of both heat generation and variations in the aspect ratio on the streamlines, isotherms and the rate of heat transfer from the walls of the enclosure are presented. Investigations are performed for water taking Prandtl number to be Pr=11.58 and the Rayleigh number to be Ra=10⁵.     

Conjugate Phase Change Heat Transfer in an Inclined Compound Cavity Partially Filled with a Porous Medium: A Deformed Mesh Approach

In this paper, the melting process of a PCM inside an inclined compound enclosure partially filled with a porous medium is theoretically addressed using a novel deformed mesh method. The sub-domain area of the compound enclosure is made of a porous layer and clear region. The right wall of the enclosure is adjacent to the clear region and is subject to a constant temperature of T_c. The left wall, which is connected to the porous layer, is thick and thermally conductive. The thick wall is partially subject to the hot temperature of T_h. The remaining borders of the enclosure are well insulated. The governing equations for flow and heat transfer, including the phase change effects and conjugate heat transfer at the thick wall, are introduced and transformed into a non-dimensional form. A deformed grid method is utilized to track the phase change front in the solid and liquid regions. The melting front movement is controlled by the Stefan condition. The finite element method, along with Arbitrary Eulerian–Lagrangian (ALE) moving grid technique, is employed to solve the non-dimensional governing equations. The modeling approach and the accuracy of the utilized numerical approach are verified by comparison of the results with several experimental and numerical studies, available in the literature. The effect of conjugate wall thickness, inclination angle, and the porous layer thickness on the phase change heat transfer of PCM is investigated. The outcomes show that the rates of melting and heat transfer are enhanced as the thickness of the porous layer increases. The melting rate is the highest when the inclination angle of the enclosure is 45°. An increase in the wall thickness improves the melting rate.

     

MHD natural convection in porous media-filled enclosures

In this work, the magnetohydrodynamics （MHD） natural convection heat transfer problem inside a porous medium filled with inclined rectangular enclosures is investigated numerically. The boundary conditions selected on the enclosure are two adiabatic and two isothermal walls. The governing equations, continuity, and Forchheimer extension of the Darcy law and energy are transformed into dimensionless forms by using a set of suitable variables, and then solved by using a finite difference scheme. The governing parameters are the magnetic influence number, the Darcy Rayleigh number, the inclination angle, and the aspect ratio of the enclosure. It is found that the magnetic influence number and the inclination angle have pronounced effects on the fluid flow and heat transfer in porous media-filled enclosures.     

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.     

Combined natural convection heat and mass transfer in an enclosure having finite thickness walls

Unsteady three-dimensional conjugate heat and mass transfer in an enclosure having finite thickness heat-conducting walls has been analyzed numerically. The governing unsteady, three-dimensional flow, energy and contaminant transport equations for the gas cavity and unsteady heat conduction equation for solid walls, written in dimensionless terms of the vector potential functions, the vorticity vector, the temperature and the concentration, have been solved using an iterative implicit finite-difference method. Main attention was paid to the effects of the Rayleigh number, buoyancy ratio and the dimensionless time on the flow structure and heat and mass transfer regimes. It should be noted that the dominant cause of the oscillations in the dimensionless time dependences of the average Nusselt number on the heat source surface and the average Sherwood number on the contaminant source surface at Ra>5?10⁵ is the mutual influence of the analyzed object geometry and the thermo-diffusivity impact on the flow. The change in the buoyancy ratio can lead to the essential modifications of the flow, temperature and concentration fields owing to the significant influence of the concentration gradient.     

An experimental study of the effect of wall temperature nonuniformity on natural convection in an enclosure heated from the side

An experimental study of natural convection in a parallelepipedal enclosure induced by a single vertical wall is described. The upper half of this wall was warm and the lower half cold. The other enclosure walls were insulated. The temperature and flow measurements were performed in the high Rayleigh number regime (10¹⁰<Ra<5×10¹⁰) by using water as the working fluid. The Rayleigh number was based on the enclosure height and the temperature difference between the warm and the cold part of the driving wall. The flow field featured two flat cells, one filled with warm fluid along the top horizontal wall, and the other filled with cold fluid along the bottom horizontal wall. Each of these cells was surrounded by an additional cell as tall as half the enclosure height. The above flow structure prohibited extensive thermal contact between warm and cold fluid, thus limiting the role of convection on the heat transfer process in the cavity. The findings of this study differ significantly from the findings of previous studies based on the ‘classical’ enclosure model possessing two isothermal vertical walls, the one warm and the other cold, and support the view that the use of ‘more realistic‘ temperature boundary conditions in enclosure natural convection needs careful examination.     

Experimental and numerical study on laminar natural convection in a cavity heated from bottom due to an inclined fin

Natural convection heat transfer in an inclined fin attached square enclosure is studied both experimentally and numerically. Bottom wall of enclosure has higher temperature than that of top wall while vertical walls are adiabatic. Inclined fin has also adiabatic boundary conditions. Numerical solutions have been done by writing a computer code in Fortran platform and results are compared with Fluent commercial code and experimental method. Governing parameters are Rayleigh numbers (8.10⁵ ≤ Ra ≤ 4 × 10⁶) and inclination angle (30° ≤ and ≤ 120°). The temperature measurements are done by using thermocouples distributed uniformly at the wall of the enclosure. Remarkably good agreement is obtained between the predicted results and experimental data. A correlation is also developed including all effective parameters on heat transfer and fluid flow. It was observed that heat transfer can be controlled by attaching an inclined fin onto wall.     

Effect of Conduction in Bottom Wall on Darcy–Bénard Convection in a Porous Enclosure

Conjugate natural convection-conduction heat transfer in a square porous enclosure with a finite-wall thickness is studied numerically in this article. The bottom wall is heated and the upper wall is cooled while the verticals walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and the COMSOL Multiphysics software is applied to solve the dimensionless governing equations. The governing parameters considered are the Rayleigh number (100 ≤ Ra ≤ 1000), the wall to porous thermal conductivity ratio (0.44 ≤ K _r ≤ 9.90) and the ratio of wall thickness to its height (0.02 ≤ D ≤ 0.4). The results are presented to show the effect of these parameters on the heat transfer and fluid flow characteristics. It is found that the number of contrarotative cells and the strength circulation of each cell can be controlled by the thickness of the bottom wall, the thermal conductivity ratio and the Rayleigh number. It is also observed that increasing either the Rayleigh number or the thermal conductivity ratio or both, and decreasing the thickness of the bounded wall can increase the average Nusselt number for the porous enclosure.     

正弦波温度边界下多孔介质方腔内非热平衡对流传热数值模拟

采用局部非热平衡模型,在方腔左侧壁面温度正弦波变化、右侧壁面温度均一的边界条件下,通过SIM-PLER算法数值研究了固体骨架发热多孔介质方腔内的稳态非达西自然对流,主要探讨了不同正弦波波动参数N及方腔的高宽比M/L对方腔内自然对流与传热的影响规律。计算结果表明:正弦波温度边界使得方腔内的流场出现了复杂的变化,流体及固体区域左侧壁面附近出现了周期性的正负变化的温度场分布,左侧壁面局部Nusselt数出现了周期性的震荡现象;存在一个最佳温度波动参数N=1,此时多孔介质方腔内的整体散热量达到最大值;增加方腔高宽比会显著地削弱方腔内的自然对流传热过程,小高宽比也会在一定的程度上削弱多孔介质方腔内的对流传热。     

Mixed convection flow in a lid-driven inclined square enclosure filled with a nanofluid

This work is focused on the numerical modeling of steady laminar mixed convection flow in a lid-driven inclined square enclosure filled with water–Al₂O₃ nanofluid. The left and right walls of the enclosure are kept insulated while the bottom and top walls are maintained at constant temperatures with the top surface being the hot wall and moving at a constant speed. The developed equations are given in terms of the stream function–vorticity formulation and are non-dimensionalized and then solved numerically subject to appropriate boundary conditions by a second-order accurate finite-volume method. Comparisons with previously published work are performed and found to be in good agreement. A parametric study is conducted and a set of graphical results is presented and discussed to illustrate the effects of the presence of nanoparticles and enclosure inclination angle on the flow and heat transfer characteristics. It is found that significant heat transfer enhancement can be obtained due to the presence of nanoparticles and that this is accentuated by inclination of the enclosure at moderate and large Richardson numbers.     

Fluid flow due to combined convection in lid-driven enclosure having a circular body

The present work is aimed to study mixed convection heat transfer characteristics for a lid-driven air flow within a square enclosure having a circular body. Flows are driven by the left lid, which slides in its own plane constant velocity. This wall is isothermal and it moves up or down in y direction while the other walls remain stationary. The horizontal walls are adiabatic. The cavity is differentially heated and the left wall is maintained at a higher temperature than the right wall. Three different temperature boundary conditions were applied for the inner cylinder as adiabatic, isothermal or conductive. The computation is carried out for wide ranges of Richardson numbers, diameter of inner cylinder and center and location of the inner cylinder. It was found that the most effective parameter on flow field and temperature distribution is the orientation of the moving lid. The circular body can be a control parameter for heat and fluid flow. An interesting obtained result that the thermal conductivity becomes insignificant for small values of diameter of the circular body.     

Natural convection in attic-shaped spaces subject to sudden and ramp heating boundary conditions

In this study, a discussion of the fluid dynamics in the attic space is reported, focusing on its transient response to sudden and linear changes of temperature along the two inclined walls. The transient behaviour of an attic space is relevant to our daily life. The instantaneous and non-instantaneous (ramp) heating boundary condition is applied on the sloping walls of the attic space. A theoretical understanding of the transient behaviour of the flow in the enclosure is performed through scaling analysis. A proper identification of the timescales, the velocity and the thickness relevant to the flow that develops inside the cavity makes it possible to predict theoretically the basic flow features that will survive once the thermal flow in the enclosure reaches a steady state. A time scale for the heating-up of the whole cavity together with the heat transfer scales through the inclined walls has also been obtained through scaling analysis. All scales are verified by the numerical simulations.     

Unsteady thermal entrance heat transfer in laminar pipe flows with step change in ambient temperature

A fully implicit upwind finite difference numerical scheme has been proposed to investigate the characteristics of thermal entrance heat transfer in laminar pipe flows subject to a step change in ambient temperature. In order to demonstrate the results more clearly, a modified Nusselt number is introduced. The unsteady axial variations of modified Nusselt number, bulk fluid temperature, and wall temperature and the transient temperature profiles at certain axial locations are presented graphically for various outside heat transfer coefficients. The effects of the outside heat transfer coefficient on the heat transport processes in the flow are examined in detail. The results can be comprehensively explained by the interaction between the upstream convective heat transfer and the diffusion heat transfer in the radial direction. Steady state is reached when the axial convection balances the radial diffusion.     

Thermal Non-equilibrium Natural Convection in a Square Enclosure with Heat-Generating Porous Layer on Inner Walls

Differentially heated enclosure with heat-generating porous layer on inner walls is studied computationally for non-Darcy flow and thermal non-equilibrium models. In this study, this problem is investigated for different internal and external Rayleigh numbers, Darcy numbers, porosity-scaled thermal conductivity ratio, solid-/fluid-scaled heat transfer coefficient and dimensionless thickness of the porous layer. The results indicate that the dimensionless thickness of the porous layer has an important effect on the heat transfer in the enclosure. It was found that the thermal non-equilibrium model is needed for small values of the porosity-scaled thermal conductivity ratio and the solid-/fluid-scaled heat transfer coefficient. It is shown that the convection of heat due to internal heat generation is increased in the enclosure when the ratio of internal Rayleigh number to external Rayleigh number is larger.