Numerical simulation of non-isothermal thin liquid film flow on inclined plane using an implicit finite difference scheme

The classical problem of the stability and dynamics of thin liquid films on solid surfaces has been studied extensively. Particularly, thin liquid films subjected to various physico-chemical effects such as thermocapillarity, solutal-Marangoni and evaporative instabilities at the film surface has been the focus of research for more than two decades. Various flow configurations of thin film such as thin film on plane, inclined, and wavy surfaces has been the subject of recent investigations. An inclined film compared to a horizontal film, also experiences the gravity force which may significantly influence the nonlinear dynamics of the film coupled with other forces. In this research, we attempt to study the stability and dynamics of thin liquid films subjected to thermocapillarity and evaporative instabilities at the free surface besides instability owing to ubiquitous van der Waals attraction, using numerical simulations. For a Newtonian liquid, flow in thin liquid film on a planar support and bounded by a passive gas, is represented by Navier-Stokes equation, equation of continuity and appropriate boundary conditions. The external effects are incorporated in the body force term of the Navier-Stokes equation. These governing equations are simplified using the so called long-wave approximation to arrive at a nonlinear partial differential equation, henceforth called equation of evolution (EOE), which describes the time evolution of the interfacial instability in the film caused by internal and/ or external effects. Efficient numerical method is required for the solution of the equation of evolution (EOE) in order to comprehend the nonlinear dynamics of the thin film. Here we present the results of our numerical simulation using Crank-Nicholson implicit finite difference scheme applied to the thin film model incorporating instabilities owing to gravity, evaporation and thermo-capillarity. Comparison of our results with those obtained from Spectral method, show remarkable agreement for most of the cases investigated.