Iterative linearisation schemes for doubly degenerate parabolic equations

From MaRDI portal
Publication:2008652

DOI10.1007/978-3-319-96415-7_3zbMATH Open1427.65237arXiv1801.00846OpenAlexW2908196371MaRDI QIDQ2008652FDOQ2008652


Authors: Jakub Wiktor Both, K. Kumar, J. M. Nordbotten, I. S. Pop, F. A. Radu Edit this on Wikidata


Publication date: 26 November 2019

Abstract: Mathematical models for flow and reactive transport in porous media often involve non-linear, degenerate parabolic equations. Their solutions have low regularity, and therefore lower order schemes are used for the numerical approximation. Here the backward Euler method is combined with a mixed finite element method scheme, which results in a stable and locally mass-conservative scheme. At the same time, at each time step one has to solve a non-linear algebraic system, for which linear iterations are needed. Finding robust and convergent ones is particularly challenging here, since both slow and fast diffusion cases are allowed. Commonly used schemes, like Newton and Picard iterations, are defined either for non-degenerate problems, or after regularising the problem in the case of degenerate ones. Convergence is guaranteed only if the initial guess is sufficiently close to the solution, which translates into severe restrictions on the time step. Here we discuss a linear iterative scheme which builds on the L-scheme, and does not employ any regularisation. We prove its rigourous convergence, which is obtained for mild restrictions on the time step. Finally, we give numerical results confirming the theoretical ones, and compare the behaviour of the scheme with other schemes.


Full work available at URL: https://arxiv.org/abs/1801.00846




Recommendations




Cited In (5)





This page was built for publication: Iterative linearisation schemes for doubly degenerate parabolic equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2008652)