An adaptive moving mesh method for thin film flow equations with surface tension

We present an adaptive moving mesh method for the numerical solution of thin liquid film spreading flows with surface tension. We follow the r-adaptive moving mesh technique which utilises a mesh density function and moving mesh partial differential equations (MMPDEs) to adapt and move the mesh coupled to the PDE(s) describing the thin film flow problem. Numerical experiments are performed on two one dimensional thin film flow equations to test the accuracy and efficiency of the method. This technique accurately resolves the multiple one-dimensional structures observed in these test problems. Moreover, it reduces the computational effort in comparison to the numerical solution using the finite difference scheme on a fixed uniform mesh. We are the first to implement r-adaptive schemes to higher order parabolic PDEs.This method accurately resolves the solution and reduces the computational effort.We introduce a new mesh density function that resolves multiple solution structures.This mesh density function is also adapted to include multiple solution components.

[1]  Tim G. Myers,et al.  Thin Films with High Surface Tension , 1998, SIAM Rev..

[2]  Richard V. Craster,et al.  Fingering phenomena created by a soluble surfactant deposition on a thin liquid film , 2004 .

[3]  Lou Kondic,et al.  Instabilities in Gravity Driven Flow of Thin Fluid Films , 2003, SIAM Rev..

[4]  Paul H. Muir,et al.  Algorithm 688: EPDCOL: a more efficient PDECOL code , 1991, TOMS.

[5]  M. Shearer,et al.  The Motion of a Thin Liquid Film Driven by Surfactant and Gravity , 2006, SIAM J. Appl. Math..

[6]  Linda R. Petzold,et al.  Observations on an adaptive moving grid method for one-dimensional systems of partial differential equations , 1987 .

[7]  S. Naire,et al.  The spreading and stability of a surfactant-laden drop on a prewetted substrate , 2006, Journal of Fluid Mechanics.

[8]  Philip H. Gaskell,et al.  The efficient and accurate solution of continuous thin film flow over surface patterning and past occlusions , 2008 .

[9]  Lou Kondic,et al.  Pattern formation in the flow of thin films down an incline: Constant flux configuration , 2001 .

[10]  Jinchao Xu,et al.  A new adaptive local mesh refinement algorithm and its application on fourth order thin film flow problem , 2007, J. Comput. Phys..

[11]  Joke Blom,et al.  An adaptive moving grid method for one-dimensional systems of partial differential equations , 1989 .

[12]  Robert D. Russell,et al.  Adaptive Moving Mesh Methods , 2010 .

[13]  Herbolzheimer,et al.  Model for the fingering instability of spreading surfactant drops. , 1990, Physical review letters.

[14]  R. Craster,et al.  Dynamics and stability of thin liquid films , 2009 .

[15]  Paul A. Zegeling,et al.  Algorithm 731: A moving-grid interface for systems of one-dimensional time-dependent partial differential equations , 1994, TOMS.

[16]  A. Bertozzi THE MATHEMATICS OF MOVING CONTACT LINES IN THIN LIQUID FILMS , 1998 .

[17]  R. Craster,et al.  Flow of surfactant-laden thin films down an inclined plane , 2004 .

[18]  R. W. Griffiths The Dynamics of Lava Flows , 2000 .

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

[20]  L. Petzold A description of dassl: a differential/algebraic system solver , 1982 .

[21]  Weizhang Huang,et al.  Moving Mesh Methods Based on Moving Mesh Partial Differential Equations , 1994 .

[22]  M. Shearer,et al.  Gravity-driven thin liquid films with insoluble surfactant: smooth traveling waves , 2007, European Journal of Applied Mathematics.

[23]  Thomas Y. Hou,et al.  An efficient dynamically adaptive mesh for potentially singular solutions , 2001 .

[24]  Robert D. Russell,et al.  Adaptivity with moving grids , 2009, Acta Numerica.

[25]  T. G. Myers,et al.  SURFACE TENSION DRIVEN THIN FILM FLOWS , 1996 .

[26]  R. Braun Dynamics of the Tear Film , 2012 .

[27]  Joke Blom,et al.  A moving grid method for one-dimensional PDEs based on the method of lines , 1988 .

[28]  J. Verwer,et al.  A numerical study of three moving-grid methods for one-dimensional partial differential equations which are based on the method of lines , 1990 .

[29]  R. Craster,et al.  Dynamics of a climbing surfactant-laden film II: stability. , 2012, Journal of colloid and interface science.

[30]  Weizhang Huang,et al.  Analysis Of Moving Mesh Partial Differential Equations With Spatial Smoothing , 1997 .

[31]  R. Craster,et al.  A note on the coating of an inclined plane in the presence of soluble surfactant. , 2006, Journal of colloid and interface science.

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

[33]  Martin Rumpf,et al.  Nonnegativity preserving convergent schemes for the thin film equation , 2000, Numerische Mathematik.

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

[35]  Chris J. Budd,et al.  Monge-Ampére based moving mesh methods for numerical weather prediction, with applications to the Eady problem , 2013, J. Comput. Phys..

[36]  S. Naire,et al.  The spreading and stability of a surfactant-laden drop on an inclined prewetted substrate , 2015, Journal of Fluid Mechanics.

[37]  Linda R. Petzold,et al.  Using Krylov Methods in the Solution of Large-Scale Differential-Algebraic Systems , 1994, SIAM J. Sci. Comput..

[38]  George Beckett,et al.  Convergence analysis of finite difference approximations on equidistributed grids to a singularly perturbed boundary value problem , 2000 .

[39]  Emily Jane Walsh Moving mesh methods for problems in meteorology , 2010 .

[40]  R. Craster,et al.  Fingering phenomena associated with insoluble surfactant spreading on thin liquid films , 2004, Journal of Fluid Mechanics.

[41]  Y. C. Lee,et al.  An efficient adaptive multigrid algorithm for predicting thin film flow on surfaces containing localised topographic features , 2007 .

[42]  Weizhang Huang,et al.  A moving collocation method for solving time dependent partial differential equations , 1996 .

[43]  R. Craster,et al.  Coating of an inclined plane in the presence of insoluble surfactant. , 2005, Journal of colloid and interface science.

[44]  R. Craster,et al.  Surfactant-induced fingering phenomena in thin film flow down an inclined plane , 2005 .

[45]  P E King-Smith,et al.  Single-equation models for the tear film in a blink cycle: realistic lid motion. , 2007, Mathematical medicine and biology : a journal of the IMA.

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

[47]  R. Craster,et al.  Surfactant-induced fingering phenomena beyond the critical micelle concentration , 2006, Journal of Fluid Mechanics.

[48]  J B Grotberg,et al.  Respiratory fluid mechanics and transport processes. , 2001, Annual review of biomedical engineering.

[49]  Yibao Li,et al.  Adaptive mesh refinement for simulation of thin film flows , 2014 .