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
Plateau–Rayleigh instability in a soft viscoelastic material
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
- Award ID(s):
- 1750208
- PAR ID:
- 10219184
- Date Published:
- Journal Name:
- Soft Matter
- ISSN:
- 1744-683X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We explore basic mechanisms for the instability of finite-amplitude interfacial gravity waves through a two-dimensional linear stability analysis of the periodic and irrotational plane motion of the interface between two unbounded homogeneous fluids of different density in the absence of background currents. The flow domains are conformally mapped into two half-planes, where the time-varying interface is always mapped onto the real axis. This unsteady conformal mapping technique with a suitable representation of the interface reduces the linear stability problem to a generalized eigenvalue problem, and allows us to accurately compute the growth rates of unstable disturbances superimposed on steady waves for a wide range of parameters. Numerical results show that the wave-induced Kelvin–Helmholtz (KH) instability due to the tangential velocity jump across the interface produced by the steady waves is the major instability mechanism. Any disturbances whose dominant wavenumbers are greater than a critical value grow exponentially. This critical wavenumber that depends on the steady wave steepness and the density ratio can be approximated by a local KH stability analysis, where the spatial variation of the wave-induced currents is neglected. It is shown, however, that the growth rates need to be found numerically with care and the successive collisions of eigenvalues result in the wave-induced KH instability. In addition, the present study extends the previous results for the small-wavenumber instability, such as modulational instability, of relatively small-amplitude steady waves to finite-amplitude ones.more » « less
-
The Moore–Saffman–Tsai–Widnall (MSTW) instability is a parametric instability that arises in strained vortex columns. The strain is assumed to be weak and perpendicular to the vortex axis. In this second part of our investigation of vortex instability including density and surface tension effects, a linear stability analysis for this situation is presented. The instability is caused by resonance between two Kelvin waves with azimuthal wavenumber separated by two. The dispersion relation for Kelvin waves and resonant modes are obtained. Results show that the stationary resonant waves for $m=\pm 1$ are more unstable when the density ratio $\rho_2/\rho_1$ , the ratio of vortex to ambient fluid density, approaches zero, whereas the growth rate is maximised near $\rho _2/\rho _1 =0.215$ for the resonance $(m,m+2)=(0,2)$ . Surface tension suppresses the instability, but its effect is less significant than that of density. As the azimuthal wavenumber $m$ increases, the MSTW instability decays, in contrast to the curvature instability examined in Part 1 (Chang & Llewellyn Smith, J. Fluid Mech. vol. 913, 2021, A14).more » « less
-
null (Ed.)The Rayleigh–Taylor instability of a linear viscoelastic fluid overlying a passive gas is considered, where, under neutral conditions, the key dimensionless groups are the Bond number and the Weissenberg number. The branching behaviour upon instability to sinusoidal disturbances is determined by weak nonlinear analysis with the Bond number advanced from its critical value at neutral stability. It is shown that the solutions emanating from the critical state either branch out supercritically to steady waves at predictable wavelengths or break up subcritically with a wavelength having a single node. The nonlinear analysis leads to the counterintuitive observation that Rayleigh–Taylor instability of a viscoelastic fluid in a laterally unbounded layer must always result in saturated steady waves. The analysis also shows that the subcritical breakup in a viscoelastic fluid can only occur if the layer is laterally bounded below a critical horizontal width. If the special case of an infinitely deep viscoelastic layer is considered, a simple expression is obtained from which the transition between steady saturated waves and subcritical behaviour can be determined in terms of the leading dimensionless groups. This expression reveals that the supercritical saturation of the free surface is due to the influence of the normal elastic stresses, while the subcritical rupture of the free surface is attributed to the influence of capillary effects. In short, depending upon the magnitude of the scaled shear modulus, there exists a wavenumber at which a transition from saturated waves to subcritical breakup occurs.more » « less
-
Inward radial Rayleigh-Be'nard-Poiseuille flow can exhibit a buoyancy-driven instability when the Rayleigh number exceeds a critical value. Furthermore, similar to plane Rayleigh-Be'nard-Poiseuille flow, a viscous Tollmien-Schlichting instability can occur when the Reynolds number is high enough. Direct numerical simulations were carried out with a compressible Navier-Stokes code in cylindrical coordinates to investigate the spatial stability of the inward radial flow inside the collector of a hypothetical solar chimney power plant. The convective terms were discretized with fifth-order-accurate upwind-biased compact finite-differences and the viscous terms were discretized with fourth-order-accurate compact finite differences. For cases with buoyancy-driven instability, steady three-dimensional waves are strongly amplified. The spatial growth rates vary significantly in the radial direction and lower azimuthal mode numbers are amplified closer to the outflow. Traveling oblique modes are amplified as well. The growth rates of the oblique modes decrease with increasing frequency. In addition to the purely radial flow, a spiral flow with swept inflow was examined. Overall lower growth rates are observed for the spiral flow compared to the radial flow. Different from the radial flow, the relative wave angles and growth rates of the left and right traveling oblique modes are not identical. A plane RBP case with viscosity-driven instability by Chung et al. was considered as well. The reported growth rate and phase speed were matched with good accuracy.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D1SM00019E
