Rheology of poly(methyl methacrylate-<Emphasis Type="Italic">co</Emphasis>-styrene) particles suspended in water: effects of electrostatic surface layer

 Linear and nonlinear viscoelastic properties were examined for aqueous suspensions of monodisperse poly(methyl methacrylate-co-styrene) (MS) particles having the radius a ₀ =45 nm and the volume fractions φ=0.428−0.448. These particles had surface charges and the resulting electrostatic surface layer (electric double layer) had a thickness of t_s=5.7 nm. At low frequencies in the linear viscoelastic regime, the MS particles behaved approximately as the Brownian hard particles having an effective radius a _eff=a ₀ + t_s, and the dependence of their zero-shear viscosity η₀ on an effective volume fraction φ_eff (={a _eff/a ₀}³φ) agreed with the φ dependence of η₀ of ideal hard-core silica suspensions. In a range of φ_eff < 0.63, this φ_eff dependence was well described by the Brady theory. However, the φ_eff dependence of the high-frequency plateau modulus was weaker and the terminal relaxation mode distribution was narrower for the MS suspensions than for the hard-core suspensions. This result suggested that the electrostatic surface layer of the MS particles was soft and penetrable (at high frequencies). In fact, this “softness” was more clearly observed in the nonlinear regime: the nonlinear damping against step strain was weaker and the thinning under steady shear was less significant for the MS suspension than for the hard-core silica suspensions having the same φ_eff. These weaker nonlinearities of the concentrated MS particles with φ_eff∼ 0.63 (maximum volume fraction for random packing) suggested that the surface layers of those particles were mutually penetrating to provide the particles with a rather large mobility. Received: 10 July 2001 Accepted: 2 November 2001     

The interchain pressure effect in shear rheology

Recently, the tube diameter relaxation time in the evolution equation of the molecular stress function (MSF) model (Wagner et al., J Rheol 49: 1317–1327, 2005) with the interchain pressure effect (Marrucci and Ianniruberto, Macromolecules 37:3934–3942, 2004) included was shown to be equal to three times the Rouse time in the limit of small chain stretch. From this result, an advanced version of the MSF model was proposed, allowing modeling of the transient and steady-state elongational viscosity data of monodisperse polystyrene melts without using any nonlinear parameter, i.e., solely based on the linear viscoelastic characterization of the melts (Wagner and Rolón-Garrido 2009a, b). In this work, the same approach is extended to model experimental data in shear flow. The shear viscosity of two polybutadiene solutions (Ravindranath and Wang, J Rheol 52(3):681–695, 2008), of four styrene-butadiene random copolymer melts (Boukany et al., J Rheol 53(3):617–629, 2009), and of four polyisoprene melts (Auhl et al., J Rheol 52(3):801–835, 2008) as well as the shear viscosity and the first and second normal stress differences of a polystyrene melt (Schweizer et al., J Rheol 48(6):1345–1363, 2004), are analyzed. The capability of the MSF model with the interchain pressure effect included in the evolution equation of the chain stretch to model shear rheology on the basis of linear viscoelastic data alone is confirmed.     

Investigation of the different base fluid effects on the nanofluids heat transfer and pressure drop

A numerical study of laminar forced convective flows of three different nanofluids through a horizontal circular tube with a constant heat flux condition has been performed. The effect of Al₂O₃ volume concentration 0 ≤ φ ≤ 0.09 in the pure water, water-ethylene glycol mixture and pure ethylene glycol as base fluids, and Reynolds number of 100 ≤ Re ≤ 2,000 for different power inputs in the range of 10 ≤ Q(W) ≤ 400 have been investigated. In this study, all of the nanofluid properties are temperature and nanoparticle volume concentration dependent. The governing equations have been solved using finite volume approach with the SIMPLER algorithm. The results indicate an increase in the averaged heat transfer coefficient with increasing the mass of ethylene glycol in the water base fluid, solid concentration and Reynolds number. From the investigations it can be inferred that, the pressure drop and pumping power in the nanofluids at low solid volumetric concentration (φ < 3%) is approximately the same as in the pure base fluid in the various Reynolds numbers, but the higher solid nanoparticle volume concentration causes a penalty drop in the pressure. Moreover, this study shows it is possible to achieve a higher heat transfer rate with lower wall shear stress with the use of proper nanofluids.     

Nonlinear rheology of concentrated spherical silica suspensions: 3. Concentration dependence

Nonlinear rheology was examined for concentrated suspensions of spherical silica particles (with radius of 40 nm) in viscous media, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture and pure ethylene glycol. The particles were randomly and isotropically dispersed in the media in the quiescent state, and their effective volume fraction φ_eff ranged from 0.36 to 0.59. For small strains, the particles exhibited linear relaxation of the Brownian stress σ_B due to their diffusion. For large step strains γ, the nonlinear relaxation modulus G(t,γ) exhibited strong damping and obeyed the time-strain separability. This damping was related to γ-insensitivity of strain-induced anisotropy in the particle distribution that resulted in decreases of σ_B/γ. The damping became stronger for larger φ_eff. This φ_eff dependence was related to a hard-core volume effect, i.e., strain-induced collision of the particles that is enhanced for larger φ_eff. Under steady/transient shear flow, the particles exhibited thinning and thickening at low and high γ˙, respectively. The thinning behavior was well described by a BKZ constitutive equation using the G(t,γ) data and attributable to decreases of a Brownian contribution, σ_B/γ˙. The thickening behavior, not described by this equation, was related to dynamic clustering of the particles and corresponding enhancement of the hydrodynamic stress at high γ˙. In this thickening regime, the viscosity growth η₊ after start-up of flow was scaled with a strain γ˙t. Specifically, critical strains γ_d and γ_s for the onset of thickening and achievement of the steadily thickened state were independent of γ˙ but decreased with increasing φ_eff. This φ_eff dependence was again related to the hard-core volume effect, flow-induced collision of the particles enhanced for larger φ_eff. Received: 26 June 1998 Accepted: 9 December 1998     

Rheological analysis of highly concentrated w/o emulsions

A series of highly concentrated lipophilic cosmetic emulsions were analysed, in order to determine their rheological and textural properties, as a function of their microstructure. The originality of this study lies in the methodology used, especially the shear-stress scanning analysis. The results of a very powerful and comprehensive dynamic rheological analysis suggest the existence of two critical volume fraction values: besides the “close-packed” value φ _c , a “slack-packed” value φ₀, close to 0.60 could be demonstrated. It has been shown that the close-packed structure is stable under shear; in constrast, the slack-packed configuration, defined as φ₀<φ<φ _c is unstable under shear. A comparison with theoretical models, especially that of Princen, showed good agreement and allowed the close-packed value φ _c to be defined more precisely as 0.67. The gap between 0.67 and 0.74 is probably indicative of a highly polydisperse distribution, as confirmed by microscopic analysis. Flow experiments confirmed the validity of Princen‘s model. Received: 20 February 1997 Accepted: 20 January 1998     

Visualization of the breakdown of dilute particles gels during steady-shear

A feature of models concerning the rheology of coagulated suspensions is the development of shear planes. An experimental set-up was developed in which we investigated, by visual analysis, whether shear planes really develop in such systems during steady-shear. A transparent coagulated PTFE dispersion was used, in which the refractive indices of the continuous and dispersed phases were matched, for the formation of a gel. Coagulation was effected by adding NaCl to a concentration of 0.5 M. Polystyrene particles were built into the gel structure as tracer particles. During steady-shear the velocities and trajectories of the tracer particles were analyzed by Particle Tracking Velocimetry (PTV). Indeed layers with approximately the same velocity were observed for a coagulated PTFE-gel, during steady-shear. Deviations from rectilinear motion of the tracer particles were observed. These observations correspond with the assumptions of the giant floc model in which the shear is not distributed homogeneously, but limited to certain shear planes. Deviations from rectilinear motion of the tracer particles in a dilute gel correspond with the results found by Folkersma et al. (1998) in which the distance by which a moving particle entrains its neighbours was assumed to be larger at low volume fractions (≈0.1) than at high volume fractions (≈0.4). Received: 7 May 1998 Accepted: 31 August 1998     

Reological behaviour of concentrated monodisperse suspensions as a function of preshear conditions and temperature: an experimental study

 The influence of preshearing on the rheological behaviour of model suspensions was investigated with a stress-controlled cone-and-plate rheometer. The used matrix fluids showed Newtonian behaviour over the whole range of applied shear stresses. Highly monodisperse spherical glass spheres with various particle diameters were used as fillers. By applying steady preshearing at a low preshear stress, where a diffusion of particles can be expected, it was found for all model suspensions investigated at volume fractions ranging from 0.20 to 0.35 that the time-temperature superposition in the steady shear and in the dynamic mode holds within the chosen temperature range. Furthermore, all presheared model suspensions displayed a high and a low frequency range which are either separated by a shoulder or by a plateau value of G′ at intermediate frequencies. It could clearly be demonstrated that the low frequency range strongly depends on the preshear conditions. Hence, the features observed in the low frequency range can be attributed to a structure formation of a particulate network. In the high frequency range a frequency-dependent behaviour was observed which obeys the classical behaviour of Newtonian fluids (G′∝ω², G′′∝ω). The resulting temperature shift factors from the dynamic and the steady shear mode are identical and independent of the volume fraction and the particle size of the filler. Received: 29 November 2000 Accepted: 12 July 2001     

Shear-banding structure orientated in the vorticity direction observed for equimolar micellar solution

 The non-monotonic shear flow of a viscoelastic equimolar aqueous surfactant solution (cetylpyridinium chloride-sodium salicylate) is investigated rheologically and optically in a transparent strain-controlled Taylor Couette flow cell. As reported before, this particular wormlike micellar solution exhibits first a shear thinning and then a pronounced shear-thickening behavior. Once this shear-thickening regime is reached, a transient phase separation/shear banding of the solution into turbid and clear ring-like patterns orientated perpendicular to the vorticity axis, i.e., stacked like pancakes, is observed (Wheeler et al. 1998; Fischer 2000). The solution exhibit several unique features as no induction period of the shear induced phase, no structural build-up at the inner rotating cylinder, jumping pancake structure of clear and turbid ringlike phases, and oscillating shear stresses appear once the pancake structure is present. According to our analysis this flow phenomenon is not purely a mechanical or rheological driven hydrodynamic instability but one has to take into account structural changes of the oriented micellar aggregates (flow induced non-equilibrium phase transition) as proposed by several authors. Although this particular flow behavior and the underlying mixture of shear induced phases and mechanical instabilities is not fully understood yet, some classification characteristics based on a recent theoretical approach by Schmitt et al. (1995) and Porte et al. (1997) where a strong coupling between the flow instability (non-homogeneous flow profile due to the bands) and the structural changes causes the observed transient phenomena can be derived. In reference to the presented model the observed orientation of the rings is typical for complex fluids that undergo a spinodal phase separation coupled with a thermodynamic flow instability. In contrast to other shear banding phenomena, this one is observed in parallel plate, cone-plate, and Couette flow cell as well as under controlled stress and controlled rate conditions. Therefore, it adds an additional aspect to the present discussion on shear banding phenomena, i.e., the coupling of hydrodynamics and phase transition of rheological complex fluids. Received: 8 January 2001 Accepted: 15 May 2001     

Structure and orientation dynamics of sepiolite fibers–poly(ethylene oxide) aqueous suspensions under extensional and shear flow,probed by in situ SAXS

The structure and orientation dynamics of sepiolite clay fibers about 1,000 nm long and 10 nm thick, suspended in an aqueous poly(ehtylene oxide) matrix of 10⁵ g/mol molecular mass, have been studied under control extensional and shear flow. A new extensional flow cell developed at the “Laboratoire de Rhéologie” and the combined rheology and small angle X-ray scattering (Rheo-SAXS) setup available at the European Synchrotron Radiation Facility have allowed access to in situ and time-resolved fiber orientations and structure properties in the volume of suspensions under flow. In the volume fractions and shear rate domains for which the suspensions exhibit shear-thinning properties, two regimes of orientation separated by a critical strain rate have been identified under extensional flow.     

Non-isobaric Marangoni boundary layer flow for Cu,Al<Subscript>2</Subscript>O<Subscript>3</Subscript> and TiO<Subscript>2</Subscript> nanoparticles in a water based fluid

In this paper, a non-isobaric Marangoni boundary layer flow that can be formed along the interface of immiscible nanofluids in surface driven flows due to an imposed temperature gradient, is considered. The solution is determined using a similarity solution for both the momentum and energy equations and assuming developing boundary layer flow along the interface of the immiscible nanofluids. The resulting system of nonlinear ordinary differential equations is solved numerically using the shooting method along with the Runge-Kutta-Fehlberg method. Numerical results are obtained for the interface velocity, the surface temperature gradient as well as the velocity and temperature profiles for some values of the governing parameters, namely the nanoparticle volume fraction φ (0≤φ≤0.2) and the constant exponent β. Three different types of nanoparticles, namely Cu, Al₂O₃ and TiO₂ are considered by using water-based fluid with Prandtl number Pr =6.2. It was found that nanoparticles with low thermal conductivity, TiO₂, have better enhancement on heat transfer compared to Al₂O₃ and Cu. The results also indicate that dual solutions exist when β<0.5. The paper complements also the work by Golia and Viviani (Meccanica 21:200–204, 1986) concerning the dual solutions in the case of adverse pressure gradient.     

An Analytical Model for Capillary Pressure–Saturation Relation for Gas–Liquid System in a Packed-Bed of Spherical Particles

Capillary pressure is considered in packed-beds of spherical particles. In the case of gas–liquid flows in packed-bed reactors, capillary pressure gradients can have a significant influence on liquid distribution and, consequently, on the overall reactor performance. In particular, capillary pressure is important for non-uniform liquid distribution, causing liquid spreading as it flows down the packing. An analytical model for capillary pressure–saturation relation is developed for the pendular and funicular regions and the factors affecting capillary pressure in the capillary region are discussed. The present model is compared to the capillary pressure models of Grosser et al. (AIChE J., 34:1850–1860, 1988) and Attou and Ferschneider (Chem. Eng. Sci., 55:491–511, 2000) and to the experiments of Dodds and Srivastava (Part Part Syst. Charact., 23:29–39, 2006) and Dullien et al. (J. Colloid Interface Sci., 127:362–372, 1989). The non-homogeneity of real packings is considered through particle size and porosity distributions. The model is based on the assumption that the particles are covered with a liquid film, which provides hydrodynamic continuity. This makes the model more suitable for porous or rough particles than for non-porous smooth particles. The main improvements of the present model are found in the pendular region, where the liquid dispersion due to capillary pressure gradients is most significant. The model can be used to improve the hydrodynamic models (e.g., CFD and cellular automata models) for packed-bed reactors, such as trickle-bed reactors, where gas, liquid, and solid phases are present. Models for such reactors have become quite common lately (Sáez and Carbonell, AIChE J., 31:52–62, 1985; Holub et al., Chem. Eng. Sci, 47, 2343–2348, 1992; Attou et al., Chem. Eng. Sci., 54:785–802, 1999; Iliuta and Larachi, Chem. Eng. Sci., 54:5039–5045, 1999, IJCRE 3:R4, 2005; Narasimhan et al., AIChE J., 48:2459–2474, 2002), but they still lack proper terms causing liquid dispersion.     

An empirical model predicting the viscosity of highly concentrated, bimodal dispersions with colloidal interactions

The relationship between particle size distribution and viscosity of concentrated dispersions is of great industrial importance, since it is the key to get high solids dispersions or suspensions. The problem is treated here experimentally as well as theoretically for the special case of strongly interacting colloidal particles. An empirical model based on a generalized Quemada equation is used to describe η as a function of volume fraction for mono- as well as multimodal dispersions. The pre-factor η˜ accounts for the shear rate dependence of η and does not affect the shape of the η vs φ curves. It is shown here for the first time that colloidal interactions do not show up in the maximum packing parameter and φ_max can be calculated from the particle size distribution without further knowledge of the interactions among the suspended particles. On the other hand, the exponent ɛ is controlled by the interactions among the particles. Starting from a limiting value of 2 for non-interacting either colloidal or non-colloidal particles, ɛ generally increases strongly with decreasing particle size. For a given particle system it then can be expressed as a function of the number average particle diameter. As a consequence, the viscosity of bimodal dispersions varies not only with the size ratio of large to small particles, but also depends on the absolute particle size going through a minimum as the size ratio increases. Furthermore, the well-known viscosity minimum for bimodal dispersions with volumetric mixing ratios of around 30/70 of small to large particles is shown to vanish if colloidal interactions contribute significantly. Received: 7 June 2000/Accepted: 12 February 2001     

A constrained theory for single crystal shape memory wires with application to restrained recovery

The theory of thin wires developed in Dret and Meunier (Comptes Rendus de l’Académie des Sciences. Série I. Mathématique 337:143–147, 2003) is adapted to phase-transforming materials with large elastic moduli in the sense discussed in James and Rizzoni (J Elast 59:399–436, 2000). The result is a one-dimensional constitutive model for shape memory wires, characterized by a small number of material constants. The model is used to analyze self-accommodated and detwinned microstructures and to study superelasticity. It also turns out that the model successfully reproduces the behavior of shape memory wires in experiments of restrained recovery (Tsoi et al. in Mater Sci Eng A 368:299–310, 2004; Tsoi in 50:3535–3544, 2002; S̆ittner et al. in Mater Sci Eng A 286:298–311, 2000; vokoun in Smart Mater Struct 12:680–685, 2003; Zheng and Cui in Intermetallics 12:1305–1309, 2004; Zheng et al. in J Mater Sci Technol 20(4):390–394, 2004). In particular, the model is able to predict the shift to higher transformation temperatures on heating. The model also captures the effect of prestraining on the evolution of the recovery stress and of the martensite volume fraction.     

Derivation of the Kazemi–Gilman–Elsharkawy Generalized Dual Porosity Shape Factor

Kazemi et al. (SPE Reserv Eng 7(2):219–227, 1992) suggested an empirical matrix-fracture transfer function, verified based on experimental data of Mattax and Kyte (Trans AIME 225(15):177–184, 1962), to model fluid flow in naturally fractured dual porosity petroleum reservoirs using a dual-porosity numerical simulator. Their generalized shape factor should be valid for all possible irregular matrix blocks. The factor is calculated based on the volume of the matrix block, the surface open to flow in all directions and the distances of these surfaces to the centre of the matrix block. The summation is done over all open surfaces of a matrix block. Kazemi et al. (1992) showed that for rectangles and cylinders the formula reduces to the well-known forms of the shape factor. By the time, many authors indicated the validity of the formula, but no theoretical proof was offered for that so far. This study derives the Kazemi et al. (1992) shape factor using control volume finite difference discretization on the fracture-matrix dual continuum. The matrix blocks are handled as Voronoi polyhedra. The derivation is given for both isotropic and tensorial matrix permeability. Based on this derivation the authors conclude that the Kazemi et al. (SPE Reserv Eng 7(2):219–227, 1992) formula is exact under pseudo-steady-state conditions within the dual continuum mathematical concept of natural fractured dual porosity systems.     

Numerical solution of singular integral equations in stress concentration problems under longitudinal shear loading

Summary The paper deals with numerical solutions of singular integral equations in stress concentration problems for longitudinal shear loading. The body force method is used to formulate the problem as a system of singular integral equations with Cauchy-type singularities, where unknown functions are densities of body forces distributed in the longitudinal direction of an infinite body. First, four kinds of fundamental density functions are introduced to satisfy completely the boundary conditions for an elliptical boundary in the range 0≤φ _k ≤2π. To explain the idea of the fundamental densities, four kinds of equivalent auxiliary body force densities are defined in the range 0≤φ _k ≤π/2, and necessary conditions that the densities must satisfy are described. Then, four kinds of fundamental density functions are explained as sample functions to satisfy the necessary conditions. Next, the unknown functions of the body force densities are approximated by a linear combination of the fundamental density functions and weight functions, which are unknown. Calculations are carried out for several arrangements of elliptical holes. It is found that the present method yields rapidly converging numerical results. The body force densities and stress distributions along the boundaries are shown in figures to demonstrate the accuracy of the present solutions. Received 26 May 1998; accepted for publication 27 November 1998     

Rheological properties and dispersion stability of magnetorheological (MR) suspensions

In the present article, the rheological responses and dispersion stability of magnetorheological (MR) fluids were investigated experimentally. Suspensions of magnetite and carbonyl iron particles were prepared as model MR fluids. Under an external magnetic field (H ₀) and a steady shear flow, the yield stress depends upon H ₀ ^3/2. The Yield stress depended on the volume fraction of the particle (φ) linearly only at low concentration and increased faster at high fraction. Rheological behavior of MR fluids subjected to a small-strain oscillatory shear flow was investigated as a function of the strain amplitude, frequency, and the external magnetic field. In order to improve the stability of MR fluid, ferromagnetic Co-γ-Fe₂O₃ and CrO₂ particles were added as the stabilizing and thickening agent in the carbonyl iron suspension. Such needle-like particles seem to play a role in the steric repulsion between the relatively large carbonyl iron particles, resulting in improved stability against rapid sedimentation of dense iron particles. Furthermore, the additive-containing MR suspensions exhibited larger yield stress, especially at higher magnetic field strength. Received: 4 April 2000 Accepted: 6 November 2000     

Water in water emulsions: phase separation and rheology of biopolymer solutions

Partially miscible polymers in solution do not separate into two macroscopic phases; in general they behave as viscoelastic fluids containing droplets of the minority phase dispersed into a continuous majority phase (emulsion type systems). Both phases contain two types of polymers and solvent in variable amounts. With time, the smaller droplets tend to merge into larger ones and eventually sediment. Provided the time stability of the emulsion is long enough and the size of the droplets does not exceed a few tens of microns, the emulsion can be characterized by conventional rheological methods as an effective medium, both in the linear regime (viscoelasticity) and under flow. We investigated a ternary system composed by two biopolymers, a protein (caseinate) and a polysaccharide (alginate) in aqueous solution and established an analogy between these phase separated solutions and immiscible blends of polymers. We first characterized the biopolymers and determined the phase diagram at room temperature that we interpreted within the framework of the Edmond and Ogston model. For the rheological investigations, starting with an initial composition of the system, we separated the two phases by centrifugation. The individual phases were then characterized through their viscoelastic and flow behaviors. By recombining variable amounts of these phases, thereby varying only their volume fractions, we were able to prepare stable emulsions with constituents having constant compositions. The effect of shear on these emulsions was investigated. After different shearing protocols, the size of the droplets was derived from the Palierne model and the flow curves were analyzed. The droplet sizes were compared to the critical capillary numbers and coalescence predictions. The flow curves and the dynamic viscosities of the emulsions were interpreted with a model recently proposed by Kroy et al. that extends earlier work of Oldroyd (1953), Schowalter et al. (1968), and Frankel and Acrivos (1970). Received: 11 September 2000/Accepted: 21 December 2000     

Influence of inclination angle on buoyancy-driven convection in triangular enclosure filled with a fluid-saturated porous medium

The present investigation deals with the numerical analysis of steady-state laminar buoyancy-driven convection in an inclined triangular enclosure filled with fluid saturated porous media using the Darcy law equation. One wall of the enclosure is isothermally heated and the other is cooled, while the remaining wall is adiabatic. The effect of inclination angle on natural convection is investigated by varying the angle of inclination (φ) between 0° and 360°. The governing transformed equations are solved numerically using a finite-difference method. Obtained results are shown in the form of streamlines, isotherms, mean Nusselt numbers and dimensionless stream function for different values of the Rayleigh number Ra in the range 100 ≤ Ra ≤ 1,000. It is found that the values of the maximum and minimum mean Nusselt number are reached for φ = 330° and φ = 210° , respectively. However, the lowest flow strength is formed at φ = 240° for all values of Ra.     

The effect of mixing particles of different size on the electrorheological response under steady shear flow

 The effect of mixing particles of different sizes on the electrorheological response of suspensions under steady shear flow was investigated. Two sizes, 15 μm and 50 μm, of monodisperse spherical sulfonated poly(styrene-co-divinylbenzene) particles were used. Several electrorheological fluids were made containing different proportions of small and large particles dispersed in silicone oil, but with constant overall particulate concentration. It was found that the mixed size system produced the highest electrorheological response under the shear rates used (10 s⁻¹ to 500 s⁻¹), which is the opposite trend to previous studies of bimodal systems with larger size ratios. Received: 21 December 2000 Accepted: 29 March 2001     

The Giesekus fluid in ω–D form for steady two-dimensional flows

Numerical results of the simulation of the Giesekus model in ω–D form, which has previously been introduced in Part I of this study, are presented. The model has been applied to the flow of a concentrated polymer solution through a planar 3.97:1 contraction. To obtain an accurate fit of the rheological properties of the fluid a four-mode model is used. The predictions of the numerical simulations are directly compared with the experimental results published by Quinzani et al. in 1994. For the velocity fields a good quantitative agreement is reached, especially in the upstream channel. Regarding the shear stress and first normal stress difference, qualitative predictions of the experimental profiles are obtained.