Convergence in C([0,T];L2(Ω)) of weak solutions to perturbed doubly degenerate parabolic equations

We study the behaviour of solutions to a class of nonlinear degenerate parabolic problems when the data are perturbed. The class includes the Richards equation, Stefan problem and the parabolic p-Laplace equation. We show that, up to a subsequence, weak solutions of the perturbed problem converge uniformly-in-time to weak solutions of the original problem as the perturbed data approach the original data. We do not assume uniqueness or regularity. When uniqueness is known, our result demonstrates that the weak solution is uniformly temporally stable to perturbations of the data. Beginning with a proof of temporally-uniform, spatially-weak convergence, we strengthen the latter by relating the unknown to an underlying convex structure that emerges naturally from energy estimates. The double degeneracy — shown to be equivalent to a maximal monotone operator framework — is handled with techniques inspired by a classical monotonicity argument and a simple variant of the compensated compactness phenomenon.