WIAS Preprint No. 2889, (2021)

An assessment of solvers for algebraically stabilized discretizations of convection-diffusion-reaction equations


Authors

  • Jha, Abhinav
  • Pártl, Ondřej
    ORCID: 0000-0002-1932-7172
  • Ahmed, Naveed
    ORCID: 0000-0002-9322-0373
  • Kuzmin, Dmitri

2020 Mathematics Subject Classification

  • 65M12 65M15 65M60

Keywords

  • finite element methods, discrete maximum principles, algebraic flux correction, flux-corrected transport, monolithic convex limiting, iterative solvers

DOI

10.20347/WIAS.PREPRINT.2889

Abstract

We consider flux-corrected finite element discretizations of 3D convection-dominated transport problems and assess the computational efficiency of algorithms based on such approximations. The methods under investigation include flux-corrected transport schemes and monolithic limiters. We discretize in space using a continuous Galerkin method and P1 or Q1 finite elements. Time integration is performed using the Crank-Nicolson method or an explicit strong stability preserving Runge-Kutta method. Nonlinear systems are solved using a fixed-point iteration method, which requires solution of large linear systems at each iteration or time step. The great variety of options in the choice of discretization methods and solver components calls for a dedicated comparative study of existing approaches. To perform such a study, we define new 3D test problems for time dependent and stationary convection-diffusion-reaction equations. The results of our numerical experiments illustrate how the limiting technique, time discretization and solver impact on the overall performance.

Appeared in

Download Documents