“…The well-posedness of (1.1) in the setting of classical solutions is addressed in [10], while nonnegative global weak solutions are constructed, by different approaches, in [3,8,9,13,14]. The rich dynamics described by the thin film Muskat problem (1.1) is very well illustrated by the laboratory experiments described in [23], but also by the numerical simulations reported in [1,2,15].…”