An Assessment of Solvers for Algebraically Stabilized Discretizations of Convection-Diffusion-Reaction Equations

by   Abhinav Jha, et al.

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 ℙ_1 or ℚ_1 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.



There are no comments yet.


page 12

page 26


Analysis of Flux Corrected Transport Schemes for Evolutionary Convection-Diffusion-Reaction Equations

We report in this paper the analysis for the linear and nonlinear versio...

Structure-Preserving and Efficient Numerical Methods for Ion Transport

Ion transport, often described by the Poisson–Nernst–Planck (PNP) equati...

Numerical simulation of multiscale fault systems with rate- and state-dependent friction

We consider the deformation of a geological structure with non-intersect...

Optimization of two-level methods for DG discretizations of reaction-diffusion equations

We analyze and optimize two-level methods applied to a symmetric interio...

Asymptotic preserving schemes for the FitzHugh-Nagumo transport equation with strong local interactions

This paper is devoted to the numerical approximation of the spatially-ex...

A note on the application of the Guermond-Pasquetti mass lumping correction technique for convection-diffusion problems

We provide a careful Fourier analysis of the Guermond-Pasquetti mass lum...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.