On exact solutions of the flow equations for Bingham visco-plastic fluids through an eccentric annular cross-section

We discuss properties of solutions of the Bingham flow equations for visco-plastic fluids through an eccentric annular cross-section. Particularly, we perform arguments which are not in favor of the well-known Szabo–Hassager’s conjecture that the rigid zone is confined by circles provided the eccentricity is small (J Non-Newtonian Fluid Mech 45:149-169, 1992).     

In situ observations of unusual drop deformation and wobbling in simple shear flow

Deformation and wobbling of a liquid drop immersed in a liquid matrix were studied under mild shear conditions for various viscosity ratios. In situ visualization experiments were conducted on a homemade transparent Couette cell incorporated to the Paar Physica MCR500 shear rheometer. The effect of drop or matrix elasticity was examined and was found to play a major role in both deformation and wobbling processes. Experimental results were compared to Jackson and Tucker (J Rheol 47:659–682, 2003), Maffettone and Minale (J Non-Newton Fluid Mech 78:227–241, 1998) and Yu and Bousmina (J Rheol 47:1011–1039, 2003) ellipsoidal models. It was found that the agreement between the Newtonian models and the experimental results required an increase in the drop viscosity. Such increment in viscosity was found to scale with the first normal stress difference.     

An incompressible smoothed particle hydrodynamics scheme for Newtonian/non-Newtonian multiphase flows including semi-analytical solutions for two-phase inelastic Poiseuille flows

An incompressible smoothed particle hydrodynamics (ISPH) method is developed for the modeling of multiphase Newtonian and inelastic non-Newtonian flows at low density ratios. This new method is the multiphase extension of Xenakis et al, J. Non-Newtonian Fluid Mech., 218, 1-15, which has been shown to be stable and accurate, with a virtually noise-free pressure field for single-phase non-Newtonian flows. For the validation of the method a semi-analytical solution of a two-phase Newtonian/non-Newtonian (inelastic) Poiseuille flow is derived. The developed method is also compared with the benchmark multiphase case of the Rayleigh Taylor instability and a submarine landslide, thereby demonstrating capability in both Newtonian/Newtonian and Newtonian/non-Newtonian two-phase applications. Comparisons with analytical solutions, experimental and previously published results are conducted and show that the proposed methodology can accurately predict the free-surface and interface profiles of complex incompressible multi-phase flows at low-density ratios relevant, for example, to geophysical environmental applications.     

A differential model for the rheological properties of concentrated suspensions with weakly viscoelastic matrices

A model for the rheological properties of a concentrated suspension in weakly viscoelastic fluid matrices is proposed. The model is derived according to the Roscoe differential procedure described in 1952. The analytical results produced recently by Greco et al. (J Non-Newton Fluid Mech 147:1–10, 2007) and Housiadas and Tanner (J Non-Newton Fluid Mech 162:88–92, 2009) for dilute suspensions of neutrally buoyant, non-Brownian rigid spheres in weakly viscoelastic matrix fluids are the key results which are used as a base to predict the properties of concentrated suspensions. The results are compared with the few available experimental data from the literature, showing promising trends for the viscometric properties of the suspensions. In particular, one sees the rapidly increasing value of −N₂/N₁ as concentration increases.     

Stability for the 3D Navier–Stokes Equations with Nonzero far Field Velocity on Exterior Domains

Concerning to the non-stationary Navier–Stokes flow with a nonzero constant velocity at infinity, just a few results have been obtained, while most of the results are for the flow with the zero velocity at infinity. The temporal stability of stationary solutions for the Navier–Stokes flow with a nonzero constant velocity at infinity has been studied by Enomoto and Shibata (J Math Fluid Mech 7:339–367, 2005), in L ^p spaces for p ≥ 3. In this article, we first extend their result to the case \frac32 < p{\frac{3}{2} < p} by modifying the method in Bae and Jin (J Math Fluid Mech 10:423–433, 2008) that was used to obtain weighted estimates for the Navier–Stokes flow with the zero velocity at infinity. Then, by using our generalized temporal estimates we obtain the weighted stability of stationary solutions for the Navier–Stokes flow with a nonzero velocity at infinity.     

A method for calculating rheological and morphological properties of constant-volume polymer blend models in inhomogeneous shear fields

We show how to formulate two-point boundary-value problems in order to compute fully-developed laminar channel and tube flow profiles for viscoelastic fluid models. The formulation is applied to Couette and pressure-driven flows separately, or a combination of both. The application of this methodology is illustrated analytically for the Upper-Convected Maxwell Model, and it is applied computationally for the Phan-Thien/Tanner and Giesekus Models. Numerical solutions exist for the last two models [J.Y. Yoo, H.C. Choi, On the steady simple shear flows of the one-mode Giesekus fluid, Rheol. Acta 28 (1989) 13–24; P.J. Oliveira, F.T. Pinho, Analytical solution for fully developed channel and pipe flow of Phan-Thien–Tanner fluids, J. Fluid Mech. 387 (1999) 271–280; M.A. Alves, F.T. Pinho, P.J. Oliveira, Study of steady pipe and channel flows of a single-mode Phan-Thien–Tanner fluid, J. Non-Newtonian Fluid Mech. 101 (2001) 55–76], allowing verification of the computational technique. Subsequently, the computational algorithm is applied to the constant-volume polymer blend models of Maffettone and Minale [P.L. Maffettone, M. Minale, Equation of change for ellipsoidal drops in viscous flow, J. Non-Newtonian Fluid Mech. 84 (1999) 105–106 (Erratum), J. Non-Newtonian Fluid Mech. 78 (1998) 227–241] and Dressler and Edwards [M. Dressler, B.J. Edwards, The influence of matrix viscoelasticity on the rheology of polymer blends, Rheol. Acta 43 (2004) 257–282; M. Dressler, B.J. Edwards, Rheology of polymer blends with matrix-phase viscoelasticity and a narrow droplet size distribution, J. Non-Newtonian Fluid Mech. 120 (2004) 189–205]. Rheological and morphological properties of the model blends are thus obtained as functions of the spatial position within the channel, applied pressure drop, and shear rate at the wall.     

Group classification of one-dimensional nonisentropic equations of fluids with internal inertia

A systematic application of the group analysis method for modeling fluids with internal inertia is presented. The equations studied include models such as the nonlinear one-velocity model of a bubbly fluid (with incompressible liquid phase) at small volume concentration of gas bubbles (Iordanski Zhurnal Prikladnoj Mekhaniki i Tekhnitheskoj Fiziki 3, 102–111, 1960; Kogarko Dokl. AS USSR 137, 1331–1333, 1961; Wijngaarden J. Fluid Mech. 33, 465–474, 1968), and the dispersive shallow water model (Green and Naghdi J. Fluid Mech. 78, 237–246, 1976; Salmon 1988). These models are obtained for special types of the potential function W(r,[(r)\dot],S){W(\rho,\dot \rho,S)} (Gavrilyuk and Teshukov Continuum Mech. Thermodyn. 13, 365–382, 2001). The main feature of the present paper is the study of the potential functions with W _S  ≠ 0. The group classification separates these models into 73 different classes.     

Viscoelasticity in inkjet printing

We investigate the effects of viscoelasticity on drop generation in inkjet printing. In drop-on-demand printing, individual ink ‘drops’ are ejected from a nozzle by imposed pressure pulses. Upon exiting the nozzle, the shape of each ‘drop’ is that of a nearly spherical bead with a long thin trailing ligament. This ligament subsequently breaks up under the Rayleigh instability, typically into several small droplets (known as satellite drops). These satellite drops can create unwanted splash on the target substrate and a reduction in printing quality. Satellite drops can potentially be eliminated by adding polymer to the ink; elastic stresses can act to contract the trailing ligament into the main drop before capillary breakup occurs. However, elasticity can also reduce the drop velocity and can delay or even prevent the break-off of the drop from the ink reservoir within the nozzle. To achieve optimal drop shape and speed, non-Newtonian parameters such as the polymer concentration and molecular weight must be chosen correctly. We explore this parameter space via numerical simulations, using the Lagrangian–Eulerian finite-element method of Harlen et al. (J Non-Newtonian Fluid Mech 60:81–104, 1995). Results are compared with experimental observations taken from real printheads.     

Large-Eddy Simulation of Dispersion from Line Sources in a Turbulent Channel Flow

Large-eddy simulations of the dispersion from scalar line sources at various locations within a fully developed turbulent channel flow at Re = uh/ν = 10,400 are presented. Both mean and fluctuating scalar quantities are compared with those from the single available set of experimental data (Lavertu and Mydlarski, J Fluid Mech 528:135–172, 2005) and differences are highlighted and discussed. The results are also discussed in the context of scalar dispersion in other kinds of turbulent flows, e.g. homogeneous shear-flow. Initial computations at a much lower Reynolds number are also reported and compared with the two available direct numerical simulation data sets.     

Comparison of a quasi-Newtonian fluid with a viscoelastic fluid in planar contraction flow

The flow of a 5.0 wt.% solution of polyisobutylene in tetradecane through a planar 4 : 1 contraction exhibiting a shear thinning viscosity is simulated using the flow-type sensitive quasi-Newtonian fluid model. The shear viscosity is fitted by the Giesekus model, which, with the chosen parameters, leads to an extension thickening elongational viscosity. The stress and velocity fields of the numerical simulations are compared with the experimental results of Quinzani et al. [J. Non-Newtonian Fluid Mech. 52 (1994) 1–36] and the numerical results of the viscoelastic simulation using the Giesekus model of Azaiez et al. [J. Non-Newtonian Fluid Mech. 62 (1996) 253–277]. It can be shown that the quasi-Newtonian fluid qualitatively predicts the essential features of the flow in the vicinity of the contraction.     

Very high Reynolds number boundary layers over 3D sparse roughness and obstacles: the mean flow

Hot-wire and oil-film interferometry measurements are taken for 3D rough wall boundary layers at very high Reynolds numbers (61,000 < Re θ < 120,000) with low blockage ratios, 10 < δ/H < 135, and high roughness, 100 < H ⁺ < 4,900. The results cover flows over both rough walls and over obstacles and are compared with and provide extension to recent lower Reynolds number results. The validity of the Townsend ‘wall similarity hypothesis’ in relation to consistently increasing 3D roughness is interrogated. In agreement with recent work, Schultz and Flack (J Fluid Mech 580:381–405, 2007) and Castro (J Fluid Mech 585:469–485, 2007) found that, for relatively low roughness, Townsend’s hypothesis holds for the mean velocity field. With increasing roughness, the equilibrium layer diminishes and gradually vanishes. The viscous component of the wall shear stress decreases, while the turbulent component increases as the roughness effects extend across the boundary layer.     

Modelling elongational and shear rheology of two LDPE melts

Experimental data of two low-density polyethylene (LDPE) melts at 200°C for both shear flow (transient and steady shear viscosity as well as transient and steady first normal stress coefficient) and elongational flow (transient and steady-state elongational viscosity) as published by Pivokonsky et al. (J Non-Newtonian Fluid Mech 135:58–67, 2006) were analysed using the molecular stress function model for broadly distributed, randomly branched molecular structures. For quantitative modelling of melt rheology in both types of flow and in a very wide range of deformation rates, only three nonlinear viscoelastic material parameters are needed: Whilst the rotational parameter, a ₂, and the structural parameter, β, are found to be equal for the two melts considered, the melts differ in the parameter describing maximum stretch of the polymer chains.     

Nonequilibrium stretching dynamics of dilute and entangled linear polymers in extensional flow

We propose an extension of the FENE-CR model for dilute polymer solutions [M.D. Chilcott, J.M. Rallison, Creeping flow of dilute polymer solutions past cylinders and spheres, J. Non-Newtonian Fluid Mech. 29 (1988) 382–432] and the Rouse-CCR tube model for linear entangled polymers [A.E. Likhtman, R.S. Graham, Simple constitutive equation for linear polymer melts derived from molecular theory: Rolie–Poly equation, J. Non-Newtonian Fluid Mech. 114 (2003) 1–12], to describe the nonequilibrium stretching dynamics of polymer chains in strong extensional flows. The resulting models, designed to capture the progressive changes in the average internal structure (kinked state) of the polymer chain, include an ‘effective’ maximum contour length that depends on local flow dynamics. The rheological behavior of the modified models is compared with various results already published in the literature for entangled polystyrene solutions, and for the Kramers chain model (dilute polymer solutions). It is shown that the FENE-CR model with an ‘effective’ maximum contour length is able to describe correctly the hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation also observed and computed by Doyle et al. [P.S. Doyle, E.S.G. Shaqfeh, G.H. McKinley, S.H. Spiegelberg, Relaxation of dilute polymer solutions following extensional flow, J. Non-Newtonian Fluid Mech. 76 (1998) 79–110] and Li and Larson [L. Li, R.G. Larson, Excluded volume effects on the birefringence and stress of dilute polymer solutions in extensional flow, Rheol. Acta 39 (2000) 419–427] using Brownian dynamics simulations of bead–spring model. The Rolie–Poly model with an ‘effective’ maximum contour length exhibits a less pronounced hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation.     

Characterizing developing adverse pressure gradient flows subject to surface roughness

An experimental study was conducted to examine the effects of surface roughness and adverse pressure gradient (APG) on the development of a turbulent boundary layer. Hot-wire anemometry measurements were carried out using single and X-wire probes in all regions of a developing APG flow in an open return wind tunnel test section. The same experimental conditions (i.e., T _∞, U _ref, and C _p) were maintained for smooth, k ⁺ = 0, and rough, k ⁺ = 41–60, surfaces with Reynolds number based on momentum thickness, 3,000 < Re _θ < 40,000. The experiment was carefully designed such that the x-dependence in the flow field was known. Despite this fact, only a very small region of the boundary layer showed a balance of the various terms in the integrated boundary layer equation. The skin friction computed from this technique showed up to a 58% increase due to the surface roughness. Various equilibrium parameters were studied and the effect of roughness was investigated. The generated flow was not in equilibrium according to the Clauser (J Aero Sci 21:91–108, 1954) definition due to its developing nature. After a development region, the flow reached the equilibrium condition as defined by Castillo and George (2001), where Λ = const, is the pressure gradient parameter. Moreover, it was found that this equilibrium condition can be used to classify developing APG flows. Furthermore, the Zagarola and Smits (J Fluid Mech 373:33–79, 1998a) scaling of the mean velocity deficit, U _∞δ*/δ, can also be used as a criteria to classify developing APG flows which supports the equilibrium condition of Castillo and George (2001). With this information a ‘full APG region’ was defined.     

Use of particle image velocimetry to study heterogeneous drag reduction

Large polymer filaments can form when drag reducing polymers are injected through wall slots. The presence of these structures enhances the performance of the drag reducing function by mechanisms which are not understood. This paper shows how particle image velocimetry (PIV) techniques can be used to study changes in the configuration of the injected polymer and in the structure of the velocity field with increasing drag reduction. The filaments are found to behave as solid bodies which break up in high shear regions close to a boundary. The breakup process provides an explanation of why the filaments are not observed close to a wall and offers the possibility of providing a heterogeneous distribution of small aggregates of polymers which could be more effective than uniformly distributed molecules as suggested by Hoyer and Gyr (J Non-Newton Fluid Mech 65:221–240, 1996; J Fluids Eng 120:818–823, 1998), Dunlop and Cox (Phys Fluids 20:203–213, 1977) and Vlachogiannis et al. (Phys Fluid 15:3786–3794, 2004). PIV measurements show dramatic qualitative changes in the velocity patterns at maximum drag reduction.     

Stability of plane Poiseuille–Couette flows of a piezo-viscous fluid

We examine stability of fully developed isothermal unidirectional plane Poiseuille–Couette flows of an incompressible fluid whose viscosity depends linearly on the pressure as previously considered in Hron et al. [J. Hron, J. Málek, K.R. Rajagopal, Simple flows of fluids with pressure-dependent viscosities, Proc. R. Soc. Lond. A 457 (2001) 1603–1622] and Suslov and Tran [S.A. Suslov, T.D. Tran, Revisiting plane Couette–Poiseuille flows of a piezo-viscous fluid, J. Non-Newtonian Fluid Mech. 154 (2008) 170–178]. Stability results for a piezo-viscous fluid are compared with those for a Newtonian fluid with constant viscosity. We show that piezo-viscous effects generally lead to stabilisation of a primary flow when the applied pressure gradient is increased. We also show that the flow becomes less stable as the pressure and therefore the fluid viscosity decrease downstream. These features drastically distinguish flows of a piezo-viscous fluid from those of its constant-viscosity counterpart. At the same time the increase in the boundary velocity results in a flow stabilisation which is similar to that observed in Newtonian fluids with constant viscosity.     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 2: single-crystal plasticity

Variational formulations are constructed for rate-independent problems in small-deformation single-crystal strain-gradient plasticity. The framework, based on that of Gurtin (J Mech Phys Solids 50: 5–32, 2002), makes use of the flow rule expressed in terms of the dissipation function. Provision is made for energetic and dissipative microstresses. Both recoverable and non-recoverable defect energies are incorporated into the variational framework. The recoverable energies include those that depend smoothly on the slip gradients, the Burgers tensor, or on the dislocation densities (Gurtin et al. J Mech Phys Solids 55:1853–1878, 2007), as well as an energy proposed by Ohno and Okumura (J Mech Phys Solids 55:1879–1898, 2007), which leads to excellent agreement with experimental results, and which is positively homogeneous and therefore not differentiable at zero slip gradient. Furthermore, the variational formulation accommodates a non-recoverable energy due to Ohno et al. (Int J Mod Phys B 22:5937–5942, 2008), which is also positively homogeneous, and a function of the accumulated dislocation density. Conditions for the existence and uniqueness of solutions are established for the various examples of defect energy, with or without the presence of hardening or slip resistance.     

μ-PIV measurements of the ensemble flow fields surrounding a migrating semi-infinite bubble

Microscale particle image velocimetry measurements of ensemble flow fields surrounding a steadily migrating semi-infinite bubble through the novel adaptation of a computer controlled linear motor flow control system. The system was programmed to generate a square wave velocity input in order to produce accurate constant bubble propagation repeatedly and effectively through a fused glass capillary tube. We present a novel technique for re-positioning of the coordinate axis to the bubble tip frame of reference in each instantaneous field through the analysis of the sudden change of standard deviation of centerline velocity profiles across the bubble interface. Ensemble averages were then computed in this bubble tip frame of reference. Combined fluid systems of water/air, glycerol/air, and glycerol/Si-oil were used to investigate flows comparable to computational simulations described in Smith and Gaver III (J Fluid Mech 601:1–23, 2008) and to past experimental observations of interfacial shape. Fluorescent particle images were also analyzed to measure the residual film thickness trailing behind the bubble. The flow fields and film thickness agree very well with the computational simulations as well as existing experimental and analytical results. Particle accumulation and migration associated with the flow patterns near the bubble tip after long experimental durations are discussed as potential sources of error in the experimental method.     

Elongational viscosity of LDPE with various structures: employing a new evolution equation in MSF theory

Molecular stress function theory with new strain energy function is used to analyze transient extensional viscosity data of seven low-density polyethylene (LDPE) melts with various molecular structures as published by Stadler et al. (Rheol Acta 48:479–490, 2009) Pivokonsky et al. (J Non Newton Fluid Mech 135:58–67, 2006) and Wagner et al. (J Rheol 47(3):779–793, 2003). The new strain energy function has three nonlinear viscoelastic material parameters and assumes that the total stored energy of a branched molecule is given by different backbone and side chains stretching. The model parameters have been fitted for each LDPE in order to correlate with the supposed macromolecular structure expected from the type of synthesis. Most probable molecular structures for these LDPEs are comb and Cayley tree structures for respectively low- and high-molecular weight parts.     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.