skip to main content

Title: Nonlinear rheology of entangled wormlike micellar solutions predicted by a micelle-slip-spring model
We examine linear and nonlinear shear and extensional rheological properties using a “micelle-slip-spring model” [T. Sato et al., J. Rheol. 64, 1045–1061 (2020)] that incorporates breakage and rejoining events into the slip-spring model originally developed by Likhtman [Macromolecules 38, 6128–6139 (2005)] for unbreakable polymers. We here employ the Fraenkel potential for main chain springs and slip-springs to address the effect of finite extensibility. Moreover, to improve extensional properties under a strong extensional flow, stress-induced micelle breakage (SIMB) is incorporated into the micelle-slip-spring model. Thus, this model is the first model that includes the entanglement constraint, Rouse modes, finite extensibility, breakage and rejoining events, and stress-induced micelle breakage. Computational expense currently limits the model to micellar solutions with moderate numbers of entanglements ([Formula: see text]), but for such solutions, nearly quantitative agreement is attained for the start-up of the shearing flow. The model in the extensional flow cannot yet be tested owing to the lack of data for this entanglement level. The transient and steady shear properties predicted by the micelle-slip-spring model for a moderate shear rate region without significant chain stretch are fit well by the Giesekus model but not by the Phan–Thien/Tanner (PTT) model, which is consistent with the ability of the Giesekus model to match experimental shear data. The extensional viscosities obtained by the micelle-slip-spring model with SIMB show thickening followed by thinning, which is in qualitative agreement with experimental trends. Additionally, the extensional rheological properties of the micelle-slip-spring model with or without SIMB are poorly predicted by both the Giesekus and the PTT models using a single nonlinear parameter. Thus, future work should seek a constitutive model able to capture the behavior of the slip-spring model in shear and extensional flows and so provide an accurate, efficient model of micellar solution rheology.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Rheology
Page Range / eLocation ID:
639 to 656
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Poole, Robert J (Ed.)
    The FENE-P (Finitely-Extensible Nonlinear Elastic) dumbbell constitutive equation is widely used in simulations and stability analyses of free and wall-bounded viscoelastic shear flows due to its relative simplicity and accuracy in predicting macroscopic properties of dilute polymer solutions. The model contains three independent material parameters, which expressed in dimensionless form correspond to a Weissenberg number (Wi), i.e., the ratio of the dumbbell relaxation time scale to a characteristic flow time scale, a finite extensibility parameter (L), corresponding to the ratio of the fully extended dumbbell length to the root mean square end-to-end separation of the polymer chain under equilibrium conditions, and a solvent viscosity ratio. An exact solution for the rheological predictions of the FENE-P model in steady simple shear flow is available [Sureshkumar et al., Phys Fluids (1997)], but the resulting nonlinear and nested set of equations do not readily reveal the key shear-thinning physics that dominates at high Wi as a result of the finite extensibility of the polymer chain. In this note we review a simple way of evaluating the steady material functions characterizing the nonlinear evolution of the polymeric contributions to the shear stress and first normal stress difference as the shear rate increases, provide asymptotic expansions as a function of Wi , and show that it is in fact possible to construct universal master curves for these two material functions as well as the corresponding stress ratio. Steady shear flow experiments on three highly elastic dilute polymer solutions of different finite extensibilities also follow the identified master curves. The governing dimensionless parameter for these master curves is Wi/L and it is only in strong shear flows exceeding Wi/L > 1 that the effects of finite extensibility of the polymer chains dominate the evolution of polymeric stresses in the flow field. We suggest that reporting the magnitude of Wi/L when performing stability analyses or simulating shear-dominated flows with the FENE-P model will help clarify finite extensibility effects. 
    more » « less
  2. We investigate aspects of the spherical squirmer model employing both large-scale numerical simulations and asymptotic methods when the squirmer is placed in weakly elastic fluids. The fluids are modelled by differential equations, including the upper-convected Maxwell (UCM)/Oldroyd-B, finite-extensibility nonlinear elastic model with Peterlin approximation (FENE-P) and Giesekus models. The squirmer model we examine is characterized by two dimensionless parameters related to the fluid velocity at the surface of the micro-swimmer: the slip parameter $\xi $ and the swirl parameter $\zeta $ . We show that, for swimming in UCM/Oldroyd-B fluids, the elastic stress becomes singular at a critical Weissenberg number, Wi , that depends only on $\xi$ . This singularity for the UCM/Oldroyd-B models is independent of the domain exterior to the swimmer, or any other forces considered in the momentum balance for the fluid – we believe that this is the first time such a singularity has been explicitly demonstrated. Moreover, we show that the behaviour of the solution at the poles is purely extensional in character and is the primary reason for the singularity in the Oldroyd-B model. When the Giesekus or the FENE-P models are utilized, the singularity is removed. We also investigate the mechanism behind the speed and rotation rate enhancement associated with the addition of swirl in the swimmer's gait. We demonstrate that, for all models, the speed is enhanced by swirl, but the mechanism of enhancement depends intrinsically on the rheological model employed. Special attention is paid to the propulsive role of the pressure and elucidated upon. We also study the region of convergence of the perturbation solutions in terms of Wi . When techniques that accelerate the convergence of series are applied, transformed solutions are derived that are in very good agreement with the results obtained by simulations. Finally, both the analytical and numerical results clearly indicate that the low- Wi region is more important than one would expect and demonstrate all the major phenomena observed when swimming with swirl in a viscoelastic fluid. 
    more » « less
  3. We investigate the effects of micellar entanglement number on the kinetics of shear banding flow formation in a Taylor–Couette flow. Three sets of wormlike micellar solutions, each set with a similar fluid elasticity and zero-shear-rate viscosity, but with varying entanglement densities, are studied under the startup of steady shear. Our experiments indicate that in the set with low fluid elasticity, the transient shear banding flow is characterized by the formation of a transient flow reversal in a range of entanglement densities. Outside of this range, the transient flow reversal is not observed. For the sets of medium and high elasticities, the transient flow reversals exist for relatively small entanglement densities and disappear for large entanglement densities. Our analysis shows that wall slip and elastic instabilities do not affect the transient flow feature. We identify a correlation between micellar entanglement number, the width of the stress plateau, and the extent of the transient flow reversal. As the micellar entanglement number increases, the width of the stress plateau first increases; then, at a higher micellar entanglement number, the plateau width decreases. Therefore, we hypothesize that the transient flow reversal is connected to the micellar entanglement number through the width of the stress plateau. 
    more » « less
  4. We combine a slip-spring model with an ‘entangled kink dynamics’ (EKD) model for strong uniaxial extensional flows (with Rouse Weissenberg number W i R ≫ 1 ) of long ( M w > 1   Mkg / mol for polystyrene) entangled polymers in solutions and melts. The slip-spring model captures the dynamics up to the formation of a ‘kinked’ or folded state, while the kink dynamics simulation tracks the dynamics from that point forward to complete extension. We show that a single-chain slip-spring model using affine motion of the slip-spring anchor points produces unrealistically high tension near the center of the chain once the Hencky strain exceeds around unity or so, exceeding the maximum tension that a chain entangled with a second chain is able to support. This unrealistic tension is alleviated by pairing the slip links on one chain with those on a second chain, and allowing some of the large tension on one of the two to be transferred to the second chain, producing non-affine motion of each. This explicit pairing of entanglements mimics the entanglement pairing also used in the EKD model, and allows the slip spring simulations to be carried out to strains high enough for the EKD model to become valid. We show that results nearly equivalent to those from paired chains are obtained in a single-chain slip-spring simulation by simply specifying that the tension in a slip spring cannot exceed the theoretical maximum value of ζ ′ ϵ ˙ L 2 / 8 where ζ ′ , ϵ ˙ and L are the friction per unit length, strain rate and contour length of the chain, respectively. The effects of constraint release (CR) and regeneration of entanglements is also studied and found to have little effect on the chain statistics up to the formation of the kinked state. The resulting hybrid model provides a fast, simple, simulation method to study the response of high molecular weight ( M w > 1   Mkg / mol ) polymers in fast flows ( W i R ≫ 1 ), where conventional simulation techniques are less applicable due to computational cost. 
    more » « less
  5. We report experiments on flow of wormlike micellar solutions past a falling sphere. By increasing the salt-to-surfactant concentration ratio, and beyond a viscosity peak, wormlike micelles experience a transition from linear to branched microstructure. Two viscoelastic wormlike micelles with salt to surfactant concentrations on each side of the viscosity peak are considered. Our results indicate three significant differences in flows of branched and linear micelles. First, while the sphere drag correction factor rapidly decreases upon increasing Weissenberg number in linear micelles, it shows an apparent local maximum at Wi ≈ 3 in branched micelles. Second, despite its high viscoelasticity, the time-averaged flow of branched micelles around the falling sphere exhibits a fore-and-aft symmetry, while a strong negative wake is observed in linear micelles at relatively weaker flows. Third, branched micelles exhibit a stronger flow-induced birefringence than linear micelles in an otherwise identical condition. Our hypothesis is that subject to strong flows around the falling sphere, branched micelles can relax much more efficiently than linear wormlike micelles through sliding of the branched junctions. This additional stress relaxation mechanism may facilitate micellar orientation, produce a marginal sphere drag reduction and a Newtonian-like flow profile around the falling sphere. Finally, unsteady flow is observed in both linear and branched micellar solutions beyond some critical thresholds of the extensional Weissenber number. Our results corroborate a recently proposed criterion for onset of instability in flow of wormlike micelles past a falling sphere, thereby, suggesting that micellar branching does not affect the mechanism of flow instability. 
    more » « less