Numerical studies of the fingering phenomena for the thin film equation

• Published in: International journal for numerical methods in fluids 2011-12, Vol.67 (11), p.1358-1372
• Main Authors: Li, Yibao, Lee, Hyun Geun, Yoon, Daeki, Hwang, Woonjae, Shin, Suyeon, Ha, Youngsoo, Kim, Junseok
• Format: Article
• Language: English
• Subjects: Applied fluid mechanics; Computational methods in fluid dynamics; Exact sciences and technology; Film thickness; fingering instability; finite difference; Finite difference method; Fluid dynamics; Fundamental areas of phenomenology (including applications); Heat treatment; Hydrodynamics, hydraulics, hydrostatics; Liquid films; Marangoni stress; Mathematical analysis; nonlinear diffusion equation; nonlinear multigrid method; Nonlinearity; Numerical analysis; Physics; thin film; Thin films
• ISSN: 0271-2091; 1097-0363; 1097-0363
• DOI: 10.1002/fld.2420

We present a new interpretation of the fingering phenomena of the thin liquid film layer through numerical investigations. The governing partial differential equation is ht + (h2−h3)x = −∇·(h3∇Δh), which arises in the context of thin liquid films driven by a thermal gradient with a counteracting gravitational force, where h = h(x, y, t) is the liquid film height. A robust and accurate finite difference method is developed for the thin liquid film equation. For the advection part (h2−h3)x, we use an implicit essentially non‐oscillatory (ENO)‐type scheme and get a good stability property. For the diffusion part −∇·(h3∇Δh), we use an implicit Euler's method. The resulting nonlinear discrete system is solved by an efficient nonlinear multigrid method. Numerical experiments indicate that higher the film thickness, the faster the film front evolves. The concave front has higher film thickness than the convex front. Therefore, the concave front has higher speed than the convex front and this leads to the fingering phenomena. Copyright © 2010 John Wiley & Sons, Ltd.