Magnetohydrodynamic (MHD) flow of a second grade fluid in a channel with porous wall

An analysis has been carried out to study the effect of magnetic field on an electrically conducting fluid of second grade in a parallel channel. The coolant fluid is injected into the porous channel through one side of the channel wall into the other heated impermeable wall. The combined effect of inertia, viscous, viscoelastic and magnetic forces are studied. The basic equations governing the flow and heat transfer are reduced to a set of ordinary differential equations by using appropriate transformations for velocity and temperature. Numerical solutions of these equations are obtained with the help of Runge-Kutta fourth order method in association with quasi-linear shooting technique. Numerical results for velocity field, temperature field, skin friction and Nusselt number are presented in terms of elastic parameter, Hartmann number, Prandtl number and Reynolds number. Special case of our results is in good agreement with earlier published work.     

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.     

Combined influence of wall conductance,pressure work,viscous dissipation and joule heating on MHD channel flow heat transfer

To investigate the combined influence of viscous dissipation, pressure work, Joule heating, arbitrary voltage ratio, unequal wall conductances and wall heat fluxes on the fully developed laminar MHD channel flow heat transfer, the exact solution of the energy equations for fluid and channel walls are derived assuming the Hartmann velocity profile. It is concluded that there can be a substantial difference, depending upon Hartmann number, electric field and Brinkman number, between the Nusselt number considering the wall conductance and that neglecting it. Representative results are presented in diagrams.     

Thermofluidynamic analysis of the flow in a sharp 180° turn channel

The aim of the present study is to obtain surface flow visualisation, as well as local and spanwise averaged heat transfer measurements near a 180° sharp turn in a rectangular channel. The channel aspect ratio (width to height ratio) varies from 1 to 5 and the ratio between the width of the channel and that of the partition wall is always equal to 5. Heat transfer measurements are performed by means of the heated-thin-foil technique, which practically corresponds to a constant heat flux boundary condition, and by using infrared (IR) thermography. Two different heating conditions, in particular heating from one side (asymmetrical), or from two sides (symmetrical), are implemented. The convective heat transfer coefficient is evaluated from the measured temperature maps and the local bulk temperature of the flow which is obtained by making a one-dimensional balance along the channel. Results are presented in terms of local, or averaged, Nusselt number which is normalised with the classical Dittus and Boelter correlation. The fluid used during the test is air and the Reynolds number, based on the flow average velocity and channel hydraulic diameter, is varied between 16,000 and 60,000.     

Analytic Approximate Solution for a Flow of a Second-Grade Viscoelastic Fluid in a Converging Porous Channel

The problem of a two-dimensional steady flow of a second-grade fluid in a converging porous channel is considered. It is assumed that the fluid is injected into the channel through one wall and sucked from the channel through the other wall at the same velocity, which is inversely proportional to the distance along the wall from the channel origin. The equations governing the flow are reduced to ordinary differential equations. The boundary-value problem described by the latter equations is solved by the homotopy perturbation method. The effects of the Reynolds and crossflow Reynolds number on the flow characteristics are examined.     

Unsteady flow and heat transfer of porous media sandwiched between viscous fluids

The problem of unsteady oscillatory flow and heat transfer of porous medin sandwiched between viscous fluids has been considered through a horizontal channel with isothermal wall temperatures. The flow in the porous medium is modeled using the Brinkman equation. The governing partial differential equations are transformed to ordinary differential equations by collecting the non-periodic and periodic terms. Closed-form solutions for each region are found after applying the boundary and interface conditions. The influence of physical parameters, such as the porous parameter, the frequency parameter, the periodic frequency parameter, the viscosity ratios, the conductivity ratios, and the Prandtl number, on the velocity and temperature fields is computed numerically and presented graphically. In addition, the numerical values of the Nusselt number at the top and bottom walls are derived and tabulated.     

A numerical investigation of wall effects in three‐dimensional,laminar flow over a backward facing step with a constant aspect and expansion ratio

The wall effects are investigated in the three‐dimensional laminar flow over a backward‐facing step. For this purpose, a numerical experiment is designed under actual laboratory conditions. The aspect ratio of the computational domain is 1:40 and the expansion ratio is 1:2. The Reynolds number ranges from 100 to 950. The governing equations are the steady state, isothermal and incompressible Navier–Stokes equations for Newtonian fluids. They are solved with a homemade Galerkin finite element code. The computations are validated with data from available laboratory and numerical experiments. The results focus on the variation of both velocity profiles and lengths of eddies along the lower and upper wall in the spanwise direction. Calculated streamlines in the streamwise and transverse direction show how the flow is distorted near the lateral wall and how it develops up to the plane of symmetry. The study of skin friction lines along the top and bottom wall of the domain reveals a flow that takes place in the spanwise direction. This spanwise component of the flow becomes more dominant with increasing Reynolds number and is impossible to be sustained at steady state for Reynolds numbers higher than 950 for this particular geometry. Copyright © 2012 John Wiley & Sons, Ltd.     

Reynolds stress modelling of rectangular open‐channel flow

A Reynolds stress model for the numerical simulation of uniform 3D turbulent open‐channel flows is described. The finite volume method is used for the numerical solution of the flow equations and transport equations of the Reynolds stress components. The overall solution strategy is the SIMPLER algorithm, and the power‐law scheme is used to discretize the convection and diffusion terms in the governing equations. The developed model is applied to a flow at a Reynolds number of 77000 in a rectangular channel with a width to depth ratio of 2. The simulated mean flow and turbulence structures are compared with measured and computed data from the literature. The computed flow vectors in the plane normal to the streamwise direction show a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface as well as the free surface and bottom vortices. This small vortex causes a significant increase in the wall shear stress in the vicinity of the free surface. A budget analysis of the streamwise vorticity is carried out. It is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation of vortical ideal fluid flow through curved channel

A numerical algorithm to study the boundary‐value problem in which the governing equations are the steady Euler equations and the vorticity is given on the inflow parts of the domain boundary is developed. The Euler equations are implemented in terms of the stream function and vorticity. An irregular physical domain is transformed into a rectangle in the computational domain and the Euler equations are rewritten with respect to a curvilinear co‐ordinate system. The convergence of the finite‐difference equations to the exact solution is shown experimentally for the test problems by comparing the computational results with the exact solutions on the sequence of grids. To find the pressure from the known vorticity and stream function, the Euler equations are utilized in the Gromeka–Lamb form. The numerical algorithm is illustrated with several examples of steady flow through a two‐dimensional channel with curved walls. The analysis of calculations shows strong dependence of the pressure field on the vorticity given at the inflow parts of the boundary. Plots of the flow structure and isobars, for different geometries of channel and for different values of vorticity on entrance, are also presented. Copyright © 2003 John Wiley & Sons, Ltd.     

Turbulent natural convection in vertical parallel-plate channels

The problem of buoyancy driven turbulent flow in parallel-plate channels is investigated. The investigation is limited to vertical channels of uniform cross-section with different modes of heating. The details of the flow and thermal fields are obtained from the solution of the conservation equations of mass, momentum, and energy in addition to equations of the low Reynolds number turbulence model. The study covers Rayleigh number ranging from 10⁵ to 10⁷ and focuses on the effect of channel geometry on the characteristic of the flow and thermal fields as well as the local and average Nusselt number variation. A Nusselt number correlation has been developed in terms of a modified Rayleigh number and channel aspect ratio for the cases of symmetrically heated isothermal and isoflux conditions.     

Incompressible flow in a rapidly rotating cylinder with a source/sink distribution in the lateral wall and a free inner boundary

The flow of an incompressible fluid in a rapidly rotating right circular cylinder is considered. A source/sink mass distribution at the lateral wall, which is azimuthally uniform and symmetric across the midplane, causes a deviation from wheel flow. The container is only partially full and the inner free surface is allowed to deviate slightly from the vertical. A finite-difference solution of the full axisymmetric, non-linear governing equations was used to obtain the flow field. A special implicit technique for the Coriolis terms which maintains geostrophy was developed and is described. The results obtained for a low Rossby number flow compare quite favourably with the linearized solution. Results are also presented for a case wherein the non-linear terms are important.     

Microchannel membrane separation applied to confined thin film desorption

The concept of a confined thin film to enhance the desorption process is based on a reduced mass diffusion resistance. A wide thin film is formed into a microchannel by using a porous membrane as one wall of the channel enabling vapor extraction along the flow. Heat added to the channel results in vapor generation and subsequent extraction through the membrane. This experimental study investigates the performance of vapor extraction as a function of confined thin film thickness, pressure difference across the membrane and inlet concentration to the microchannel. In addition, heat added to the system was varied and results are presented in terms of the wall superheat temperature relative to the inlet saturated conditions of the binary fluid. The test section was equipped with a transparent window to observe bubble formation and vapor extraction. Results show that the performance, measured by the vapor release rate, increases for reduced channel thickness, for increased pressure difference across the membrane, and for lower inlet concentration. Results show that lower wall superheat correspond to higher heat transfer coefficients. Trends of Nusselt number and Sherwood number versus both channel Reynolds number and the product of the Reynolds number and Schmidt number are presented. Bubble formation in the channel does not degrade overall performance provided a critical heat flux condition does not occur.     

Low-Reynolds-number motion of a deformable drop between two parallel plane walls

The motion of a three-dimensional deformable drop between two parallel plane walls in a low-Reynolds-number Poiseuille flow is examined using a boundary-integral algorithm that employs the Green’s function for the domain between two infinite plane walls, which incorporates the wall effects without discretization of the walls. We have developed an economical calculation scheme that allows long-time dynamical simulations, so that both transient and steady-state shapes and velocities are obtained. Results are presented for neutrally buoyant drops having various viscosity, size, deformability, and channel position. For nearly spherical drops, the decrease in translational velocity relative to the undisturbed fluid velocity at the drop center increases with drop size, proximity of the drop to one or both walls, and drop-to-medium viscosity ratio. When deformable drops are initially placed off the centerline of flow, lateral migration towards the channel center is observed, where the drops obtain steady shapes and translational velocities for subcritical capillary numbers. With increasing capillary number, the drops become more deformed and have larger steady velocities due to larger drop-to-wall clearances. Non-monotonic behavior for the lateral migration velocities with increasing viscosity ratio is observed. Simulation results for large drops with non-deformed spherical diameters exceeding the channel height are also presented.     

Mixed Convection in a Vertical Porous Channel

A numerical study of mixed convection in a vertical channel filled with a porous medium including the effect of inertial forces is studied by taking into account the effect of viscous and Darcy dissipations. The flow is modeled using the Brinkman–Forchheimer-extended Darcy equations. The two boundaries are considered as isothermal–isothermal, isoflux–isothermal and isothermal–isoflux for the left and right walls of the channel and kept either at equal or at different temperatures. The governing equations are solved numerically by finite difference method with Southwell–Over–Relaxation technique for extended Darcy model and analytically using perturbation series method for Darcian model. The velocity and temperature fields are obtained for various porous parameter, inertia effect, product of Brinkman number and Grashof number and the ratio of Grashof number and Reynolds number for equal and different wall temperatures. Nusselt number at the walls is also determined for three types of thermal boundary conditions. The viscous dissipation enhances the flow reversal in the case of downward flow while it counters the flow in the case of upward flow. The Darcy and inertial drag terms suppress the flow. It is found that analytical and numerical solutions agree very well for the Darcian model. An erratum to this article is available at .     

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.     

Unsteady mixed convective heat transfer of two immiscible fluids confined between long vertical wavy wall and parallel flat wall

The combined effects of thermal and mass convection of viscous incompressible and immiscible fluids through a vertical wavy wall and a smooth flat wall are analyzed. The dimensionless governing equations are perturbed into a mean part (the zeroth-order) and a perturbed part (the first-order). The first-order quantities are obtained by the perturbation series expansion for short wavelength, in which the terms of the exponential order arise. The analytical expressions for the zeroth-order, the first-order, and the total solutions are obtained. The numerical computations are presented graphically to show the salient features of the fluid flow and the heat transfer characteristics. Separate solutions are matched at the interface by using suitable matching conditions. The shear stress and the Nusselt number are also analyzed for variations of the governing parameters. It is observed that the Grashof number, the viscosity ratio, the width ratio, and the conductivity ratio promote the velocity parallel to the flow direction. A reversal effect is observed for the velocity perpendicular to the flow direction.     

The effects of surface roughness geometry of flow undergoing Poiseuille flow by molecular dynamics simulation

Numerical simulation of Poiseuille flow of liquid Argon in a rough nano-channel using the non-equilibrium molecular dynamics simulation is performed. Density and velocity profiles across the channel are investigated in which roughness is implemented only on the lower wall. The Lennard–Jones potential is used to model the interactions between all particles. The effects of surface roughness geometry, gap between roughness elements (or roughness periodicity), surface roughness height and surface attraction energy on the behavior of the flow undergoing Poiseuille flow are presented. Results show that surface shape and roughness height have a decisive role on the flow behaviors. In fact, by increasing the roughness ratio (height to base ratio), the slip velocity and the maximum velocity in the channel cross section are reduced, and the density fluctuations near the wall increases. Results also show that the maximum density near the wall for a rough surface is less than a smooth wall. Moreover, the simulation results show that the effect of triangle roughness surface on the flow behavior is more than the cylindrical ones.     

Thermophoretic augmentation of particle deposition in natural convection flow through a parallel plate channel

Present paper deals with temperature driven mass deposition rate of particles known as thermophoretic wall flux when a hot flue gas in natural convection flow through a cooled isothermal vertical parallel plate channel. Present study finds application in particle filters used to trap soot particles from post combustion gases issuing out of small furnaces with low technical implications. Governing equations are solved using finite difference marching technique with channel inlet values as initial values. Channel heights required to regain hydrostatic pressure at the exit are estimated for various entry velocities. Effect of temperature ratio between wall and gas on thermophoretic wall flux is analysed and wall flux found to increase with decrease in temperature ratio. Results are compared with published works wherever possible and can be used to predict particle deposition rate as well as the conditions favourable for maximum particle deposition rate.     

Numerical Simulation of Turbulent Flow and Heat Transfer in a Three-Dimensional Channel Coupled with Flow Through Porous Structures

This study investigates numerically the turbulent flow and heat transfer characteristics of a T-junction mixing, where a porous media flow is vertically discharged in a 3D fully developed channel flow. The fluid equations for the porous medium are solved in a pore structure level using an Speziale, Sarkar and Gatski turbulence model and validated with open literature data. Overall, two types of porous structures, consisted of square pores, are investigated over a wide range of Reynolds numbers: an in-line and a staggered pore structure arrangement. The flow patterns, including the reattachment length in the channel, the velocity field inside the porous medium as well as the fluctuation velocity at the interface, are found to be strongly affected by the velocity ratio between the transversely interacting flow streams. In addition, the heat transfer examination of the flow domain reveals that the temperature distribution in the porous structure is more uniform for the staggered array. The local heat transfer distributions inside the porous structure are also studied, and the general heat transfer rates are correlated in terms of area-averaged Nusselt number accounting for the effects of Reynolds number, velocity ratio as well as the geometrical arrangement of the porous structures.     

Effect of span length and temperature on the 3-D confined flow around a vortex promoter

This article presents a numerical study on the influence of span length and wall temperature on the 3-D flow pattern around a square section vortex promoter located inside a micro-channel in the low Reynolds number regime. The first objective of the work is to quantify the critical Reynolds number that defines the onset of vortex shedding and to identify the different regimes that appear as a function of the channel aspect ratio (span to height ratio). We found that the critical Reynolds number for the onset of the Karman street regime increases as the aspect ratio decreases. In particular, for the aspect ratio of 1/2 the critical Reynolds number is nearly six times the critical Reynolds number of the 2-D problem. An intermediate oscillating regime between the steady and the Karman street solutions was also found to exist within a rather wide range of Reynolds numbers for small channel aspect ratios. The second objective was to investigate the influence of the vortex promoter wall temperature on both vortex shedding and flow pattern. This has practical engineering implications because the working fluid considered in the article is water that has a viscosity that depends significantly on temperature and promotes a strong coupling between the momentum and energy equations that influences the system behaviour. Results indicate that high surface temperature on the prism promotes the onset of the Karman street, suggesting design guidelines for micro-channel based heat sinks that make use of vortex promoters.