BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Uwe Thiele

Uwe Thiele appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softphysics.flu-dynnlin.PSphysics.chem-phcond-mat.stat-mechBiological Physicscond-mat.mes-hallcond-mat.mtrl-scicond-mat.othernlin.CDpatt-sol

Catalog footprint

What is connected

41works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2025arXiv

Dynamic Models for Two Nonreciprocally Coupled Fields: A Microscopic Derivation for Zero, One, and Two Conservation Laws

We construct dynamic models governing two nonreciprocally coupled fields for several cases with zero, one, and two conservation laws. Starting from two microscopic nonreciprocally coupled Ising models, and using the mean-field approximation, we obtain closed-form evolution equations for the spatially resolved magnetization in each lattice. Only allowing for single spin-flip dynamics, the macroscopic equations in the thermodynamic limit are closely related to the nonreciprocal Allen-Cahn equations, i.e. conservation laws are absent. Likewise, only accounting for spin-exchange dynamics within each lattice, the thermodynamic limit yields equations similar to the nonreciprocal Cahn-Hilliard model, i.e. with two conservation laws. In the case of spin-exchange dynamics within and between the two lattices, we obtain two nonreciprocally coupled equations that add up to one conservation law. For each of these cases, we systematically map out the linear instabilities that can arise. Moreover, combining the different dynamics gives a large number of further models. Our results provide a microscopic foundation for a broad class of nonreciprocal field theories, establishing a direct link between nonequilibrium statistical mechanics and macroscopic continuum descriptions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bifurcations of front motion in passive and active Allen-Cahn-type equations

The well-known cubic Allen-Cahn (AC) equation is a simple gradient dynamics (or variational) model for a nonconserved order parameter field. After revising main literature results for the occuring different types of moving fronts, we employ path continuation to determine their bifurcation diagram in dependence of the external field strength or chemical potential. We then employ the same methodology to systematically analyse fronts for more involved AC-type models. In particular, we consider a cubic-quintic variational AC model and two different nonvariational generalisations. We determine and compare the bifurcation diagrams of front solutions in the four considered models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient calculation of phase coexistence and phase diagrams: application to a binary phase-field crystal model

We show that one can employ well-established numerical continuation methods to efficiently calculate the phase diagram for thermodynamic systems. In particular, this involves the determination of lines of phase coexistence related to first order phase transitions and the continuation of triple points. To illustrate the method we apply it to a binary Phase-Field-Crystal model for the crystallisation of a mixture of two types of particles. The resulting phase diagram is determined for one- and two-dimensional domains. In the former case it is compared to the diagram obtained from a one-mode approximation. The various observed liquid and crystalline phases and their stable and metastable coexistence are discussed as well as the temperature-dependence of the phase diagrams. This includes the (dis)appearance of critical points and triple points. We also relate bifurcation diagrams for finite-size systems to the thermodynamics of phase transitions in the infinite-size limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

First order phase transitions and the thermodynamic limit

We consider simple mean field continuum models for first order liquid-liquid demixing and solid-liquid phase transitions and show how the Maxwell construction at phase coexistence emerges on going from finite-size closed systems to the thermodynamic limit. The theories considered are the Cahn-Hilliard model of phase separation, which is also a model for the liquid-gas transition, and the phase field crystal model of the solid-liquid transition. Our results show that states comprising the Maxwell line depend strongly on the mean density with spatially localized structures playing a key role in the approach to the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase-field crystal description of active crystallites: elastic and inelastic collisions

The active Phase-Field-Crystal (aPFC) model combines elements of the Toner-Tu theory for self-propelled particles and the classical Phase-Field-Crystal (PFC) model that describes the transition between liquid to crystalline phases. In the liquid-crystal coexistence region of the PFC model, crystalline clusters exist in the form of localized states that coexist with the homogeneous background. At sufficiently strong self-propulsion strength (or activity) they start to travel. We employ numerical path continuation and direct time simulations to first investigate the existence regions of different types of localized states for a one-dimensional system. The results are summarized in morphological phase diagrams in the parameter plane spanned by activity and mean concentration. Then we focus on the interaction of traveling localized states studying their collision behavior. As a result we distinguish 'elastic' and 'inelastic' collisions. In the former, localized states recover their properties after a collision while in the latter they may annihilate or form resting or traveling bound states. In passing, we describe oscillating and modulated traveling localized states that have as the steadily traveling localized states no counterpart in the classical PFC model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-dimensional localized states in an active phase-field-crystal model

The active phase-field-crystal (active PFC) model provides a simple microscopic mean field description of crystallization in active systems. It combines the PFC model (or conserved Swift-Hohenberg equation) of colloidal crystallization and aspects of the Toner-Tu theory for self-propelled particles. We employ the active PFC model to study the occurrence of localized and periodic active crystals in two spatial dimensions. Due to the activity, crystalline states can undergo a drift instability and start to travel while keeping their spatial structure. Based on linear stability analyses, time simulations and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of resting and traveling states. We explore, for instance, how the slanted homoclinic snaking of steady localized states found for the passive PFC model is modified by activity. The analysis is carried out for the model in two spatial dimensions. Morphological phase diagrams showing the regions of existence of various solution types are presented merging the results from all the analysis tools employed. We also study how activity influences the crystal structure with transitions from hexagons to rhombic and stripe patterns. This in-depth analysis of a simple PFC model for active crystals and swarm formation provides a clear general understanding of the observed multistability and associated hysteresis effects, and identifies thresholds for qualitative changes in behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-dimensional patterns in dip coating -- first steps on the continuation path

We present a brief comparative investigation of the bifurcation structure related to the formation of two-dimensional deposition patterns as described by continuum models of Cahn-Hilliard type. These are, on the one hand a driven Cahn-Hilliard model for Langmuir-Blodgett transfer of a surfactant layer from the surface of a bath onto a moving plate and on the other hand a driven thin-film equation modelling the surface acoustic wave-driven coating of a plate by a simple liquid. In both cases, we present selected two-dimensional steady states corresponding to deposition patterns and discuss the main structure of the corresponding bifurcation diagrams.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dip-coating with prestructured substrates: transfer of simple liquids and Langmuir-Blodgett monolayers

When a plate is withdrawn from a liquid bath, either a static meniscus forms in the transition region between the bath and the substrate or a liquid film of finite thickness (a Landau-Levich film) is transferred onto the moving substrate. If the substrate is inhomogeneous, e.g., has a prestructure consisting of stripes of different wettabilities, the meniscus can be deformed or show a complex dynamic behavior. Here we study the free surface shape and dynamics of a dragged meniscus occurring for striped prestructures with two orientations, parallel and perpendicular to the transfer direction. A thin film model is employed that accounts for capillarity through a Laplace pressure and for the spatially varying wettability through a Derjaguin (or disjoining) pressure. Numerical continuation is used to obtain steady free surface profiles and corresponding bifurcation diagrams in the case of substrates with different homogeneous wettabilities. Direct numerical simulations are employed in the case of the various striped prestructures. The final part illustrates the importance of our findings for particular applications that involve complex liquids by modeling a Langmuir-Blodgett transfer experiment. There, one transfers a monolayer of an insoluble surfactant that covers the surface of the bath onto the moving substrate. The resulting pattern formation phenomena can be crucially influenced by the hydrodynamics of the liquid meniscus that itself depends on the prestructure on the substrate. In particular, we show how prestructure stripes parallel to the transfer direction lead to the formation of bent stripes in the surfactant coverage after transfer and present similar experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From a thin film model for passive suspensions towards the description of osmotic biofilm spreading

Biofilms are ubiquitous macro-colonies of bacteria that develop at various interfaces (solid-liquid, solid-gas or liquid-gas). The formation of biofilms starts with the attachment of individual bacteria to an interface, where they proliferate and produce a slimy polymeric matrix - two processes that result in colony growth and spreading. Recent experiments on the growth of biofilms on agar substrates under air have shown that for certain bacterial strains, the production of the extracellular matrix and the resulting osmotic influx of nutrient-rich water from the agar into the biofilm are more crucial for the spreading behaviour of a biofilm than the motility of individual bacteria. We present a model which describes the biofilm evolution and the advancing biofilm edge for this spreading mechanism. The model is based on a gradient dynamics formulation for thin films of biologically passive liquid mixtures and suspensions, supplemented by bioactive processes which play a decisive role in the osmotic spreading of biofilms. It explicitly includes the wetting properties of the biofilm on the agar substrate via a disjoining pressure and can therefore give insight into the interplay between passive surface forces and bioactive growth processes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gradient dynamics models for liquid films with soluble surfactant

In this paper we propose equations of motion for the dynamics of liquid films of surfactant suspensions that consist of a general gradient dynamics framework based on an underlying energy functional. This extends the gradient dynamics approach to dissipative non-equilibrium thin film systems with several variables, and casts their dynamic equations into a form that reproduces Onsager's reciprocity relations. We first discuss the general form of gradient dynamics models for an arbitrary number of fields and discuss simple well-known examples with one or two fields. Next, we develop the gradient dynamics (three field) model for a thin liquid film covered by soluble surfactant and discuss how it automatically results in consistent convective (driven by pressure gradients, Marangoni forces and Korteweg stresses), diffusive, adsorption/desorption, and evaporation fluxes. We then show that in the dilute limit, the model reduces to the well-known hydrodynamic form that includes Marangoni fluxes due to a linear equation of state. In this case the energy functional incorporates wetting energy, surface energy of the free interface (constant contribution plus an entropic term) and bulk mixing entropy. Subsequently, as an example, we show how various extensions of the energy functional result in consistent dynamical models that account for nonlinear equations of state, concentration-dependent wettability and surfactant and film bulk decomposition phase transitions. We conclude with a discussion of further possible extensions towards systems with micelles, surfactant adsorption at the solid substrate and bioactive behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Morphological transitions of sliding drops: Dynamics and bifurcations

We study fully three-dimensional droplets that slide down an incline by employing a thin-film equation that accounts for capillarity, wettability, and a lateral driving force in small-gradient (or long-wave) approximation. In particular, we focus on qualitative changes in the morphology and behavior of stationary sliding drops. We employ the inclination angle of the substrate as control parameter and use continuation techniques to analyze for several fixed droplet sizes the bifurcation diagram of stationary droplets, their linear stability, and relevant eigenmodes. The obtained predictions on existence ranges and instabilities are tested via direct numerical simulations that are also used to investigate a branch of time-periodic behavior (corresponding to repeated breakup-coalescence cycles, where the breakup is also denoted as pearling) which emerges at a global instability, the related hysteresis in behavior, and a period-doubling cascade. The nontrivial oscillatory behavior close to a Hopf bifurcation of drops with a finite-length tail is also studied. Finally, it is shown that the main features of the bifurcation diagram follow scaling laws over several decades of the droplet size.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Liquid drops on a surface: using density functional theory to calculate the binding potential and drop profiles and comparing with results from mesoscopic modelling

The contribution to the free energy for a film of liquid of thickness $h$ on a solid surface, due to the interactions between the solid-liquid and liquid-gas interfaces is given by the binding potential, $g(h)$. The precise form of $g(h)$ determines whether or not the liquid wets the surface. Note that differentiating $g(h)$ gives the Derjaguin or disjoining pressure. We develop a microscopic density functional theory (DFT) based method for calculating $g(h)$, allowing us to relate the form of $g(h)$ to the nature of the molecular interactions in the system. We present results based on using a simple lattice gas model, to demonstrate the procedure. In order to describe the static and dynamic behaviour of non-uniform liquid films and drops on surfaces, a mesoscopic free energy based on $g(h)$ is often used. We calculate such equilibrium film height profiles and also directly calculate using DFT the corresponding density profiles for liquid drops on surfaces. Comparing quantities such as the contact angle and also the shape of the drops, we find good agreement between the two methods. We also study in detail the effect on $g(h)$ of truncating the range of the dispersion forces, both those between the fluid molecules and those between the fluid and wall. We find that truncating can have a significant effect on $g(h)$ and the associated wetting behaviour of the fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modelling Pattern Formation in Dip-Coating Experiments

We briefly review selected mathematical models that describe the dynamics of pattern formation phenomena in dip-coating and Langmuir-Blodgett transfer experiments, where solutions or suspensions are transferred onto a substrate producing patterned deposit layers with structure length from hundreds of nanometres to tens of micrometres. The models are presented with a focus on their gradient dynamics formulations that clearly shows how the dynamics is governed by particular free energy functionals and facilitates the comparison of the models. In particular, we include a discussion of models based on long-wave hydrodynamics as well as of more phenomenological models that focus on the pattern formation processes in such systems. The models and their relations are elucidated and examples of resulting patterns are discussed before we conclude with a discussion of implications of the gradient dynamics formulation and of some related open issues.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coarsening modes of clusters of aggregating particles

There are two modes by which clusters of aggregating particles can coalesce: The clusters can merge either (i) by the Ostwald ripening process in which particles diffuse from one cluster to the other whilst the cluster centres remain stationary, or (ii) by means of a cluster translation mode, in which the clusters move towards each other and join. To understand in detail the interplay between these different modes, we study a model system of hard particles with an additional attraction between them. The particles diffuse along narrow channels with smooth or periodically corrugated walls, so that the system may be treated as one-dimensional. When the attraction between the particles is strong enough, they aggregate to form clusters. The channel potential influences whether clusters can move easily or not through the system and can prevent cluster motion. We use Dynamical Density Functional theory to study the dynamics of the aggregation process, focusing in particular on the coalescence of two equal size clusters. As long as the particle hard-core diameter is non-zero, we find that the coalescence process can be halted by a sufficiently strong corrugation potential. The period of the potential determines the size of the final stable clusters. For the case of smooth channel walls, we demonstrate that there is a cross-over in the dominance of the two different coarsening modes, that depends on the strength of the attraction between particles, the cluster sizes and the separation distance between clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collapsed heteroclinic snaking near a heteroclinic chain in dragged meniscus problems

We study a liquid film that is deposited onto a flat plate that is inclined at a constant angle to the horizontal and is extracted from a liquid bath at a constant speed. We additionally assume that there is a constant temperature gradient along the plate that induces a Marangoni shear stress. We analyse steady-state solutions of a long-wave evolution equation for the film thickness. Using centre manifold theory, we first obtain an asymptotic expansion of solutions in the bath region. The presence of the temperature gradient significantly changes these expansions and leads to the presence of logarithmic terms that are absent otherwise. Next, we obtain numerical solutions of the steady-state equation and analyse the behaviour of the solutions as the plate velocity is changed. We observe that the bifurcation curve exhibits snaking behaviour when the plate inclination angle is beyond a certain critical value. Otherwise, the bifurcation curve is monotonic. The solutions along these curves are characterised by a foot-like structure that is formed close to the meniscus and is preceded by a thin precursor film further up the plate. The length of the foot increases along the bifurcation curve. Finally, we explain that the snaking behaviour of the bifurcation curves is caused by the existence of an infinite number of heteroclinic orbits close to a heteroclinic chain that connects in an appropriate three-dimensional phase space the fixed point corresponding to the precursor film with the fixed point corresponding to the foot and then with the fixed point corresponding to the bath.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Equilibrium Gibbs' Criterion for Completely Wetting Volatile Liquids

During the spreading of a liquid over a solid substrate, the contact line can stay pinned at sharp edges until the contact angle exceeds a critical value. At (or sufficiently near) equilibrium, this is known as Gibbs' criterion. Here, we show both experimentally and theoretically that for completely wetting volatile liquids there also exists a dynamically-produced critical angle for depinning, which increases with the evaporation rate. This suggests that one may introduce a simple modification of the Gibbs' criterion for (de)pinning, that accounts for the non-equilibrium effect of evaporation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stability of liquid films covered by a carpet of self-propelled surfactant particles

We consider a carpet of self-propelled particles at the liquid-gas interface of a liquid film on a solid substrate. The particles excert an excess pressure on the interface and also move along the interface while the swimming direction performs rotational diffusion. We study the intricate influence of these self-propelled insoluble surfactants on the stability of the film surface and show that depending on the strength of in-surface rotational diffusion and the absolute value of the in-surface swimming velocity several characteristic instability modes can occur. In particular, rotational diffusion can either stabilize the film or induce instabilities of different character.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Introduction to Microscopic Theories for Inhomogeneous Liquids: Getting Started with Density Functional Theory

Classical density functional theory (DFT) is a statistical mechanical theory for calculating the density profiles of the molecules in a liquid. It is widely used, for example. to calculate the density distribution of the molecules in the vicinity of a confining wall, the interfacial tension, the wetting behaviour and many other properties of nonuniform liquids. DFT can however be somewhat daunting to students entering the field, because of the many connections to other areas of liquid-state science that are required and used to develop the theories. Here we give an introduction to some of the key ideas, based on a lattice-gas (Ising) model fluid. This builds on knowledge covered in most undergraduate statistical mechanics and thermodynamics courses and so students can quickly get to the stage of calculating density profiles, etc for themselves. We derive a simple DFT for the lattice-gas and present some typical results that can readily be calculated using the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gradient dynamics description for films of mixtures and suspensions - the case of dewetting triggered by coupled film height and concentration fluctuations

A gradient dynamics model based on an extended interface Hamiltonian is presented that is able to describe the dynamics of structuring processes in thin films of liquid mixtures, solutions and suspensions on solid substrates including coupled dewetting and decomposition. After discussing known limiting cases the model is employed to investigate the dewetting of thin films of liquid mixtures and suspensions under the influence of effective long-range van der Waals forces that depend on solute concentration. The occurring fluxes are discussed and it is shown that spinodal dewetting may be triggered through the coupling of film height and concentration fluctuations. Fully nonlinear calculations provide the time evolution and resulting steady film height and concentration profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Localized states in the conserved Swift-Hohenberg equation with cubic nonlinearity

The conserved Swift-Hohenberg equation with cubic nonlinearity provides the simplest microscopic description of the thermodynamic transition from a fluid state to a crystalline state. The resulting phase field crystal model describes a variety of spatially localized structures, in addition to different spatially extended periodic structures. The location of these structures in the temperature versus mean order parameter plane is determined using a combination of numerical continuation in one dimension and direct numerical simulation in two and three dimensions. Localized states are found in the region of thermodynamic coexistence between the homogeneous and structured phases, and may lie outside of the binodal for these states. The results are related to the phenomenon of slanted snaking but take the form of standard homoclinic snaking when the mean order parameter is plotted as a function of the chemical potential, and are expected to carry over to related models with a conserved order parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modelling spreading dynamics of liquid crystals in three spatial dimensions

We study spreading dynamics of nematic liquid crystal droplets within the framework of the long-wave approximation. A fourth order nonlinear parabolic partial differential equation governing the free surface evolution is derived. The influence of elastic distortion energy and of imposed anchoring variations at the substrate are explored through linear stability analysis and scaling arguments, which yield useful insight and predictions for the behaviour of spreading droplets. This behaviour is captured by fully nonlinear time-dependent simulations of three dimensional droplets spreading in the presence of anchoring variations that model simple defects in the nematic orientation at the substrate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Note on the hydrodynamic description of thin nematic films: strong anchoring model

We discuss the long-wave hydrodynamic model for a thin film of nematic liquid crystal in the limit of strong anchoring at the free surface and at the substrate. We rigorously clarify how the elastic energy enters the evolution equation for the film thickness in order to provide a solid basis for further investigation: several conflicting models exist in the literature that predict qualitatively different behaviour. We consolidate the various approaches and show that the long-wave model derived through an asymptotic expansion of the full nemato-hydrodynamic equations with consistent boundary conditions agrees with the model one obtains by employing a thermodynamically motivated gradient dynamics formulation based on an underlying free energy functional. As a result, we find that in the case of strong anchoring the elastic distortion energy is always stabilising. To support the discussion in the main part of the paper, an appendix gives the full derivation of the evolution equation for the film thickness via asymptotic expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Parameter passing between Molecular Dynamics and continuum models for droplets on solid substrates - The static case

We study equilibrium properties of polymer films and droplets on a solid substrate employing particle-based simulation techniques (Molecular Dynamics) and a continuum description. Parameter-passing techniques are explored that facilitate a detailed comparison of the two models. In particular, the liquid-vapor, solid-liquid and solid-vapor interface tensions, and the Derjaguin or disjoining pressure are determined by Molecular Dynamics simulations. This information is then introduced into continuum descriptions accounting for (i) the full curvature and (ii) a long-wave approximation of the curvature (thin film model). A comparison of the dependence of the contact angle on droplet size indicates that the theories agree well if the contact angles are defined in a compatible manner.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Patterned deposition at moving contact lines

When a simple or complex liquid recedes from a smooth solid substrate it often leaves a homogeneous or structured deposit behind. In the case of a receding non-volatile pure liquid the deposit might be a liquid film or an arrangement of droplets depending on the receding speed of the meniscus and the wetting properties of the system. For complex liquids with volatile components as, e.g., polymer solutions and particle or surfactant suspensions, the deposit might be a homogeneous or structured layer of solute - with structures ranging from line patterns that can be orthogonal or parallel to the receding contact line via hexagonal or square arrangements of drops to complicated hierarchical structures. We review a number of recent experiments and modelling approaches with a particular focus on mesoscopic hydrodynamic long-wave models. The conclusion highlights open question and speculates about future developments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modelling the formation of structured deposits at receding contact lines of evaporating solutions and suspensions

When a film of a liquid suspension of nanoparticles or a polymer solution is deposited on a surface, it may dewet from the surface and as the solvent evaporates the solute particles/polymer can be deposited on the surface in regular line patterns. In this paper we explore a hydrodynamic model for the process that is based on a long-wave approximation that predicts the deposition of irregular and regular line patterns. This is due to a self-organised pinning-depinning cycle that resembles a stick-slip motion of the contact line. We present a detailed analysis of how the line pattern properties depend on quantities such as the evaporation rate, the solute concentration, the Péclet number, the chemical potential of the ambient vapour, the disjoining pressure, and the intrinsic viscosity. The results are related to several experiments and to depinning transitions in other soft matter systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The relation of steady evaporating drops fed by an influx and freely evaporating drops

We discuss a thin film evolution equation for a wetting evaporating liquid on a smooth solid substrate. The model is valid for slowly evaporating small sessile droplets when thermal effects are insignificant, while wettability and capillarity play a major role. The model is first employed to study steady evaporating drops that are fed locally through the substrate. An asymptotic analysis focuses on the precursor film and the transition region towards the bulk drop and a numerical continuation of steady drops determines their fully non-linear profiles. Following this, we study the time evolution of freely evaporating drops without influx for several initial drop shapes. As a result we find that drops initially spread if their initial contact angle is larger than the apparent contact angle of large steady evaporating drops with influx. Otherwise they recede right from the beginning.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermodynamically consistent description of the hydrodynamics of free surfaces covered by insoluble surfactants of high concentration

In this paper we propose several models that describe the dynamics of liquid films which are covered by a high concentration layer of insoluble surfactant. First, we briefly review the 'classical' hydrodynamic form of the coupled evolution equations for the film height and surfactant concentration that are well established for small concentrations. Then we re-formulate the basic model as a gradient dynamics based on an underlying free energy functional that accounts for wettability and capillarity. Based on this re-formulation in the framework of nonequilibrium thermodynamics, we propose extensions of the basic hydrodynamic model that account for (i) nonlinear equations of state, (ii) surfactant-dependent wettability, (iii) surfactant phase transitions, and (iv) substrate-mediated condensation. In passing, we discuss important differences to most of the models found in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modelling fluids and crystals using a two-component modified phase field crystal model

A modified phase field crystal model in which the free energy may be minimised by an order parameter profile having isolated bumps is investigated. The phase diagram is calculated in one and two dimensions and we locate the regions where modulated and uniform phases are formed and also regions where localised states are formed. We investigate the effectiveness of the phase field crystal model for describing fluids and crystals with defects. We further consider a two component model and elucidate how the structure transforms from hexagonal crystalline ordering to square ordering as the concentration changes. Our conclusion contains a discussion of possible interpretations of the order parameter field.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ratcheting of driven attracting colloidal particles: Temporal density oscillations and current multiplicity

We consider the unidirectional particle transport in a suspension of colloidal particles which interact with each other via a pair potential having a hard-core repulsion plus an attractive tail. The colloids are confined within a long narrow channel and are driven along by a DC or an AC external potential. In addition, the walls of the channel interact with the particles via a ratchet-like periodic potential. We use dynamical density functional theory to compute the average particle current. In the case of DC drive, we show that as the attraction strength between the colloids is increased beyond a critical value, the stationary density distribution of the particles loses its stability leading to depinning and a time dependent density profile. Attraction induced symmetry breaking gives rise to the coexistence of stable stationary density profiles with different spatial periods and time-periodic density profiles, each characterized by different values for the particle current.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rayleigh and depinning instabilities of forced liquid ridges on heterogeneous substrates

Depinning of two-dimensional liquid ridges and three-dimensional drops on an inclined substrate is studied within the lubrication approximation. The structures are pinned to wetting heterogeneities arising from variations of the strength of the short-range polar contribution to the disjoining pressure. The case of a periodic array of hydrophobic stripes transverse to the slope is studied in detail using a combination of direct numerical simulation and branch-following techniques. Under appropriate conditions the ridges may either depin and slide downslope as the slope is increased, or first breakup into drops via a transverse instability, prior to depinning. The different transition scenarios are examined together with the stability properties of the different possible states of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Decomposition driven interface evolution for layers of binary mixtures: {II}. Influence of convective transport on linear stability

We study the linear stability with respect to lateral perturbations of free surface films of polymer mixtures on solid substrates. The study focuses on the stability properties of the stratified and homogeneous steady film states studied in Part I [U. Thiele, S. Madruga and L. Frastia, Phys. Fluids 19, 122106 (2007)]. To this aim, the linearized bulk equations and boundary equations are solved using continuation techniques for several different cases of energetic bias at the surfaces, corresponding to linear and quadratic solutal Marangoni effects. For purely diffusive transport, an increase in film thickness either exponentially decreases the lateral instability or entirely stabilizes the film. Including convective transport leads to a further destabilization as compared to the purely diffusive case. In some cases the inclusion of convective transport and the related widening of the range of available film configurations (it is then able to change its surface profile) change the stability behavior qualitatively. We furthermore present results regarding the dependence of the instability on several other parameters, namely, the Reynolds number, the Surface tension number and the ratio of the typical velocities of convective and diffusive transport.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical model for the formation of patterned deposits at receding contact lines

We describe the formation of deposition patterns that are observed in many different experiments where a three-phase contact line of a volatile nanoparticle suspension or polymer solution recedes. A dynamical model based on a long-wave approximation predicts the deposition of irregular and regular line patterns due to self-organised pinning-depinning cycles corresponding to a stick-slip motion of the contact line. We analyze how the line pattern properties depend on the evaporation rate and solute concentration.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the depinning of a drop of partially wetting liquid on a rotating cylinder

We discuss the analogy of the behaviour of films and drops of liquid on a rotating horizontal cylinder on the one hand and substrates with regular one-dimensional wettability patterns on the other hand. Based on the similarity between the respective governing long-wave equations we show that a drop of partially wetting liquid on a rotating cylinder undergoes a depinning transition when the rotation speed is increased. The transition occurs via a sniper bifurcation as in a recently described scenario for drops depinning on heterogeneous substrates.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thin film evolution equations from (evaporating) dewetting liquid layers to epitaxial growth

In the present contribution we review basic mathematical results for three physical systems involving self-organising solid or liquid films at solid surfaces. The films may undergo a structuring process by dewetting, evaporation/condensation or epitaxial growth, respectively. We highlight similarities and differences of the three systems based on the observation that in certain limits all of them may be described using models of similar form, i.e., time evolution equations for the film thickness profile. Those equations represent gradient dynamics characterized by mobility functions and an underlying energy functional. Two basic steps of mathematical analysis are used to compare the different system. First, we discuss the linear stability of homogeneous steady states, i.e., flat films; and second the systematics of non-trivial steady states, i.e., drop/hole states for dewetting films and quantum dot states in epitaxial growth, respectively. Our aim is to illustrate that the underlying solution structure might be very complex as in the case of epitaxial growth but can be better understood when comparing to the much simpler results for the dewetting liquid film. We furthermore show that the numerical continuation techniques employed can shed some light on this structure in a more convenient way than time-stepping methods. Finally we discuss that the usage of the employed general formulation does not only relate seemingly not related physical systems mathematically, but does as well allow to discuss model extensions in a more unified way.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Transport of free surface liquid films and drops by external ratchets and self-ratcheting mechanisms

We discuss the usage of ratchet mechanisms to transport a continuous phase in several micro-fluidic settings. In particular, we study the transport of a dielectric liquid in a heterogeneous ratchet capacitor that is periodically switched on and off. The second system consists of drops on a solid substrate that are transported by different types of harmonic substrate vibrations. We argue that the latter can be seen as a self-ratcheting process and discuss analogies between the employed class of thin film equations and Fokker-Planck equations for transport of discrete objects in a 'particle ratchet'.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Self-ratcheting Stokes drops driven by oblique vibrations

We develop and analyze a minimal hydrodynamic model in the overdamped limit to understand why a drop climbs a smooth homogeneous incline that is harmonically vibrated at an angle different from the substrate normal [Brunet, Eggers and Deegan, Phys. Rev. Lett. 99, 144501 (2007)]. We find that the vibration component orthogonal to the substrate induces a nonlinear (anharmonic) response in the drop shape. This results in an asymmetric response of the drop to the parallel vibration and, in consequence, in the observed net motion. Beside establishing the basic mechanism, we identify scaling laws valid in a broad frequency range and a flow reversal at high frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Decomposition driven interface evolution for layers of binary mixtures: I. Model derivation and stratified base states

A dynamical model is proposed to describe the coupled decomposition and profile evolution of a free surface film of a binary mixture. An example is a thin film of a polymer blend on a solid substrate undergoing simultaneous phase separation and dewetting. The model is based on model-H describing the coupled transport of the mass of one component (convective Cahn-Hilliard equation) and momentum (Navier-Stokes-Korteweg equations) supplemented by appropriate boundary conditions at the solid substrate and the free surface. General transport equations are derived using phenomenological non-equilibrium thermodynamics for a general non-isothermal setting taking into account Soret and Dufour effects and interfacial viscosity for the internal diffuse interface between the two components. Focusing on an isothermal setting the resulting model is compared to literature results and its base states corresponding to homogeneous or vertically stratified flat layers are analysed.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Drops on an inclined heterogeneous substrate: onset of sliding motion

Pinning and depinning of drops on an inclined heterogeneous substrate is studied as a function of the inclination and heterogeneity amplitude. Two types of heterogeneity are considered: a hydrophobic defect that blocks the droplet in front, and a hydrophilic one that holds it at the back. Two different types of depinning leading to sliding motion are identified, and the resulting stick-slip motion is studied numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Liquid transport generated by a flashing field-induced wettability ratchet

We develop and analyze a model for ratchet-driven macroscopic transport of a continuous phase. The transport relies on a field-induced dewetting-spreading cycle of a liquid film with a free surface based on a switchable, spatially asymmetric, periodic interaction of the liquid-gas interface and the substrate. The concept is exemplified using an evolution equation for a dielectric liquid film under an inhomogeneous voltage. We analyse the influence of the various phases of the ratchet cycle on the transport properties. Conditions for maximal transport and the efficiency of transport under load are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Asymptotic theory for a moving droplet driven by a wettability gradient

An asymptotic theory is developed for a moving drop driven by a wettability gradient. We distinguish the mesoscale where an exact solution is known for the properly simplified problem. This solution is matched at both -- the advancing and the receding side -- to respective solutions of the problem on the microscale. On the microscale the velocity of movement is used as the small parameter of an asymptotic expansion. Matching gives the droplet shape, velocity of movement as a function of the imposed wettability gradient and droplet volume.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Stratified spatiotemporal chaos in anisotropic reaction-diffusion systems

Numerical simulations of two dimensional pattern formation in an anisotropic bistable reaction-diffusion medium reveal a new dynamical state, stratified spatiotemporal chaos, characterized by strong correlations along one of the principal axes. Equations that describe the dependence of front motion on the angle illustrate the mechanism leading to stratified chaos.

0 citations0 reviewsNo code linkedOpen access