An adaptive, space-time discretized linear iterative scheme for doubly-degenerate parabolic problems
For researchers simulating degenerate diffusion problems with free boundaries, this work provides a fully adaptive solver with rigorous error control, though the approach is incremental over existing linearization and adaptivity techniques.
The paper presents an adaptive space-time solver for doubly-degenerate parabolic problems that uses an iterative linearization scheme reducing the problem to a sequence of heat equations. Numerical experiments show efficient allocation of computational resources with rapid error decay in terms of total degrees of freedom.
Degenerate diffusion problems, where the governing parabolic equation can change type to either an ordinary differential equation or an elliptic equation, model many real life applications. Due to the presence of free-boundaries, accurate numerical simulation of such problems require extremely small mesh and time step sizes locally. To remediate this issue, in this work, we consider a space-time formulation of the problem based on an efficient splitting of the nonlinearities. First, an iterative linearization scheme is proposed to resolve the nonlinearities that effectively reduces to solving a sequence of heat equations. Unconditional convergence of the scheme is proven even for double degenerate cases with linear convergence achieved if the problem is non-degenerate. Next, the dual norm of the nonlinear residual is decomposed into a linearization error component and a discretization error component corresponding to the heat equation. This leads to reliable and fully computable a posteriori estimates for the problem that are robust with respect to the nonlinearities/degeneracies. These estimates are used then in a fully adaptive (discretization + linearization) space-time solver. Numerical experiments for multiple test cases (one and two dimensions in space) demonstrate that this solver efficiently allocates the computational resources in the space-time domain, resulting in a rapid decay of error in terms of total degrees of freedom spent.