Finite element simulation of dynamic wetting flows as an interface formation process
KAUST Grant NumberKUK-C1-013-04
Permanent link to this recordhttp://hdl.handle.net/10754/598329
MetadataShow full item record
AbstractA mathematically challenging model of dynamic wetting as a process of interface formation has been, for the first time, fully incorporated into a numerical code based on the finite element method and applied, as a test case, to the problem of capillary rise. The motivation for this work comes from the fact that, as discovered experimentally more than a decade ago, the key variable in dynamic wetting flows - the dynamic contact angle - depends not just on the velocity of the three-phase contact line but on the entire flow field/geometry. Hence, to describe this effect, it becomes necessary to use the mathematical model that has this dependence as its integral part. A new physical effect, termed the 'hydrodynamic resist to dynamic wetting', is discovered where the influence of the capillary's radius on the dynamic contact angle, and hence on the global flow, is computed. The capabilities of the numerical framework are then demonstrated by comparing the results to experiments on the unsteady capillary rise, where excellent agreement is obtained. Practical recommendations on the spatial resolution required by the numerical scheme for a given set of non-dimensional similarity parameters are provided, and a comparison to asymptotic results available in limiting cases confirms that the code is converging to the correct solution. The appendix gives a user-friendly step-by-step guide specifying the entire implementation and allowing the reader to easily reproduce all presented results, including the benchmark calculations. © 2012 Elsevier Inc.
CitationSprittles JE, Shikhmurzaev YD (2013) Finite element simulation of dynamic wetting flows as an interface formation process. Journal of Computational Physics 233: 34–65. Available: http://dx.doi.org/10.1016/j.jcp.2012.07.018.
SponsorsThe authors would like to thank Dr Mark Wilson, Dr Paul Suckling and Dr Alex Lukyanov for many stimulating discussions about the FEM implementation of dynamic wetting phenomena and Jonathan Simmons for carefully proof reading the manuscript. JES kindly acknowledges the financial support of EPSRC via a Postdoctoral Fellowship in Mathematics.This publication was based on work supported in part by Award No KUK-C1-013-04, made by King Abdullah University of Science and Technology (KAUST).
JournalJournal of Computational Physics