Search results

A fast false implicit transient scheme FITS is developed to predict the two‐dimensional steady‐state solutions of buoyancy‐assisted laminar internal flows. This new scheme uses the control volume based on power law technique in conjugation with the alternating direction implicit (ADI) and the successive grid refinement (SGR) procedures to solve the transient vorticity and energy transport equations. The ADI procedure allows the power law, which gives an excellent approximation to the exact 1‐D solution, to be applied locally in one‐dimensional sense for each sweep in the co‐ordinates’ directions. This in turn increased the solution accuracy and hence permits the use of a larger time increment. As a result a remarkable increase in the convergence rate to steady‐state is achieved. The final solution is obtained by successively refining the grid as the solution advances in time. The efficiency of FITS is verified by comparing the present predictions with three steady‐state benchmark solutions: natural convection of a heat generating fluid in a rectangular enclosure, natural convection inside a cavity with two isothermal walls, and a vertical buoyancy‐assisted laminar backward‐facing step flow.

The purpose of this paper is to consider a two‐dimensional, steady, mixed convection flow in an enclosure filled with a fluid‐saturated uniform porous medium. The left wall is moving down and the right wall is moving up and are maintained at cold and hot constant temperatures, respectively. The top and the bottom walls are fixed and are thermally insulated.

The governing equations are normalized and solved numerically with appropriate boundary conditions by finite‐volume approach using third‐order accurate scheme (deferred QUICK).

The study is conducted by varying the key parameters, i.e. Richardson number (Ri=Gr/Re²), Darcy number (Da=κ /H²) and Grashof number (Gr=gβ H³ Δ T/ν ²) and fixing Prandtl number at (Pr=0.71). A parametric study is conducted and a set of streamline and isotherm plots are presented. The average Nusselt number reaches a value of 1 asymptotically when the Ri is gradually increased for Gr up to 10³. The asymptotic value is 1.5 for Gr=10⁴. A heat transfer correlation is also presented.

The study of the mixed convection problem with lid‐driven flows in enclosures will be useful in cooling of electronic devices, lubrication technologies, chemical processing equipment, etc.

This study aims to perform a numerical study of air convection in a rectangular enclosure with two isothermal blocks and oscillating bottom wall temperature under laminar flow conditions. The geometry of the enclosure contains two isothermal blocks placed equidistant along the streamwise direction. The top wall is assumed to be cold (low temperature). The bottom wall temperature is either kept as constant or sinusoidally varied with time. The vertical walls are considered as adiabatic. The flow is diagonally upwards and assisted by the buoyancy force. The inlet is positioned at the bottom of the left wall, and the outlet is placed at the top of the right wall. The parameters considered in this paper are Rayleigh number (104-106), Prantdl number (0.71), amplitude of temperature oscillation (0-0.5) and the period (0.2). The effects of these parameters on heat transfer and fluid flow inside the open cavity are studied. The periodic results of fluid flow are illustrated with streamlines and the heat transfer is represented by isotherms and time-averaged Nusselt number. By virtue of increasing buoyancy, the heat transfer accelerates with an increase in the Rayleigh number. Also, the heat transfer is intensive with an increase in the bottom wall temperature.

The momentum and energy equations are solved simultaneously. The energy equation (3) is initially solved using the alternating direction implicit (ADI) method. The results of the energy equation are updated into the vorticity equation. The unsteady vorticity transport equation is also solved using the ADI method. Dimensionless time step equal to 0.01 is used for high Ra (10⁵ and 10⁶) and 0.001 is used for low Ra (10⁴). Convergence criteria of 10⁻⁵ is used during the vorticity, stream function and temperature calculations, as the sum of error should be very small.

Numerical study of air convection in a rectangular enclosure with two isothermal blocks and oscillating bottom wall temperature is performed under laminar flow condition. The effect of the isothermal blocks on the heat transfer is analyzed for different Rayleigh numbers and the following conclusions are arrived. The hydrodynamic blockage effect is subdued by the isothermal heating of square blocks. Based on the streamline diagrams, it is found that the formation of vortices is greatly influenced by the Rayleigh number when all the walls are exposed to a constant wall temperature. The influence of amplitude on the heat transfer is remarkable on the wall exposed to oscillating temperature and is subtle on the opposite static cold wall. The heat transfer increases with an increase in the Rayleigh number and temperature.

Flow is assumed to be two-dimensional and laminar subject to oscillatory boundary condition. The present investigation aims to study natural convection inside the cavity filled with air whose bottom wall is subject to time-variant temperature. The buoyancy is further intensified through two isothermal square blocks placed equidistant along the streamwise direction at mid-height.

The authors have developed a CFD solver to simulate the situation. Effect of Rayleigh number subject to oscillatory thermal boundary condition is simulated. Streamline contour and isotherm contour are presented. Local and average Nusselt numbers are presented.

Turbulent fluid flow and heat transfer characteristics are analyzed numerically for fluids flowing through a rotating periodical two‐pass square channel. The smooth walls of this two‐pass channel are subject to a constant heat flux. A two‐equation k‐ε turbulence model with modified terms for Coriolis and rotational buoyancy is employed to resolve this elliptic problem. The duct through‐flow rate and rotating speed are fixed constantly; while the wall heat flux into the fluid is varied to examine the rotating buoyancy effect on the heat transfer and fluid flow characteristics. It is disclosed that the changes in local heat transfer due to the rotational buoyancy in the radially outward flow are more significant than those in the radially inward flow. However, the channel averaged heat transfer is altered slightly due to the rotational buoyancy in the both ducts. Whenever the buoyancy effects are sufficiently strong, the flow reversal appears over the leading face of the radially outward‐flow channel, and the radial distance for initiation of flow separation decreases with increasing the buoyancy parameter. A comparison of the present numerical results with the available experimental data by taking buoyancy into consideration is also presented.

The aim of this article is to present the results of a parametric analysis of the entropy generation due to mixed convection in the entry‐developing region between two differentially heated isothermal vertical plates.

The entropy generation was estimated via a numerical solution of the mass, momentum and energy conservation equations governing the flow and heat transfer in the vertical channel between the two parallel plates. The resultant temperature and velocity profiles were used to estimate the entropy generation and other heat transfer parameters over a wide range of the operating parameters. The investigated parameters include the buoyancy parameter (Gr/Re), Eckert number (Ec), Reynolds number (Re), Prandtl number (Pr) and the ratio of the dimensionless temperature of the two plates (θ_T).

The optimum values of the buoyancy parameter (Gr/Re) optimum at which the entropy generation assumes its minimum for the problem under consideration have been obtained numerically and presented over a wide range of the other operating parameters. The effect of the other operating parameters on the entropy generation is presented and discussed as well.

The results of this investigation are limited to the geometry of vertical channel parallel plates under isothermal boundary conditions. However, the concept of minimization of entropy generation via controlling the buoyancy parameter is applicable for any other geometry under any other thermal boundary conditions.

The results presented in this paper can be used for optimum designs of heat transfer equipment based on the principle of entropy generation minimization with particular focus on the optimum design of plate and frame heat exchanger and the optimization of electronic packages and stacked packaging of laminar‐convection‐cooled printed circuits.

This paper introduces the entropy generation minimization via controlling the operating parameters and clearly identifies the optimum buoyancy parameter (Gr/Re) at which entropy generation assumes its minimum under different operating conditions.

This paper's aim is to examine flow and heat transfer through vertical channels between parallel plates, which is of prime importance in the design of cooling systems for electronic equipment such as that of finned cold plates in general, plate‐and‐frame heat exchangers, etc.

Numerical and analytical solutions are presented to investigate the heat transfer enhancement and the pressure drop reduction due to buoyancy effects (for buoyancy‐aided flow) for the developing laminar mixed convection in vertical channel between parallel plates in the vicinity of the critical values of the buoyancy parameter (Gr/Re)crt that are obtained analytically. The numerical solutions are presented for a wide range of the buoyancy parameters Gr/Re that cover both of buoyancy‐opposed and buoyancy‐aided flow situations under each of the isothermal boundary conditions under investigation.

Buoyancy parameters greater than the critical values result in building‐up the pressure downstream of the entrance such that the vertical channel might act as a thermal diffuser with possible incipient flow reversal. Locations at which the pressure gradient vanishes and the locations at which the pressure‐buildup starts have been numerically obtained and presented for all the investigated cases.

The study is limited to the laminar flow situation.

The results clearly show that for buoyancy‐aided flow, the increase of the buoyancy parameter enhances the heat transfer and reduces the pressure drop across the vertical channel. These findings are very useful for cooling channel or chimney designs.

The study is original and presents new findings, since none of the previous studies reported the conditions for which pressure buildup might take place due to mixed convection in vertical channels between parallel plates.

The purpose of this paper is to study numerically the fluid flow and heat transfer in an inclined channel provided with heated porous blocks on its lower plate.

The Brinkman‐Forchheimer extended Darcy model with the Boussinesq approximation is adopted for the flow in the porous regions. The governing equations with the appropriate boundary conditions are solved by the control volume method. The effect of some pertinent parameters such as the buoyancy force intensity, the porous blocks shape and height, the porous medium permeability and the Reynolds number are analyzed for various inclination angles ranging from −90° to +90°.

The results reveal, essentially, that the inclination angle of the channel can alter substantially the fluid flow and heat transfer mechanisms, especially at high Richardson and Darcy numbers. In this case, the maximum and minimum global Nusselt numbers are reached for α=+90° and α=−90°, respectively.

The results obtained in this work are valid for an inclined channel with porous blocks attached on the heated parts of the lower plate, whereas the upper wall is thermally insulated.

The results obtained in this worky can be used in the thermal control of electronic components. The use of porous blocks mounted on the heat sources will increase the rate of heat removal in order to maintain the electronic components at an acceptable operating temperature.

The paper provides an interesting method to improve the cooling of electronic devices by use of a porous medium.

The purpose of this paper is to provide a solution for natural convection in a cavity with a partial heater in case of volumetric heating and analysis of the entropy generation.

The control volume method based on three-dimensional (3D) vorticity-potential vector was applied to solve governing equations of natural convection in a 3D cavity with a fin for different governing parameters as external Rayleigh numbers (103=Ra_E=106), internal Rayleigh numbers 103=Ra_I=106, partition height (0.25=h=0.75) and partition location (0.25=c=0.75). A code was written by using Fortran platform.

The edge of the fin becomes important on entropy generation. The ratio of the Ra_I/Ra_E plays the important role on natural convection and entropy generation. The variation of external Rayleigh number becomes insignificant for the Ra_I/Ra_E>1.

The originality of this work is to analyze the entropy generation and natural convection in a cubical cavity with volumetrically heating.

The purpose of this paper is to propose an efficient/robust numerical algorithm for solving the two‐dimensional laminar mixed‐convection in a lid‐driven cavity using the mixed finite element (FE) technique.

A numerical algorithm was based on the so‐called consistent splitting scheme, which improved the numerical accuracy of the primitive variables. In order to obtain a stable solution, two choices of mixed FEs, the Taylor‐Hood and Crouzeix‐Raviart types, were used. Two mesh layouts were considered; uniform and non‐uniform.

To verify that the proposed scheme had a second‐order accuracy, some numerical results are presented and compared with the known solution. The answer was confirmative. Numerically accurate solutions were obtained for a fixed Prandtl number, Pr=0.71, for a range of the Reynolds number, Re from 100 to 3,000, and for a range of the Richardson number, Ri from 0.001 to 100. The results from these calculations, using the mixed FE consistent splitting scheme, agreed with the existing ones.

Further extensions of this work could include the influence of various choices of Reynolds numbers, Prandtl numbers and Richardson numbers, and the effect of aspect ratio.

The present work was the first to apply a mixed FE in association with the consistent splitting scheme to the mixed convection problem.

The purpose of this paper is to study thermal (natural) convection in nine different containers involving the same area (area= 1 sq. unit) and identical heat input at the bottom wall (isothermal/sinusoidal heating). Containers are categorized into three classes based on geometric configurations [Class 1 (square, tilted square and parallelogram), Class 2 (trapezoidal type 1, trapezoidal type 2 and triangle) and Class 3 (convex, concave and triangle with curved hypotenuse)].

The governing equations are solved by using the Galerkin finite element method for various processing fluids (Pr = 0.025 and 155) and Rayleigh numbers (10³ ≤ Ra ≤ 10⁵) involving nine different containers. Finite element-based heat flow visualization via heatlines has been adopted to study heat distribution at various sections. Average Nusselt number at the bottom wall ( Nub¯) and spatially average temperature (θ^) have also been calculated based on finite element basis functions.

Based on enhanced heating criteria (higher Nub¯ and higher θ^), the containers are preferred as follows, Class 1: square and parallelogram, Class 2: trapezoidal type 1 and trapezoidal type 2 and Class 3: convex (higher θ^) and concave (higher Nub¯).

The comparison of heat flow distributions and isotherms in nine containers gives a clear perspective for choosing appropriate containers at various process parameters (Pr and Ra). The results for current work may be useful to obtain enhancement of the thermal processing rate in various process industries.

Heatlines provide a complete understanding of heat flow path and heat distribution within nine containers. Various cold zones and thermal mixing zones have been highlighted and these zones are found to be altered with various shapes of containers. The importance of containers with curved walls for enhanced thermal processing rate is clearly established.