Development and application of a parallel multigrid solver for the simulation of spreading droplets

We consider the parallel application of an efficient solver developed for the accurate solution of a range of droplet spreading flows modelled as a coupled set of nonlinear lubrication equations. The underlying numerical scheme is based upon a second order finite difference discretization in space and a second order, fully implicit, adaptive scheme in time. At each time step this leads to the need to solve a large system of nonlinear algebraic equations, for which the full approximation storage (FAS) multigrid algorithm is employed. The motion of the contact line between the three phases (liquid, air and the solid substrate) is based upon the assumption of a thin precursor film, with a corresponding disjoining pressure term in the governing equations. It is the inclusion of this precursor film in the model that motivates the need for a parallel solution method. This is because the thickness of such a film must be very small in order to yield realistic predictions whilst the finite difference grid must be correspondingly fine in order to obtain accurate numerical solutions. Results are presented which demonstrate that the parallel implementation is sufficiently efficient and robust to allow reliable numerical solutions to be obtained for a level of mesh resolution that is an order of magnitude finer than is possible using a single processor. Copyright c 2007 John Wiley & Sons, Ltd.

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

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

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

[4]  L. Schwartz,et al.  Hysteretic Effects in Droplet Motions on Heterogeneous Substrates: Direct Numerical Simulation , 1998 .

[5]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[6]  Long Chen.  INTRODUCTION TO MULTIGRID METHODS , 2005 .

[7]  Martin Berzins,et al.  Parallelization and scalability issues of a multilevel elastohydrodynamic lubrication solver , 2007, Concurr. Comput. Pract. Exp..

[8]  Randolph E. Bank,et al.  A new parallel domain decomposition method for the adaptive finite element solution of elliptic partial differential equations , 2001, Concurr. Comput. Pract. Exp..

[9]  D. B. Graves,et al.  Spin coating over topography , 1993 .

[10]  P. Wesseling.  An Introduction to Multigrid Methods , 1992 .

[11]  Peter K. Jimack,et al.  Efficient and accurate time adaptive multigrid simulations of droplet spreading , 2004 .

[12]  Victor Starov,et al.  Spreading of liquid drops over dry surfaces , 1994 .

[13]  Peter K. Jimack,et al.  A fully implicit, fully adaptive time and space discretisation method for phase-field simulation of binary alloy solidification , 2007, J. Comput. Phys..

[14]  L. Schwartz,et al.  Simulation of Droplet Motion on Low-Energy and Heterogeneous Surfaces , 1998 .

[15]  Mathieu Sellier,et al.  The numerical simulation of thin film flow over heterogeneous substrates , 2003 .

[16]  A. Brandt.  Guide to multigrid development , 1982 .

[17]  S. Herminghaus,et al.  Wetting: Statics and dynamics , 1997 .