Search results

The purpose of this paper is to develop a combined method for three-dimensional incompressible flow and heat transfer by the spectral collocation method (SCM) and the artificial compressibility method (ACM), and further to study the performance of the combined method SCM-ACM for three-dimensional incompressible flow and heat transfer.

The partial differentials in space are discretized by the SCM with Chebyshev polynomial and Chebyshev–Gauss–Lobbatto collocation points. The unsteady artificial compressibility equations are solved to obtain the steady results by the ACM. Three-dimensional exact solutions with trigonometric function form and exponential function form are constructed to test the accuracy of the combined method.

The SCM-ACM is developed successfully for three-dimensional incompressible flow and heat transfer with high accuracy that the minimum value of variance can reach. The accuracy increases exponentially along with time marching steps. The accuracy is also improved exponentially with the increasing of nodes before stable accuracy is achieved, while it keeps stably with the increasing of the time step. The central processing unit time increases exponentially with the increasing of nodes and decreasing of the time step.

It is difficult for the implementation of the implicit scheme by the developed SCM-ACM. The SCM-ACM can be used for solving unsteady impressible fluid flow and heat transfer.

The SCM-ACM is applied for two classic cases of lid-driven cavity flow and natural convection in cubic cavities. The present results show good agreement with the published results with much fewer nodes.

The combined method SCM-ACM is developed, firstly, for solving three-dimensional incompressible fluid flow and heat transfer by the SCM and ACM. The performance of SCM-ACM is investigated. This combined method provides a new choice for solving three-dimensional fluid flow and heat transfer with high accuracy.

This paper aims to present briefly a unified fractional step method for fluid dynamics, incompressible solid mechanics and heat transfer calculations. The proposed method is demonstrated by solving compressible and incompressible flows, solid mechanics and conjugate heat transfer problems.

The finite element method is used for the spatial discretization of the equations. The fluid dynamics algorithm used is often referred to as the characteristic‐based split scheme.

The proposed method can be employed as a unified approach to fluid dynamics, heat transfer and solid mechanics problems.

The idea of using a unified approach to fluid dynamics and incompressible solid mechanics problems is proposed. The proposed approach will be valuable in complicated engineering problems such as fluid‐structure interaction and problems involving conjugate heat transfer and thermal stresses.

The purpose of this paper is to consider the unsteady natural convection flow of a viscous incompressible fluid, which is induced by differential heating on the solid vertical boundary of an open‐ended rectangular cavity with the two horizontal surfaces which are permeable and maintained at the temperature of ambient fluid. Attention is focused on how the flow and heat transfer is affected by variations of the buoyancy force, as well as by the permeability of the surfaces.

An upwind finite‐difference method in conjunction with a successive over‐relaxation iteration technique is used to solve the governing boundary layer equations. To do this, the first and second derivatives were approximated by central differences and were used in the vorticity, energy and Poisson equations. To preserve the conservative property, the finite‐difference forms of the vorticity and energy equations were written in conservative form for the convective terms.

Local rate of heat transfer from the heated surface increases owing to an increase in the value of Ra. In the region near the bottom surface, the heat transfer from the left vertical surface decreases, but that increases in the region near the upper surface. Due to blowing of fluid through the permeable surfaces, the rate of heat transfer is higher than the situation where fluid is being withdrawn. This difference was found to be higher in the case of larger value of Ra.

The analysis is valid for unsteady, two‐dimensional natural convection flow of a viscous fluid filled in an open‐ended rectangular enclosure. An extension to three‐dimensional flow case is left for future work.

The method is very useful to analyze solar receiver systems, fire research, electronic cooling, brake housing of an aircraft and many environmental geothermal processes.

The results of this study may be of some interest to engineers interested in heat transfer in ventilated rooms or enclosures.

The conception of a heat function, just like the stream function used in a laminar two dimensional incompressible flow field visualization, has been introduced to visualize the convective heat transfer or the flow of energy around a sphere when the sphere is either being cooled or heated by a stream of fluid flowing around it. The heat function is developed in a spherical polar coordinate and is used to generate the heat lines around the sphere. The heat lines essentially show the magnitude and direction of energy transfer around the sphere with and without the existence of a finite radial velocity at the surface. The steady state hydrodynamic field around the sphere is numerically obtained up to a maximum Reynolds number of 100 and the corresponding thermal field has been obtained by solving the steady state energy equation. The field properties thus obtained are utilized to form the heat function, which becomes an effective tool for visualization of convective heat transfer.

In part I of this study, a three‐dimensional finite difference iterative solver capable of handling the coupled Navier‐Stokes and energy equations for incompressible viscous flows was described and validated with two‐ and three‐dimensional benchmarks. Part II describes the results of the computational study of two distinct complex geometries: 1) two‐dimensional and three‐dimensional natural convection in cavity whose surface is cooled while two internal blocks are heated; 2) two‐dimensional and three‐dimensional natural convection in the region defined by two interconnected cavities of different sizes which are differentially heated. All computations have been performed for a Prandtl number of 1.0, and different values of the Rayleigh number ranging between 10³ and 10⁶ depending on the problem. In the first problem, three‐dimensional effects in the top region of the cavity trap fluid in vortices near the top of the heated blocks adversely affecting heat transfer in the region while enhancing it in the region between the two heated blocks. In the second problem, the sudden expansion of fluid as it leaves the top cavity and enters the bottom one generates three‐dimensional wakes in the bottom cavity that enhance the convective heat transfer across the system walls near them. These studies tend to suggest that three‐dimensional effects play a very important role in the enhancement of convective heat transfer in complex geometries, especially at higher Rayleigh numbers.

A computational fluid dynamics (CFD) analysis has been carried out on the aerodynamic and thermal behavior of an incompressible Newtonian fluid having a constant property and flowing turbulently through a two-dimensional horizontal high-performance heat transfer channel with a rectangular cross section. The top surface of the channel was kept at a constant temperature, while it was made sure to maintain the adiabatic condition of the bottom surface. Two obstacles, with different shapes, i.e. flat rectangular and V-shaped, were inserted into the channel; they were fixed to the top and bottom surfaces of the channel in a periodically staggered manner to force vortices to improve the mixing and consequently the heat transfer. The first fin-type obstacle is placed on the heated top channel surface, and the second baffle-type one is placed on the insulated bottom surface. Five different obstacle situations were considered in this study, which are referred as cases FF (flat fin and flat baffle), FVD (flat fin and V-downstream baffle), FVU (flat fin and V-upstream baffle), VVD (V-downstream fin and V-downstream baffle) and VVU (V-Upstream fin and V-upstream baffle).

The flow model is governed by Reynolds-averaged Navier–Stokes equations with the k-epsilon turbulence model and the energy equation. These governing equations are discretized by the finite volume method, in two dimensions, using the commercial CFD software FLUENT software with the Semi Implicit Method for Pressure Linked Equations (SIMPLE) algorithm for handling the pressure-velocity coupling. Air is the test fluid with the flow rate in terms of Reynolds numbers ranging from 12,000 to 32,000.

Important deformations and large recirculation regions were observed in the flow field. A vortex causes a rotary motion inside the flow field, which enhances the mixing by bringing the packets of fluid from the near-wall region of the channel to the bulk and the other way around. The largest value of the axial variations of the Nusselt number and skin friction coefficient is found in the region facing the baffle, while the smallest value is in the region near the fin, for all cases. The thermal enhancement factor (TEF) was also introduced and discussed to assess the performance of the channel for various obstacle situations. It is found that the TEF values are 1.273-1.368, 1.377-1.573, 1.444-1.833, 1.398-1.565 and 1.348-1.592 for FF, FVD, FVU, VVD and VVU respectively, depending on the Re values. In all cases, the TEF was found to be much larger than unity; its maximum value was around 1.833 for FVU at the highest Reynolds number. Therefore, the FVU may be considered as the best geometrical configuration when using the obstacles to improve the heat transfer efficiency inside the channel.

This study can be a real application in the field of shell-and-tube heat exchangers and flat plate solar air collectors.

The purpose of this paper is to carry out a numerical study on the dynamic and thermal behavior of a fluid with a constant property and flowing turbulently through a two-dimensional horizontal rectangular channel. The upper surface was put in a constant temperature condition, while the lower one was thermally insulated. Two transverse, solid-type obstacles, having different shapes, i.e. flat rectangular and V-shaped, were inserted into the channel and fixed to the top and bottom walls of the channel, in a periodically staggered manner to force vortices to improve the mixing, and consequently the heat transfer. The flat rectangular obstacle was put in the first position and was placed on the hot top wall of the channel. However, the second V-shaped obstacle was placed on the insulated bottom wall, at an attack angle of 45°; its position was varied to find the optimum configuration for optimal heat transfer.

The fluid is considered Newtonian, incompressible with constant properties. The Reynolds averaged Navier–Stokes equations, along with the standard k-epsilon turbulence model and the energy equation, are used to control the channel flow model. The finite volume method is used to integrate all the equations in two-dimensions; the commercial CFD software FLUENT along with the SIMPLE-algorithm is used for pressure-velocity coupling. Various values of the Reynolds number and obstacle spacing were selected to perform the numerical runs, using air as the working medium.

The channel containing the flat fin and the 45° V-shaped baffle with a large Reynolds number gave higher heat transfer and friction loss than the one with a smaller Reynolds number. Also, short separation distances between obstacles provided higher values of the ratios Nu/Nu₀ and f/f₀ and a larger thermal enhancement factor (TEF) than do larger distances.

This is an original work, as it uses a novel method for the improvement of heat transfer in completely new flow geometry.

The purpose of this paper is to present a novel numerical method to solve incompressible flows with natural and mixed convections using pseudo‐compressibility formulation.

The present method is derived using the framework of Harten Lax and van Leer with contact (HLLC) method of Toro, Spruce and Spears, that was originally developed for compressible gas dynamics equations. This work generalizes the algorithm described in the previous paper to the case where heat transfer is involved. Here, the solution of the Riemann problem is approximated by a three‐wave system.

A few test cases involving incompressible laminar flows inside 2D square cavity for various Rayleigh and Reynolds numbers are considered for validating the present method. The computed results from the present method are found to be quite promising.

Although pseudo‐compressibility formulation has been found to have superior performance and has the potential to have numerical treatments similar to compressible flow equations, only two numerical methods have been applied so far; namely Jameson method and Roes flux difference splitting method. A new sophisticated numerical method, following the framework of HLLC method, is derived and implemented for solving pseudo‐compressibility‐based incompressible flow equations with heat transfer.

The aim of this paper is to investigate the convective heat transfer in swirl tubes, which are obtained by roto‐translating a circular section eccentric with respect to the rotation axis. The geometry is numerically investigated with the aim of evaluating the convective heat transfer enhancement effect due to the secondary flow induced by the centrifugal force.

The governing equations, i.e. continuity, momentum and energy equations, are integrated numerically within Comsol Multiphysics^® environment, under the assumption of incompressible Newtonian and constant properties fluid and of periodically fully developed laminar flow for what concerns both the hydrodynamic and the thermal problem under the uniform wall heat flux thermal boundary condition.

The heat transfer performance of the geometry is discussed in relation to the flow pattern. In particular, the numerical results show that two different stable flow regimes may exist, according to the ratio of the Reynolds number to the dimensionless helix pitch. The Nusselt number augmentation becomes significant for high Prandtl number fluids when a critical Re/P* value, corresponding to the onset of the centrifugal forces induced secondary flow, is reached.

The geometry here investigated represents an interesting solution to enhance the convective heat transfer in situations in which the flow, although disturbed, persists in the laminar regime. This type of enhanced tubes shows then interesting heat transfer performances (which becomes particularly significant for high Prandtl number values) by thus suggesting convenient applications also for highly viscous fluids which are often treated under the laminar flow regime.

The paper aims to consider heat transfer in incompressible flow in a rotating flat microchannel with allowance for boundary slip conditions of the first and second order. The novelty of the paper encompasses analytical and numerical solutions of the problem, with the latter based on the lattice Boltzmann method (LBM). The analytical solution of the problem includes relations for the velocity and temperature profiles and for the Nusselt number depending on the rotation rate of the microchannel and slip velocity. It was demonstrated that the velocity profiles at high rotation rates transform from parabolic to M-shaped with a minimum at the channel axis. The temperature profiles tend to become uniform (i.e. almost constant). An increase in the channel rotation rate contributes to the increase in the Nusselt number. An increase in the Prandtl number causes a similar effect. The trend caused by the effect of the second-order slip boundary conditions depends on the closure hypothesis. It is shown that heat transfer in a flat microchannel can be successfully modeled using the LBM methodology, which takes into account the second-order boundary conditions.

The paper is based on the comparisons of an analytical solution and a numerical solution, which employs the lattice Boltzmann method. Both mathematical approaches used the first-order and second-order slip boundary conditions. The results obtained using both methods agree well with each other.

The analytical solution of the problem includes relations for the velocity and temperature profiles and for the Nusselt number depending on the rotation rate of the microchannel and slip velocity. It was demonstrated that the velocity profiles at high rotation rates transform from parabolic to M-shaped with a minimum at the channel axis. The temperature profiles tend to become uniform (i.e. almost constant). The increase in the channel rotation rate contributes to the increase in the Nusselt number. An increase in the Prandtl number causes the similar effect. The trend caused by the effect of the second-order slip boundary conditions depends on the closure hypothesis. It is shown that heat transfer in a flat microchannel can be successfully modeled using the LBM methodology, which considers the second-order boundary conditions.

The novelty of the paper encompasses analytical and numerical solutions of the problem, whereas the latter are based on the LBM.