[1] 
G. Grün and M. Rumpf.
Nonnegativity preserving convergent schemes for the thin film
equation.
Numerische Mathematik, 87:113152, 2000. [ bib  .pdf 1 ] We present numerical schemes for fourth order degenerate parabolic equations that arise e.g. in lubrication theory for time evolution of thin films of viscous fluids. We prove convergence and nonnegativity results in arbitrary space dimensions. A proper choice of the discrete mobility enables us to establish discrete counterparts of the essential integral estimates known from the continuous setting. Hence, the numerical cost in each time step reduces to the solution of a linear system involving a sparse matrix. Furthermore, by introducing a time step control that makes use of an explicit formula for the normal velocity of the free boundary we keep the numerical cost for tracing the free boundary low.
