Normal and tangential stress fluctuations during jamming

Under the application of high shear rates, concentrated suspensions of colloidal particles jam. The stress necessary to induce their flow abruptly increases and exhibits very large fluctuations. At the transition, the first normal stress difference increases and develops fluctuations. The sign of the first normal stress difference changes with the volume fraction of the suspension, which becomes dilatant at high volume fraction. I review an experimental study of this jamming transition and experimental evidence of the dilatancy of the suspension.     

Numerical simulation and visualization of fiber suspension in a turbulent round jet

The turbulent properties of the fiber suspension in a turbulent round jet are numerically simulated and visualized, and some of the results are compared with the experimental data. The effects of the Reynolds number, fiber volume fraction, and aspect ratio are analyzed. The results show that the fiber injection in the flow has a delay effect on the streamwise velocity decay along the jet axis, and such an effect becomes more obvious with the increases in the fiber volume fraction and aspect ratio and the decrease in the Reynolds number. The flow with fibers shows an increase in the streamwise velocity along the radial direction, and the increase magnitude is directly proportional to the fiber volume fraction and aspect ratio and inversely proportional to the Reynolds number. The presence of fibers makes the turbulent kinetic energy and Reynolds stress increase, and the extent increases with the fiber volume fraction, Reynolds number, and fiber aspect ratio.     

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.     

Rheology of an emulsion of viscoelastic drops in steady shear

Steady shear rheology of a dilute emulsion with viscoelastic inclusions is numerically investigated using direct numerical simulations. Batchelor's formulation for rheology of a viscous emulsion is extended for a viscoelastic system. Viscoelasticity is modeled using the Oldroyd-B constitutive equation. A front-tracking finite difference code is used to numerically determine the drop shape, and solve for the velocity and stress fields. The effective stress of the viscoelastic emulsion has three different components due to interfacial tension, viscosity difference (not considered here) and the drop phase viscoelasticity. The interfacial contributions – first and second normal stress differences and shear stresses – vary with Capillary number in a manner similar to those of a Newtonian system. However the shear viscosity decreases with viscoelasticity at low Capillary numbers, and increases at high Capillary numbers. The first normal stress difference due to interfacial contribution decreases with increasing drop phase viscoelasticity. The first normal stress difference due to the drop phase viscoelasticity is found to have a complex dependence on Capillary and Deborah numbers, in contrast with the linear mixing rule. Drop phase viscoelasticity does not contribute significantly to effective shear viscosity of the emulsion. The total first normal stress difference shows an increase with drop phase viscoelasticity at high Capillary numbers. However at low Capillary numbers, a non-monotonic behavior is observed. The results are explained by examining the stress field and the drop shape.     

Investigation of the stability of a plane-channel suspension flow with account for finite particle volume fraction

Within the framework of the two-fluid approach, a variant of a heterogeneous-medium model which takes into account a finite volume fraction of the inclusions and a small but finite phase velocity slip is proposed. The interphase momentum exchange is described by the Stokes force with the Brinkman correction for the finite particle volume fraction. The suspension viscosity depends on the particle volume fraction in accordance with the Einstein formula. Within the framework of the model constructed, a formulation of the problem of linear stability of plane-parallel two-phase flows is proposed. As an example, the stability of a channel suspension flow is considered. The system of equations for small disturbances with the boundary conditions is reduced to an eigenvalue problem for a fourth-order ordinary differential equation. Using the orthogonalization method, the dependence of the critical Reynolds number on the governing nondimensional parameters of the problem is studied numerically. It is shown that taking a finite volume fraction of the inclusions into account significantly affects the laminar-turbulent transition limit.     

Effect of viscosity ratio on the motion of drops flowing on an inclined surface

The flow of two-dimensional drops on an inclined channel is studied by numerical simulations at finite Reynolds numbers. The effect of viscosity ratio on the behaviour of the two-phase medium is examined. The flow is driven by the acceleration due to gravity, and there is no pressure gradient along the flow direction. An implicit version of the finite difference/front-tracking method was developed and used in the present study. The lateral migration of a drop is studied first. It is found that the equilibrium position of a drop moves away from the channel floor as the viscosity ratio increases. However, the trend reverses beyond a certain viscosity ratio. Simulations with 40 drops in a relatively large channel show that there exists a limiting viscosity ratio where the drops behave like solid particles, and the effect of internal circulation of drops becomes negligible. This limiting condition resembles the granular flow regime except that the effect of interstitial fluid is present. The limiting viscosity ratio depends on the flow conditions (80 for \(Re=10\), and 200 for \(Re=20\)). There are two peaks in the areal fraction distribution of drops across the channel which is different from granular flow regime. It is also found that the peak in areal fraction distribution of drops moves away from the channel floor as the inclination angle of the channel increases.     

Comparative study of lattice‐Boltzmann and finite volume methods for the simulation of laminar flow through a 4:1 planar contraction

In the present paper, a comparative study of numerical solutions for Newtonian fluids based on the lattice‐Boltzmann method (LBM) and the classical finite volume method (FVM) is presented for the laminar flow through a 4:1 planar contraction at a Reynolds number of value one, Re=1. In this study, the stress field for LBM is directly obtained from the distribution function. The calculations of the stress based on the FVM‐data use the evaluations of velocity gradients with finite differences. The stress field for both LBM and FVM is expressed in the present study in terms of the shear stress and the first normal stress difference. The lateral and axial profiles of the velocity, the shear stress and the first normal stress difference for both methods are investigated. It is shown that the LBM results for the velocity and the stresses are in excellent agreement with the FVM results. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical study of particulate suspension flow through wavy‐walled channels

The particulate suspension flow in a channel whose walls describe a travelling wave motion is examined numerically. A perturbation method is employed and the primitive variables are expanded in a series with the wall amplitude as the perturbation parameter. The boundary conditions are applied at the mean surface of the channel and the first‐order perturbation quantities are numerically determined by solving the governing system of ordinary differential equations by shooting technique. The present approach does not impose any restriction on the Reynolds number of the flow and the wave number and frequency of the wavy‐walled channel, although it is limited by the linear analysis. The wall shear stress and the positions of flow separation and reattachment points are computed and the influence of the volume fraction density of the particles is examined. The variations of velocity and pressure of the particulate suspension flow with frequency of excitation are also presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of turbulent MHD Couette nanofluid flow and heat transfer using hybrid DTM–FDM

Unsteady turbulent magnetohydrodynamic nanofluid hydrothermal treatment is studied. The zero-equation turbulence model is used to simulate turbulent flow. The modeling results obtained by applying the hybrid differential transformation method–finite difference method to solve this problem confirm its viability. An analytical procedure is used for finding the effects of the problem parameters. Results indicate that the average Nusselt number over the lower plate depends linearly on volume fraction of nanofluid, Hall parameter, turbulent Eckert number, and Reynolds number whereas it is inversely proportional on the Hartmann number and the turbulent parameter.     

Orientation distribution of fibers and rheological property in fiber suspensions flowing in a turbulent boundary layer

A model relating the translational and rotational transport of orientation distribution function （ODF） of fibers to the gradient of mean ODF and the dispersion coefficients is proposed to derive the mean equation for the ODE Then the ODF of fibers is predicted by numerically solving the mean equation for the ODF together with the equations of turbulent boundary layer flow. Finally the shear stress and first normal stress difference of fiber suspensions are obtained. The results, some of which agree with the available relevant experimental data, show that the most fibers tend to orient to the flow direction. The fiber aspect ratio and Reynolds number have significant and negligible effects on the orientation dis- tribution of fibers, respectively. The additional normal stress due to the presence of fibers is anisotropic. The shear stress of fiber suspension is larger than that of Newtonian solvent, and the first normal stress difference is much less than the shear stress. Both the additional shear stress and the first normal stress difference increase with increasing the fiber concentration and decreasing fiber aspect ratio.     

Three-dimensional numerical simulation of sedimenting drops inside a vertical channel

Numerical simulation of sedimenting deformable drops inside a vertical channel has been performed at finite Reynolds numbers. The channel is confined by two vertical walls in x-direction and is periodic in y- and z-directions. Results are obtained using a finite difference/front-tracking method. The main dimensionless parameters are the Reynolds number, the Bond number and ratio of the length of channel to the diameter of drops. The effect of these parameters on lateral migration of drops is investigated. It is found that the wall repulsion is the main mechanism of the lateral migration of the drop, and drop migrates toward the channel axis. When the Reynolds number is relatively low, two different lateral migration regimes are observed: migration with monotonic approach and migration with damped oscillations. These regimes are affected by the dimensionless parameters. When the Bond number increases, the oscillations of drop around the centerline of channel are stronger and drop reaches the channel centerline in a larger period. Results of lateral migration of one drop are consistent with perturbation theory, and two-dimensional numerical simulations performed by Feng et al. (1994). The drag coefficient has also been calculated, and effect of various parameters has been discussed. Two drops interaction is similar to that observed by Feng et al. (1994) for two-dimensional circular cylinders. Results are consistent with experiments performed by Wu and Manasseh (1998). Simulations of four sedimenting drops show that depending on the relative size of drops, they either fall in two rows or they form a single horizontal layer and settle with a unique velocity.     

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.     

Analysis of a meshless solver for high Reynolds number flow

In this paper, the extension of an upwind least‐square based meshless solver to high Reynolds number flow is explored, and the properties of the meshless solver are analyzed both theoretically and numerically. Existing works have verified the meshless solver mostly with inviscid flows and low Reynolds number flows, and in this work, we are interested in the behavior of the meshless solver for high Reynolds number flow, especially in the near‐wall region. With both theoretical and numerical analysis, the effects of two parameters on the meshless solver are identified. The first one is the misalignment effect caused by the significantly skewed supporting points, and it is found that the meshless solver still yields accurate prediction. It is a very interesting property and is opposite to the median‐dual control volume based vertex‐centered finite volume method, which is known to give degraded result with stretched triangular/tetrahedral cells in the near‐wall region. The second parameter is the curvature, and according to theoretical analysis, it is found in the region with both large aspect ratio and curvature, and the streamwise residual is less affected; however, the wall‐normal counterpart suffers from accuracy degradation. In this paper, an improved method that uses a meshless solver for the streamwise residual and finite difference for wall‐normal residual is developed. This method is proved to be less sensitive to the curvature and provides improved accuracy. This work presents an understanding of the meshless solver for high Reynolds number flow computation, and the analysis in this paper is verified with a series of numerical experiments. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Dissipative particle dynamics simulation of polymer drops in a periodic shear flow

The steady-state and transient shear flow dynamics of polymer drops in a microchannel are investigated using the dissipative particle dynamics (DPD) method. The polymer drop is made up of 10% DPD solvent particles and 90% finite extensible non-linear elastic (FENE) bead spring chains, with each chain consisting of 16 beads. The channel’s upper and lower walls are made up of three layers of DPD particles, respectively, perpendicular to Z-axis, and moving in opposite directions to generate the shear flow field. Periodic boundary conditions are implemented in the X and Y directions. With FENE chains, shear thinning and normal stress difference effects are observed. The “colour” method is employed to model immiscible fluids according to Rothman–Keller method; the χ-parameters in Flory–Huggins-type models are also analysed accordingly. The interfacial tension is computed using the Irving–Kirkwood equation. For polymer drops in a steady-state shear field, the relationship between the deformation parameter (D_def) and the capillary number (Ca) can be delineated into a linear and nonlinear regime, in qualitative agreement with experimental results of Guido et al. [J. Rheol. 42 (2) (1998) 395]. In the present study, Ca<0.22, in the linear regime. As the shear rate increases further, the drop elongates; a sufficiently deformed drop will break up; and a possible coalescence may occur for two neighbouring drops. Dynamical equilibrium between break-up and coalescence results in a steady-state average droplet-size distribution. In a shear reversal flow, an elongated and oriented polymer drop retracts towards a roughly spherical shape, with a decrease in the first normal stress difference. The polymer drop is found to undergo a tumbling mode at high Schmidt numbers. A stress analysis shows that the stress response is different from that of a suspension of solid spheres. An overshoot in the strain is observed for the polymer drop under extension due to the memory of the FENE chains.     

Finite difference and finite element methods for mhd channel flows

A Galerkin finite element method and two finite difference techniques of the control volume variety have been used to study magnetohydrodynamic channel flows as a function of the Reynolds number, interaction parameter, electrode length and wall conductivity. The finite element and finite difference formulations use unequally spaced grids to accurately resolve the flow field near the channel wall and electrode edges where steep flow gradients are expected. It is shown that the axial velocity profiles are distorted into M-shapes by the applied electromagnetic field and that the distortion increases as the Reynolds number, interaction parameter and electrode length are increased. It is also shown that the finite element method predicts larger electromagnetic pinch effects at the electrode entrance and exit and larger pressure rises along the electrodes than the primitive-variable and streamfunction–vorticity finite difference formulations. However, the primitive-variable formulation predicts steeper axial velocity gradients at the channel walls and lower axial velocities at the channel centreline than the streamfunction–vorticity finite difference and the finite element methods. The differences between the results of the finite difference and finite element methods are attributed to the different grids used in the calculations and to the methods used to evaluate the pressure field. In particular, the computation of the velocity field from the streamfunction–vorticity formulation introduces computational noise, which is somewhat smoothed out when the pressure field is calculated by integrating the Navier–Stokes equations. It is also shown that the wall electric potential increases as the wall conductivity increases and that, at sufficiently high interaction parameters, recirculation zones may be created at the channel centreline, whereas the flow near the wall may show jet-like characteristics.     

DRAG FORCE IN DENSE GAS—PARTICLE TWO—PHASE FLOW

Numerical simulations of flow over a stationary particle in a dense gas-particle two-phase flow have been carried out for small Reynolds numbers (less than 100). In order to study the influence of the particles interaction on the drag force, three particle arrangements have been tested: a single particle, two particles placed in the flow direction and many particles located regularly in the flow field. The Navier-Stokes equations are discretized in the three-dimensional space using finite volume method. For the first and second cases, the numerical results agree reasonably well with the data in literature. For the third case, i.e., the multiparticle case, the influence of the particle volume fraction and Reynolds numbers on the drag force has been investigated. The results show that the computational values of the drag ratio agree approximately with the published results at higher Reynolds numbers (from 34.2 to 68.4), but there is a large difference between them at small Reynolds numbers. The project supported by the Special Funds for Major State Basis Research Projects in China (G19990222).     

The effect of real viscosity on the heat transfer of water based Al2O3 nanofluids in a two-sided lid-driven differentially heated rectangular cavity

The heat transfer and fluid flow behavior of water based Al₂O₃ nanofluids are numerically investigated inside a two-sided lid-driven differentially heated rectangular cavity. Physical properties which have major effects on the heat transfer of nanofluids such as viscosity and thermal conductivity are experimentally investigated and correlated and subsequently used as input data in the numerical simulation. Transport equations are numerically solved with finite volume approach using SIMPLEC algorithm. It was found that not only the thermal conductivity but also the viscosity of nanofluids has a key role in the heat transfer of nanofluids. The results show that at low Reynolds number, increasing the volume fraction of nanoparticles increases the viscosity and has a deteriorating effect on the heat transfer of nanofluids. At high Reynolds number, the increase in the viscosity is compensated by force convection and the increase in the volume fraction of nanoparticles which results in an increase in heat transfer is in coincidence with experimental results.     

Investigation of Coulomb force effects on ethylene glycol based nanofluid laminar flow in a porous enclosure

Forced convection heat transfer of ethylene glycol based nanofluid with Fe_3O_4 inside a porous medium is studied using the electric field. The control volume based finite element method(CVFEM) is selected for numerical simulation. The impact of the radiation parameter(R_d), the supplied voltage(?φ), the volume fraction of nanofluid(?), the Darcy number(Da), and the Reynolds number(Re) on nanofluid treatment is demonstrated. Results prove that thermal radiation increases the temperature gradient near the positive electrode. Distortion of isotherms increases with the enhance of the Darcy number and the Coulomb force.     

Peristaltic flow of MHD Jeffrey fluid through finite length cylindrical tube

The present investigation studies the peristaltic flow of the Jeffrey fluid through a tube of finite length. The fluid is electrically conducting in the presence of an applied magnetic field. Analysis is carried out under the assumption of long wavelength and low Reynolds number approximations. Expressions of the pressure gradient, volume flow rate, average volume flow rate, and local wall shear stress are obtained. The effects of relaxation time, retardation time, Hartman number on pressure, local wall shear stress, and mechanical efficiency of peristaltic pump are studied. The reflux phenomenon is also investigated. The case of propagation of a non-integral number of waves along the tube walls, which are inherent characteristics of finite length vessels, is also examined.