Conditional full stability of positivity-preserving finite difference scheme for diffusion-advection-reaction models

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:341

Article

Company, Rafael¹  Egorova, Vera N.²  Jodar, Lucas¹ 

[1] Univ Politecn Valencia, Camino Vera S-N, E-46022 Valencia, Spain

[2] BCAM, Al Mazarredo 14, Bilbao 48009, Basque Country, Spain

关键词: Diffusion-advection-reaction;    Semi-discretization;    Exponential time differencing;    Finite difference;    Numerical analysis;   

DOI  :  10.1016/j.cam.2018.02.031

来源: Elsevier

【 摘 要 】

The matter of the stability for multidimensional diffusion-advection-reaction problems treated with the semi-discretization method is remaining challenge because when all the stepsizes tend simultaneously to zero the involved size of the problem grows without bounds. Solution of such problems is constructed by starting with a semi-discretization approach followed by a full discretization using exponential time differencing and matrix quadrature rules. Analysis of the time variation of the numerical solution with respect to previous time level together with the use of logarithmic norm of matrices is the basis of the stability result. Sufficient stability conditions on stepsizes, that also guarantee positivity and boundedness of the solution, are found. Numerical examples in different fields prove its competitiveness with other relevant methods. (C) 2018 Elsevier B.V. All rights reserved.

【 授权许可】

Free