A soft viscoelastic drop has dynamics governed by the balance between surface tension, viscosity, and elasticity, with the material rheology often being frequency dependent, which are utilized in bioprinting technologies for tissue engineering and drop-deposition processes for splash suppression. We study the free and forced oscillations of a soft viscoelastic drop deriving (1) the dispersion relationship for free oscillations, and (2) the frequency response for forced oscillations, of a soft material with arbitrary rheology. We then restrict our analysis to the classical cases of a Kelvin–Voigt and Maxwell model, which are relevant to soft gels and polymer fluids, respectively. We compute the complex frequencies, which are characterized by an oscillation frequency and decay rate, as they depend upon the dimensionless elastocapillary and Deborah numbers and map the boundary between regions of underdamped and overdamped motions. We conclude by illustrating how our theoretical predictions for the frequency-response diagram could be used in conjunction with drop-oscillation experiments as a “drop vibration rheometer”, suggesting future experiments using either ultrasonic levitation or a microgravity environment.
- Award ID(s):
- 1750208
- NSF-PAR ID:
- 10139814
- Date Published:
- Journal Name:
- Soft Matter
- Volume:
- 15
- Issue:
- 45
- ISSN:
- 1744-683X
- Page Range / eLocation ID:
- 9244 to 9252
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
null (Ed.)A soft cylindrical interface endowed with surface tension can be unstable to wavy undulations. This is known as the Plateau–Rayleigh instability (PRI) and for solids the instability is governed by the competition between elasticity and capillarity. A dynamic stability analysis is performed for the cases of a soft (i) cylinder and (ii) cylindrical cavity assuming the material is viscoelastic with power-law rheology. The governing equations are made time-independent through the Laplace transform from which a solution is constructed using displacement potentials. The dispersion relationships are then derived, which depend upon the dimensionless elastocapillary number, solid Deborah number, and compressibility number, and the static stability limit, critical disturbance, and maximum growth rate are computed. This dynamic analysis recovers previous literature results in the appropriate limits. Elasticity stabilizes and compressibility destabilizes the PRI. For an incompressible material, viscoelasticity does not affect stability but does decrease the growth rate and shift the critical wavenumber to lower values. The critical wavenumber shows a more complex dependence upon compressibility for the cylinder but exhibits a predictable trend for the cylindrical cavity.more » « less
-
Abstract Very long period (VLP, 2–100 s) seismic signals at basaltic volcanoes like Kilauea, Hawai‘i, and Mount Erebus, Antarctica are likely from resonant oscillations of magma within the shallow plumbing system. The system consists of conduits connected to cracks (dikes and sills) or reservoirs of other shapes. A quantitative understanding of wave propagation and resonance in a coupled conduit‐crack system is required to interpret observations. In this work, we idealize the system as an axisymmetric conduit coupled to a tabular crack, accounting for fluid inertia, compressibility, and viscosity, buoyancy, and crack wall elasticity. We perform time domain simulations and eigenmode analyses of the governing equations, linearized about a rest state. The fundamental mode or conduit‐reservoir mode reflects the balance of conduit magma inertia with buoyancy (and, for small cracks, crack wall elasticity). Magma oscillates in an effectively incompressible manner within the conduit, deflating and inflating the crack, which couples to the surrounding solid to produce observable surface displacements. For sufficiently low viscosity magmas, viscous effects are confined to boundary layers. Shorter period modes are primarily reverberating crack waves with negligible coupling to the conduit. Finally, we introduce an approximate reduced model for the conduit‐reservoir mode, which can also handle more general reservoir geometries (e.g., spherical chambers). The reduced model connects the observable VLP period and quality factor to two uniquely constrained parameters: the inviscid oscillation period
T 0and the viscous diffusion timeτ v i s across the conduit radius. Our models can be extended to study the seismic response of more complex magmatic systems. -
Surface waves are excited at the boundary of a mechanically vibrated cylindrical container and are referred to as edge waves. Resonant waves are considered, which are formed by a travelling wave formed at the edge and constructively interfering with its centre reflection. These waves exhibit an axisymmetric spatial structure defined by the mode number $n$ . Viscoelastic effects are investigated using two materials with tunable properties; (i) glycerol/water mixtures (viscosity) and (ii) agarose gels (elasticity). Long-exposure white-light imaging is used to quantify the magnitude of the wave slope from which frequency-response diagrams are obtained via frequency sweeps. Resonance peaks and bandwidths are identified. These results show that for a given $n$ , the resonance frequency decreases with viscosity and increases with elasticity. The amplitude of the resonance peaks are much lower for gels and decrease further with mode number, indicating that much larger driving amplitudes are needed to overcome the elasticity and excite edge waves. The natural frequencies for a viscoelastic fluid in a cylindrical container with a pinned contact-line are computed from a theoretical model that depends upon the dimensionless Ohnesorge number ${Oh}$ , elastocapillary number ${Ec}$ and Bond number ${Bo}$ . All show good agreement with experimental observations. The eigenvalue problem is equivalent to the classic damped-driven oscillator model on linear operators with viscosity appearing as a damping force and elasticity and surface tension as restorative forces, consistent with our physical interpretation of these viscoelastic effects.more » « less
-
We investigate the modes of deformation of an initially spherical bubble immersed in a homogeneous and isotropic turbulent background flow. We perform direct numerical simulations of the two-phase incompressible Navier–Stokes equations, considering a low-density bubble in the high-density turbulent flow at various Weber numbers (the ratio of turbulent and surface tension forces) using the air–water density ratio. We discuss a theoretical framework for the bubble deformation in a turbulent flow using a spherical harmonic decomposition. We propose, for each mode of bubble deformation, a forcing term given by the statistics of velocity and pressure fluctuations, evaluated on a sphere of the same radius. This approach formally relates the bubble deformation and the background turbulent velocity fluctuations, in the limit of small deformations. The growth of the total surface deformation and of each individual mode is computed from the direct numerical simulations using an appropriate Voronoi decomposition of the bubble surface. We show that two successive temporal regimes occur: the first regime corresponds to deformations driven only by inertial forces, with the interface deformation growing linearly in time, in agreement with the model predictions, whereas the second regime results from a balance between inertial forces and surface tension. The transition time between the two regimes is given by the period of the first Rayleigh mode of bubble oscillation. We discuss how our approach can be used to relate the bubble lifetime to the turbulence statistics and eventually show that at high Weber numbers, bubble lifetime can be deduced from the statistics of turbulent fluctuations at the bubble scale.more » « less