Abstract
In this paper a variant of nonlinear exponential Euler scheme is proposed for solving nonlinear heat conduction problems. The method is based on nonlinear iterations where at each iteration a linear initial-value problem has to be solved. We compare this method to the backward Euler method combined with nonlinear iterations. For both methods we show monotonicity and boundedness of the solutions and give sufficient conditions for convergence of the nonlinear iterations. Numerical tests are presented to examine performance of the two schemes. The presented exponential Euler scheme is implemented based on restarted Krylov subspace methods and, hence, is essentially explicit (involves only matrix-vector products).
Publisher
Keldysh Institute of Applied Mathematics
Subject
General Materials Science
Reference9 articles.
1. Hundsdorfer W., Verwer J. G. Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations. –– Springer Verlag, 2003.
2. Тихонов, А. Н., Самарский А. А. Уравнения математической физики. –– 6е изд. –– Изд-во МГУ, 1999.
3. Horn R. A., Johnson C. R. Topics in Matrix Analysis. –– Cambridge University Press, 1991.
4. Hochbruck M., Ostermann A. Exponential integrators // Acta Numer. –– 2010. –– Vol. 19. –– P. 209–286.
5. Dekker K., Verwer J. G. Stability of Runge–Kutta methods for stiff non-linear differential equations. –– North-Holland Elsevier Science Publishers, 1984.