Simulations of deformable systems in fluids under shear flow using an arbitrary Lagrangian Eulerian technique
Introduction
Microfluidics, i.e. the processing of fluids in channels which have a characteristic length in the order of microns, has gained an increasing importance over the last 15 years. In several microfluidic systems, as cytometric, sorting and diagnostic devices, suspensions of soft particles flow in miniaturized channels [1], [2]: in such a situation, a wide range of phenomena can arise due to the deformability and the elasticity of the suspended objects, the geometry of the channels, the interactions between the objects and the walls of the flow cells, and eventually the complexity of the suspending media.
Among deformable systems, drops have been extensively investigated since the 1930s, when Taylor experimentally studied the deformation of Newtonian drops in a Newtonian fluid under simple shear flow in the low capillary number regime [3], [4]. For drops, the capillary number, which is the ratio between viscous and interfacial forces, is defined as , where is the ambient fluid viscosity, is the radius of the undeformed drop, is the shear rate and is the interfacial tension between the two fluids. The drops are found to deform and assume an ellipsoidal shape with a fixed orientation with respect to the flow direction, the deformation at the steady state linearly increasing with Ca according to the law:where D is the Taylor deformation parameter, defined as the ratio of the difference and the sum of the major and minor semiaxes of the ellipsoid in the shear plane, and is the ratio between the viscosities of the drop and the suspending matrix. During both the transient and steady deformation, the fluid bulk of the drops circulates around the center of mass. Taylor’s results have been afterwards confirmed theoretically, experimentally and computationally [5], [6], [7]. Kennedy and co-authors [8] showed that if the viscosity ratio is less than about 4, when Ca reaches a critical value, the interfacial tension is no more capable of contrasting the deformation and the drop breaks up. On the contrary, above such critical -value, the drop does not break up even at high Ca-values. Several papers focus on identifying break-up conditions for drops in shear flow [6], [7], [9]. The presence of solid walls in the vicinity of a drop can substantially influence its deformation: in confined shear flow, very elongated non-ellipsoidal deformed shapes have been detected [10]. If one considers fluid complexity and in particular viscoelasticity, a new parameter has to be taken into account: the Deborah number De, which is the ratio between the fluid and the flow characteristic times. On whatever side viscoelasticity acts (i.e. drop, matrix or both), in shear flow the stationary orientation angle of a drop is lower than in the corresponding Newtonian/Newtonian case [11]. A phenomenological model is available for drops under all flow conditions in the small deformation regime, which can handle viscoelasticity in both the matrix and the dispersed phase [12].
Elastic particles are quite interesting from both a scientific and a technological point of view, since they are present, for example, in filled polymers and are also models for more complicated systems. Moreover, not much attention has been devoted to such objects and a wide comprehension of their behaviors in flow is still lacking. In 1946, Fröhlich and Sack [13] were the first to investigate the behavior of a suspension of elastic spheres in a Newtonian fluid undergoing extensional flow, linking the extensional stress of the suspension to the strain rate. In 1967, Roscoe [14] studied the behavior of a dilute suspension of viscoelastic spheres in a Newtonian fluid subjected to shear flow, obtaining from a theoretical point of view the deformation of the objects, the stress and the viscosity of the suspension. He showed that the particles attain an ellipsoidal steady deformed shape with fixed orientation, and that the material within the ellipsoids continuously rotates. In the same year, Goddard and Miller [15] derived a constitutive equation for dilute suspensions of slightly deformed elastic spheres. In 1981, Murata [16] studied the deformation of a spherical elastic particle in an arbitrary flow field by using a perturbation analysis. Since then, almost nothing has been done on elastic particles until the end of last decade, when Gao and Hu [17] performed a 2D perturbation analysis from which they derived the relationship between the steady elliptical deformation of an elastic particle, expressed by the Taylor deformation parameter D, and the dimensionless parameter which measures the relative weights of the viscous and elastic forces in the system (the elastic capillary number Ca). The authors showed that for small Ca-values, at the stationary , and validated such result by means of 2D Arbitrary Lagrangian Eulerian Finite Element Method (ALE FEM) simulations. In 2011, Gao et al. [18] studied the behavior of an initially spherical elastic particle suspended in a Newtonian fluid in shear flow through a polarization method [19], [20], coming to a validation of Roscoe’s results.
More complicated soft systems are vesicles and capsules, that are constituted by liquid drops wrapped in thin elastic membranes. Such objects, which are often non-spherical at rest, are widely used in the biotechnological field as drug deliverers and models for cells (e.g. Red Blood Cells). Several theoretical, experimental, and numerical works are available on the deformation of vesicles and capsules in Newtonian fluids under unconfined shear flow, as reported in their review by Finken and co-authors [21]. Depending on the physical parameters, different regimes of motion are detected, such as stationary tank-treading (TT), where the orientation angle with respect to the flow direction is fixed and the elastic membrane rotates around the internal fluid like the treads of a tank, tumbling (TU), where the whole object periodically rotates around the vorticity axis, vacillating–breathing (VB) (also known as oscillating (OS)), where the orientation angle oscillates around a mean value, and some intermittent regimes, that consist in an evolution of the motion regime in time (e.g., ). Stationary TT and TU being the two extreme behaviors, an increase in is found to promote TU, whereas an increase in the flow intensity promotes TT.
In this paper, an arbitrary Lagrangian Eulerian finite element method based code for viscoelastic fluids using well-known stabilization techniques such as SUPG [22], DEVSS [23], [24], and log-conformation [25], [26] is adapted and extended to study the behavior of drops and elastic particles suspended in Newtonian and viscoelastic fluids under unbounded shear flow in 3D. To the best of our knowledge, we show the deformation of elastic particles in unbounded shear flow of Upper Convected Maxwell (UCM) and Giesekus (Gsk) fluids for the first time.
The paper is organized as follows: in Section 2, the governing equations for the systems considered are presented; in Section 3, the numerical technique adopted to solve the model equations is explained in detail; in Section 4, numerical convergence tests are shown; in Section 5, results are presented; finally, in Section 6, some conclusions are drawn.
Section snippets
A liquid drop in a fluid under shear flow
In Fig. 1, a schematic drawing is reported of an initially spherical liquid drop suspended in a fluid under simple shear flow. For both the suspended and the suspending phase it is assumed that inertia can be neglected and that the volume is constant (incompressible). Therefore, the mass and momentum balance for both the drop and the suspending phase reduce towhere and are the velocity vector and the stress tensor, respectively. can in turn be expressed as:where p is the
Volume discretization
Both the matrix fluid and the droplet fluid (or solid particle) domain are discretized using the finite element method employing a mesh of quadratic tetrahedra. The interface mesh aligns with element faces (quadratic triangles), which are the same on the matrix and droplet side (conforming geometry). In Fig. 6, the mesh employed in the simulations presented in this paper is shown. In the box on the right, a zoom of the interface mesh of the initially spherical suspended object is displayed.
We
Convergence tests
Before running simulations, convergence tests have to be performed in space and time. This also helps in optimizing the mesh resolution and time-step such that the results are sufficiently reliable without wasting too much computational time.
Deformation of a Newtonian drop in a Newtonian fluid
A Newtonian drop is suspended in a Newtonian fluid under unconfined shear flow. As known from Taylor [4], in the moderate capillary number regime, the drop deforms into an ellipsoid and attains a fixed orientation with respect to the direction of the flow. In Fig. 10a, the steady state values of the Taylor deformation parameter emerging from our numerical simulations are plotted vs. the capillary number Ca as black circles. The viscosity ratio is . Such results are compared with the
Conclusions
In this paper, an arbitrary Lagrangian Eulerian finite element method based numerical code for viscoelastic fluids using well-known stabilization techniques (SUPG, DEVSS, log-conformation) is adapted and extended to study the behavior of soft systems such as liquid drops and elastic particles suspended in fluids under unbounded shear flow in 3D. To the best of our knowledge, in this paper, 3D numerical simulations of elastic particles in a viscoelastic fluid are presented for the first time.
References (38)
- et al.
Motion and deformation of liquid drops, and the rheology of dilute emulsions in simple shear flow
Comput Fluids
(1994) - et al.
Deformation of a Newtonian drop in a viscoelastic matrix under steady shear flow experimental validation of slow flow theory
J Non-Newton Fluid Mech
(2003) - et al.
Deformation of elastic particles in viscous shear flow
J Comput Phys
(2009) - et al.
Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier–Stokes equations
Comp Meth Appl Mech Eng
(1982) - et al.
A new mixed finite element method for computing viscoelastic flows
J Non-Newton Fluid Mech
(1995) - et al.
Time dependent finite element analysis of the linear stability of viscoelastic flows with interfaces
J Non-Newton Fluid Mech
(2003) - et al.
Constitutive laws for the matrix-logarithm of the conformation tensor
J Non-Newton Fluid Mech
(2004) - et al.
Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithms
J Non-Newton Fluid Mech
(2005) - et al.
A finite-volume/boundary-element method for flow past interfaces in the presence of surfactants, with application to shear flow past a viscous drop
Comput Fluids
(1998) - et al.
Fictitious boundary and moving mesh methods for the numerical simulation of rigid particulate flows
J Comput Phys
(2007)
Stability analysis of polymer shear flows using the extended pom–pom constitutive equations
J Non-Newton Fluid Mech
Decoupled second-order transient schemes for the flow of viscoelastic fluids without a viscous solvent contribution
J Non-Newton Fluid Mech
An algorithm for the simulation of transient viscoelastic flows with free surfaces
J Comput Phys
Direct numerical simulations of fluid–solid systems using the arbitrary Lagrangian–Eulerian technique
J Comput Phys
A second-order theory for shear deformation of drops
J Colloid Interface Sci
Lab-on-a-chip: microfluidics in drug discovery
Nat Rev Drug Discov
Continuous flow separations in microfluidic devices
Lab Chip
The viscosity of a fluid containing small drops of another fluid
Proc R Soc Lond A
The formation of emulsions in definable fields of flow
Proc R Soc Lond A
Cited by (63)
Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method
2023, ParticuologyCitation Excerpt :Jaensson et al. (2015) increased the efficiency of simulations of up to 64 particles suspended in viscoelastic shear flows. Villone et al. (2014) extended this method to simulate a 3D deformed particle suspended in a viscoelastic fluid. Choi et al. (2010) proposed an extended finite element method (FEM)–temporary ALE method, which dramatically increased the computing efficiency, but it still simulated only two moving particles.
Dynamics and rheology of soft particles
2023, Advances in Chemical EngineeringOn the hydrodynamic mutual interactions among cells for high-throughput microfluidic holographic cyto-tomography
2022, Optics and Lasers in EngineeringAnisotropic Oldroyd-type models for non-colloidal suspensions of viscoelastic particles in Newtonian and yield-stress fluids via homogenization
2021, Journal of Non-Newtonian Fluid MechanicsCitation Excerpt :However, thus far, this has only been attempted for small strains and, therefore, this approach could only be used for small oscillatory shear conditions. On the other hand, numerical simulations have already demonstrated the effect of viscoelasticity in the suspending fluid for dilute suspensions [78]. This work has focused on the use of the Hashin–Shtrikman–Willis-type estimates of Ponte Castañeda and Willis [21] for the suspensions, as a compromise between accuracy and complexity.
Computational interfacial rheology
2021, Journal of Non-Newtonian Fluid MechanicsBenchmark solutions for flows with rheologically complex interfaces
2020, Journal of Non-Newtonian Fluid Mechanics