An experimental and numerical investigation of the dynamics of microconfined droplets in systems with one viscoelastic phase

The dynamics of single droplets in a bounded shear flow is experimentally and numerically investigated for blends that contain one viscoelastic component. Results are presented for systems with a viscosity ratio of 1.5 and a Deborah number for the viscoelastic phase of 1. The numerical algorithm is a volume-of-fluid method for tracking the placement of the two liquids. First, we demonstrate the validation of the code with an existing boundary integral method and with experimental data for confined systems containing Newtonian components. This is followed by numerical simulations and experimental data for the combined effect of geometrical confinement and component viscoelasticity on the droplet dynamics after startup of shear flow at a moderate capillary number. The viscoelastic liquids are Boger fluids, which are modeled with the Oldroyd-B constitutive model and the Giesekus model. Confinement substantially increases the viscoelastic stresses and the elongation rates in and around the droplet. We show that the latter can be dramatic for the use of the Oldroyd-B model in confined systems with viscoelastic components. A sensitivity analysis for the choice of the model parameters in the Giesekus constitutive equation is presented.     

Development and implementation of VOF-PROST for 3D viscoelastic liquid–liquid simulations

We implement a volume-of-fluid algorithm with a parabolic re-construction of the interface for the calculation of the surface tension force (VOF-PROST). This achieves higher accuracy for drop deformation simulations in comparison with existing VOF methods based on a piecewise linear interface re-construction. The algorithm is formulated for the Giesekus constitutive law. The evolution of a drop suspended in a second liquid and undergoing simple shear is simulated. Numerical results are first checked against two cases in the literature: the small deformation theory for second-order liquids, and an Oldroyd-B extensional flow simulation. We then address the experimental data of Guido et al. (2003) for a Newtonian drop in a viscoelastic matrix liquid. The data deviate from existing theories as the capillary number increases, and reasons for this are explored here with the Oldroyd-B and Giesekus models.     

Forward roll coating flows of viscoelastic liquids

Roll coating is distinguished by the use of one or more gaps between rotating cylinders to meter and apply a liquid layer to a substrate. Except at low speed, the two-dimensional film splitting flow that occurs in forward roll coating is unstable; a three-dimensional steady flow sets in, resulting in more or less regular stripes in the machine direction. For Newtonian liquids, the stability of the two-dimensional flow is determined by the competition of capillary and viscous forces: the onset of meniscus nonuniformity is marked by a critical value of the capillary number. Although most of the liquids coated industrially are non-Newtonian polymeric solutions and dispersions, most of the theoretical analyses of film splitting flows relied on the Newtonian model. Non-Newtonian behavior can drastically change the nature of the flow near the free surface; when minute amounts of flexible polymer are present, the onset of the three-dimensional instability occurs at much lower speeds than in the Newtonian case.Forward roll coating flow is analyzed here with two differential constitutive models, the Oldroyd-B and the FENE-P equations. The results show that the elastic stresses change the flow near the film splitting meniscus by reducing and eventually eliminating the recirculation present at low capillary number. When the recirculation disappears, the difference of the tangential and normal stresses (i.e., the hoop stress) at the free surface becomes positive and grows dramatically with fluid elasticity, which explains how viscoelasticity destabilizes the flow in terms of the analysis of Graham [M.D. Graham, Interfacial hoop stress and instability of viscoelastic free surface flows, Phys. Fluids 15 (2003) 1702–1710].     

Dispersion of linear gravity waves on a viscoelastic fluid in an horizontal canal

A characteristic equation is derived that describes the spatial decay of linear surface gravity waves on Maxwell fluids. Except at small frequencies, the derived dispersion relation is different from the temporal decay dispersion relation which is normally studied within fluid mechanics. The implications for waves on viscous Newtonian fluids using the two different dispersion relations is briefly discussed. The wave number is measured experimentally as function of the frequency in a horizontal canal. Seven Newtonian fluids and four viscoelastic liquids with constant viscosity have been used in the experiments. The spatial decay theory for Newtonian fluids fits well to the experimental data. The model and experiments are used to determine limits for the Maxwell fluid time numbers for the four viscoelastic liquids. As a result of low viscosity it was not possible within this study to obtain these time numbers from oscillatory experiments. Therefore, a comparison of surface gravity wave experiments with theory is applicable as a method to evaluate memory times of low viscosity viscoelastic fluids.     

Selective withdrawal of polymer solutions: Computations

This paper reports numerical simulations of selective withdrawal of Newtonian and polymeric liquids, and complements the experimental study reported in the accompanying paper (Zhou and Feng [2]). We use finite elements to solve the Navier–Stokes and constitutive equations in the liquid on an adaptively refined unstructured grid, with an arbitrary Lagrangian–Eulerian scheme to track its free surface. The rheology of the viscoelastic liquids are modeled by the Oldroyd-B and Giesekus equations, and the physical and geometric parameters are matched with those in the experiments. The computed interfacial deformation is in general agreement with the experimental observations. In particular, the critical condition for interfacial rupture is predicted to quantitative accuracy. Furthermore, we combine the numerical and experimental data to explore the potential of selective withdrawal as an extensional rheometer. For Newtonian fluids, the measured steady elongational viscosity is within 47% of the actual value, apparently with better accuracy than other methods applicable to low-viscosity liquids. For polymer solutions, an estimated maximum error of 300% compares favorably with prior measurements.     

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.     

Effects of viscoelasticity on the retraction of a sheared drop

Effects of drop and matrix viscoelasticity on the retraction of a sheared drop are numerically investigated. Retraction of an Oldroyd-B drop in a Newtonian matrix is initially faster and later slower with increasing drop Deborah number. The observed behavior is explained using an ordinary differential equation model representing the dominant balance between various forces during retraction. The initial faster relaxation of viscoelastic drops is due to viscoelastic stresses pulling the drop interface at the tips inward. The later slower retraction is due to the slowly-relaxing viscoelastic forces at the equator, where they act against the capillary force. The drop inclination decreases substantially during retraction unlike in a Newtonian case. Matrix viscoelasticity slows the relaxation of a Newtonian drop because of the increasingly slow relaxation of highly stretched polymers near the drop tip with increasing Deborah number. Increasing the ratio of polymeric to total viscosity further accentuates the viscoelastic effects in both cases. For an Oldroyd-B drop in an Oldroyd-B matrix, a competition between the dispersed and the continuous phase elasticities, represented by their ratio, determines the dynamics; larger values of the ratio leads again to initial faster and later slower retraction.     

Two-dimensional study of drop deformation under simple shear for Oldroyd-B liquids

The effect of viscoelasticity on the deformation of a circular drop suspended in a second liquid in shear is investigated with direct numerical simulations. A numerical algorithm based on the volume-of-fluid method for interface tracking is implemented in two dimensions with the Oldroyd-B constitutive model for viscoelastic liquids. The code is verified against a normal mode analysis for the stability of two-layer flow in a channel; theoretical growth rates are reproduced for the interface height, velocity and stress components. Drop simulations are performed for drop and matrix liquids of different viscosities and elasticities. A new feature is found for the case of equal viscosity, when the matrix liquid is highly elastic and surface tension is low; hook-like structures form at the drop tips. This is due to the growth of first normal stress differences that occur slightly above the front tip and below the back tip as the matrix elasticity increases above a threshold value.     

纤维悬浮液搅拌流动的数值模拟

由于缺乏适当的本构方程，对纤维悬浮液流动的研究一直局限于纤维的牛顿流体悬浮液。本文采用ＭＵＣＭ模型对作者最近提出的纤维Ｏｌｄｒｏｙｄ－Ｂ流体悬浮液的本构方程作了改进，并对锚式桨搅拌槽的二维Ｏｌｄｒｏｙｄ－Ｂ流体和牛顿流体纤维悬浮液搅拌流动作了数值模拟。模拟的结果表明，本文所用的模型和方法能有效地抑制过大局部应力的影响并合理地处理流体的记忆效应。     

Anomalous rolling of spheres down an inclined plane

A sphere in air will roll down a plane that is tilted away from the vertical. The only couple acting about the point of contact between the sphere and the plane is due to the component of the weight of the sphere along the plane, provided that air friction is negligible. If on the other hand the sphere is immersed in a liquid, hydrodynamic forces will enter into the couples that turn the sphere, and the rotation of the sphere can be anomalous, i.e., as if rolling up the plane while it falls. In this paper we shall show that anomalous rolling is a characteristic phenomenon that can be observed in every viscoelastic liquid tested so far. Anomalous rolling is normal for hydrodynamically levitated spheres, both in Newtonian and viscoelastic liquids. Normal and anomalous rolling are different names for dry and hydrodynamic rolling. Spheres dropped at a vertical wall in Newtonian liquids are forced into anomalous rotation and are pushed away from the wall while in viscoelastic liquids, they are forced into anomalous rotation, but are pushed toward the wall. If the wall is inclined and the fluid is Newtonian, the spheres will rotate normally for dry rolling, but the same spheres rotate anomalously in viscoelastic liquids when the angle of inclination from the vertical is less than some critical value. The hydrodynamic mechanisms underway in the settling of circular particles in a Newtonian fluid at a vertical wall are revealed by an exact numerical simulation based on a finite-element solution of the Navier-Stokes equations and Newton's equations of motion for a rigid body.     

Effect of time-periodic vertical oscillations of the Rayleigh–Bénard system on nonlinear convection in viscoelastic liquids

A study of heat transport in Rayleigh–Bénard convection in viscoelastic liquids with/without gravity modulation is made using a most minimal representation of Fourier series and a representation with higher modes. The Oldroyd-B constitutive relation is considered. The resulting non-autonomous Lorenz model (generalized Khayat–Lorenz model of four modes and seven modes) is solved numerically using the adaptive-grid Runge–Kutta–Fehlberg45 method to quantify the heat transport. The effect of gravity modulation is shown to be stabilizing there by leading to a situation of reduced heat transfer. The Deborah number is shown to have an antagonistic influence on convection compared to the stabilizing effect of modulation amplitude and elastic ratio. The results in respect of Maxwell, Rivlin–Ericksen and Newtonian liquids are obtained as particular cases of the present study. A transformation of the momentum equations illustrates the equivalence of present approach and the one due to Khayat that uses normal stresses explicitly.     

The influence of Oldroyd-B and PTT lubricants on moving journal bearing systems

The influence of viscoelasticity on the performance of statically and dynamically loaded journal bearings is considered. The lubricant in the system is modelled using either the Oldroyd-B or linear PTT models. Significant viscoelastic effects are presented for both moderate and narrow gap journal bearing configurations. The dynamical behaviour of the journal bearing system is shown to be dependent on the fluid model, the relaxation time and also the gap size.     

Squeezing a viscoelastic liquid from a wedge: an exact solution

The paper reports an exact solution for the squeezing flow from a wedge of a general viscoelastic liquid. To obtain numerical values for the field variables, a network model that allows stress overshoot and shear-thinning in the start-up of a shear flow is adopted. It is found that both these features are important in this transient flow; stress overshoot is responsible for a stiffer response of the fluid (compared to the inelastic case) at moderate time —at large time, shear-thinning dominates and the fluid behaves like an inelastic fluid. On the other hand, the Oldroyd-B fluid always predicts a softer response than the Newtonian one. Furthermore, there is a limiting Weissenberg number above which one component of the stresses of the Oldroyd-B fluid increases unboundedly with time. This limiting Weissenberg number is approximately $\text{sol2}\text{3}$.     

Viscoelastic effects in double‐pipe single‐pass counterflow heat exchangers

The conducto‐convective heat loss from a viscoelastic liquid, in the core of a double‐pipe heat exchanger arrangement, to a cooler Newtonian fluid flowing in the outer annulus is investigated with direct numerical simulations. A numerical algorithm based on the finite difference method is implemented in time and space with the Giesekus constitutive model for the viscoelastic liquids. The flow of both the annulus and core‐fluids is considered to be Poiseuille flow, driven by respective pressure gradients. In general, the results show that a viscoelastic core‐fluid leads to slightly lower (albeit comparable) attainable temperatures in the core‐fluid stream as compared with a corresponding Newtonian fluid. Copyright © 2008 John Wiley & Sons, Ltd.     

Surfactant spreading on a thin weakly viscoelastic film

A mathematical model is presented for surfactant-driven thin weakly viscoelastic film flows on a flat, impermeable plane. The Oldroyd-B constitutive relation is used to model the viscoelastic fluid. Lubrication theory and a perturbation expansion in powers of the Weissenberg number (We) are employed, which give rise to non-linear coupled evolution equations governing the transport of insoluble surfactant and thin liquid film thickness. Spreading on a Newtonian film is recovered to leading order and corrections to viscoelasticity are obtained at order We. These equations are solved numerically over a wide range of viscosity ratio (ratio of solvent viscosity to the sum of solvent and polymeric viscosities), pre-existing surfactant level and Peclet number (Pe). The effect of viscoelasticity on surfactant transport and fluid flow is investigated and the mechanisms underlying this effect are explored. Shear stress, streamwise normal stress and the temporal rate of change of extra shear stress generated from gradients in surfactant concentration dominate thin viscoelastic film flows whereas only shear stresses play a role in Newtonian thin film flows. Our results also reveal that, for weak viscoelasticity, the influence of viscosity ratio on the evolution of surfactant concentration and film thickness can be significant and varies considerably, depending on the concentration of pre-existing surfactant and surfactant surface diffusivity.     

Uniform flow of viscoelastic fluids past a confined falling cylinder

Uniform steady flow of viscoelastic fluids past a cylinder placed between two moving parallel plates is investigated numerically with a finite-volume method. This configuration is equivalent to the steady settling of a cylinder in a viscoelastic fluid, and here, a 50% blockage ratio is considered. Five constitutive models are employed (UCM, Oldroyd-B, FENE-CR, PTT and Giesekus) to assess the effect of rheological properties on the flow kinematics and wake patterns. Simulations were carried out under creeping flow conditions, using very fine meshes, especially in the wake of the cylinder where large normal stresses are observed at high Deborah numbers. Some of the results are compared with numerical data from the literature, mainly in terms of a drag coefficient, and significant discrepancies are found, especially for the constant-viscosity constitutive models. Accurate solutions could be obtained up to maximum Deborah numbers clearly in excess of those reported in the literature, especially with the PTT and FENE-CR models. The existence or not of a negative wake is identified for each set of model parameters.     

Drop oscillations under simple shear in a highly viscoelastic matrix

Recent two-dimensional numerical simulations and experiments have shown that, when a drop undergoes shear in a viscoelastic matrix liquid, the deformation can undergo an overshoot. I implement a volume-of-fluid algorithm with a paraboloid reconstruction of the interface for the calculation of the surface tension force for three-dimensional direct numerical simulations for a Newtonian drop in an Oldroyd-B liquid near criticalities. Weissenberg numbers up to 1 at viscosity ratio 1 and retardation parameter 0.5 are examined. Critical capillary numbers rise with the Weissenberg number. Just below criticality, drop deformation begins to undergo an overshoot when the Weissenberg number is sufficiently high. The overshoot becomes more pronounced, and at higher matrix Weissenberg numbers, such as 0.8, drop deformation undergoes novel oscillations before settling to a stationary shape. Breakup simulations are also described.     

Stokes＇ first problem for a viscoelastic fluid with the generalized Oldroyd-B model

The flow near a wall suddenly set in motion for a viscoelastic fluid with the generalized Oldroyd-B model is studied. The fractional calculus approach is used in the constitutive relationship of fluid model. Exact analytical solutions of velocity and stress are obtained by using the discrete Laplace transform of the sequential fractional derivative and the Fox H-function. The obtained results indicate that some well known solutions for the Newtonian fluid, the generalized second grade fluid as well as the ordinary Oldroyd-B fluid, as limiting cases, are included in our solutions. The project supported by the National Natural Science Foundation of China (10272067), the Doctoral Program Foundation of the Education Ministry of China (20030422046), the Natural Science Foundation of Shandong Province, China (Y2006A14) and the Research Foundation of Shandong University at Weihai. The English text was polished by Keren Wang.     

Stability of small radial oscillations of a spherical layer of liquid

An isothermal spherical layer of a viscoelastic liquid described by the one-parameter Maxwell model is considered. When the model parameter is taken equal to zero, a model of a purely viscous Newtonian fluid is obtained. The stability of the spherical layer of liquid with respect to small radial perturbations of the velocity and pressure is investigated for both types of liquids. Leningrad. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 170–171, September–October, 1988.     

Phosphate Removal Through Nano-Zero-Valent Iron Permeable Reactive Barrier; Column Experiment and Reactive Solute Transport Modeling

The stability of natural penetrative convection arising due to a uniform internal heat source in a vertical porous layer saturated with an Oldroyd-B fluid is investigated. The vertical walls of the porous layer are impermeable and maintained at different uniform temperatures. The energy stability analysis performed reveals that the system is unconditionally stable even in the presence of internal heating in the case of Newtonian fluids, while for viscoelastic fluids the base flow is found to be unstable. As the energy stability analysis of Gill type is unable to decide the stability of the system, the Galerkin method is used to solve the complex eigenvalue problem. The internal heating introduces asymmetry in the basic flow and amounts to the existence of different set of onset modes. The internal heating and stress relaxation parameter facilitates instability of the system while increasing strain retardation parameter discloses stabilizing effect on the system. Moreover, the critical Darcy–Rayleigh number, wave number and wave speed become invariant as Ns becomes large. The streamlines and isotherms presented herein demonstrate the development of complex dynamics at the critical state.