The present invention is directed towards systems, methods and a computer-readable medium for simulating the evolution of a height of an evaporating droplet. The simulation includes a simulation space with boundary conditions. The simulation includes generating a height function that is representative of the height of the droplet at a first point in time at a plurality of points in the simulation space based upon a lubrication equation that is a differential function describing variation of the height function over time. The simulation determines the height function at a second point in time by finding an approximate solution that satisfies the lubrication equations and boundary conditions.