A 2D boundary element method for simulating the deformation of axisymmetric compound non‐Newtonian drops

The boundary integral formulation of the solution to the Stokes equations is used to describe the deformation of small compound non‐Newtonian axisymmetric drops suspended in a Newtonian fluid that is subjected to an axisymmetric flow field. The non‐Newtonian stress is treated as a source term in the Stokes equations, which yields an extra integral over the domains containing non‐Newtonian material. By transforming the integral representation for the velocity to cylindrical co‐ordinates and performing the integration over the azimuthal direction analytically, the dimension of the problem can be reduced from three to two. A boundary element method for the remaining two‐dimensional problem aimed at the simulation of the deformation of such axisymmetric compound non‐Newtonian drops is developed. Apart from a numerical validation of the method, simulation results for a drop consisting of an Oldroyd‐B fluid and a viscoelastic material are presented. Moreover, the method is extended to compound drops that are composed of a viscous inner core encapsulated by a viscoelastic material. The simulation results for these drops are verified against theoretical results from literature. Moreover, it is shown that the method can be used to identify the dominant break‐up mechanism of compound drops in relation to the specific non‐Newtonian character of the membrane. Copyright © 1999 John Wiley & Sons, Ltd.     

An experimental study of drop migration in shear flow between concentric cylinders

The phenomenon of migration of liquid drops in Couette flow between concentric cylinders due to non-Newtonian fluid properties and shape deformation has been studied experimentally. The results agree very well with the theory of Chan and Leal, which included the effect of hydrodynamic interaction with the bounding walls, and that of velocity profile curvature in a Couette device. Significant observations that were not reported in previous studies include the migration of a deformable Newtonian drop to an equilibrium position between the centerline and the inner rotor, and the competition between normal stresses and shape deformation effects for the case of a Newtonian drop in a non-Newtonian fluid.     

Deformation and breakup of viscoelastic drops in planar extensional flows

A computer-controlled four-roll mill was used to investigate the deformation and break-up of polymeric drops in the well-characterized flow of an immiscible Newtonian fluid. Aqueous polymer solutions ranging in concentration from 160 ppm to 3% by weight were examined. For zero-shear-rate viscosity ratios greater than order 1, the deformation of the drops closely followed that of Newtonian fluids, irrespective of the droplet material. However, drops with viscosity ratios less than order 1 had significantly smaller critical deformations and the critical capillary number was found to be substantially smaller. Two modes of drop break-up were discovered that differed substantially from that observed for Newtonian drops in the inclusion of cusped ends and tip streaming.     

A boundary‐only approach to the deformation of a shear‐thinning drop in extensional Newtonian flow

The influence of shear thinning on drop deformation is examined through a numerical simulation. A two‐dimensional formulation within the scope of the boundary element method (BEM) is proposed for a drop driven by the ambient flow inside a channel of a general shape, with emphasis on a convergent–divergent channel. The drop is assumed to be shear thinning, obeying the Carreau–Bird model and the suspending fluid is Newtonian. The viscosity of the drop at any time is estimated on the basis of a rate‐of‐strain averaged over the region occupied by the drop. The viscosity thus changes from one time step to the next, and it is strongly influenced by drop deformation. It is found that small drops, flowing on the axis, elongate in the convergent part of the channel, then regain their spherical form in the divergent part; thus confirming experimental observations. Newtonian drops placed off‐axis are found to rotate during the flow with the period related to the initial extension, i.e. to the drop aspect ratio. This rotation is strongly prohibited by shear thinning. The formulation is validated by monitoring the local change of viscosity along the interface between the drop and the suspending fluid. It is found that the viscosity averaged over the drop compares, generally to within a few per cent, with the exact viscosity along the interface.     

Deformation and breakup of a non-Newtonian slender drop in an extensional flow

The deformation and breakup of a non-Newtonian slender drop in a Newtonian liquid in a simple extensional and creeping flow has been theoretically studied. The power-law was chosen for the fluid inside the drop, and the deformation of the drop is described by a single ordinary differential equation, which was numerically solved. Asymptotic analytical expressions for the local radius were derived near the center and close to the end of the drop. The results for the shape of the drop and the breakup criterion are presented as a function of the capillary number, the viscosity ratio and type of non-Newtonian fluid inside the drop. An approximate analytical solution is also suggested which is in good agreement with the numerical results.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

Models for the deformation of a single ellipsoidal drop: a review

Dilute polymer blends and immiscible liquid emulsions are characterized by a globular morphology. The dynamics of a single drop subjected to an imposed flow field has been considered to be a valuable model system to get information on dilute blends. This problem has been studied either theoretically by developing exact theories for small drop deformations or by developing simplified models often based on phenomenological assumptions. In this paper, a critical overview of the available models for the dynamics of a single drop is presented, discussing four different systems, namely the Newtonian system, where a single Newtonian drop is immersed in an infinite Newtonian matrix; the non-Newtonian system, where at least one of the components, the drop fluid or the matrix one, is non-Newtonian; the confined Newtonian system, where the matrix is confined and wall effects alter the drop dynamics; and the confined non-Newtonian system.     

Unsteady flow of a fourth-grade fluid due to an oscillating plate

The governing non-linear high-order, sixth-order in space and third-order in time, differential equation is constructed for the unsteady flow of an incompressible conducting fourth-grade fluid in a semi-infinite domain. The unsteady flow is induced by a periodically oscillating two-dimensional infinite porous plate with suction/blowing, located in a uniform magnetic field. It is shown that by augmenting additional boundary conditions at infinity based on asymptotic structures and transforming the semi-infinite physical space to a bounded computational domain by means of a coordinate transformation, it is possible to obtain numerical solutions of the non-linear magnetohydrodynamic equation. In particular, due to the unsymmetry of the boundary conditions, in numerical simulations non-central difference schemes are constructed and employed to approximate the emerging higher-order spatial derivatives. Effects of material parameters, uniform suction or blowing past the porous plate, exerted magnetic field and oscillation frequency of the plate on the time-dependent flow, especially on the boundary layer structure near the plate, are numerically analysed and discussed. The flow behaviour of the fourth-grade non-Newtonian fluid is also compared with those of the Newtonian fluid.     

Three‐dimensional boundary element analysis of drop deformation in confined flow for Newtonian and viscoelastic systems

An adaptive (Lagrangian) boundary element approach is proposed for the general three‐dimensional drop deformation in confined flow. The adaptive method is stable as it includes remeshing capabilities of the deforming interface between drop and suspending fluid, and thus can handle large deformations. Both drop and surrounding fluid are viscous incompressible and can be Newtonian or viscoelastic. A boundary‐only formulation is implemented for fluids obeying the linear Jeffrey's constitutive equation. Similarly to the formulation for two‐dimensional Newtonian fluids (Khayat RE, Luciani A, Utracki LA. Boundary element analysis of planar drop deformation in confined flow. Part I. Newtonian fluids. Engineering Analysis of Boundary Elements 1997; 19 : 279), the method requires the solution of two simultaneous integral equations on the interface between the two fluids and the confining solid boundary. Although the problem is formulated for any confining geometry, the method is illustrated for a deforming drop as it is driven by the ambient flow inside a cylindrical tube. The accuracy of the method is assessed by comparison with the analytical solution for two‐phase radial spherical flow, leading to good agreement. The influence of mesh refinement is examined for a drop in simple shear flow. Copyright © 2000 John Wiley & Sons, Ltd.     

Derivation of a spectral pressureless formulation for direct numerical simulation of incompressible channel flows based on a functional formalism

In this article, a new computational spectral algorithm is developed for simulation of general three-dimensional, time-dependent, incompressible channel flow. The development is based on a general functional formalism of non-equilibrium thermodynamics, and, although it is illustrated here for a Newtonian fluid, it is easily adapted to non-Newtonian fluids. The advantage of this algorithm is that the scalar pressure is eliminated from the discrete spectral analog to the equations of motion, which are expressed solely in terms of the spectral coefficients of the velocity vector field. This alleviates the need for the application of boundary conditions on the pressure, the specification of which can be a major source of difficulty in direct numerical simulations. At the same time, the velocity spectrum is quite general, and not subject to any a priori constraints. Thus, it is anticipated that the ideas exposed in the present algorithm can lead to the development of better numerical simulation techniques for complicated three-dimensional and turbulent flows.     

Motion and shape of an axisymmetric viscoplastic drop slowly falling through a viscous fluid

Slow sedimentation of a deformable drop of Bingham fluid in an unbounded Newtonian medium is studied using a variation of the integral equation method (Toose et al., J Eng Math 30:131–150, 1996, Int J Numer Methods Fluids 30:653–674, 1999). The Green function for the Stokes equation is used, and the non-Newtonian stress is treated as a source term. The computations are performed for a range of physical parameters of the system. It is demonstrated that initially deformed drop similar to Newtonian ones breaks up for high capillary number, Ca, and stabilizes to steady shapes at low Ca. Estimations of critical capillary number for specific initial deformations demonstrated its growth (increase in the stability of the drop) with the yield stress magnitude both for prolate and oblate initial shapes. Prolate initial shapes become more stable with the increase of the plastic viscosity. In contrast to this, for low yield stress, oblate shapes are destabilized with the growth of the plastic viscosity. This effect is similar to the effect of the viscosity of a Newtonian drop on its stability. However, at higher yield stress, the effect of plastic viscosity is reversed.     

Displacement of a Newtonian fluid by a non-Newtonian fluid in a porous medium

This paper presents an analytical Buckley-Leverett-type solution for one-dimensibnal immiscible displacement of a Newtonian fluid by a non-Newtonian fluid in porous media. The non-Newtonian fluid viscosity is assumed to be a function of the flow potential gradient and the non-Newtonian phase saturation. To apply this method to field problems a practical procedure has been developed which is based on the analytical solution and is similar to the graphic technique of Welge. Our solution can be regarded as an extension of the Buckley-Leverett method to Non-Newtonian fluids. The analytical result reveals how the saturation profile and the displacement efficiency are controlled not only by the relative permeabilities, as in the Buckley-Leverett solution, but also by the inherent complexities of the non-Newtonian fluid. Two examples of the application of the solution are given. One application is the verification of a numerical model, which has been developed for simulation of flow of immiscible non-Newtonian and Newtonian fluids in porous media. Excellent agreement between the numerical and analytical results has been obtained using a power-law non-Newtonian fluid. Another application is to examine the effects of non-Newtonian behavior on immiscible displacement of a Newtonian fluid by a power-law non-Newtonian fluid.     

A new mesh relaxation approach and automatic time‐step control method for boundary integral simulations of a viscous drop

We present a numerical methodology for the simulation of a viscous drop under simple shear flows by using the boundary integral method. The present work treats only a single drop in an unbounded fluid‐flow, but the results can be directly applied to studies on the rheology of dilute emulsions, in which the hydrodynamic interactions between two or more drops can be neglected. Singular and non‐singular integral representations of the velocity field are considered. Several aspects of the method are presented, including a new mesh relaxation approach and an automatic time‐step control method. The relaxation strategy is used in order to contain the distortion of the mesh and is performed by using relaxation iterations in a virtual temporal march between each physical time step of the simulation and monitoring the standard deviation of the areas of the elements. The automatic time‐step control method uses a global quantity related to the drop deformation in order to automatically set the temporal integration time step. It is carried out in a way to keep the local integration error less than a given tolerance. This strategy reduces the computational cost of the simulation by dramatically reducing the number of time steps in the temporal integration process. Copyright © 2016 John Wiley & Sons, Ltd.     

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

Direct numerical simulation of the turbulent energy balance and the shear stresses in power-law fluid flows in pipes

The results of direct numerical simulation of turbulent flows of non-Newtonian pseudoplastic fluids in a straight pipe are presented. The data on the distributions of the turbulent stress tensor components and the shear stress and turbulent kinetic energy balances are obtained for steady turbulent flows at the Reynolds numbers of 10⁴ and 2×10⁴. As distinct from Newtonian fluid flows, the viscous shear stresses turn out to be significant even far from the wall. In power-law fluid flows the mechanism of the energy transport from axial to transverse component fluctuations is suppressed. It is shown that with decrease in the fluid index the turbulent transfer of the momentum and the velocity fluctuations between the wall layer and the flow core reduces, while the turbulent energy flux toward the wall increases. The earlier-proposed models for the average viscosity and the non-Newtonian one-point correlations are in good agreement with the data of direct numerical simulation.     

Experimental investigation of viscoelastic drop deformation in Newtonian matrix at high capillary number under simple shear flow

The steady state deformation of a viscoelastic drop (Boger fluid) in a Newtonian liquid at high capillary number under simple shear flow is investigated by direct visualization using a specially designed Couette apparatus which enables visualization from two perpendicular directions. Two drop deformation modes are found: (1) Mode I – drop deformation in the flow direction and (2) Mode II – drop deformation in the vorticity direction. The drop deformation mode depends on the relative strength of the elastic contribution to viscous contribution. If the elastic contribution is weak compared to the viscous contribution, the drop elongates in the flow direction via Mode I. If the elastic contribution is strong, the drop elongates in the vorticity direction via Mode II. The drop size also affects the drop deformation. At the same capillary number, bigger drops have larger deformations than smaller drops.     

Numerical simulation of the non-Newtonian blood flow through a mechanical aortic valve

This work focuses on the comparison between Newtonian and non-Newtonian blood flows through a bileaflet mechanical heart valve in the aortic root. The blood, in fact, is a concentrated suspension of cells, mainly red blood cells, in a Newtonian matrix, the plasma, and consequently its overall behavior is that of a non-Newtonian fluid owing to the action of the cells’ membrane on the fluid part. The common practice, however, assumes the blood in large vessels as a Newtonian fluid since the shear rate is generally high and the effective viscosity becomes independent of the former. In this paper, we show that this is not always the case even in the aorta, the largest artery of the systemic circulation, owing to the pulsatile and transitional nature of the flow. Unexpectedly, for most of the pulsating cycle and in a large part of the fluid volume, the shear rate is smaller than the threshold level for the blood to display a constant effective viscosity and its shear thinning character might affect the system dynamics. A direct inspection of the various flow features has shown that the valve dynamics, the transvalvular pressure drop and the large-scale features of the flow are very similar for the Newtonian and non-Newtonian fluid models. On the other hand, the mechanical damage of the red blood cells (hemolysis), induced by the altered stress values in the flow, is larger for the non-Newtonian fluid model than for the Newtonian one.     

Spreading of inkjet droplet of non-Newtonian fluid on solid surface with controlled contact angle at low Weber and Reynolds numbers

The dynamics of inkjet droplet of non-Newtonian fluid on glass substrates was investigated experimentally and compared with that of Newtonian fluid. The non-Newtonian fluids used here were 100 ppm solutions of polyethylene oxide (300k, 600k and 900k) dissolved in the 1:1 mixture of water and glycerin. Weber number (We) was 2–35 and Ohnesorge number was fixed at 0.057 ± 0.003. The wettability of solid substrate was also varied. The diameter of inkjet droplets in the present study was about 50 μm and was much smaller than the size of the previous studies on drop impact. Due to the development of a thin and long thread at the rear of the main drop the jetting window of polymer solution was much narrower than that of Newtonian fluid, and hence the experimental range of Weber number was restricted. The impact scenarios of non-Newtonian inkjet droplets were found to be qualitatively different from those of Newtonian droplets during the receding phase while they were almost the same as the Newtonian fluid case during the kinematic phase. The spreading diameter at the equilibrium was well correlated with the modified Weber number (We′ = We/(1 − cos θ_eq)) as in the case of Newtonian fluid, where θ_eq is the equilibrium contact angle. The similarity or disparity between the Newtonian and non-Newtonian cases was discussed considering the conformation of polymer chains during each stage of drop deformation.