Bioprinting technologies rely on the formation of soft gel drops for printing tissue scaffolds and the dynamics of these drops can affect the process. A model is developed to describe the oscillations of a spherical gel drop with finite shear modulus, whose interface is held by surface tension. The governing elastodynamic equations are derived and a solution is constructed using displacement potentials decomposed into a spherical harmonic basis. The resulting nonlinear characteristic equation depends upon two dimensionless numbers, elastocapillary and compressibility, and admits two types of solutions, (i) spheroidal (or shape change) modes and (ii) torsional (rotational) modes. The torsional modes are unaffected by capillarity, whereas the frequency of shape oscillations depend upon both the elastocapillary and compressibility numbers. Two asymptotic dispersion relationships are derived and the limiting cases of the inviscid Rayleigh drop and elastic globe are recovered. For a fixed polar wavenumber, there exists an infinity of radial modes that each transition from an elasticity wave to a capillary wave upon increasing the elastocapillary number. At the transition, there is a qualitative change in the deformation field and a set of recirculation vortices develop at the free surface. Two special modes that concern volume oscillations and translational motion are characterized. A new instability is documented that reflects the balance between surface tension and compressibility effects due to the elasticity of the drop.
more »
« less
Pair interactions between viscous drops in a viscoelastic matrix in free shear: Transition from passing to tumbling trajectories
Shear-induced pair interactions between viscous drops suspended in a viscoelastic matrix are numerically investigated examining the effects of elasticity and drop deformability on their post-collision trajectory. Two different trajectory types are identified depending on the Weissenberg number Wi and capillary number Ca. Drops suspended in a Newtonian matrix ( Wi = 0.0) show a passing trajectory where drops slide past each other and separate in the stream-wise direction. However, when increasing the Weissenberg number above a critical value, a tumbling/doublet trajectory is observed where two drops rotate around the midpoint of the line joining their centers, as was also seen previously for rigid particles. The tumbling trajectory is explained by investigating the flow around a single drop in shear. Elasticity generates a larger region of spiraling streamlines around a drop, which, during a pair interaction, traps the second drop giving rise to the tumbling pair. Decreasing deformability (lower Ca) and increasing viscoelasticity (higher Wi) favor a tumbling trajectory. With simulations sweeping the parameter space, we obtain a phase plot of the two different trajectories as functions of Ca and Wi. Treating the tension along the curved streamlines due to the non-zero first normal stress difference in the viscoelastic medium as an enhancement to the interfacial tension, we have developed an approximate force balance model for the zone of spiraling streamlines. It qualitatively captures the observed scaling of the critical Ca and Wi values at the phase boundary. The effects of unequal size, initial configuration, and non-unity viscosity ratio are briefly investigated.
more » « less- Award ID(s):
- 2019507
- NSF-PAR ID:
- 10364437
- Publisher / Repository:
- Society of Rheology
- Date Published:
- Journal Name:
- Journal of Rheology
- Volume:
- 66
- Issue:
- 3
- ISSN:
- 0148-6055
- Page Range / eLocation ID:
- p. 571-584
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
We analyse the electrophoresis of a weakly charged particle with a thin double layer in a dilute polymer solution. The particle velocity in polymer solutions modelled with different constitutive equations is calculated using a regular perturbation in the polymer concentration and the generalized reciprocal theorem. The analysis shows that the polymer is strongly stretched in two regions, the birefringent strand and the high-shear region inside the double layer. The electrophoretic velocity of the particle always decreases with the addition of polymers due to both increased viscosity and fluid elasticity. At a small Weissenberg number ( $Wi$ ), which is the product of the polymer relaxation time and the shear rate, the polymers inside the double layer contribute to most of the velocity reduction by increasing the fluid viscosity. With increasing $Wi$ , viscoelasticity decreases and shear thinning increases the particle velocity. Polymer elasticity alters the fluid velocity disturbance outside the double layer from that of a neutral squirmer to a puller-type squirmer. At high $Wi$ , the strong extensional stress inside the birefringent strand downstream of the particle dominates the velocity reduction. The scaling of the birefringent strand is used to estimate the particle velocity.more » « less
-
Growth of the microfluidics field has triggered numerous advances in focusing and separating microparticles, with such systems rapidly finding applications in biomedical, chemical, and environmental fields. The use of shear-thinning viscoelastic fluids in microfluidic channels is leading to evolution of elasto-inertial focusing. Herein, we showed that the interplay between the elastic and shear-gradient lift forces, as well as the secondary flow transversal drag force that is caused by the non-zero second normal stress difference, lead to different particle focusing patterns in the elasto-inertial regime. Experiments and 3D simulations were performed to study the effects of flowrate, particle size, and the shear-thinning extent of the fluid on the focusing patterns. The Giesekus constitutive equation was used in the simulations to capture the shear-thinning and viscoelastic behaviors of the solution used in the experiments. At low flowrate, with Weissenberg number Wi ~ O(1), both the elastic force and secondary flow effects push particles towards the channel center. However, at a high flowrate, Wi ~ O(10), the elastic force direction is reversed in the central regions. This remarkable behavior of the elastic force, combined with the enhanced shear-gradient lift at the high flowrate, pushes particles away from the channel center. Additionally, a precise prediction of the focusing position can only be made when the shear-thinning extent of the fluid is correctly estimated in the modeling. The shear-thinning also gives rise to the unique behavior of the inertial forces near the channel walls which is linked with the ‘warped’ velocity profile in such fluids.more » « less
-
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
