Instability and dripping of electrified liquid films flowing down inverted substrates

We consider the gravity-driven flow of a perfect dielectric, viscous, thin liquid film, wetting a flat substrate inclined at a nonzero angle to the horizontal. The dynamics of the thin film is influenced by an electric field which is set up parallel to the substrate surface—this nonlocal physical mechanism has a linearly stabilizing effect on the interfacial dynamics. Our particular interest is in fluid films that are hanging from the underside of the substrate; these films may drip depending on physical parameters, and we investigate whether a sufficiently strong electric field can suppress such nonlinear phenomena. For a non-electrified flow, it was observed by Brun et al. [Phys. Fluids 27, 084107 (2015)] that the thresholds of linear absolute instability and dripping are reasonably close. In the present study, we incorporate an electric field and analyze the absolute and convective instabilities of a hierarchy of reduced-order models to predict the dripping limit in parameter space. The spatial stability results for the reduced-order models are verified by performing an impulse-response analysis with direct numerical simulations (DNS) of the Navier–Stokes equations coupled to the appropriate electrical equations. Guided by the results of the linear theory, we perform DNS on extended domains with inflow and outflow conditions (mimicking an experimental setup) to investigate the dripping limit for both non-electrified and electrified liquid films. For the latter, we find that the absolute instability threshold provides an order-of-magnitude estimate for the electric-field strength required to suppress dripping; the linear theory may thus be used to determine the feasibility of dripping suppression given a set of geometrical, fluid, and electrical parameters

Artificial viscosity model to mitigate numerical artefacts at fluid interfaces with surface tension

The numerical onset of parasitic and spurious artefacts in the vicinity of uid interfaces with surface tension is an important and well-recognised problem with respect to the accuracy and numerical stability of interfacial ow simulations. Issues of particular interest are spurious capillary waves, which are spatially underresolved by the computational mesh yet impose very restrictive time-step requirements, as well as parasitic currents, typically the result of a numerically unbalanced curvature evaluation. We present an arti cial viscosity model to mitigate numerical artefacts at surface-tension-dominated interfaces without adversely a ecting the accuracy of the physical solution. The proposed methodology computes an additional interfacial shear stress term, including an interface viscosity, based on the local ow data and uid properties that reduces the impact of numerical artefacts and dissipates underresolved small scale interface movements. Furthermore, the presented methodology can be readily applied to model surface shear viscosity, for instance to simulate the dissipative e ect of surface-active substances adsorbed at the interface. The presented analysis of numerical test cases demonstrates the e cacy of the proposed methodology in diminishing the adverse impact of parasitic and spurious interfacial artefacts on the convergence and stability of the numerical solution algorithm as well as on the overall accuracy of the simulation results

Stability of a bi-layer free film: simultaneous or individual rupture events?

We consider the stability of a long free film of liquid composed of two immiscible layers of differing viscosities, where each layer experiences a van der Waals force between its interfaces. We analyse the different ways the system can exhibit interfacial instability when the liquid layers are sufficiently thin. For an excess of surfactant on one gas–liquid interface the coupling between the layers is relatively weak and the instability manifests as temporally separated rupture events in each layer. Conversely, in the absence of surfactant the coupling between the layers is much stronger and the instability manifests as rupture of both layers simultaneously. These features are consistent with recent experimental observations

Electrohydrodynamic deformation and rotation of a particle-coated drop

A dielectric drop suspended in conducting liquid and subjected to an uniform electric field deforms into an ellipsoid whose major axis is either perpendicular or tilted (due to Quincke rotation effect) relative to the applied field. We experimentally study the effect of surface-adsorbed colloidal particles on these classic electrohydrodynamic phenomena. We observe that at high surface coverage (>90%), the electrohydrodynamic flow is suppressed, oblate drop deformation is enhanced, and the threshold for tilt is decreased compared to the particle-free drop. The deformation data are well explained by a capsule model, which assumes that the particle monolayer acts as an elastic interface. The reduction of the threshold field for rotation is likely related to drop asphericity

Localized Instabilities and Spinodal Decomposition in Driven Systems in the Presence of Elasticity

We study numerically and analytically the instabilities associated with phase separation in a solid layer on which an external material ux is imposed. The first instability is localized within a boundary layer at the exposed free surface by a process akin to spinodal decomposition. In the limiting static case, when there is no material ux, the coherent spinodal decomposition is recovered. In the present problem stability analysis of the time-dependent and non-uniform base states as well as numerical simulations of the full governing equations are used to establish the dependence of the wavelength and onset of the instability on parameter settings and its transient nature as the patterns eventually coarsen into a at moving front. The second instability is related to the Mullins- Sekerka instability in the presence of elasticity and arises at the moving front between the two phases when the ux is reversed. Stability analyses of the full model and the corresponding sharp-interface model are carried out and compared. Our results demonstrate how interface and bulk instabilities can be analysed within the same framework which allows to identify and distinguish each of them clearly. The relevance for a detailed understanding of both instabilities and their interconnections in a realistic setting are demonstrated for a system of equations modelling the lithiation/delithiation processes within the context of Lithium ion batteries.Comment: 8 figures, 19 page