Search

The transition to dripping in the gravity-driven flow of a liquid film under an inclined plate is investigated at zero Reynolds number. Computations are carried out on a periodic domain assuming either a fixed fluid volume or a fixed flow rate for a hierarchy of models: two lubrication models with either linearised curvature or full curvature (the LCM and FCM, respectively), and the full equations of Stokes flow. Of particular interest is the breakdown of travelling-wave solutions as the plate inclination angle is increased. For any fixed volume the LCM reaches the horizontal state where it attains a cosine-shaped profile. For sufficiently small volume, the FCM and Stokes solutions attain a weak Young-Laplace equilibrium profile, the approach to which is described by an asymptotic analysis generalising that of Kalliadasis & Chang (1994) for the LCM. For large volumes, the bifurcation curves for the FCM and Stokes model have a turning point so that the fully inverted state is never reached. For fixed flow rate the LCM blows up at a critical angle that is well predicted by asymptotic analysis. The bifurcation curve for the FCM either has a turning point or else reaches a point at which the surface profile has an infinite slope singularity, indicating the onset of multi-valuedness. The latter is confirmed by the Stokes model which can be continued to obtain overturning surface profiles. Overall the thin-film models either provide an accurate prediction for dripping onset or else supply an upper bound on the critical inclination angle.

Locomotion is essential for living cells. It enables bacteria and algae to explore space for food, cancer to spread, and immune system to fight infections. Motile cells display trajectories of intriguing complexity, from regular (e.g. circular, helical, and so on) to irregular motions (run-tumble), the origin of which has remained elusive for over a century. This dynamics versatility is conventionally attributed to the shape asymmetry of the motile entity, to the suspending media, and/or to stochastic regulation. We propose here a universal approach highlighting that these movements are generic, occurring for a large class of cells and artificial microswimmers, without the need of invoking shape asymmetry nor stochasticity, but are encoded in their inherent nonlinear evolution. We show, in particular, that for a circular and spherical particle moving in a simple fluid, circular, helical and chaotic motions (akin to a persistent random walk) emerge naturally in different regions of parameter space. This establishes the operating principles for complex trajectories manifestation of motile systems, and offers a new vision with minimal ingredients. The reduced evolution equations based on symmetries are consistent with those derived for a model of an autophoretic particle including diffusion, emission/absorption at the particle surface and hydrodynamics, and provide qualitative and quantitative agreement., Comment: 11 pages, 3 figures

Autonomous locomotion is a ubiquitous phenomenon in biology and in physics of active systems at microscopic scale. This includes prokaryotic, eukaryotic cells (crawling and swimming) and artificial swimmers. An outstanding feature is the ability of these entities to follow complex trajectories, ranging from straight, curved (circular, helical...), to random-like ones. The non-straight nature of these trajectories is often explained as a consequence of the asymmetry of the particle or the medium in which it moves, or due to the presence of bounding walls, etc... Here, we show that straight, circular and helical trajectories emerge naturally in the absence of asymmetry of the swimmer or that of suspending medium. Our first proof is based on general considerations, without referring to an explicit form of a model. We show that these three trajectories correspond to self-congruent solutions. Self-congruency means that the states of the system at different moments of time can be made identical by an appropriate combination of rotation and translation of the coordinate space. We show that these solutions are exhibited by spherically symmetric particles as a result of a series of pitchfork bifurcations as the activity is increased. Self-congruent dynamics in one and two dimensions are analyzed as well. Finally, we present a simple explicit nonlinear exactly solvable model of fully isotropic phoretic particle that shows the transitions from a non-motile state to straight motion to circular motion to helical motion as a series of spontaneous symmetry-breaking bifurcations. Whether a system exhibits or not a given trajectory only depends on the numerical values of parameters entering the model, while asymmetry of swimmer shape, or anisotropy of the suspending medium , or influence of bounding walls are not necessary., Comment: 23 pages, 6 figures

We consider a 2D model of an autophoretic particle in which the particle has a circular shape and emits/absorbs a solute that diffuses and is advected by the suspending fluid. Beyond a certain emission/absorption rate (characterized by a dimensionless P\'eclet number, $Pe$) the particle is known to undergo a bifurcation from a non motile to a motile state, with different trajectories, going from a straight to circular and to a chaotic motion by progressively increasing $Pe$. From the full model involving solute diffusion and advection, we derive a reduced closed model which involves only two time-dependent amplitudes $C_1(t)$ and $C_2(t)$ corresponding to the first two Fourier modes of the solute concentration field. This model consists of two coupled nonlinear ordinary differential equations for $C_1$ and $C_2$ and presents several great advantages:(i) the straight and circular motions can be handled fully analytically, (ii) complex motions such as chaos can be analyzed numerically very efficiently in comparison to the numerically expensive full model involving partial differential equations, (iii) the reduced model has a universal form dictated only by symmetries, (iv) the model can be extended to higher Fourier modes. The derivation method is exemplified for a 2D model, for simplicity, but can easily be extended to 3D, not only for the presently selected phoretic model, but also for any model in which chemical activity triggers locomotion. A typical example can be found, for example, in the field of cell motility involving acto-myosin kinetics. This strategy offers an interesting way to cope with swimmers on the basis of ordinary differential equations, allowing for analytical tractability and efficient numerical treatment., Comment: 13 pages, 2 figures

We explore flow of a completely wetting fluid in a funnel, with particular focus on contact line instabilities at the fluid front. While the flow in a funnel may be related to a number of other flow configurations as limiting cases, understanding its stability is complicated due to the presence of additional azimuthal curvature, as well as due to convergent flow effects imposed by the geometry. Convergent nature of the flow leads to thickening of the film, therefore influencing its stability properties. In this work, we analyze these stability properties by combining physical experiments, asymptotic modeling, self-similar type of analysis and numerical simulations. We show that appropriate long-wave based model supported by the input from experiments, simulations and linear stability analysis origination from the flow down an incline plane provides a basic insight allowing to understand the development of contact line instability and emerging lengthscales.

We consider the Cahn-Hilliard (CH) equation with a Burgers-type convective term that is used as a model of coarsening dynamics in laterally driven phase-separating systems. In the absence of driving, it is known that solutions to the standard CH equation are characterized by an initial stage of phase separation into regions of one phase surrounded by the other phase (i.e., clusters or drops/holes or islands are obtained) followed by the coarsening process, where the average size of the structures grows in time and their number decreases. Moreover, two main coarsening modes have been identified in the literature, namely, coarsening due to volume transfer and due to translation. In the opposite limit of strong driving, the well-known Kuramoto-Sivashinsky (KS) equation is recovered, which may produce complicated chaotic spatio-temporal oscillations. The primary aim of the present work is to perform a detailed and systematic investigation of the transitions in the solutions of the convective CH (cCH) equation for a wide range of parameter values, and, in particular, to understand in detail how the coarsening dynamics is affected by an increase of the strength of the lateral driving force. Considering symmetric two-drop states, we find that one of the coarsening modes is stabilized at relatively weak driving, and the type of the remaining mode may change as driving increases. Furthermore, there exist intervals in the driving strength where coarsening is completely stabilized. In the intervals where the symmetric two-drop states are unstable they can evolve, for example, into one-drop states, two-drop states of broken symmetry or even time-periodic two-drop states that consist of two traveling drops that periodically exchange mass. We present detailed stability diagrams for symmetric two-drop states in various parameter planes and corroborate our findings by selected time simulations.

The flow of an electrified liquid film down an inclined plane wall is investigated with the focus on coherent structures in the form of travelling waves on the film surface, in particular, single-hump solitary pulses and their interactions. The flow structures are analysed first using a long-wave model, which is valid in the presence of weak inertia, and second using the Stokes equations. For obtuse angles, gravity is destablising and solitary pulses exist even in the absence of an electric field. For acute angles, spatially non-uniform solutions exist only beyond a critical value of the electric field strength; moreover, solitary-pulse solutions are present only at sufficiently high supercritical elec- tric field strengths. The electric field increases the amplitude of the pulses, can generate recirculation zones in the humps, and alters the far-field decay of the pulse tails from ex- ponential to algebraic with a significant impact on pulse interactions. A weak-interaction theory predicts an infinite sequence of bound-state solutions for non-electrified flow, and a finite set for electrified flow. The existence of single-hump pulse solutions and two-pulse bound states is confirmed for the Stokes equations via boundary-element computations. In addition, the electric field is shown to trigger a switch from absolute instability to convective instability, thereby regularising the dynamics, and this is confirmed by time- dependent simulations of the long-wave model.

We study contact line induced instabilities for a thin film of fluid under destabilizing gravitational force in three dimensional setting. In the previous work (Phys. Fluids, {\bf 22}, 052105 (2010)), we considered two dimensional flow, finding formation of surface waves whose properties within the implemented long wave model depend on a single parameter, $D=(3Ca)^{1/3}\cot\alpha$, where $Ca$ is the capillary number and $\alpha$ is the inclination angle. In the present work we consider fully 3D setting and discuss the influence of the additional dimension on stability properties of the flow. In particular, we concentrate on the coupling between the surface instability and the transverse (fingering) instabilities of the film front. We furthermore consider these instabilities in the setting where fluid viscosity varies in the transverse direction. It is found that the flow pattern strongly depends on the inclination angle and the viscosity gradient.