On the Hydrodynamics and Heat Convection of an Impinging External Flow Upon a Cylinder with Transpiration and Embedded in a Porous Medium

This paper extends the existing studies of heat convection by an external flow impinging upon a flat porous insert to that on a circular cylinder inside a porous medium. The surface of the cylinder is subject to constant temperature and can include uniform or non-uniform transpiration. These cylindrical configurations are introduced in the analyses of stagnation-point flows in porous media for the first time. The equations governing steady transport of momentum and thermal energy in porous media are reduced to simpler nonlinear differential equations and subsequently solved numerically. This reveals the dimensionless velocity and temperature fields of the stagnation-point flow, as well as the Nusselt number and shear stress on the surface of the cylinder. The results show that transpiration on the surface of the cylinder and Reynolds number of the external flow dominate the fluid dynamics and heat transfer problems. In particular, non-uniform transpiration is shown to significantly affect the thermal and hydrodynamic responses of the system in the circumferential direction. However, the permeability and porosity of the porous medium are found to have relatively smaller influences.     

Flow structures and heat transfer on dimples in a staggered arrangement

Vortex structures and heat transfer enhancement mechanism of turbulent flow over a staggered array of dimples in a narrow channel have been investigated using Large Eddy Simulation (LES), Laser Doppler Velocimetry (LDV) and pressure measurements for Reynolds numbers Re_H = 6521 and Re_H = 13,042.The flow and temperature fields are calculated by LES using dynamic mixed model applied both for the velocity and temperature. Simulations have been validated with experimental data obtained for smooth and dimpled channels and empiric correlations. The flow structures determined by LES inside the dimple are chaotic and consist of small eddies with a broad range of scales where coherent structures are hardly to detect. Proper Orthogonal Decomposition (POD) method is applied on resolved LES fields of pressure and velocity to identify spatial–temporal structures hidden in the random fluctuations. For both Reynolds numbers it was found that the dimple package with a depth h to diameter D ratio of h/D = 0.26 provides the maximum thermo-hydraulic performance. The heat transfer rate could be enhanced up to 201% compared to a smooth channel.     

Analysis of Thermal Flow in a Rotating Porous U-Turn Duct Using Lattice Boltzmann Method

The lattice Boltzmann method is carried out to investigate the heat transfer enhancement in a U-turn duct which is partially filled with a porous media. The porous layer is inserted at the core of the duct and is modeled using the Brinkman–Forchheimer assumptions. In order to validate the results, first a channel flow problem without any porous layer is compared with available data. Second, the porous Couette flow and partially porous channel flow are successfully compared with the studies of other researchers. Then, fluid flow in a clear U-turn duct is studied looking carefully at the velocity, curvature and rotation effects. Finally, the effects of porous layer thickness on the rate of heat transfer and pressure drop are investigated. Parametric studies are conducted to evaluate the effects of various parameters (i.e., Reynolds number, Darcy number, rotation number), highlighting their influences on the thermo-hydrodynamics behavior of the flow. The optimum values of porous layer thickness are presented for specific flow parameters.     

Numerical Study of Heat Transfer in a Rectangular Channel with Porous Covering Obstacles

A numerical study is performed to analyze steady laminar forced convection in a channel in which discrete heat sources covered with porous material are placed on the bottom wall. Hydrodynamic and heat transfer results are reported. The flow in the porous medium is modeled using the Darcy–Brinkman–Forchheimer model. A computer program based on control volume method with appropriate averaging for diffusion coefficient is developed to solve the coupling between solid, fluid, and porous region. The effects of parameters such as Reynolds number, Prandtl number, inertia coefficient, and thermal conductivity ratio are considered. The results reveal that the porous cover with high thermal conductivity enhances the heat transfer from the solid blocks significantly and decreases the maximum temperature on the heated solid blocks. The mean Nusselt number increases with increase of Reynolds number and Prandtl number, and decrease of inertia coefficient. The pressure drop along the channel increases rapidly with the increase of Reynolds number.     

A pore scale study on turbulent combustion in porous media

This paper presents pore scale simulation of turbulent combustion of air/methane mixture in porous media to investigate the effects of multidimensionality and turbulence on the flame within the pores of porous media. In order to investigate combustion in the pores of porous medium, a simple but often used porous medium consisting of a staggered arrangement of square cylinders is considered in the present study. Results of turbulent kinetic energy, turbulent viscosity ratio, temperature, flame speed, convective heat transfer and thermal conductivity are presented and compared for laminar and turbulent simulations. It is shown that the turbulent kinetic energy increases from the inlet of burner, because of turbulence created by the solid matrix with a sudden jump or reduction at the flame front due to increase in temperature and velocity. Also, the pore scale simulation revealed that the laminarization of flow occurs after flame front in the combustion zone and turbulence effects are important mainly in the preheat zone. It is shown that turbulence enhances the diffusion processes in the preheat zone, but it is not enough to affect the maximum flame speed, temperature distribution and convective heat transfer in the porous burner. The dimensionless parameters associated with the Borghi–Peters diagram of turbulent combustion have been analyzed for the case of combustion in porous media and it is found that the combustion in the porous burner considered in the present study concerns the range of well stirred reactor very close to the laminar flame region.     

Generalized Plain Couette Flow and Heat Transfer in a Composite Channel

An analytical study of fluid flow and heat transfer in a composite channel is presented. The channel walls are maintained at different constant temperatures in such a way that the temperatures do not allow for free convection. The upper plate is considered to be moving and the lower plate is fixed. The flow is modeled using Darcy–Lapwood–Brinkman equation. The viscous and Darcy dissipation terms are included in the energy equation. By applying suitable matching and boundary conditions, an exact solution has been obtained for the velocity and temperature distributions in the two regions of the composite channel. The effects of various parameters such as the porous medium parameter, viscosity ratio, height ratio, conductivity ratio, Eckert number, and Prandtl number on the velocity and temperature fields are presented graphically and discussed.     

Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel

We analyse the convection flow of a viscous fluid through a horizontal channel enclosing a fully saturated porous medium. The Galerkin finite element analysis is used to discuss the flow and heat transfer through the porous medium using serendipity elements. The velocity, the temperature distributions and the rate of heat transfer are analysed for variations in the governing parameters. The profiles at different vertical levels are asymmetric curves, exhibiting reversal flow everywhere except on the midplane. In a given porous medium, for fixed G or N, the temperature in the fluid region at any position in fluids with a higher Prandtl number, is much higher than in fluids with a lower Prandtl number. Likewise, other parameters being fixed, lesser the permeability of the medium, lower the temperature in the flow field. Nu reduces across the flow at all axial positions, while it enhances along the axial direction of the channel. Nu reduces with decrease in the Darcy parameter D, and thus lesser the permeability of the medium, lesser the rate of heat transfer across the boundary at any axial position of the channel.     

Generalised Couette flow of two immiscible viscous fluids with heat transfer using Brinkman model

The present study deals with generalised Couette flow of two viscous, incompressible, immiscible fluids with heat transfer in presence of heat source through two straight parallel horizontal walls. The lower wall is bounded below, by a naturally permeable material of high porosity and the flow inside the porous medium is assumed to be of moderate permeability, modelled by Brinkman equation. The flow domain is divided into three zones to obtain analytical solutions of the momentum and energy equations. To link various flow regions, appropriate matching conditions have been used. The effects of permeability parameter, Reynolds number and viscous parameter on velocity field and the effects of Reynolds number, viscous parameter, permeability parameter, constant heat source and Brinkman number on temperature distribution in different zones are discussed graphically. The mass flow rate, skin-friction factor and rates of heat transfer at the upper boundary and porous interface are discussed with the help of tables.     

Conjugate heat transfer inside a porous channel

Analytical and numerical analyses have been performed for fully developed forced convection in a fluid-saturated porous medium channel bounded by two parallel plates. The channel walls are assumed to be finite in thickness. Conduction heat transfer inside the channel wall is also accounted and the full problem is treated as a conjugate heat transfer problem. The flow in the porous material is described by the Darcy–Brinkman momentum equation. The outer surfaces of the solid walls are treated as isothermal. A temperature dependent volumetric heat generation is considered inside the solid wall only. Analytical expressions for velocity, temperature, and Nusselt number are obtained after simplifying and solving the governing differential equations with reasonable approximations. Subsequent results obtained by numerical calculations show an excellent agreement with the analytical results.     

Fully Developed Mixed Convection Flow in a Vertical Channel Containing Porous and Fluid Layer with Isothermal or Isoflux Boundaries

An analysis of fully developed combined free and forced convective flow in a fluid saturated porous medium channel bounded by two vertical parallel plates is presented. The flow is modeled using Brinkman equation model. The viscous and Darcy dissipation terms are also included in the energy equation. Three types of thermal boundary conditions such as isothermal–isothermal, isoflux–isothermal, and isothermal–isoflux for the left–right walls of the channel are considered. Analytical solutions for the governing ordinary differential equations are obtained by perturbation series method. In addition, closed form expressions for the Nusselt number at both the left and right channel walls are derived. Results have been presented for a wide range of governing parameters such as porous parameter, ratio of Grashof number and Reynolds number, viscosity ratio, width ratio, and conductivity ratio on velocity, and temperature fields. It is found that the presence of porous matrix in one of the region reduces the velocity and temperature.     

Response to Comments on the Article Titled ‘The Starting Flow in Ducts Filled with a Darcy–Brinkman Medium’

The fully developed flow and constant flux heat transfer in super-elliptic ducts filled with a porous (Darcy–Brinkman) medium are studied. Super-elliptic ducts resemble rectangular ducts with rounded corners. An efficient Ritz method is used to determine the velocity and temperature fields. Extensive tables for friction factor–Reynolds number product and Nusselt number are given.     

Developing and fully developed turbulent flow in ribbed channels

Wall-mounted roughness features, such as ribs, are often placed along the walls of a channel to increase the convective surface area and to augment heat transfer and mixing by increasing turbulence. Depending on the relative roughness size and orientation, the ribs also have varying degrees of increased pressure losses. Designs that use ribs to promote heat transfer encompass the full range of having only a few streamwise ribs, which do not allow fully developed flow conditions, to multiple streamwise ribs, which do allow the flow to become fully developed. The majority of previous studies have focused on perturbing the geometry of the rib with little attention to the spatially and temporally varying flow characteristics and their dependence on the Reynolds number. A staggered rib-roughened channel study was performed using time-resolved digital particle image velocimetry (TRDPIV). Both the developing (entry region) and a fully developed region were interrogated for three Reynolds numbers of 2,500, 10,000, and 20,000. The results indicate that the flow was more sensitive to Reynolds number at the inlet than within the fully developed region. Despite having a similar mean-averaged flowfield structure over the full Reynolds number range investigated, the population and distribution of coherent structures and turbulent dissipation within the fully developed region were also found to be Reynolds number dependent. Exploring the time-accurate flow characteristics revealed that in addition to vortices shed from the rib shear layer, the region of the rib wake was governed by a periodic process of bursting of the wake vortices resulting in the intermittent ejection of the inter-rib recirculation region into the core flow. This periodic process was the driving mechanism resulting in mixing and heat transfer augmentation. A quadrant-splitting burst analysis was also performed to determine the characteristic frequency and duration of inter-rib bursting as well as the wake shedding frequency, both of which were determined to be Reynolds number dependent.     

Steady and Unsteady Heat Transfer in a Channel Partially Filled with Porous Media Under Thermal Non-Equilibrium Condition

Steady and pulsatile flow and heat transfer in a channel lined with two porous layers subject to constant wall heat flux under local thermal non-equilibrium (LTNE) condition is numerically investigated. To do this, a physical boundary condition in the interface of porous media and clear region of the channel is derived. The objective of this work is, first, to assess the effects of local solid-to-fluid heat transfer (a criterion indicating on departure from local thermal equilibrium (LTE) condition), solid-to-fluid thermal conductivity ratio and porous layer thickness on convective heat transfer in steady condition inside a channel partially filled with porous media; second, to examine the impact of pulsatile flow on heat transfer in the same channel. The effects of LTNE condition and thermal conductivity ratio in pulsatile flow are also briefly discussed. It is observed that Nusselt number inside the channel increases when the problem is tending to LTE condition. Therefore, careless consideration of LTE may lead to overestimation of heat transfer. Solid-to-fluid thermal conductivity ratio is also shown to enhance heat transfer in constant porous media thickness. It is also revealed that an increase in the amplitude of pulsation may result in enhancement of Nusselt number, while Nusselt number has a minimum in a certain frequency for each value of amplitude.     

Direct Numerical Simulation of Flow in a Ribbed Channel

The issue of turbine lifetime is an important one, particularly for modern turbines operating at high temperature regimes. A cooling design such as ribs may achieve an improved lifetime and complex mechanisms of heat transfer need to be well studied. In this paper, a Direct Numerical Simulation (DNS) is presented for a 3-D channel flow with two square ribs on the lower wall. The full unsteady compressible Navier-Stokes equations are solved with an original hybrid finite difference/finite element scheme. The Reynolds number of the simulation is 7 000 based on the bulk velocity at the inlet and the channel height. The present study is mainly devoted to understand the mechanism of heat transfer at the wall through the topological analysis of the flow and the temperature flux. Results show that the large-scale structures generated by obstacles splash onto the lower surface and induce longitudinal vortices which enhance heat transfer at the wall. A comprehensive data base including 56 correlations was set up for testing and improving turbulence models for this complex, separated flow.     

Heat transfer and flow behavior around four staggered elliptic cylinders in cross flow

Experimental and numerical studies were carried out to investigate forced convection heat transfer and flow features around the downstream elliptic cylinder in four staggered cylinders in cross flow. The elliptic cylinders examined had an axis ratio (b/c) of 1:2, and they were arranged with zero angle of attack to the upstream flow. The present heat transfer measurements were obtained by heating only the downstream elliptic cylinder (test cylinder) under the condition of constant heat flux. The testing fluid was air and the Reynolds number based on the major axis length (c) was ranged from 4,000 to 45,570. The tested longitudinal spacing ratio (S_x/c) and the transversal spacing ratio (S_y/b) were in the ranges of 1.5 ≤ S_x/c ≤ 4.0 and 1.5 ≤ S_y/b ≤ 4.0, respectively. The air flow pattern and temperature fields around the four staggered elliptic cylinders were predicted by using CFD software package. Also, a flow visualization study was made to show the flow features around the elliptic cylinders. It was observed that Nu_m of the downstream elliptic cylinder in four staggered cylinders was higher than that of three in-line cylinders for all tested spacing ratios and Reynolds numbers except for Re = 4,000. It was clear that, at lower Reynolds number values (Re < 14,100), the average Nusselt number of the downstream elliptic cylinder in three staggered arrangement was higher than that of the downstream cylinder in four staggered arrangement for all tested spacing ratios. On the other hand, at Re > 14,100, the tested elliptic cylinder in four staggered arrangement had the higher values of the average Nusselt number. Moreover, in four staggered arrangement, the maximum average Nusselt number enhancement ratio (average Nusselt number of the tested downstream cylinder/average Nusselt number of a single elliptic cylinder) was found to be about 2.0, and was obtained for spacing ratios of S_x/c = 2.5, S_y/b = 2.5 and at Re = 32,000. Finally, the average Nusselt number of the tested cylinder in four staggered arrangement was correlated in terms of Reynolds number and cylinder spacing ratios.     

Experimental investigation of water laminar mixed-convection flow with sub-millimeter bubbles in a vertical channel

This paper describes flow and heat transfer characteristics of laminar mixed-convection flows of water with sub-millimeter bubbles in a vertical channel. We use thermocouples and a particle tracking velocimetry technique for the temperature and velocity measurements. The working fluid used is tap water, and hydrogen bubbles generated by electrolysis of the water are used as the sub-millimeter bubbles. The Reynolds number of the main flow ranges from 100 to 200. The ratio of the heat transfer coefficient with sub-millimeter-bubble injection to that without injection (the heat transfer coefficient ratio) ranges from 1.24 to 1.38. The heat transfer coefficient ratio decreases with the increase in the Reynolds number. We conclude from velocity measurements that this decrease is mainly caused by a decrease in the advection effect due to sub-millimeter bubbles.     

Flow and heat transfer characteristics behind vortex generators - A benchmark dataset

The velocity field and wall heat transfer distributions for internal flows in the presence of longitudinal vortices have been experimentally investigated. A transient method based on temperature measurement with thermochromic liquid crystals was used to obtain the heat transfer distribution behind a tetrahedral, full-body vortex generator. With the focus on the longitudinal vortices, the flow field was captured by a three-component particle image velocimetry system. Mean values as well as velocity fluctuations have been assessed. The combined investigation of heat transfer and flow field describes in detail the physical conditions. For a channel Reynolds number of 300,000 a dataset has been obtained, which can be used for validation of numerical models.     

Forced Convection Heat Transfer of Nanofluids in a Porous Channel

This article is concerned with the effects of flow and migration of nanoparticles on heat transfer in a straight channel occupied with a porous medium. Investigation of force convective heat transfer of nanofluids in a porous channel has not been considered completely in the literature and this challenge is generally considered to be an open research topic that may require more study. The fully developed flow and steady Darcy?CBrinkman?CForchheimer equation is employed in porous channel. The thermal equilibrium model is assumed between nanofluid and solid phases. It is assumed that the nanoparticles are distributed non-uniformly inside the channel. As a result the volume fraction distribution equation is also coupled with governing equations. The effects of parameters such as Lewis number, Schmidt number, Brownian diffusion, and thermophoresis on the heat transfer are completely studied. The results show that the local Nusselt number is decreased when the Lewis number is increased. It is observed that as the Schmidt number is increased, the wall temperature gradient is decreased and as a consequence the local Nusselt number is decreased. The effects of Lewis number, Schmidt number, and modified diffusivity ratio on the volume fraction distribution are also studied and discussed.     

Friction and Heat Transfer in a Laminar Separated Flow behind a Rectangular Step with Porous Injection or Suction

Results of a numerical study of a laminar separated flow behind a rectangular step on a porous surface with uniform injection or suction are described. Two cases are considered: an unconfined flow past a step and flow evolution in a confined channel (duct). It is shown that mass transfer on the surface causes strong changes in the flow structure and substantially affects the position of the reattachment point, as well as friction and heat transfer. More intense injection leads first to an increase in the separation-zone length and then to its rapid vanishing due to boundary-layer displacement. Vice versa, suction at high Reynolds numbers Re _s > 100 reduces the separation-zone length. The duct flow has a complicated distribution of friction and heat-transfer coefficients along the porous surface owing to the coupled effect of the transverse flow of the substance and changes in the main flow velocity due to mass transfer. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 18–28, January–February, 2006.     

Modeling Fluid Flow in Fractured Porous Media with the Interfacial Conditions Between Porous Medium and Fracture

One of the most popular models that has been applied to predict the fluid velocity inside the fracture with impermeable walls is the cubic law. It highlights that the mean flux along the fracture is proportional to the cubic of fracture aperture. However, for a fractured porous medium, the normal and tangential interface conditions between the fracture and porous matrix can change the velocity profile inside the fracture. In this paper, a correction factor is introduced for flow equation along the fracture by imposing the continuity of normal and tangential components of velocity at the interface between the fracture and porous matrix. As a result, the mean velocity inside the fracture depends not only on the fracture aperture, but also on a set of non-dimensional numbers, including the matrix porosity, the ratio of intrinsic permeability of fracture to that of matrix, the wall Reynolds number, and the ratio of normal velocity on one wall to the other. Finally, the introduced correction factor is employed within the extended finite element method, which is widely used for numerical simulation of fluid flow within the fractured porous media. Several numerical results are presented for the fluid flow through a specimen containing single fracture, in order to investigate the deviation from the cubic law in different case studies.