A lagged diffusivity method for reaction-convection-diffusion equations with Dirichlet boundary conditions

In this paper we solve a 2D nonlinear, non-steady reaction–convection–diffusion equation subject to Dirichlet boundary conditions by an iterative procedure consisting in lagging the diffusion term. First, we analyze the procedure, which we call Lagged Diffusivity Method. In particular, we provide a proof of the uniqueness of the solution and of the convergence of the lagged iteration when some assumptions are satisfied. We also describe outer and inner solvers, with special regard to how to link the stopping criteria in an efficient way. Numerical experiments are then introduced in order to evaluate the role of different linear solvers and of other components of the solution procedure, considering also the effect of the discretization.