On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method

Abstract

We present a robust computational framework for advective-diffusive-reactive systems that satisfies maximum principles, the non-negative constraint, and element-wise species balance property. The proposed methodology is valid on general computational grids, can handle heterogeneous anisotropic media, and provides accurate numerical solutions even for very high Péclet numbers. The significant contribution of this paper is to incorporate advection (which makes the spatial part of the differential operator non-self-adjoint) into the non-negative computational framework, and overcome numerical challenges associated with advection. We employ low-order mixed finite element formulations based on least-squares formalism, and enforce explicit constraints on the discrete problem to meet the desired properties. The resulting constrained discrete problem belongs to convex quadratic programming for which a unique solution exists. Maximum principles and the non-negative constraint give rise to bound constraints while element-wise species balance gives rise to equality constraints. The resulting convex quadratic programming problems are solved using an interior-point algorithm. Several numerical results pertaining to advection-dominated problems are presented to illustrate the robustness, convergence, and the overall performance of the proposed computational framework.

References in corpus (5)

• A numerical framework for diffusion-controlled bimolecular-reactive systems to enforce maximum principles and non-negative constraint
• A Note on why Enforcing Discrete Maximum Principles by a simple a Posteriori Cutoff is a Good Idea
• Numerical methods for convection-diffusion problems or The 30 years war
• Sign-preserving of principal eigenfunctions in P1 finite element approximation of eigenvalue problems of second-order elliptic operators
• A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations