Time Integration and Steady-State Continuation for 2d Lubrication Equations

Lubrication equations describe many structuring processes of thin liquid films. We develop and apply a numerical framework 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. Both algorithms employ a Cayley transform to overcome numerical problems resulting from scale separation in space and time. An adaptive time-step allows one to study the dynamics close to hetero- or homoclinic connections. The developed framework is employed, on the one hand, to analyze 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 analyze 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...

[1]  Ioannis G. Kevrekidis,et al.  Bifurcations and pattern formation in the "regularized" Kuramoto-Sivashinsky equation , 1992 .

[2]  I. Kevrekidis,et al.  Back in the saddle again: a computer assisted study of the Kuramoto-Sivashinsky equation , 1990 .

[3]  P. Rosenau,et al.  On a nonlinear thermocapillary effect in thin liquid layers , 1994, Journal of Fluid Mechanics.

[4]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[5]  S. Chou,et al.  Novel polymer patterns formed by lithographically induced self-assembly (LISA). , 2005, Langmuir : the ACS journal of surfaces and colloids.

[6]  U. Thiele,et al.  Long-wave theory of bounded two-layer films with a free liquid–liquid interface: Short- and long-time evolution , 2004, nlin/0412040.

[7]  J. Vega,et al.  Control of Rayleigh-Taylor instability by vertical vibration in large aspect ratio containers. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  Y. Saad.  Analysis of some Krylov subspace approximations to the matrix exponential operator , 1992 .

[9]  Ashutosh Sharma,et al.  Instability of thin polymer films on coated substrates : Rupture, dewetting, and drop formation , 1996 .

[10]  U. Thiele,et al.  Morphology changes in the evolution of liquid two-layer films. , 2004, The Journal of chemical physics.

[11]  L. Tuckerman,et al.  A method for exponential propagation of large systems of stiff nonlinear differential equations , 1989 .

[12]  Self-propelled running droplets on solid substrates driven by chemical reactions , 2005, The European physical journal. E, Soft matter.

[13]  L. Tuckerman,et al.  Bifurcation Analysis for Timesteppers , 2000 .

[14]  M. Golubitsky,et al.  Bifurcation on the hexagonal lattice and the planar Bénard problem , 1983, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[15]  S. Troian,et al.  A Theoretical Study of Instabilities at the Advancing Front of Thermally Driven Coating Films , 1997, Journal of colloid and interface science.

[16]  Andrea L. Bertozzi,et al.  CONTACT LINE STABILITY AND UNDERCOMPRESSIVE SHOCKS IN DRIVEN THIN FILM FLOW , 1998 .

[17]  T. P. Witelskia,et al.  ADI schemes for higher-order nonlinear diffusion equations , 2003 .

[18]  Unther Gr¨un.  ON THE CONVERGENCE OF ENTROPY CONSISTENT SCHEMES FOR LUBRICATION TYPE EQUATIONS IN MULTIPLE SPACE DIMENSIONS , 2000 .

[19]  Herbert E. Huppert,et al.  Flow and instability of a viscous current down a slope , 1982, Nature.

[20]  A. Münch.  Pinch-off transition in Marangoni-driven thin films. , 2003, Physical review letters.

[21]  E. Knobloch,et al.  Symmetry and Symmetry-Breaking Bifurcations in Fluid Dynamics , 1991 .

[22]  Serafim Kalliadasis.  Nonlinear instability of a contact line driven by gravity , 2000, Journal of Fluid Mechanics.

[23]  Stephan Herminghaus,et al.  Dynamics and structure formation in thin polymer melt films , 2005 .

[24]  M. J. Robbins,et al.  Modelling approaches to the dewetting of evaporating thin films of nanoparticle suspensions , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[25]  U. Thiele,et al.  Dynamics of a horizontal thin liquid film in the presence of reactive surfactants , 2007 .

[26]  Driven drops on heterogeneous substrates: onset of sliding motion. , 2006, Physical review letters.

[27]  P. Hänggi,et al.  Depinning of three-dimensional drops from wettability defects , 2008, 0811.2918.

[28]  George M. Homsy,et al.  Stability of Newtonian and viscoelastic dynamic contact lines , 1996 .

[29]  Ashutosh Sharma.  Relationship of thin film stability and morphology to macroscopic parameters of wetting in the apolar and polar systems , 1993 .

[30]  Eli Ruckenstein,et al.  Spontaneous rupture of thin liquid films , 1974 .

[31]  Eusebius J. Doedel,et al.  Numerical Methods for Bifurcation Problems and Large-Scale Dynamical Systems , 2012 .

[32]  Andrea L. Bertozzi,et al.  Undercompressive shocks in thin film flows , 1999 .

[33]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[34]  P. Hänggi,et al.  PAPER www.rsc.org/softmatter | Soft Matter Ratchet-driven fluid transport in bounded two-layer films of immiscible liquids , 2007 .

[35]  Ashutosh Sharma,et al.  PATTERN FORMATION IN UNSTABLE THIN LIQUID FILMS , 1998 .

[36]  Paul Wheeler,et al.  Computation of Spiral Spectra , 2005, SIAM J. Appl. Dyn. Syst..

[37]  Thomas Podgorski,et al.  Corners, cusps, and pearls in running drops. , 2001, Physical review letters.

[38]  Ralf Blossey,et al.  Complex dewetting scenarios captured by thin-film models , 2003, Nature materials.

[39]  S. Bankoff,et al.  Long-scale evolution of thin liquid films , 1997 .

[40]  U. Thiele.  Open questions and promising new fields in dewetting , 2003, The European physical journal. E, Soft matter.

[41]  D. Bonn,et al.  Wetting and Spreading , 2009 .

[42]  M. Bestehorn,et al.  Surface patterns of laterally extended thin liquid films in three dimensions. , 2001, Physical review letters.

[43]  T. Manteuffel.  The Tchebychev iteration for nonsymmetric linear systems , 1977 .

[44]  R. Seydel.  Practical Bifurcation and Stability Analysis , 1994 .

[45]  Y. Pomeau,et al.  Film rupture in the diffuse interface model coupled to hydrodynamics. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[46]  Andrea L. Bertozzi,et al.  Linear stability and transient growth in driven contact lines , 1997 .

[47]  Andreas Münch,et al.  Dewetting rates of thin liquid films , 2005 .

[48]  U. Thiele.  Structure Formation in Thin Liquid Films , 2007 .

[49]  Diem Ho.  Tchebychev acceleration technique for large scale nonsymmetric matrices , 1989 .

[50]  L. Pismen,et al.  Marangoni waves induced by a multistable chemical reaction on thin liquid films , 1984 .

[51]  L. Brusch,et al.  Modelling Thin Film Dewetting on Structured Surfaces and Templates: Bifurcation Analysis and Numerical Simulations , 2022 .

[52]  M. Sadkane.  A block Arnoldi-Chebyshev method for computing the leading eigenpairs of large sparse unsymmetric matrices , 1993 .

[53]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[54]  E. Knobloch,et al.  On the depinning of a driven drop on a heterogeneous substrate , 2006 .

[55]  Lou Kondic,et al.  Computing three-dimensional thin film flows including contact lines , 2002 .

[56]  P. Rosenau,et al.  Formation of patterns induced by thermocapillarity and gravity , 1992 .

[57]  É. Raphaël,et al.  Dewetting of thin polymer films , 2006, The European physical journal. E, Soft matter.

[58]  Andrea L. Bertozzi,et al.  Shock Solutions for Particle-Laden Thin Films , 2005, SIAM J. Appl. Math..

[59]  U. Thiele.  Thin film evolution equations from (evaporating) dewetting liquid layers to epitaxial growth , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[60]  L. Limat,et al.  Shape and motion of drops sliding down an inclined plane , 2005, Journal of Fluid Mechanics.

[61]  Martin Rumpf,et al.  Simulation of singularities and instabilities arising in thin film flow , 2001, European Journal of Applied Mathematics.

[62]  Oliver E. Jensen,et al.  Insoluble surfactant spreading on a thin viscous film: shock evolution and film rupture , 1992, Journal of Fluid Mechanics.

[63]  F. Heslot,et al.  Fingering instability of thin spreading films driven by temperature gradients , 1990, Nature.

[64]  Ashutosh Sharma,et al.  Electric-field-induced interfacial instabilities and morphologies of thin viscous and elastic bilayers. , 2009, Langmuir : the ACS journal of surfaces and colloids.

[65]  Oron,et al.  Stable localized patterns in thin liquid films. , 1992, Physical review letters.

[66]  Dynamical model for chemically driven running droplets. , 2004, Physical review letters.

[67]  U. Thiele,et al.  Thin Films of Soft Matter , 2007 .

[68]  U. Thiele,et al.  3D Large scale Marangoni convection in liquid films , 2003 .

[69]  A. Sharma,et al.  Auto-optimization of dewetting rates by rim instabilities in slipping polymer films. , 2001, Physical review letters.

[70]  E. Knobloch,et al.  Long-wave Marangoni instability with vibration , 2005, Journal of Fluid Mechanics.

[71]  Ashutosh Sharma,et al.  The ABC of pattern evolution in self-destruction of thin polymer films , 2000 .

[72]  Edgar Knobloch,et al.  Thin liquid films on a slightly inclined heated plate , 2004 .

[73]  L. Schwartz,et al.  Fingering phenomena for driven coating films , 2000 .

[74]  Reiter,et al.  Dewetting of thin polymer films. , 1992, Physical review letters.

[75]  Alexander Oron,et al.  Nonlinear dynamics of three-dimensional long-wave Marangoni instability in thin liquid films , 2000 .

[76]  R. Craster,et al.  Models for Marangoni drying , 2001 .

[77]  J. H. Conway,et al.  Groups, Combinatorics & Geometry: The orbifold notation for surface groups , 1992 .

[78]  M. Velarde,et al.  Dewetting: film rupture by nucleation in the spinodal regime. , 2001, Physical review letters.

[79]  V. Mitlin.  Dewetting of Solid Surface: Analogy with Spinodal Decomposition , 1993 .

[80]  Ashutosh Sharma,et al.  Electric field induced instability and pattern formation in thin liquid films. , 2005, Langmuir : the ACS journal of surfaces and colloids.

[81]  Jörg Liesen,et al.  Convergence analysis of Krylov subspace methods , 2004 .

[82]  Andrea L. Bertozzi,et al.  Dewetting films: bifurcations and concentrations , 2001 .

[83]  Antonio Palacios,et al.  Noise-Induced intermittency in Cellular Pattern-Forming Systems , 2007, Int. J. Bifurc. Chaos.

[84]  J. Swift,et al.  Suppression of dripping from a ceiling. , 2001, Physical review letters.

[85]  J. Wesfreid,et al.  Two-dimensional patterns in Rayleigh-Taylor instability of a thin layer , 1992, Journal of Fluid Mechanics.

[86]  A. Oron,et al.  Three-dimensional nonlinear dynamics of thin liquid films , 2000, Physical review letters.

[87]  K. Meerbergen,et al.  Matrix transformations for computing rightmost eigenvalues of large sparse non-symmetric eigenvalue problems , 1996 .

[88]  Andrea L. Bertozzi,et al.  Positivity-Preserving Numerical Schemes for Lubrication-Type Equations , 1999, SIAM J. Numer. Anal..

[89]  Zhiqun Lin,et al.  Electric field induced instabilities at liquid'liquid interfaces , 2001 .

[90]  L. Limat,et al.  Self-similar flow and contact line geometry at the rear of cornered drops , 2005 .

[91]  E. Knobloch,et al.  Front and back instability of a liquid film on a slightly inclined plate , 2003 .

[92]  Thomas P. Witelski,et al.  Coarsening dynamics of dewetting films. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[93]  A. Nepomnyashchy,et al.  Convective Cahn-Hilliard models: from coarsening to roughening. , 2001, Physical review letters.

[94]  R. Behringer,et al.  Reverse undercompressive shock structures in driven thin film flow. , 2002, Physical review letters.

[95]  M. Fermigier,et al.  Rayleigh-Taylor Instability of a Thin Layer , 1991 .

[96]  Y. Pomeau,et al.  Sliding drops in the diffuse interface model coupled to hydrodynamics. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.