Oscillations of a ring-constrained charged drop

Free drops of uncharged and charged inviscid, conducting fluids subjected to small-amplitude perturbations undergo linear oscillations (Rayleigh, Proc. R. Soc. London, vol. 29, no. 196–199, 1879, pp. 71–97; Rayleigh, Philos. Mag., vol. 14, no. 87, 1882, pp. 184–186). There exist a countably infinite number of oscillation modes, n=2,3,… role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3en=2,3,…n=2,3,…, each of which has a characteristic frequency and mode shape. Presence of charge (Q role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eQQ) lowers modal frequencies and leads to instability when Q\u3eQR role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eQ\u3eQRQ\u3eQR (Rayleigh limit). The n=0 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3en=0n=0 and n=1 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3en=1n=1 modes are disallowed because they violate volume conservation and cause centre of mass (COM) motion. Thus, the first mode to become unstable is the n=2 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3en=2n=2 prolate–oblate mode. For free drops, there is a one-to-one correspondence between mode number and shape (Legendre polynomial Pn role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3ePnPn). Recent research has shifted to studying oscillations of spherical drops constrained by solid rings. Pinning the drop introduces a new low-frequency mode of oscillation (n=1 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3en=1n=1), one associated primarily with COM translation of the constrained drop. We analyse theoretically the effect of charge on oscillations of constrained drops. Using normal modes and solving a linear operator eigenvalue problem, we determine the frequency of each oscillation mode. Results demonstrate that for ring-constrained charged drops (RCCDs), the association between mode number and shape is lost. For certain pinning locations, oscillations exhibit eigenvalue veering as Q role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eQQ increases. While slightly charged RCCDs pinned at zeros of P2 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eP2P2 have a first mode that involves COM motion and a second mode that entails prolate–oblate oscillations, the modes flip as Q role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eQQ increases. Thereafter, prolate–oblate oscillations of RCCDs adopt the role of being the first mode because they exhibit the lowest vibration frequency. At the Rayleigh limit, the first eigenmode – prolate–oblate oscillations – loses stability while the second – involving COM motion – remains stable

Pinch-off of a surfactant-covered jet

Surfactants at fluid interfaces not only lower and cause gradients in surface tension but can induce additional surface rheological effects in response to dilatational and shear deformations. Surface tension and surface viscosities are both functions of surfactant concentration. Measurement of surface tension and determination of its effects on interfacial flows are now well established. Measurement of surface viscosities, however, is notoriously difficult. Consequently, quantitative characterization of their effects in interfacial flows has proven challenging. One reason behind this difficulty is that, with most existing methods of measurement, it is often impossible to isolate the effects of surface viscous stresses from those due to Marangoni stresses. Here, a combined asymptotic and numerical analysis is presented of the pinch-off of a surfactant-covered Newtonian liquid jet. Similarity solutions obtained from slender-jet theory and numerical solutions are presented for jets with and without surface rheological effects. Near pinch-off, it is demonstrated that Marangoni stresses become negligible compared to other forces. The rate of jet thinning is shown to be significantly lowered by surface viscous effects. From analysis of the dynamics near the pinch-off singularity, a simple analytical formula is derived for inferring surface viscosities. Three-dimensional, axisymmetric simulations confirm the validity of the asymptotic analyses but also demonstrate that a thinning jet traverses a number of intermediate regimes before eventually entering the final asymptotic regime

Self-similar rupture of thin films of power-law fluids on a substrate

Thinning and rupture of a thin film of a power-law fluid on a solid substrate under the balance between destabilizing van der Waals pressure and stabilizing capillary pressure is analysed. In a power-law fluid, viscosity is not constant but is proportional to the deformation rate raised to the n−1 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3en−1n−1 power, where 00n=1n=1 for a Newtonian fluid). In the first part of the paper, use is made of the slenderness of the film and the lubrication approximation is applied to the equations of motion to derive a spatially one-dimensional nonlinear evolution equation for film thickness. The variation with time remaining until rupture of the film thickness, the lateral length scale, fluid velocity and viscosity is determined analytically and confirmed by numerical simulations for both line rupture and point rupture. The self-similarity of the numerically computed film profiles in the vicinity of the location where the film thickness is a minimum is demonstrated by rescaling of the transient profiles with the scales deduced from theory. It is then shown that, in contrast to films of Newtonian fluids undergoing rupture for which inertia is always negligible, inertia can become important during thinning of films of power-law fluids in certain situations. The critical conditions for which inertia becomes important and the lubrication approximation is no longer valid are determined analytically. In the second part of the paper, thinning and rupture of thin films of power-law fluids in situations when inertia is important are simulated by solving numerically the spatially two-dimensional, transient Cauchy momentum and continuity equations. It is shown that as such films continue to thin, a change of scaling occurs from a regime in which van der Waals, capillary and viscous forces are important to one where the dominant balance of forces is between van der Waals, capillary and inertial forces while viscous force is negligible

Inertial impedance of coalescence during collision of liquid drops

The fluid dynamics of the collision and coalescence of liquid drops has intrigued scientists and engineers for more than a century owing to its ubiquitousness in nature, e.g. raindrop coalescence, and industrial applications, e.g. breaking of emulsions in the oil and gas industry. The complexity of the underlying dynamics, which includes occurrence of hydrodynamic singularities, has required study of the problem at different scales – macroscopic, mesoscopic and molecular – using stochastic and deterministic methods. In this work, a multi-scale, deterministic method is adopted to simulate the approach, collision, and eventual coalescence of two drops where the drops as well as the ambient fluid are incompressible, Newtonian fluids. The free boundary problem governing the dynamics consists of the Navier–Stokes system and associated initial and boundary conditions that have been augmented to account for the effects of disjoining pressure as the separation between the drops becomes of the order of a few hundred nanometres. This free boundary problem is solved by a Galerkin finite element-based algorithm. The interplay of inertial, viscous, capillary and van der Waals forces on the coalescence dynamics is investigated. It is shown that, in certain situations, because of inertia two drops that are driven together can first bounce before ultimately coalescing. This bounce delays coalescence and can result in the computed value of the film drainage time departing significantly from that predicted from existing scaling theories

Dynamics of Thin Films Near Singularities Under the Influence of Non-Newtonian Rheology

Free surface flows where the shape of the interface separating two fluids is unknown apriori are an important area of interest in fluid dynamics. The study of free surface flows such as the breakup and coalescence of drops, and thinning and rupture of films lends itself to a diverse range of industrial applications, such as inkjet printing, crop spraying, foam and emulsion stability, and nanolithography, and helps develop an understanding of natural phenomena such as sea spray generation in oceans, or the dynamics of tear films in our eyes. In free surface flows, singularities are commonly observed in finite time, such as when the radius of a thread goes to zero upon pinchoff or when the thickness of a film becomes zero upon rupture. Dynamics in the vicinity of singularities usually lack a length scale and exhibit self-similarity. In such cases, universal scaling laws that govern the temporal behavior of measurable physical quantities such as the thickness of a film can be determined from asymptotic analysis and verified by high-resolution experiments and numerical simulations. These scaling laws provide deep insight into the underlying physics, and help delineate the regions of parameter space in which certain forces are dominant, while others are negligible. While the majority of previous works on singularities in free-surface flows deal with Newtonian fluids, many fluids in daily use and industry exhibit non-Newtonian rheology, such as polymer-laden, emulsion, foam, and suspension flows. The primary goal of this thesis is to investigate the thinning and rupture of thin films of non-Newtonian fluids exhibiting deformation-rate-thinning (power-law) rheology due to attractive intermolecular van der Waals forces. This is accomplished by means of intermediate asymptotic analysis and numerical simulations which utilize a robust Arbitrary Eulerian-Lagrangian (ALE) method that employs the Galerkin/Finite-Element Method for spatial discretization. For thinning of sheets of power-law fluids, a significant finding is the discovery of a previously undiscovered scaling regime where capillary, viscous and van der Waals forces due to attraction between the surfaces of the sheet, are in balance. For thinning of supported thin films, the breakdown of the lubrication approximation used almost exclusively in the past to study such systems, is shown to occur for films of power-law fluids through theory and confirmed by two dimensional simulations. The universality of scaling laws determined for rupture of supported films is shown by studying the impact of a bubble immersed in a power-law fluid with a solid wall. Emulsions, which are fine dispersions of drops of one liquid in another immiscible liquid, are commonly encountered in a variety of industries such as food, oil and gas, pharmaceuticals, and chemicals. Stability over a specified time frame is desirable in some applications, such as the shelf life of food products, while rapid separation into its constituent phases is required in others, such as when separating out brine from crude oil. The timescale over which coalescence of two drops of the dispersed phase occurs is crucial in determining emulsion stability

Local dynamics during thinning and rupture of liquid sheets of power-law fluids

Rupture of liquid sheets of power-law fluids surrounded by a gas is analysed under the competing influences of pressure due to van der Waals attraction, inertia, viscous stress and capillary pressure due to surface tension. Results of a combined theoretical and computational study are presented over the entire range of parameters governing the thinning of a power-law fluid of power-law exponent 00n=1n=1: Newtonian fluid) and Ohnesorge number 0≤Oh0≤Oh\u3c∞0≤OhOh≡μ0/ρh0σ−−−−√Oh≡μ0/ρh0σ, and μ0,ρ,h0 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eμ0,ρ,h0μ0,ρ,h0 and σ role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eσσ stand for the zero-deformation-rate viscosity, density, the initial sheet half-thickness and surface tension, respectively. The dynamics in the vicinity of the space–time singularity where the sheet ruptures is asymptotically self-similar, and thus the variation with time remaining until rupture τ≡tR−t role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eτ≡tR−tτ≡tR−t, where tR role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3etRtR is the time instant t role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3ett at which the sheet ruptures, of sheet half-thickness, lateral length scale and lateral velocity is determined analytically and confirmed by simulations. For sheets for which inertia is negligible (Oh−1=0 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eOh−1=0Oh−1=0), two distinct viscous scaling regimes are found, one for 0.580.58n≤0.58n≤0.58. The thinning dynamics of inviscid sheets (Oh=0 role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; border: 0px; font-variant: inherit; font-stretch: inherit; line-height: normal; font-family: inherit; vertical-align: baseline; display: inline-table; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; position: relative; \u3eOh=0Oh=0) is identical to that of Newtonian ones. For real fluids for which neither viscosity nor inertia is negligible, it is shown that the aforementioned creeping and inertial flow regimes are transitory and the thinning of power-law sheets exhibits a remarkably richer set of scaling transitions compared with Newtonian sheets