On constitutive equations for electrorheological materials

Constitutive equations for electrorheological (ER) fluids have been based on experimental results for steady shearing flows and constant electric fields. The fluids have been modeled as being rigid until a yield stress is reached. Additional stress is then proportional to the shear rate. Recent experimental results indicate that ER materials have a regime of solid-like response when deformed from a rest state. They behave in a viscoelastic-like manner under sinusoidal shearing and exhibit time-dependent response under sudden changes in shear rate or electric field. In this work, a constitutive theory for ER materials is presented which accounts for these recent experimental observations. The stress is given by a functional of the deformation gradient history and the electric field vector. Using the methods of continuum mechanics, a general three-dimensional constitutive equation is obtained. A sample constitutive equation is introduced which is then used to determine the response of an ER material for different shear histories. The calculated shear response is shown to be qualitatively similar to that observed experimentally.     

Four-field Galerkin/least-squares formulation for viscoelastic fluids

A new Galerkin/Least-Squares (GLS) stabilized finite element method is presented for computing viscoelastic flows of complex fluids described by the conformation tensor; it extends the well-established GLS method for computing flows of incompressible Newtonian fluids. GLS methods are attractive for large-scale computations because they yield linear systems that can be solved easily with iterative solvers (e.g., the Generalized Minimum Residual method) and because they allow simple combinations of interpolation functions that can be conveniently and efficiently implemented on modern distributed-memory cache-based clusters.Like other state-of-the-art methods for computing viscoelastic flows (e.g., DEVSS-TG/SUPG), the new GLS method introduces a separate variable to represent the velocity gradient; with the aid of this variable, the conservation equations of mass, momentum, conformation, and the definition of velocity gradient are converted into a set of first-order partial differential equations in four unknown fields—pressure, velocity, conformation, and velocity gradient. The unknown fields are represented by low-order (continuous piecewise linear or bilinear) finite element basis functions.The method is applied to the Oldroyd-B constitutive equation and is tested in two benchmark problems—flow in a planar channel and flow past a cylinder in a channel. Results show that (1) the mesh-convergence rate of GLS is comparable to the DEVSS-TG/SUPG method; (2) the LS stabilization permits using equal-order basis functions for all fields; (3) GLS handles effectively the advective terms in the evolution equation of the conformation tensor; and (4) GLS yields accurate results at lower computational costs than DEVSS-type methods.     

Direct velocity gradient measurements using a dual-point laser—Doppler anemometer

The present paper describes a modified laser—Doppler anemometer set up to permit local velocity measurements at two separated points in space. In this way direct measurements of local shear and elongational rates become possible, providing information relevant in many flow studies of non-Newtonian fluids. Details of the optical and electronic sub-units are described and applications of the velocity gradient anemometer in channel flow are given to demonstrate its potential for measurement of rheological fluid properties and to provide insight into flows of non-Newtonian fluids.     

Laminar flow in three-dimensional square–square expansions

In this work we investigate the three-dimensional laminar flow of Newtonian and viscoelastic fluids through square–square expansions. The experimental results obtained in this simple geometry provide useful data for benchmarking purposes in complex three-dimensional flows. Visualizations of the flow patterns were performed using streak photography, the velocity field of the flow was measured in detail using particle image velocimetry and additionally, pressure drop measurements were carried out. The Newtonian fluid flow was investigated for the expansion ratios of 1:2.4, 1:4 and 1:8 and the experimental results were compared with numerical predictions. For all expansion ratios studied, a corner vortex is observed downstream of the expansion and an increase of the flow inertia leads to an enhancement of the vortex size. Good agreement is found between experimental and numerical results. The flow of the two non-Newtonian fluids was investigated experimentally for expansion ratios of 1:2.4, 1:4, 1:8 and 1:12, and compared with numerical simulations using the Oldroyd-B, FENE-MCR and sPTT constitutive equations. For both the Boger and shear-thinning viscoelastic fluids, a corner vortex appears downstream of the expansion, which decreases in size and strength when the elasticity of the flow is increased. For all fluids and expansion ratios studied, the recirculations that are formed downstream of the square–square expansion exhibit a three-dimensional structure evidenced by a helical flow, which is also predicted in the numerical simulations.     

A spectral/hp least‐squares finite element analysis of the Carreau–Yasuda fluids

A least‐squares finite element model with spectral/hp approximations was developed for steady, two‐dimensional flows of non‐Newtonian fluids obeying the Carreau–Yasuda constitutive model. The finite element model consists of velocity, pressure, and stress fields as independent variables (hence, called a mixed model). Least‐squares models offer an alternative variational setting to the conventional weak‐form Galerkin models for the Navier–Stokes equations, and no compatibility conditions on the approximation spaces used for the velocity, pressure, and stress fields are necessary when the polynomial order (p) used is sufficiently high (say, p > 3, as determined numerically). Also, the use of the spectral/hp elements in conjunction with the least‐squares formulation with high p alleviates various forms of locking, which often appear in low‐order least‐squares finite element models for incompressible viscous fluids, and accurate results can be obtained with exponential convergence. To verify and validate, benchmark problems of Kovasznay flow, backward‐facing step flow, and lid‐driven square cavity flow are used. Then the effect of different parameters of the Carreau–Yasuda constitutive model on the flow characteristics is studied parametrically. Copyright © 2016 John Wiley & Sons, Ltd.     

Two-dimensional viscoelastic flows: experimentation and modeling

Extensive experimental data on the birefringence in converging and diverging flows of a polymeric melt have been obtained. The birefringence and pressure drop measurements were carried out in working cells of planar geometry having different contraction angles and contraction ratios. For investigation of diverging or abrupt expansion flow, the direction of flow in the cells was reversed. The theoretical predictions are based upon the Leonov constitutive equation and a finite element scheme with streamwise integration.In contrast to Newtonian and second-order fluids, viscoelastic fluids at high shear rates show significant differences in pressure drop and birefringence (i.e. stresses) in converging and diverging flows. For a constant flow rate, the pressure drop is higher and the birefringence smaller in diverging flows than in converging flows. This difference increases with increasing flow rate. Further, for the same contraction ratio but different contraction angles, the birefringence maximum increases considerably with contraction angle. In addition, an increase in contraction ratio has the same effect.The viscoelastic constitutive equation of Leonov has been shown to describe all the above viscoelastic effects observed in the experiments. In general, a reasonable agreement between theory and experiment has been obtained, which shows the usefulness of the Leonov model in describing actual flows.     

GPU accelerated simulations of three-dimensional flow of power-law fluids in a driven cube

Newtonian fluid flow in two- and three-dimensional cavities with a moving wall has been studied extensively in a number of previous works. However, relatively a fewer number of studies have considered the motion of non-Newtonian fluids such as shear thinning and shear thickening power law fluids. In this paper, we have simulated the three-dimensional, non-Newtonian flow of a power law fluid in a cubic cavity driven by shear from the top wall. We have used an in-house developed fractional step code, implemented on a Graphics Processor Unit. Three Reynolds numbers have been studied with power law index set to 0.5, 1.0 and 1.5. The flow patterns, viscosity distributions and velocity profiles are presented for Reynolds numbers of 100, 400 and 1000. All three Reynolds numbers are found to yield steady state flows. Tabulated values of velocity are given for the nine cases studied, including the Newtonian cases.     

Flow characteristics in hydraulically equilibrium two-phase flows in a vertical 2 × 3 rod bundle channel

In order to increase data on two-phase flow distribution in a multi-subchannel system, being similar to a rod bundle, experiments have been carried out using water and air at ambient pressure and temperature as the working fluids and a newly constructed 2 × 3 rod bundle channel as the test channel. The channel contained six rods in rectangular array and two-kinds of six subchannels, simulating a BWR fuel rod bundle. Experimental data on flow distribution and pressure drop along each subchannel axis were obtained in various single- and two-phase flows under a hydraulic equilibrium flow condition. From the measured pressure drop in the single-phase flow, friction factor data in each subchannel were obtained. The two-phase pressure drop data were compared with calculations by a simple, one-dimensional, one-pressure two-fluid model. In addition, Taylor bubble velocity in each subchannel in slug-churn flows was measured with a double needle contact probe. Using the bubble velocity data, we obtained a subchannel void fraction in each subchannel, and discussed a relationship of the subchannel void fractions between two different subchannels. Results of such experiments and discussions are presented in this paper.     

Flow of wormlike micelle solutions through a periodic array of cylinders

Solutions of self-assembled wormlike micelles are used with ever increasing frequency in a multitude of consumer products ranging from cosmetic to industrial applications. Owing to the wide range of applications, flows of interest are often complex in nature; exhibiting both extensional and shear regions that can make modeling and prediction both challenging and valuable. Adding to the complexity, the micellar dynamics are continually changing, resulting in a number of interesting phenomena, such as shear banding and extensional flow instabilities. In this paper, we present the results of our investigation into the flow fields generated by a controllable and idealized porous media: a periodic array of cylinders. Our test channel geometry consists of six equally spaced cylinders, arranged perpendicular to the flow. By systematically varying the Deborah number, the flow kinematics, stability and pressure drop were measured. A combination of particle image velocimetry in conjunction with flush mount pressure transducers were used to characterize the flow, while flow induced birefringence measurements were used to determine micelle deformation and alignment. The pressure drop was found to decrease initially due to the shear thinning of the test fluid, and then exhibit a dramatic upturn as other elastic effects begin to dominate. We present evidence of the onset of an elastic instability in one of the test fluids above a critical Deborah number manifest in fluctuating transient pressure drop measurements and asymmetric streamlines. We argue that this disparity in the two test fluids can be attributed to the measurable differences in their extensional rheology.     

On the consequences of material frame-indifference in algebraic stress models

The principle of material frame-indifference (MFI) is a fundamental and controversial principle of continuum mechanics that has been invoked recently to derive nonlinear algebraic models for stresses of viscoelastic liquids. The purpose of the present study is to identify regions of a flow field where MFI should be considered. Such regions are identified by computing the angular velocity of the principal directions of the rate-of-deformation tensor in order to obtain an Euclidean objective vorticity tensor. An analysis is carried out for uniform shear and extensional flows, and for a Couette flow. The method is then applied to the planar flow through an abrupt 4:1 contraction and to the two-dimensional stream past a circular cylinder. The main results are: (1) MFI should be taken into account in regions characterized by the transition between two different kinematics and a significant velocity magnitude, and (2) MFI can be safely ignored in regions of pure viscometric behaviour as well as in recirculation regions. The consequences of MFI being taken into account are then examined upon using the Euclidean objective vorticity tensor in a simple algebraic constitutive law for viscoelastic fluids.     

Extension of Bernoulli's theorem on steady flows of inviscid fluids to steady flows of plastic solids

A differential equation is derived which is valid along any streamline in a steady flow of a general continuum where the streamline is also a trajectory of principal stress. For special materials and flow conditions this equation can be integrated to give algebraic relations between variables along the streamline. For inviscid fluids this leads to Bernoulli's famous theorem relating pressure and velocity along any streamline. For three-dimensional ideal flows of Tresca plastic solids and planar ideal flows of general rigid/perfectly plastic solids, it also leads to known results along any streamline. For other special constitutive materials in rigid/plastic solids additional streamline relations are obtained.     

A note on a derivative scheme for the finite volume method applied to incompressible viscous fluid

In the last two decades, there have been significant advances in application of the finite volume method to a wide spectrum of physical phenomena, ranging from heat transfer and compressible/incompressible fluid flow to solid mechanics. The finite volume method requires for different applications similar approximations of derivatives at the control surface. Within this framework, this study discusses a derivative scheme used in the finite volume method for incompressible viscous fluids. The numerical scheme is based on an implicit technique associated with the SIMPLE method to attain pressure–velocity coupling. The present work addresses simulation of two-dimensional flows in plane channels with and without contractions. The results show velocities and pressure fields with good agreement when compared to analytical results.     

A qualitative investigation of the equations of quasi-one-dimensional magnetohydrodynamic channel flow

A qualitative investigation of the system of differential equations describing the quasi-one-dimensional flow of an electrically conducting medium at small magnetic Reynolds numbers gives an idea of the different possible flow patterns occuring when the electromagnetic field and channel shape are given in different ways. Such a treatment is essential for the calculation of one-dimensional flows, and also for the solution of variational problems [1].In the literature devoted to this question studies have been made of flow in a one-dimensional electromagnetic field and a channel of constant cross section [2], as well as of the flow when the magnetic field is described by specially given functions of the flow velocity [3]. These cases reduce to the analysis of integral curves in a plane.In the present paper the investigation is carried out for an arbitrary distribution of the electric and magnetic fields and channel shape, which leads to a consideration of the behavior of integral curves in three-dimensional space. The qualitative results are illustrated by examples.     

Effect of viscoelasticity on the rotation of a sphere in shear flow

When particles are dispersed in viscoelastic rather than Newtonian media, the hydrodynamics will be changed entailing differences in suspension rheology. The disturbance velocity profiles and stress distributions around the particle will depend on the viscoelastic material functions. Even in inertialess flows, changes in particle rotation and migration will occur. The problem of the rotation of a single spherical particle in simple shear flow in viscoelastic fluids was recently studied to understand the effects of changes in the rheological properties with both numerical simulations [D’Avino et al., J. Rheol. 52 (2008) 1331–1346] and experiments [Snijkers et al., J. Rheol. 53 (2009) 459–480]. In the simulations, different constitutive models were used to demonstrate the effects of different rheological behavior. In the experiments, fluids with different constitutive properties were chosen. In both studies a slowing down of the rotation speed of the particles was found, when compared to the Newtonian case, as elasticity increases. Surprisingly, the extent of the slowing down of the rotation rate did not depend strongly on the details of the fluid rheology, but primarily on the Weissenberg number defined as the ratio between the first normal stress difference and the shear stress.In the present work, a quantitative comparison between the experimental measurements and novel simulation results is made by considering more realistic constitutive equations as compared to the model fluids used in previous numerical simulations [D’Avino et al., J. Rheol. 52 (2008) 1331–1346]. A multimode Giesekus model with Newtonian solvent as constitutive equation is fitted to the experimentally obtained linear and nonlinear fluid properties and used to simulate the rotation of a torque-free sphere in a range of Weissenberg numbers similar to those in the experiments. A good agreement between the experimental and numerical results is obtained. The local torque and pressure distributions on the particle surface calculated by simulations are shown.     

Modification of near-wall turbulence structure in a shear-driven three-dimensional turbulent boundary layer

 Most high Reynolds number flows of engineering interest are three-dimensional in nature. Key features of three-dimensional turbulent boundary layers (3DTBLs) include: non-colateral shear stress and strain rate vectors, and decreasing ratio of the shear stresses to the turbulent kinetic energy with increasing three-dimensionality. These are indicators that the skewing has a significant effect on the structure of turbulence. In order to further investigate the flow physics and turbulence structure of these complex flows, an innovative method for generating a planar shear-driven 3DTBL was developed. A specialized facility incorporating a relatively simple geometry and allowing for varying strengths of crossflow was constructed to facilitate studies where the skewing is decoupled from the confounding effects of streamwise pressure gradient and curvature. On-line planar particle image velocimetry (PIV) measurements and flow visualization results indicate that the experimental configuration generates the desired complex flow, which exhibits typical characteristics associated with 3DTBLs. Furthermore, spanwise shear results in modification of the near-wall turbulence structure. Analysis of near-wall flow visualization photographs revealed a reduction of mean streak length with increasing spanwise shear, while streak spacing remained relatively constant. In the most strongly sheared case, where the belt velocity is twice that of the freestream velocity, the mean streak length was reduced by approximately 50%. Received: 28 October 1997/Accepted: 4 February 1998     

Three dimensional finite element analysis of the flow of polymer melts

The finite element simulation of a selection of two- and three-dimensional flow problems is presented, based upon the use of four different constitutive models for polymer melts (Oldroyd-B, Rolie-Poly, Pom-Pom and XPP). The mathematical and computational models are first introduced, before their application to a range of visco-elastic flows is described. Results demonstrate that the finite element models used here are able to re-produce predictions made by other published numerical simulations and, significantly, by carefully conducted physical experiments using a commercial-grade polystyrene melt in a three-dimensional contraction geometry. The paper also presents a systematic comparison and evaluation of the differences between two- and three-dimensional simulations of two different flow regimes: flow of an Oldroyd-B fluid around a cylinder and flow of a Rolie-Poly fluid into the contraction geometry. This comparison allows new observations to be made concerning the relatively poor quality of two-dimensional simulations for flows in even quite deep channels.     

Pure axial flow of viscoelastic fluids in rectangular microchannels under combined effects of electro-osmosis and hydrodynamics

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien–Tanner (PTT) model with the Gordon–Schowalter convected derivative which covers the upper convected Maxwell, Johnson–Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson–Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid system and is verified against 1D analytical solution of the velocity profile with less than 0.06% relative error. Also, a parametric study is carried out to investigate the effect of channel aspect ratio (width to height), wall zeta potential and the Debye–Hückel parameter on 2D velocity profile, volumetric flow rate and the Poiseuille number in the mixed EO/PD flows of viscoelastic fluids with different Weissenberg numbers. Our results show that, for low channel aspect ratios, the previous 1D analytical models underestimate the velocity profile at the channel half-width centerline in the case of favorable pressure gradients and overestimate it in the case of adverse pressure gradients. The results reveal that the inapplicability of the Debye–Hückel approximation at high zeta potentials is more significant for higher Weissenberg number fluids. Also, it is found that, under the specified values of electrokinetic parameters, there is a threshold for velocity scale ratio in which the Poiseuille number is approximately independent of channel aspect ratio.     

Particle imaging velocimetry of a wet aqueous foam with an underlying liquid film

 We investigate the utility of particle imaging velocimetry (PIV) for performing kinematic measurements in wet aqueous foam with a liquid film beneath it. The flow velocities are measured near the walls of a square cross-section horizontal duct. The flow velocities are useful for validating the rheological models. We show that there is a discrepancy between the velocity profiles in the wet foam and the Bingham plastic model of flow. The velocity measurements reveal a more complex flow pattern, which may be analysed following three different regimes: a plug flow, a shear flow in a vertical plane and a three-dimensional shear flow. The transition between the plug flow and the shear flows may be explained by a shear-induced migration of bubbles. Received: 25 April 2000 / Accepted: 26 February 2001     

Velocity-profile deformation in EHD flow in a two-dimensional channel

The study of unipolar-charged fluids in the presence of external and induced electric fields has recently taken on great importance. The characteristics of one-dimensional EGD flows [1, 2] and developed laminar flows of a viscous fluid [3] have been clarified in several studies made in this field. However, the study of three-dimensional flows of such media is actually just beginning. Here, along with the analysis of three-dimensional boundary layers and jets [4], there is considerable interest in the study of spatial (two-dimensional and three-dimensional) EHD flows of an inviscid fluid, since in many engineering devices the zone of interaction of the flow with the electric fields does not exceed a few channel diameters, which makes it possible to neglect viscous effects.In this paper we examine some aspects of two-dimensional EHD flows of a viscous incompressible medium for infinitely large electric Reynolds numbers. The perturbations of the hydrodynamic parameters of the flow downstream from the zone of action of the electrostatic forces are determined. It is shown that in many cases the flow parameters outside this zone may be determined without solving the complete system of EHD partial differential equations.