Numerical Investigation of the Steady State of a Driven Thin Film Equation

A third-order ordinary differential equation with application in the flow of a thin liquid film is considered. The boundary conditions come from Tanner's problem for the surface tension driven flow of a thin film. Symmetric and nonsymmetric finite difference schemes are implemented in order to obtain steady state solutions. We show that a central difference approximation to the third derivative in the model equation produces a solution curve with oscillations. A difference scheme based on a combination of forward and backward differences produces a smooth accurate solution curve. The stability of these schemes is analysed through the use of a von Neumann stability analysis

Instabilities in free-surface electroosmotic flows

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.With the recent development of novel microfluidic devices electroosmotic flows with fluid/fluid interfaces have emerged as very important subjects of investigation. Two immiscible fluids may need to be transported in a microchannel, or one side of a channel may be open to air for various purposes, including adsorption of airborne molecules to liquid for high-sensitivity substance detection. The liquid/liquid or liquid/gas interface in these cases can deform, resulting in significant corrugations followed sometimes by incipient rupture of liquid layers. For electroosmotic flow the rupture, leading to shortcircuit, can cause overall failure of the device. It is thus imperative to know the conditions for the rupture as well as the initial interfacial instability. Studies based on the Debye-Huckle approximation reveal that all free-surface electroosmotic flows of thickness larger than the Debye screening length are unstable and selectively lead to rupture. Layers of the order of Debye screening length, however, are not properly described by the Debye-Huckle approximation. Even for micro-scale layers, the rupture phenomenon can make local layer thickness to be nanoscale. A fully coupled system of hydrodynamics, electric field, and ionic distribution need to be analyzed. In this paper linear instability and subsequent nonlinear developments of a nanoscale free-surface electroosmotic flow are reported.This study is sponsored by the Ministry of Education, Science and Technology of Korea through the World Class University Grant

Effect of disjoining pressure in a thin film equation with\ud non-uniform forcing

We explore the effect of disjoining pressure on a thin film equation in the presence of a non-uniform body force, motivated by a model describing the reverse draining of a magnetic film. To this end, we use a combination of numerical investigations and analytical considerations. The disjoining pressure has a regularizing influence on the evolution of the system and appears to select a single steady-state solution for fixed height boundary conditions; this is in contrast with the existence of a continuum of locally attracting solutions that exist in the absence of disjoining pressure for the same boundary conditions. We numerically implement matched asymptotics expansions to construct equilibrium solutions and also investigate how they behave as the disjoining pressure is sent to zero. Finally, we consider the effect of the competition between forcing and disjoining pressure on the coarsening dynamics of the thin film for fixed contact angle boundary conditions

Time integration and steady-state continuation for 2d lubrication equations

Lubrication equations allow to describe many structurin processes of thin liquid films. We develop and apply numerical tools suitable for their analysis employing a dynamical systems approach. In particular, we present a time integration algorithm based on exponential propagation and an algorithm for steady-state continuation. In both algorithms a Cayley transform is employed to overcome numerical problems resulting from scale separation in space and time. An adaptive time-step allows to study the dynamics close to hetero- or homoclinic connections. The developed framework is employed on the one hand to analyse different phases of the dewetting of a liquid film on a horizontal homogeneous substrate. On the other hand, we consider the depinning of drops pinned by a wettability defect. Time-stepping and path-following are used in both cases to analyse steady-state solutions and their bifurcations as well as dynamic processes on short and long time-scales. Both examples are treated for two- and three-dimensional physical settings and prove that the developed algorithms are reliable and efficient for 1d and 2d lubrication equations, respectively.Comment: 33 pages, 16 figure

Bifurcation analysis of the behavior of partially wetting liquids on a rotating cylinder

We discuss the behavior of partially wetting liquids on a rotating cylinder using a model that takes into account the effects of gravity, viscosity, rotation, surface tension and wettability. Such a system can be considered as a prototype for many other systems where the interplay of spatial heterogeneity and a lateral driving force in the proximity of a first- or second-order phase transition results in intricate behavior. So does a partially wetting drop on a rotating cylinder undergo a depinning transition as the rotation speed is increased, whereas for ideally wetting liquids the behavior \bfuwe{only changes quantitatively. We analyze the bifurcations that occur when the rotation speed is increased for several values of the equilibrium contact angle of the partially wetting liquids. This allows us to discuss how the entire bifurcation structure and the flow behavior it encodes changes with changing wettability. We employ various numerical continuation techniques that allow us to track stable/unstable steady and time-periodic film and drop thickness profiles. We support our findings by time-dependent numerical simulations and asymptotic analyses of steady and time-periodic profiles for large rotation numbers

Generalized linear stability of noninertial coating flows over topographical features

The transient evolution of perturbations to steady lubrication flow over a topographically patterned surface is investigated via a nonmodal linear stability analysis of the non-normal disturbance operator. In contrast to the capillary ridges that form near moving contact lines, the stationary capillary ridges near trenches or elevations have only stable eigenvalues. Minimal transient amplification of perturbations occurs, regardless of the magnitude or steepness of the topographical features. The absence of transient amplification and the stability of the ridge are explained on physical grounds. By comparison to unstable ridge formation on smooth, flat, and homogeneous surfaces, the lack of closed, recirculating streamlines beneath the capillary ridge is linked to the linear stability