Linear stability analysis of electrohydrodynamic instabilities at fluid interfaces in the “small feature” limit

Abstract

The endeavour to effectively harness interfacial electrohydrodynamic instabilities, to create small patterns, involves reducing the wavelength of the instability. This can be accomplished by decreasing the separation between the electrodes which may not always be possible. One may therefore have to reduce the surface tension or increase the applied voltage at a fixed electrode spacing. This can result in the wavelength of the pattern becoming of the same order as the electrode separation. Pease and Russel (J. Chem. Phys. 118, 3790 (2003)) were the first to argue that the commonly used Thin-Film Approximation (TFA) that involves an asymptotic expansion in the small parameter δ = (ε(0) φ(0)(2)/(γh(0)))(1/2) (where ε (0) is the permittivity of vacuum, φ (0) is the root mean square value of the applied potential, γ is the surface tension and h (0) is the thickness of the thin film) need not always be valid and γ may not be small in experiments. Higher-order corrections to the TFA might therefore be necessary. We extend the Direct Current (DC) field analysis of Pease and Russel to an Alternating Current (AC) field. AC field has been suggested as an effective way of controlling the wavelength of electrohydrodynamic instabilities at fluid-fluid interfaces. Infact, the perfect and leaky dielectric limits can be realised in the same fluid at very high and very low electric field frequencies, respectively. Recently, Roberts and Kumar (J. Fluid Mech. 631, 255 (2009)) carried out an analysis using TFA to investigate AC-field-induced instabilities at air-polymer interfaces. We propose a Generalized Model (GM), without the lubrication approximation, and carry out detailed comparison with the TFA. We consider the top fluid to be air, a perfect dielectric, and the bottom fluid to be a perfect or a leaky dielectric. The analysis is carried out for both DC and AC fields, and the deviation from TFA is expressed in terms of the parameter B = γh(0)/(ε(0) φ(0)(2) ) = δ(-2). We discuss variation of the wavelength of the fastest growing mode with frequency of the applied field for any arbitrary value of B, unlike the analysis of Roberts and Kumar which is restricted to B ≫ 1(δ ≪ 1) . We also revisit the analysis of Pease and Russel for instabilities under DC field and present the results in terms of the single parameter, B.