“…This model includes both Richards' model (with ζ (s) = s), which describes the flow of water in an underground medium, and Stefan's model (with β (s) = s), which arises in the study of the heat diffusion in a melting medium. The numerical approximation of both Richards' and Stefan's models has been extensively studied in the literature (see the fundamental work on the Stefan problem [13], and [14] for a review of some numerical approximations, and see [12] for the Richards problem), but the convergence analysis of the considered schemes received a much reduced coverage and consists mostly in establishing space-time averaged (e.g. in L 2 (Ω × (0, T ))) results (in the case of finite volume schemes, see for example [6,9]).…”