Viscoplastic flow in slightly varying channels with wall slip pertaining to a magnetorheological (MR) polishing process

In this paper, the first analytical endeavor into the fluid dynamic modeling of an MR polishing process is reported. The velocity and shear stress fields of an MR fluid running through a thin slippery channel with a slightly varying height are analytically solved using a bi-viscosity constitutive model and a Navier slip model. Estimations of the mechanical power density and the total power per unit depth applied onto the channel surfaces are also presented. Analytical solutions for the Couette–Poiseuille flow behavior of a bi-viscous fluid flowing through either parallel or non-uniform channels are obtained, and the associated necessary and sufficient conditions characterizing a total of 5 types of flow are derived.The behaviors of the fluid are examined through the use of a parametric diagram of Bingham number (Bn) and Couette number (Co), i.e., Bn–Co or 1/Bn–1/Co diagram, by changing the geometric and operating conditions. Using these diagrams, variations in the rheological characteristics of the flow are investigated in great detail, with a special focus on the movement of the pseudo-core region. Finally, the mechanical power density field obtained for the flow in a converging–diverging channel is used to explain the wear mechanism in the MR polishing process. The effects of the power density field and the total power on the material removal rate (MRR) and the within-workpiece nonuniformity (WIWNU) with respect to various geometric and operating conditions are evaluated.     

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.     

Determination of yield stress fluid behaviour from inclined plane test

The aim of this paper is to determine precisely under which conditions an inclined plane can be used as a rheometer, which could represent a practical and rapid technique for various types of industrial or natural viscoplastic coarse suspensions. We first examine its efficiency and relevancy for determining fluid yield stress in a straight way by measuring the deepest fluid layer able to stay on the inclined plane. We have made experiments with different materials (clay-water suspensions) whose yield stress ranged from 35 to 90 Pa, using 1 m long open rectangular channels with a slope ranging from 10 to 30° and a width ranging from 5 to 25 cm. Our procedure involved measuring the final fluid depth far from edges a long time after the end of the slow gravity-induced emptying of a dam placed upstream. The fluid yield stress was also estimated independently by fitting a Herschel-Bulkley model to simple shear rheometry data obtained within a relatively wide shear rate range. A good agreement between inclined rectangular channel tests and independent usual rheometrical tests is obtained even for aspect ratios (flow depth to channel width ratio) as large as 1 when one assumes that, when the fluid has stopped, the side and bottom wall shear stresses are equal to the fluid yield stress. These results prove the efficiency of the inclined plane test for determining yield stress when appropriate experimental precautions are taken for both tests. In addition we examine the possibility of determining the simple shear flow curve of a mud suspension from fluid depth, velocity and discharge measurements of different steady flows in a wide open channel (8 m long; 60 cm wide) equipped with a recirculating system. The results obtained from inclined plane tests are in good agreement with independent rheometrical data (with torsional geometries). However it is technically difficult to cover a wide shear rate range from the inclined plane technique since this requires a rather wide channel flow rate range.     

Computation of Stokes Flow in a Channel with a Collapsible Segment

The numerical solution of Stokes flow in two-dimensional channel in which a segment of one wall is formed by an elastic membrane under longitudinal tension and the remaining channel boundary is rigid is considered. This model problem is being used to gain an understanding of the complex interactions that occurs between the fluid flow and the wall mechanics when fluid flows through a collapsible tube, examples of which are widespread in physiology. Previous work by Pedley considered a similar system using lubrication theory in which the wall slopes are assumed small. The results showed that as the longitudinal wall tension is reduce, the downstream end of the collapsible segment becomes ever steeper, thus violating the assumptions. Here, lubrication theory is abandoned and a numerical solution of the full governing equations, including the complete expression for wall curvature, is sought using an iterative scheme. The effect of the variation in wall tension due to the fluid shear stresses at the compliant boundary is also included.Results are presented for a range of transmural (internal minus external) pressures and wall tensions. It is found, however, that as the wall tension is reduced, the iterative scheme considered fails to converge. This similar behaviour to that seen by Silliman & Scriven in viscous free-surface flows. Possible reasons for this breakdown together with alternative solution strategies are discussed.     

Thin shear boundary layers in flow of hot aluminium

The flow of hot aluminium in channels is investigated. The constitutive relation considered for the flow stress is the inverse hyperbolic sine function of the Zener–Hollomon parameter. Analytical solutions for the flow are derived. At high strain rates, an exponential velocity profile close to the channel walls is predicted indicating the existence of a thin shear boundary layer characterized by strong shear. The characteristic length scale for the exponential velocity profile is a function of material parameters in the constitutive relation and the inverse of the local pressure gradient. The analytical prediction of a thin shear boundary layer close to the channel walls for large strain rates is consistent with the observed microstructure in an extruded section.     

Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

Interstitial fluid flow: simulation of mechanical environment of cells in the interosseous membrane

In vitro experiments have shown that subtle fluid flow environment plays a significant role in living biological tissues,while there is no in vivo practical dynamical measurement of the interstitial fluid flow velocity.On the basis of a new finding that capillaries and collagen fibrils in the interosseous membrane form a parallel array,we set up a porous media model simulating the flow field with FLUENT software,studied the shear stress on interstitial cells’ surface due to the interstitial fluid flow,and analyzed the effect of flow on protein space distribution around the cells.The numerical simulation results show that the parallel nature of capillaries could lead to directional interstitial fluid flow in the direction of capillaries.Interstitial fluid flow would induce shear stress on the membrane of interstitial cells,up to 30 Pa or so,which reaches or exceeds the threshold values of cells’ biological response observed in vitro.Interstitial fluid flow would induce nonuniform spacial distribution of secretion protein of mast cells.Shear tress on cells could be affected by capillary parameters such as the distance between the adjacent capillaries,blood pressure and the permeability coefficient of capillary’s wall.The interstitial pressure and the interstitial porosity could also affect the shear stress on cells.In conclusion,numerical simulation provides an effective way for in vivo dynamic interstitial velocity research,helps to set up the vivid subtle interstitial flow environment of cells,and is beneficial to understanding the physiological functions of interstitial fluid flow.     

Stable and accurate schemes for smoothed dissipative particle dynamics

Smoothed dissipative particle dynamics (SDPD) is a mesoscopic particle method that allows to select the level of resolution at which a fluid is simulated. The numerical integration of its equations of motion still suffers from the lack of numerical schemes satisfying all the desired properties such as energy conservation and stability. Similarities between SDPD and dissipative particle dynamics with energy (DPDE) conservation, which is another coarse-grained model, enable adaptation of recent numerical schemes developed for DPDE to the SDPD setting. In this article, a Metropolis step in the integration of the fluctuation/dissipation part of SDPD is introduced to improve its stability.     

Viscoelastic flow in a two-dimensional collapsible channel

We compute the flow of three viscoelastic fluids (Oldroyd-B, FENE-P, and Owens blood model) in a two-dimensional channel partly bounded by a tensioned membrane, a benchmark geometry for fluid–structure interactions. The predicted flow patterns are compared to those of a Newtonian liquid. We find that computations fail beyond a limiting Weissenberg number. Flow fields and membrane shape differ significantly because of the different degree of shear thinning and molecular extensibility underlying the three different microstructural models.     

A Two-Layer Mathematical Model of Blood Flow in Porous Constricted Blood Vessels

In this paper, we discussed a mathematical model for two-layered non-Newtonian blood flow through porous constricted blood vessels. The core region of blood flow contains the suspension of erythrocytes as non-Newtonian Casson fluid and the peripheral region contains the plasma flow as Newtonian fluid. The wall of porous constricted blood vessel configured as thin transition Brinkman layer over layered by Darcy region. The boundary of fluid layer is defined as stress jump condition of Ocha-Tapiya and Beavers–Joseph. In this paper, we obtained an analytic expression for velocity, flow rate, wall shear stress. The effect of permeability, plasma layer thickness, yield stress and shape of the constriction on velocity in core & peripheral region, wall shear stress and flow rate is discussed graphically. This is found throughout the discussion that permeability and plasma layer thickness have accountable effect on various flow parameters which gives an important observation for diseased blood vessels.     

The rheological behavior of flexible-chain polymers in the region of high shear rates and stresses,the critical process of spurting,and supercritical conditions of their movement at T > Tg

Methods of capillary viscometry were used in studying the rheological properties and behavior of a broad range of rubbers, including polymers with narrow and wide molecular-weight-distribution as well as commercial rubber grades, at widely varying shear rates and stresses. As is shown, in full conformity with the previously conducted experiments, during transition from a fluid to highelastic (quasi-cross-linked) state, they are chracterized by spurting followed by sliding over the channel walls. This relaxation transition is characterized by a critical shear stress value invariant with respect to the molecular weight, molecularweight distribution and temperature. The parameters defining spurting of polymer flow as a function of molecular-weight characteristics, temperature, and channel geometry have been investigated in detail. It is shown for the first time that under supercritical conditions the rate of polymer flow through channels does not depend, in the first approximation, on the molecular weight of the polymer, its molecularweight distribution, temperature, and filling, but is determined only by the shear stress.     

Study of wall properties on peristaltic transport of a couple stress fluid

In this paper, we investigate the peristaltic transport of a couple stress fluid in a channel with compliant walls. Perturbation method has been used to get the solution. The flow is induced by sinusoidal traveling waves along the channel walls. The effects of wall damping, wall elastance, wall tension and couple stress parameter on the flow are investigated using the equations of fluid as well as deformable boundaries. It is found that the mean velocity at boundaries decreases with increasing couple-stress parameter and wall damping and increases with increasing wall tension and wall elastance, while the mean axial velocity increases with increasing wall tension and wall elastance and decreases with couple-stress parameter and wall damping.     

Pressure driven flow of dilute polymer solutions in a channel: Analytic results for very small and for very large channels

Kinetic theory of dilute macromolecular solutions is applied to pressure driven flow in a small channel where wall- (and interfacial) layers have to be reckoned with. The complete rheology is studied. It turns out that for very small channels both the shear stress and the normal stress are an order of magnitude larger than corresponding quantities in simple shear. On the other hand, when the channel is so wide that the wall layers are very thin in comparison, agreement with results appropriate for simple shear is found. The volume flow rate-pressure difference relation is derived and compared to the prediction which utilizes the slip velocity concept. For very small channels, this concept is five orders of magnitude off, but reproduces asymptotically correct results for very large channels.     

Stress relaxation behavior of co-continuous PS/PMMA blends after step shear strain

Stress relaxation probing on the immiscible blends is an attractive route to reveal the time-dependent morphology–viscoelasticity correlations under/after flow. However, a comprehensive understanding on the stress relaxation of co-continuous blends, especially after subjected to a shear strain, is still lacking. In this work, the stress relaxation behavior of co-continuous polystyrene/poly(methyl methacrylate) (50/50) blends with different annealing times, strain levels, and temperatures was examined under step shear strain and was correlated with the development of their morphologies. It was found that co-continuous blends display a fast relaxation process which corresponded to the relaxation of bulk polymer and a second slower relaxation process due to the recovery of co-continuous morphology. The stress relaxation rates of co-continuous blends tend to decrease due to the coarsening of instable co-continuous structure during annealing. Furthermore, the stress relaxation of the co-continuous blends is strongly affected by the change of viscosity and interfacial tension caused by the temperature. The contribution of morphological coarsening, viscosity, and interfacial tension variation on the stress relaxation behavior of co-continuous blends was discussed based on the Lee–Park model and time–temperature superposition principle, respectively.     

Experimental validation of a dynamic resonant wall shear stress sensor

Experimental measurements are used to validate a numerical model of a dynamic resonant wall shear stress sensor. The numerical model consists of an unsteady two-dimensional boundary-layer model for the flow and a simple mechanical model for the sensor itself. The sensor’s sensitivity to wall shear stress is experimentally determined in a flat-plate boundary layer, and the results agree closely with those from the numerical simulations. Using the validated model, it is determined that the energy lost in each sensor oscillation due to the interaction between the sensor and fluid increases with increasing mean wall shear stress.     

Analytical investigation of the fluid flow in the interface region between a porous medium and a clear fluid in channels partially filled with a porous medium

In this paper analytical solutions for the steady fully developed laminar fluid flow in the parallel-plate and cylindrical channels partially filled with a porous medium and partially with a clear fluid are presented. The Brinkman-extended Darcy equation is utilized to model the flow in a porous region. The solutions account for the boundary effects and for the stress jump boundary condition at the interface recently suggested by Ochoa-Tapia and Whitaker. The dependence of the velocity on the Darcy number and on the adjustable coefficient in the stress jump boundary condition is investigated. It is shown that accounting for a jump in the shear stress at the interface essentially influences velocity profiles.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

Bingham fluid simulation with the incompressible lattice Boltzmann model

The Bingham fluid flow is numerically studied using the lattice Boltzmann method by incorporating the Papanastasiou exponential modification approach. The He–Luo incompressible lattice Boltzmann model is employed to avoid numerical instability usually encountered in non-Newtonian fluid simulations due to a strong non-linear relationship between the shear rate tensor and the rate-of-strain tensor. First, the value of the regularization parameter in Bingham fluid mimicking is analyzed and a method to determine the value is proposed. Then, the model is validated by pressure-driven planar channel flow and planar sudden expansion flow. The velocity profiles for the pressure-driven planar channel flow are in good agreement with analytical solutions. The calculated reattachment lengths for a 2:1 planar sudden expansion flow also agree well with the available data. Finally, the Bingham flow over a cavity is studied, and the streamlines and yielded/unyielded regions are discussed.     

Hydrodynamics study in a plane ultrafiltration module using an electrochemical method and particle image velocimetry visualization

 The wall shear stress is determined at the surface of a plane sheet of Plexiglas, taking the place of a membrane, using an electrochemical method. Several microelectrodes are mounted flush to this plane plate, and maps of shear stress are determined for two inlet and outlet configurations and three channel heights. The heterogeneity of the wall shear stress is observed for both configurations. Furthermore, the study of the turbulence features of the flow shows a decreasing fluctuating rate of velocity gradient when the channel height is decreased. The wall velocity gradients and turbulent intensity rates analysis are confirmed by a flow visualization using the particle image velocimetry method. Received: 25 September 2000 / Accepted: 23 April 2001     

Numerical modeling of turbulent compound channel flow using the lattice Boltzmann method

The flow of water in a straight compound channel with prismatic cross section is investigated with a relatively new tool, the lattice Boltzmann method. The large eddy simulation model is added in the lattice Boltzmann model for nonlinear shallow water equations (LABSWE^TM) so that the turbulence, caused by lateral exchange of momentum in the shear layer between the main channel and floodplain, can be taken into account and modeled efficiently. To validate the numerical model, a symmetrical compound channel with trapezoidal main channel and flat floodplain is tested. Similar to most natural watercourses, the floodplain has higher roughness values than the main channel. Different relative depths, D_r (the ratio of the depth of flow on the floodplain to that in the main channel), are considered. The Reynolds number is set at 30 000 in the main channel. The lateral distributions of the longitudinal velocity, the boundary shear stress, the Reynolds stress and the apparent shear stress across the channel are obtained after the large eddy simulation is performed. The results of numerical simulations are compared with the available experiment data, which show that the LABSWE^TM is capable of modeling the features of flow turbulence in compound channels and is sufficiently accurate for practical applications in engineering. Copyright © 2008 John Wiley & Sons, Ltd.