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

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