The comparative efficiency of time-splitting and multigrid schemes for the solution of lubrication models of three-dimensional (3D), thin film flow is demonstrated via detailed comparisons for benchmark gravity-driven continuous film and droplet spreading problems. Data are presented which show the effect of (1) problem formulation (either as a single fourth-order partial differential equation for the film thickness or two coupled second-order equations for film thickness and pressure) and (2) grid density on the choice of fixed time-step, CPU time per time-step and overall efficiency of each scheme. On the basis of these findings, recommendations are given as to the most efficient combination of problem formulation and numerical scheme for solving 3D thin film free surface flows over practical, engineering surfaces.
CrasterR. V.MatarO. K.Dynamics and stability of thin liquid films. Rev. Mod. Phys., 2009, 81(3), 1131–1198.
2.
HelbigK.NasarekR.Gambaryan-RoismanT.StephanP.Effect of longitudinal minigrooves on flow stability and wave characteristics of falling liquid films. J. Heat Transfer, 2009, 131(1), 011601–011601. (011608 pp.).
3.
HelbigK.AlexeevA.Gambaryan-RoismanT.StephanP.Evaporation of falling and shear-driven thin films on smooth and grooved surfaces. Flow Turbul. Combust., 2005, 75(1–4), 85–104.
4.
EtsionI.State of the art in laser surface texturing. J. Tribol., 2005, 127(1), 248–253.
5.
Kistler, S. F. and Schweizer, P. M. (eds) Liquid film coating, 1997 (Chapman and Hall, London).
DecréM.BaretJ.-C.Gravity-driven flows of viscous liquids over two-dimensional topographies. J. Fluid Mecha., 2003, 487, 147–166.
8.
PozrikidisC.ThoroddsenS. T.The deformation of a liquid film flowing down an inclined plane wall over a small particle arrested on the wall. Phys. Fluids A: Fluid Dynam., 1991, 3(11), 2546–2558. .
9.
BlythM. G.PozrikidisC.Film flow down an inclined plane over a three dimensional obstacle. Phys. Fluids, 2006, 18(5), 052104–052104.
10.
BaxterS. J.PowerH.CliffeK. A.HibberdS.Three-dimensional thin film flow over and around an obstacle on an inclined plane. Phys. Fluids, 2009, 21(3), 032102–032123.
11.
GaskellP. H.JimackP. K.SellierM.ThompsonH. M.Flow of evaporating, gravity-driven thin liquid films over topography. Phys. Fluids, 2006, 18, 1–1.
12.
SultanE.BoudaoudA.Ben AmarM.Evaporation of a thin film: diffusion of the vapour and Marangoni instabilities. J. Fluid Mech., 2005, 543, 183–202.
13.
SchwartzL. W.RoyR. V.EleyR. R.PrincenH. M.Surfactant-driven motion and splitting of droplets on a substrate. J. Eng. Math., 2004, 50(2–3), 157–175.
14.
VeremieievS.ThompsonH. M.LeeY. C.GaskellP. H.Computing inertial two- and three-dimensional thin film flow on planar surfaces featuring topography. Comput. Fluids, 2010, 39, 431–450.
15.
SaprykinS.KoopmansR. J.KalliadasisS. C.Free-surface thin-film flows over topography: influence of inertia and viscoelasticity. J. Fluid Mech., 2007, 578, 271–293.
16.
LenzR. D.KumarS.Competitive displacement of thin liquid films on chemically patterned substrates. J. Fluid Mech., 2007, 571, 33–57.
17.
LeeY. C.ThompsonH. M.GaskellP. H.Thin film flow over flexible membranes containing surface texturing: bio-inspired solutions. Proc. IMechE Part J: J. Engineering Tribology, 2009, 223, JET552–JET552.
18.
GaskellP. H.SavageM. D.SummersJ. L.ThompsonH. M.Modelling and analysis of meniscus roll coating. J. Fluid Mech., 1995, 298, 113–127.
19.
SummersJ. L.ThompsonH. M.GaskellP. H.Flow structure and transfer jets in a contra-rotating rigid-roll coating system. Theoret. Comput. Fluid Dynam., 2004, 17(3), 189–212.
SellierM.LeeY. C.ThompsonH. M.GaskellP. H.Thin film flow on surfaces containing localized arbitrary occlusions. Comput. Fluids, 2009, 38(1), 171–182.
22.
ChristovC. I.PontesJ.WalgraefD.VelardeM. G.Implicit time splitting for fourth-order parabolic equations. Comput. Methods Appl. Mech. Eng., 1997, 148(3–4), 209–224. .
23.
GaskellP. H.JimackP. K.SellierM.ThompsonH. M.Efficient and accurate time adaptive multigrid simulations of droplet spreading. Int. J. Numer. Methods Fluids, 2004, 45(11), 1161–1186. .
24.
GaskellP. H.JimackP. K.SellierM.ThompsonH. M.WilsonM. C. T.Gravity-driven flow of continuous thin liquid films on non-porous substrates with topography. J. Fluid Mech., 2004, 509, 253–280.
25.
KohY. Y.LeeY. C.GaskellP. H.JimackP. K.ThompsonH. M.Droplet spreading: quantitative comparisons with experiment. Eur. Phys. J., Special Topics, 2009, 166(1), 117–120.
26.
TannerL. H.The spreading of silicone oil drops on horizontal surfaces. J. Phys. D, 1979, 12, 1473–1484.
27.
StillwagonL. E.LarsonR. G.Leveling of thin films over uneven substrates during spin coating. Phys. Fluids A (Fluid Dynamics), 1990, 2(11), 1937–1944. .
28.
LeeY. C.ThompsonH. M.GaskellP. H.An efficient adaptive multigrid algorithm for predicting thin film flow on surfaces containing localised topographic features. Comput. Fluids, 2007, 36(5), 838–855. .
AltasI.DymJ.GuptaM. M.ManoharR. P.Multigrid solution of automatically generated high-order discretizations for the biharmonic equation. SIAM J. Sci. Comput., 1998, 19(5), 1575–1585.
32.
HennS.A Multigrid method for a fourth-order diffusion equation with application to image processing. SIAM J. Sci. Comput., 2005, 27(3), 831–849.
33.
LeeY. C.ThompsonH. M.GaskellP. H.FILMPAR: A parallel algorithm designed for the efficient and accurate computation of thin film flow on functional surfaces containing micro-structure. Comput. Physi. Commun., 2009, 180, 2634–2649.
34.
DouglasJ.Jr.A note on the alternating direction implicit method for the numerical solution of heat flow problems. Proc. Am. Math. Soc., 1955, 3(1), 28–41.
35.
DouglasJ.Jr.On the numerical integration of d2u/dx2+d2u/dy2=du/dt by implicit methods. J. Soc. Ind. Appl. Math., 1955, 3(1), 42–65.
36.
Daniels, N., Ehret, P., Gaskell, P. H., Thompson, H. M., and Decré, M. Multigrid methods for thin liquid film spreading flows. In Proceedings of the 1st International Conference on Computational Fluid Dynamics, 29 August 2001, pp. 279–284 (Springer, Kyoto, Japan).
