Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

Heat transfer effects on the stability of low speed plane Couette-Poiseuille flow

The stability problem of low-speed plane Couette-Poiseuille flow of air under heat transfer effects is solved numerically using the linear stability theory. Stability equations obtained from two-dimensional equations of motion and their boundary conditions result in an eigenvalue problem that is solved using an efficient shoot-search technique. Variable fluid properties are accounted for both in the basic flow and the perturbation (stability) equations. A parametric study is performed in order to assess the roles of moving wall velocity and heat transfer. It is found that the moving wall velocity and the location of the critical layers play decisive roles in the instability mechanism. The flow becomes unconditionally stable whenever the moving wall velocity exceeds half of the maximum velocity in the channel. With wall heating and Mach number effects included, the flow is stabilized.     

Transient natural convection flow of a particulate suspension through a vertical channel

Continuum equations governing transient, laminar, fully-developed natural convection flow of a particulate suspension through an infinitely long vertical channel are developed. The equations account for particulate viscous effects which are absent from the original dusty-gas model. The walls of the channel are maintained at constant but different temperatures. No-slip boundary conditions are employed for the particle phase at the channel walls. The general transient problem is solved analytically using trigonometric Fourier series and the Laplace transform method. A parametric study of some physical parameters involved in the problem is performed to illustrate the influence of these parameters on the flow and thermal aspects of the problem.     

A particle suspension model: an unstructured finite-volume implementation

In this paper, a modified particle temperature model for concentrated suspensions is proposed, which allows for the shear-induced migration of particles. The migration is modelled by a convection–diffusion equation, derived from the particle mass and momentum conservation. The model is implemented in an unstructured finite volume method and is utilized to investigate the shear-induced particle migration in channel flow. The profiles and the evolution of the velocity, concentration and particle temperature along the channel are presented. The entrance lengths needed to reach a fully developed profile of the corresponding field variables are also checked against different averaged concentrations and different relative particle radii. Comparison with available experimental data is made whenever possible.     

THREE-DIMENSIONAL VISCOUS FLOW THROUGH A ROTATING CHANNEL: A PSEUDOSPECTRAL MATRIX METHOD APPROACH

A Fourier–Chebyshev pseudospectral method is used for the numerical simulation of incompressible flows in a three-dimensio nal channel of square cross-section with rotation. Realistic, non-periodic boundary conditions that impose no-slip conditions in two directions (spanwis e and vertical directions) are used. The Navier–Stokes equations are integrated in time using a fractional step method. The Poisson equations for pressure and the Helmholtz equation for velocity are solved using a matrix diagonalization (eigenfunction decomposition) method, through which we are able to reduce a three-dimensional matrix problem to a simple algebraic vector equation. This results in signficant savings in computer storage requirement, particularly for large-scale computations. Verification of the numerical algorithm and code is carried out by comparing with a limiting case of an exact steady state solution for a one-dimensional channel flow and also with a two-dimensional rotating channel case. Two-cell and four-cell two-dimensional flow patterns are observed in the numerical experiment. It is found that the four-cell flow pattern is stable to symmetri cal disturbances but unstable to asymmetrical disturbances.     

A confrontation of 1D and 2D RKDG numerical simulation of transitional flow at open‐channel junction

In this study, a comparison between the 1D and 2D numerical simulation of transitional flow in open‐channel networks is presented and completely described allowing for a full comprehension of the modeling water flow. For flow in an open‐channel network, mutual effects exist among the channel branches joining at a junction. Therefore, for the 1D study, the whole system (branches and junction) cannot be treated individually. The 1D Saint Venant equations calculating the flow in the branches are then supplemented by various equations used at the junction: a discharge flow conservation equation between the branches arriving and leaving the junction, and a momentum or energy conservation equation. The disadvantages of the 1D study are that the equations used at the junction are of empirical nature due to certain parameters given by experimental results and moreover they often present a reduced field of validity. On the contrary, for the 2D study, the entire network is considered as a single unit and the flow in all the branches and junctions is solved simultaneously. Therefore, we simply apply the 2D Saint Venant equations, which are solved by a second‐order Runge–Kutta discontinuous Galerkin method. Finally, the experimental results obtained by Hager are used to validate and to compare the two approaches 1D and 2D. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Entropy generation in flow field subjected to a porous block in a vertical channel

Entropy generation in the flow field subjected to a porous block situated in a vertical channel is examined. The effects of channel inlet port height (vertical height between channel inlet port and the block center), porosity, and block aspect ratio on the entropy generation rate due to fluid friction and heat transfer in the fluid are examined. The governing equations of flow, heat transfer, and entropy are solved numerically using a control volume approach. Air is used as the flowing fluid in the channel. A uniform heat flux is considered in the block and natural convection is accommodated in the analysis. It is found that entropy generation rate due to fluid friction increases with increasing inlet port height, while this increase becomes gradual for entropy generation rate due to heat transfer for the inlet port height exceeding 0.03 m. The porosity lowers entropy generation rate due to fluid friction and heat transfer. The effect of block aspect ratio on entropy generation rate is notable; in which case, entropy generation rate increases for the block aspect ratio of 1:2.     

Droplet break-up through an oblique shock wave

The interaction of a two-phase flow with a wedge where a stationary shock wave is initially settled is studied in a two-dimensional configuration. Before the introduction of the dispersed phase, the flow around the wedge is a supersonic one phase flow such as an attached stationary shock wave is present. Then, the dispersed phase is introduced upstream the initial position of the stationary shock wave. The purpose of this study is to point out two-phase and droplets break-up effects on the oblique shock wave. The two-dimensional equations are solved by a TVD scheme where fluxes are computed by using Riemann solver for the gas phase equations and also for the dispersed phase equations wich is an original approach due to the authors (Saurel et al. 1994). In addition to drag forces and heat and mass transfers, the process of droplets fragmentation based on the particle oscillation is considered. Accepted April 28, 1995     

Heat Transfer Characteristics of Reciprocating Flows in Channels Partially Filled with Porous Medium

This paper presents the effect of introducing a porous medium on the flow regime and heat transfer of a two-dimensional channel through which the flow is reciprocating. The channel is discretely heated from above and is insulated in the bottom which can simulate a cooling mechanism for compact circuit boards. In this ideal geometry, a fully developed reciprocating flow is established via oscillating pressure gradient. In side boundaries, velocity and temperature are assumed to be periodic. A certain volume of this channel is occupied by a porous medium which is shown to be an effecting tool for augmentation of heat transfer. At first, Momentum equations of the domain are solved analytically (Brinkman-extended Darcy model is used for porous region) and then the energy equation is solved numerically using alternating direction implicit (ADI) method. Finally a case study is investigated for a high-porous and high-conductive medium (Aluminum alloy T-6201) and the enhancing effect and optimization criteria are discussed in the result section.     

Scaling Transformation for the Effect of Temperature-Dependent Fluid Viscosity with Thermophoresis Particle Deposition on MHD-Free Convective Heat and Mass Transfer Over a Porous Stretching Surface

This article concerns with a steady two-dimensional flow of an electrically conducting incompressible fluid over a vertical stretching sheet. The flow is permeated by a uniform transverse magnetic field. The fluid viscosity is assumed to vary as a linear function of temperature. A scaling group of transformations is applied to the governing equations. The system remains invariant due to some relations among the parameters of the transformations. After finding three absolute invariants, a third-order ordinary differential equation corresponding to the momentum equation, and two second-order ordinary differential equations corresponding to energy and diffusion equations are derived. The equations along with the boundary conditions are solved numerically. It is found that the decrease in the temperature-dependent fluid viscosity makes the velocity to decrease with the increasing distance of the stretching sheet. At a particular point of the sheet, the fluid velocity decreases with the decreasing viscosity but the temperature increases in this case. Impact of thermophoresis particle deposition in the presence of temperature-dependent fluid viscosity plays an important role on the concentration boundary layer. The results, thus, obtained are presented graphically and discussed.     

Motion of a circular cylinder in a viscous liquid between parallel plates

The two-dimensional motion of a cylinder in a viscous fluid between two parallel walls of a vertical channel is studied. It is found that when the cylinder moves very closely along one of the channel walls, it always rotates in the direction opposite to that of contact rolling along the nearest wall. When the cylinder is away from the walls, its rotation depends on the Reynolds number of the flow. In this study two numerical methods were used. One is for the unsteady motion of a sedimenting cylinder initially released from a position close to one of the channel walls, where the Navier-Stokes equations are solved for the fluid and Newton's equations of motion are solved for the rigid cylinder. The other method is for the steady flow in which a cylinder is fixed in a uniform flow field where the channel walls are sliding past the cylinder at the speed of the approaching flow, or equivalently a cylinder is moving with a constant velocity in a quiescent fluid. The flow field, the drag, the side force (lift), and the torque experienced by the cylinder are studied in detail. The effects of the cylinder location in the channel, the size of the channel relative to the cylinder diameter, and the Reynolds number of the flow are examined. In the limit when the cylinder is translating very closely along one of the walls, the flow in the gap between the cylinder and the wall is solved analytically using lubrication theory, and the numerical solution in the other region is used to piece together the whole flow field.This research was supported by NSF DMR91-20668 through the Laboratory for Research on the Structure of Matter at the University of Pennsylvania and from the Research Foundation of the University of Pennsylvania.     

SIMULATION OF THE FLOW OVER ELLIPTIC AIRFOILS OSCILLATING AT LARGE ANGLES OF ATTACK

The incompressible, viscous flow over two-dimensional elliptic airfoils oscillating in pitch at large angles of attack, such that flow separation occurs, has been simulated numerically for a Reynolds number of 3000. A vortex method is used to solve the two-dimensional Navier–Stokes equations in vorticity/stream-function form using a time-marching approach. Using an operator-splitting method the convection and diffusion equations are solved sequentially at each time step. The convection equation is solved using a vortex-in-cell method, and the diffusion equation using a second-order ADI finite-difference scheme. Elliptic profiles are obtained by mapping a circle in a computational domain into the physical domain using a Joukowski transformation. The effects of several parameters on the flow field are considered, such as: frequency of oscillation, mean angle of attack, location of pitch-axis and the thickness ratio of the ellipse. The results obtained are shown to compare favourably with available experimental results.     

Particle suspension in a simple shear flow

A hydrodynamic model describing the particle distribution over the cross-section of a finely dispersed flow is proposed. The model is constructed on the basis of notions concerning the diffusion of particles induced by their random displacements in the process of relative motion of neighboring layers at constant shear velocity. It is shown that the suspension capacity of the flow is large for small particles due to thermal fluctuations and for relatively large particles due to shear-induced particle pulsations. There are critical particle sizes for which the particles are suspended and transported by the flow less effectively than larger or smaller particles.Ekaterinburg. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 112–121, January–February, 1995.     

Rheology of a suspension of particles in viscoelastic fluids

In this paper we study the bulk stress of a suspension of rigid particles in viscoelastic fluids. We first apply the theoretical framework provided by Batchelor [J. Fluid Mech. 41 (1970) 545] to derive an analytical expression for the bulk stress of a suspension of rigid particles in a second-order fluid under the limit of dilute and creeping flow conditions. The application of the suspension balance model using this analytical expression leads to the prediction of the migration of particles towards the centerline of the channel in pressure-driven flows. This is in agreement with experimental observations. We next examine the effects of inertia (or flow Reynolds number) on the rheology of dilute suspensions in Oldroyd-B fluids by two-dimensional direct numerical simulations. Simulation results are verified by comparing them with the analytical expression in the creeping flow limit. It is seen that the particle contribution to the first normal stress difference is positive and increases with the elasticity of the fluid and the Reynolds number. The ratio of the first normal stress coefficient of the suspension and the suspending fluid decreases as the Reynolds number is increased. The effective viscosity of the suspension shows a shear-thinning behavior (in spite of a non-shear-thinning suspending fluid) which becomes more pronounced as the fluid elasticity increases.     

Flow-driven orientation dynamics of semiflexible fiber systems

We study the flow-induced orientation dynamics of semiflexible fibers in dilute fiber suspensions. Starting from the equations of motion for a two-rod model of flexible fibers in Stokes flow, the Smoluchowski equation for a connected monomer orientation distribution function is derived. We then obtain a set of equations for the time dependence of the first and second moments of the orientation distribution function, thus extending the Folgar Tucker equations for short rigid fiber suspensions to flexible fiber suspensions. The resulting generalized equations for the orientation dynamics of a suspension of flexible fibers are solved for simple channel flow. It is shown that all qualitative effects of bending and straightening of fibers and their influence on the orientation of flexible fibers are captured within our model. A scalar measure for the distribution of bending in a flow is introduced, which allows to detect the degree of bending of fibers. Paper was presented at the 3rd Annual Rheology Conference, AERC 2006, April 27–29, 2006, Crete, Greece.     

Simulation of dynamic stall for a NACA 0012 airfoil using a vortex method

The unsteady, incompressible, viscous laminar flow over a NACA 0012 airfoil is simulated, and the effects of several parameters investigated. A vortex method is used to solve the two-dimensional Navier–Stokes equations in the vorticity/stream-function form. By applying an operator-splitting method, the “convection” and “diffusion” equations are solved sequentially at each time step. The convection equation is solved using the vortex-in-cell method, and the diffusion equation using a second-order ADI finite difference scheme. The airfoil profile is obtained by mapping a circle in the computational domain into the physical domain through a Joukowski transformation. The effects of several parameters are investigated, such as the reduced frequency, mean angle of attack, location of pitch axis, and the Reynolds number. It is observed that the reduced frequency has the most influence on the flow field.     

Heat transfer analysis of a particle-containing channel flow

This paper describes an analytical model of heat transfer in a two-dimensional, steady, nonreacting particle-containing channel flow. An idealized gas flow of specified uniform velocity between insulated parallel plates is assumed and the nonvaporizing particles are conceptualized as contained within an thin sheet injected at the symmetry plane. Two dimensionless parameters that affect the solution are described. These are the effective gas diffusivityK and the dimensionless particle number densityP. The linear, coupled differential equations governing the energy exchange between the gas and liquid phases are solved by means of the Green's function technique. This procedure yields a Volterra integral-series equation as the solution of the gas-phase energy equation. A series solution of this integral equation is obtained by the method of successive substitutions and terms up to second order are calculated.     

SPH-based numerical treatment of the interfacial interaction of flow with porous media

In this paper, the macroscopic equations of mass and momentum are developed and discretized based on the smoothed particle hydrodynamics (SPH) formulation for the interaction at an interface of flow with porous media. The theoretical background of flow through porous media is investigated to highlight the key constraints that should be satisfied, particularly at the interface between the porous media flow and the overlying free flow. The study aims to investigate the derivation of the porous flow equations, computation of the porosity, and treatment of the interfacial boundary layer. It addresses weak assumptions that are commonly adopted for interfacial flow simulation in particle-based methods. As support to the theoretical analysis, a two-dimensional weakly compressible SPH model is developed based on the proposed interfacial treatment. The equations in this model are written in terms of the intrinsic averages and in the Lagrangian form. The effect of particle volume change due to the spatial change of porosity is taken into account, and the extra stress terms in the momentum equation are approximated by using Ergun's equation and the subparticle scale model to represent the drag and turbulence effects, respectively. Four benchmark test cases covering a range of flow scenarios are simulated to examine the influence of the porous boundary on the internal, interface, and external flows. The capacity of the modified SPH model to predict velocity distributions and water surface behavior is fully examined with a focus on the flow conditions at the interfacial boundary between the overlying free flow and the underlying porous media.     

A Numerical Method for 2-D Bed Morphology Calculations

A fully coupled two-dimensional sub-critical and/or supercritical, free-surface flow numerical model is developed to calculate bed variations in alluvial channels. Vertically averaged free-surface flow equations in conjunction with sediment transport equation are numerically solved using an explicit finite-volume scheme in integral form. The capabilities of the proposed method are first demonstrated by analyzing supercritical flow in an expansion channel. Thereafter, one and two-dimensional applications referring to aggradation and scouring are reported. For each of these test cases, computed results compare satisfactorily with measurements as well as with other numerical solutions. The method is stable, reliable and accurate, although time consuming, handling a variety of sediment transport equations with rapid changes of sediment transport at the boundaries.